×

Combining discrete equations method and upwind downwind-controlled splitting for non-reacting and reacting two-fluid computations: two dimensional case. (English) Zbl 1391.76582

Summary: This work deals with impermeable and permeable interfaces and the design of numerical strategies allowing multi-dimensional propagation of these interfaces on general unstructured grids. The numerical context is the (reactive) discrete equations method (DEM/RDEM) for the Baer-Nunziato type non-equilibrium multiphase model allowing a diffused interface, and meanwhile preserving the global conservation, which is of fundamental importance for studying long term combustion phenomena in large-scale geometries. Another advantage of RDEM for combustion lies in its ability to compute both deflagration and detonation, provided an appropriate reactive Riemann solver is inserted within the method. The present paper is a sequel to the recent publication [K. Tang et al., Comput. Fluids 93, 74–90 (2014; Zbl 1391.76581)] where an anti-diffusive approach and an original upwind downwind-controlled splitting method (UDCS) were combined with the 1D formulation of the DEM and RDEM. The method successfully developed in 1D for computing inert interfaces (e.g., impermeable water gas shock tube problem) and flame interfaces (e.g. Chapman-Jouguet deflagration and strong detonation wave) with excellent robustness and accuracy properties is extended here to two dimensional problems. The proposed low- and anti-diffusive versions of the multi-D UDCS strategy form an original contribution to the modeling of multifluid flows on unstructured grids. This multi-D extension relies on a general derivation of the downwind factors involved in the formulation of UDCS. In particular, the proposed UDCS anti-diffusive algorithm represents a new alternative to the “Extended-Vofire” solver [V. Faucher and S. Kokh, “Extended Vofire algorithm for fast transient fluid-structure dynamics with liquid-gas flows and interfaces”, J. Fluids Struct. 39, 102–125 (2013; doi:10.1016/j.jfluidstructs.2013.02.014)] for unstructured meshes. Numerical experiments performed for non-reacting gas-gas and liquid-gas shock bubble interactions as well as for a model combustion problem demonstrate the combination of DEM/RDEM with UDCS yields excellent robustness/accuracy properties. Some remaining issues linked to the modeling of flame propagation in multi-dimensional cases are eventually discussed.

MSC:

76M25 Other numerical methods (fluid mechanics) (MSC2010)
76V05 Reaction effects in flows
80A25 Combustion

Citations:

Zbl 1391.76581

Software:

HLLC
Full Text: DOI

References:

