A reduced model for compressible viscous heat-conducting multicomponent flows

In the present paper we propose a reduced temperature non-equilibrium model for simulating multicomponent flows with inter-phase heat transfer, diffusion processes (including the viscosity and the heat conduction) and external energy sources. We derive three equivalent formulations for the proposed model. The first formulation consists of balance equations for partial densities, the mixture momentum, the mixture total energy, and phase volume fractions. The second formulation is symmetric and obtained by replacing the equations for the mixture total energy and volume fractions in the first formulation with balance equations for the phase total energy. Replacing one of the phase total energy equation of the second formulation with the mixture total energy equation gives the third formulation. All the three formulations assume velocity and pressure equilibrium across the material interface. These equivalent forms provide different physical perspectives and numerical conveniences. Temperature equilibration and continuity across the material interfaces are achieved with the instantaneous thermal relaxation. Temperature equilibrium is maintained during the heat conduction process. The proposed models are proved to respect the thermodynamical laws. For numerical solution, the model is split into a hyperbolic partial differential equation (PDE) system and parabolic PDE systems. The former is solved with the high-order Godunov finite volume method that ensures the pressure-velocitytemperature (PVT) equilibrium conduction. The parabolic PDEs are solved with both the implicit and the explicit locally iterative method (LIM) based on Chebyshev parameters. Numerical results are presented for several multicomponent flow problems with diffusion processes. Furthermore, we apply the proposed model to simulate the target ablation problem that is of significance to inertial confinement fusion. Comparisons with one-temperature models in literature demonstrate the ability to maintain the PVT property and superior convergence performance of the proposed model in solving multicomponent problems with diffusions.

[1]  Richard Saurel,et al.  A multiphase model for compressible flows with interfaces, shocks, detonation waves and cavitation , 2001, Journal of Fluid Mechanics.

[2]  Eric Johnsen,et al.  Maintaining interface equilibrium conditions in compressible multiphase flows using interface capturing , 2015, J. Comput. Phys..

[3]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[4]  Tim Colonius,et al.  Finite-volume WENO scheme for viscous compressible multicomponent flows , 2014, J. Comput. Phys..

[5]  Peter D. M. Spelt,et al.  Simulations of viscous and compressible gas-gas flows using high-order finite difference schemes , 2018, J. Comput. Phys..

[6]  Robert I. Nigmatulin,et al.  Dynamics of multiphase media , 1991 .

[7]  Rao V. Garimella,et al.  A comparative study of interface reconstruction methods for multi-material ALE simulations , 2010 .

[8]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[9]  Eleuterio F. Toro,et al.  Conservative hyperbolic formulation for compressible two-phase flow with different phase pressures and temperatures , 2007 .

[10]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[11]  Chi-Wang Shu Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[12]  Boniface Nkonga,et al.  Sharpening diffuse interfaces with compressible fluids on unstructured meshes , 2017, J. Comput. Phys..

[13]  Grégoire Allaire,et al.  A five-equation model for the simulation of interfaces between compressible fluids , 2002 .

[14]  Halvor Lund,et al.  A Hierarchy of Relaxation Models for Two-Phase Flow , 2012, SIAM J. Appl. Math..

[15]  Jean-Marc Hérard,et al.  A three-phase flow model , 2007, Math. Comput. Model..

[16]  M. Baer,et al.  A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials , 1986 .

[17]  R. Abgrall,et al.  A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows , 1999 .

[18]  D. Stewart,et al.  Two-phase modeling of deflagration-to-detonation transition in granular materials: Reduced equations , 2001 .

[19]  Richard Saurel,et al.  A compressible flow model with capillary effects , 2005 .

[20]  Eric Johnsen,et al.  Preventing numerical errors generated by interface-capturing schemes in compressible multi-material flows , 2012, J. Comput. Phys..

[21]  Ye Wen-hu NUMERICAL STUDY OF LASER ABLATIVE RAYLEIGH-TAYLOR INSTABILITY , 1999 .

[22]  W. Sutherland LII. The viscosity of gases and molecular force , 1893 .

[23]  Carlos Pantano,et al.  Diffuse-Interface Capturing Methods for Compressible Two-Phase Flows , 2018 .

[24]  Igor Menshov,et al.  Eulerian Model for Simulating Multi-Fluid Flows with an Arbitrary Number of Immiscible Compressible Components , 2020, J. Sci. Comput..

[25]  V. Zhukov Explicit methods of numerical integration for parabolic equations , 2011 .

[26]  L. Spitzer,et al.  TRANSPORT PHENOMENA IN A COMPLETELY IONIZED GAS , 1953 .

[27]  Hervé Guillard,et al.  A five equation reduced model for compressible two phase flow problems , 2005 .

[28]  V. Zhukov,et al.  Explicit time integration of the Navier-Stokes equations using the local iteration method , 2019, Keldysh Institute Preprints.

[29]  R. Bonazza,et al.  Shock-Bubble Interactions , 2011 .

[30]  Jiaquan Gao,et al.  How to prevent pressure oscillations in multicomponent flow calculations , 2000, Proceedings Fourth International Conference/Exhibition on High Performance Computing in the Asia-Pacific Region.

[31]  Eric Johnsen,et al.  Implementation of WENO schemes in compressible multicomponent flow problems , 2005, J. Comput. Phys..

[32]  Ben Thornber,et al.  A five-equation model for the simulation of miscible and viscous compressible fluids , 2018, J. Comput. Phys..

[33]  R. P. Drake,et al.  The time scale for the transition to turbulence in a high Reynolds number, accelerated flow , 2003 .

[34]  Richard Saurel,et al.  Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures , 2009, J. Comput. Phys..

[35]  Boniface Nkonga,et al.  Towards the direct numerical simulation of nucleate boiling flows , 2014 .

[36]  Barry Koren,et al.  A new formulation of Kapila's five-equation model for compressible two-fluid flow, and its numerical treatment , 2010, J. Comput. Phys..

[37]  Stefano Atzeni,et al.  The Physics of Inertial Fusion: Beam Plasma Interaction, Hydrodynamics, Hot Dense Matter , 2004 .

[38]  J. Haas,et al.  Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities , 1987, Journal of Fluid Mechanics.