third try...
|
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 |
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 |
Quote:
|
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. |
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"? |
Hi
'chemistrySh' is the conventional 'chemistry source term' of your energy equation. you can see its definition in createFields.H -Nakul |
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! |
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. |
Thank you,nakul!It is kind of you!
|
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 |
help regarding reactionFoam solver .
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 PHP Code:
Regards , sonu . |
Nothing. "$FOAM_CASE" should be resolved by the solver at runtime.
|
greeting adhiraj ,
thanks for replying . sori but can you elaborate a bit more i am not getting what you want to say . |
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.
|
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 |
What error do you get when you retain "$FOAM_CASE/..." ?
Can you post the error message in its entirety? |
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // |
1 Attachment(s)
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 |
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:
Thanks , Regard sonu . |
All times are GMT -4. The time now is 04:30. |