zdz0309 / QSSAchemistrySolver

The quasi-steady-state approximation chemistry solver for OpenFOAM 4 and 7

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

QSSAchemistrySolver

The quasi-steady-state approximation chemistry solver for OpenFOAM 4 and 7

This repository is migrated from adhiraj-dasgupta/unsupportedContribOF23x. For OpenFOAM-2.3, please refer to that repository.

compile the code and test it

  • chose a branch to download. This repository works for OpenFOAM-4.x (master branch) and OpenFOAM-7 (of7 branch).
  • come into the CODE folder and wmake to compile
    git clone git@github.com:ZmengXu/QSSAchemistrySolver.git
    cd CODE
    wmake
  • come into the CASE folder and run the demo case
    cd CASE
    chmod 777 Allrun Allclean
    ./Allrun

run QSS in your own case

Here is the instructions for your own case:

  • In $FOAM_CASE/system/controlDict, compile to get a QSS mechanism library, in shell typing the following command to compile mechanism
    gfortran -shared -fPIC -o libMechnism19sCH4.so ckwyp.f
  • In $FOAM_CASE/system/controlDict
    libs
    (
        "./chemkin/libMechnism19sCH4.so"
        "libQSSchemistryModel.so"
    );
  • In $FOAM_CASE/constant/chemistryProperties
    chemistryType
    {
        chemistrySolver QSS;
        chemistryThermo psi;
    }
    QSSCoeffs
    {
        solver          seulex;
        absTol          1e-12;
        relTol          0.1;
    }

Here's the temperature trace results for 19 species QSS and 30 species mechanism.

TemperatureTrace

About

The quasi-steady-state approximation chemistry solver for OpenFOAM 4 and 7

License:MIT License


Languages

Language:Fortran 75.8%Language:C 11.7%Language:C++ 10.7%Language:Shell 1.9%