Numerical solution of three-dimensional heterogeneous solid propellants

We have described, for the first time, a fully coupled low Mach number numerical algorithm which can be used to investigate the combustion of composite propellants. The code uses a body fitted grid along the moving interface by means of a mapping technique. Specifically, the surface is assumed to be single valued, thus allowing the use of a mapping function that effectively maps the propagating corrugated surface into a stationary flat surface. The transformed connection conditions are then applied along the stationary flat surface, allowing second-order one-sided derivatives to be used. This has proven to be a very robust and efficient way to treat the surface and the connection conditions. Numerical tests are performed and the scheme is shown to be second-order accurate in the spatial directions as well as in time. Selected results are presented for heterogeneous propellants.

[1]  C. F. Price,et al.  A model of composite solid-propellant combustion based on multiple flames , 1970 .

[2]  John B. Bell,et al.  An Adaptive Projection Method for Unsteady, Low-Mach Number Combustion , 1998 .

[3]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[4]  C HindmarshAlan.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980 .

[5]  R. Miller.  Effects of particle size on reduced smoke propellant ballistics , 1982 .

[6]  J. Buckmaster,et al.  A numerical study of periodic sandwich propellants with oxygenated binders , 2003 .

[7]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[8]  Thomas L. Jackson,et al.  Random Packing of Heterogeneous Propellants , 2001 .

[9]  John C. Strikwerda,et al.  Finite Difference Methods for the Stokes and Navier–Stokes Equations , 1984 .

[10]  S. Son,et al.  Unsteady combustion modeling of energetic solids , 1995 .

[11]  Donald Scott Stewart,et al.  Nonsteady burning of periodic sandwich propellants with complete coupling between the solid and gas phases , 2001 .

[12]  Xu Zhou,et al.  Exploring propellant science using a fully-coupled combustion code , 2003 .

[13]  Ananias G. Tomboulides,et al.  A Quasi-Two-Dimensional Benchmark Problem for Low Mach Number Compressible Codes , 1998 .

[14]  J. Buckmaster,et al.  Three-dimensional heterogeneous propellant combustion , 2002 .

[15]  Thomas L. Jackson,et al.  Random Packs and Their Use in Modeling Heterogeneous Solid Propellant Combustion , 2001 .

[16]  L. Massa,et al.  Unsteady combustion modelling of energetic solids, revisited , 2004 .

[17]  Steven A. Orszag,et al.  Numerical Simulation of Low Mach Number Reactive Flows , 1997 .

[18]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[19]  R. Sani,et al.  On pressure boundary conditions for the incompressible Navier‐Stokes equations , 1987 .

[20]  David E. Keyes,et al.  Numerical Solution of Two-Dimensional Axisymmetric Laminar Diffusion Flames , 1986 .

[21]  J. Buckmaster,et al.  Heterogeneous Propellant Combustion , 2002 .

[22]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[23]  John B. Bell,et al.  An Adaptive Mesh Projection Method for Viscous Incompressible Flow , 1997, SIAM J. Sci. Comput..

[24]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[25]  A. Hindmarsh.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

[26]  J. Kan.  A second-order accurate pressure correction scheme for viscous incompressible flow , 1986 .

[27]  H. S. Mukunda,et al.  Sandwich propellant combustion: Modeling and experimental comparison , 2002 .