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

LES/RANS on Channel Flow Velocity Profile

Register Blogs Community New Posts Updated Threads Search

Like Tree8Likes
  • 2 Post By FMDenaro
  • 1 Post By sbaffini
  • 1 Post By FMDenaro
  • 1 Post By sbaffini
  • 2 Post By flotus1
  • 1 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 25, 2017, 04:48
Default LES/RANS on Channel Flow Velocity Profile
  #1
New Member
 
Join Date: Aug 2013
Posts: 7
Rep Power: 13
hx293 is on a distinguished road
Hello everyone,

I just started learning CFD with the channel flow case. I tried different RANS models and LES SGS models on predicting the velocity profile at the channel outlet, and I found most RANS results agreed well with former experiments and DNS results, but LES results did not.

Is there any idea why the LES results deviate the exp. and DNS results so much?

Attached figs contain results and illustrations of the BC settings. The Reynolds number of the problem (defined by channel height and velocity at the channel center) is around 10,000, and the simulated channel length is enough to reach fully developed turbulent flow at channel outlet. The y+ value of the near wall mesh is around 0.2. The delta x+ (streamwise) value of the near wall mesh is around 6, and the delta y+ (pitchwise) value is around 4. The solver I used is Fluent.

Thank you for your attention.
Attached Images
File Type: jpg fig1.jpg (38.4 KB, 80 views)
File Type: png fig2.png (73.3 KB, 87 views)
File Type: png fig3.png (60.6 KB, 74 views)
hx293 is offline   Reply With Quote

Old   August 25, 2017, 05:17
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
First, one problem is that you are simulating a spatial developing BL and that requires much care in prescribing the LES inlet profile. Then, the velocity profile has a typical effect of a too diffusive impact of the SGS model (or also discretization).
1) please, give details about the inlet velocity and the discretization
2) try to simulate first the periodical channel flow at Re_tau=590 and compare with available DNS
sbaffini and juliom like this.
FMDenaro is offline   Reply With Quote

Old   August 25, 2017, 06:07
Default
  #3
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,190
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
You can't really do it that way. A developing BL is quite an order of magnitude more difficult to get right with respect to a channel flow. You can use the latter to find suitable discretization settings while taking the bc out of the picture. Then you can work on a suitable set of bc for your BL setting. Also, the unstructured grid part is certainly not helping. You might find my phd thesis helpful in setting a channel flow with Fluent:

https://www.cfd-online.com/Forums/bl...hesis-les.html
FMDenaro likes this.
sbaffini is offline   Reply With Quote

Old   August 27, 2017, 00:51
Default
  #4
New Member
 
Join Date: Aug 2013
Posts: 7
Rep Power: 13
hx293 is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
You can't really do it that way. A developing BL is quite an order of magnitude more difficult to get right with respect to a channel flow. You can use the latter to find suitable discretization settings while taking the bc out of the picture. Then you can work on a suitable set of bc for your BL setting. Also, the unstructured grid part is certainly not helping. You might find my phd thesis helpful in setting a channel flow with Fluent:

https://www.cfd-online.com/Forums/bl...hesis-les.html
Hi Paolo,

Thank you for your comment. I just have few more questions:

1. In this model, why steady RANS method matched the exp./DNS results while unsteady LES did not? Is it caused by steady/unsteady settings, or caused by RANS/LES methods?

2. Why unstructured grid is not recommended compared to structured one? Will increasing the grid density of unstructured mesh be helpful? Or we must use structured mesh in LES?

I am looking forward to your reply.
hx293 is offline   Reply With Quote

Old   August 27, 2017, 01:09
Default
  #5
New Member
 
Join Date: Aug 2013
Posts: 7
Rep Power: 13
hx293 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
First, one problem is that you are simulating a spatial developing BL and that requires much care in prescribing the LES inlet profile. Then, the velocity profile has a typical effect of a too diffusive impact of the SGS model (or also discretization).
1) please, give details about the inlet velocity and the discretization
2) try to simulate first the periodical channel flow at Re_tau=590 and compare with available DNS
Hi Filippo,

Thank you for your comments.

I am using uniform velocity at the inlet, with the turbulent intensity of 5% induced by the Spectral Synthesizer method. Attached gif illustrates the inlet vorticity fluctuation induced by this method. The near-wall mesh of the case is structured, its sizing has been mentioned before. The mesh at the channel center is unstructured.

I agree with you that the periodical channel flow case is a better idea, and I found it quite popular in some recent JFM papers. However, I am not curious about the transition process, but only interested in the fully developed velocity profile. Is this model accurate enough for this purpose? Maybe with further refinement of the mesh?

I am looking forward to your reply.
Attached Images
File Type: jpg Vorticity_pit.jpg (46.7 KB, 51 views)
hx293 is offline   Reply With Quote

