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

Boundary conditions of nut in LES

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

Like Tree15Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 21, 2019, 04:37
Default Boundary conditions of nut in LES
  #1
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
Hello Foamers,

From OpenFOAM tutorials, for most LES cases, zeroGradient is used at inlets, outlets and walls for nut (the turbulent viscosity). Why does nut behave like zeroGradient at the walls? Theoretically, I would think it should equal to 0.

I ask this question because when I simulate air flow inside a circular pipe (bulk velocity 66 m/s, so subsonic flow. Re is around 15,000), by using zeroGradient for nut, I always end up with a velocity profile that has large velocity at the center, which looks like a laminar flow. By using fixedValue of 0 for nut, I get a flatter velocity profile, which makes more sense because it is turbulence.

I'm guessing that this should also have something to do with the mesh resolution. Maybe zeroGradient works fine if the mesh is very fine near the wall?

One thing I've noticed is that for RANS, calculated is used for nut, and a wall function is often specified. What if I use these wall functions in LES? Is that physically right?

Thanks in advance,
Ruiyan
cryabroad is offline   Reply With Quote

Old   October 21, 2019, 06:08
Default
  #2
New Member
 
Adam
Join Date: Jan 2019
Posts: 21
Rep Power: 7
boundary93 is on a distinguished road
zeroGradient with LES makes only sense if you y+ value is below 1, otherwise you should use a wallmodel.
boundary93 is offline   Reply With Quote

Old   October 21, 2019, 23:42
Default
  #3
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
Quote:
Originally Posted by boundary93 View Post
zeroGradient with LES makes only sense if you y+ value is below 1, otherwise you should use a wallmodel.
I agree with you, but what about fixedValue of 0? If my mesh is fine enough (y+<1 for the first cell center), it should be fixedValue or zeroGradient? In my view, physically, no turbulence is present at the wall, so nut should be 0 there.

What confuses me is that why nut needs boundary values? How are they used? I'm always under the impression that nut only has internal values, which should be enough for calculating velocities.

Still waiting for other people to shed some light on this topic!
cryabroad is offline   Reply With Quote

Old   October 22, 2019, 04:43
Default
  #4
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 16
Santiago is on a distinguished road
Quote:
In my view, physically, no turbulence is present at the wall, so nut should be 0 there.
It's not your view, it is like that. The thing is that OPENFOAM implements integral wall models in a rather esoteric way: by changing viscosity instead of velocities at the wall hence one of the reasons why you need boundary values for nut.

Quote:
What confuses me is that why nut needs boundary values?
This depends on what model you are using for turbulence. For instance, the SA model proposes a transport equation for nu (or a variant of nu), thus you need to specify boundary conditions for said equation.

On the other hand, the Smagorinsky family of models does not require boundary values for nut, since this is obtained by algebraic operations. In that case, your BCs must be 'calculated', unless you use wall models.
Santiago is offline   Reply With Quote

Old   October 22, 2019, 23:29
Default
  #5
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
Thank you for your input, it is very helpful!

