Diverging result for Temperature field in interFoam

ala · December 6, 2010, 12:09

Hey,

and here come: my 1dCode and OF solver and test case.

Have fun!
Greetings, Alicja

xiazhy · January 23, 2011, 01:50

Hi chiven

i am a fresher in OF..I want to simulate the high speed Non_newtonian jet in air. i want to know which solver i should use, and i hope i can get some adivce from you about LES+VOF+Non_newtonian. Thanks a lot ! any help whould be much appreciated!

zhenyan

chiven · January 23, 2011, 02:04

Quote:

Originally Posted by xiazhy

Hi chiven

i am a fresher in OF..I want to simulate the high speed Non_newtonian jet in air. i want to know which solver i should use, and i hope i can get some adivce from you about LES+VOF+Non_newtonian. Thanks a lot ! any help whould be much appreciated!

zhenyan

Hi, zhenyan, cheering. If you are trying to use the VOF model, maybe interFoam solver is enough for you. Best regards, Chiven

xiazhy · January 23, 2011, 02:38

Hi chiven

Thanks a lot for your so quick reply. i will have a try, and if i meet question i will ask you for advice. May i have your email address? My email address is liping_he@sina.cn Thanks!

Goodluck

zhenyan

linch · February 3, 2011, 09:32

Hi foamers,

I went on and tried to implement temperature dependent density (just temperature dependent, without pressure dependency).

Because the velocity isn't divergence-free any more, I had to add source terms into the pEqn.H and alphaEqn.H.

My p_rghEqn now looks like:

Code:

fvm::laplacian(rAUf, p_rgh) == fvc::div(phi)+ (alpha1*fvc::ddt(rho1)+(scalar(1)-alpha1)*fvc::ddt(rho2))/rho

The problem is that in this case my pressure field starts to oscillate and I get chessboard looking like pressure regions alternating every time step. The mean pressure in these regions is ok.

Do you have any suggestions to help it? Discretisation schemes?

Regads,
Ilya

P.S. The oscillation occur more intensively in near-wall regions and not in regions with the highest velocity divergence (== highest temperature gradients )

linch · February 3, 2011, 10:42

Here is an example:
div(phi) due to density change and p-field oscillations in the domain corner (zoom)

The scales of both pics are different, only one corner of the computational domain is shown in the right one but it looks similar for all four corners.

Near wall regions are affected stronger. Oscillations occur not immediately and the their amplitude increases with time.

ala · March 4, 2011, 05:54

Hi Foamers,

I have solved my problem in 1D (please see here: http://www.cfd-online.com/Forums/ope...terfoam-4.html). If I have a moving interface I also need a convection term, but i'm a student ; ).

Unfortunately, I still have problems with the energy equation in 2D, especially using the interFoam solver, like Sascha612.

Some hints? It would be great to know some literature about two phase flow with energy equation. My mail address is <alicja.mura at googlemail.com>.

I also think, that I have problems because ofthe gamma equation discretisation. Do you have further information (literature, papers, and so on) about it? I have seen that some people cite technical reports of OpenFOAM .

greetings,
alicja

siavash_abadeh · September 4, 2011, 02:32

hi dears
i have a case that is 2 phase with heat transfer.
i want to add energy eq to interfoam.
is there anyone who can help me.
thanks in advance

siavash_abadeh · September 4, 2011, 08:09

hi dears
i have a case that i want to solve energy equation with interfoam(vof)
how can i add this equation.
is there anyone who can help me.
thanks in advance

siavash_abadeh · September 4, 2011, 08:26

hi dears
i have a case that i want to solve energy equation to interfoam.
is there anyone who can help me(is there ready code).
thanks in advance

linch · September 5, 2011, 05:43

What's wrong with this one?

vonboett · January 25, 2012, 17:25

Ovies code looks promising to me. I want to adress a different problem in the same way. Simulating shallow landslides, the different fluid mixtures have a viscosity that strongly depends on the water content. Instead of temperature exchange, I want to model an exchange of water content between the phases. Anyone tried with multiphaseInterFoam yet?

ovie · February 10, 2012, 16:46

Quote:

Originally Posted by vonboett

Ovies code looks promising to me. I want to adress a different problem in the same way. Simulating shallow landslides, the different fluid mixtures have a viscosity that strongly depends on the water content. Instead of temperature exchange, I want to model an exchange of water content between the phases. Anyone tried with multiphaseInterFoam yet?

How do you intend to model the exchange of water between phases? Is it a concentration driven diffusion process for water mass fraction? Interfoam assumes the phases are incompressible i.e. no mixing, which means the water has to be converted to the new phase its been transported/converted into. I really dont understand the problem definition or the physics of the problem. Maybe you should kindly elaborate some more.

Thanks.

vonboett · February 13, 2012, 04:16

Hi Ovie,

yes, I want to adress the problem as a concentration driven diffusion process, although that is not quite right. The change of water concentration in a phase of shallow landslide material does not affect the density of that material as much as its rheology. So in a first approach, I would leave the density of the phase same as it is but only change the local viscosity in dependency to the local water concentration. This water concentration however has to follow the flow of the phase aswell as a diffusion driven by the water concnetration gradient. It is still at the stage of an idea, but maybe there will be some experimental data available soon which makes the "try out" an interesting opportunity. It would turn into a bigger research project that treates the physics in an appropriate way if the results of this test would turn out to be promising.

So, this is still the brainstorm stage, any Idea is welcome :-)

