CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Modeling Evaporation at a Heated Wall

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

LinkBack Thread Tools Display Modes
Old   February 10, 2011, 10:20
Arrow Modeling Evaporation at a Heated Wall
Join Date: Sep 2010
Posts: 35
Rep Power: 7
Edy is on a distinguished road
Hi Foamers,

I want to model nucleate boiling in cylindrical pipes, and one important part in nucleate boiling is ... boiling at the wall!

I use the Kurul & Podowski heat flux partitioning. In my solver I calculate the three contributions. The evaporation heat flux density q_e is defined as volScalarField which is null everywhere except for near wall cells where:

q_e=Pi/6*d_dep*freq*NSD*rho_v*h_lv (where d_dep is the departure diameter of bubbles,freq the detachment frequency, NSD the nucleation site density and h_lv is the latent heat).

Once this is known I can determine the evaporation mass flow rate at the wall, in [kg/m/s], which is defined as:

m_e = q_e/h_lv

But all the equations solved in my solver (alphaEqn, UEqns, ...) are dealing with cells, meaning finite volumes, therefore I need a volume evaporation rate, whose unit is in [kg/m/s]. I thought about calculating it as :

Gamma_e = m_e*(area of the cell in contact with the heated wall)/(volume of the cell)

But if x is the direction of the flow:

Gamma_e = m_e*(x_cell*z_cell)/(x_cell*y_cell*z_cell)
= m_e/y_cell

And my evaporation rate depends on the mesh refinement which is of course non physical...

Has anyone an idea how to solve this problem? Cause I've been thinking about it for a couple of days unsuccessfully...

Thanks for reading, giving a tip or helping!


Edy is offline   Reply With Quote

Old   February 10, 2011, 10:28
Join Date: Sep 2010
Posts: 35
Rep Power: 7
Edy is on a distinguished road
I wrote "mass flow rate" for m_e but it is indeed rather a evaporation mass flux!
Edy is offline   Reply With Quote

Old   February 11, 2011, 05:32
Senior Member
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,147
Blog Entries: 1
Rep Power: 15
nimasam is on a distinguished road
in openFoam there is no interface reconstruction so you can not convert mass flux into mass evaporation per volume! like :
massflux * interfaceArea/ volume

you can only use mag(grad(alpha)) which seems its not compeletly true for local interface area estimation
how ever i eager in this subject, can u share ur efforts (aricles and source code)?
nimasam65 <@>
nimasam is offline   Reply With Quote


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
evaporation modeling - have problem with temperature field MOhaj Main CFD Forum 0 April 25, 2010 00:21
problem with temperature field in modeling evaporation with fluent MOhaj FLUENT 0 April 24, 2010 23:54
evaporation of drop on a heated plate ztdep Main CFD Forum 0 March 23, 2010 23:08
Question about bcdefw.f for wall temperature bc. Jimmy CD-adapco 10 March 18, 2008 16:28
Convergence difficulty in modeling evaporation Lawrence J CFX 0 November 30, 2004 10:08

All times are GMT -4. The time now is 23:19.