CFD Online Logo CFD Online URL
Home > Forums > FLUENT

why enthalpy residual is great ???

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

LinkBack Thread Tools Display Modes
Old   January 9, 2007, 12:43
Default why enthalpy residual is great ???
Posts: n/a
Hi CFD specialists ;

I am working on two phase flow with mass transfer process;I need to check following expression for computing residual of energy equation without energy source term in steady state. ( div.[ rho U (E+p/rho)] -div.[k (grad (T))]=0.0 (k is effective conductivity)

Fluent uses finite volume method for discretization of left hand side of upper equation according to the following form :

S= [(sigma on faces) {rho U.dA (E+p/rho)}-(sigma on faces ) { k grad(T).dA}]/ cell volume

In single phase flow without energy source term , S should be zero. I use from a udf in order to certify this object , Udf computes right hand side of upper equation ( S); The part of udf which calculates S is presented at following:

thread_loop_c (t,domain)




int i=1;








convection=(F_UDMI(f,tf,1)+F_UDMI(f,tf,2)/C_R(c,t))*F_FLUX(f,tf) ;

diffusion=C_UDMI(c,t,11)*A_by_es/ds*(F_T(f,tf)-C_T(c,t)); } else{

INTERIOR_FACE_GEOMETRY(f,tf,A,ds,es,a_by_es,dr0,dr 1);

if(F_C0(f,tf)==c) {



diffusion=0.5*(C_UDMI(c,t,11)+C_UDMI(F_C1(f,tf),t, 11))*A_by_es/ds*(C_T(F_C1(f,tf),t)-C_T(c,t));

if(F_FLUX(f,tf)>0) {

convection=C_UDMI(c,t,12)*F_FLUX(f,tf); /*method is first order upwind*/

else {


} }

else {



diffusion=0.5*(C_UDMI(c,t,11)+C_UDMI(F_C0(f,tf),t, 11))*A_by_es/ds*(C_T(F_C0(f,tf),t)-C_T(c,t));

if(F_FLUX(f,tf)<0) {


} else {


} } }

C_UDMI(c,t,4)=C_UDMI(c,t,4)-diffusion+convection ;


C_UDMI(c,t,4)=C_UDMI(c,t,4)/C_VOLUME(c,t) ;


end_c_loop(c,t) }

C_UDMI(c,t,11) indicates to effective conductivity in last iterration;

C_UDMI(c,t,12) indicates to [e+P/rho] in last iteration;

F_UDMI(f,tf,1)indicates to enthalpy in face;

F_UDMI(f,tf,2)indicates to pressure in face;

C_UDMI(c,t,4) is equal to S in context .

But in any iteration I faced to S in range of -250000 to 350000 ; whereas in any iteration S should be zero . Why it happended ?

Have you a suggestion in order that I could remedy this problem ? Does fluent use another method for calculating of S ?

I am appreciated for everybody who help me in solving this problem.



  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
transsonic nozzle with rhoSimpleFoam Unseen OpenFOAM Running, Solving & CFD 7 April 16, 2014 03:38
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM 6 April 12, 2011 11:24
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16

All times are GMT -4. The time now is 20:00.