CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Verification of this phase change modeling

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 15, 2011, 02:20
Default Verification of this phase change modeling
  #1
New Member
 
kawamura
Join Date: May 2011
Location: Ishikawa,Japan
Posts: 7
Rep Power: 14
kawamura is on a distinguished road
Hi everyone.

I'm trying to two phase flow with phase change using VOF method.
There are some approach for mass transfer modeling.
I'm interested in the following way.
mDot = K/L * grad(T)・grad(F)
Note that K refers to the mixture thermal conductivity, F refers to the volume fraction and L is the latent heat.

But I couldn't find that verification.
If you understand this correctry, please let me know the verification!
Thanks a lot.
kawamura is offline   Reply With Quote

Old   December 20, 2011, 02:54
Default
  #2
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 74
Rep Power: 16
Ohbuchi is on a distinguished road
Hi, Kawamura-san,

I can not make a verification for this notation.
But I think grad(F) should be grad(F)/mag(grad(F)).
This means unit normal vector of interface.
Ohbuchi is offline   Reply With Quote

Old   December 20, 2011, 04:46
Default
  #3
New Member
 
kawamura
Join Date: May 2011
Location: Ishikawa,Japan
Posts: 7
Rep Power: 14
kawamura is on a distinguished road
Hi Ohbuchi-san,

thanks for your response.
but FLUENT and some people show this expression in their paper .
there is non-existence of normal vector.

I think the most accurate expression is...
mDot = 1/L *(qg-ql)・n *Aint/V (kg/m^3 /s) in evaporating
where qg and ql is heat flux in gas and liquid phase respectively,Aint is interface area,V is volume of cell and n is unit normal vector.

but getting the correct interface area is so difficult in VOF.
so I think that grad(F) may works as n・Aint/V .
but i don't get a correct understanding about this
kawamura is offline   Reply With Quote

Old   December 21, 2011, 00:59
Default
  #4
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 74
Rep Power: 16
Ohbuchi is on a distinguished road
Hi. Kawamura-san,
The value of grad(F) will be very large on sharp interface, so I think some limiter function will be needed to restrict the value within 0-1.
And, I think it's impossible to determine interface area from grad(F) without assuming specific flow regime(such as stratified flow or bubbly flow).


BTW, the expression of mas transfer by heat transfer;
mDot = (k/L) grad(T) grad(F)/mag(grad(F)),
was originally introduced in following paper.

Son,G. and Dhir,V.K. A level set method for analysis of film boiling on an immersed solid surface. Nummer. Heat Transfer B,vol.52,pp.153-177.(2007)

This expression is only valid on stratified flow regime.
Ohbuchi is offline   Reply With Quote

Old   December 21, 2011, 01:14
Default
  #5
New Member
 
kawamura
Join Date: May 2011
Location: Ishikawa,Japan
Posts: 7
Rep Power: 14
kawamura is on a distinguished road
Hi Ohbuchi-san,

It's useful information for me!
Now it's starting to make sense.

I consider about assuming specific flow.

Thank you for your kindness.

Kawamura
kawamura is offline   Reply With Quote

Reply

Tags
phase change, vof

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
UDF of "anisotropic phase function" in radiation modeling whiz Fluent UDF and Scheme Programming 0 September 1, 2009 04:26
thermal phase change question CFDflying CFX 1 February 18, 2004 05:10
compressible two phase flow in CFX4.4 youngan CFX 0 July 2, 2003 00:32
phase change rookie FLUENT 3 June 12, 2003 21:47
Phase change modelling Sebastien DANIEL FLUENT 1 July 6, 2000 04:35


All times are GMT -4. The time now is 07:35.