PDE Simulation

กก3/6/2003

1.      Flow control equations:
Axis-symmetric Euler equations + one species equation.

2.      Numerical method: CESE method

3.      Computational Domain: [0, 56.9 ft] x [-30, 30 ft]
(Because of symmetry, only up half of the domain is calculated)

4.      Mesh: 6,400,000 cells
 
(Inside and near thrust tube, the mesh is uniform; the others are stretching)

5.      Boundary conditions:
-- On solid wall and symmetric line, reflective boundary condition is used.
-- On outlet and open surface, non-reflective boundary condition is used.

6.      Flow conditions:
Chamber temperature T = 510 R;   Chamber pressure P = 15.0 Psia;
Equivalence ratio F = 1.3;  Freq. = 60 Hz
C-J detonation velocity VC-J = 6200 ft/sec.;  Pressure jump P2/P1 = 16.667.
Diameter of tube = 5 inches, and its length =26.138inches.

7.   Pressure contour:
At time 40.15 ms
At time 99.36 ms
M
ovie (up to 6 cycles, 100ms): (17 MB GIF format)

8.      Domain and checking  points distribution:

9.      Pressure history (up to 6 cycles, 100ms) :

Pressure history on the wall at the dead end of the thrust tube.

5 ft 10 ft 15 ft 20 ft 25 ft 30 ft 35 ft 40 ft 45 ft 50 ft
0o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
3o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
6o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
9o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
12o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
15o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
18o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
21o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
24o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
27o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
30o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
33o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
36o P-t P-t P-t P-t P-t P-t P-t P-t P-t P-t
39o P-t P-t P-t P-t P-t P-t P-t P-t P-t --
42o P-t P-t P-t P-t P-t P-t P-t P-t P-t --
45o P-t P-t P-t P-t P-t P-t P-t P-t -- --