Study of the EEqn.H in rhoPimpleDyMFoam.

Horacio Aguerre · May 27, 2013, 09:03

Hi,

I'm working with the solver "rhoPimpleDyMFoam" from the OpenFOAM 2-2-x version. In specific i was trying to replicate the polytropic process in an adiabatic compression.
The case i studied is a pistone-cylinder sistem where air is compressed by the pistone patch which is a moving wall. The pistone velocity is adjusted to 0,1(m/s) to aproximate a reversible process.

When i choose energy "he" as internal energy i could reproduce the polytropic process in an adiabatic compression.

However when i choose energy "he" as entalphy I couldn't achieve the theoretical expression $pV^{1,4}=constant$ for air. Instead the system looks as it would be non adiabatic or like the conversion between work and internal energy was not made properly.

Therefore i started to study the EEqn.H.

I compare the theoretical balance of total energy in terms of entalphy with what is written in the EEqn.H.

Balance of a tensorial property in an arbritary moving domain.

$\dfrac{d}{dt} \int_{V}\rho \phi dV + \int_{S}\rho \phi\left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}=-\int_{S}\rho \mathbf{Q_{\phi}} \cdot \mathbf{ds}+\int_{V}S_{\phi}dV$

Balance of mechanical energy.

$\phi \equiv \frac{\mathbf{v} \cdot \mathbf{v}}{2} = K$

$\dfrac{d}{dt} \int_{V}\rho K dV + \int_{S}\rho K \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S}\left(p\mathbf{v}+ \boldsymbol{\tau} \cdot \mathbf{v}\right) \cdot \mathbf{ds} + \int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV+ \int_{V}\rho\left(\mathbf{v} \cdot \mathbf{g}\right) dV$

Balance of internal energy.

$\phi \equiv \widehat{u}$

$\dfrac{d}{dt} \int_{V}\rho \widehat{u} dV + \int_{S}\rho \widehat{u} \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV$

The balance of internal energy can be expressed in terms of entalphy.

$\widehat{h}=\widehat{u}+\frac{p}{\rho}\implies \widehat{u}=\widehat{h}-\frac{p}{\rho}$

Replacing in the balance of internal energy.

$\dfrac{d}{dt} \int_{V}\rho \left(\widehat{h}-\frac{p}{\rho}\right) dV + \int_{S}\rho \left(\widehat{h}-\frac{p}{\rho}\right) \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV$

$\dfrac{d}{dt} \int_{V}\rho \widehat{h} dV+ \int_{S}\rho \widehat{h} \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}-\dfrac{d}{dt} \int_{V}p dV - \int_{S} p \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV$

Balance of total energy.

The sum of the balances of mechanical and internal energy results in the balance of total energy.

$\dfrac{d}{dt} \int_{V}\rho \left(\widehat{h}+ K \right) dV + \int_{S}\rho \left(\widehat{h}+ K\right) \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds} -\dfrac{d}{dt} \int_{V}p dV - \int_{S} p \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{S}\left(p\mathbf{v}+ \boldsymbol{\tau} \cdot \mathbf{v}\right) \cdot \mathbf{ds} + \int_{V}\rho\left(\mathbf{v} \cdot \mathbf{g}\right) dV$

If the terms of work of viscous forces and the gravity are neglected, the balance is:

$\dfrac{d}{dt} \int_{V}\rho \left(\widehat{h}+ K \right) dV + \int_{S}\rho \left(\widehat{h}+ K\right) \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds} -\dfrac{d}{dt} \int_{V}p dV + \int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ $= -\int_{S} \mathbf{q} \cdot \mathbf{ds}$

rhoPimpleDyMFoam: EEqn.H

The balance of total energy in the EEqn.H file with

considered as entalphy $\widehat{h}$ is:

Code:

fvm::ddt(rho, he) + fvm::div(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K)
- dpdt - fvm::laplacian(turbulence ->alphaEff(), he) == fvOptions(rho, he)

Where, $\phi=\rho_{f} \left(\mathbf{v}-\mathbf{v_{b}}\right)_{f} \cdot \mathbf{S_{f}}$ , is relative to the mesh.

If fvOptions(rho, he) are considered as zero, the term $\int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ from the final balance is no considered in EEqn.H.

Then i try adding in the EEqn.H the therm $\int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ .

Code:

fvm::ddt(rho, he) + fvm::div(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K)
- dpdt + fvc::div(fvc::meshPhi(U),p) - fvm::laplacian(turbulence ->alphaEff(), he) == fvOptions(rho, he)

The result was correct. As in the first case (he=internal energy) the adiabatic compression was reproduced with the relation $pV^{1,4}=constant$ .

