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

LES of turbulent channel flows

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

Like Tree6Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   December 2, 2008, 11:35
Default Hi Philippe, Yes, you can d
  #21
Senior Member
 
Cedric DUPRAT
Join Date: Mar 2009
Location: Belgium
Posts: 169
Rep Power: 7
cedric_duprat is on a distinguished road
Hi Philippe,

Yes, you can do that, just telling that the shear stress in your periodic channel is compensated by the favorable pressure gradient.
So you get u_tau to get U+ and y+.
Do you plan to let your results avaiable ? I'm quite interresting in that just to compare with my results.

Notice that with the value of u_tau, you can calculate a Re_tau which is better to compare with DNS results.

Cedric
cedric_duprat is offline   Reply With Quote

Old   December 2, 2008, 17:48
Default Good. I'm glad it's the good w
  #22
Member
 
Philippe B. Vincent
Join Date: Mar 2009
Location: Quebec, Canada
Posts: 32
Rep Power: 7
philippebv is on a distinguished road
Good. I'm glad it's the good way to calculate u_tau. I don't know to what extend I will be working on this test case, but sure I can share the results. I'll keep you posted.

thanks,

Philippe
philippebv is offline   Reply With Quote

Old   December 9, 2008, 08:09
Default Hi Cedric, if you have used
  #23
Senior Member
 
Fabian Braennstroem
Join Date: Mar 2009
Posts: 407
Rep Power: 9
braennstroem is on a distinguished road
Hi Cedric,

if you have used any scripts for plotting the above profiles, it would be nice, to share these for the turbulence SIG :-) If there is some kind of 'standard', others might provide there results to some kind of database as well!?

Fabian
braennstroem is offline   Reply With Quote

Old   December 9, 2008, 09:40
Default Hi Fabian, Well, I used the
  #24
Senior Member
 
Cedric DUPRAT
Join Date: Mar 2009
Location: Belgium
Posts: 169
Rep Power: 7
cedric_duprat is on a distinguished road
Hi Fabian,

Well, I used the postChannel tool to get lot's of two colums files. but after that, it's quite .... home made :o)
I'll try to clean one of these gnuplot files to do a standard format.
Do you plan to do a LES channel tutorial ? :-)

Cedric
cedric_duprat is offline   Reply With Quote

Old   December 10, 2008, 04:22
Default Hi Cedric, yes, such a tuto
  #25
Senior Member
 
Fabian Braennstroem
Join Date: Mar 2009
Posts: 407
Rep Power: 9
braennstroem is on a distinguished road
Hi Cedric,

yes, such a tutorial would be a first step!? The base is probably already in the channelOodle tutorial, though comparisons with other filters and models would be nice. Even comparison with RANS models and their y+ sensitivity would be interesting.

Fabian
braennstroem is offline   Reply With Quote

Old   January 4, 2009, 03:34
Default Hi OpenFOAM I want to simulat
  #26
New Member
 
Jianying Jiao
Join Date: Mar 2009
Posts: 6
Rep Power: 7
jiao is on a distinguished road
Hi OpenFOAM
I want to simulation channel flow with dynSmagorinsky(channelOodel channel395), but i can't get "good" results for a channel flow with OpenFOAM.
I want to validate OpenFOAM from Moser DNS (Retau = 395).
the size of my channel : 4*2*2,80*50*60
the schemes of my channel:

ddtSchemes
{
default backward;
}

gradSchemes
{
default Gauss linear;
grad Gauss linear;
grad Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,k) Gauss limitedLinear 1;
div(phi,B) Gauss limitedLinear 1;
div Gauss linear;
div(phi,nuTilda) Gauss limitedLinear 1;
div((nuEff*dev(grad .T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A ),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DBEff,B) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}

it is in openFoam 1.4.1.
my result:



than you for your view.
Jianying
jiao is offline   Reply With Quote

Old   January 4, 2009, 03:38
Default Hi OpenFOAM I want to simulat
  #27
New Member
 
Jianying Jiao
Join Date: Mar 2009
Posts: 6
Rep Power: 7
jiao is on a distinguished road
Hi OpenFOAM
I want to simulation channel flow with dynSmagorinsky(channelOodel channel395), but i can't get "good" results for a channel flow with OpenFOAM.
I want to validate OpenFOAM from Moser DNS (Retau = 395).
the size of my channel : 4*2*2,80*50*60
the schemes of my channel:

ddtSchemes
{
default backward;
}