Old   August 27, 2017, 06:47
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
The case you want to simulate is not so simple as it would seem...

Starting from a uniform profile, you need to solve the transition that, even if you supplied a superimposed fluctuation, is not physical but numerical.
I see also that many RANS solutions do not match enough the DNS solution.

I have some further question

1) how do you compute the u_tau velocity? (I see you have Re_tau=300, right?). What is the criterion you use for the streamwise lenght?
2) did you use the static or dynamic Smagorinsky model?
3) are you using the central bounded scheme?
4) what about the spanwise lenght of the domain where you set periodic conditions?
5) the first part of the domain has "symmetric" condition. Why?
FMDenaro is offline   Reply With Quote

Old   August 27, 2017, 12:37
Default
  #7
New Member
 
Join Date: Aug 2013
Posts: 7
Rep Power: 13
hx293 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The case you want to simulate is not so simple as it would seem...

Starting from a uniform profile, you need to solve the transition that, even if you supplied a superimposed fluctuation, is not physical but numerical.
I see also that many RANS solutions do not match enough the DNS solution.

I have some further question

1) how do you compute the u_tau velocity? (I see you have Re_tau=300, right?). What is the criterion you use for the streamwise lenght?
2) did you use the static or dynamic Smagorinsky model?
3) are you using the central bounded scheme?
4) what about the spanwise lenght of the domain where you set periodic conditions?
5) the first part of the domain has "symmetric" condition. Why?

Before heading to the details, I would like to highlight my motivation that I am only interested in the velocity profile of the fully-developed turbulent channel flow.

Although turbulent fluctuation exist, the time-averaged velocity profile of the fully-developed channel flow should stay the same through streamwise direction. My hypothesis (might be a stupid one) is that given enough channel length, the non-physical phenomena near the channel inlet (i.e. unappropriated velocity profile as inlet BC, etc.) has no influence on the downstream fully-developed velocity profile. My RANS results actually fits well with the exp./DNS, which shows that the hypothesis is not that bad.

Back to the details:

(1) u_tao = sqrt(tau/density), and tau is calculated by viscosity times velocity gradient; Re_tau ~ 300 is correct. The streamwise length is chose as the minimum length that required to achieve fully-developed flow at outlet, which is determined by trial and error.
(2) I am using static Smagorinsky model.
(3) I am using the central bounded scheme.
(4) The calculation domain size is 100mm*10mm*3mm in streamwise, spanwise(vertical to the channel plane) and pitchwise(vertical to the periodic plane) directions, respectively.
(5) The symmetric BC is setted because I think it is one of the physically-right scenario.

I only had basic CFD courses, but not studied turbulence modeling. Any recommended readings (for beginners) are welcomed.

Thank you for helping.
hx293 is offline   Reply With Quote

Old   August 27, 2017, 12:48
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
First, I do not reccomend the static Smagorinsgy model. Try setting the dynamic model.

Then, RANS and LES are very different in the meaning of the variables you solve, the former being statistically steady, the latter being time dependent.
That means that from the LES solutions, you need to apply ensemble averaging to get a meaningful steady profile. The lenght to have a RANS solution well developed is not as same as the lenght required for LES.

I strongly suggest also to check the spanwise energy spectra to see if the LES solution is quite correct.

There are many textbooks that should be read to go in deep in the simulation of turbulence.
I suggest to my students this list:
https://www.researchgate.net/publica..._Part_I_and_II
juliom likes this.
FMDenaro is offline   Reply With Quote

Old   August 27, 2017, 13:07
Default
  #9
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Just to add that I suggest also to do a run without any SGS model (LES no-model or unresolved DNS) to check the real effect of the model onto the solution.
FMDenaro is offline   Reply With Quote

Old   August 28, 2017, 05:44
Default
  #10
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,190
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
As stated by Filippo, RANS and LES/DNS are two really different matters.

The way you approached the problem in RANS is ok. It can be steady and you don't even need the third spanwise dimension. Also, the grid is of no great concern.

The difference in LES for this case is that you actually want to reproduce all the main mechanisms leading to the transition or, for fully turbulent scenarios, the BL growth.

Basically, this requires accuracy down to a certain flow scale. Fluent by itself, the grid you choose in particular, not to mention the specific scheme and SGS model you are using, are not well suited for this problem. Eventually, for very fine grids, they are going to be ok, but we are probably talking order 10^8 cells for your problem.

While there is probably room for improvement, I have no experience on this (LES of BL in Fluent).

You can try with these modifications:

- no sgs model at all (laminar or static smagorinsky with the constant equal to 0)
- unbounded central scheme (to be activated via TUI if using laminar)
- fully structured grid (should not be difficult in your case)
- second order time advancement with a dt leading to a courant around 0.1

