k-omega-SST: what happens when nut is set to 'calculated' at the wall for y+<1

sapujapu · May 11, 2020, 17:09

Hello everyone,

I'm trying to model a heat transfer problem and therefore want to pay special attention to my boundary layer treatment.

The plan is to resolve the boundary layer with prisms so that y+<1 and therefore not rely on wall functions. I want to use the k-omega SST turbulence model as the literature shows that it performs well in my problem and because it's capable of fully resolving the boundary layer.

My original plan was to set the following BC at the wall:

k: kLowReWallFunction

omega: omegaWallFunction

nut: nutLowReWallFunctions

However, I am quite unsure about the BC for nut. From the OF BC Guide, I know that nutLowReWallFunction simply sets nut to zero. Intuitively, this makes sense to me, as there are no turbulent fluctuations at the wall and therefore zero turbulent viscosity. Please correct me though if I'm wrong here.

A colleague of mine told me that he usually uses 'calculated' for the nut wall BC. I don't really understand the point behind this though: if this BC is used, how exactly is nut calculated and how is this result going to look like? Is it simply calculated as $\nu_t = a_1 \frac{k}{max(a_1 \omega , b_1 , F_{23}S)}$ ? (Equation for nut fromk-omega-SST model)

If yes, it looks like it does go to zero because the turbulent kinetic energy k should be zero at the wall (at least according to my intuition).

I hope my question is somewhat clear, it would be really cool if someone could help me out there. Thanks!

geth03 · June 3, 2020, 06:05

i also use the k-omega-sst turbulence model.
for walls, i use nutkWallFunction.
nut takes really small values at walls, like you assumed, but in my case
never is 0. (e-5 to e-10)

your equation for nut is correct.
here is the link to the model:
https://www.openfoam.com/documentati...omega-sst.html

cryabroad · June 4, 2020, 04:44

"calculated" means to calculate the quantity based on the way it should be calculated, instead of doing something to it. So you are right, if you use "calculated" on nut, it means you calculate nut based on how it should be calculated, i.e., from the model itself.

I'm not so sure about having such a fine mesh for RANS though. To me, the whole idea of using RANS is to reduce computational resources. So, the mesh doesn't have to be very fine, and wall models are used to give reasonable results for a coarse mesh (especially towards the wall). Maybe you are having troubles getting reasonable results when using wall models? If that's the case do you mind share the problems? I'm kind of interested.

Ship Designer · June 13, 2020, 08:41

Hello,

In the formula for calculating nut, what are the symbols F_23 and S? Thanks.

May 11, 2020, 17:09	k-omega-SST: what happens when nut is set to 'calculated' at the wall for y+<1	#1
sapujapu New Member Join Date: Feb 2020 Posts: 10 Rep Power: 6	Hello everyone, I'm trying to model a heat transfer problem and therefore want to pay special attention to my boundary layer treatment. The plan is to resolve the boundary layer with prisms so that y+<1 and therefore not rely on wall functions. I want to use the k-omega SST turbulence model as the literature shows that it performs well in my problem and because it's capable of fully resolving the boundary layer. My original plan was to set the following BC at the wall: k: kLowReWallFunction omega: omegaWallFunction nut: nutLowReWallFunctions However, I am quite unsure about the BC for nut. From the OF BC Guide, I know that nutLowReWallFunction simply sets nut to zero. Intuitively, this makes sense to me, as there are no turbulent fluctuations at the wall and therefore zero turbulent viscosity. Please correct me though if I'm wrong here. A colleague of mine told me that he usually uses 'calculated' for the nut wall BC. I don't really understand the point behind this though: if this BC is used, how exactly is nut calculated and how is this result going to look like? Is it simply calculated as $\nu_t = a_1 \frac{k}{max(a_1 \omega , b_1 , F_{23}S)}$ ? (Equation for nut fromk-omega-SST model) If yes, it looks like it does go to zero because the turbulent kinetic energy k should be zero at the wall (at least according to my intuition). I hope my question is somewhat clear, it would be really cool if someone could help me out there. Thanks!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
does Hyperthreading affect the application of UDF?	SJSW	Fluent UDF and Scheme Programming	11	October 10, 2018 22:28
Implementation of new omega and nut wall functions	zordiack	OpenFOAM Programming & Development	4	December 6, 2016 15:17
Boundary Conditions for k omega SST	dancfd	OpenFOAM Pre-Processing	0	June 9, 2011 23:25
Env variable not set	gruber2	OpenFOAM Installation	5	December 30, 2005 04:27
Multicomponent fluid	Andrea	CFX	2	October 11, 2004 05:12

June 3, 2020, 06:05		#2
geth03 Senior Member Join Date: Dec 2019 Location: Cologne, Germany Posts: 355 Rep Power: 8	i also use the k-omega-sst turbulence model. for walls, i use nutkWallFunction. nut takes really small values at walls, like you assumed, but in my case never is 0. (e-5 to e-10) your equation for nut is correct. here is the link to the model: https://www.openfoam.com/documentati...omega-sst.html

June 4, 2020, 04:44		#3
cryabroad Senior Member Ruiyan Chen Join Date: Jul 2016 Location: Hangzhou, China Posts: 162 Rep Power: 9	"calculated" means to calculate the quantity based on the way it should be calculated, instead of doing something to it. So you are right, if you use "calculated" on nut, it means you calculate nut based on how it should be calculated, i.e., from the model itself. I'm not so sure about having such a fine mesh for RANS though. To me, the whole idea of using RANS is to reduce computational resources. So, the mesh doesn't have to be very fine, and wall models are used to give reasonable results for a coarse mesh (especially towards the wall). Maybe you are having troubles getting reasonable results when using wall models? If that's the case do you mind share the problems? I'm kind of interested.

June 13, 2020, 08:41		#4
Ship Designer Senior Member Claudio Boezio Join Date: May 2020 Location: Europe Posts: 137 Rep Power: 6	Hello, In the formula for calculating nut, what are the symbols F_23 and S? Thanks.