Есть также второй набор файлов, используемых чтобы контролировать химию в решателе.
Вот они:
· chemistryProperties (химические свойства)
· combustionProperties (свойства горения)
· thermophysicalProperties (теплофизические (термодинамические?) свойства)
В файле chemistryProperties химия и турбулентная реакция (turbulentReaction) должны быть включены (turned on), как химический решатель использует ODE. Оставьте настройки в combustionProperties как есть. В thermophysicalPropertiesт убедитесь, что пути к вашим файлам хемкина например chem.inp и therm.dat верны.
Эффективная диффузия
Эффективная диффузия принемается как эффективная вязкость что подразумевает турболентное число Шмидта один. Ламинарная вязкость в свою очередь рассчитывается с использованием закона Сатерленда, где константы жестко закодированы в фактический ридер хемкин как
и 
Настройка системы и запуск проекта
В controlDict используйте следующие настройки:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application reactingFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 1;
deltaT 1e-02;
writeControl adjustableRunTime;
writeInterval 1e-01;
purgeWrite 0;
writeFormat binary;
writePrecision 6;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
adjustTimeStep yes;
maxCo 0.1;
runTimeModifiable yes;
// ************************************************************************* //
Заускайте проект с помощью
blockMesh. reactingFoamTutorial
reactingFoam. reactingFoamTutorial
3 Adaption to OpenFOAM-1.7
It is possible to use the tutorial in OpenFOAM-1.7.1 with some minor changes:
· a "g"-file is needed in the /constant directory, which can be used from the new combustion/dieselFoam/aachenBomb-tutorial
· file "constant/chemistryProperties"
o put the header (from "version" to "object") into FoamFile{... }
o chemistrySolver should be ode instead of ODE
o there should be a odeCoeffs-dictionary instead of ODECoeffs
o added the line psyChemistryModel ODEChemistryModel<gasThermoPhysics>; after the header
· file "constant/thermophysicalProperties"
o replace line thermoType hMixtureThermo<reactingMixture>; by thermoType hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;
|
|
o check for correct places of CHEMKINFile and CHEMKINThermoFile, it works with CHEMKINFile "../chemkin/chem.inp"; and CHEMKINThermoFile "../chemkin/therm.dat";
· file "cons*-tant/turbulenceProperties"
o replace the solver section with one from a similar file from the new aachenBomb-tutorial
· file "0/T"
o in boundary condition "outlet", replace type fixedValue; value uniform 800; with type zeroGradient;
· a "constant/RASProperties"-file could be needed, examples can be found e.g. in the new aachenBomb-tutorial
· file "system/controlDict"
o set startFrom to startTime instead of latestTime