[linch]

third try...

[Scot]

Hi Valerio-
What's the latest on the multispecies diffusion libraries you were working on (Fick model, Maxwell-Stefan, etc)? Anything released yet? This would be helpful! With so many solvers in openFoam, it seems there's not enough attention given to chemistry and more is always welcome. Thank you and best wishes,
Scot

[novyno]

Hi Scot,

have a look on http://www.extend-project.de/user-gr...ewgroup/groups

Send me an e-mail and I can send you the code.

Regards

Valerio

[skarnani]

Originally Posted by nakul

Hi Valerio,

Correct me if I am wrong but I think that if assumption of Le=1 is there then we should have alpha instead of mu in the species equation.

In hs-equ. also for Le=1 they have used alphaEff().

According to me muEff() is there to account for turbulent diffusion. Please clarify me if I am wrong!!!

I agree. I'm a little confused by the unit Lewis number note. By equating alpha and mu, doesn't this imply that Pr = 1. Laminar simulations run with this case resulted in artificially high flame temperatures (~2100K). I was able to get more reasonable flame temperatures by multiplying the mu*1.42, which would give Pr = 0.7. Any thoughts or similar experiences would be helpful.

[Mostafa]

Dear All

I am going to run a case includes some liquid species that react with each other. I am going to use the reactingFoam as a solver (OF 2.0.1).
All species have constant proprties and I want to use the following thermo type:
thermoType hsPsiThermo<reactingMixture<constTransport<specieT hermo<eConstThermo<icoPolynomial>>>>>>

but when I set it in the constat/thermophysicalProperties file and run the case, I get this error:

------------------------------------------------------------------------------------
--> FOAM FATAL ERROR:
Inconsistent thermo package selected:

hsPsiThermo<reactingMixture<constTransport<specieT hermo<eConstThermo<icoPolynomial>>>>>>

Please select a thermo package based on gasThermoPhysics. Valid options include:

3
(
hsPsiMixtureThermo<singleStepReactingMixture<gasTh ermoPhysics>>
hsPsiMixtureThermo<multiComponentMixture<gasThermo Physics>>
hsPsiMixtureThermo<reactingMixture<gasThermoPhysic s>>
)



From function autoPtr<hsCombustionThermo> hsCombustionThermo::NewType(const fvMesh&, const word&)
in file combustionThermo/hsCombustionThermo/hsCombustionThermoNew.C at line 116.

FOAM exiting
-------------------------------------------------------------------------------------

What should I do to use the reactingFoam for constant propertie species?

Thanks in advance for your help.

[zqlhzx]

Hi nakul:
would you tell me what "chemistrySh" represents in hEqn.H in reactingFoam?Please!Is there any theroy can corresponding to the term "chemistrySh"?

[nakul]

Hi

'chemistrySh' is the conventional 'chemistry source term' of your energy equation. you can see its definition in createFields.H

-Nakul

[zqlhzx]

thanks,nakul,I find the define of chemistrySh in the chemisty.H,the following is codes:

if (chemistry.chemistry())
{
Info<< "Solving chemistry" << endl;
chemistry.solve
(
runTime.value() - runTime.deltaTValue(),

runTime.deltaTValue()


);

// turbulent time scale

if (turbulentReaction)
{
tmp<volScalarField> tepsilon(turbulence->epsilon());
const volScalarField& epsilon = tepsilon();
tmp<volScalarField> tmuEff(turbulence->muEff());
const volScalarField& muEff = tmuEff();
tmp<volScalarField> ttc(chemistry.tc());
const volScalarField& tc = ttc();

forAll(epsilon, i)
{
if (epsilon[i] > 0)

{
// Chalmers PaSR model
scalar tk = Cmix.value()*Foam::sqrt(muEff[i]/rho[i]/epsilon[i]);
kappa[i] =
(runTime.deltaTValue() + tc[i])
/(runTime.deltaTValue() + tc[i] + tk);
}
else
{
// Return to laminar combustion

kappa[i] = 1.0;
}
}
}
else
{
kappa = 1.0;

}
chemistrySh = kappa*chemistry.Sh()();
}
But I really don not know what formula expression can corresponding to the term "chemistrySh",would you tell me its formula expression about the term?Or if you have some paper about the term ,could you let me take look?please!

[nakul]

Hi

Read this paper - http://powerlab.fsb.hr/ped/kturbo/Op...olmPhD2008.pdf

It has details on how reaction source term is modeled in OF.

As I said in previous post, chemistrySH is chemical source term multiplied by a factor 'kappa'. kappa takes into account the effect of turbulence. It is so because reactingFoam solves chemistry based upon turbulence of the flow.

The philosophy is similar to Eddy Dissipation but formulation is different.

[zqlhzx]

Thank you,nakul!It is kind of you!

[gregor]

Hi guys,

