[chenyr10]

Hello:

I want to use reduced mechanism(with subroutine called ckwyp.f) to determine the reaction source term in reactingFoam.

Is there any openfoam extensions or simple way to use this reduced mechanism instead of calculating from chemistry solver?

Thanks!

[Slanth]

Hi, have you find a way to solve this problem, I am also working on that now. I think maybe I should replace the original omega function (in ODEchemistrymode) with the fortran soubrountine?

[arvindpj]

You can check out ISAT-CK7 + Cantera + OpenFOAM and also then new TDAC models.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      chemistryProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

chemistryType
{
    chemistrySolver   ISAT; // ISAT
    chemistryThermo   psi;
}

chemistry           on;

ISATCoeffs
{
    saveISATtree on;
    externalCKWYP off;
    constantPressure off;
    ISATABabsErr 1e-5;
    ISATABrelErr 0;
    DDASACabsErr 1e-8;
    DDASACrelErr 1e-9;
    ISATABsizeMB 500;    
}

initialChemicalTimeStep 1e-7;

EulerImplicitCoeffs
{
    cTauChem        1;
    equilibriumRateLimiter off;
}

odeCoeffs
{
    solver          SIBS;
    absTol          1e-12;
    relTol          0.01;
}

// ************************************************************************* //

Cheers,