Stabilizing Transient Source term
 

I am solving a energy transport equation of the form

drhoH/dt + del.(rho V H) = -dP/dt+...+.

in a collocated arrangement. I see that as the simulation progresses, the sources due to dp/dt is highly unstable.

How can I stabilize my computation? The computation is for an incompressible flow but with heat transfer.

--

CFDtoy

Re: Stabilizing Transient Source term
 

What kind of time-marching scheme do you use?.For multi-stage R-K it might be unclear how to discretize dp/dt on RHS.

Re: Stabilizing Transient Source term
 

I am using a simple backward euler (implicit for stability!) Does any time splitting work?

CFDtoy

Re: Stabilizing Transient Source term
 

The stability analysis is based on linearized theories and hence you cannot necessarily guarantee stability for non linear problems.Did you check if explicit time marching worked for this particular case?

Re: Stabilizing Transient Source term
 

explicit time stepping posed a lot of problems..since i hd problems in converging them.

normally when the source terms are not functions of the flw variables and cause such kind of stability problem..is there a way to reduce their propagation !>!? I am not aware of time splitting schmes for these purposes?

anythin related?

Thanks for your cmments harish.

-- CFDtoy

Re: Stabilizing Transient Source term
 

I have never come across anything either.Is there any reason for using the equation for enthalpy instead of total energy or internal energy.Those equations do not have a time dependent pressure term and your implicit euler implementation might be easier.

Re: Stabilizing Transient Source term
 

In case of total enthalpy, you do require the dp/dt term on the RHS. Also, finding temp using total enthalpy seemed easier.

CFDtoy

Re: Stabilizing Transient Source term
 

I think it would be the same even if you use the equation for internal energy.