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

Energy Equation Solution Algorithm

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 23, 2006, 04:49
Default Energy Equation Solution Algorithm
  #1
Chandra
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 "non-steady" 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 = time-counter;

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

  Reply With Quote

Old   June 23, 2006, 05:03
Default Re: Energy Equation Solution Algorithm
  #2
O.
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))?

  Reply With Quote

Old   June 23, 2006, 06:11
Default Re: Energy Equation Solution Algorithm
  #3
Chandra
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 cell-center as below: u.dT/dx = 0.5*[u(i,j,k,n)+u(i-1,j,k,n)]*[T(i+1,j,k,n)-T(i-1,j,k,n)]/2dx. Here, n: time-step.

All the terms in the equation is discretized like the above on time-step '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 face-1 and a thin circulat slit near to the peripheral of the same face-1. A water jet is entering through the central hole, strikes on the opposite face-2 and the water goes out through the narrow-circular-slit. A heater having constant heat-flux is attached to the striking face-2. 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 central-hole, the temperature profile is showing a band of high and low temperature fields...like a pressure distribution in a wave-field.

If it is convinient to you, may I email the model and the result?

Thank you, Chandra
  Reply With Quote

Old   June 23, 2006, 06:25
Default Re: Energy Equation Solution Algorithm
  #4
O.
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 odd-even decoupling ;-)

Btw. What was your motivation for this:
0.5*[u(i,j,k,n)+u(i-1,j,k,n)]*[T(i+1,j,k,n)-T(i-1,j,k,n)]/2dx
rather than
u(i,j,k,n)*[T(i+1,j,k,n)-T(i-1,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(i-1,j,k,n)+u(i,j,k,n)]*[T(i,j,k,n)-T(i-1,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!
  Reply With Quote

Old   June 23, 2006, 07:09
Default Re: Energy Equation Solution Algorithm
  #5
Chandra
Guest
 
Posts: n/a
Hello, thank you very much for pointing out the mistake. I was knowing the odd-even 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(central-difference) the momemtum equations on correponsing "FACE-centers" rather than the "CELL-center"....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(i-1,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(i-1,j,k,n)*[T(i+1,j,k,n)-T(i-1,j,k,n)]/2dx ??

Motivation behind using: u.dT/dx = 0.5*[u(i,j,k,n)+u(i-1,j,k,n)]*[T(i+1,j,k,n)-T(i-1,j,k,n)]/2dx was that T(i,j,k) is defined on cell-center on staggerred grid and so discretize the temperature equation simply at the point-of-definition i.e. "cell-center" using central-difference. Point-of-definitions was used to discretize the momentum equations.

Thank you very much again for your help.

Chandra
  Reply With Quote

Old   June 23, 2006, 07:21
Default Re: Energy Equation Solution Algorithm
  #6
O.
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 ...

  Reply With Quote

Old   June 23, 2006, 08:05
Default Re: Energy Equation Solution Algorithm
  #7
Chandra
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

  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
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
k-e 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 high-Re flows? Zhong Lei Main CFD Forum 23 May 14, 1999 13:22


All times are GMT -4. The time now is 19:56.