The question is,

should the term $\int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ be introduced in any way in fvOptions(rho, he) or it must be explicit in the EEqn.H?

Thanks,

Horacio

santiagomarquezd · June 12, 2013, 23:15

Any clues?. It seems we have to upload the question to the bug tracker.

Regards.

santiagomarquezd · July 8, 2013, 19:43

http://www.openfoam.org/mantisbt/view.php?id=909

ALU · August 20, 2013, 09:33

Hi,

I'm trying to simulate exactly the same case, but I get very unphysical solutions even using the internal energy formulation. I also tried your proposal concerning the EEqn.H file without any success.

You said you were able to reproduce the adiabatic compression, so I would appreciate it if you had a look at my settings or if you shared me your case.

Horacio Aguerre · September 4, 2013, 13:37

Quote:

Originally Posted by ALU

Hi,

I'm trying to simulate exactly the same case, but I get very unphysical solutions even using the internal energy formulation. I also tried your proposal concerning the EEqn.H file without any success.

You said you were able to reproduce the adiabatic compression, so I would appreciate it if you had a look at my settings or if you shared me your case.

Hi ALU,

I tried your case without modifying it and i get bad results too.

I made the following modifications,

File: /0$ U

If you want to reproduce an expansion, you should set for the "top" velocity the same velocity of the "patch", in this case, 0.1 m/s for the third coordinate.

Remember that the type "fixedValue" indicates absolute velocities. If you want to set relative to domain velocities you can use type "movingWallVelocity".

File: /system$ fvSolution

I added in the subdictionary PIMPLE the flag "correctPhi" equal to "no".

The correctPhi equation makes the linear solver to diverge.

For compressible fluids, i think that this equation is no necessary as the non conservative effects due to the mesh change in the flux are absorbed by the density equation.

I run the case with a time step of 0.001s from time 0 to time 1s and the adiabatic expansion was obtained successfully.

Quote:

Originally Posted by ALU

I also tried your proposal concerning the EEqn.H file without any success

The modification of the EEqn.H is necessary when you set, in the dictionary "thermophysicalProperties" (in thermoType{}), the energy to sensibleEnthalpy, but we have here sensibleInternalEnergy.

I hope this can help you,

Regards,

Horacio

ALU · September 5, 2013, 03:20

Hi Horacio Aguerre,

thank you very much for the correction of my case and the detailed explanation, this really helps a lot!

Regards

zqlhzx · September 11, 2013, 01:32

Hi ,Horacio Aguerre
could you see my problem and tell some thing?I think you can answer my question.I moves fireFoam of 2.0.1 to 2.2.1,then rewrite hEqn according to the EEqn.H in fireFoam of 2.2.1 and other necessary codes.But ,when I wmake my rewrited fireFoam,it gives me errors:undefine div(phi,h) in fvscheme.
The following is codes of EEqn.H:
{
volScalarField& he = thermo.he();
fvScalarMatrix EEqn
(
fvm::ddt(rho, he) + fvm::div(phi, he)
+ fvc::ddt(rho, K) + fvc::div(phi, K)
+ (
he.name() == "e"
? fvc::div
(
fvc::absolute(phi/fvc::interpolate(rho), U),
p,
"div(phiv,p)"
)
: -DpDt
)
- fvm::laplacian(turbulence->alphaEff(), he)
==
fvOptions(rho, he)
);
EEqn.relax();
fvOptions.constrain(EEqn);
EEqn.solve();
fvOptions.correct(he);
thermo.correct();
Info<< "T min/max = " << min(T).value() << ", " << max(T).value() << endl;
}
There is div(phi,he) insteading of div(phi,h),So It should be no need to define (phi,h) in fvshceme,but define div(phi,he).However,I see the fvsheme in smallpoolfire2D of 2.2.1,there indeed is div(phi,h).In fact,there is the term of "div(phi,he)" in EEqn.H,why define div(phi,h) in fvsheme?I think I have some thing not understand.!

Horacio Aguerre · September 11, 2013, 10:11

Quote:

Originally Posted by zqlhzx

