|
[Sponsors] |
Low y+ value for enhanced wall treatement. |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 26, 2019, 04:24 |
Low y+ value for enhanced wall treatement.
|
#1 |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
Hello,
y+~1 is offered for enhanced wall treatment, I wonder using very low y+ value (for example 0.2) would be a problem for accuracy? Thanks in advance. |
|
August 26, 2019, 04:41 |
|
#2 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
The criterium is y+ up to ~1. So y+ values lower than 1 are usually not a problem for wall-resolving boundary layer treatment in turbulence models.
You just have to pay attention with very low wall distances or very high aspect ratio cells near the wall. These can result in problems with numerical accuracy, and some solvers (e.g. Fluent) offer special numerical formulations for these cases. Last edited by flotus1; August 26, 2019 at 11:26. |
|
August 26, 2019, 10:05 |
|
#3 | |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
Quote:
|
||
August 26, 2019, 11:24 |
|
#4 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
It has less to do with the value of y+ and more with the absolute wall distance. There is no clear threshold for either of them.
When you import a mesh in fluent with very small cells at the boundaries (compared to the domain size?) or very high aspect ratio cells, you will get a warning and a hint to the counter-measures. Some solvers (e.g. CCM+) just have a limiter for the smallest -dimensional- wall distance. |
|
August 27, 2019, 00:22 |
|
#5 |
Senior Member
Ashwani
Join Date: Sep 2013
Location: Hyderabad
Posts: 154
Rep Power: 12 |
You can use y+=0.2. It would not be much of a problem. However, as mentioned just have to be careful at very low y+ and high aspect ratio cell.
Also, you should check your y+ during the solution evolution and in the final solution, as y+ is a post-processing parameter. |
|
September 2, 2019, 12:01 |
|
#6 |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
Reduced the number of elements along wall normal direction. This caused increase in y+ (still lower than 1) and aspect ratio is halved but my heat transfer results didn't changed. So neither aspect ratio nor y+ was my problem. I am still digging.
|
|
September 2, 2019, 12:53 |
|
#7 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
||
September 2, 2019, 13:47 |
|
#8 | |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
Quote:
Last edited by CaptainCombo; September 2, 2019 at 16:01. |
||
September 6, 2019, 09:25 |
|
#9 |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
If I swtich near wall tratement from "enhanced wall tratemnet" to "scaleable wall function" without changing any other thing. I am getting true results. Can anyone give me a clue, what is the reason could be?
|
|
September 6, 2019, 10:34 |
|
#10 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
Two or more errors cancelling each other out?
|
|
September 7, 2019, 17:30 |
|
#11 |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
I tried litterally everything to get true result by using enhanced wall treatement but all my efforts didn't change anything.
- Tried 40 layer for inflation, in order to put enough number of elements to boundary layer. - Used elements with the aspect ratio 1 and y+ was ~1.1, this cost me 6M elements such a simple problem but it didn't work. - Investigated the effect of turbulence intensity at the inlet by changing it between 1-5%, I saw that turbulence intensity does not have big influence. - Used enhanced wall treatement options for pressure gradient and thermal effects but they didn't cause any difference. I had other attemps but it does not worth to mention here. As a result my numerical simulation gives local Nu number near the end of the channel 30% higher than Gnielinski correlation. I checked that thermal and hydraulic boundary layers becomes fully developed before the end of the channel. I shared case and data of one of my tries. Maby someone see something I don't see. Case: https://drive.google.com/open?id=1ul...DSl8AUmZ2byfNe Data:https://drive.google.com/open?id=1yA...28Vp02BLtqdRUv |
|
September 7, 2019, 21:24 |
|
#12 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
You are really not giving us much to work with here.
1) What exactly are you trying to simulate? 2) How exactly did you do it - describe in a few words, add images of computational domain and mesh. Not everybody has Ansys Fluent readily available and has the time to go through the setup. And some might even be reluctant to download files from unknown sources. |
|
September 8, 2019, 07:24 |
|
#13 |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
I am trying to calculate convective heat transfer from the walls of straight channel under turbulent flow condition (Re = 10000). I am using relizable k-e turbulence model. The channel is 1.5m long and the problem is assumed to be two dimensional. Constant heat flux is applied at walls, distance between upper and lower walls is 0.01m. Uniform velocity and temperature profile is defined at inlet however I checked thermal and hydraulic boundary layers reach fully developed state where I calculate Nusselt number. I am sure about my calculation method of Nusselt number because I solved similar problem for laminar flow condition; my numerical simulation gives almost same results with analytical methods.
Here are mesh and y+ for one of my attemps; |
|
September 9, 2019, 06:46 |
|
#14 | |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
Quote:
Is there any particular reason why you did not use periodic boundary conditions for this study? It would eliminate any chance of the flow not being fully developed, while also cutting down computational cost drastically. A mapped mesh would also make it easier on the eye. I am not too firm on which turbulence model is the right one for heat transfer, but IIRC k-omega SST is the default recommendation by Ansys. |
||
September 10, 2019, 03:17 |
|
#15 |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
In this (https://www.sciencedirect.com/scienc...09250995000354) paper they compared various low Re turbulence models. So, I ran simulations for turbulent pipe flow for the same Reynolds number as in the paper (Re = 7000). I applied the streamwise periodic boundary condition to be sure about flow development. In addition, I added the constant heat flux boundary condition to the same simulations in order to calculate Nu. Here is the velocity profile in the viscous sublayer;
post a picture Since bulk temperature does not change for all turbulence models, I think I need to focus on the viscous sublayer, is that true? My simulation with using Lam & Bremhorst gives closest heat transfer predictions. Nu is only higher 1.96% than empirical calculation. However, Realizable k-e overestimates 17.91% and SST k-w intermittency model overestimates 24%. @flotus1 What is IIRC, I assumed you mean intermitency transition SST k-w. |
|
September 10, 2019, 04:56 |
|
#16 | |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,399
Rep Power: 46 |
Sorry for the confusion. IIRC means "if I recall correctly", a common abbreviation to save precious time when communicating on the interwebs. Didn't really work this time
Quote:
Heat transfer isn't really my strong suit, maybe someone else can chime in? |
||
September 10, 2019, 09:56 |
|
#17 |
Senior Member
Blanco
Join Date: Mar 2009
Location: Torino, Italy
Posts: 193
Rep Power: 17 |
Hi,
are you simulating a constant density flow? Because that can change significantly the problem in my experience. In addition, are you using the friction factor calculated directly using the Moody's diagram, w/ a known pipe roughness? |
|
September 11, 2019, 06:42 |
|
#18 | |
Member
Join Date: May 2012
Posts: 30
Rep Power: 13 |
Quote:
Hello, I am using pressure-based solver and constant fluid properties. I didn't specify surface roughness, Fluent uses smooth surface by default. After I saw your post I calculated friction factor (f) for low Reynolds number model of Lam & Bremhorst and SST k-omega just as it is happening in the heat transfer low-re-ke model of LB predicts friction factor better than SST k-w. I think the reason behind the wrong estimation of heat transfer also causing wrong surface friction. By the way are you author of this paper; https://www.sciencedirect.com/scienc...90072910000803 |
||
September 12, 2019, 09:03 |
|
#19 | |
Senior Member
Blanco
Join Date: Mar 2009
Location: Torino, Italy
Posts: 193
Rep Power: 17 |
Ok, in addition there's also the turbulent Prandtl number. It's typically set at 0.9 which is "reasonable" for air (and maybe water if I remember correctly, or at least it seems reasonable for most flows of these two fluids), but if you're simulating some exotic fluid I would expect different values. Unfortunately this as it is still a research area.
As far as the article is concerned, no I'm not in that pubblication, but you're guess on my name is still very close to reality Quote:
|
||
September 13, 2019, 03:55 |
|
#20 |
Senior Member
Andrea
Join Date: Feb 2012
Location: Leeds, UK
Posts: 179
Rep Power: 16 |
Just to add that not only the value of the turbulent Prandtl number is tricky to evaluate, but the assumption of a constant Pr_t is often a major simplification with little physical justification. Actually, the SGDH assumption for the turbulent heat flux is indeed a major simplification with no sound physical basis in most cases - we just use it because it's simple and robust and often it is the only option that we have in commercial codes
It's interesting to observe that the SGDH with Pr_t~1 breaks down even with non-exotic unity Prandtl number fluids in very basic flows in natural and mixed convection, such as Rayleigh-Benard convection and unstable stratified mixed convection in channels, as reported here: https://www.sciencedirect.com/scienc...29549319302407 One quick workaround to relax the constant Pr_t assumption in the SGDH context is to use a local correlation - see, for instance, section 3 here: https://core.ac.uk/download/pdf/82053898.pdf |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[mesh manipulation] mergeMeshes problem | Attesz | OpenFOAM Meshing & Mesh Conversion | 3 | July 29, 2015 04:15 |
Wall boundary conditions for low Re models | callumso | OpenFOAM Running, Solving & CFD | 10 | June 12, 2013 17:16 |
Boundary Layer roughness/ low reynolds wall treatment | Luigi_ | STAR-CCM+ | 1 | March 14, 2012 08:40 |
__ wall function for low yPlus __ | maysmech | OpenFOAM | 14 | January 18, 2011 09:18 |
wall function for low Re | M. Essuri | FLUENT | 2 | December 18, 2006 03:11 |