However, I strongly suggest to go back to a simpler case to start LES while also reading some book on it. Otherwise your main risk is wasting a huge amount of time and resources without ever getting any meaningful result.
FMDenaro likes this.
sbaffini is offline   Reply With Quote

Old   August 28, 2017, 06:54
Default
  #11
Super Moderator
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,427
Rep Power: 49
flotus1 has a spectacular aura aboutflotus1 has a spectacular aura about
Quote:
Originally Posted by hx293 View Post
Before heading to the details, I would like to highlight my motivation that I am only interested in the velocity profile of the fully-developed turbulent channel flow.
Then I would recommend to model a fully developed channel flow instead of a developing channel flow in a long channel and assuming that it is long enough.
This is usually done with periodic boundary conditions between "inlet" and "outlet" and the symmetry boundaries respectively.
FMDenaro and RobertHB like this.
flotus1 is offline   Reply With Quote

Old   August 29, 2017, 01:09
Default
  #12
New Member
 
Join Date: Aug 2013
Posts: 7
Rep Power: 13
hx293 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Just to add that I suggest also to do a run without any SGS model (LES no-model or unresolved DNS) to check the real effect of the model onto the solution.
Now I kind of understand the limitation of my model. I tried your suggestion to run an unresolved DNS case, and the result is very similar to my former LES results. Thus, it must be the numerical model.

I will try to run a simple case to avoid the transition process. Thank you for your help.
hx293 is offline   Reply With Quote

Old   August 29, 2017, 01:17
Default
  #13
New Member
 
Join Date: Aug 2013
Posts: 7
Rep Power: 13
hx293 is on a distinguished road
Quote:
Originally Posted by flotus1 View Post
Then I would recommend to model a fully developed channel flow instead of a developing channel flow in a long channel and assuming that it is long enough.
This is usually done with periodic boundary conditions between "inlet" and "outlet" and the symmetry boundaries respectively.
Thank you for your comments. Your suggestion sounds clever!

But I am not sure if I could set the inlet flow velocity (in other words, the Reynolds number) as I want. In my RANS model, this can be easily achieved by the inlet condition directly. In your suggested model, how to set this velocity?
hx293 is offline   Reply With Quote

Old   August 29, 2017, 04:26
Default
  #14
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by hx293 View Post
Now I kind of understand the limitation of my model. I tried your suggestion to run an unresolved DNS case, and the result is very similar to my former LES results. Thus, it must be the numerical model.

I will try to run a simple case to avoid the transition process. Thank you for your help.
This is a typical problem due to the numerical viscosity introduced by the scheme
randolph likes this.
FMDenaro is offline   Reply With Quote

Old   August 29, 2017, 12:17
Default
  #15
Super Moderator
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,427
Rep Power: 49
flotus1 has a spectacular aura aboutflotus1 has a spectacular aura about
Quote:
Originally Posted by hx293 View Post
Thank you for your comments. Your suggestion sounds clever!

But I am not sure if I could set the inlet flow velocity (in other words, the Reynolds number) as I want. In my RANS model, this can be easily achieved by the inlet condition directly. In your suggested model, how to set this velocity?
Well thanks a lot, but this is not really my idea. It is just the commonly accepted method to simulate fully developed channel flows.
Implementations may vary depending on the software you use, but usually a body force is used to drive the flow. You might have to make a few adjustments to the magnitude of the force until you get the desired flow rate.
flotus1 is offline   Reply With Quote

Old   August 29, 2017, 12:34
Default
  #16
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I agree, setting periodic streamwise and spanwise conditions is a way to study a temporal evolving boundary layer, no velocity inlet/outlet must be prescribed but only cyclic conditions. The driving force is the constant pressure gradient or the mass-prescribed constraint.

The thesis of Paolo describes the setting, if you want to se a comparison with other LES codes, you can have a look here:
https://www.researchgate.net/publica...mulation_codes
FMDenaro is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to calculate mean velocity profile for channel flow michel1988 OpenFOAM Post-Processing 3 November 20, 2019 10:09
atmBoundaryLayerInletVelocity - Velocity Profile not continuous through domain sdfij6354 OpenFOAM Running, Solving & CFD 3 July 26, 2017 17:16
UDF parbolic inlet velocity profile for 3D channel flow. srv537 FLUENT 0 August 6, 2016 03:21
InterFoam - Validation for velocity profile in simple channel me.ouda OpenFOAM Running, Solving & CFD 0 October 19, 2015 07:42
How to make a UDF to have a velocity profile in a square channel Gigis Fluent UDF and Scheme Programming 8 January 13, 2013 23:20


All times are GMT -4. The time now is 13:58.