1,728 research outputs found
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
All speed scheme for the low mach number limit of the Isentropic Euler equation
An all speed scheme for the Isentropic Euler equation is presented in this
paper. When the Mach number tends to zero, the compressible Euler equation
converges to its incompressible counterpart, in which the density becomes a
constant. Increasing approximation errors and severe stability constraints are
the main difficulty in the low Mach regime. The key idea of our all speed
scheme is the special semi-implicit time discretization, in which the low Mach
number stiff term is divided into two parts, one being treated explicitly and
the other one implicitly. Moreover, the flux of the density equation is also
treated implicitly and an elliptic type equation is derived to obtain the
density. In this way, the correct limit can be captured without requesting the
mesh size and time step to be smaller than the Mach number. Compared with
previous semi-implicit methods, nonphysical oscillations can be suppressed. We
develop this semi-implicit time discretization in the framework of a first
order local Lax-Friedrich (LLF) scheme and numerical tests are displayed to
demonstrate its performances
A hierarchy for modeling high speed propulsion systems
General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery, such as inlets, ramjets, and scramjets. The discussion is separated into four areas: (1) computational fluid dynamics models for the entire nonlinear system or high order nonlinear models; (2) high order linearized models derived from fundamental physics; (3) low order linear models obtained from the other high order models; and (4) low order nonlinear models (order here refers to the number of dynamic states). Included in the discussion are any special considerations based on the relevant control system designs. The methods discussed are for the quasi-one-dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, including moving normal shocks, hammershocks, simple subsonic combustion via heat addition, temperature dependent gases, detonations, and thermal choking. The report also contains a comprehensive list of papers and theses generated by this grant
A Space-time Smooth Artificial Viscosity Method For Nonlinear Conservation Laws
We introduce a new methodology for adding localized, space-time smooth,
artificial viscosity to nonlinear systems of conservation laws which propagate
shock waves, rarefactions, and contact discontinuities, which we call the
-method. We shall focus our attention on the compressible Euler equations in
one space dimension. The novel feature of our approach involves the coupling of
a linear scalar reaction-diffusion equation to our system of conservation laws,
whose solution is the coefficient to an additional (and artificial)
term added to the flux, which determines the location, localization, and
strength of the artificial viscosity. Near shock discontinuities, is
large and localized, and transitions smoothly in space-time to zero away from
discontinuities. Our approach is a provably convergent, spacetime-regularized
variant of the original idea of Richtmeyer and Von Neumann, and is provided at
the level of the PDE, thus allowing a host of numerical discretization schemes
to be employed. We demonstrate the effectiveness of the -method with three
different numerical implementations and apply these to a collection of
classical problems: the Sod shock-tube, the Osher-Shu shock-tube, the
Woodward-Colella blast wave and the Leblanc shock-tube. First, we use a
classical continuous finite-element implementation using second-order
discretization in both space and time, FEM-C. Second, we use a simplified WENO
scheme within our -method framework, WENO-C. Third, we use WENO with the
Lax-Friedrichs flux together with the -equation, and call this WENO-LF-C.
All three schemes yield higher-order discretization strategies, which provide
sharp shock resolution with minimal overshoot and noise, and compare well with
higher-order WENO schemes that employ approximate Riemann solvers,
outperforming them for the difficult Leblanc shock tube experiment.Comment: 34 pages, 27 figure
Flux Splitting for stiff equations: A notion on stability
For low Mach number flows, there is a strong recent interest in the
development and analysis of IMEX (implicit/explicit) schemes, which rely on a
splitting of the convective flux into stiff and nonstiff parts. A key
ingredient of the analysis is the so-called Asymptotic Preserving (AP)
property, which guarantees uniform consistency and stability as the Mach number
goes to zero. While many authors have focussed on asymptotic consistency, we
study asymptotic stability in this paper: does an IMEX scheme allow for a CFL
number which is independent of the Mach number? We derive a stability criterion
for a general linear hyperbolic system. In the decisive eigenvalue analysis,
the advective term, the upwind diffusion and a quadratic term stemming from the
truncation in time all interact in a subtle way. As an application, we show
that a new class of splittings based on characteristic decomposition, for which
the commutator vanishes, avoids the deterioration of the time step which has
sometimes been observed in the literature
- …