Numerical solution of inviscid and viscous flows using modern schemes and quadrilateral or triangular mesh
On a hybrid finite-volume-particle method
We present a hybrid finite-volume-particle numerical method for computing the transport of a passive pollutant by a flow. The flow is modeled by the one- and two-dimensional Saint-Venant system of shallow water equations and the pollutant propagation is described by a transport equation. This paper is an extension of our previous work [Chertock, Kurganov and Petrova, J. Sci. Comput. (to appear)], where the one-dimensional finite-volume-particle method has been proposed. The core idea behind the...
On a hybrid finite-volume-particle method
We present a hybrid finite-volume-particle numerical method for computing the transport of a passive pollutant by a flow. The flow is modeled by the one- and two-dimensional Saint-Venant system of shallow water equations and the pollutant propagation is described by a transport equation. This paper is an extension of our previous work [Chertock, Kurganov and Petrova, J. Sci. Comput. (to appear)], where the one-dimensional finite-volume-particle method has been proposed. The core idea behind the...
On a stabilized colocated Finite Volume scheme for the Stokes problem
We present and analyse in this paper a novel colocated Finite Volume scheme for the solution of the Stokes problem. It has been developed following two main ideas. On one hand, the discretization of the pressure gradient term is built as the discrete transposed of the velocity divergence term, the latter being evaluated using a natural finite volume approximation; this leads to a non-standard interpolation formula for the expression of the pressure on the edges of the control volumes. On the other...
On discontinuous Galerkin methods for nonlinear convection-diffusion problems and compressible flow
On discontinuous Galerkin methods for nonlinear convection-diffusion problems and compressible flow
The paper is concerned with the discontinuous Galerkin finite element method for the numerical solution of nonlinear conservation laws and nonlinear convection-diffusion problems with emphasis on applications to the simulation of compressible flows. We discuss two versions of this method: (a) Finite volume discontinuous Galerkin method, which is a generalization of the combined finite volume—finite element method. Its advantage is the use of only one mesh (in contrast to the combined finite volume—finite...
On numerical solution of the incompressible Navier-Stokes equations with static or total pressure specified on boundaries.
On some implicit and semi-implicit staggered schemes for the shallow water and Euler equations
In this paper, we propose implicit and semi-implicit in time finite volume schemes for the barotropic Euler equations (hence, as a particular case, for the shallow water equations) and for the full Euler equations, based on staggered discretizations. For structured meshes, we use the MAC finite volume scheme, and, for general mixed quadrangular/hexahedral and simplicial meshes, we use the discrete unknowns of the Rannacher−Turek or Crouzeix−Raviart finite elements. We first show that a solution...
On the Charney Conjecture of Data Assimilation Employing Temperature Measurements Alone: The Paradigm of 3D Planetary Geostrophic Model
Analyzing the validity and success of a data assimilation algorithmwhen some state variable observations are not available is an important problem in meteorology and engineering. We present an improved data assimilation algorithm for recovering the exact full reference solution (i.e. the velocity and temperature) of the 3D Planetary Geostrophic model, at an exponential rate in time, by employing coarse spatial mesh observations of the temperature alone. This provides, in the case of this paradigm,...
On the connection between some Riemann-solver free approaches to the approximation of multi-dimensional systems of hyperbolic conservation laws
In this paper, we present some interesting connections between a number of Riemann-solver free approaches to the numerical solution of multi-dimensional systems of conservation laws. As a main part, we present a new and elementary derivation of Fey’s Method of Transport (MoT) (respectively the second author’s ICE version of the scheme) and the state decompositions which form the basis of it. The only tools that we use are quadrature rules applied to the moment integral used in the gas kinetic derivation...
On the connection between some Riemann-solver free approaches to the approximation of multi-dimensional systems of hyperbolic conservation laws
In this paper, we present some interesting connections between a number of Riemann-solver free approaches to the numerical solution of multi-dimensional systems of conservation laws. As a main part, we present a new and elementary derivation of Fey's Method of Transport (MoT) (respectively the second author's ICE version of the scheme) and the state decompositions which form the basis of it. The only tools that we use are quadrature rules applied to the moment integral used in the...
Reduced basis method for finite volume approximations of parametrized linear evolution equations
The model order reduction methodology of reduced basis (RB) techniques offers efficient treatment of parametrized partial differential equations (P2DEs) by providing both approximate solution procedures and efficient error estimates. RB-methods have so far mainly been applied to finite element schemes for elliptic and parabolic problems. In the current study we extend the methodology to general linear evolution schemes such as finite volume schemes for parabolic and hyperbolic evolution equations....
Relaxation and numerical approximation of a two-fluid two-pressure diphasic model
This paper is concerned with the numerical approximation of the solutions of a two-fluid two-pressure model used in the modelling of two-phase flows. We present a relaxation strategy for easily dealing with both the nonlinearities associated with the pressure laws and the nonconservative terms that are inherently present in the set of convective equations and that couple the two phases. In particular, the proposed approximate Riemann solver is given by explicit formulas, preserves the natural...
Relaxation models of phase transition flows
In this work, we propose a general framework for the construction of pressure law for phase transition. These equations of state are particularly suitable for a use in a relaxation finite volume scheme. The approach is based on a constrained convex optimization problem on the mixture entropy. It is valid for both miscible and immiscible mixtures. We also propose a rough pressure law for modelling a super-critical fluid.
Relaxation schemes for the multicomponent Euler system
We show that it is possible to construct a class of entropic schemes for the multicomponent Euler system describing a gas or fluid homogeneous mixture at thermodynamic equilibrium by applying a relaxation technique. A first order Chapman–Enskog expansion shows that the relaxed system formally converges when the relaxation frequencies go to the infinity toward a multicomponent Navier–Stokes system with the classical Fick and Newton laws, with a thermal diffusion which can be assimilated to a Soret...
Relaxation schemes for the multicomponent Euler system
We show that it is possible to construct a class of entropic schemes for the multicomponent Euler system describing a gas or fluid homogeneous mixture at thermodynamic equilibrium by applying a relaxation technique. A first order Chapman–Enskog expansion shows that the relaxed system formally converges when the relaxation frequencies go to the infinity toward a multicomponent Navier–Stokes system with the classical Fick and Newton laws, with a thermal diffusion which can be assimilated to a Soret...
Residual and hierarchical a posteriori error estimates for nonconforming mixed finite element methods
We analyze residual and hierarchical a posteriori error estimates for nonconforming finite element approximations of elliptic problems with variable coefficients. We consider a finite volume box scheme equivalent to a nonconforming mixed finite element method in a Petrov–Galerkin setting. We prove that all the estimators yield global upper and local lower bounds for the discretization error. Finally, we present results illustrating the efficiency of the estimators, for instance, in the simulation...
Residual and hierarchical a posteriori error estimates for nonconforming mixed finite element methods
We analyze residual and hierarchical a posteriori error estimates for nonconforming finite element approximations of elliptic problems with variable coefficients. We consider a finite volume box scheme equivalent to a nonconforming mixed finite element method in a Petrov–Galerkin setting. We prove that all the estimators yield global upper and local lower bounds for the discretization error. Finally, we present results illustrating the efficiency of the estimators, for instance, in the simulation...
Second-order MUSCL schemes based on Dual Mesh Gradient Reconstruction (DMGR)
We discuss new MUSCL reconstructions to approximate the solutions of hyperbolic systems of conservations laws on 2D unstructured meshes. To address such an issue, we write two MUSCL schemes on two overlapping meshes. A gradient reconstruction procedure is next defined by involving both approximations coming from each MUSCL scheme. This process increases the number of numerical unknowns, but it allows to reconstruct very accurate gradients. Moreover a particular attention is paid on the limitation...
Simulation of the Three-Dimensional Flow of Blood Using a Shear-Thinning Viscoelastic Fluid Model
This paper is concerned with the numerical simulation of a thermodynamically compatible viscoelastic shear-thinning fluid model, particularly well suited to describe the rheological response of blood, under physiological conditions. Numerical simulations are performed in two idealized three-dimensional geometries, a stenosis and a curved vessel, to investigate the combined effects of flow inertia, viscosity and viscoelasticity in these geometries....