Troels Dyhr Pedersen.indd - Solid Mechanics
Troels Dyhr Pedersen.indd - Solid Mechanics
Troels Dyhr Pedersen.indd - Solid Mechanics
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
- 73 - -<br />
Chemical reaction mechanism implementation<br />
STAR-CD offers two choices of solvers for version 4.10: The “coupled complex<br />
chemistry” solver and DARS-CFD. While DARS-CFD would be the preferred option,<br />
licenses were not available and could not be acquired in time for this investigation.<br />
The coupled complex chemistry solver allows up the 30 species in a mechanism. These<br />
species are defined as scalars in the software. The reaction mechanism is implemented as<br />
a list of forward reactions. Reverse reactions are calculated by STAR assuming chemical<br />
equilibrium. The reduced mechanism for DME combustion developed earlier includes<br />
only 27 species and is thus applicable to the solver.<br />
Formatting of the input file is critical and must be done according to the rules dictated by<br />
the manual. These are slightly different from the format used by CHEMKIN. It should be<br />
mentioned that it is absolutely critical to check the contents of the cpl.inp01-echo file<br />
created by STAR, which holds the information that STAR-CD uses. It was found that<br />
minor deviations from the required procedure in making the input file resulted in severe<br />
misinterpretations by STAR, such as several orders of magnitude in reaction parameters.<br />
The transport properties are defined in the file tran.dat which is identical to the one used<br />
by CHEMKIN.<br />
12.4.2 Simulation results<br />
Although a number of variations were made to the setup of the test case, not much<br />
difference was found in the transient solution in those cases that successfully resulted in<br />
detonations. The general tendency for development in temperature, pressure and gas<br />
velocity is shown in figures 29-31. Each figure is composed of 10 frames recorded at 5<br />
microsecond intervals, with the first frame at top. The scale for temperature is Kelvin, for<br />
pressure it is Pascal and for gas velocity it is m/s.<br />
To create a disturbance that will ultimately trigger the detonation, a small section at the<br />
end of the volume was subjected to an elevated temperature of 1500 K. This establishes a<br />
flame front (left of the pictures) as well as inducing a small pressure wave that travels<br />
through the volume. This method is widely used to trigger detonation in experiments as<br />
well, typically by igniting the mixture with a spark plug.<br />
While the flame is travelling it is accelerated into the unburned mixture by the expansion<br />
of the gas in the burned fraction behind the flame. This movement compresses the<br />
mixture ahead, which is also starting to react. At some point, typically when that flame<br />
front has moved less than a centimetre through the volume, a minor section of the<br />
mixture auto ignites some distance ahead of the flame. When this happens, the rapid<br />
expansion triggers the development of the detonation.