Best wishes,

Albi

Gildeh · February 15, 2012, 16:31

Hi Ovie and every body,

However, that is a long time you posted these notes here, I hope you can help me. I am doing the same job as you with interFoam solver. I am pretty new to OF and trying to solve my problem with that. I did your procedures you mentioned here. But 2 problems:
1. You compiled the files twoPhaseMixture.H and twoPhaseMixture.C in the src file where you installed OF. I copied incompressibleTwoPhaseMixture file in another file and changed those two files there. What should I do if I want to compile those files in this new place and then use them in my new solver that is intended to work instead of interFoam.
2. I could not compile the previous step and just jumped to the second step (step 3 up to 7 in your procedure) and tried to compile my solver that again I faced errors that I think I have to do the first step firstly. Am I right?

I am pretty confused with this problem and am not sure that I could implement your idea properly.

I appreciate if you or any one else can help.

Gildeh

Sasy · September 7, 2013, 08:23

Hi Ovie and every body,
i want simulate evaporation and condenstion in vertical tube,i use interFoam solver and i konw i should modify this solver by adding Temp.eqn and source term.i can implemented the temprature equation in interFoam As said ovie..how add mass source in interFoam??

ovie · September 24, 2013, 20:38

Quote:

Originally Posted by Sasy

Hi Ovie and every body,
i want simulate evaporation and condenstion in vertical tube,i use interFoam solver and i konw i should modify this solver by adding Temp.eqn and source term.i can implemented the temprature equation in interFoam As said ovie..how add mass source in interFoam??

Hi:

I think there are a few papers where evaporation is modeled as source term for a continuum e.g. Hardt, S. & Wondra, F. 2008 Evaporation model for interfacial flows based on a continuum-field representation of the source terms. This is a very good point to start. The paper also cites a number of helpful references as well.

Good luck.

vishal_s · September 30, 2013, 04:38

hello, I am interested to add temperature in interFoam. But in OpenFOAM 2.2.1 in the multiphase ----> interFoam directory, there is no twoPhaseMixture.H or twoPhaseMixture.C file. Do i need them to incorporate them or the temperature model can work without it. I am working on two phase Jet flow.

Kindly reply

thejaraju · December 19, 2013, 05:06

Quote:

Originally Posted by ovie

No worries and good luck!

Hi ovie

This is thejaraj from bangalore....
from the past 2 weeks i am trying to install the viscoelasticInterfoam for die swell but i am getting the following......

In file included from viscoelasticInterFoam.C:75:0:
alphaEqnSubCycle.H: In function ‘int main(int, char**)’:
alphaEqnSubCycle.H:3:15: error: ‘piso’ was not declared in this scope
/home/theja/OpenFOAM/OpenFOAM-2.2.x/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nOuterCorr’ [-Wunused-variable]
/home/theja/OpenFOAM/OpenFOAM-2.2.x/src/finiteVolume/lnInclude/readPISOControls.H:15:16: warning: unused variable ‘transonic’ [-Wunused-variable]
/home/theja/OpenFOAM/OpenFOAM-2.2.x/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nOuterCorr’ [-Wunused-variable]
/home/theja/OpenFOAM/OpenFOAM-2.2.x/src/finiteVolume/lnInclude/readPISOControls.H:6:15: warning: unused variable ‘nCorr’ [-Wunused-variable]
/home/theja/OpenFOAM/OpenFOAM-2.2.x/src/finiteVolume/lnInclude/readPISOControls.H:12:16: warning: unused variable ‘momentumPredictor’ [-Wunused-variable]
/home/theja/OpenFOAM/OpenFOAM-2.2.x/src/finiteVolume/lnInclude/readPISOControls.H:15:16: warning: unused variable ‘transonic’ [-Wunused-variable]
/home/theja/OpenFOAM/OpenFOAM-2.2.x/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable]
make: *** [Make/linux64GccDPOpt/viscoelasticInterFoam.o] Error 1
theja@theja:~/OpenFOAM/OpenFOAM-2.2.x/applications/solvers/viscoelasticInterFoam$ ls
alphaEqn.H alphaEqnSubCycle.H correctPhi.H createFields.H Make pEqn.H viscoelasticInterFoam.C viscoelasticInterFoam.dep viscoelasticModel.H
alphaEqn.H~ alphaEqnSubCycle.H~ correctPhi.H~ createFields.H~ newViscoelasticLaw.C UEqn.H viscoelasticInterFoam.C~ viscoelasticLaw.H
theja@theja:~/OpenFOAM/OpenFOAM-2.2.x/applications/solvers/viscoelasticInterFoam$