gradSchemes
{
default Gauss linear;
grad Gauss linear;
grad Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,k) Gauss limitedLinear 1;
div(phi,B) Gauss limitedLinear 1;
div Gauss linear;
div(phi,nuTilda) Gauss limitedLinear 1;
div((nuEff*dev(grad .T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A ),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DBEff,B) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}

it is in openFoam 1.4.1.
my result:






than you for your view.
Jianying
jiao is offline   Reply With Quote

Old   January 16, 2009, 13:45
Default Hello In the channel solver
  #28
Member
 
mohd mojab
Join Date: Mar 2009
Posts: 31
Rep Power: 7
mou_mi is on a distinguished road
Hello

In the channel solver for having periodic pressure, it use;

// Extract the velocity in the flow direction
dimensionedScalar magUbarStar =
(flowDirection & U)().weightedAverage(mesh.V());

// Calculate the pressure gradient increment needed to
// adjust the average flow-rate to the correct value
dimensionedScalar gradPplus =
(magUbar - magUbarStar)/rUA.weightedAverage(mesh.V());

Would anyone explain for me what is the meaning of this terms and how they work that can remain the mass flux constant?

Thank you
mou
mou_mi is offline   Reply With Quote

Old   January 17, 2009, 13:26
Default Hi Mou, the weightedAverage
  #29
Senior Member
 
Cedric DUPRAT
Join Date: Mar 2009
Location: Belgium
Posts: 169
Rep Power: 7
cedric_duprat is on a distinguished road
Hi Mou,

the weightedAverage commande calculate the volume average of a volSclalarField.
So if you have a volScalarField "toto", toto.weightedAverage(mesh.V()); gives you the average value of "toto" in the all domain.

With that, the line (flowDirection & U)().weightedAverage(mesh.V()); give the average velocity in the flowDirection.

In the transportProperties you fixed Ubar (and magUbar) so the term (magUbar - magUbarStar) gives the difference between the mass flow average you have in your calculation and the mass flow average you want.

You can write this error as a pressure drop in you periodic channel which explain the "gradPplus" name.

Then, you add a streamwise body force to the Navier-Stokes equation to correct this error.
The idea is to get (magUbar - magUbarStar) as small as possible, then the mass flow average you wanted is the mass flow averae you get in your calculation.

I hope it helps you.

Cedric
cedric_duprat is offline   Reply With Quote

Old   January 19, 2009, 10:35
Default Thank you Cedric, It was so h
  #30
Member
 
mohd mojab
Join Date: Mar 2009
Posts: 31
Rep Power: 7
mou_mi is on a distinguished road
Thank you Cedric,
It was so helpful.

Mou
mou_mi is offline   Reply With Quote

Old   February 13, 2009, 11:39
Default Hi Jianying: Have you found
  #31
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
Hi Jianying:

Have you found a way to improve your results? I did a similar simulation using channel395, and here is my results:




nzy102 is offline   Reply With Quote

Old   February 13, 2009, 11:50
Default Hi Jianying: Have you found
  #32
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
Hi Jianying:

Have you found a way to improve your results? I did a similar simulation using channel395, and here is my results:





Neither Umean or rms looks good to me.

Here is my blockmeshDict

================================================== =========================================
convertToMeters 1;

vertices
(
// (0 0 0)
// (4 0 0)
// (0 1 0)
// (4 1 0)
// (0 2 0)
// (4 2 0)
// (0 0 2)
// (4 0 2)
// (0 1 2)
// (4 1 2)
// (0 2 2)
// (4 2 2)

(0 0 0)
(6.28 0 0)
(0 1 0)
(6.28 1 0)
(0 2 0)
(6.28 2 0)
(0 0 3.14)
(6.28 0 3.14)
(0 1 3.14)
(6.28 1 3.14)
(0 2 3.14)
(6.28 2 3.14)

);

blocks
(
// hex (0 1 3 2 6 7 9 8) (40 25 30) simpleGrading (1 10.7028 1)
// hex (2 3 5 4 8 9 11 10) (40 25 30) simpleGrading (1 0.0984 1)
hex (0 1 3 2 6 7 9 8) (80 50 60) simpleGrading (1 10 1)
hex (2 3 5 4 8 9 11 10) (80 50 60) simpleGrading (1 0.1 1)
);

edges
(
);

