A new formulation of Kapila's five-equation model for compressible two-fluid flow, and its numerical treatment

被引:74
作者
Kreeft, Jasper J. [1 ]
Koren, Barry [2 ,3 ]
机构
[1] Delft Univ Technol, Fac Aerosp Engn, NL-2600 GB Delft, Netherlands
[2] Ctr Wiskunde & Informat, NL-1090 GB Amsterdam, Netherlands
[3] Leiden Univ, Math Inst, NL-2300 RA Leiden, Netherlands
关键词
Compressible two-fluid flow; Two-fluid mixture flow; Exchange terms; Momentum and energy exchange; Interface capturing; Approximate Riemann solver; Shock tube; Shock-bubble-interaction; RELAXATION-PROJECTION METHOD; TO-DETONATION TRANSITION; 2-PHASE FLOW; CONSERVATION-LAWS; EULER EQUATIONS; SCHEMES; FLUIDS; SIMULATION; INTERFACES; DYNAMICS;
D O I
10.1016/j.jcp.2010.04.025
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
A new formulation of Kapila's five-equation model for inviscid, non-heat-conducting, compressible two-fluid flow is derived, together with an appropriate numerical method. The new formulation uses flow equations based on conservation laws and exchange laws only. The two fluids exchange momentum and energy, for which exchange terms are derived from physical laws. All equations are written as a single system of equations in integral form. No equation is used to describe the topology of the two-fluid flow. Relations for the Riemann invariants of the governing equations are derived, and used in the construction of an Osher-type approximate Riemann solver. A consistent finite-volume discretization of the exchange terms is proposed. The exchange terms have distinct contributions in the cell interior and at the cell faces. For the exchange-term evaluation at the cell faces, the same Riemann solver as used for the flux evaluation is exploited. Numerical results are presented for two-fluid shock-tube and shock-bubble-interaction problems, the former also for a two-fluid mixture case. All results show good resemblance with reference results. (C) 2010 Elsevier Inc. All rights reserved.
引用
收藏
页码:6220 / 6242
页数:23
相关论文
共 29 条
[1]   How to prevent pressure oscillations in multicomponent flow calculations: A quasi conservative approach [J].
Abgrall, R .
JOURNAL OF COMPUTATIONAL PHYSICS, 1996, 125 (01) :150-160
[2]   Asymptotic expansion of a multiscale numerical scheme for compressible multiphase flow [J].
Abgrall, R ;
Perrier, V .
MULTISCALE MODELING & SIMULATION, 2006, 5 (01) :84-115
[3]   A five-equation model for the simulation of interfaces between compressible fluids [J].
Allaire, G ;
Clerc, S ;
Kokh, S .
JOURNAL OF COMPUTATIONAL PHYSICS, 2002, 181 (02) :577-616
[4]  
[Anonymous], 1999, Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science
[5]  
[Anonymous], 2002, Level Set Methods and Dynamic Implicit Surfaces
[6]   A 2-PHASE MIXTURE THEORY FOR THE DEFLAGRATION-TO-DETONATION TRANSITION (DDT) IN REACTIVE ANTIGRANULOCYTES-MATERIALS [J].
BAER, MR ;
NUNZIATO, JW .
INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 1986, 12 (06) :861-889
[7]   Numerical simulation of the homogeneous equilibrium model for two-phase flows [J].
Clerc, S .
JOURNAL OF COMPUTATIONAL PHYSICS, 2000, 161 (01) :354-375
[8]  
Guillard H., 2006, ECCOMAS CFD
[9]   INTERACTION OF WEAK SHOCK-WAVES WITH CYLINDRICAL AND SPHERICAL GAS INHOMOGENEITIES [J].
HAAS, JF ;
STURTEVANT, B .
JOURNAL OF FLUID MECHANICS, 1987, 181 :41-76
[10]   MULTIPLE GRID AND OSHER SCHEME FOR THE EFFICIENT SOLUTION OF THE STEADY EULER EQUATIONS [J].
HEMKER, PW ;
SPEKREIJSE, SP .
APPLIED NUMERICAL MATHEMATICS, 1986, 2 (06) :475-493