This work presents a mathematical model for the two-phase flows in the mortar systems and demonstrates the application of approximate Riemann solver on such model. The mathematical model for the two-phase gas-dynamical processes in the mortar tube consists of a system of first-order, nonlinear coupled partial differential equations with inhomogeneous terms. The model poses an initial value problem with discontinuous initial and boundary conditions that arise due to the design complexity and nonuniformity of granular propellant distribution in the mortar tube. The governing equations in this model possess characteristics of the Riemann problem. Therefore, a high-resolution Godunov-type shock-capturing approach was used to address the formation of flow structure such as shock waves, contact discontinuities, and rarefaction waves. A linearized approximate Riemann solver based on the Roe–Pike method was modified for the two-phase flows to compute fully nonlinear wave interactions and to directly provide upwinding properties in the scheme. An entropy fix based on Harten–Heyman method was used with van Leer flux limiter for total variation diminishing. The three-dimensional effects were simulated by incorporating an unsplit multidimensional wave propagation method, which accounted for discontinuities traveling in both normal and oblique coordinate directions. A mesh generation algorithm was developed to account for the projectile motion and coupled with the approximate Riemann solver. The numerical method was verified by using exact solutions of three test problems. The specific system considered in this work is a 120 mm mortar system, which contains an ignition cartridge that discharges hot gas-phase products and unburned granular propellants into the mortar tube through multiple vent-holes on its surface. The model for the mortar system was coupled with the solution of the transient gas-dynamic behavior in the ignition cartridge. The numerical results were validated with experimental data. Based on the close comparison between the calculated results and test data, it was found that the approximate Riemann solver is a suitable method for studying the two-phase combustion processes in mortar systems.

1.
Deledicque
,
V.
, and
Papalexandris
,
M. V.
, 2007, “
An Exact Riemann Solver for Compressible Two-Phase Flow Models Containing Non-Conservative Products
,”
J. Comput. Phys.
0021-9991,
222
, pp.
217
245
.
2.
Ishii
,
M.
, 1975,
Thermo-Fluid Dynamic Theory of Two-Phase Flows
,
Eyrolles
,
Paris
.
3.
Drew
,
D. A.
, and
Passman
,
S. L.
, 1998,
Theory of Multicomponent Fluids
,
Springer
,
New York
.
4.
Enwald
,
H.
,
Peirano
,
E.
, and
Almstedt
,
A. E.
, 1986, “
Eulerian Two-Phase Flow Theory Applied to Fluidization
,”
Int. J. Multiphase Flow
0301-9322,
222
, pp.
21
66
.
5.
Gough
,
P. S.
, and
Zwarts
,
F. J.
, 1979, “
Modeling Heterogeneous Two-Phase Reacting Flow
,”
AIAA J.
0001-1452,
17
(
1
), pp.
17
25
.
6.
Baer
,
M. R.
, and
Nunziato
,
J. W.
, 1986, “
A Two-phase Mixture Theory for the Deflagration-to-Detonation Transition (DDT) in Reactive Granular Materials
,”
Int. J. Multiphase Flow
0301-9322,
12
(
6
), pp.
861
889
.
7.
Powers
,
J. M.
,
Stewart
,
D. S.
, and
Krier
,
H. K.
, 1990, “
Theory of Two-Phase Detonation—Part II: Modeling
,”
Combust. Flame
0010-2180,
80
, pp.
264
279
.
8.
Saurel
,
R.
, and
Abgrall
,
R.
, 1999, “
A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows
,”
J. Comput. Phys.
0021-9991,
150
, pp.
425
467
.
9.
Papalexandris
,
M. V.
, 2004, “
Numerical Simulations of Detonations in Mixtures of Gases and Solid Particles
,”
J. Fluid Mech.
0022-1120,
507
, pp.
95
142
.
10.
Kuo
,
K. K.
, 1986,
Principles of Combustion
,
Wiley
,
New York
, Chap. 8.
11.
Godunov
,
S. K.
, 1959, “
A Finite Difference Method for the Computation of Discontinuous Solutions of the Equations of Fluid Dynamics
,”
Mat. Sb.
0368-8666,
47
, pp.
357
393
.
12.
Roe
,
P. L.
, and
Pike
,
J.
, 1984, “
Efficient Construction and Utilization of Approximate Riemann Solutions
,”
Computing Methods in Applied Science and Engineering
,
North-Holland
,
Amsterdam
.
13.
Acharya
,
R.
, 2009, “
Modeling and Simulation of Interior Ballistic Processes in 120 mm Mortar System
,” Ph.D. thesis, Pennsylvania State University.
14.
Roe
,
P. L.
, 1981, “
Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes
,”
J. Comput. Phys.
0021-9991,
43
, pp.
357
372
.
15.
Harten
,
A.
, and
Hyman
,
J. M.
, 1983, “
Self Adjusting Grid Methods for One-Dimensional Hyperbolic Conservation Laws
,”
J. Comput. Phys.
0021-9991,
50
, pp.
235
269
.
16.
Sweby
,
P. K.
, 1982, “
Shock Capturing Schemes
,” Ph.D. thesis, Department of Mathematics, University of Reading, UK.
17.
Harten
,
A.
, 1983, “
High Resolution Schemes for Hyperbolic Conservation Laws
,”
J. Comput. Phys.
0021-9991,
49
, pp.
357
393
.
18.
van Leer
,
B.
, 1992, “
Progress in Multi-Dimensional Upwind Differencing
,”
NASA
Technical Report No. CR-l89708/ICASE 92-43.
19.
Chakravarthy
,
S. R.
, and
Osher
,
S.
, 1983, “
Numerical Experiments With the Osher Upwind Scheme for the Euler Equations
,”
AIAA J.
0001-1452,
21
(
9
), pp.
1241
1248
.
20.
LeVeque
,
R. J.
, 1997, “
Wave Propagation Algorithms for Multi-Dimensional Hyperbolic Systems
,”
J. Comput. Phys.
0021-9991,
131
, pp.
327
353
.
21.
Langseth
,
J. O.
, and
LeVeque
,
R. J.
, 2000, “
A Wave-Propagation Method for Three-Dimensional Hyperbolic Conservation Laws
,”
J. Comput. Phys.
0021-9991,
165
, pp.
126
166
.
22.
Toro
,
E. F.
, 1997,
Riemann Solvers and Numerical Methods for Fluid Dynamics
, 1st ed.,
Springer
,
Heidelberg
.
23.
Woodward
,
P.
, and
Colella
,
P.
, 1984, “
The Numerical Simulation of Two-Dimensional Fluid Flow With Strong Shocks
,”
J. Comput. Phys.
0021-9991,
54
, pp.
115
173
.
You do not currently have access to this content.