1,229 research outputs found
Augmented resolution of linear hyperbolic systems under nonconservative form
Hyperbolic systems under nonconservative form arise in numerous applications
modeling physical processes, for example from the relaxation of more general
equations (e.g. with dissipative terms). This paper reviews an existing class
of augmented Roe schemes and discusses their application to linear
nonconservative hyperbolic systems with source terms. We extend existing
augmented methods by redefining them within a common framework which uses a
geometric reinterpretation of source terms. This results in intrinsically
well-balanced numerical discretizations. We discuss two equivalent
formulations: (1) a nonconservative approach and (2) a conservative
reformulation of the problem. The equilibrium properties of the schemes are
examined and the conditions for the preservation of the well-balanced property
are provided. Transient and steady state test cases for linear acoustics and
hyperbolic heat equations are presented. A complete set of benchmark problems
with analytical solution, including transient and steady situations with
discontinuities in the medium properties, are presented and used to assess the
equilibrium properties of the schemes. It is shown that the proposed schemes
satisfy the expected equilibrium and convergence properties
A high-order nonconservative approach for hyperbolic equations in fluid dynamics
It is well known, thanks to Lax-Wendroff theorem, that the local conservation
of a numerical scheme for a conservative hyperbolic system is a simple and
systematic way to guarantee that, if stable, a scheme will provide a sequence
of solutions that will converge to a weak solution of the continuous problem.
In [1], it is shown that a nonconservative scheme will not provide a good
solution. The question of using, nevertheless, a nonconservative formulation of
the system and getting the correct solution has been a long-standing debate. In
this paper, we show how get a relevant weak solution from a pressure-based
formulation of the Euler equations of fluid mechanics. This is useful when
dealing with nonlinear equations of state because it is easier to compute the
internal energy from the pressure than the opposite. This makes it possible to
get oscillation free solutions, contrarily to classical conservative methods.
An extension to multiphase flows is also discussed, as well as a
multidimensional extension
Well-balanced -adaptive and moving mesh space-time discontinuous Galerkin method for the shallow water equations
In this article we introduce a well-balanced discontinuous Galerkin method for the shallow water equations on moving meshes. Particular emphasis will be given on -adaptation in which mesh points of an initially uniform mesh move to concentrate in regions where interesting behaviour of the solution is observed. Obtaining well-balanced numerical schemes for the shallow water equations on fixed meshes is nontrivial and has been a topic of much research. In [S. Rhebergen, O. Bokhove, J.J.W. van der Vegt, Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations, J. Comput. Phys. 227 (2008) 1887–1922] we introduced a well-balanced discontinuous Galerkin method using the theory of weak solutions for nonconservative products introduced in [G. Dal Maso, P.G. LeFloch, F. Murat, Definition and weak stability of nonconservative products, J. Math. Pures Appl. 74 (1995) 483–548]. In this article we continue this approach and prove well-balancedness of a discontinuous Galerkin method for the shallow water equations on moving meshes. Numerical simulations are then performed to verify the -adaptive method in combination with the space-time discontinuous Galerkin method against analytical solutions and showing its robustness on more complex problems
A water wave model with horizontal circulation and accurate dispersion
We describe a new water wave model which is variational, and combines a depth-averaged vertical (component of) vorticity with depth-dependent potential flow. The model facilitates the further restriction of the vertical profile of the velocity potential to n-th order polynomials or a finite element profile with a small number of elements (say), leading to a framework for efficient modelling of the interaction of steepening and breaking waves near the shore with a large-scale horizontal flow. The equations are derived from a constrained variational formulation which leads to conservation laws for energy, mass, momentum and vertical vorticity (or circulation). We show that the potential flow water wave equations and the shallow-water equations are recovered in the relevant limits, and provide approximate shock relations for the model which can be used in numerical schemes to model breaking waves
Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations
The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme
Entropy stable DGSEM for nonlinear hyperbolic systems in nonconservative form with application to two-phase flows
In this work, we consider the discretization of nonlinear hyperbolic systems
in nonconservative form with the high-order discontinuous Galerkin spectral
element method (DGSEM) based on collocation of quadrature and interpolation
points (Kopriva and Gassner, J. Sci. Comput., 44 (2010), pp.136--155; Carpenter
et al., SIAM J. Sci. Comput., 36 (2014), pp.~B835-B867). We present a general
framework for the design of such schemes that satisfy a semi-discrete entropy
inequality for a given convex entropy function at any approximation order. The
framework is closely related to the one introduced for conservation laws by
Chen and Shu (J. Comput. Phys., 345 (2017), pp.~427--461) and relies on the
modification of the integral over discretization elements where we replace the
physical fluxes by entropy conservative numerical fluxes from Castro et al.
(SIAM J. Numer. Anal., 51 (2013), pp.~1371--1391), while entropy stable
numerical fluxes are used at element interfaces. Time discretization is
performed with strong-stability preserving Runge-Kutta schemes. We use this
framework for the discretization of two systems in one space-dimension: a
system with a nonconservative product associated to a
linearly-degenerate field for which the DGSEM fails to capture the physically
relevant solution, and the isentropic Baer-Nunziato model. For the latter, we
derive conditions on the numerical parameters of the discrete scheme to further
keep positivity of the partial densities and a maximum principle on the void
fractions. Numerical experiments support the conclusions of the present
analysis and highlight stability and robustness of the present schemes
- …