|
[Sponsors] |
October 20, 2011, 09:36 |
Initial condition in OpenFoam
|
#1 |
Member
Gitesh
Join Date: Jan 2010
Location: Finland
Posts: 73
Rep Power: 16 |
Hello,
I want to initialize hydrostatic pressure in my simulation. I have tried it to do it by 'funkySetFields -time 0'. It works for zero time but at first iteration it suddenly changes. I have pressure inlet and pressure outlet as a BC's. In 0/p I have given fixed values for pressure uniform 101320(for outlet) and uniform 118600 (for inlet). Can anyone help me in this problem ? Best regards, Gitesh |
|
October 20, 2011, 22:56 |
|
#2 |
Senior Member
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16 |
What solver are you using?
Can you provide a little more information about what you have tried? Step by step? Did you use setFields to initialize water surface? |
|
October 21, 2011, 04:44 |
|
#3 |
Member
Gitesh
Join Date: Jan 2010
Location: Finland
Posts: 73
Rep Power: 16 |
Hello mgdenno,
OK. So my system is air is feeding by pipe into water tank. I am using twoPhaseEulerFoam solver. I am using 'funkySetFields' for case initialization. It works ok. But at after first iteration it calculate different things. You can see in pictures (pressure) in attached file. Best regards, Gitesh |
|
October 21, 2011, 16:00 |
|
#4 |
Senior Member
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16 |
Hi Gitesh,
I am not at all familiar with twoPhaseEulerFoam, but which variables did you initialize? pressure? phase? Looks like maybe you only initialized the pressure but not the phase? Also which way is gravity acting? MD Last edited by mgdenno; October 21, 2011 at 16:38. Reason: Elaborate |
|
October 22, 2011, 03:21 |
|
#5 |
Member
Gitesh
Join Date: Jan 2010
Location: Finland
Posts: 73
Rep Power: 16 |
Hello,
I also initialized phase. It works ok but the problem in pressure. Gravity is in +ve x direction. BR, Gitesh |
|
October 23, 2011, 03:06 |
|
#6 | |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
Best,
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
||
October 23, 2011, 03:46 |
|
#7 |
Member
Gitesh
Join Date: Jan 2010
Location: Finland
Posts: 73
Rep Power: 16 |
Hello alberto,
I want to initialize hydrostatic pressure for water. Also I have to give pressure inlet value because air is coming from there. Moreover I have pressure outlet there from where air will be out. So, I have used funkySetFields for hydrostatic pressure initialization. So, is there any other bc from where we can directly calculate pressure gradient with out initialization ? Thank you !! BR, Gitesh |
|
October 23, 2011, 03:53 |
|
#8 | |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
Best,
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
||
October 24, 2011, 02:08 |
|
#9 |
Member
Gitesh
Join Date: Jan 2010
Location: Finland
Posts: 73
Rep Power: 16 |
Hello Albarto,
Thank you !! So, you mean I have to give velocity of air at inlet and to initialize pressure gradient by funkySetFields ? BR, Gitesh |
|
October 24, 2011, 02:11 |
|
#10 | |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
1) Specify velocity at the inlet, and set the condition on the pressure to zeroGradient there. 2) Fix the value of the pressure at the outlet, and set the velocity to zeroGradient or inletOutlet, depending on your case. 3) You can use a uniform initialization for the pressure in this case, since it does not really matter. The solver will find the correct pressure filed at the first time step. Best,
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
||
May 27, 2013, 10:08 |
|
#11 | |
Senior Member
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 15 |
Quote:
i am working on 2D multi element airfoils, i have a Question about 2th suggestion. "Fix the value of the pressure at the outlet, and set the velocity to zeroGradient or inletOutlet, depending on your case." what do you mean with depending on your case,in my case which one is better? inletOutlet or fixedValue ? |
||
May 27, 2013, 19:08 |
|
#12 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
hi
inletOutlet is an "outlet" BC not inlet.this is the same zeroGradient only with a difference. When there is a backflow in outlet,inletOutlet BC sets a value on the cells that have backflow by the value we have specified. So if you think may you have backflow use inletOulet otherwise zeroGradient suffices. Hope it helps.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
May 28, 2013, 03:32 |
|
#13 | |
Senior Member
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 15 |
Quote:
i understand what you said, thank you very much. if my domain is large enough around my airfoil, it doesn’t need to use the inletOutlet boundary condition, is it right? |
||
May 28, 2013, 09:15 |
|
#14 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
yes,inletOutlet is used in internal flows in common.
I didn't received your case you said that had sent.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
May 28, 2013, 11:44 |
|
#15 | |
Senior Member
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 15 |
Quote:
i put my multi element airfoil, thank you very much for your help. |
||
May 19, 2016, 16:40 |
|
#16 |
Senior Member
CFD_Lovers
Join Date: Mar 2015
Posts: 168
Rep Power: 11 |
hi friends! i have the problem of hydrostatic pressure initializing in a tank full of water like Gitesh P!
and what is funkysetField? how can i set hydrostatic pressure condition in my domain ? i shoul say my inlet is not the whole left side of my geometry , it is a little nozzle in there (left side of geometry) so what can i do? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
icoLagrangianFoam OF1.6 myNewParticleSolver | heavy_user | OpenFOAM | 23 | June 2, 2020 03:18 |
SimpleFoam k and epsilon bounded | nedved | OpenFOAM Running, Solving & CFD | 16 | March 4, 2017 09:30 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 16:33 |
lift and drag on ship superstructures | vaina74 | OpenFOAM Running, Solving & CFD | 3 | June 8, 2010 13:30 |
SimpleFoam k and epsilon bounded | nedved | OpenFOAM Running, Solving & CFD | 1 | November 25, 2008 21:21 |