
[Sponsors] 
November 17, 2000, 12:28 
Wall BC's for a NS solver

#1 
Guest
Posts: n/a

Sponsored Links
I am currently writing a 3D unsteady NavierStokes solver using the FVM and LaxWendroff scheme (1 step) for the time integration. I have a problem to implement the noslip BC's, at least to correct the viscous flux. The algorithm is the type of nodebased so that I have the state variables stored in the nodes of each cells. The shear stress is then evaluate in the midcell edges and contribute to each surrounding nodes. Concerning the wall boundary treatment, I put a dummy cell row inside the wall where each nodes contains a negative velocity corresponding to the real one inside the real flow domain. This dummy cell row contribute to the correction of fluxes for the next time step state values in the vertex nodes of the cells attach to the wall. When the code is run, the divergence appears quite quickly. Actually, It is seen that there is oscillation of pressure then velocity at the wall (like oddeven case ...) Has anyone of you an idea why I am facing this problem (I have already a numerical smoothing) and also how handle the wall BC's with the method I use. Cheers. 

Sponsored Links 
November 20, 2000, 12:12 
Re: Wall BC's for a NS solver

#2 
Guest
Posts: n/a

hi there,
I guess that you are using a collocated grid arrangement! If that's the case, to avoid some pressure velocity decorolation you must use the RhieChow model. Rhie, C.M. and Chow, W.L. "A numerical study of the turbulent flow past an isolated airfoil with trailing edge separation." AIAA Journal, Vol 21, 1983, pp 15251532. If you are actually using this model, well i'm sure that you could find some very useful, and clearly explained, information in this following paper: Morinishi, Lund, Vasilyev and Moin "Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow." J. Comp. Phys 143, 90124, 1998. You can download this paper online at the following address: http://mae.uta.edu/~lund/ I hope this helps. Sincerely, Frederic Felten CFD Laboratory UT Arlington. http://utacfdb.uta.edu/~felten/ 

November 20, 2000, 13:26 
Re: Wall BC's for a NS solver

#3 
Guest
Posts: n/a

You seem to have different problems with your solver, the nondivergence of the velocity field and the oddeven spatial oscillations.
The oddeven oscillation of pressure and velocity are due to their mutual decoupling. If even order discretizations are used to compute both the source term of the Poisson equation and the pressure gradient term(s) in the momentum equation(s), then the velocity and pressure at a given point are insensitive to each other (i.e., do not depend on each other). So you need to use some kind of biased (oddorder) discretization for one of these terms. This is a standard practice when using nonstaggered (collocated) grids. Rhie and Chow method can be used but a more recent method by Zang et al. (Journal of Comp. Phys., 1994) might be more helpful. One of the other authors is Prof. Street from Stanford Univ. Even with oscillations, you still should have a divergence free (but incorrect) velocity field when using nonstaggered grids. Divergence becomes nonzero eventually after the oscillation have grown significantly in amplitude. You were saying that the divergence starts appearing close to the wall which suggests that your pressure BC is not accurate. Spatial oscillations in the wall normal velocity component are produced if an appropriate BC is not used. See Gresho (Annual Review of Fluid Mech., year = 199?) for details. Here is a possible source of error at the walls. For each boundary cell, you have ghost cell. The velocity components in the ghost cell are negative of the components in the boundary cell, the cell face (i.e., the wall) has zero velocity and hence there is not convective flux. However, this works only for Euler flows. In incompressible viscous flows, velocity divergence condition demands that the wall normal component of velocity has zero gradient along the wall normal. If your wall is aligned to the xz plane then, dv/dy = 0 at the wall (i.e. the cell face where you compute the viscous flux). So there is no normal viscous stress at the wall. For this condition to be valid on the wall, v(ghost cell) = v(boundary cell). So you have to enforce this condition while computing viscous stresses. If you use the condition v(ghost cell) = v(boundary cell) that is required for the noslip condition, you create nonzero dv/dy and the near wall flow will not be accurate. 

November 20, 2000, 18:25 
Re: Wall BC's for a NS solver

#4 
Guest
Posts: n/a

Hi Kalyan,
A small correction. The method of Zang, Street, and Koseff is the same as the RhieChow method as far as the elimination of pv decoupling is concerned! chidu... 

November 23, 2000, 07:13 
help

#5 
Guest
Posts: n/a

je suis un chercheur marocain dans le domaine CFD hydraulic option les ecoulements a surface libre je cherche des articles ou livres qui parlent de la modelisation mathematique et numerique des écoulements surfaciques merci


November 23, 2000, 07:55 
Re: help

#6 
Guest
Posts: n/a

Pour obtenir de l'aide, je vous suggere de poster a nouveau votre message. Mais comme un nouveau message (presentement il est vu comme une reponse a une autre question). De plus, je crois qu'il serait preferable de le faire en englais. Vous pourrez ainsi rejoindre un plus grand nombre de chercheurs.
Bonne Chance. 

December 8, 2000, 22:45 
Re: Wall BC's for a NS solver

#7 
Guest
Posts: n/a

>For this condition to be valid on the wall, v(ghost >cell) = v(boundary cell). So you have to enforce this >condition while computing viscous stresses. If you
:use the condition v(ghost cell) = v(boundary cell) :that is required for the noslip condition, you :create nonzero dv/dy and the near wall flow will not :be accurate. I have never been able to stomach this image method on bounded domains...I can see that it is OK when the wall is infinite in the tangential directions, but I cannot see an image w.r.t. the wall of a (say) rectangular domain. Can anyone shed some light ? Best, Peter 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Very technical question about solving wall boundary layer ...  jlb001  FLUENT  6  December 27, 2014 06:56 
Define Two UDF Wall Temperature BC's at once?  Thomas Eire  Main CFD Forum  0  June 17, 2008 05:06 
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug  unoder  OpenFOAM Installation  11  January 30, 2008 21:30 
Multicomponent fluid  Andrea  CFX  2  October 11, 2004 05:12 
what the result is negatif pressure at inlet  chong chee nan  FLUENT  0  December 29, 2001 06:13 
Sponsored Links 