
[Sponsors] 
March 8, 1999, 07:38 
A question about Peclet number

#1 
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 convectivediffusive 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? 

March 8, 1999, 10:05 
Re: A question about Peclet number

#2 
Guest
Posts: n/a

Remember that at each point there are three velocity components, so, there are three such numbers. In a more general convectiondiffusion equation, they are related to the coefficient of the firstorder differential terms. The coefficient does not have to be the velocity components. But there are always three firstorder terms. The numbers associated with each direction ( coordinates) are not the same. For 1D, there is one cell Reynolds number, for 2D, there are two cell Reynolds number, and for 3D, there are three cell Reynolds number,etc... A good exercise is : try to solve a standard 2D cavity flow problem. " try it and you will like it !"


March 8, 1999, 21:36 
Re: A question about Peclet number

#3 
Guest
Posts: n/a

Let's discuss the most simple case, 1D convectiondiffusion equation of temperature,
u ( dT/dx )  a ( d^2 T/dx^2 ) = 0, where u has boundarylayerlike property, i.e., u(0) = 0 and u(1) = 1. By using 2ndorder 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? 

March 9, 1999, 05:41 
Re: A question about Peclet number

#4 
Guest
Posts: n/a

In fact for a more general problem of convectiondiffusion , 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 2exponential scheme 3upwind scheme: not appropriate for low Pe, for high values of Pe ,the diffusion is over estimated. 4hybrid 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 

March 9, 1999, 06:24 
Re: A question about Peclet number

#5 
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. 

March 9, 1999, 10:13 
Re: A question about Peclet number

#6 
Guest
Posts: n/a

(1). specify the complete uprofile from x=0 to x=1,(2). specify T at x=0, and T at x=1,(3). write the finitedifference equation or the finitevolume equation from your 1D convectiondiffusion equation,using central difference for convection term and diffusion term,(4).write a Fortran code,(5). run the program, solve the 1D 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.


March 10, 1999, 01:37 
Re: A question about Peclet number

#7 
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 1D AdvectionDiffusion equation for a constant advection velocity(U) with UDS and CDS constant and variable grid spacing and compare to the exact solution......................................Duan e


March 11, 1999, 11:55 
Re: A question about Peclet number

#8 
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 (CouranFriedrichLewy) 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 

March 12, 1999, 09:26 
Re: continue my message

#9 
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 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
DecomposePar unequal number of shared faces  maka  OpenFOAM PreProcessing  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 