CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

LES channel flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 5, 2010, 11:52
Default LES channel flow
  #1
New Member
 
Ilia Popov
Join Date: Feb 2010
Posts: 2
Rep Power: 0
ilia is on a distinguished road
Helllo all!

I am trying to simulate channel turbulent flow with LES.

I am using a box of size 6x2x4 ( x is flowwise, z is spanwise, y is normal to walls).

Boundary conditions are cyclic in x and z directions. Noslip conditions at walls.

I set mu = 1/180, pressure gradient gradP = 1, and delta t = 1e-3. (This gives Re_tau = 180).

The mesh size is 32*32*32

I am using channelFoam solver modified to keep gradP constant.

As initial conditions I use boxTurb utility.

LES model i use is Smagorinsky with vanDriest damping.

The question is: What discretisation schemes shall I use?

I tried to use MUSCL scheme and CrankNicholson time scheme but there was some standing waves in flowwise with wavelength of 4 cells.

Then I tried to use some simplier schemes, like Euler and linear, the solution becomes laminar.

Could you suggest what schemes shall i use to both avoid standing waves and resolve turbulence?
ilia is offline   Reply With Quote

Old   February 9, 2010, 05:10
Default
  #2
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
For a channel flow, I suggest backward for time and linear for U.
eugene is offline   Reply With Quote

Old   February 9, 2010, 06:00
Default
  #3
New Member
 
Ilia Popov
Join Date: Feb 2010
Posts: 2
Rep Power: 0
ilia is on a distinguished road
Thank you!

So, you recommend using the same schemes as in the channelFoam tutorial. But why MUSCL gives not good results in this case and you do not recommend using it?
ilia is offline   Reply With Quote

Old   February 9, 2010, 12:42
Default
  #4
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
I have never tried MUSCL for LES, but it includes upwind components, so it will dissipate a lot of energy from the system which makes it unsuitable for conventional LES.
eugene is offline   Reply With Quote

Old   June 16, 2010, 16:46
Default
  #5
New Member
 
Matt James
Join Date: Jun 2010
Location: Marinette,WI, USA
Posts: 25
Rep Power: 7
mdjames is on a distinguished road
Hi Ilia,

How did you go about modifying channelFoam with regards to maintaining a constant pressure gradient?
mdjames is offline   Reply With Quote

Old   June 17, 2010, 06:12
Default
  #6
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
The current channelFoam modifies the imposed gradP in a goal-seek fashion to try and maintain a fixed flow rate. To apply a fixed gradP, you just have to change the code to use a constant gradP. Very simple.
eugene is offline   Reply With Quote

Old   April 8, 2011, 16:39
Default
  #7
New Member
 
David Thériault
Join Date: Nov 2010
Posts: 8
Rep Power: 6
dav1dt is on a distinguished road
Hello everyone

I think I'm at the right place.

I try to maintain a constant pressure gradient in a solver derived from channelFoam. I encounter some difficulties and I want to know if it comes from my modification, maybe someone could help me.

The pertinent part of this code is:

while (runTime.run())
{
#include "readTimeControls.H"
#include "readPISOControls.H"
#include "CourantNo.H"
#include "setDeltaT.H"

runTime++;

Info<< "Time = " << runTime.timeName() << nl << endl;
Info<< "Coriolis parameter = " << fc.value() << nl << endl;
Info<< "pressure gradient = " << gradP.value() << nl << endl;

sgsModel->correct();

fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
+ sgsModel->divDevBeff(U)
+ (fc^U)
==
-gradP
);

if (momentumPredictor)
{
solve(UEqn == -fvc::grad(p));
}

...................

Where graP is declared as a constant vector. When I use this solver, several timeStep after reaching a turbulent flow, the third component of fvc::grad(p) accelerates as a formula 1 even if I had declared a zero in the third component of my constant gradP.

My question is: knowing that I'm pretty sure that problem doesn't come from the other source term (fc^U), my modification appears to be legit or not? A hint could be very useful.

Thanks
dav1dt is offline   Reply With Quote

Old   April 18, 2011, 16:58
Default
  #8
New Member
 
David Thériault
Join Date: Nov 2010
Posts: 8
Rep Power: 6
dav1dt is on a distinguished road
I found it, beginner's problem!
dav1dt 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
LES of turbulent channel flows cedric_duprat OpenFOAM Running, Solving & CFD 118 July 6, 2015 19:51
setup problems - LES pipe flow with cyclic BC (1) and direct mapped inlet (2) florian_krause OpenFOAM 22 June 13, 2013 21:25
Please help with LES in a channel flow felipe Main CFD Forum 0 January 14, 2010 12:43
Fluent LES Channel Flow santosh FLUENT 0 January 7, 2010 22:11
Signal decay in les channel flow pankaj saha Main CFD Forum 0 July 25, 2008 12:34


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