I hope you can help me in solving this error

thanking you in advance

Thejaraju

Bernhard · December 19, 2013, 06:22

You are working in OpenFOAM 2.2.x. The solver you are trying to compare has probably been developed for some older of version of OpenFOAM. Please check what has been changed between this previous version and the one you are using. Incorporate these changes in the solver you are trying to compile. Then you're done

February 3, 2011, 09:32		#65
linch Senior Member Illya Shevchuk Join Date: Aug 2009 Location: Darmstadt, Germany Posts: 176 Rep Power: 17	Hi foamers, I went on and tried to implement temperature dependent density (just temperature dependent, without pressure dependency). Because the velocity isn't divergence-free any more, I had to add source terms into the pEqn.H and alphaEqn.H. My p_rghEqn now looks like: Code: fvm::laplacian(rAUf, p_rgh) == fvc::div(phi)+ (alpha1fvc::ddt(rho1)+(scalar(1)-alpha1)fvc::ddt(rho2))/rho The problem is that in this case my pressure field starts to oscillate and I get chessboard looking like pressure regions alternating every time step. The mean pressure in these regions is ok. Do you have any suggestions to help it? Discretisation schemes? Regads, Ilya P.S. The oscillation occur more intensively in near-wall regions and not in regions with the highest velocity divergence (== highest temperature gradients ) Last edited by linch; February 4, 2011 at 08:03.

February 3, 2011, 10:42		#66
linch Senior Member Illya Shevchuk Join Date: Aug 2009 Location: Darmstadt, Germany Posts: 176 Rep Power: 17	Here is an example: div(phi) due to density change and p-field oscillations in the domain corner (zoom) The scales of both pics are different, only one corner of the computational domain is shown in the right one but it looks similar for all four corners. Near wall regions are affected stronger. Oscillations occur not immediately and the their amplitude increases with time. Last edited by linch; February 4, 2011 at 09:53.

February 15, 2012, 16:31	temperature in interFoam	#75
Gildeh New Member Hossein Join Date: Feb 2012 Posts: 13 Rep Power: 14	Hi Ovie and every body, However, that is a long time you posted these notes here, I hope you can help me. I am doing the same job as you with interFoam solver. I am pretty new to OF and trying to solve my problem with that. I did your procedures you mentioned here. But 2 problems: 1. You compiled the files twoPhaseMixture.H and twoPhaseMixture.C in the src file where you installed OF. I copied incompressibleTwoPhaseMixture file in another file and changed those two files there. What should I do if I want to compile those files in this new place and then use them in my new solver that is intended to work instead of interFoam. 2. I could not compile the previous step and just jumped to the second step (step 3 up to 7 in your procedure) and tried to compile my solver that again I faced errors that I think I have to do the first step firstly. Am I right? I am pretty confused with this problem and am not sure that I could implement your idea properly. I appreciate if you or any one else can help. Gildeh

September 7, 2013, 08:23	how add source term to inetrFoam	#76
Sasy New Member sasan Join Date: Sep 2013 Posts: 28 Rep Power: 13	Hi Ovie and every body, i want simulate evaporation and condenstion in vertical tube,i use interFoam solver and i konw i should modify this solver by adding Temp.eqn and source term.i can implemented the temprature equation in interFoam As said ovie..how add mass source in interFoam??

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Getting a concentration field around a bubble in InterFoam	azman	OpenFOAM Running, Solving & CFD	3	June 7, 2022 05:21
Moving mesh	Niklas Wikstrom (Wikstrom)	OpenFOAM Running, Solving & CFD	122	June 15, 2014 07:20
Adding temperature field to InterFoam	yapalparvi	OpenFOAM Running, Solving & CFD	8	October 14, 2009 21:18
Problem with rhoSimpleFoam	matteo_gautero	OpenFOAM Running, Solving & CFD	0	February 28, 2008 07:51
Problems calculating field gh with interFoam	cricke	OpenFOAM Running, Solving & CFD	0	December 10, 2007 08:17

