CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Problem with phase change - Relating to TEMPERATURE

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

LinkBack Thread Tools Search this Thread Display Modes
Old   August 5, 2019, 01:12
Default Problem with phase change - Relating to TEMPERATURE
Senior Member
Join Date: Oct 2016
Posts: 349
Rep Power: 11
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Hello Foamers, In my phase change problem, I want to set the liquid evaporation equal to interfacial temperature (saturation temperature).
My conditions at the interface (boiling plane) are as follows:
T= 100 (saturation temperature when phase change occurs)
@T=100 --> [w(moisture content)->at high (-) temperature side =0]

Based on it, mass and heat balances are calculated at interface as,
Mass balance (r) = rho * alpha * (dx/dt) [kg/m2.s] (calculated at -> low (+) temperature side)
Heat balance (-K dT/dx) = (r) * latent heat [W/m2]
forAll(mesh.cells(), celli) {
if (alpha[celli] <= scalar(0.1)) //10% of moisture
T[celli] = 373.0; //Tsat = 100deg - T[celli] = Tsat.value();
alpha[celli]= scalar(0); //Moisture content = 0
volScalarField VX = V.component(vector::X);
// Compute mass source field (mass balance) (rho*alpha1*interface_VELOCITY)
volScalarField mDot = -(rho*alpha*VX);
// Compute energy source field using enthalpy of evaporation (energy balance)
//hDot = -evapEnthalpy*mDot;
{ fvScalarMatrix TEqn
fvm::ddt(rho*cp, T) // fvm::ddt(cp, T)
- mDot * cp * gradTx
- fvm::laplacian(K, T) // - fvm::laplacian(lambda/rho, T)
Hereby, attaching the solver. and attached figure gives the equations and explanations about boiling and convection.
ENERGY EQU: rho*Cp*(delT/delt) = del/delx[K*delT/delx ] + r*cp*delT/delx + rho*delQ/delT (W/m3)
REF - Clear details here in this MANUSCRIPT ~ I

I have partially set the condition of boiling in the solver, However, compilation is success - but there found no changes in results of alpha and temperature.
Some assistance will be highly helpful. Thank you
Attached Images
File Type: jpg ATTACHMENT_1_Boiling.jpg (110.9 KB, 13 views)
File Type: png ATTACHMENT_2_Convection.png (153.7 KB, 10 views)
Attached Files
File Type: gz myboilingFoam_solver.tar.gz (112.8 KB, 8 views)
Kummi is offline   Reply With Quote


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
Phase Change >> which solver ? >> met a problem I can't solve, please help me with it Kummi OpenFOAM Programming & Development 6 September 8, 2021 04:58
Phase Change Problem karthikeyas Main CFD Forum 1 January 30, 2018 01:47
Fluent solving phase change problem krasus FLUENT 0 May 23, 2017 15:30
Material change during phase transformation (FePO4 > LiFePO4) concentration distribut kimsw0103 CFX 0 January 15, 2016 22:05
Geometry setup for phase change problem Fred G. Kang Main CFD Forum 1 October 14, 1998 11:41

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