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

evaporation of R134a in a channel

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

Like Tree4Likes
  • 2 Post By peterliu
  • 2 Post By peterliu

Reply
 
LinkBack Thread Tools Display Modes
Old   May 19, 2011, 10:48
Default evaporation of R134a in a channel
  #1
New Member
 
qingming.liu
Join Date: May 2011
Posts: 4
Blog Entries: 1
Rep Power: 5
peterliu is on a distinguished road
Hi
I am doing simulation of evaporation of R134a in a channel and I did not observed phase change from R 134a to R134a gas. Does anyone can help me ?

FLUENT
Version: 3d, dp, pbns, vof, lam, transient (3d, double precision, pressure-based, VOF, laminar, transient)
Release: 12.1.4
Models
------

Model Settings
---------------------------------------------------------
Space 3D
Time Unsteady, 1st-Order Implicit
Viscous Laminar
Heat Transfer Enabled


Material Properties
-------------------

Material: r134a-gas (fluid)

Property Units Method Value(s)
----------------------------------------------------------------
Density kg/m3 constant 32.203999
Cp (Specific Heat) j/kg-k constant 1030.98
Thermal Conductivity w/m-k constant 0.01381
Viscosity kg/m-s constant 1.1687e-05
Molecular Weight kg/kgmol constant 102.03
Standard State Enthalpy j/kgmol constant 42063000
Reference Temperature k constant 298.14999

Material: r134a-liquid (fluid)

Property Units Method Value(s)
-----------------------------------------------------------------
Density kg/m3 constant 1207.3
Cp (Specific Heat) j/kg-k constant 1424.41
Thermal Conductivity w/m-k constant 0.081198998
Viscosity kg/m-s constant 0.00019525
Molecular Weight kg/kgmol constant 102.03
Standard State Enthalpy j/kgmol constant 23909000
Reference Temperature k constant 298.14999
Thermal Expansion Coefficient 1/k constant 0
Speed of Sound m/s none #f

Boundary Conditions
-------------------

Zones

name id type
------------------------------------
wall-solid 10009 wall
inlet 6 velocity-inlet
outlet 7 pressure-outlet

wall-solid

Condition Value
-------------------------------------------------------------
Thermal BC Type 1
Temperature (k) 300
Heat Flux (w/m2) 15000
Convective Heat Transfer Coefficient (w/m2-k) 0
Free Stream Temperature (k) 300


inlet

Condition Value
--------------------------------------------
Velocity Magnitude (m/s) 0.2
Coordinate System 0
Temperature (k) 295.15

Solver Settings
---------------
---------------------------------------------
Time Step (s) 0.00050000002
Max. Iterations Per Time Step 100

Relaxation

Variable Relaxation Factor
-----------------------------------
Pressure 0.3
Density 1
Body Forces 1
Momentum 0.7
Volume Fraction 0.5
Energy 1

Pressure-Velocity Coupling

Parameter Value
------------------
Type SIMPLE

Discretization Scheme

Variable Scheme
------------------------------------
Pressure PRESTO!
Momentum First Order Upwind
Volume Fraction First Order Upwind
Energy First Order Upwind
Attached Images
File Type: jpg flowboiling_out_vel.jpg (70.1 KB, 33 views)
File Type: jpg flowboiling_out_tem.jpg (67.1 KB, 28 views)
File Type: jpg flowboiling_out_phase.jpg (54.4 KB, 25 views)
jin5370 and prakashpvl like this.
peterliu is offline   Reply With Quote

Old   May 19, 2011, 10:53
Default
  #2
New Member
 
qingming.liu
Join Date: May 2011
Posts: 4
Blog Entries: 1
Rep Power: 5
peterliu is on a distinguished road
Forgot to mentioned the UDF, I setted the saturated temperature to 298.15K for R 134a in 5.6bar pressure.

/* UDF to define a simple mass transfer based on Saturation
Temperature. The "from" phase is the gas and the "to" phase is the
liquid phase */

#include "udf.h"

DEFINE_MASS_TRANSFER(liq_gas_source, cell, thread, from_index,
from_species_index, to_index, to_species_index)
{
real m_lg;
real T_SAT = 298.15;
Thread *gas = THREAD_SUB_THREAD(thread, from_index);
Thread *liq = THREAD_SUB_THREAD(thread, to_index);

m_lg = 0.;
if (C_T(cell, liq) >= T_SAT)
{
m_lg = -0.1*C_VOF(cell,liq)*C_R(cell,liq)*
fabs(C_T(cell,liq)-T_SAT)/T_SAT;
}
if ((m_lg == 0. ) && (C_T(cell, gas) <= T_SAT))
{
m_lg = 0.1*C_VOF(cell,gas)*C_R(cell,gas)*
fabs(T_SAT-C_T(cell,gas))/T_SAT;
}

return (m_lg);
}
ddmdddghj and prakashpvl like this.
peterliu is offline   Reply With Quote

Old   August 2, 2012, 08:37
Default evaporation of R134a in a channel
  #3
New Member
 
Vinod
Join Date: Aug 2012
Posts: 1
Rep Power: 0
vinod500ind is on a distinguished road
Hi,
I'm having similar kind of problem. That's unable to observe the phase change.

Did you found any solution on this?

Please do let me know.

Regards,

Vinod
vinod500ind is offline   Reply With Quote

Old   July 8, 2013, 01:27
Default
  #4
New Member
 
mahdad
Join Date: Jun 2013
Posts: 2
Rep Power: 0
ma_ta65 is on a distinguished road
hi
i want the standard state enthalpy of r404a(liquid and gas) for evaporation modeling.
ma_ta65 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
Open Channel Flow ElanMorin FLUENT 3 August 19, 2011 14:21
Open channel flow motaba Main CFD Forum 4 March 26, 2011 04:22
tecplot 3D velocity contours inside a channel vetnav Tecplot 4 July 14, 2010 20:03
About Cooling channel alefem FLOW-3D 1 May 28, 2010 11:15
evaporation of liquid water in a horizonal channel Ahmad Al-Zoubi CFX 2 February 2, 2007 04:55


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