Numerical Solution of Two-Dimensional Eulerian Reactive Flow

Numerical Solution of Two-Dimensional Eulerian Reactive Flow 

The Eulerian equations of motion are more useful for numerical solution of highly distorted fluid flow than are Lagrangian equations of motion. Multicomponent Eulerian calculations require equations of state for mixed cells and methods for moving mass and its associated state values into and out of mixed cells. These complications are avoided by Lagrangian calculations. Harlow s particle-in-cell (PIC) method uses particles for the mass movement. The first reactive Eulerian hydrodynamic code EIC (Explosive-in-cell) used the PIC method, and it is described in reference 2. The discrete nature of the mass movement introduced pressure and temperature variations from cycle to cycle of the calculation that were unacceptable for many reactive fluid dynamic problems. A one-component continuous mass transport Eulerian code developed in 1966 proved useful for solving many one-component problems of interest in reactive fluid dynamics. The need for a multicomponent Eulerian code resulted in a second 2DE code, described in reference 4. Elastic-plastic flow and real viscosity were added in 1976. The technique was extended to three dimensions in the 1970 s and the resulting 3DE code is described in Appendix D. 

FORTRAN and executable versions of the 2DE code are on the CD-ROM that operate under the Windows XP or VISTA operating systems. The codes use ABSOFT FORTRAN with PLPLOT GRAPHICS. 

Bulk or Artificial Viscosity XX Viscosity Deviator XZ Viscosity Deviator Gas Constant ZZ Viscosity Deviator XX Elastic Stress Deviator XZ Elastic Stress Deviator ZZ Elastic Stress Deviator Temperature K 

Particle Velocity in X or R direction (i) Particle Velocity in Z direction (j) 

---

The Eulerian (bottom-up) approach is to start with the convective-diffusion equation and through Reynolds averaging, obtain time-smoothed transport equations that describe micromixing effectively. Several schemes have been proposed to close the two terms in the time-smoothed equations, namely, scalar turbulent flux in reactive mixing, and the mean reaction rate (Bourne and Toor, 1977 Brodkey and Lewalle, 1985 Dutta and Tarbell, 1989 Fox, 1992 Li and Toor, 1986). However, numerical solution of the three-dimensional transport equations for reacting flows using CFD codes are prohibitive in terms of the numerical effort required, especially for the case of multiple reactions with... 

---