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

Initial conditionfor turbulent channel flow in LES

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 30, 2007, 13:04
Default Initial conditionfor turbulent channel flow in LES
  #1
pankaj saha
Guest
 
Posts: n/a
hi, i am solving fully developed plane channel flow by LES. I am using periodic boundary condition in steramwise(x-direction),and spanwise(z-direction).A no slip b.c for top and bottom wall(y-direction).

For initial condition i am taking a fully developed mean profile with random perturbation for the u-component (streamwise)of velocity. and v=0,w=0. like as follows...

For the whole domain---- u(x,y,z)=(1-y**8)+0.005*(RAND(x,y,z)-0.5 v(x,y,z)=0.0 w(x,y,z)=0.0

And in fortran programming language---

IND=0

do 55 i=1,iim

do 55 j=1,jim

do 55 k=1,kim

u(i,j,k)=(1-y(j)**8)+0.005*(RAND(IND)-0.5)

v(i,j,k)=0.0

w(i,j,k)=0.0

p(i,j,k)=1.0

IND=IND+1

55 continue

i,j,k--is spatial grid in three directions

iim=maximum grid no in x-direction

jim=maximum grid no in y-direction

kim=maximum grid no in z-direction

RAND(IND)------is a inbuilt function in fortran. As, the value of IND , change in the loop--then RAND(IND), will give a value , between--0.0-1.0

my question is---------

1. Is my approach correct.

2.Do i need random values for v and w velocity component.

3.Do i need to change the random function.

Please help me. A warm thanks in advance..

  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
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM 6 April 12, 2011 11:24
MRFSimpleFOAM goes divergenced! renyun0511 OpenFOAM Running, Solving & CFD 0 November 19, 2009 03:11
Negative value of k causing simulation to stop velan OpenFOAM Running, Solving & CFD 1 October 17, 2008 05:36


All times are GMT -4. The time now is 14:25.