Propellant grains burning in a simple internal ballistics demonstration model.
*UNIT_SYSTEM
SI
*PARAMETER
rho = 2200.0, "powder density"
Cv = 700.0, "heat capacity"
gamma = 1.4, "Cp/Cv ratio"
e0 = 1.2e6, "combustion energy per unit mass"
cvol = 0.3, "initial co-volume"
a = 1.0e-5, "burn rate parameter 1"
n = 0.6, "burn rate parameter 2"
c = 0.0001, "burn rate parameter 3"
patm = 1.0e5, "burn rate parameter 4"
Ti = 500.0, "ignition temperature"
pi = 1.0e6, "ignition pressure"
xi = 0.0, "ignition point X"
yi = 0.0, "ignition point Y"
zi = 0.0, "ignition point Z"
ti = 0.0, "ignition time"
Ri = 0.02, "ignition radius"
Np = 800000, "total number of particles"
#
# --- Propellant material properties ---
#
*MAT_POWDER_BURN
2, [%rho]
[%Cv], [%gamma], [%e0], [%cvol], 123, [%Ti], [%pi]
*FUNCTION
"burn rate"
123
%a*max(%patm,p)^%n + %c
#
# --- Ignition location ---
#
*POWDER_BURN_IGNITE
"ignition location"
1
[%xi], [%yi], [%zi], [%ti], [%Ri]
#
# --- Define number of particles in PART command ---
#
*PART
"propellant"
2, 2, 0, 0, 0, 3.0e-8, 0, 0, [%Np]
#
# --- Define region in space where the discrete particles are active ---
#
*PARTICLE_DOMAIN
ALL, 0, 0
-0.06, -0.06, -0.01, 0.06, 0.06, 1.0
*END