|
[Sponsors] |
March 12, 2004, 16:17 |
Quick Question - Wall Function
|
#1 |
Guest
Posts: n/a
|
Dear CFD people:
Could you correct me if I am wrong: Is it true that, the way wall function is implemented is by: 1. Prescribed stress at the wall: If log wall is used then we have to do some iteration since it is impossible to express "wall stress" explicitly (Put it in right hand side of equation completely). 2. To calculate wall stress we need to know; Tangential velocity (U)at the "next-to-wall" cells; The "Y" => distance between the cell and the wall; Density at the "next-to-wall" cells; and Fluid viscosity. The question is, we need (U) to calculate wall strees, but at the same time we need the wall stress to calculate (U), because wall stress should decelerate (U) ??? How does it actually work ?/? Where can I get some references about "Prescribed Wall Stress" boundary condition. Thank's for the help D.T |
|
March 13, 2004, 15:08 |
Re: Quick Question - Wall Function
|
#2 |
Guest
Posts: n/a
|
See Lars Davidson: Numerical Methods for Turbulent Flow
http://www.tfd.chalmers.se/gr-kurs/MTF071 page 93 Best Regards |
|
March 16, 2004, 05:29 |
Re: Quick Question - Wall Function
|
#3 |
Guest
Posts: n/a
|
Dear D.T.
You are correct that the wall stress need to be computed using an iteratively method, as the log-law is: U/u_tau=LN(y+)/kappa (y+=u_tau*y*rho/mu) and hence both the RHS and the LHS need u_tau and thus the wall-stress from: tau_w=rho*u_tau^2 You can solve this by using an iterative scheme (guessing, computing and then updating). This is not that bad, as the whole CFD-process is an iterative process (steady-state), for time-accurate computations one should possibly be a bit more careful. The wall stress is then imposed in the code by either a source term: S=tau_w*A or via a modified turbulent viscosity from: mu_eff*U/y=tau_w |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
channelFoam for a 3D pipe | AlmostSurelyRob | OpenFOAM | 3 | June 24, 2011 14:06 |
Face merging, patch definition, wall function | Sandu69 | OpenFOAM | 2 | March 24, 2010 12:48 |
Question about bcdefw.f for wall temperature bc. | Jimmy | Siemens | 10 | March 18, 2008 16:28 |
Wall function problem in Fluent | mefpz | FLUENT | 1 | October 10, 2007 14:43 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |