i'm simulating heat and mass transfer inside two phase closed thermosyphon using 2-d axysymmetric geometry.i've selected pressure based solver and transient.I am using ANSYS /FLUENT 19 student version .VOF method with Lee Model is used . Water is used as secondary phase and Vapor is used as primary phase. SIMPLE algorithm scheme for pressure velocity coupling and a first order upwind scheme for the determination of momentum and energy is included in the model. Geo –Reconstruct and PRESTO discretization for the volume fraction and pressure interpolation scheme, respectively, are also performed in these simulations. but,after running the calculation,contour is showing only in the evaporator region as shown in attached figure.I've also patched the liquid water fill ratio with volume fraction=1 in the evaporator region.Can anyone help me out?

# simulation of two phase heat and mass transfer inside two phase closed thermosyphon

- 240 Views
- Last Post 20 November 2018

First of all: Please avoid generating Entropy by creating several topics regarding the same questions.

1/2D Axis Symmetric is only for the cases which are axis-symmetric. Can you confirm that from your side?

2/The contour shows only that some vapor has been generated in the evaporation section. Try to clip the contour to some lower values to detected condensation regions. You can even check the Mass transfer under Phase Interaction for post-processing.

3/With the Lee Model you have coefficients (frequency) and those would require a high amount of tuning and this tuning won't make any sense without exp. data or references. Generally the condensation frequency is set to lower value than the one of evaporation but it is case dependent.

Best regards,

Amine

sir.i am simulating the work done by B.Fadhl https://www.sciencedirect.com/science/article/pii/S1359431113004699.i've selected an axisymmetric domain as attached .The line at x=0 is the centerline.i've further splitted the right wall into three parts namely evaporator,adiabatic and condenser in ANSYS SPACECLAIM 19 student version.Frequency is 0.1 in Lee model as given in the article attached. The problem that i'm facing is that Contour shows the phase change only in the evaporator section.

I have said a couple of times: we do not have access to all journals across the world. From the abstract: the authors are using an UDF.

I wanted just to let you check if the case if rotational about the x-axis. If this is not fulfilled either you go to planar 2D or fully 3D.

You can clip the contour lines to something smaller to detect the condensation zones. Please check too the phase interaction post-processing variable.

Best regards,

Amine

it is rotional about the y-axis. Thermosyphon is cylindrical in shape.So,can i go for planar 2D.

No you do not need to do that. You can stick to axis-symmetric but we need to rotate it. So the "Y-direction" of the model will be your X-Axis. You can then take only the half of the diameter and model the section to be aligned with the X-Axis.

Best regards,

Amine

Hello Prashant,

Just to add to Amine's comments:

- Please read the paper and make sure you are using the same set of source terms as the authors. If I am to speculate, they are using sources terms responsible of evaporation and condensation in the mass and energy equations.
- In addition to this, please make sure you are using consistent initial conditions, similar to what the authors have done.
- Another important factor is the mesh. If there is a way to verify the mesh, please make sure you are using the same resolution of mesh.

Replicating the results from a paper is always a challenging task. The most important suggestion I would like to give you is this - when in doubt about the exact steps used by the authors, you can always contact the authors of the paper via email. This is important as they might be able to clarify the exact model for you. Please let us know and we might be able to help you with questions regarding best usage and practices of the modeling tools.

Thank you.

Best Regards,

Karthik

So,is there any problem having vertical direction aligned with y-axis not with x-axis?

In Fluent if you want to use the 2D Axis-symmetric solver then the symmetry/ rotational axis needs to be the x-axis. You can then just apply gravity in the -X direction.

Best regards,

Amine

sir, i that case,what will be the sign of g.Will it be g=9.8m/s^2?

Example if the height goes from 0.0 to +HMAX in the positive x-direction. Then g will be -9.81 m/s^2.

Best regards,

Amine

in ANSYS FLUENT19.0.Is is important to add source term UDF,if i have already selected the phase interaction as Lee model for evaporation-condensation in Fluent.Since,there is no momentum source term in the VOF model,when selected Lee model,will Ansys Fluent will take care of the source terms?

All the stuff described in the paper is already implemented in FLUENT. You do not need to add any UDF. All required secondary source terms are accounted for.

Best regards,

Amine

for the energy source term,since latent heat of vaporization is 2455 kj/kg.I've incorporated the following udf for energy source term.Please,check it sir.

#include "udf.h"

#include "sg_mphase.h"

#define T_SAT 373

#define LAT_HT 2455.1345e3

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)

{

Thread *pri_th, *sec_th;

real m_dot;

pri_th = THREAD_SUB_THREAD(mix_th, 0);

sec_th = THREAD_SUB_THREAD(mix_th, 1);

if(C_T(cell, mix_th)>=T_SAT)

{

m_dot = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;

dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*LAT_HT/T_SAT;

}

if(C_T(cell, mix_th)<=T_SAT)

{

m_dot = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(T_SAT-C_T(cell,pri_th))/T_SAT;

dS[eqn] = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*LAT_HT/T_SAT;

}

return LAT_HT*m_dot; }

since latent heat of vaporisation is 42455kj/kg.I'm Using the attached udf for energy source

sorry,latent heat=2455kj/kg

You do not need require any UDF. Furthermore if UDF then I will go for DEFINE_LINEARIZED_MASS_TRANSFER.

Best regards,

Amine

sir, in paper they have considered latent heat of vaporisation as 2455kj/kg,how to incorporate that value in ANSYS FLUENT 19 STUDENT VERSION.

You use your working boundary conditions. Regarding mass transfer just use the implemented model and just provide the saturation temperature and change the frequencies.

Latent Heat in Fluent is in J/kgmol and you have the molar weight so you can incoprarte theier setting in Fluent.

You need to carry out some tutorials in order to get really started...

Best regards,

Amine

sir,i've followed all your instructions ,simulation is started in 2-D domain but the coontinuity equation is diverging and the error appears.I've attached two pics

Same as usual: smaller time steps, lower URF's.. etc.

If it does not help then perhaps other community members might help you because as ANSYS Stuff I am not allowed to give more that I have already done.

Best regards,

Amine

Hello Prashant,

first of all,

Gram atomic mass of water (H2O) = (1×2+16) gm = (2+16) gm = 18 gm

Number of molecules in 18 gm of water = 1 mol

Number of the molecule in 1 gm of water = 1/18 mol

Number of molecule in 3.6 gm of water = (1/18 x 3.6)mol = 0.2 mol

Now go to report and get the mass of water, please refer below fig.

in my case, it is approx. 2.5 kg (initially)

Number of molecule in 2500 gm of water = (1/18 x 2500)mol = 138.89 mol

latent heat 2455 KJ/Kg = 2455*e3 J/Kg

the standard latent heat of enthalpy = (2455*E3/138.89) = 17675.85 J/(Kg.mol)

Or just take the molar weight in kg/kmol and multiply by the standard state enthalpy in j/kg.

Best regards,

Amine

Hello, I`m working on the same model. Can we discuss about it if you want.

Thank you

##### Search

##### Change Language

##### Categories

##### This Weeks High Earners

- abenhadj 133
- peteroznewman 132
- rwoolhou 127
- kkanade 89
- jj77 47
- FlaviaGelatti 44
- tsiriaks 44
- Max4 27
- EngineerMech 24
- SandeepMedikonda 14