CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

compressible flow with phase change

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

Like Tree7Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   September 17, 2013, 06:52
Default
  #21
New Member
 
Join Date: Sep 2013
Posts: 21
Rep Power: 4
ma-tri-x is on a distinguished road
Ah I see. Would be very great to use it! I also got now the meaning of beta and I've rewritten the simpleLiquid to the tait-Equation of state:

rho=rho_0*((p+B)/(p_0 +B))^(1/gamma); gamma=7.15

B = prefactorA + prefactorB*T[°C] + prefactorC*T²[°C²]

psi = d rho/d p

rhoHat not used (set to scalar(1.0)). B is calculated in thermoCorrection(), rho and psi are calculated in psiCorrection()

But it's not very stable and I don't know why... it does not seem to exceed the courant numbers...

Anyway, for sufficiently small starting values for my bubble (r=160mum, alpha=0, T_0<600K, p_rgh_0<4e6Pa) the picture really looks like something expected. More than with any solver before. So it would be really great to get your update soon!

Are you going to change the solver routines or the equations?
ma-tri-x is offline   Reply With Quote

Old   October 30, 2013, 09:48
Default
  #22
New Member
 
weiyi xue
Join Date: Oct 2013
Posts: 10
Rep Power: 4
simon95 is on a distinguished road
Hi Matvej,

Thanks for your great job.

I am looking for the solver that is compatible with openfoam2.2.x. If you complete a porting of version. please publish it on this site.

Thanks

Simon
simon95 is offline   Reply With Quote

Old   November 1, 2013, 05:55
Default
  #23
New Member
 
weiyi xue
Join Date: Oct 2013
Posts: 10
Rep Power: 4
simon95 is on a distinguished road
Sorry, I had a very simple question. Regarding tutorial water-flashing, I think that phase2 is atmosphere and rho of phase2 should be set to 1 in file transportProperties. I could not understand why rho of phase2 is set to 0.1 in transportProperties.

Thanks in advance.

Simon
simon95 is offline   Reply With Quote

Old   January 22, 2014, 18:13
Default
  #24
Member
 
Jace
Join Date: Oct 2012
Posts: 77
Rep Power: 6
zhengzh5 is on a distinguished road
Hi,

just out of curiosity...

you mentioned that the governing equations have been rewritten along with other significant changes to the solver. Do you mind briefly touched upon what exactly are the improvements, and what are the motivations behind them. I assume there's some flaws in the previous set of governing equations?

Thanks!

Jason
zhengzh5 is offline   Reply With Quote

Old   January 24, 2014, 05:13
Default
  #25
Senior Member
 
Phoevos
Join Date: Mar 2009
Posts: 103
Rep Power: 9
fivos is on a distinguished road
Hi,

I tried to test the new solver (successfully compiled from provided source code for OpenFOAM 2.1.x) on one of the tutorial cases, but it does not work.

I tried to run the hemispherical body (http://openfoamwiki.net/index.php/Co...haseChangeFoam, case 1), but it always diverges at some point at 70ms. I tried to run with constant time step and at 69ms Courant No blows up. I switched to auto-timestepping with max. Co. 0.05, it goes a little bit further, but still crashes.

Has anyone managed to run this particular case for 20s (or more), as it is in the controlDict?

Generally, when running I notice that there are significant oscillations in pressure field and I also get phase change (cavitation) at the slip wall, which does not seem reasonable...

I have not modified any of the tutorial settings or the solver itself...

Any ideas on that?

Thanks in advance.
fivos is offline   Reply With Quote

Old   March 19, 2014, 11:48
Default
  #26
New Member
 
Juan Carlos
Join Date: Feb 2010
Posts: 9
Rep Power: 8
Juan Carlos is on a distinguished road
Hi Matvej,
I would like to use your solver for compressible two-phase flow with phase change, and I would like to know if the new version of this solver for version OF2.2.X it's available
Thank you in advance for sharing your work
Best Regards
Juan C.
Juan Carlos is offline   Reply With Quote

Old   April 24, 2014, 14:50
Default Follow-up
  #27
New Member
 
Joseph Prince
Join Date: Apr 2014
Posts: 1
Rep Power: 0
jfp6 is on a distinguished road
Were you ever able to finish this solver?

Thanks,

Quote:
Originally Posted by mkraposhin View Post
Hello, ma-tri-x!

I'm sorry for the late response!

Today i'm working on the new version of this solver. The system of equations has been rewritten totally, so, please, wait about 2-4 weeks.
Of course, i will made versions for OF2.1.X and OF2.2.X

Changes are very significant,
herlonyoshio likes this.
jfp6 is offline   Reply With Quote

Old   January 5, 2015, 11:17
Default
  #28
Member
 
Christian
Join Date: Apr 2014
Posts: 69
Rep Power: 4
crixman is on a distinguished road
Hi Marco,
any news on the hydraulic resistance test case?
I would like to test it as well! Also, having the solver for OF2.2.x would be great


Quote:
Originally Posted by sfigato View Post

Hi,

I already sent to you an mail!

Regarding the test case, I have also an hydraulic resistence with mineral oil. Moreover, I am preparing an hydraulic valve with cavitation (the geometry looks like really cool).. but it takes time to be ready!

Regards
Marco
crixman is offline   Reply With Quote

Old   April 1, 2015, 18:49
Default
  #29
New Member
 
zhouhoucun
Join Date: Dec 2014
Posts: 11
Rep Power: 3
zhouhoucun is on a distinguished road
Quote:
Originally Posted by sfigato View Post
Hi,

here is the link, where you can download the description of the model:

http://fluidsengineering.asmedigital...icleid=1429576

The model includes also the gas cavitation! You cannot find the full description of the model with the gas mass transfer (one more transport equation)! If you give to me your mail I will send to you it!

thank you for sharing your work
Best Regards

Marco
Dear Macro, could share me your code plesse? My email address is 470861844@qq.com, Thank you in advance.

Yours,
Zhou
zhouhoucun is offline   Reply With Quote

Old   July 21, 2015, 10:04
Default
  #30
Senior Member
 
rkhr
Join Date: May 2011
Posts: 231
Rep Power: 8
Kanarya is on a distinguished road
Hi All,

I would like to use the code as well but I could not find it.
Is the code still available?

Thanks!
Kanarya is offline   Reply With Quote

Old   October 16, 2015, 02:50
Default
  #31
New Member
 
Federica
Join Date: Oct 2015
Posts: 13
Rep Power: 2
fkika is on a distinguished road
Hi matvey.

I am not able to reach the page for the download.
The solver is still available? Can I implemeted it in OpenFOAM-2.3?

Thank you

Federica
fkika is offline   Reply With Quote

Old   January 14, 2016, 09:39
Default Link for download
  #32
bhh
Member
 
Bjorn H. Hjertager
Join Date: Mar 2009
Posts: 71
Rep Power: 9
bhh is on a distinguished road
Quote:
Originally Posted by mkraposhin View Post
Hi, i made update to the solver it is available at the same address

Test cases for the hemispherical body with the Zwart and the Singhal models are still running.

Some notes on the implementation of volume fraction transport in OpenFOAM.
1)
It is formulated in volume fluxes, but it derived from mass fluxes:
ddt(alpha_l * rho_l) + div(alpha_l * rho_l * U) = m_l