would you mind running my enthalpy diffusion test case (https://github.com/greole/ValidationCases and http://www.cfd-online.com/Forums/openfoam/123227-enthalpy-diffusion-testcase.htm) with your reactingFoam version to see if the mean enthalpy ist conserved ?

Help would be much apprechiated

[yash.aesi]

greeting oll ,

i am relatively new to OF. nowadays i am trying to simulate my case of different species with the reactingFoam solver and m using openFoam 2.2-x . fist i solved the tutorial of counterFlow Flame after that nw i tried to run my case . but i am not able to understand what to write in place of bold marked text in thermophysical properties file shown below :

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

chemistryReader foamChemistryReader;

foamChemistryFile "$FOAM_CASE/constant/reactions";

foamChemistryThermoFile "$FOAM_CASE/constant/thermo.compressibleGas";

reactions and thermo.compressiblegas files are there in constant folder . but if i change CASE to my problem name then it will say :

PHP Code:

--> FOAM FATAL ERROR: 
Unknown variable name 'FOAM_counterFlow2Dmy'

    From function stringOps::inplaceExpand(string&, const bool)
    in file primitives/strings/stringOps/stringOps.C at line 749.

FOAM exiting 


thanks in advance ,
Regards ,
sonu .

[adhiraj]

Nothing. "$FOAM_CASE" should be resolved by the solver at runtime.

[yash.aesi]

greeting adhiraj ,

thanks for replying . sori but can you elaborate a bit more i am not getting what you want to say .

[adhiraj]

What I meant was that if you have the phrase "$FOAM_CASE" in the file, during running, the solver will be able to resolve the path. You do not need to explicitly put in your case directory name.

[yash.aesi]

thanks adhiraj for quick replying .
but if i keep it as it is means "$FOAM_CASE" then it will end in big error of print stack and if i change it like i said above then again end up with diffrent error .

one more doubt i am using OF 2.2-X. so there in constant folder we are having reaction and thermo. file so do we need to separately keep chemkin folder alone with our case files .


thanks again ,
Regards ,
sonu

[adhiraj]

What error do you get when you retain "$FOAM_CASE/..." ?
Can you post the error message in its entirety?

[yash.aesi]

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Creating reaction model

Selecting combustion model PaSR<psiChemistryCombustion>
Selecting chemistry type 
{
    chemistrySolver ode;
    chemistryThermo psi;
}

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

Selecting chemistryReader foamChemistryReader
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::fvPatchField<double>::operator/=(Foam::fvPatchField<double> const&) at ??:?
#4  Foam::FieldField<Foam::fvPatchField, double>::operator/=(Foam::FieldField<Foam::fvPatchField, double> const&) at ??:?
#5  Foam::multiComponentMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::correctMassFractions() at ??:?
#6  Foam::multiComponentMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::multiComponentMixture(Foam::dictionary const&, Foam::List<Foam::word> const&, Foam::HashPtrTable<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >, Foam::word, Foam::string::hash> const&, Foam::fvMesh const&) at ??:?
#7  Foam::reactingMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::reactingMixture(Foam::dictionary const&, Foam::fvMesh const&) at ??:?
#8  Foam::SpecieMixture<Foam::reactingMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::SpecieMixture(Foam::dictionary const&, Foam::fvMesh const&) at ??:?
#9  Foam::heThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#10  Foam::hePsiThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::hePsiThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#11  Foam::psiReactionThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#12  Foam::autoPtr<Foam::psiReactionThermo> Foam::basicThermo::New<Foam::psiReactionThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#13  Foam::psiReactionThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#14  Foam::psiChemistryModel::psiChemistryModel(Foam::fvMesh const&) at ??:?
#15  Foam::chemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > >::chemistryModel(Foam::fvMesh const&) at ??:?
#16  Foam::ode<Foam::chemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::ode(Foam::fvMesh const&) at ??:?
#17  Foam::psiChemistryModel::addfvMeshConstructorToTable<Foam::ode<Foam::chemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh const&) at ??:?
#18  Foam::autoPtr<Foam::psiChemistryModel> Foam::basicChemistryModel::New<Foam::psiChemistryModel>(Foam::fvMesh const&) at ??:?
#19  Foam::psiChemistryModel::New(Foam::fvMesh const&) at ??:?
#20  Foam::combustionModels::psiChemistryCombustion::psiChemistryCombustion(Foam::word const&, Foam::fvMesh const&) at ??:?
#21  Foam::combustionModels::PaSR<Foam::combustionModels::psiChemistryCombustion>::PaSR(Foam::word const&, Foam::fvMesh const&) at ??:?
#22  Foam::combustionModels::psiCombustionModel::adddictionaryConstructorToTable<Foam::combustionModels::PaSR<Foam::combustionModels::psiChemistryCombustion> >::New(Foam::word const&, Foam::fvMesh const&) at ??:?
#23  Foam::combustionModels::psiCombustionModel::New(Foam::fvMesh const&) at ??:?
#24  
 at ??:?
#25  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#26  
 at ??:?
Floating point exception (core dumped)

[yash.aesi]

Greeting adhiraj ,
The above posted error comes if i retain it as it is and i dont have any idea what kind of error it is .
if i run the tutorial of counterflowflame2D there its working fine . i am attaching my constant folder .
please have a look and guide me .
thank you .

Regards ,
sonu

[yash.aesi]

greeting oll ,

Above error is solved its just bcz of the wrong BC's , nw my case is running
but there is no combustion ( combustion is on ) . The BC's for temp for fuel inlet is 300 k and air inlet is 500 k . but in this case there combustion is not happening why ?

also if i try to change the BC's to some higher temperature like 600 k then following error is coming :

PHP Code:

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
--> FOAM Warning : 
    From function janafThermo<EquationOfState>::limit(const scalar T) const
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000;  T = 5006.76 


so can anyone guide me where i am doing wrong ??

Thanks ,
Regard
sonu .