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

Diverging result for Temperature field in interFoam

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

Like Tree10Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   December 6, 2010, 12:09
Default
  #61
ala
New Member
 
Alicja M
Join Date: Mar 2009
Location: Erlangen, DE
Posts: 26
Rep Power: 7
ala is on a distinguished road
Hey,

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

Have fun!
Greetings, Alicja
Attached Files
File Type: zip 1d_heat_diff.zip (14.8 KB, 72 views)
File Type: zip disDiffusion.zip (7.0 KB, 62 views)
File Type: zip 1dCode.zip (1.5 KB, 72 views)
ala is offline   Reply With Quote

Old   January 23, 2011, 01:50
Default
  #62
New Member
 
zhenyan xia
Join Date: Jan 2011
Location: Tianjin,China
Posts: 6
Rep Power: 5
xiazhy is on a distinguished road
Send a message via MSN to 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
xiazhy is offline   Reply With Quote

Old   January 23, 2011, 02:04
Default
  #63
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 7
chiven is on a distinguished road
Quote:
Originally Posted by xiazhy View Post
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
chiven is offline   Reply With Quote

Old   January 23, 2011, 02:38
Default
  #64
New Member
 
zhenyan xia
Join Date: Jan 2011
Location: Tianjin,China
Posts: 6
Rep Power: 5
xiazhy is on a distinguished road
Send a message via MSN to xiazhy
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
xiazhy is offline   Reply With Quote

Old   February 3, 2011, 09:32
Default
  #65
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 168
Rep Power: 7
linch is on a distinguished road
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 )

Last edited by linch; February 4, 2011 at 08:03.
linch is offline   Reply With Quote

Old   February 3, 2011, 10:42
Default
  #66
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 168
Rep Power: 7
linch is on a distinguished road
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.
linch is offline   Reply With Quote

Old   March 4, 2011, 05:54
Default
  #67
ala
New Member
 
Alicja M
Join Date: Mar 2009
Location: Erlangen, DE
Posts: 26
Rep Power: 7
ala is on a distinguished road
Hi Foamers,
I have solved my problem in 1D (please see here: Diverging result for Temperature field in interFoam). 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 of
the 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
ala is offline   Reply With Quote

Old   September 4, 2011, 02:32
Default
  #68
New Member
 
Join Date: Nov 2010
Posts: 10
Rep Power: 6
siavash_abadeh is on a distinguished road
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 is offline   Reply With Quote

Old   September 4, 2011, 08:09
Default
  #69
New Member
 
Join Date: Nov 2010
Posts: 10
Rep Power: 6
siavash_abadeh is on a distinguished road
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 is offline   Reply With Quote

Old   September 4, 2011, 08:26
Default
  #70
New Member
 
Join Date: Nov 2010
Posts: 10
Rep Power: 6
siavash_abadeh is on a distinguished road
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
siavash_abadeh is offline   Reply With Quote

Old   September 5, 2011, 05:43
Default
  #71
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 168
Rep Power: 7
linch is on a distinguished road
What's wrong with this one?
linch is offline   Reply With Quote

Old   January 25, 2012, 17:25
Default
  #72
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 175
Rep Power: 6
vonboett is on a distinguished road
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?
vonboett is offline   Reply With Quote

Old   February 10, 2012, 16:46
Default
  #73
Member
 
Ovie Doro
Join Date: Jul 2009
Posts: 99
Rep Power: 7
ovie is on a distinguished road
Quote:
Originally Posted by vonboett View Post
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.
ovie is offline   Reply With Quote

Old   February 13, 2012, 04:16
Default
  #74
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 175
Rep Power: 6
vonboett is on a distinguished road
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
vonboett is offline   Reply With Quote

Old   February 15, 2012, 16:31
Default temperature in interFoam
  #75
New Member
 
Hossein
Join Date: Feb 2012
Posts: 13
Rep Power: 4
Gildeh is on a distinguished road
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
Gildeh is offline   Reply With Quote

Old   September 7, 2013, 08:23
Default how add source term to inetrFoam
  #76
New Member
 
sasan
Join Date: Sep 2013
Posts: 28
Rep Power: 3
Sasy is on a distinguished road
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??
Sasy is offline   Reply With Quote

Old   September 24, 2013, 20:38
Default
  #77
Member
 
Ovie Doro
Join Date: Jul 2009
Posts: 99
Rep Power: 7
ovie is on a distinguished road
Quote:
Originally Posted by Sasy View Post
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.
ovie is offline   Reply With Quote

Old   September 30, 2013, 04:38
Default
  #78
Member
 
vishal
Join Date: Mar 2013
Posts: 54
Rep Power: 3
vishal_s is on a distinguished road
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
vishal_s is offline   Reply With Quote

Old   December 19, 2013, 05:06
Default
  #79
New Member
 
thejaraju
Join Date: Dec 2010
Posts: 8
Rep Power: 5
thejaraju is on a distinguished road
Quote:
Originally Posted by ovie View Post
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
thejaraju is offline   Reply With Quote

Old   December 19, 2013, 06:22
Default
  #80
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
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
Bernhard is offline   Reply With Quote

Reply

Tags
interfoam energy, temperature field

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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
Getting a concentration field around a bubble in InterFoam azman OpenFOAM Running, Solving & CFD 2 January 7, 2013 14:47
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


All times are GMT -4. The time now is 15:43.