scholarly journals All Speed Scheme for the Low Mach Number Limit of the Isentropic Euler Equations

2011 ◽  
Vol 10 (1) ◽  
pp. 1-31 ◽  
Author(s):  
Pierre Degond ◽  
Min Tang
Keyword(s):  
Mach Number ◽  
Euler Equations ◽  
Computational Cost ◽  
Time Discretization ◽  
Implicit Time ◽  
Elliptic Type ◽  
Implicit Methods ◽  
Time Step ◽  
Low Mach Number ◽  
Isentropic Euler Equations

AbstractAn all speed scheme for the Isentropic Euler equations is presented in this paper. When the Mach number tends to zero, the compressible Euler equations converge to their 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, firstly, nonphysical oscillations can be suppressed by choosing proper parameter, besides, only a linear elliptic equation needs to be solved implicitly which reduces much computational cost. We develop this semi-implicit time discretization in the framework of a first order Local Lax-Friedrichs (or Rusanov) scheme and numerical tests are displayed to demonstrate its performances.

GLOBAL EXISTENCE FOR PHASE TRANSITION PROBLEMS VIA A VARIATIONAL SCHEME

2004 ◽  
Vol 01 (04) ◽  
pp. 747-768
Author(s):  
CHRISTIAN ROHDE ◽  
MAI DUC THANH
Keyword(s):  
Phase Transition ◽  
Surface Tension ◽  
Initial Data ◽  
Time Discretization ◽  
Approximate Solutions ◽  
Implicit Time ◽  
Time Step ◽  
Dynamical Phase Transition ◽  
Transition Dynamics ◽  
Variational Scheme

We construct approximate solutions of the initial value problem for dynamical phase transition problems via a variational scheme in one space dimension. First, we deal with a local model of phase transition dynamics which contains second and third order spatial derivatives modeling the effects of viscosity and surface tension. Assuming that the initial data are periodic, we prove the convergence of approximate solutions to a weak solution which satisfies the natural dissipation inequality. We note that this result still holds for non-periodic initial data. Second, we consider a model of phase transition dynamics with only Lipschitz continuous stress–strain function which contains a non-local convolution term to take account of surface tension. We also establish the existence of weak solutions. In both cases the proof relies on implicit time discretization and the analysis of a minimization problem at each time step.

scholarly journals Exact solution and the multidimensional Godunov scheme for the acoustic equations

2022 ◽  
Author(s):  
Wasilij Barsukow ◽  
Christian Klingenberg
Keyword(s):  
Exact Solution ◽  
Mach Number ◽  
Euler Equations ◽  
Two Dimensions ◽  
Finite Volume Scheme ◽  
Godunov Scheme ◽  
Low Mach Number ◽  
Low Mach Number Limit ◽  
Acoustic Equations ◽  
Spatial Dimensions

The acoustic equations derived as a linearization of the Euler equations are a valuable system for studies of multi-dimensional solutions. Additionally they possess a low Mach number limit analogous to that of the Euler equations. Aiming at understanding the behaviour of the multi-dimensional Godunov scheme in this limit, first the exact solution of the corresponding Cauchy problem in three spatial dimensions is derived. The appearance of logarithmic singularities in the exact solution of the 4-quadrant Riemann Problem in two dimensions is discussed. The solution formulae are then used to obtain the multidimensional Godunov finite volume scheme in two dimensions. It is shown to be superior to the dimensionally split upwind/Roe scheme concerning its domain of stability and ability to resolve multi-dimensional Riemann problems. It is shown experimentally and theoretically that despite taking into account multi-dimensional information it is, however, not able to resolve the low Mach number limit.

A FULLY IMPLICIT SOLVER OF 3-D EULER EQUATIONS ON MULTIBLOCK CURVILINEAR GRIDS

2005 ◽  
Vol 19 (28n29) ◽  
pp. 1483-1486 ◽  
Author(s):  
HAI-QING SI ◽  
TONG-GUANG WANG ◽  
XIAO-YUN LUO
Keyword(s):  
Euler Equations ◽  
Three Dimensional ◽  
Implicit Time ◽  
Tvd Scheme ◽  
Time Step ◽  
Vector Operation ◽  
Fully Implicit ◽  
Curvilinear Grids ◽  
The Matrix ◽  
Implicit Solver

A fully implicit unfactored algorithm for three-dimensional Euler equations is developed and tested on multi-block curvilinear meshes. The convective terms are discretized using an upwind TVD scheme. The large sparse linear system generated at each implicit time step is solved by GMRES* method combined with the block incomplete lower-upper preconditioner. In order to reduce the memory requirements and the matrix-vector operation counts, an approximate method is used to derive the Jacobian matrix, which only costs half of the computational efforts of the exact Jacobian calculation. The comparison between the numerical results and the experimental data shows good agreement, which demonstrates that the implicit algorithm presented is effective and efficient.

A Fast Semi-Implicit Level Set Method for Curvature Dependent Flows with an Application to Limit Cycles Extraction in Dynamical Systems

2015 ◽  
Vol 18 (1) ◽  
pp. 203-229 ◽  
Author(s):  
Guoqiao You ◽  
Shingyu Leung
Keyword(s):  
Dynamical Systems ◽  
Level Set Method ◽  
Level Set ◽  
Limit Cycles ◽  
Time Discretization ◽  
Numerical Approach ◽  
Initial Guess ◽  
Implicit Time ◽  
Time Step ◽  
Image Regularization

AbstractWe propose a new semi-implicit level set approach to a class of curvature dependent flows. The method generalizes a recent algorithm proposed for the motion by mean curvature where the interface is updated by solving the Rudin-Osher-Fatemi (ROF) model for image regularization. Our proposal is general enough so that one can easily extend and apply the method to other curvature dependent motions. Since the derivation is based on a semi-implicit time discretization, this suggests that the numerical scheme is stable even using a time-step significantly larger than that of the corresponding explicit method. As an interesting application of the numerical approach, we propose a new variational approach for extracting limit cycles in dynamical systems. The resulting algorithm can automatically detect multiple limit cycles staying inside the initial guess with no condition imposed on the number nor the location of the limit cycles. Further, we also propose in this work an Eulerian approach based on the level set method to test if the limit cycles are stable or unstable.

scholarly journals A Hybrid Adaptive Low-Mach-Number/Compressible Method for the Euler Equations

2017 ◽  
Author(s):  
Emmanuel Motheau ◽  
Max Duarte ◽  
Ann S. Almgren ◽  
John B. Bell
Keyword(s):  
Mach Number ◽  
Euler Equations ◽  
Low Mach Number

scholarly journals Contrast-Independent, Partially-Explicit Time Discretizations for Nonlinear Multiscale Problems

Mathematics ◽  
2021 ◽  
Vol 9 (23) ◽  
pp. 3000
Author(s):  
Eric T. Chung ◽  
Yalchin Efendiev ◽  
Wing Tat Leung ◽  
Wenyuan Li
Keyword(s):  
Linear Equations ◽  
Nonlinear Problems ◽  
Time Discretization ◽  
Numerical Results ◽  
Explicit Methods ◽  
Implicit Methods ◽  
Time Step ◽  
Multiscale Problems ◽  
Time Stepping ◽  
Space Decomposition

This work continues a line of work on developing partially explicit methods for multiscale problems. In our previous works, we considered linear multiscale problems where the spatial heterogeneities are at the subgrid level and are not resolved. In these works, we have introduced contrast-independent, partially explicit time discretizations for linear equations. The contrast-independent, partially explicit time discretization divides the spatial space into two components: contrast dependent (fast) and contrast independent (slow) spaces defined via multiscale space decomposition. Following this decomposition, temporal splitting was proposed, which treats fast components implicitly and slow components explicitly. The space decomposition and temporal splitting are chosen such that they guarantees stability, and we formulated a condition for the time stepping. This condition was formulated as a condition on slow spaces. In this paper, we extend this approach to nonlinear problems. We propose a splitting approach and derive a condition that guarantees stability. This condition requires some type of contrast-independent spaces for slow components of the solution. We present numerical results and show that the proposed methods provide results similar to implicit methods with a time step that is independent of the contrast.

Sign in / Sign up

Export Citation Format

Share Document