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

A question about Peclet number

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 8, 1999, 07:38
Default A question about Peclet number
  #1
Zhong Lei
Guest
 
Posts: n/a
I have a simple problem which I have not understood.

It is well known that the Peclet number ( Peclet = u dx / \nu) must be smaller than 2.0 to maintain numerical stability in a linear convective-diffusive problem, where the convective velocity remains constant. However for a nonlinear problem such as boundary layer flows, as the flow goes away from the wall, the local velocity is increased. To satisfy the condition of Peclet number < 2, does it mean that the spatial size dx should be decreased? Certaninly, it does not. But why?
  Reply With Quote

Old   March 8, 1999, 10:05
Default Re: A question about Peclet number
  #2
John C. Chien
Guest
 
Posts: n/a
Remember that at each point there are three velocity components, so, there are three such numbers. In a more general convection-diffusion equation, they are related to the coefficient of the first-order differential terms. The coefficient does not have to be the velocity components. But there are always three first-order terms. The numbers associated with each direction ( coordinates) are not the same. For 1-D, there is one cell Reynolds number, for 2-D, there are two cell Reynolds number, and for 3-D, there are three cell Reynolds number,etc... A good exercise is : try to solve a standard 2-D cavity flow problem. " try it and you will like it !"
  Reply With Quote

Old   March 8, 1999, 21:36
Default Re: A question about Peclet number
  #3
Zhong Lei
Guest
 
Posts: n/a
Let's discuss the most simple case, 1-D convection-diffusion equation of temperature,

u ( dT/dx ) - a ( d^2 T/dx^2 ) = 0,

where u has boundary-layer-like property, i.e., u(0) = 0 and u(1) = 1. By using 2nd-order central differential scheme, we get the Peclet number = u dx/a. To satisfy the condition Pe < 2, dx should be decreased as u is increased. How to explain?
  Reply With Quote

Old   March 9, 1999, 05:41
Default Re: A question about Peclet number
  #4
Essemiani
Guest
 
Posts: n/a
In fact for a more general problem of convection-diffusion , your pde has the following form:

d/dx(ro*u*T)=d/dx(nu*dT/dx), u and T are your flowfield variables.

then you can integrate this pde on a control volume .

Thev problem you'll have is the correct estimation of the variables values at the interfaces of your control volume.

to do this you can use different discretisation schemes:

1- central difference scheme: limitation Pe<2 2-exponential scheme 3-upwind scheme: not appropriate for low Pe, for high values of Pe ,the diffusion is over estimated. 4-hybrid scheme: maximum error for Pe=2

and finally the power law scheme, which is recommended ( no dependance to Pe)

for more details , check the book of Patankar " "numerical heat transfer and fluid flow"

good luck karim
  Reply With Quote

Old   March 9, 1999, 06:24
Default Re: A question about Peclet number
  #5
Zhong Lei
Guest
 
Posts: n/a
Patankar's blending scheme is very effective to maintain numerical stability, so it has been widely used in industrial applications. But it often results in losing accuracy of computation when Peclet number > 2. In this case, the accuracy is blended between first and second order.

Maybe the answer of my question is that, as the distance is increased away from the wall, the diffusion terms become smaller and smaller, thus the constraint of Peclet number becomes less important. Consquently a large spatial size is allowed.

But, this answer is incomplete in mathematics.
  Reply With Quote

Old   March 9, 1999, 10:13
Default Re: A question about Peclet number
  #6
John C. Chien
Guest
 
Posts: n/a
(1). specify the complete u-profile from x=0 to x=1,(2). specify T at x=0, and T at x=1,(3). write the finite-difference equation or the finite-volume equation from your 1-D convection-diffusion equation,using central difference for convection term and diffusion term,(4).write a Fortran code,(5). run the program, solve the 1-D problem,(6). plot the temperature distribution, examine the temperature profile, (7). change the mesh size based on the condition derived when necessary. (8). tell us your result about the temperature profile.
  Reply With Quote

Old   March 10, 1999, 01:37
Default Re: A question about Peclet number
  #7
Duane Baker
Guest
 
Posts: n/a
Hi, This is a very good suggestion from John C. and along the way, it may be useful to look at Example section 3.11 in Ferziger and Peric' "Computational Methods for Fluid Dynamics" which solves the 1-D Advection-Diffusion equation for a constant advection velocity(U) with UDS and CDS constant and variable grid spacing and compare to the exact solution......................................Duan e
  Reply With Quote

Old   March 11, 1999, 11:55
Default Re: A question about Peclet number
  #8
X. Ye
Guest
 
Posts: n/a
Dear Zhong Lei,

As I know, Fromm and Roache derived incorrectly this stability criteria (Pe<2) from the first Von Neumann stability condition (see Charls Hirsch "Numerical Computaions of Internal and External flows". The correct criteria of Hirt is: Pe<2/CFL, where CFL (Couran-Friedrich-Lewy) number is a product of eigen value and time step. For example, for the explicit method, CFL should be smaller than 1, that means, if CFL is sufficient small, Pe can be very large. But the Peclet number can give the impression about the accuracy.

X. Ye
  Reply With Quote

Old   March 12, 1999, 09:26
Default Re: continue my message
  #9
X. Ye
Guest
 
Posts: n/a
The CFL number is proportional to time step, the sum of sound and gas velocity but oppositional to the grid size. So, if you reduce the grid size, the stability will be reduced.

X. Ye
  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
DecomposePar unequal number of shared faces maka OpenFOAM Pre-Processing 6 August 12, 2010 09:01
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
Unaligned accesses on IA64 andre OpenFOAM 5 June 23, 2008 10:37
Question on Schmidt number ap Main CFD Forum 0 March 26, 2008 19:58
A question about the low Reynolds number k-¦Åmodel Haribort CFX 0 December 27, 2001 02:49


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