28.08.2013 Views

Troels Dyhr Pedersen.indd - Solid Mechanics

Troels Dyhr Pedersen.indd - Solid Mechanics

Troels Dyhr Pedersen.indd - Solid Mechanics

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!