[1] Abgrall, R.; Saurel, R., Discrete equations for physical and numerical compressible multiphase mixtures, J Comput Phys, 186, 361-396, (2003) · Zbl 1072.76594
[2] Baer, M.; Nunziato, J., A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials, Int J Multiphase Flow, 12, 6, 861-889, (1986) · Zbl 0609.76114
[3] Banks, J. W.; Schwendeman, D. W.; Kapila, A. K.; Henshaw, W. D., A high-resolution Godunov method for compressible multi-material flow on overlapping grids, J Comput Phys, 223, 1, 262-297, (2007) · Zbl 1163.76032
[4] Beccantini, A.; Studer, E., The reactive Riemann problem for thermally perfect gases at all combustion regimes, Int J Numer Methods Fluids, 64, 269-313, (2010) · Zbl 1272.80010
[5] Borghi R, Destriau M. La combustion et les flames. Editions Technip; 1995.
[6] Chertock, A.; Karni, S.; Kurganov, A., Interface tracking method for compressible multifluids, ESAIM: Math Modell Numer Anal, 42, 991-1019, (2008) · Zbl 1232.76030
[7] Chinnayya, A.; Daniel, E.; Saurel, R., Modelling detonation waves in heterogeneous energetic materials, J Comput Phys, 196, 2, 490-538, (2004) · Zbl 1109.76335
[8] Delanaye M. Polynomial reconstruction finite volume schemes for the compressible Euler and Navier-Stokes equations on unstructured adaptive grids. PhD thesis, Université de Liège; 1998.
[9] Delanaye M, Essers JA. An accurate finite volume scheme for Euler and Navier-Stokes equations on unstructured adaptive grids. In: 12th computational fluid dynamics conference (AIAA CFD Conf.), San Diego, June 1995. American Institute of Aeronautics and Astronautics; 1995.
[10] Després, B.; Lagoutière, F.; Labourasse, E.; Marmajou, I., An antidissipative transport scheme on unstructured meshes for multicomponent flows, Int J Finite Vol, 7, 30-65, (2010)
[11] Drew, D. A.; Passman, S. L., Theory of multicomponent fluids, (Applied mathematical sciences, vol. 135, (1998), Springer New York) · Zbl 0919.76003
[12] EUROPLEXUS. A computer program for the finite element simulation of fluid-structure systems under transient dynamic loading, 28 May; 2014. <http://www-epx.cea.fr>.
[13] EUROPLEXUS. Numerical simulation of hydrogen explosion in a reactor with Europleuxs, 28 May; 2014. <http://www-epx.cea.fr/index.php/what-is-epx/27-videos#Explo_H2>.
[14] Faucher V, Galon P, Beccantini A, Crouzet F, Debaud F, Gautier T. Hybrid parallel strategy for the simulation of fast transient accidental situations at reactor scale. In: SNA & MC 2013 - joint international conference on supercomputing in nuclear applications + Monte Carlo, October 27-31, 2013, La Cité des Sciences et de l’Industrie, Paris, France; 2013.
[15] Faucher, V.; Kokh, S., Extended vofire algorithm for fast transient fluid-structure dynamics with liquid-gas flows and interfaces, J Fluids Struct, 39, 102-125, (2013)
[16] Franquet, E.; Perrier, V., Runge-Kutta discontinuous Galerkin method for interface flows with a maximum preserving limiter, Comput Fluids, 65, 2-7, (2012) · Zbl 1365.76116
[17] Franquet, E.; Perrier, V., Runge-Kutta discontinuous Galerkin method for the approximation of Baer and Nunziato type multiphase models, J Comput Phys, 231, 11, 4096-4141, (2012) · Zbl 1426.76530
[18] Franquet, E.; Perrier, V., Runge-Kutta discontinuous Galerkin method for reactive multiphase flows, Comput Fluids, 83, 157-163, (2013) · Zbl 1290.76159
[19] Haas, J.-F.; Sturtevant, B., Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities, J Fluid Mech, 181, 41-76, (1987)
[20] Ishii, M.; Hibiki, T., Thermo-fluid dynamics of two-phase flow, (2011), Springer New York, NY · Zbl 1209.76001
[21] Jameson, A., A perspective on computational algorithms for aerodynamic analysis and design, Progr Aerospace Sci, 37, 197-243, (2001)
[22] Kokh, S.; Lagoutière, F., An anti-diffusive numerical scheme for the simulation of interfaces between compressible fluids by means of a five-equation model, J Comput Phys, 229, 8, 2773-2809, (2010) · Zbl 1302.76129
[23] Kreeft, J. J.; Koren, B., A new formulation of kapila’s five-equation model for compressible two-fluid flow, and its numerical treatment, J Comput Phys, 229, 18, 6220-6242, (2010) · Zbl 1425.76216
[24] Kuhl, A. L.; Kamel, M. M.; Oppenheim, A. K., Pressure waves generated by steady flames, Sympos (Int) Combust, 14, 1, 1201-1215, (1973)
[25] Le Métayer, O.; Massoni, J.; Saurel, R., Modelling evaporation fronts with reactive Riemann solvers, J Comput Phys, 205, 567-610, (2005) · Zbl 1088.76051
[26] Nonomura, T.; Kitamura, K.; Fujii, K., A simple interface sharpening technique with a hyperbolic tangent function applied to compressible two-fluid modeling, J Comput Phys, 258, 95-117, (2014) · Zbl 1349.76514
[27] Nourgaliev, R. R.; Dinh, T. N.; Theofanous, T. G., Adaptive characteristics-based matching for compressible multifluid dynamics, J Comput Phys, 213, 2, 500-529, (2006) · Zbl 1136.76396
[28] Quirk, J. J.; Karni, S., On the dynamics of a shock-bubble interaction, J Fluid Mech, 318, 129-163, (1996) · Zbl 0877.76046
[29] Saurel, R.; Abgrall, R., A multiphase Godunov method for compressible multifluid and multiphase flows, J Comput Phys, 150, 425-467, (1999) · Zbl 0937.76053
[30] Saurel, R.; Abgrall, R., A simple method for compressible multifluid flows, SIAM J Sci Comput, 21, 3, 1115, (1999) · Zbl 0957.76057
[31] Shyue, K.-M., A wave-propagation based volume tracking method for compressible multicomponent flow in two space dimensions, J Comput Phys, 215, 219-244, (2006) · Zbl 1140.76401
[32] Studer E, Beccantini A, Kudriakov S, Velikorodny A. Hydrogen combustion modelling in large scale geometries. In: ICONE 21 - international conference on nuclear engineering, July 29 - August 2, 2013, Chengdu, China; 2013.
[33] Tang K. Combining discrete equations method and upwind downwind-controlled splitting for non-reacting and reacting two-fluid computations. PhD thesis, Grenoble Institute of Technology; 2012. · Zbl 1391.76582
[34] Tang K, Beccantini A, Corre C. Combining discrete equations method and upwind downwind-controlled splitting for non-reacting and reacting two-fluid computations. In: Seventh international conference on computational fluid dynamics (ICCFD7), Big Island, Hawaii, July 9-13; 2012. · Zbl 1391.76582
[35] Tang, K.; Beccantini, A.; Corre, C., Combining discrete equations method and upwind downwind-controlled splitting for non-reacting and reacting two-fluid computations: one dimensional case, Comput Fluids, 93, 74-90, (2014) · Zbl 1391.76581
[36] Tokareva, S. A.; Toro, E. F., HLLC-type Riemann solver for the Baer-Nunziato equations of compressible two-phase flow, J Comput Phys, 229, 10, 3573-3604, (2010) · Zbl 1391.76440
[37] Ullah, M. A.; Gao, W.; Mao, D.-K., Towards front-tracking based on conservation in two space dimensions III, tracking interfaces, J Comput Phys, 242, 268-303, (2013) · Zbl 1314.65113
[38] Velikorodny A, Studer E, Kudriakov S, Beccantini A. Combustion modeling in large scale volumes. In: ICHS 2013 - international conference on hydrogen safety, September 9-11, 2013, Brussels - Belgium; 2013.
[39] Zein, A.; Hantke, M.; Warnecke, G., Modeling phase transition for compressible two-phase flows applied to metastable liquids, J Comput Phys, 229, 8, 2964-2998, (2010) · Zbl 1307.76079
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.