Hi ,Horacio Aguerre
could you see my problem and tell some thing?I think you can answer my question.I moves fireFoam of 2.0.1 to 2.2.1,then rewrite hEqn according to the EEqn.H in fireFoam of 2.2.1 and other necessary codes.But ,when I wmake my rewrited fireFoam,it gives me errors:undefine div(phi,h) in fvscheme.
The following is codes of EEqn.H:
{
volScalarField& he = thermo.he();
fvScalarMatrix EEqn
(
fvm::ddt(rho, he) + fvm::div(phi, he)
+ fvc::ddt(rho, K) + fvc::div(phi, K)
+ (
he.name() == "e"
? fvc::div
(
fvc::absolute(phi/fvc::interpolate(rho), U),
p,
"div(phiv,p)"
)
: -DpDt
)
- fvm::laplacian(turbulence->alphaEff(), he)
==
fvOptions(rho, he)
);
EEqn.relax();
fvOptions.constrain(EEqn);
EEqn.solve();
fvOptions.correct(he);
thermo.correct();
Info<< "T min/max = " << min(T).value() << ", " << max(T).value() << endl;
}
There is div(phi,he) insteading of div(phi,h),So It should be no need to define (phi,h) in fvshceme,but define div(phi,he).However,I see the fvsheme in smallpoolfire2D of 2.2.1,there indeed is div(phi,h).In fact,there is the term of "div(phi,he)" in EEqn.H,why define div(phi,h) in fvsheme?I think I have some thing not understand.!

Hi zqlhzx,

The EEqn is able to solve the energy balance in terms of enthalpy ("h") or internal energy ("e"). Both balances have common expresions. For this expresions it is used the generic energy variable "he". However the expresions differ. Then appears the conditional statement that selects the appropiate source term for the respective form of energy.

Code:

 he.name() == "e"

Then "he" is like as a generic energy that is used in the definition of the EEqn, but in fact it is "h" or "e".

In fvSchemes and in fvSolution you must set for the energy parameters the name "he.name()"

If you choose enthalpy "he.name()=h" you must define "div(phi,h)" in fvSchemes.
If you choose internal energy "he.name()=e" you must define "div(phi,e)" in fvSchemes.

In the OpenFOAM official web site there is an explanation on how to set up the thermophysical properties.

HTML Code:

http://www.openfoam.org/docs/user/thermophysical.php

Regards,

Horacio

zqlhzx · September 12, 2013, 05:05

Thanks ,Horacio Aguerre!It is kind of you!

zqlhzx · September 13, 2013, 10:34

Dear Horacio Aguerre:
Thank you for your last reply!I have another problem I would like to consult you.In my thermophysicalPropertise,I want to define thermo type like this:
thermoType
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
energy absoluteEnthalpy;
equationOfState perfectGas;
specie specie;
}I want use "absoluteEnthalpy;",but when I run my case,it tells me have no this thermo combination.I found just "psiuReactionThermo" Class can use absoluteEnthalpy. "psiReactionThermo"and "rhoReactionThermo"can not use absoluteEnthalpy.They only have sensibleEnthalpy.My case have chemical reaction and I must use absoluteEnthalpy,can you give me advise how to change my thermo type?Thanks again!

santiagomarquezd · October 11, 2013, 14:43

Henry changed the code using the snippet suggested by Horacio, but in pEqn.H

https://github.com/OpenFOAM/OpenFOAM...d6c6406b07d7bb

Regards.

Saicharan · August 19, 2022, 03:47

Quote:

Originally Posted by Horacio Aguerre

Balance of a tensorial property in an arbritary moving domain.

$\dfrac{d}{dt} \int_{V}\rho \phi dV + \int_{S}\rho \phi\left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}=-\int_{S}\rho \mathbf{Q_{\phi}} \cdot \mathbf{ds}+\int_{V}S_{\phi}dV$

Balance of mechanical energy.

$\phi \equiv \frac{\mathbf{v} \cdot \mathbf{v}}{2} = K$

$\dfrac{d}{dt} \int_{V}\rho K dV + \int_{S}\rho K \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S}\left(p\mathbf{v}+ \boldsymbol{\tau} \cdot \mathbf{v}\right) \cdot \mathbf{ds} + \int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV+ \int_{V}\rho\left(\mathbf{v} \cdot \mathbf{g}\right) dV$

Hi Horacio,

"Balance of a tensorial property in an arbitrary moving domain" in the equation left hand side looks like reynolds transport theorem. Can I know what is on RHS.
And the term on RHS of blacance of mechanical energy [p(napla.V) +t:napla.V] physically mean?
Thanks and regards