ddt(alpha_l) + div(alpha_l U_l) = m_l/rho_l + (alpha_l/rho_l)*DDt(rho_l)

when fluid is incompressible and flow is isothermal, last term is zero, then we have:
ddt(alpha_l) + div(alpha_l U_l) = m_l

m_l - is the sum of m^+ and m^- of the Singhal model. In OpenFOAM m_l is represented by method mDotAlphal. That's why i think that we need to scale terms in Singhal model when switching from the mass fraction to the volume fraction

2)
In Singhal model pressure difference is calculated as: dP=pSat + c*k - p, but OpenFOAM interPhaseChangeFoam solver do not contain term "c*k" in the equation for pressure. And there is no such term in my model. That's why i don't include 0.5*0.39*rho_l*k when i estimate mass evaporation term.
If you want to insert this term in the model (interPhaseChangeFoam or compressibleMixingPhaseChangeFoam), then you need to change equation for pressure to obtain mass and momentum continuity

About NCG.
I think, i can add transport equations for the non-condensible gas in the model in three step:
1) Simple transport of the volume gas fraction without accounting for it's influence on the total mass of the system (alpha_g <= 1.0e-5 or another small value)
2) Transport of the volume gas fraction without desorption or absorption
3) And if two previous models will work, i can try to implement model with desorption and absorption

If you are still interested, it will be nice to run test cases, which you had proposed by e-mail

Hi,
It seems that the link to downloading the code is wrong. Is it possible to get an updated link for downloading?
Kind regards
Bjorn
bhh is offline   Reply With Quote

Old   January 15, 2016, 10:35
Default Link for downloading the code from Matvey Kraposhin
  #33
bhh
Member
 
Bjorn H. Hjertager
Join Date: Mar 2009
Posts: 71
Rep Power: 9
bhh is on a distinguished road
Hi,
It seems that the link for downloading this code is inactive. Does anybody in this thread have any infomation on how to download the code from Matvey Kraposhin mentioned in his post?
Kind regards
Bjørn
bhh is offline   Reply With Quote

Old   January 17, 2016, 21:29
Default Link for downloading the code from Matvey Kraposhin
  #34
New Member
 
popovich
Join Date: Jan 2016
Posts: 9
Rep Power: 2
popovich is on a distinguished road
Hi,
It seems that the link to downloading the code is wrong. Is it possible to get an updated link for downloading?
Kind regards
popovich
popovich is offline   Reply With Quote

Reply

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
Modelling Supersonic TwoPhase Flow with Phase Change wes OpenFOAM Running, Solving & CFD 8 April 26, 2016 07:21
Multiphase flow and Phase change due to heat transferevaporation haghajani OpenFOAM Running, Solving & CFD 86 April 9, 2013 02:49
phase change modeling Danial Q Main CFD Forum 0 April 5, 2012 01:14
low speed compressible flow lily CFX 2 November 16, 2005 06:15
Two phase flow models Atholl Main CFD Forum 2 May 7, 2002 03:49


All times are GMT -4. The time now is 17:27.