patches
(
wall bottomWall
(
(0 1 7 6)
)
wall topWall
(
(4 10 11 5)
)
cyclic sides1
(
(0 2 3 1)
(6 7 9 8)
)
cyclic sides2
(
(2 4 5 3)
(8 9 11 10)
)
cyclic inout1
(
(1 3 9 7)
(0 6 8 2)
)
cyclic inout2
(
(3 5 11 9)
(2 8 10 4)
)
);

mergePatchPairs
(
);
================================================== =========================

SGS model I used:

LESmodel oneEqEddy;

delta vanDriest;
================================================== =========================

The initial condition of flow was obtained using perturbU. The flow was run for about 200 flow through time to get a steady state then sampled for another 200 flow through times. The problem I am having is:

I calculated utau based on mean pressure gradient. The utau based on 200 flow through times is 0.0073. The actual value is 0.0079. Is there anyway to improve this?

Thank you.

Ning
nzy102 is offline   Reply With Quote

Old   February 13, 2009, 11:59
Default Sorry, my bad. here is rms:
  #33
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
Sorry, my bad. here is rms:

/image{umean}
nzy102 is offline   Reply With Quote

Old   February 13, 2009, 11:59
Default Sorry, my bad. here is rms:
  #34
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
Sorry, my bad. here is rms:

/image{rms}
nzy102 is offline   Reply With Quote

Old   February 13, 2009, 12:03
Default Geeze, can't post an image:
  #35
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
Geeze, can't post an image:


nzy102 is offline   Reply With Quote

Old   February 18, 2009, 17:16
Default I missed contribution of SGS m
  #36
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
I missed contribution of SGS model from my previous calculation. Adding that (T = Bmean + R) gives the following plot:



The top two plots are taken from Eugene's thesis. Looks like that Van Driest model significantly overpredicts Vrms near the wall. So my question is how to improve my results? What is the best model in the openfoam I should use? Thank you.

Ning
nzy102 is offline   Reply With Quote

Old   February 18, 2009, 17:44
Default I missed contribution of SGS m
  #37
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
I missed contribution of SGS model from my previous calculation. Adding that (T = Bmean + R) gives the following plot:

Eugene's data:





my data:







The top two plots are taken from Eugene's thesis. Looks like that Van Driest model significantly overpredicts Vrms near the wall. So my question is how to improve my results? What is the best model in the openfoam I should use? Thank you.

Ning
nzy102 is offline   Reply With Quote

Old   February 18, 2009, 17:49
Default I missed contribution of SGS m
  #38
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
I missed contribution of SGS model from my previous calculation. Adding that (T = Bmean + R) gives the following plot:

Eugene's data:





my data:







The top two plots are taken from Eugene's thesis. Looks like that Van Driest model significantly overpredicts Vrms near the wall. So my question is how to improve my results? What is the best model in the openfoam I should use? Thank you.

Ning
nzy102 is offline   Reply With Quote

Old   February 19, 2009, 06:24
Default Hi Ning Thank you for your an
  #39
New Member
 
Jianying Jiao
Join Date: Mar 2009
Posts: 6
Rep Power: 7
jiao is on a distinguished road
Hi Ning
Thank you for your answer to my question. Your result is better than my result, but my mesh is poorer than yours.
I read some paper about channel flow, their result is better, their mesh (64*64*64), temporal term(RK3 or RK4). I think the reason which lead to my poor result is caused the temporal term. the temporal discretization term is second-order.
Give me some idear. Thank you.
Jianying
jiao is offline   Reply With Quote

Old   February 19, 2009, 12:11
Default Hi Jianying: Have you tried
  #40
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 7
nzy102 is on a distinguished road
Hi Jianying:

Have you tried one equation model? It should give you better results than smagrinsky model. Read Eugene's thesis. You will have better sense of what grid resolution you need. Regarding the schemes, I am using second-order schemes for both time and space. Hope this will help.

Ning
nzy102 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
Pressure inlet boundary conditions for open channel flows jack2000 OpenFOAM Running, Solving & CFD 3 October 21, 2012 14:10
LES In Turbulent in channel flow pankaj saha Main CFD Forum 8 April 15, 2009 11:34
LES In Turbulent in channel flow pankaj saha Main CFD Forum 13 August 9, 2008 02:51
Turbulent channel flow roberthino OpenFOAM Running, Solving & CFD 5 August 15, 2007 08:35
Bc for turbulent channel flow roberthino OpenFOAM Running, Solving & CFD 0 August 13, 2007 08:12


All times are GMT -4. The time now is 21:42.