Search

We propose a one dimensional fully nonlinear model for surface waves in channels accounting for dispersive effects strictly related to multi-dimensional hydrostatic processes related to geometrical variations, namely to the bathymetry. Our model generalizes and improves those proposed by Chassagne et al. 2018 and 2019, as well as Quezada de Luna et al JFM 2021. The new model is obtained as an approximation of the hyperbolic shallow water equations under some hypotheses on the anisotropy of the propagation. For historical reasons the new model is referred to as geometrical Green-Naghdi type (gGN) system. The gGN system is fully nonlinear and weakly dispersive. It is endowed with an exact energy conservation law, and admits several types of exact travelling wave solutions. We analyze the propagation properties compare them to those of Chassagne et al. JFM 2019, including a characterization of undular bores via the analogy by Lemoine. Finally the system is recast in two useful forms appropriate for its numerical approximations, whose properties are discussed. Our results allow to further confirm the existence of these dispersive processes (highlighted in Chassagne et al JFM 2019) having a fully hydrostatic nature and modelled by the 2D hyperbolic nonlinear shallow water simulations. Finally, the numerical fully validate the new gGN model proposed.

We develop structure-preserving numerical methods for the Serre-Green-Naghdi equations, a model for weakly dispersive free-surface waves. We consider both the classical form, requiring the inversion of a non-linear elliptic operator, and a hyperbolic approximation of the equations, allowing fully explicit time stepping. Systems for both flat and variable topography are studied. Our novel numerical methods conserve both the total water mass and the total energy. In addition, the methods for the original Serre-Green-Naghdi equations conserve the total momentum for flat bathymetry. For variable topography, all the methods proposed are well-balanced for the lake-at-rest state. We provide a theoretical setting allowing us to construct schemes of any kind (finite difference, finite element, discontinuous Galerkin, spectral, etc.) as long as summation-by-parts operators are available in the chosen setting. Energy-stable variants are proposed by adding a consistent high-order artificial viscosity term. The proposed methods are validated through a large set of benchmarks to verify all the theoretical properties. Whenever possible, comparisons with exact, reference numerical, or experimental data are carried out. The impressive advantage of structure preservation, and in particular energy preservation, to resolve accurately dispersive wave propagation on very coarse meshes is demonstrated by several of the tests.

Numerical methods for hyperbolic PDEs require stabilization. For linear acoustics, divergence-free vector fields should remain stationary, but classical Finite Difference methods add incompatible diffusion that dramatically restricts the set of discrete stationary states of the numerical method. Compatible diffusion should vanish on stationary states, e.g. should be a gradient of the divergence. Some Finite Element methods allow to naturally embed this grad-div structure, e.g. the SUPG method or OSS. We prove here that the particular discretization associated to them still fails to be constraint preserving. We then introduce a new framework on Cartesian grids based on surface (volume in 3D) integrated operators inspired by Global Flux quadrature and related to mimetic approaches. We are able to construct constraint-compatible stabilization operators (e.g. of SUPG-type) and show that the resulting methods are vorticity-preserving. We show that the Global Flux approach is even super-convergent on stationary states, we characterize the kernels of the discrete operators and we provide projections onto them.

We present a novel quasi-conservative arbitrary high order accurate ADER discontinuous Galerkin (DG) method allowing to efficiently use a non-conservative form of the considered partial differential system, so that the governing equations can be solved directly in the most physically relevant set of variables. This is particularly interesting for multi-material flows with moving interfaces and steep, large magnitude contact discontinuities, as well as in presence of highly non-linear thermodynamics. However, the non-conservative formulation of course introduces a conservation error which would normally lead to a wrong approximation of shock waves. Hence, from the theoretical point of view, we give a formal definition of the conservation defect of non-conservative schemes and we analyze this defect providing a local quasi-conservation condition, which allows us to prove a modified Lax-Wendroff theorem. Then, to deal with shock waves in practice, we exploit the framework of the so-called a posteriori subcell finite volume (FV) limiter, so that, in troubled cells appropriately detected, we can incorporate a local conservation correction. Our corrected FV update entirely removes the local conservation defect, allowing, at least formally, to fit in the hypotheses of the proposed modified Lax-Wendroff theorem. Here, the shock-triggered troubled cells are detected by combining physical admissibility criteria, a discrete maximum principle and a shock sensor inspired by Lagrangian hydrodynamics. To prove the capabilities of our novel approach, first, we show that we are able to recover the same results given by conservative schemes on classical benchmarks for the single-fluid Euler equations. We then conclude by showing the improved reliability of our scheme on the multi-fluid Euler system on examples like the interaction of a shock with a helium bubble.

We propose a novel approach to approximate numerically shock waves. The method combines the unstructured shock-fitting approach developed in the last decade by some of the authors, with ideas coming from embedded boundary techniques. The numerical method obtained allows avoiding the re-meshing phase required by the unstructured fitting method, while guaranteeing accuracy properties very close to those of the fitting approach. This new method has many similarities with front tracking approaches, and paves the way to shock-tracking techniques truly independent on the data and mesh structure used by the flow solver. The approach is tested on several problems showing accuracy properties very close to those of more expensive fitting methods, with a considerable gain in flexibility and generality.

