Implicit finite element schemes for stationary compressible particle-laden gas flows

The derivation of macroscopic models for particle-laden gas flows is reviewed. Semi-implicit and Newton-like finite element methods are developed for the stationary two-fluid model governing compressible particle-laden gas flows. The Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable. To suppress numerical oscillations, the spatial discretization is performed by a high-resolution finite element scheme based on algebraic flux correction. A multidimensional limiter of TVD type is employed. An important goal is the efficient computation of stationary solutions in a wide range of Mach numbers. This is a challenging task due to oscillatory correction factors associated with TVD-type flux limiters and the additional strong nonlinearity caused by interfacial coupling terms. A semi-implicit scheme is derived by a time-lagged linearization of the nonlinear residual, and a Newton-like method is obtained in the limit of infinite CFL numbers. The original Jacobian is replaced by a low-order approximation. Special emphasis is laid on the numerical treatment of weakly imposed boundary conditions. It is shown that the proposed approach offers unconditional stability and faster convergence rates for increasing CFL numbers. The strongly coupled solver is compared to operator splitting techniques, which are shown to be less robust.     

Capabilities of a quasi-gasdynamic algorithm as applied to inviscid gas flow simulation

Ten well-known one-dimensional test problems reflecting the characteristic features of unsteady inviscid gas flows are successfully solved by a unified numerical algorithm based on the quasigasdynamic system of equations. In all the cases, the numerical solution converges to a self-similar one as the spatial grid is refined.     

The inviscid limit for an inflow problem of compressible viscous gas in presence of both shocks and boundary layers

In this paper, we study the inviscid limit problem for the Navier-Stokes equations of one-dimensional compressible viscous gas on half plane. We prove that if the solution of the inviscid Euler system on half plane is piecewise smooth with a single shock satisfying the entropy condition, then there exist solutions to Navier-Stokes equations which converge to the inviscid solution away from the shock discontinuity and the boundary at an optimal rate of ε¹ as the viscosity ε tends to zero.     

Finite element LES and VMS methods on tetrahedral meshes

Finite element methods for problems given in complex domains are often based on tetrahedral meshes. This paper demonstrates that the so-called rational Large Eddy Simulation model and a projection-based Variational Multiscale method can be extended in a straightforward way to tetrahedral meshes. Numerical studies are performed with an inf-sup stable second order pair of finite elements with discontinuous pressure approximation.     

On a regularization of the compressible Euler equations for an isothermal gas

We consider a Leray-type regularization of the compressible Euler equations for an isothermal gas. The regularized system depends on a small parameter α>0. Using Riemann invariants, we prove the existence of smooth solutions for the regularized system for every α>0. The regularization mechanism is a non-linear bending of characteristics that prevents their finite-time crossing. We prove that, in the α→0 limit, the regularized solutions converge strongly. However, based on our analysis and numerical simulations, the limit is not the unique entropy solution of the Euler equations. The numerical method used to support this claim is derived from the Riemann invariants for the regularized system. This method is guaranteed to preserve the monotonicity of characteristics.     

A large deviations principle for stochastic flows of viscous fluids

We study the well-posedness of a stochastic differential equation on the two dimensional torus ${\mathbb{T}}^2$, driven by an infinite dimensional Wiener process with drift in the Sobolev space $L^2(0,T;H^1\left({\mathbb{T}}^2\right))$. The solution corresponds to a stochastic Lagrangian flow in the sense of DiPerna Lions. By taking into account that the motion of a viscous incompressible fluid on the torus can be described through a suitable stochastic differential equation of the previous type, we study the inviscid limit. By establishing a large deviations principle, we show that, as the viscosity goes to zero, the Lagrangian stochastic Navier–Stokes flow approaches the Euler deterministic Lagrangian flow with an exponential rate function.     

Vortical and self-similar flows of 2D compressible Euler equations

This paper presents the vortical and self-similar solutions for 2D compressible Euler equations using the separation method. These solutions complement Makino’s solutions in radial symmetry without rotation. The rotational solutions provide new information that furthers our understanding of ocean vortices and reference examples for numerical methods. In addition, the corresponding blowup, time-periodic or global existence conditions are classified through an analysis of the new Emden equation. A conjecture regarding rotational solutions in 3D is also made.     

Note on the uniqueness of subsonic Euler flows in an axisymmetric nozzle

In this paper, we consider the uniqueness of globally subsonic compressible flows through an infinitely long axisymmetric nozzle. The flow is governed by the steady Euler equations and satisfies no-flow boundary conditions on the nozzle walls. We will show that for given mass flux and Bernoulli’s function in the upstream, the subsonic flow is unique in the class of all axisymmetric solutions, which possess the asymptotic behaviors at the far fields. This result extends the uniqueness of solutions in the previous paper Du and Duan (2011) [1].     

Finite element analysis of laminar and turbulent flows using LES and subgrid-scale models

Numerical simulations of laminar and turbulent flows in a lid driven cavity and over a backward-facing step are presented in this work. The main objectives of this research are to know more about the structure of turbulent flows, to identify their three-dimensional characteristic and to study physical effects due to heat transfer. The filtered Navier–Stokes equations are used to simulate large scales, however they are supplemented by subgrid-scale (SGS) models to simulate the energy transfer from large scales toward subgrid-scales, where this energy will be dissipated by molecular viscosity. Two SGS models are applied: the classical Smagorinsky’s model and the Dynamic model for large eddy simulation (LES). Both models are implemented in a three-dimensional finite element code using linear tetrahedral elements. Qualitative and quantitative aspects of two and three-dimensional flows in a lid-driven cavity and over a backward-facing step, using LES, are analyzed comparing numerical and experimental results obtained by other authors.     

