CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

having trouble using reactingFoam with reactions turned off

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree4Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 16, 2011, 07:13
Default
  #21
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 174
Rep Power: 8
linch is on a distinguished road
third try...
linch is offline   Reply With Quote

Old   June 21, 2011, 20:02
Default
  #22
New Member
 
Scot Johnson
Join Date: Mar 2011
Posts: 25
Rep Power: 6
Scot is on a distinguished road
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
Scot is offline   Reply With Quote

Old   June 29, 2011, 15:44
Default
  #23
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 8
novyno is on a distinguished road
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
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Old   September 7, 2011, 19:15
Default
  #24
New Member
 
Sunny Karnani
Join Date: Apr 2010
Posts: 22
Rep Power: 7
skarnani is on a distinguished road
Quote:
Originally Posted by nakul View Post
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.
skarnani is offline   Reply With Quote

Old   April 6, 2012, 16:23
Default
  #25
New Member
 
Mostafa Moghaddami
Join Date: Oct 2009
Posts: 13
Rep Power: 7
Mostafa is on a distinguished road
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.
Mostafa is offline   Reply With Quote

Old   September 3, 2013, 21:26
Default
  #26
Member
 
赵庆良
Join Date: Aug 2013
Posts: 56
Rep Power: 4
zqlhzx is on a distinguished road
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"?
zqlhzx is offline   Reply With Quote

Old   September 4, 2013, 08:53
Default
  #27
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi

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

-Nakul
nakul is offline   Reply With Quote

Old   September 4, 2013, 09:21
Default
  #28
Member
 
赵庆良
Join Date: Aug 2013
Posts: 56
Rep Power: 4
zqlhzx is on a distinguished road
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!
zqlhzx is offline   Reply With Quote

Old   September 5, 2013, 03:13
Default
  #29
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
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.
nakul is offline   Reply With Quote

Old   September 5, 2013, 03:20
Default
  #30
Member
 
赵庆良
Join Date: Aug 2013
Posts: 56
Rep Power: 4
zqlhzx is on a distinguished road
Thank you,nakul!It is kind of you!
zqlhzx is offline   Reply With Quote

Old   September 8, 2013, 11:36
Default
  #31
Member
 
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 78
Rep Power: 8
gregor is on a distinguished road
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
gregor is offline   Reply With Quote

Old   September 11, 2013, 08:39
Default help regarding reactionFoam solver .
  #32
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
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 :
Code:
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 .
yash.aesi is offline   Reply With Quote

Old   September 11, 2013, 12:34
Default
  #33
Senior Member
 
Adhiraj
Join Date: Sep 2010
Location: Pennsylvania, United States
Posts: 101
Rep Power: 6
adhiraj is on a distinguished road
Nothing. "$FOAM_CASE" should be resolved by the solver at runtime.
adhiraj is offline   Reply With Quote

Old   September 11, 2013, 13:27
Default
  #34
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
greeting adhiraj ,

thanks for replying . sori but can you elaborate a bit more i am not getting what you want to say .
yash.aesi is offline   Reply With Quote

Old   September 11, 2013, 14:48
Default
  #35
Senior Member
 
Adhiraj
Join Date: Sep 2010
Location: Pennsylvania, United States
Posts: 101
Rep Power: 6
adhiraj is on a distinguished road
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.
adhiraj is offline   Reply With Quote

Old   September 11, 2013, 15:16
Default
  #36
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
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
yash.aesi is offline   Reply With Quote

Old   September 12, 2013, 06:30
Default
  #37
Senior Member
 
Adhiraj
Join Date: Sep 2010
Location: Pennsylvania, United States
Posts: 101
Rep Power: 6
adhiraj is on a distinguished road
What error do you get when you retain "$FOAM_CASE/..." ?
Can you post the error message in its entirety?
adhiraj is offline   Reply With Quote

Old   September 12, 2013, 06:47
Default
  #38
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 is offline   Reply With Quote

Old   September 12, 2013, 07:00
Default
  #39
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
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
Attached Files
File Type: zip constant.zip (4.4 KB, 17 views)
yash.aesi is offline   Reply With Quote

Old   September 14, 2013, 12:07
Default
  #40
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 4
yash.aesi is on a distinguished road
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 rhoInitial residual 0, Final residual 0No 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<EquationOfStateout of temperature range 200 -> 5000;  5006.76 
so can anyone guide me where i am doing wrong ??

Thanks ,
Regard
sonu .
yash.aesi is offline   Reply With Quote

Reply

Tags
hseqn, reactingfoam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
ReactingFoam with surface reactions robert_mornhinweg OpenFOAM Running, Solving & CFD 8 April 21, 2015 06:20
No reactions in reactingFoam 2.1 OMN OpenFOAM Running, Solving & CFD 16 April 7, 2015 12:14
PaSR + infinite reaction rate in reactingFoam --> no reactions occurring tatu OpenFOAM Running, Solving & CFD 1 November 21, 2012 06:02
reactingFoam wedge handling wrong U dhondupant OpenFOAM Bugs 1 December 9, 2010 08:34
ReactingFoam without reactions lasb OpenFOAM Running, Solving & CFD 5 June 10, 2008 08:50


All times are GMT -4. The time now is 20:56.