In this paper we present a novel approach for the design of high order general boundary conditions when approximating solutions of the Euler equations on domains with curved boundaries, using meshes which may not be boundary conformal. When dealing with curved boundaries and/or unfitted discretizations, the consistency of boundary conditions is a well-known challenge, especially in the context of high order schemes. In order to tackle such consistency problems, the so-called Reconstruction for Off-site Data (ROD) method has been recently introduced in the finite volume framework: it is based on performing a boundary polynomial reconstruction that embeds the considered boundary treatment thanks to the implementation of a constrained minimization problem. This work is devoted to the development of the ROD approach in the context of discontinuous finite elements. We use the genuine space-time nature of the local ADER predictors to reformulate the ROD as a single space-time reconstruction procedure. This allows us to avoid a new reconstruction (linear system inversion) at each sub-time node and retrieve a single space-time polynomial that embeds the considered boundary conditions for the entire space-time element. Several numerical experiments are presented proving the consistency of the new approach for all kinds of boundary conditions. Computations involving the interaction of shocks with embedded curved boundaries are made possible through an a posteriori limiting technique.

We present a new high-order accurate spectral element solution to the two-dimensional scalar Poisson equation subject to a general Robin boundary condition. The solution is based on a simplified version of the shifted boundary method employing a continuous arbitrary order $hp$-Galerkin spectral element method as the numerical discretization procedure. The simplification relies on a polynomial correction to avoid explicitly evaluating high-order partial derivatives from the Taylor series expansion, which traditionally have been used within the shifted boundary method. In this setting, we apply an extrapolation and novel interpolation approach to project the basis functions from the true domain onto the approximate surrogate domain. The resulting solution provides a method that naturally incorporates curved geometrical features of the domain, overcomes complex and cumbersome mesh generation, and avoids problems with small-cut-cells. Dirichlet, Neumann, and general Robin boundary conditions are enforced weakly through: i) a generalized Nitsche's method and ii) a generalized Aubin's method. For this, a consistent asymptotic preserving formulation of the embedded Robin formulations is presented. We present several numerical experiments and analysis of the algorithmic properties of the different weak formulations. With this, we include convergence studies under polynomial, $p$, increase of the basis functions, mesh, $h$, refinement, and matrix conditioning to highlight the spectral and algebraic convergence features, respectively. This is done to assess the influence of errors across variational formulations, polynomial order, mesh size, and mappings between the true and surrogate boundaries., Comment: 1 table, 18 figures

In this work, the high order accuracy and the well-balanced (WB) properties of some novel continuous interior penalty (CIP) stabilizations for the Shallow Water (SW) equations are investigated. The underlying arbitrary high order numerical framework is given by a Residual Distribution (RD)/continuous Galerkin (CG) finite element method (FEM) setting for the space discretization coupled with a Deferred Correction (DeC) time integration, to have a fully-explicit scheme. If, on the one hand, the introduced CIP stabilizations are all specifically designed to guarantee the exact preservation of the lake at rest steady state, on the other hand, some of them make use of general structures to tackle the preservation of general steady states, whose explicit analytical expression is not known. Several basis functions have been considered in the numerical experiments and, in all cases, the numerical results confirm the high order accuracy and the ability of the novel stabilizations to exactly preserve the lake at rest steady state and to capture small perturbations of such equilibrium. Moreover, some of them, based on the notions of space residual and global flux, have shown very good performances and superconvergences in the context of general steady solutions not known in closed-form. Many elements introduced here can be extended to other hyperbolic systems, e.g., to the Euler equations with gravity.

We present a novel approach for solving the shallow water equations using a discontinuous Galerkin spectral element method. The method we propose has three main features. First, it enjoys a discrete well-balanced property, in a spirit similar to the one of e.g. [20]. As in the reference, our scheme does not require any a-priori knowledge of the steady equilibrium, moreover it does not involve the explicit solution of any local auxiliary problem to approximate such equilibrium. The scheme is also arbitrarily high order, and verifies a continuous in time cell entropy equality. The latter becomes an inequality as soon as additional dissipation is added to the method. The method is constructed starting from a global flux approach in which an additional flux term is constructed as the primitive of the source. We show that, in the context of nodal spectral finite elements, this can be translated into a simple modification of the integral of the source term. We prove that, when using Gauss-Lobatto nodal finite elements this modified integration is equivalent at steady state to a high order Gauss collocation method applied to an ODE for the flux. This method is superconvergent at the collocation points, thus providing a discrete well-balanced property very similar in spirit to the one proposed in [20], albeit not needing the explicit computation of a local approximation of the steady state. To control the entropy production, we introduce artificial viscosity corrections at the cell level and incorporate them into the scheme. We provide theoretical and numerical characterizations of the accuracy and equilibrium preservation of these corrections. Through extensive numerical benchmarking, we validate our theoretical predictions, with considerable improvements in accuracy for steady states, as well as enhanced robustness for more complex scenarios

In this work we propose a simple but effective high order polynomial correction allowing to enhance the consistency of all kind of boundary conditions for the Euler equations (Dirichlet, characteristic far-field and slip-wall), both in 2D and 3D, preserving a high order of accuracy without the need of curved meshes. The method proposed is a simplified reformulation of the Shifted Boundary Method (SBM) and relies on a correction based on the extrapolated value of the in cell polynomial to the true geometry, thus not requiring the explicit evaluation of high order Taylor series. Moreover, this strategy could be easily implemented into any already existing finite element and finite volume code. Several validation tests are presented to prove the convergence properties up to order four for 2D and 3D simulations with curved boundaries, as well as an effective extension to flows with shocks.