May 27, 2013, 09:03	Study of the EEqn.H in rhoPimpleDyMFoam.	#1
Horacio Aguerre New Member Join Date: May 2013 Location: Santa Fé, Argentina Posts: 5 Rep Power: 13	Hi, I'm working with the solver "rhoPimpleDyMFoam" from the OpenFOAM 2-2-x version. In specific i was trying to replicate the polytropic process in an adiabatic compression. The case i studied is a pistone-cylinder sistem where air is compressed by the pistone patch which is a moving wall. The pistone velocity is adjusted to 0,1(m/s) to aproximate a reversible process. When i choose energy "he" as internal energy i could reproduce the polytropic process in an adiabatic compression. However when i choose energy "he" as entalphy I couldn't achieve the theoretical expression $pV^{1,4}=constant$ for air. Instead the system looks as it would be non adiabatic or like the conversion between work and internal energy was not made properly. Therefore i started to study the EEqn.H. I compare the theoretical balance of total energy in terms of entalphy with what is written in the EEqn.H. Balance of a tensorial property in an arbritary moving domain. $\dfrac{d}{dt} \int_{V}\rho \phi dV + \int_{S}\rho \phi\left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}=-\int_{S}\rho \mathbf{Q_{\phi}} \cdot \mathbf{ds}+\int_{V}S_{\phi}dV$ Balance of mechanical energy. $\phi \equiv \frac{\mathbf{v} \cdot \mathbf{v}}{2} = K$ $\dfrac{d}{dt} \int_{V}\rho K dV + \int_{S}\rho K \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S}\left(p\mathbf{v}+ \boldsymbol{\tau} \cdot \mathbf{v}\right) \cdot \mathbf{ds} + \int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV+ \int_{V}\rho\left(\mathbf{v} \cdot \mathbf{g}\right) dV$ Balance of internal energy. $\phi \equiv \widehat{u}$ $\dfrac{d}{dt} \int_{V}\rho \widehat{u} dV + \int_{S}\rho \widehat{u} \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV$ The balance of internal energy can be expressed in terms of entalphy. $\widehat{h}=\widehat{u}+\frac{p}{\rho}\implies \widehat{u}=\widehat{h}-\frac{p}{\rho}$ Replacing in the balance of internal energy. $\dfrac{d}{dt} \int_{V}\rho \left(\widehat{h}-\frac{p}{\rho}\right) dV + \int_{S}\rho \left(\widehat{h}-\frac{p}{\rho}\right) \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV$ $\dfrac{d}{dt} \int_{V}\rho \widehat{h} dV+ \int_{S}\rho \widehat{h} \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}-\dfrac{d}{dt} \int_{V}p dV - \int_{S} p \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{V} \left[p \left(\boldsymbol{\nabla} \cdot \mathbf{v}\right)+\boldsymbol{\tau}:\boldsymbol{\nabla v}\right] dV$ Balance of total energy. The sum of the balances of mechanical and internal energy results in the balance of total energy. $\dfrac{d}{dt} \int_{V}\rho \left(\widehat{h}+ K \right) dV + \int_{S}\rho \left(\widehat{h}+ K\right) \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds} -\dfrac{d}{dt} \int_{V}p dV - \int_{S} p \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds}$ $=-\int_{S} \mathbf{q} \cdot \mathbf{ds}-\int_{S}\left(p\mathbf{v}+ \boldsymbol{\tau} \cdot \mathbf{v}\right) \cdot \mathbf{ds} + \int_{V}\rho\left(\mathbf{v} \cdot \mathbf{g}\right) dV$ If the terms of work of viscous forces and the gravity are neglected, the balance is: $\dfrac{d}{dt} \int_{V}\rho \left(\widehat{h}+ K \right) dV + \int_{S}\rho \left(\widehat{h}+ K\right) \left(\mathbf{v}-\mathbf{v_{b}}\right) \cdot \mathbf{ds} -\dfrac{d}{dt} \int_{V}p dV + \int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ $= -\int_{S} \mathbf{q} \cdot \mathbf{ds}$ rhoPimpleDyMFoam: EEqn.H The balance of total energy in the EEqn.H file with considered as entalphy $\widehat{h}$ is: Code: fvm::ddt(rho, he) + fvm::div(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K) - dpdt - fvm::laplacian(turbulence ->alphaEff(), he) == fvOptions(rho, he) Where, $\phi=\rho_{f} \left(\mathbf{v}-\mathbf{v_{b}}\right)_{f} \cdot \mathbf{S_{f}}$ , is relative to the mesh. If fvOptions(rho, he) are considered as zero, the term $\int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ from the final balance is no considered in EEqn.H. Then i try adding in the EEqn.H the therm $\int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ . Code: fvm::ddt(rho, he) + fvm::div(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K) - dpdt + fvc::div(fvc::meshPhi(U),p) - fvm::laplacian(turbulence ->alphaEff(), he) == fvOptions(rho, he) The result was correct. As in the first case (he=internal energy) the adiabatic compression was reproduced with the relation $pV^{1,4}=constant$ . The question is, should the term $\int_{S} p \mathbf{v_{b}} \cdot \mathbf{ds}$ be introduced in any way in fvOptions(rho, he) or it must be explicit in the EEqn.H? Thanks, Horacio mturcios777, fredo490, mgg and 2 others like this.