January 23, 2011, 01:50		#62
xiazhy New Member zhenyan xia Join Date: Jan 2011 Location: Tianjin,China Posts: 6 Rep Power: 15	Hi chiven i am a fresher in OF..I want to simulate the high speed Non_newtonian jet in air. i want to know which solver i should use, and i hope i can get some adivce from you about LES+VOF+Non_newtonian. Thanks a lot ! any help whould be much appreciated! zhenyan

January 23, 2011, 02:38		#64
xiazhy New Member zhenyan xia Join Date: Jan 2011 Location: Tianjin,China Posts: 6 Rep Power: 15	Hi chiven Thanks a lot for your so quick reply. i will have a try, and if i meet question i will ask you for advice. May i have your email address? My email address is liping_he@sina.cn Thanks! Goodluck zhenyan

March 4, 2011, 05:54		#67
ala New Member Alicja M Join Date: Mar 2009 Location: Erlangen, DE Posts: 26 Rep Power: 17	Hi Foamers, I have solved my problem in 1D (please see here: http://www.cfd-online.com/Forums/ope...terfoam-4.html). If I have a moving interface I also need a convection term, but i'm a student ; ). Unfortunately, I still have problems with the energy equation in 2D, especially using the interFoam solver, like Sascha612. Some hints? It would be great to know some literature about two phase flow with energy equation. My mail address is <alicja.mura at googlemail.com>. I also think, that I have problems because ofthe gamma equation discretisation. Do you have further information (literature, papers, and so on) about it? I have seen that some people cite technical reports of OpenFOAM . greetings, alicja

September 4, 2011, 02:32		#68
siavash_abadeh New Member Join Date: Nov 2010 Posts: 10 Rep Power: 16	hi dears i have a case that is 2 phase with heat transfer. i want to add energy eq to interfoam. is there anyone who can help me. thanks in advance

September 4, 2011, 08:09		#69
siavash_abadeh New Member Join Date: Nov 2010 Posts: 10 Rep Power: 16	hi dears i have a case that i want to solve energy equation with interfoam(vof) how can i add this equation. is there anyone who can help me. thanks in advance

September 4, 2011, 08:26		#70
siavash_abadeh New Member Join Date: Nov 2010 Posts: 10 Rep Power: 16	hi dears i have a case that i want to solve energy equation to interfoam. is there anyone who can help me(is there ready code). thanks in advance

September 5, 2011, 05:43		#71
linch Senior Member Illya Shevchuk Join Date: Aug 2009 Location: Darmstadt, Germany Posts: 176 Rep Power: 17	What's wrong with this one?

January 25, 2012, 17:25		#72
vonboett Senior Member Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 240 Rep Power: 17	Ovies code looks promising to me. I want to adress a different problem in the same way. Simulating shallow landslides, the different fluid mixtures have a viscosity that strongly depends on the water content. Instead of temperature exchange, I want to model an exchange of water content between the phases. Anyone tried with multiphaseInterFoam yet?

February 13, 2012, 04:16		#74
vonboett Senior Member Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 240 Rep Power: 17	Hi Ovie, yes, I want to adress the problem as a concentration driven diffusion process, although that is not quite right. The change of water concentration in a phase of shallow landslide material does not affect the density of that material as much as its rheology. So in a first approach, I would leave the density of the phase same as it is but only change the local viscosity in dependency to the local water concentration. This water concentration however has to follow the flow of the phase aswell as a diffusion driven by the water concnetration gradient. It is still at the stage of an idea, but maybe there will be some experimental data available soon which makes the "try out" an interesting opportunity. It would turn into a bigger research project that treates the physics in an appropriate way if the results of this test would turn out to be promising. So, this is still the brainstorm stage, any Idea is welcome :-) Best wishes, Albi

September 30, 2013, 04:38		#78
vishal_s Member vishal Join Date: Mar 2013 Posts: 73 Rep Power: 13	hello, I am interested to add temperature in interFoam. But in OpenFOAM 2.2.1 in the multiphase ----> interFoam directory, there is no twoPhaseMixture.H or twoPhaseMixture.C file. Do i need them to incorporate them or the temperature model can work without it. I am working on two phase Jet flow. Kindly reply

December 19, 2013, 06:22		#80
Bernhard Senior Member Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 22	You are working in OpenFOAM 2.2.x. The solver you are trying to compare has probably been developed for some older of version of OpenFOAM. Please check what has been changed between this previous version and the one you are using. Incorporate these changes in the solver you are trying to compile. Then you're done