# WALL FUNCTION: Implement process

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

 May 21, 2004, 01:58 WALL FUNCTION: Implement process #1 boling2002 Guest   Posts: n/a I am currently working on implementing the wall function to Wilcox wall function.Now I am trying the case for flat plate. My code is using JST scheme. The main steps I am using are: 1) Compute u_{\tau} by using y^{+}=1/k ln( u^{+} + B with given velocity from the first interior cell. 2) Compute \tau_{w} using u_{\tau} and \rho. 3) Using the integral x-component of \tau_{w} multiply the maginitude of surface vector. 4) Set the engergy item to zero since I am working for steady adibatic case. Is there anything wrong in my implementation, I could not get the resonable results. Could anyone give me your hand ? Thanks

 June 8, 2004, 09:32 Re: WALL FUNCTION: Implement process #3 Jonas Bredberg Guest   Posts: n/a Hi there, I guess you simply forgot to tell us that you have included the following functions for k and omega: k=u_tau/sqrt(C_mu) omega=u_tau/(sqrt(C_mu)*\kappa*y) Furthermore it seems to be a typo for the log-law: u+=log(y+)/\kappa+C, \kappa=0.41, C=5 When computing u_tau (from the log-law and u+=U/u_tau) you could maybe iterative some times, as the log-law needs y+ which involves u_tau itself as: y+=u_tau*y*\rho/\mu What you mean by multiplying tau_w with the surface vector, and don't know. Though I guess you do it rightly. Tau_w is the surface tension which is introduced in FVM as a source term by multiplying with the appropriate surface-area.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11 yka8150 Main CFD Forum 0 September 21, 2009 23:08 matteo_gautero OpenFOAM Running, Solving & CFD 0 February 28, 2008 07:51 mefpz FLUENT 1 October 10, 2007 13:43 Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00

All times are GMT -4. The time now is 15:44.