A general Riemann solver for Euler equations

In this paper, we present a general Riemann solver which is applied successfully to compute the Euler equations in fluid dynamics with many complex equations of state (EOS). The solver is based on a splitting method introduced by the authors. We add a linear advection term to the Euler equations in the first step, to make the numerical flux between cells easy to compute. The added linear advection term is thrown off in the second step. It does not need an iterative technique and characteristic wave decomposition for computation. This new solver is designed to permit the construction of high‐order approximations to obtain high‐order Godunov‐type schemes. A number of numerical results show its robustness. Copyright © 2007 John Wiley & Sons, Ltd.     

Strategies to cure numerical shock instability in the HLLEM Riemann solver

The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver that is known to be plagued by various forms of numerical shock instability. In this paper, we clarify that the shock instability exhibited by this scheme is primarily triggered by the spurious activation of the antidiffusive terms present in the first and third Riemann flux components on the transverse interfaces adjoining the shock front due to numerical perturbations. These erroneously activated terms are shown to counteract the favorable damping mechanism provided by its inherent HLL-type diffusive terms, causing an unphysical variation of the conserved quantity ρu both along and across the numerical shock. To prevent this, two distinct strategies are proposed termed as S elective W ave M odification and A nti D iffusion C ontrol. The former focuses on enhancing the quantity of the favorable HLL-type dissipation available on these critical flux components by carefully increasing the magnitudes of certain nonlinear wave speed estimates, while the latter focuses on directly controlling the magnitude of these critical antidiffusive terms. A linear perturbation analysis is performed to gauge the effectiveness of these cures and to estimate a von Neumann–type stability bounds on the CFL number associated with their use. Results from a variety of classic shock instability test cases show that the proposed strategies are able to provide excellent shock stable solutions even on grids that are highly elongated across the shock front without compromising the accuracy on inviscid contact or shear dominated viscous flows.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

An Approximate Roe-Type Riemann Solver for a Class of Realizable Second Order Closures

A realizable, objective second-moment turbulence closure, allowing for an entropy characterisation, is analyzed with respect to its convective subset. The distinct characteristic wave system of these equations in non-conservation form is exposed. An approximate solution to Ihe associated one-dimensional Riemann problem is constructed making use of approximate jump conditions obtained by assuming a linear path across shock waves. A numerical integration method based on a new approximate Riemann solver (flux-difference-splitting) is proposed for use in conjunction with either unstructured or structured grids. Test calculations of quasi one-dimensional flow cases demonstrate the feasibility of the current technique even where Euler-based approaches fail.     

Asymptotic solution of the thin viscous shock layer equations at low Reynolds numbers for a cold surface

Hypersonic three-dimensional viscous rarefied gas flow past blunt bodies in the neighborhood of the stagnation line is considered. The question of the applicability of the gasdynamic thin viscous shock layer model [1] is investigated for the transition flow regime from continuum to free-molecular flow. It is shown that for a power-law temperature dependence of the viscosity coefficient T the quantity (Re)^1/(1+), where = ( – 1)/2 and is the specific heat ratio, is an important determining parameter of the hypersonic flow at low Reynolds numbers. In the case of a cold surface approximate asymptotic solutions of the thin viscous shock layer equations are obtained for noslip conditions on the surface and generalized Rankine-Hugoniot relations on the shock wave at low Reynolds numbers. These solutions give simple analytic expressions for the thermal conductivity and friction coefficients as functions of the determining flow parameters. As the Reynolds number tends to zero, the values of the thermal conductivity and friction coefficients determined by this solution tend to their values in free-molecular flow for an accommodation coefficient equal to unity. This tending of the thermal conductivity and friction coefficients to the free-molecular limit takes place for both two-and three-dimensional flows. The asymptotic solutions are compared with numerical calculations and experimental data.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 159–170. Original Russian Text Copyright © 2004 by Brykina.     

Multiphase Thermohaline Convection in the Earth’s Crust: I. A New Finite Element – Finite Volume Solution Technique Combined With a New Equation of State for NaCl–H2O

We present a new finite element – finite volume (FEFV) method combined with a realistic equation of state for NaCl–H₂O to model fluid convection driven by temperature and salinity gradients. This method can deal with the nonlinear variations in fluid properties, separation of a saline fluid into a high-density, high-salinity brine phase and low-density, low-salinity vapor phase well above the critical point of pure H₂O, and geometrically complex geological structures. Similar to the well-known implicit pressure explicit saturation formulation, this approach decouples the governing equations. We formulate a fluid pressure equation that is solved using an implicit finite element method. We derive the fluid velocities from the updated pressure field and employ them in a higher-order, mass conserving finite volume formulation to solve hyperbolic parts of the conservation laws. The parabolic parts are solved by finite element methods. This FEFV method provides for geometric flexibility and numerical efficiency. The equation of state for NaCl–H₂O is valid from 0 to 750°C, 0 to 4000 bar, and 0–100 wt.% NaCl. This allows the simulation of thermohaline convection in high-temperature and high-pressure environments, such as continental or oceanic hydrothermal systems where phase separation is common.     

Multiphase Thermohaline Convection in the Earth’s Crust: II. Benchmarking and Application of a Finite Element – Finite Volume Solution Technique with a NaCl–H2O Equation of State

We present the benchmarking of a new finite element – finite volume (FEFV) solution technique capable of modeling transient multiphase thermohaline convection for geological realistic p-T-X conditions. The algorithm embeds a new and accurate equation of state for the NaCl–H₂O system. Benchmarks are carried out to compare the numerical results for the various component-processes of multiphase thermohaline convection. They include simulations of (i) convection driven by temperature and/or concentration gradients in a single-phase fluid (i.e., the Elder problem, thermal convection at different Rayleigh numbers, and a free thermohaline convection example), (ii) multiphase flow (i.e., the Buckley–Leverett problem), and (iii) energy transport in a pure H₂O fluid at liquid, vapor, supercritical, and two-phase conditions (i.e., comparison to the U.S. Geological Survey Code HYDROTHERM). The results produced with the new FEFV technique are in good agreement with the reference solutions. We further present the application of the FEFV technique to the simulation of thermohaline convection of a 400°C hot and 10 wt.% saline fluid rising from 4 km depth. During the buoyant rise, the fluid boils and separates into a high-density, high-salinity liquid phase and a low-density, low-salinity vapor phase.