2,254 research outputs found
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
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 Numerical Study of Methods for Moist Atmospheric Flows: Compressible Equations
We investigate two common numerical techniques for integrating reversible
moist processes in atmospheric flows in the context of solving the fully
compressible Euler equations. The first is a one-step, coupled technique based
on using appropriate invariant variables such that terms resulting from phase
change are eliminated in the governing equations. In the second approach, which
is a two-step scheme, separate transport equations for liquid water and vapor
water are used, and no conversion between water vapor and liquid water is
allowed in the first step, while in the second step a saturation adjustment
procedure is performed that correctly allocates the water into its two phases
based on the Clausius-Clapeyron formula. The numerical techniques we describe
are first validated by comparing to a well-established benchmark problem.
Particular attention is then paid to the effect of changing the time scale at
which the moist variables are adjusted to the saturation requirements in two
different variations of the two-step scheme. This study is motivated by the
fact that when acoustic modes are integrated separately in time (neglecting
phase change related phenomena), or when sound-proof equations are integrated,
the time scale for imposing saturation adjustment is typically much larger than
the numerical one related to the acoustics
Approximation of a compressible Navier-Stokes system by non-linear acoustical models
We analyse the existing derivation of the models of non-linear acoustics such
as the Kuznetsov equation, the NPE equation and the KZK equation. The technique
of introducing a corrector in the derivation ansatz allows to consider the
solutions of these equations as approximations of the solution of the initial
system (a com-pressible Navier-Stokes/Euler system). The validation of the
approximation ansatz is given for the KZK equation case
- …