
[Sponsors] 
June 23, 2006, 04:49 
Energy Equation Solution Algorithm

#1 
Guest
Posts: n/a

Hi everybody,
I have solved the N.S.Momentum equations and now i have to solve the energy equation for a "nonsteady" flow problem i.e. del_T/del_t + V.Grad_T = Alpha*Div_(Grad_T). I used explicit method: T(n+1) = T(n) + dt*[Alpha*Div_(Grad_T)  V.Grad_T](n); Here, n = timecounter; But, the result doesnt seem to be appropriate. I think something is wrong in the algorithm itself e.g. such kind of explicit scheme cant be used. Please suggest me something in this regard. I think the people who already have solved the CFD problems before would be havign much idea about it. Thanks a lot in anticipation, regards, Chandra 

June 23, 2006, 05:03 
Re: Energy Equation Solution Algorithm

#2 
Guest
Posts: n/a

A simple explicit scheme should work.
How did you discretize the convective part? Do you have any coupling between your energy equation (which is a simple scalar convection/diffusion equation in this case) and the momentum and continuity equations (e.g. density = f(T))? 

June 23, 2006, 06:11 
Re: Energy Equation Solution Algorithm

#3 
Guest
Posts: n/a

Hello, Thanks for your reply. I dont have any coupling. Fluid properties are assummed to be independent of temperature. Also, incompressible flow (water) and staggerred grid.
I discretized the convective part at cellcenter as below: u.dT/dx = 0.5*[u(i,j,k,n)+u(i1,j,k,n)]*[T(i+1,j,k,n)T(i1,j,k,n)]/2dx. Here, n: timestep. All the terms in the equation is discretized like the above on timestep 'n'. And T(n+1) is evaluated simply at all the cells. The problem: there is a cylinder(two faces and one round wall) with a small hole at face1 and a thin circulat slit near to the peripheral of the same face1. A water jet is entering through the central hole, strikes on the opposite face2 and the water goes out through the narrowcircularslit. A heater having constant heatflux is attached to the striking face2. I solve the equations (continuity+momentum & energy). Velocity profile is coming fine. But for temperature profile, when the initial developing vortex strikes the surface after starting from the centralhole, the temperature profile is showing a band of high and low temperature fields...like a pressure distribution in a wavefield. If it is convinient to you, may I email the model and the result? Thank you, Chandra 

June 23, 2006, 06:25 
Re: Energy Equation Solution Algorithm

#4 
Guest
Posts: n/a

You have a central discretization. This is not stable without artificial dissipation. Try to upwind the temperature. For the velocity you can try both, upwind or central.
I think you are having a classic oddeven decoupling ;) Btw. What was your motivation for this: 0.5*[u(i,j,k,n)+u(i1,j,k,n)]*[T(i+1,j,k,n)T(i1,j,k,n)]/2dx rather than u(i,j,k,n)*[T(i+1,j,k,n)T(i1,j,k,n)]/2dx ?? I would suggest this (first order only for a start) u(i,j,k,n) >= 0 : u.dT/dx = 0.5*[u(i1,j,k,n)+u(i,j,k,n)]*[T(i,j,k,n)T(i1,j,k,n)]/dx u(i,j,k,n) < 0 : u.dT/dx = 0.5*[u(i,j,k,n)+u(i+1,j,k,n)]*[T(i+1,j,k,n)T(i,j,k,n)]/dx Good luck! 

June 23, 2006, 07:09 
Re: Energy Equation Solution Algorithm

#5 
Guest
Posts: n/a

Hello, thank you very much for pointing out the mistake. I was knowing the oddeven coupling in case of momentum equations but didnt know about the temperature equation. Thanks! I had read in some of the papers that if we are using staggerred grid, we actually dont need to employ upwinding specifically, rather, staggerred grid takes care of itself once we discretize(centraldifference) the momemtum equations on correponsing "FACEcenters" rather than the "CELLcenter"....and that is what i did while solving the momentum equations.
I have a simple question, however. You suggested: u.dT/dx = u(i,j,k,n)*[T(i+1,j,k,n)T(i1,j,k,n)]/2dx. Is this for "u(i,j,k,n) > 0" ? and for "u(i,j,k,n) < 0", u.dT/dx = u(i1,j,k,n)*[T(i+1,j,k,n)T(i1,j,k,n)]/2dx ?? Motivation behind using: u.dT/dx = 0.5*[u(i,j,k,n)+u(i1,j,k,n)]*[T(i+1,j,k,n)T(i1,j,k,n)]/2dx was that T(i,j,k) is defined on cellcenter on staggerred grid and so discretize the temperature equation simply at the pointofdefinition i.e. "cellcenter" using centraldifference. Pointofdefinitions was used to discretize the momentum equations. Thank you very much again for your help. Chandra 

June 23, 2006, 07:21 
Re: Energy Equation Solution Algorithm

#6 
Guest
Posts: n/a

You don't really have a staggered grid for the temperature, hence you get the stability problem. I would go ahead with an upwind for that. Since you already have the velocities on the faces, use them and upwind the temperature. For second order you can use a MUSCL approach (extrapolate the values for T from the upwind direction). You might or might not need a limiter ...


June 23, 2006, 08:05 
Re: Energy Equation Solution Algorithm

#7 
Guest
Posts: n/a

Hello,
Thank you very much for your help. I got very useful information from you and it is going to very helpful indeed. Thank you, Chandra 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Source term energy equation for reactive flows  DaIN  Main CFD Forum  0  October 6, 2011 15:11 
energy equation in rhoCentralFoam  nakul  OpenFOAM  0  October 10, 2010 15:07 
ke turbulence model and energy equation  Blob  Main CFD Forum  0  May 29, 2009 08:35 
question about energy equation  zhou  FLUENT  0  February 24, 2004 00:55 
Why FVM for highRe flows?  Zhong Lei  Main CFD Forum  23  May 14, 1999 13:22 