View Issue Details

IDProjectCategoryView StatusLast Update
0003951OpenFOAMContributionpublic2023-01-25 15:59
Reporterevamaria Assigned Tohenry  
PrioritylowSeverityminorReproducibilityalways
Status closedResolutionsuspended 
OSUbuntuOS Version22.04 
Product Version9 
Summary0003951: Contribution for PLOG Chemistry
DescriptionAs mentioned by other reports, the new reaction format from chemkin with pressure dependent reaction rates (PLOG) is not usable in OpenFOAM.

Therefore, I modified the chemkinReader to be able to read those reaction types and added a pressureDependentArrheniusReactionRate similar to the PLOG reaction in chemkin. My current implementation is in OpenFOAM-9, but if you consider to add this feature, I would be happy to contribute it for OpenFOAM-dev.
TagsNo tags attached.

Activities

evamaria

2023-01-23 13:47

reporter  

chemkinIItoFoam.zip (322,975 bytes)
makePressureDependentReactions.C (1,549 bytes)   
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Copyright (C) 2011-2021 OpenFOAM Foundation
     \\/     M anipulation  |
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.

    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.

\*---------------------------------------------------------------------------*/

#include "makeReaction.H"

#include "pressureDependentArrheniusReactionRate.H"

#include "forGases.H"

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

namespace Foam
{
    forCoeffGases(makeIRReactions, pressureDependentArrheniusReactionRate);
}

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

henry

2023-01-25 10:42

manager   ~0012937

Your implementation of the PLOG reaction rate does not appear to be correct, see e.g.

https://cantera.org/documentation/docs-2.4/doxygen/html/dd/d89/classCantera_1_1Plog.html

could you provide the reference you used for the implementation and the validation case you used?

evamaria

2023-01-25 11:18

reporter   ~0012938

No, you are correct, there is something not correct in my PLOG reaction rate. I will try to fix it and provide a reference case with ignition delay times. But I guess, in the mean time you can close this issue.

Issue History

Date Modified Username Field Change
2023-01-23 13:47 evamaria New Issue
2023-01-23 13:47 evamaria File Added: chemkinIItoFoam.zip
2023-01-23 13:47 evamaria File Added: makePressureDependentReactions.C
2023-01-23 13:47 evamaria File Added: pressureDependentArrheniusReactionRate.zip
2023-01-25 10:42 henry Note Added: 0012937
2023-01-25 11:18 evamaria Note Added: 0012938
2023-01-25 15:59 henry Assigned To => henry
2023-01-25 15:59 henry Status new => closed
2023-01-25 15:59 henry Resolution open => suspended