Talk:Tut reactingFoam firstTutorial

From OpenFOAMWiki

1 adaption to OpenFOAM-1.6

It is possible to use the tutorial in OpenFOAM-1.6 with some minor changes:

  • a "g"-file is needed in the /constant directory, which can be used from the new aachenBomb-tutorial
  • file "constant/chemistryProperties"
    • chemistrySolver should be ode instead of ODE
    • there should be a odeCoeffs-dictionary instead of ODECoeffs
    • added the line psyChemistryModel ODEChemistryModel<gasThermoPhysics>; after the header
  • file "constant/thermophysicalProperties"
    • replace line thermoType hMixtureThermo<reactingMixture>; by thermoType hPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;
    • check for correct places of CHEMKINFile and CHEMKINThermoFile, it works with CHEMKINFile "../chemkin/chem.inp"; and CHEMKINThermoFile "../chemkin/therm.dat";
  • file "0/T"
    • 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


--Aortwein 14:01, 12 May 2010 (UTC)

Now: adaptation to OpenFOAM-1.7.1 --Aortwein 13:23, 14 December 2010 (CET)

2 older question

Thanks for the tutorial, however, I am getting the following error in the third time step, and have not been able to find the problem:

--> FOAM FATAL ERROR : attempt to use janafThermo<equationOfState> out of temperature range 200 -> 5000; T = 5000.92#0

Any suggestions or a fix would be very appreciated. Thanks again. ~E

  • a possibility could be to change the oxygen in the lowerInlet (in the the 0/O2 file) from 1.0 to (say) 0.9 or less and the nitrogen in the lowerInlet ( 0/N2) from 0.0 to 0.1 to cool down the reaction
  • note: my suggestion is not validated, since I did not have that problem yet, but the problem is also discussed here: at cfd-online.com --Aortwein 15:29, 4 December 2009 (UTC)
  • as far as I can see, the main problem lies in the "T" file in the "0"-Directory. Here, the boundary condition for "outlet" is set to "type fixedValue" with a value of 800. That should be changed to "type zeroGradient", afterwards it should work. If that change is not done, obviously that boundary condition is making problems as soon as a temperature higher than 800 K reaches the outlet. --Aortwein 10:31, 16 December 2009 (UTC)
  • One suggestion is to increase upper limit of the temperature in thermo.dat file. I have fixed this problem using this technique.

You cannot post new threads to this discussion page because it has been protected from new threads.

There are no threads on this page yet.