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

LES In Turbulent in channel flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 17, 2008, 20:09
Default LES In Turbulent in channel flow
  #1
pankaj saha
Guest
 
Posts: n/a
Hi,

I am solving Turbulent channel flow by LES. But, i am facing problem as follows.---

I am not using any software. i am using my own code in fortran language.

For the channel----

i am taking periodic b.c in streamwise and spanwise direction.

dimension--4pi*2*2pi and grid 50*56*66

Re=205, corresponding friction velocity and half channel height.

So, i am using mean driving pressure gradient=1.0

for initial condition: u(x,y,z)=(1-y**8)+0.005*(rand(iseed)-0.5). /////Rand gives value , between 0,1, y is wall normal coordinate.

also, v(x,y,z)=0.0

w(x,y,z)=0.0

using---2nd order cd for convective and diffusive discretization. and adam-bashforth for time .

time step--0.0002

The problem is , i saw the mean velocity increase and reaches a very high value during the simulation.

also, along the channel centerline, i check the instantaneous u, v, w velocity----it does not show any turbulence.

what is the remedy....can anybody share his experience , during channel flow simulation.

thanks in advance.

  Reply With Quote

Old   July 18, 2008, 11:40
Default Re: LES In Turbulent in channel flow
  #2
underGroundMan
Guest
 
Posts: n/a
How long have you run it for? Are you sure its mean velocity that is increasing? How is instantaneous velocity behaving, increasing, decreasing or staying constant?

Regards
  Reply With Quote

Old   July 18, 2008, 18:58
Default Re: LES In Turbulent in channel flow
  #3
agg
Guest
 
Posts: n/a
1) I am assuming you are solving invompressible NS equations. Are you making sure your initial condition satisfies continuity? Check Kim, Moin and Moser's 1987 JFM paper for good initial condition.

2) Is your grid resolution enough to resolve the wall shear stress region? I had solved LES (Re=180) and had finer mesh than yours.

3) Make sure your LES model is working fine, unless you have already tested this in some other simulation.

Maybe you could start with a lower Re (180) as in the JFM paper and try to figure out what is going on (just a suggestion)

-agg
  Reply With Quote

Old   July 19, 2008, 03:03
Default Re: LES In Turbulent in channel flow
  #4
Tom
Guest
 
Posts: n/a
The resolution shouldn't be a problem for this Reynolds number. Probably the flow become laminar which implies that you have to change the initial conditions. What are using? Just some random small-scale disturbances? Take some initial disturbances that have a larger length scale and make sure that the disturbances are large enough. In my DNS I usually just take some random sine disturbances. Works fine.
  Reply With Quote

Old   July 19, 2008, 08:17
Default Re: LES In Turbulent in channel flow
  #5
pankaj saha
Guest
 
Posts: n/a
Hi,

First of all , thanks to agg and tom. My, Re is 205. And , my resolution near the wall is 0.68 wall unit. Still, according to agg, i will try, with Re=180.

How, could i understand , my intial condition is divergence free or not?

As, i am using a profile like below for i.c----

u(x,y,z)=(1-y^8)+10*(rand(iseed)-0.5) v(x,y,z)=0.0 w(x,y,z)=0.0

rand(iseed) will give value between 0,1.

So, my disturbence is -5,+5 at random values.

Could, you please, tell me , whethere, this field is satisfying continuity or not? ------------------------------------------------------

Tom, is my random disturbence fullfill the necessary large scale and amplitude.

tom, could you please, tell me what is the exact expression for i.c condition you , used.

Do, you give v=0.0, w=0.0 or anything else.

--------------------------------------------------------

Thanks pankaj

  Reply With Quote

Old   July 19, 2008, 09:05
Default Re: LES In Turbulent in channel flow
  #6
Tom
Guest
 
Posts: n/a
The disturbances that you impose have a very short wave length and viscosity or the subgrid model will almost immediately dissipate/damp most of these these disturbances. Therefore, I suggested to use disturbances with a longer wave length because these will no be damped so easily. Take just a couple sine/cosine disturbances, the fine details don't matter so much, just try something. It doesn't have to be divergence free (at least with my code).
  Reply With Quote

Old   July 19, 2008, 12:03
Default Re: LES In Turbulent in channel flow
  #7
pankaj saha
Guest
 
Posts: n/a
Thanks. Tom.---

Could you please, send me the expression you used.

one more thing---how you specify the mean pressure gradient?

I have used dp/dx=1.0, beccause, my code is nondimesionalized by half channel height and friction velocity. am i correct?

or what else you are using?
  Reply With Quote

Old   August 8, 2008, 18:27
Default Re: LES In Turbulent in channel flow
  #8
pankaj saha
Guest
 
Posts: n/a
Hi,

First of all , thanks to agg and tom. My, Re is 205. And , my resolution near the wall is 0.68 wall unit. Still, according to agg, i will try, with Re=180.

How, could i understand , my intial condition is divergence free or not?

As, i am using a profile like below for i.c----

u(x,y,z)=(1-y^8)+10*(rand(iseed)-0.5) v(x,y,z)=0.0 w(x,y,z)=0.0

rand(iseed) will give value between 0,1.

So, my disturbence is -5,+5 at random values.

Could, you please, tell me , whethere, this field is satisfying continuity or not? ------------------------------------------------------

Tom, is my random disturbence fullfill the necessary large scale and amplitude.

tom, could you please, tell me what is the exact expression for i.c condition you , used.

Do, you give v=0.0, w=0.0 or anything else.

-------------------------------------------------

Could you please, send me the expression you used.

one more thing---how you specify the mean pressure gradient?

I have used dp/dx=1.0, beccause, my code is nondimesionalized by half channel height and friction velocity. am i correct?

or what else you are using?

  Reply With Quote

Old   April 15, 2009, 11:34
Default
  #9
New Member
 
rowsan alam
Join Date: Apr 2009
Posts: 2
Rep Power: 0
rowsanalam is on a distinguished road
hi pankaj myself rowsan --plz help me with your code i am also facing with same problem as you--plz send me the code--in my id --rowsanalam@gmail.com


Quote:
Originally Posted by pankaj saha
;55894
Hi,

First of all , thanks to agg and tom. My, Re is 205. And , my resolution near the wall is 0.68 wall unit. Still, according to agg, i will try, with Re=180.

How, could i understand , my intial condition is divergence free or not?

As, i am using a profile like below for i.c----

u(x,y,z)=(1-y^8)+10*(rand(iseed)-0.5) v(x,y,z)=0.0 w(x,y,z)=0.0

rand(iseed) will give value between 0,1.

So, my disturbence is -5,+5 at random values.

Could, you please, tell me , whethere, this field is satisfying continuity or not? ------------------------------------------------------

Tom, is my random disturbence fullfill the necessary large scale and amplitude.

tom, could you please, tell me what is the exact expression for i.c condition you , used.

Do, you give v=0.0, w=0.0 or anything else.

--------------------------------------------------------

Thanks pankaj
rowsanalam is offline   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
Open Channel Flow ElanMorin FLUENT 4 February 25, 2015 17:26
turbulent developed lenght in oppen channel flow zandi Main CFD Forum 0 January 9, 2010 15:52
Initial conditionfor turbulent channel flow in LES pankaj saha Main CFD Forum 0 November 30, 2007 13:04
Turbulent channel flow roberthino OpenFOAM Running, Solving & CFD 5 August 15, 2007 08:35
Pressure drop across a turbulent channel flow shiv Main CFD Forum 0 September 21, 2005 00:52


All times are GMT -4. The time now is 10:08.