June 12, 2013, 23:15		#2
santiagomarquezd Senior Member Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 452 Rep Power: 24	Any clues?. It seems we have to upload the question to the bug tracker. Regards. __________________ Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC) - CONICET/UNL Tel: 54-342-4511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe - Argentina. http://www.cimec.org.ar

July 8, 2013, 19:43		#3
santiagomarquezd Senior Member Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 452 Rep Power: 24	http://www.openfoam.org/mantisbt/view.php?id=909 __________________ Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC) - CONICET/UNL Tel: 54-342-4511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe - Argentina. http://www.cimec.org.ar

October 11, 2013, 14:43		#11
santiagomarquezd Senior Member Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 452 Rep Power: 24	Henry changed the code using the snippet suggested by Horacio, but in pEqn.H https://github.com/OpenFOAM/OpenFOAM...d6c6406b07d7bb Regards. __________________ Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC) - CONICET/UNL Tel: 54-342-4511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe - Argentina. http://www.cimec.org.ar

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
grid study, important question	hamid1	FLUENT	1	August 4, 2013 01:14
[ANSYS Meshing] grid study, important questions	hamid1	ANSYS Meshing & Geometry	2	February 10, 2012 14:28
Loop outlet data to inlet of new study	pggow	FLUENT	0	December 21, 2011 15:31
dpm study of solid-liquid study.	bhanwar	Main CFD Forum	0	September 29, 2011 09:21
Study	Mario	Main CFD Forum	0	August 24, 2009 10:59

September 5, 2013, 03:20		#6
ALU New Member Join Date: Jun 2012 Posts: 11 Rep Power: 14	Hi Horacio Aguerre, thank you very much for the correction of my case and the detailed explanation, this really helps a lot! Regards

September 11, 2013, 01:32		#7
zqlhzx Member 赵庆良 Join Date: Aug 2013 Posts: 56 Rep Power: 13	Hi ,Horacio Aguerre could you see my problem and tell some thing?I think you can answer my question.I moves fireFoam of 2.0.1 to 2.2.1,then rewrite hEqn according to the EEqn.H in fireFoam of 2.2.1 and other necessary codes.But ,when I wmake my rewrited fireFoam,it gives me errors:undefine div(phi,h) in fvscheme. The following is codes of EEqn.H: { volScalarField& he = thermo.he(); fvScalarMatrix EEqn ( fvm::ddt(rho, he) + fvm::div(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K) + ( he.name() == "e" ? fvc::div ( fvc::absolute(phi/fvc::interpolate(rho), U), p, "div(phiv,p)" ) : -DpDt ) - fvm::laplacian(turbulence->alphaEff(), he) == fvOptions(rho, he) ); EEqn.relax(); fvOptions.constrain(EEqn); EEqn.solve(); fvOptions.correct(he); thermo.correct(); Info<< "T min/max = " << min(T).value() << ", " << max(T).value() << endl; } There is div(phi,he) insteading of div(phi,h),So It should be no need to define (phi,h) in fvshceme,but define div(phi,he).However,I see the fvsheme in smallpoolfire2D of 2.2.1,there indeed is div(phi,h).In fact,there is the term of "div(phi,he)" in EEqn.H,why define div(phi,h) in fvsheme?I think I have some thing not understand.!

September 12, 2013, 05:05		#9
zqlhzx Member 赵庆良 Join Date: Aug 2013 Posts: 56 Rep Power: 13	Thanks ,Horacio Aguerre!It is kind of you!

September 13, 2013, 10:34		#10
zqlhzx Member 赵庆良 Join Date: Aug 2013 Posts: 56 Rep Power: 13	Dear Horacio Aguerre: Thank you for your last reply!I have another problem I would like to consult you.In my thermophysicalPropertise,I want to define thermo type like this: thermoType { type hePsiThermo; mixture reactingMixture; transport sutherland; thermo janaf; energy absoluteEnthalpy; equationOfState perfectGas; specie specie; }I want use "absoluteEnthalpy;",but when I run my case,it tells me have no this thermo combination.I found just "psiuReactionThermo" Class can use absoluteEnthalpy. "psiReactionThermo"and "rhoReactionThermo"can not use absoluteEnthalpy.They only have sensibleEnthalpy.My case have chemical reaction and I must use absoluteEnthalpy,can you give me advise how to change my thermo type?Thanks again!