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

transient term treatment

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 17, 2003, 22:54
Default transient term treatment
  #1
Mike
Guest
 
Posts: n/a
I have difficulty in the calculation of velocity and pressure fields for the transient flow. I am using the SIMPLE method to solve the velocity and pressure field for compressible flow. When the density is only function of pressure (isothermal) there is no problem, while when I solve the energy equation and the density is function of both pressure and temperature I get negative pressures. In a general form for the transient term in the countinuity I write:

d/dt (rho)=d/dP (rho) . dP/dt+ d/dT (rho). dT/dt

Now for ideal gas we have: rho=P/(kT)

d/dP (rho)= 1/(kT) (k is molar gas constant)

d/dT (rho)= -P/(kT^2)

The flow is laminar and the mach number is very small, however the density variation is big. Is there any trick to handle the transient part in the pressure equation?

Thanks for the help.
  Reply With Quote

Old   February 17, 2003, 23:10
Default Re: transient term treatment
  #2
mukhopadhyay
Guest
 
Posts: n/a
although i am not fully clear about your problem, it appears that the source term needs attention - i mean the linearisation.
  Reply With Quote

Old   February 18, 2003, 00:04
Default Re: transient term treatment
  #3
Mike
Guest
 
Posts: n/a
How should I do that?
  Reply With Quote

Old   February 18, 2003, 02:24
Default Re: transient term treatment
  #4
mukhopadhyay
Guest
 
Posts: n/a
suggest consult the book by Patankar.
  Reply With Quote

Old   February 18, 2003, 03:11
Default Re: transient term treatment
  #5
Mike
Guest
 
Posts: n/a
Thanks. I guess I have already done that. for example as I said I have written:

d/dt (rho)= d/dP (rho) .dP/dt + d/dT (rho) . dT/dt

thus, for the ideal gas we get:

d/dt (rho)= 1/(kT).dP/dt- P/(kT^2) . dT/dt

Now I have no problem with term 1/(kT) dP/dt, but the problem apears to be with P/(kT^2) .dT/dt. I can either used it in the central coefficient or put it in the source term using the perivious pressure. In either case it gives me negative pressure and wrong velocities!

Can you give a specific solution? Thanks again
  Reply With Quote

Old   February 18, 2003, 03:39
Default Re: transient term treatment
  #6
Rami
Guest
 
Posts: n/a
Hi Mike,

Have a look in the following paper:

K.C. Karki and S.V. Patankar, Pressure Based Calculation Procedure for Viscous Flows at All Speeds in Arbitrary Configurations, AIAA J V27 N9, 1989, pp 1167-1174.

Although steady-state is treated, it might be helpful.

Rami
  Reply With Quote

Old   February 18, 2003, 07:47
Default Re: transient term treatment
  #7
andy
Guest
 
Posts: n/a
At low Mach numbers the term in the momentum equation you are having problems with makes the momentum equation very stiff. There are several ways of handling variable density low Mach number flows and I would suggest you perform (or look up) a low Mach number asymptotic analysis of your set of equations. Subsequent to this you can differentiate the resulting equation of state in order to derive a reasonably well behaved RHS for whichever form of the pressure/continuity equation you adopt.
  Reply With Quote

Old   February 18, 2003, 10:32
Default Re: transient term treatment
  #8
gorka
Guest
 
Posts: n/a
Andy:

Your suggestion seems very interesting. What do you mean with RHS? Do you have any reference (Paper or book) regarding this issue?

Many thanks.
  Reply With Quote

Old   February 18, 2003, 11:22
Default Re: transient term treatment
  #9
andy
Guest
 
Posts: n/a
By RHS I mean Right Hand Side or the source term of the elliptic equation which is solved to enforce continuity in a typical low speed code.

I suspect most texts describing the details of a split/segrated/projection/approximate factorization type numerical method for low speed variable density flow is likely to discuss or refer to the point. I guess combustion is probably the most widely researched low speed variable density flow and I would suggest starting with a few papers in this area.

  Reply With Quote

Old   February 18, 2003, 14:26
Default Re: transient term treatment
  #10
Mike
Guest
 
Posts: n/a
Andy,

Thanks a lot. As I said it seems the main problem appears from the countinuity transient term. I will try to find some papers in the area you said. Do you have any specific paper in mind?
  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
transient source term strakakl OpenFOAM 38 November 19, 2013 02:18
pressure transient term? Romek CFX 4 April 19, 2011 21:19
How to define a transient species source term? aleisia FLUENT 3 March 7, 2011 23:58
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 15:45
[Q] Convection term treatment in UDS Ryan, Lee FLUENT 4 October 18, 2004 10:20


All times are GMT -4. The time now is 13:55.