G. Farag, S. Zhao, T. Coratger, Pierre Boivin, G. Chiavassa, et al.. A pressure-based regularized lattice-Boltzmann method for the simulation of compressible flows. Physics of Fluids, American Institute of Physics, 2020, 32 (6), pp.066106. ⟨10.1063/5.0011839⟩. ⟨hal-02885427⟩ Plus de détails...
A new pressure-based Lattice-Boltzmann method (HRR-p) is proposed for the simulation of flows for Mach numbers ranging from 0 to 1.5. Compatible with nearest neighbor lattices (e.g. D3Q19), the model consists of a predictor step comparable to classical athermal Lattice-Boltzmann methods, appended with a fully local and explicit correction step for the pressure. Energy conservation-for which the Hermi-tian quadrature is not accurate enough on such lattice-is solved via a classical finite volume MUSCL-Hancock scheme based on the entropy equation. The Euler part of the model is then validated for the transport of three canonical modes (vortex, en-tropy, and acoustic propagation), while its diffusive/viscous properties are assessed via thermal Couette flow simulations. All results match the analytical solutions, with very limited dissipation. Lastly, the robustness of the method is tested in a one dimensional shock tube and a two-dimensional shock-vortex interaction.
G. Farag, S. Zhao, T. Coratger, Pierre Boivin, G. Chiavassa, et al.. A pressure-based regularized lattice-Boltzmann method for the simulation of compressible flows. Physics of Fluids, American Institute of Physics, 2020, 32 (6), pp.066106. ⟨10.1063/5.0011839⟩. ⟨hal-02885427⟩
Xi Deng, Pierre Boivin. Diffuse interface modelling of reactive multi-phase flows applied to a sub-critical cryogenic jet. Applied Mathematical Modelling, Elsevier, 2020, ⟨10.1016/j.apm.2020.04.011⟩. ⟨hal-02561937⟩ Plus de détails...
In order to simulate cryogenic H 2 − O 2 jets under subcritical condition, a numerical model is constructed to solve compressible reactive multi-component flows which involve complex multi-physics processes such as moving material interfaces, shock waves, phase transition and combustion. The liquid and reactive gaseous mixture are described by a homogeneous mixture model with diffusion transport for heat, momentum and species. A hybrid thermodynamic closure strategy is proposed to construct an equation of state (EOS) for the mixture. The phase transition process is modeled by a recent fast relaxation method which gradually reaches the thermo-chemical equilibrium without iterative process. A simplified transport model is also implemented to ensure the accurate behavior in the limit of pure fluids and maintain computational efficiency. Last, a 12-step chemistry model is included to account for hydrogen combustion. Then the developed numerical model is solved with the finite volume method where a low dissipation AUSM (advection upstream splitting method) Riemann solver is extended for multi-component flows. A homogeneous reconstruction strategy compatible with the homogeneous mixture model is adopted to prevent numerical oscillations across material interfaces. Having included these elements, the model is validated on a number of canonical configurations, first for multi-phase flows, and second for reactive flows. These tests allow recovery of the expected behavior in both the multiphase and reactive limits, and the model capability is further demonstrated on a 2D burning cryogenic H 2 − O 2 jet, in a configuration reminiscent of rocket engine ignition.
Xi Deng, Pierre Boivin. Diffuse interface modelling of reactive multi-phase flows applied to a sub-critical cryogenic jet. Applied Mathematical Modelling, Elsevier, 2020, ⟨10.1016/j.apm.2020.04.011⟩. ⟨hal-02561937⟩
M. Tayyab, S. Zhao, Y. Feng, Pierre Boivin. Hybrid regularized Lattice-Boltzmann modelling of premixed and non-premixed combustion processes. Combustion and Flame, Elsevier, 2020, 211, pp.173-184. ⟨10.1016/j.combustflame.2019.09.029⟩. ⟨hal-02346556⟩ Plus de détails...
A Lattice-Boltzmann model for low-Mach reactive flows is presented, built upon our recently published model (Comb & Flame, 196, 2018). The approach is hybrid and couples a Lattice-Boltzmann solver for the resolution of mass and momentum conservation and a finite difference solver for the energy and species conservation. Having lifted the constant thermodynamic and transport properties assumptions, the model presented now fully accounts for the classical reactive flow thermodynamic closure: each component is assigned NASA coefficients for calculating its thermodynamic properties. A temperature-dependent viscosity is considered, from which are deduced thermo-diffusive properties via specification of Prandtl and component-specific Schmidt numbers. Another major improvement from our previous contribution is the derivation of an advanced collision kernel compatible of multi-component reactive flows stable in high shear flows. Validation is carried out first on premixed configurations, through simulation of the planar freely propagating flame, the growth of the associated Darrieus-Landau instability and three regimes of flame-vortex interaction. A double shear layer test case including a flow-stabilized diffusion flame is then presented and results are compared with DNS simulations, showing excellent agreement.
M. Tayyab, S. Zhao, Y. Feng, Pierre Boivin. Hybrid regularized Lattice-Boltzmann modelling of premixed and non-premixed combustion processes. Combustion and Flame, Elsevier, 2020, 211, pp.173-184. ⟨10.1016/j.combustflame.2019.09.029⟩. ⟨hal-02346556⟩
Yongliang Feng, Pierre Boivin, Jérome Jacob, Pierre Sagaut. Hybrid recursive regularized lattice Boltzmann simulation of humid air with application to meteorological flows. Physical Review E , American Physical Society (APS), 2019. ⟨hal-02265484⟩ Plus de détails...
An extended version of the hybrid recursive regularized Lattice-Boltzmann model which incorporates external force is developed to simulate humid air flows with phase change mechanisms under the Boussinesq approximation. Mass and momentum conservation equations are solved by a regu-larized lattice Boltzmann approach well suited for high Reynolds number flows, whereas the energy and humidity related equations are solved by a finite volume approach. Two options are investigated to account for cloud formation in atmospheric flow simulations. The first option considers a single conservation equation for total water and an appropriate invariant variable of temperature. In the other approach, liquid and vapor are considered via two separated equations, and phase transition is accounted for via a relaxation procedure. The obtained models are then systematically validated on four well-established benchmark problems including a double diffusive Rayleigh Bénard convection of humid air, 2D and 3D thermal moist rising bubble under convective atmospheric environment as well as a shallow cumulus convection in framework of large-eddy simulation.
Yongliang Feng, Pierre Boivin, Jérome Jacob, Pierre Sagaut. Hybrid recursive regularized lattice Boltzmann simulation of humid air with application to meteorological flows. Physical Review E , American Physical Society (APS), 2019. ⟨hal-02265484⟩
G. Farag, Pierre Boivin, P. Sagaut. Interaction of two-dimensional spots with a heat releasing/absorbing shock wave: linear interaction approximation results. Journal of Fluid Mechanics, Cambridge University Press (CUP), 2019, 871, pp.865-895. ⟨10.1017/jfm.2019.324⟩. ⟨hal-02142649⟩ Plus de détails...
The canonical interaction between a two-dimensional weak Gaussian disturbance (en-tropy spot, density spot, weak vortex) with an exothermic/endothermic planar shock wave is studied via the Linear Interaction Approximation. To this end, a unified framework based on an extended Kovasznay decomposition that simultaneously accounts for non-acoustic density disturbances along with a poloidal-toroidal splitting of the vorticity mode and for heat-release is proposed. An extended version of Chu's definition for the energy of disturbances in compressible flows encompassing multi-component mixtures of gases is also proposed. This new definition precludes spurious non-normal phenomena when computing the total energy of extended Kovasznay modes. Detailed results are provided for three cases, along with fully general expressions for mixed solutions that combine incoming vortical, entropy and density disturbances.
G. Farag, Pierre Boivin, P. Sagaut. Interaction of two-dimensional spots with a heat releasing/absorbing shock wave: linear interaction approximation results. Journal of Fluid Mechanics, Cambridge University Press (CUP), 2019, 871, pp.865-895. ⟨10.1017/jfm.2019.324⟩. ⟨hal-02142649⟩
Xi Deng, Pierre Boivin, Feng Xiao. A new formulation for two-wave Riemann solver accurate at contact interfaces. Physics of Fluids, American Institute of Physics, 2019, 31 (4), pp.046102. ⟨10.1063/1.5083888⟩. ⟨hal-02100764⟩ Plus de détails...
This study proposes a new formulation for Harten, Lax and van Leer (HLL) type Riemann solver which is capable of solving contact discontinuities accurately but with robustness for strong shock. It is well known that the original HLL, which has incomplete wave structures, is too dissipative to capture contact disconti-nuities accurately. On the other side, contact-capturing approximate Riemann solvers such as Harten, Lax and van Leer with Contact (HLLC) usually suffer from spurious solutions, also called carbuncle phenomenon, for strong shock. In this work a new accurate and robust HLL-type formulation, so-called HLL-BVD (HLL Riemann solver with BVD) is proposed by modifying the original HLL with BVD (boundary variation dimin-shing) algorithm. Instead of explicitly recovering the complete wave structures like the way of HLLC, the proposed method restores the missing contact with a jump-like function. The capability of solving contact discontinuities is further improved by minimizing the inherent dissipation term in HLL. Without modifying the original incomplete wave structures of HLL, the robustness for strong shock has been reserved. Thus the proposed method is free from shock instability problem. The accuracy and robustness of the new method are demonstrated through solving several one-and two-dimensional tests. Results indicate that the new formulation based on two-wave HLL-type Riemann solver is not only capable of capturing contact waves more accurately than the original HLL or HLLC but, most importantly, is free form carbuncle instability for strong shock.
Xi Deng, Pierre Boivin, Feng Xiao. A new formulation for two-wave Riemann solver accurate at contact interfaces. Physics of Fluids, American Institute of Physics, 2019, 31 (4), pp.046102. ⟨10.1063/1.5083888⟩. ⟨hal-02100764⟩
Pierre Boivin, M.A. Cannac, O. Le Metayer. A thermodynamic closure for the simulation of multiphase reactive flows. International Journal of Thermal Sciences, Elsevier, 2019, 137, pp.640-649. ⟨hal-01981954⟩ Plus de détails...
A simple thermodynamic closure for the simulation of multiphase reactive flows is presented. It combines a fully explicit thermodynamic closure appropriate for weakly thermal multiphase flow simulations, with the classical variable heat capacity ideal gas thermodynamic closure, commonly used for reactive flows simulations. Each liquid and gas component is assumed to follow the recent Noble-Abel Stiffened Gas equation of state, fully described by a set of five parameters. A new method for setting these parameters is presented and validated through comparisons with NIST references. Comparisons with a well-known cubic equation of state, Soave-Redlich-Kwong, are also included. The Noble-Abel Stiffened-Gas equation of state is then extended as to cope with variable heat capacity, to make the mixture ther-modynamic closure appropriate for multiphase reactive flows.
Pierre Boivin, M.A. Cannac, O. Le Metayer. A thermodynamic closure for the simulation of multiphase reactive flows. International Journal of Thermal Sciences, Elsevier, 2019, 137, pp.640-649. ⟨hal-01981954⟩
Journal: International Journal of Thermal Sciences
Yongliang Feng, Pierre Boivin, Jérome Jacob, Pierre Sagaut. Hybrid recursive regularized thermal lattice Boltzmann model for high subsonic compressible flows. Journal of Computational Physics, Elsevier, 2019, 394, pp.82-99. ⟨hal-02142837⟩ Plus de détails...
A thermal lattice Boltzmann model with a hybrid recursive regularization (HRR) collision operator is developed on standard lattices for simulation of subsonic and sonic compressible flows without shock. The approach is hybrid: mass and momentum conservation equations are solved using a lattice Boltzmann solver, while the energy conservation is solved under entropy form with a finite volume solver. The defect of Galilean invariance related to Mach number is corrected by the third order equilibrium distribution function , supplemented by an additional correcting term and hybrid recursive regularization. The proposed approach is assessed considering the simulation of i) an isentropic vortex convection, ii) a two dimensional acoustic pulse and iii) non-isothermal Gaussian pulse with Ma number in range of 0 to 1. Numerical simulations demonstrate that the flaw in Galilean invari-ance is effectively eliminated by the compressible HRR model. At last, the compressible laminar flows over flat plate at Ma number of 0.3 and 0.87, Reynolds number of 10 5 are considered to validate the capture of viscous and diffusive effects.
Yongliang Feng, Pierre Boivin, Jérome Jacob, Pierre Sagaut. Hybrid recursive regularized thermal lattice Boltzmann model for high subsonic compressible flows. Journal of Computational Physics, Elsevier, 2019, 394, pp.82-99. ⟨hal-02142837⟩
Yongliang Feng, Muhammad Tayyab, Pierre Boivin. A Lattice-Boltzmann model for low-Mach reactive flows. Combustion and Flame, Elsevier, 2018, 196, pp.249 - 254. ⟨10.1016/j.combustflame.2018.06.027⟩. ⟨hal-01832640⟩ Plus de détails...
A new Lattice-Boltzmann model for low-Mach reactive flows is presented. Based on standard lattices, the model is easy to implement, and is the first, to the authors' knowledge, to pass the classical freely propagating flame test case as well as the counterflow diffusion flame, with strains up to extinction. For this presentation, simplified transport properties are considered, each species being assigned a separate Lewis number. In addition, the gas mixture is assumed to be calorically perfect. Comparisons with reference solutions show excellent agreement for mass fraction profiles, flame speed in premixed mixtures, as well as maximum temperature dependence with strain rate in counterflow diffusion flames.
Yongliang Feng, Muhammad Tayyab, Pierre Boivin. A Lattice-Boltzmann model for low-Mach reactive flows. Combustion and Flame, Elsevier, 2018, 196, pp.249 - 254. ⟨10.1016/j.combustflame.2018.06.027⟩. ⟨hal-01832640⟩
Pierre Boivin, Forman Williams. Extension of a wide-range three-step hydrogen mechanism to syngas. Combustion and Flame, Elsevier, 2018, 196, pp.85-87. ⟨10.1016/j.combustflame.2018.05.034⟩. ⟨hal-02112081⟩ Plus de détails...
Previously we have shown how a single species X can be introduced, representing either HO 2 for high-temperature ignition or H 2 O 2 for low-temperature ignition, to develop an algorithm that covers the entire range of ignition, flame-propagation, and combustion conditions, without a significant degradation of accuracy, for hydrogen-air systems. By adding relevant CO chemistry to the hydrogen chemistry, this same approach can be applied to derive a comparably useful four-step reduced-chemistry description for syngas blends that have small enough concentrations of methane, other hydrocarbons , or other reactive species to be dominated by the elementary steps of the H 2 /CO system. The present communication reports the resulting extended algorithm. This work begins with the elementary steps of the detailed chemistry as listed in Table 1. We shall employ the numbering of the steps as given in the table, which identifies the 8 steps that are considered to be reversible and gives fitted parameters for the reverses of those steps.
Pierre Boivin, Forman Williams. Extension of a wide-range three-step hydrogen mechanism to syngas. Combustion and Flame, Elsevier, 2018, 196, pp.85-87. ⟨10.1016/j.combustflame.2018.05.034⟩. ⟨hal-02112081⟩
Alexandre Chiapolino, Pierre Boivin, Richard Saurel. A simple phase transition relaxation solver for liquid-vapor flows. International Journal for Numerical Methods in Fluids, Wiley, 2017, 83 (7), pp.583-605. ⟨10.1002/fld.4282⟩. ⟨hal-01359203⟩ Plus de détails...
Determining liquid-vapor phase equilibrium is often required in multiphase flow computations. Existing equilibrium solvers are either accurate but computationally expensive, or cheap but inaccurate. The present paper aims at building a fast and accurate specific phase equilibrium solver, specifically devoted to unsteady multiphase flow computations. Moreover, the solver is efficient at phase diagram bounds, where non-equilibrium pure liquid and pure gas are present. It is systematically validated against solutions based on an accurate (but expensive) solver. Its capability to deal with cavitating, evaporating and condensing two-phase flows is highlighted on severe test problems both 1D and 2D.
Alexandre Chiapolino, Pierre Boivin, Richard Saurel. A simple phase transition relaxation solver for liquid-vapor flows. International Journal for Numerical Methods in Fluids, Wiley, 2017, 83 (7), pp.583-605. ⟨10.1002/fld.4282⟩. ⟨hal-01359203⟩
Journal: International Journal for Numerical Methods in Fluids
Alexandre Chiapolino, Pierre Boivin, Richard Saurel. A simple and fast phase transition relaxation solver for compressible multicomponent two-phase flows. Computers and Fluids, Elsevier, 2017, 150, pp.31 - 45. ⟨10.1016/j.compfluid.2017.03.022⟩. ⟨hal-01502389⟩ Plus de détails...
The present paper aims at building a fast and accurate phase transition solver dedicated to unsteady multiphase flow computations. In a previous contribution (Chiapolino et al. 2017), such a solver was successfully developed to compute thermodynamic equilibrium between a liquid phase and its corresponding vapor phase. The present work extends the solver's range of application by considering a multicomponent gas phase instead of pure vapor, a necessary improvement in most practical applications. The solver proves easy to implement compared to common iterative procedures, and allows systematic CPU savings over 50%, at no cost in terms of accuracy. It is validated against solutions based on an accurate but expensive iterative solver. Its capability to deal with cavitating, evaporating and condensing two-phase flows is highlighted on severe test problems both 1D and 2D.
Alexandre Chiapolino, Pierre Boivin, Richard Saurel. A simple and fast phase transition relaxation solver for compressible multicomponent two-phase flows. Computers and Fluids, Elsevier, 2017, 150, pp.31 - 45. ⟨10.1016/j.compfluid.2017.03.022⟩. ⟨hal-01502389⟩
Richard Saurel, Olivier Le Métayer, Pierre Boivin. From Cavitating to Boiling Flows. d'Agostino L., Salvetti M.; CISM International Centre for Mechanical Sciences (Courses and Lectures). Cavitation Instabilities and Rotordynamic Effects in Turbopumps and Hydroturbines , 575, Springer pp.259-282 2017, 978-3-319-49717-4. ⟨hal-01678361⟩ Plus de détails...
A flow model is derived for the numerical simulation of interfacial flows with phase transition. The model arises from the classical multi-component Euler equations, but is associated to a non-classical thermodynamic closure: each phase is compressible and evolves in its own subvolume, with phases sharing common pressure, velocity and temperature, leading to non-trivial thermodynamic relations for the mixture. Phase transition is made possible through the introduction of Gibbs free energy relaxation terms in the equations. Capillary effects and heat conduction—essential in boiling flows—are introduced as well. The resulting multi-phase flow model is hyperbolic, valid for arbitrary density jumps at interfaces as well as arbitrary flow speeds. Its capabilities are illustrated successively through examples of nozzle induced cavitation and heated wall induced boiling.
Richard Saurel, Olivier Le Métayer, Pierre Boivin. From Cavitating to Boiling Flows. d'Agostino L., Salvetti M.; CISM International Centre for Mechanical Sciences (Courses and Lectures). Cavitation Instabilities and Rotordynamic Effects in Turbopumps and Hydroturbines , 575, Springer pp.259-282 2017, 978-3-319-49717-4. ⟨hal-01678361⟩
Richard Saurel, Pierre Boivin, Olivier Le Métayer. A general formulation for cavitating, boiling and evaporating flows. Computers and Fluids, Elsevier, 2016, 128, pp.53-64. ⟨10.1016/j.compfluid.2016.01.004⟩. ⟨hal-01277179⟩ Plus de détails...
A flow model is derived for the numerical simulation of multi-phase flows with phase transition. The model arises from the classical multi-component Euler equations, but is associated to a non-classical thermodynamic closure: each phase is compressible and evolves in its own subvolume, with phases sharing common pressure, velocity and temperature, leading to non-trivial thermodynamic relations for the mixture. Phase transition is made possible through the introduction of Gibbs free energy relaxation terms in the equations. Capillary effects and heat conduction – essential in boiling flows – are introduced as well. The resulting multi-phase flow model is hyperbolic, valid for arbitrary density jumps at interfaces as well as arbitrary flow speeds. Its capabilities are illustrated successively through examples of nozzle induced cavitation, a high-speed evaporating liquid jet, and heated wall induced boiling.
Richard Saurel, Pierre Boivin, Olivier Le Métayer. A general formulation for cavitating, boiling and evaporating flows. Computers and Fluids, Elsevier, 2016, 128, pp.53-64. ⟨10.1016/j.compfluid.2016.01.004⟩. ⟨hal-01277179⟩