Blow-up of viscous compressible reactive self-gravitating gas

In this paper,we prove some results concerning blow-up of viscous compressible reactive (selfgravitating) flows when the initial density is compactly supported and the other initial value satisfy proper conditions.It extends the work of Xin and Cho to the case of viscous compressible reactive self-gravitating flows equations.We control the lower bound of second moment by total energy and obtain the precise relationship between the size of the support of initial density and the existence time.     

An implicit finite volume method for the solution of 3D low Mach number viscous flows using a local preconditioning technique

This paper presents a cell-centered high order finite volume scheme for the solution of the three-dimensional (3D) Navier–Stokes equations with low Mach number. The system of non-linear equations is solved by means of a fully implicit pseudo-transient scheme. Each pseudo-time step is solved by a Newton-GMRes procedure. A local preconditioning technique is used to scale the speed of sound and to improve the system condition number for low Mach number and low cell Reynolds number. This preconditioning is applied to the AUSM+up flux vector splitting function. The method is tested on 2D and 3D low Mach number laminar flows.     

Existence and uniqueness of steady solutions to the magnetohydrodynamic equations of compressible flows

We establish the existence and uniqueness of a strong solution to the steady magnetohydrodynamic equations for the compressible barotropic fluids in a bounded smooth domain with a perfectly conducting boundary, under the assumption that the external force field is small.     

Finite volume element approximations of nonlocal reactive flows in porous media

In this article, we study finite volume element approximations for two‐dimensional parabolic integro‐differential equations, arising in the modeling of nonlocal reactive flows in porous media. These types of flows are also called NonFickian flows and exhibit mixing length growth. For simplicity, we consider only linear finite volume element methods, although higher‐order volume elements can be considered as well under this framework. It is proved that the finite volume element approximations derived are convergent with optimal order in H¹‐ and L²‐norm and are superconvergent in a discrete H¹‐norm. By examining the relationship between finite volume element and finite element approximations, we prove convergence in L^∞‐ and W^1,∞‐norms. These results are also new for finite volume element methods for elliptic and parabolic equations. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 285–311, 2000     

Direct numerical simulation of one type of compressible turbulence interacting with a shock wave

Direct numerical simulations of compressible turbulence interacting with an initially plane shock wave are presented. The underlying model is based on the numerical solution of the Euler equations combined with direct statistical simulation. Steady-state isentropic isotropic turbulence is considered. The amplification factors for fluctuations of the thermodynamic variables, velocity, vorticity, and kinetic energy of fluctuations are analyzed; and the correlation coefficients between flow variables are studied for Mach numbers ranging from 1.2 to 3.     

Analysis of preconditioning and multigrid for Euler flows with low-subsonic regions

For subsonic flows and upwind-discretized, linearized 1-D Euler equations, the smoothing behavior of multigrid-accelerated point Gauss-Seidel relaxation is analyzed. Error decay by convection across domain boundaries is also discussed. A fix to poor convergence rates at low Mach numbers is sought in replacing the point relaxation applied to unconditioned Euler equations, by locally implicit time-stepping applied to preconditioned Euler equations. The locally implicit iteration step is optimized for good damping of high-frequency errors. Numerical inaccuracy at low Mach numbers is also addressed.The work reported was performed in the framework of the BRITE-EURAM Aeronautics R&D Programme of the European Communities (Contract No. AER2-CT92-0040). The work was started during the second author's visit to CWI in 1993.     

Analysis of expanded mixed finite element methods for the generalized forchheimer flows of slightly compressible fluids

The nonlinear Forchheimer equations are used to describe the dynamics of fluid flows in porous media when Darcy's law is not applicable. In this article, we consider the generalized Forchheimer flows for slightly compressible fluids, and then study the expanded mixed finite element method applied to the initial boundary value problem for the resulting degenerate parabolic equation for pressure. The bounds for the solutions, time derivative, and gradient of solutions are established. Utilizing the monotonicity properties of Forchheimer equation and boundedness of solutions, a priori error estimates for solution are obtained in ‐norm, ‐norm as well as for its gradient in ‐norm for all . Optimal ‐error estimates are shown for solutions under some additional regularity assumptions. Numerical results using the lowest order Raviart–Thomas mixed element confirm the theoretical analysis regarding convergence rates. Published 2015. Numer Methods Partial Differential Eq 32: 60–85, 2016     

Error analysis of the L2 least‐squares finite element method for incompressible inviscid rotational flows

In this article we analyze the L² least‐squares finite element approximations to the incompressible inviscid rotational flow problem, which is recast into the velocity‐vorticity‐pressure formulation. The least‐squares functional is defined in terms of the sum of the squared L² norms of the residual equations over a suitable product function space. We first derive a coercivity type a priori estimate for the first‐order system problem that will play the crucial role in the error analysis. We then show that the method exhibits an optimal rate of convergence in the H¹ norm for velocity and pressure and a suboptimal rate of convergence in the L² norm for vorticity. A numerical example in two dimensions is presented, which confirms the theoretical error estimates. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004