For your first point about nut being 0 at the wall, I guess you are mainly talking about RANS? I come across this article (http://www.tfd.chalmers.se/~hani/kur...nfoamFinal.pdf) that explains the idea of changing viscosity instead of velocities in OpenFOAM, which is exactly what you mentioned. I think one needs to be very cautious when using these BCs in RANS. Most of them seem to work properly only for a certain mesh.

As for your second point, I am doing LES, so it points to me the right direction I think. Right now I'm testing results by using zeroGradient or fixedValue of 0. I'm also thinking about using nutUSpaldingWallFunction. From the article I mentioned above, it is said that this BC is kind of a "generic" boundary condition, in the sense that it automatically detects the y+ of the first cell center and operates in different modes based on your 1st cell y+. What I should test is actually calculated!
cryabroad is offline   Reply With Quote

Old   October 24, 2019, 05:47
Default
  #6
Senior Member
 
Ehsan Asgari
Join Date: Apr 2010
Posts: 473
Rep Power: 18
syavash is on a distinguished road
Quote:
Originally Posted by cryabroad View Post
Thank you for your input, it is very helpful!

For your first point about nut being 0 at the wall, I guess you are mainly talking about RANS? I come across this article (http://www.tfd.chalmers.se/~hani/kur...nfoamFinal.pdf) that explains the idea of changing viscosity instead of velocities in OpenFOAM, which is exactly what you mentioned. I think one needs to be very cautious when using these BCs in RANS. Most of them seem to work properly only for a certain mesh.

As for your second point, I am doing LES, so it points to me the right direction I think. Right now I'm testing results by using zeroGradient or fixedValue of 0. I'm also thinking about using nutUSpaldingWallFunction. From the article I mentioned above, it is said that this BC is kind of a "generic" boundary condition, in the sense that it automatically detects the y+ of the first cell center and operates in different modes based on your 1st cell y+. What I should test is actually calculated!
Hi,

As it was previously mentioned, calculated is a reasonable choice for algebraic eddy-viscosity models. As a suggestion, you can test the nut BCs on the channel395 tutorial. Mesh is small and the results are quick to obtain.
Nevertheless, I would expect to see identical results from zeroGradient and calculated.

Regards,
Syavash
syavash is offline   Reply With Quote

Old   October 24, 2019, 13:15
Default
  #7
gu1
Senior Member
 
Guilherme
Join Date: Apr 2017
Posts: 245
Rep Power: 10
gu1 is on a distinguished road
For y+<1 I use 'calculated'...
gu1 is offline   Reply With Quote

Old   October 25, 2019, 06:18
Default
  #8
New Member
 
Adam
Join Date: Jan 2019
Posts: 21
Rep Power: 7
boundary93 is on a distinguished road
When you have y+ lower than 1, its not important because the value of nut goes towards zero.

At least that's my experience.
boundary93 is offline   Reply With Quote

Old   October 28, 2019, 23:54
Default
  #9
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
I can confirm this behavior, at least for my little circular pipe calculations. I used a grid that has y+ around 0.5, and fixedValue of 0 and zeroGradient give me exactly the same results. See attached figure of u+ vs. y+. The velocity profile differs from the expected theoretical line though, and I'm trying to figure out the reason.

One quick question though, how do we get the friction velocity from OpenFOAM? My way of doing it is to use the post processing tool to calculate wallShearStress, and based on it the friction velocity can be calculated, the square root of the wallShearStress at the wall in this case (incompressible). Any better ways?
Attached Images
File Type: png uplus.png (43.9 KB, 235 views)
cryabroad is offline   Reply With Quote

Old   October 29, 2019, 00:03
Default
  #10
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
Thank you syavash, I end up of doing several tests on my simple circular pipe instead of the channel case. The results seem to show that, for not very fine mesh (y+ ~ 15), calculated does give almost the same u+ vs. y+ line compared to fixedValue of 0, but not zeroGradient. I'll post the results after I run some additional simulations to confirm this.

I'm getting a feeling that for a coarse mesh, maybe zeroGradient is not a good choice. Any other BCs (calculated, the various wall functions for nut in OpenFOAM) may work better. It makes sense though, otherwise there won't be people developing all these wall functions to make the velocity profile more realistic. Any comments are welcome!
cryabroad is offline   Reply With Quote

Old   October 29, 2019, 03:35
Default
  #11
Senior Member
 
Ehsan Asgari
Join Date: Apr 2010
Posts: 473
Rep Power: 18
syavash is on a distinguished road
Quote:
Originally Posted by cryabroad View Post
Thank you syavash, I end up of doing several tests on my simple circular pipe instead of the channel case. The results seem to show that, for not very fine mesh (y+ ~ 15), calculated does give almost the same u+ vs. y+ line compared to fixedValue of 0, but not zeroGradient. I'll post the results after I run some additional simulations to confirm this.

I'm getting a feeling that for a coarse mesh, maybe zeroGradient is not a good choice. Any other BCs (calculated, the various wall functions for nut in OpenFOAM) may work better. It makes sense though, otherwise there won't be people developing all these wall functions to make the velocity profile more realistic. Any comments are welcome!
Hi,

Thanks for sharing your results. I never tried such a coarse mesh with y+~15 at the wall. Regarding your results, are you certain that the flow has become turbulent? Can you provide some snapshots of instantaneous velocity contours?
The way you calculate u_tau is correct. However, I use the averaged velocity field for calculating tau_wall.

Regards,
Syavash
syavash is offline   Reply With Quote

Old   October 29, 2019, 22:57
Default
  #12
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
That's the thing I'm thinking about as well! Please see the attached instantaneous velocity contour, and I don't think that signals strong turbulence? Flow (Re around 15000) is from left to right, top and bottom are walls, with nut set to calculated. In the axial direction I used 125 cells for a length of 50mm, and about 30 cells in the radial direction. I applied the 1/7 power law velocity profile at the inlet, and used the LEMOS inflowGenerator from Rostock. It basically injects Lagrangian vortons (as they called) into the domain with prescribed length scale and Reynolds stress.

The reason I'm using this coarse mesh near the wall is that this short circular pipe will be served as a short inlet tube for a combustion chamber behind, and we are mainly interested in mixing of things inside the chamber. The short inlet pipe is needed because swirl flows will be introduced as well (in the long term) so it would be great if we have this pipe and let the flow develop before entering into the chamber.

Basically, the output of the pipe calculation is to provide a reasonable (at least not totally wrong) input for the chamber. So, we want to use as less cells as possible to achieve that. Maybe in the end this will be found impossible and I might as well just use as many cells as possible for this inlet pipe, but it would be ideal to use less cells there.
Attached Images
File Type: png ins_v.png (18.9 KB, 187 views)
cryabroad is offline   Reply With Quote

Old   October 30, 2019, 11:50
Default
  #13
Senior Member
 
Ehsan Asgari
Join Date: Apr 2010
Posts: 473
Rep Power: 18
syavash is on a distinguished road
Quote:
Originally Posted by cryabroad View Post
That's the thing I'm thinking about as well! Please see the attached instantaneous velocity contour, and I don't think that signals strong turbulence? Flow (Re around 15000) is from left to right, top and bottom are walls, with nut set to calculated. In the axial direction I used 125 cells for a length of 50mm, and about 30 cells in the radial direction. I applied the 1/7 power law velocity profile at the inlet, and used the LEMOS inflowGenerator from Rostock. It basically injects Lagrangian vortons (as they called) into the domain with prescribed length scale and Reynolds stress.

The reason I'm using this coarse mesh near the wall is that this short circular pipe will be served as a short inlet tube for a combustion chamber behind, and we are mainly interested in mixing of things inside the chamber. The short inlet pipe is needed because swirl flows will be introduced as well (in the long term) so it would be great if we have this pipe and let the flow develop before entering into the chamber.

Basically, the output of the pipe calculation is to provide a reasonable (at least not totally wrong) input for the chamber. So, we want to use as less cells as possible to achieve that. Maybe in the end this will be found impossible and I might as well just use as many cells as possible for this inlet pipe, but it would be ideal to use less cells there.

I guess it has not become fully turbulent yet. Why use LeMOS though? A better and more efficient way is using mapped boundary condition to develop turbulence for this geometry. See pitzDaily tutorial.

Regards,
Syavash
syavash is offline   Reply With Quote

Old   October 30, 2019, 22:16
Default
  #14
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
Thanks, I will definitely look into the mapping method!

I use LeMOS because it is similar to the vortex method in FLUENT. I have used FLUENT before and it works well, so I decided to use a similar inlet method.

Regards,
Ruiyan
cryabroad is offline   Reply With Quote

Old   November 4, 2019, 04:58
Default
  #15
New Member
 
Adam
Join Date: Jan 2019
Posts: 21
Rep Power: 7
boundary93 is on a distinguished road
For this kind of geometry I would recommend cyclic instead of mapped.
boundary93 is offline   Reply With Quote

Old   March 6, 2020, 09:12
Default
  #16
Member
 
bany
Join Date: Nov 2019
Posts: 50
Rep Power: 8
bany is on a distinguished road
Hi, Ruiyan. I think i have same problem that how to obtain a fully developed turbulence. And my domain is a cylinder too. Did you obtain a fully developed turbulence? mapped or cyclic? And what was your nut BC?
Thank you very much.
bany is offline   Reply With Quote

Old   April 12, 2020, 23:19
Default
  #17
Member
 
bany
Join Date: Nov 2019
Posts: 50
Rep Power: 8
bany is on a distinguished road
Quote:
Originally Posted by syavash View Post
I guess it has not become fully turbulent yet. Why use LeMOS though? A better and more efficient way is using mapped boundary condition to develop turbulence for this geometry. See pitzDaily tutorial.

Regards,
Syavash
Hi, when i use the mapped BC, i cannot get a fully developed pipe flow. Because i do not want to use too many cells to ensure y+<1. Can i use the wall functions or which wall functions can perform well in LES?
You can get more details in a fully developed pipe flow in LES with mapped BC
bany is offline   Reply With Quote

Old   April 24, 2020, 07:39
Default
  #18
Senior Member
 
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 10
cryabroad is on a distinguished road
I end up using the inflow generator I mentioned in previous posts, the one that comes with LEMOS. Here is the link: https://github.com/LEMOS-Rostock/LEM...nflowGenerator

For nut I used nutUSpaldingWallFunction on solid walls, and calculated on other boundaries. I have to say though, that I'm not doing very rigorous confined flow simulation. I'm more interested in the large chamber behind the feeding pipes, so as long as I have a more or less correct velocity profile at the end of the pipe (which connects to the chamber) together with some turbulence I'm fine. It may not be your case.
XJ_Wang and guanjiang.chen like this.
cryabroad is offline   Reply With Quote

Old   April 24, 2020, 09:00
Default
  #19
Member
 
bany
Join Date: Nov 2019
Posts: 50
Rep Power: 8
bany is on a distinguished road
Quote:
Originally Posted by cryabroad View Post
I end up using the inflow generator I mentioned in previous posts, the one that comes with LEMOS. Here is the link: https://github.com/LEMOS-Rostock/LEM...nflowGenerator

For nut I used nutUSpaldingWallFunction on solid walls, and calculated on other boundaries. I have to say though, that I'm not doing very rigorous confined flow simulation. I'm more interested in the large chamber behind the feeding pipes, so as long as I have a more or less correct velocity profile at the end of the pipe (which connects to the chamber) together with some turbulence I'm fine. It may not be your case.
OK, thanks for your replying! I will have a try.
bany is offline   Reply With Quote

Old   May 25, 2020, 10:12
Default
  #20
Member
 
Guanjiang Chen
Join Date: Apr 2020
Location: Bristol, United Kingdom
Posts: 54
Rep Power: 6
guanjiang.chen is on a distinguished road
Quote:
Originally Posted by cryabroad View Post
I end up using the inflow generator I mentioned in previous posts, the one that comes with LEMOS. Here is the link: https://github.com/LEMOS-Rostock/LEM...nflowGenerator

For nut I used nutUSpaldingWallFunction on solid walls, and calculated on other boundaries. I have to say though, that I'm not doing very rigorous confined flow simulation. I'm more interested in the large chamber behind the feeding pipes, so as long as I have a more or less correct velocity profile at the end of the pipe (which connects to the chamber) together with some turbulence I'm fine. It may not be your case.
Hi,

You use nutwall function in the LES calculation. What model do you use?

Sincerely,
Guanjiang
guanjiang.chen is offline   Reply With Quote

Reply

Tags
les, nut, wall function

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
LES inlet boundary conditions hz283 OpenFOAM 17 December 29, 2019 17:59
My radial inflow turbine Abo Anas CFX 27 May 11, 2018 02:44
CFD analaysis of Pelton turbine amodpanthee CFX 31 April 19, 2018 19:02
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 09:44
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32


All times are GMT -4. The time now is 19:39.