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

Query regarding law of wall for temperature?

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By aviral2000

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 28, 2019, 11:46
Default Query regarding law of wall for temperature?
  #1
New Member
 
Delhi
Join Date: Jul 2017
Posts: 10
Rep Power: 9
aviral2000 is on a distinguished road
Hello everyone!

I am slightly confused about the validity of the law of temperature. I tried to go over literature and I couldn't find a concrete answer. I was wondering if the law of wall for temperature shifts when we can change the wall temperature. I ran RANS simulations in an in-house code for a flat plate case with different wall temperature (isothermal boundary condition). I selected Twall/Tfreestream as 1.03, 1.04, 1.05, 1.1. I observed that T+ vs y+ was varying for both SA and k-w SST turbulence models. I wondering if anyone ran a similar simulation in any other code. I want to check if this behavior is expected or not. It would be great to have the opinion of the CFD community about it.

Thank you very much
aviral2000 is offline   Reply With Quote

Old   March 28, 2019, 12:59
Default
  #2
New Member
 
Delhi
Join Date: Jul 2017
Posts: 10
Rep Power: 9
aviral2000 is on a distinguished road
It would be great if someone could give me their opinion regarding it from their experience or relevant literature.
aviral2000 is offline   Reply With Quote

Old   March 29, 2019, 01:41
Default
  #3
Senior Member
 
Lane Carasik
Join Date: Aug 2014
Posts: 692
Rep Power: 15
lcarasik is on a distinguished road
Quote:
Originally Posted by aviral2000 View Post
It would be great if someone could give me their opinion regarding it from their experience or relevant literature.
You should probably give people more than a few hours to respond to your topic.
lcarasik is offline   Reply With Quote

Old   March 29, 2019, 04:55
Default
  #4
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,191
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
Also, show us what you mean by varying (some plots would be useful), how did you define your T+ and y+, your grid parameters, what other bc you used for your case, the properties you used (are you sure the Pr number is the same?).

As it is not expected in general, the first thing I can think of is that your grid is in a gray area, and as you changed the temperature difference, the y+ range chenaged as well. You should be aware that these WF are not perfect, they typically have some variability if you change your y+.
sbaffini is offline   Reply With Quote

Old   March 29, 2019, 08:53
Default
  #5
New Member
 
Delhi
Join Date: Jul 2017
Posts: 10
Rep Power: 9
aviral2000 is on a distinguished road
Quote:
Originally Posted by lcarasik View Post
You should probably give people more than a few hours to respond to your topic.
Sorry if it caused any inconvenience. In my second reply to the thread, I was appending my original question. My previous message might have focussed only to those who ran this kind of simulation. Anyway, next time, I will just edit the main question.

Keeping this un-necessary discussion aside. Can you please give any technical advice to my question?
Thank you!
aviral2000 is offline   Reply With Quote

Old   March 29, 2019, 09:04
Default
  #6
New Member
 
Delhi
Join Date: Jul 2017
Posts: 10
Rep Power: 9
aviral2000 is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Also, show us what you mean by varying (some plots would be useful), how did you define your T+ and y+, your grid parameters, what other bc you used for your case, the properties you used (are you sure the Pr number is the same?).

As it is not expected in general, the first thing I can think of is that your grid is in a gray area, and as you changed the temperature difference, the y+ range changed as well. You should be aware that these WF are not perfect, they typically have some variability if you change your y+.
Hi! Thanks for your reply.

I defined my variables as:

y^+ = \frac{y u_{\tau}}{\nu}
where, I take \nu is the freestream value of viscosity.

T^+ = \frac{T_w - T}{T_{\tau}}

where {T_{\tau}} = \frac{q_w}{\rho u c_p}

and T_w is the wall temperature,
\rho and u are freestream density and velocity respectively. c_p is specific heat at constant pressure and q_w is wall heat flux.

I should mention that this dependence on wall temperature is not due to the absolute value of Twall but the difference with the freestream temperature (T_w - T_{\infty}). For low temperature difference, this behavior is obtained which I find weird as most of the experimental data were obtained with low temperature difference so as to do a constant parametric (with respect to T) measurements.

The results are shown as follows and can be found here: https://drive.google.com/file/d/1wma...ew?usp=sharing

Note that here Twall = 1.03 in the legend is basically Twall = 1.03*T_freestream. Where T_freestream is 300 K for all the simulations.

The prandtl number (=0.71) and c_p are constant with temperature. These are not wall function simulations, they are low-Reynolds number turbulence model simulations with y+ < 0.1. I did mesh convergence study and the result is not varying. It is important to tell here that the U+ vs y+ profile is almost identical to the analytic one for all wall temperatures.

I understand your point about grid variability, however I am not running wall functions. It's the turbulence model solved to the wall with y+ < 0.1. The BC I used for the wall is isothermal wall. The case is exactly the 2-D flat plate case in NASA turbulence modeling resource. Unfortunately, they don't provide heat transfer data for validation. The velocity profile shows good agreement with CFL3D data.
aviral2000 is offline   Reply With Quote

Old   March 29, 2019, 09:47
Default
  #7
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,191
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
Let us assume that you are computing a constant density flow, otherwise things get much more involved (or, in any case, let us assume that density variations are negligible).

Let us also assume that the viscosity variation with temperature is negligible, as this is the case under which classical law of the walls actually hold.

In this scenario, your definition of T_\tau seems wrong to me, as it requires the local u_\tau and not freestream value of u at the denominator.

Second, how do you compute q_w?

In any case, my experience is with the channel flow where, using a temperature bc instead of a fixed flux one, I also experienced non collapsing profiles, but that was for different Re numbers, which might be a different scenario.
sbaffini is offline   Reply With Quote

Old   March 29, 2019, 09:56
Default
  #8
New Member
 
Delhi
Join Date: Jul 2017
Posts: 10
Rep Power: 9
aviral2000 is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Let us assume that you are computing a constant density flow, otherwise things get much more involved (or, in any case, let us assume that density variations are negligible).

Let us also assume that the viscosity variation with temperature is negligible, as this is the case under which classical law of the walls actually hold.

In this scenario, your definition of T_\tau seems wrong to me, as it requires the local u_\tau and not freestream value of u at the denominator.

Second, how do you compute q_w?
Really sorry! That was a typo. The definition of T_{\tau} is indeed as follows,
{T_{\tau}} = \frac{q_w}{\rho u_{\tau} c_p}

I computed q_w as q_w = c_p ( \frac{\mu_w}{Pr} + \frac{\mu_{t,w}}{Pr_t} ) \frac{d T}{d y} from the CFD solution where \frac{d T}{d y} = \frac{T_w - T_1}{\Delta y} where T_1 is the temperature at the first cell and \Delta y is the wall distance. I checked my Stanton number results and they match with the empirical solutions. So, I guess q_w computation (even though it is first order), is not the reason. I also tried green-gauss to get the value of temperature gradient, however, results are still the same. Hmm, I realised that I am using viscosity at the wall for q_w computation. This seems correct with the physics!

Last edited by aviral2000; March 29, 2019 at 12:39.
aviral2000 is offline   Reply With Quote

Old   March 29, 2019, 14:53
Default
  #9
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,191
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
Everything seems legit to me. As I said, I found a similar behavior when changing the Re number, for fixed wall temperature, in a periodic channel flow. Maybe the two things are related.

Yet, the fact that it gets worst with decreasing temperature difference, which is at your T+ numerator, is suspicious... maybe there is a utau dependence from the temperature difference?
sbaffini is offline   Reply With Quote

Old   March 29, 2019, 15:37
Default
  #10
New Member
 
Delhi
Join Date: Jul 2017
Posts: 10
Rep Power: 9
aviral2000 is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Everything seems legit to me. As I said, I found a similar behavior when changing the Re number, for fixed wall temperature, in a periodic channel flow. Maybe the two things are related.

Yet, the fact that it gets worst with decreasing temperature difference, which is at your T+ numerator, is suspicious... maybe there is a utau dependence from the temperature difference?
Yes, that's suspicious. The problem is that these parametric studies are not published usually. With the change in Reynolds number, I get collapsing profile in log region and then the profile diverges. This behavior is expected as the length of log region depends on the Reynolds number of the flow. Do you observe similar behavior with your simulations?

Hmm, I don't think that there is u_{\tau} dependence as the velocity profile and skin-friction results matches perfectly. As of the current moment, my guess is the computation of q_w is not accurate or probably for the specified conditions i.e. Mach number (0.2 in my case), this behavior is expected. I am testing other mach numbers to understand this behaviour. I expect the q_w to change slightly with Mach number. Do you think this is expected?

Thanks for suggestions.
jadidi.cfd likes this.
aviral2000 is offline   Reply With Quote

Old   March 30, 2019, 14:04
Default
  #11
New Member
 
Delhi
Join Date: Jul 2017
Posts: 10
Rep Power: 9
aviral2000 is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Everything seems legit to me. As I said, I found a similar behavior when changing the Re number, for fixed wall temperature, in a periodic channel flow. Maybe the two things are related.

Yet, the fact that it gets worst with decreasing temperature difference, which is at your T+ numerator, is suspicious... maybe there is a utau dependence from the temperature difference?
I think I found the reason. There is no error in my extraction of plots and temperature profile. I believe the issue is with the viscous dissipation term in the energy equation. As I am using a compressible solver, at higher mach number and higher values of \Delta T = T_w - T_{\infty} i.e. at higher Eckert number, the viscous dissipation is significant enough to affect the temperature distribution in the boundary layer. I initially found it strange as I am still within the subsonic regime (M = 0.2) when I observe the deviation in the temperature profile. The reason my Stanton number distribution matched nicely with empirical results is that I was using the following definition,

St = \frac{q_w}{\rho u c_p (T_w - T_r)} where T_r is the adiabatic wall temperature. The use of this definition makes my Stanton number distribution quite independent of Mach number, however, the temperature distribution in the boundary layer is still affected.

At lower mach number (M = 0.1), the T+ vs y+ curve for different wall temperature were matching perfectly. I also tried removing the viscous dissipation term from the RANS solver and observed that even though my St number result is different now (expected behaviour), the temperature profile for different mach numbers are matching nicely.

Most of the T+ vs y+ curves and studies are for incompressible solver which does not model viscous dissipation term. Also, for higher compressible flows, Van-Driest or other similar transformations are used to convert to velocity and temperature profile to analytic ones. The case I was doing was in sort of a gray area as I was using a compressible solver for a low-speed flow without using these transformations for the velocity profile.

That's my opinion regarding the issue. It will be great to have an opinion from you and other members of CFD community regarding it. Feel free to point out any mistakes in my reasoning or suggestions for additional verification tests.

Thank you!
aviral2000 is offline   Reply With Quote

Old   March 30, 2019, 14:19
Default
  #12
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,191
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
Quote:
Originally Posted by aviral2000 View Post
At lower mach number (M = 0.1), the T+ vs y+ curve for different wall temperature were matching perfectly. I also tried removing the viscous dissipation term from the RANS solver and observed that even though my St number result is different now (expected behaviour), the temperature profile for different mach numbers are matching nicely.
Regardless of my opinion, I think this is pretty much your answer to the question.

In more general terms, I would conclude that the main mistake was in applying an incompressible law to a compressible case. Thus, besides the viscous dissipation term, a varying density and, in case, viscosity, might invalidate it as well.
sbaffini is offline   Reply With Quote

Reply

Tags
heat transfer, law of wall, temperature, turbulence modeling


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
Wall treatment with OpenFOAM roby OpenFOAM Running, Solving & CFD 48 May 28, 2021 12:38
[Commercial meshers] Fluent3DMeshToFoam simvun OpenFOAM Meshing & Mesh Conversion 50 January 19, 2020 16:33
Centrifugal fan j0hnny CFX 13 October 1, 2019 14:55
velocity BC for pourous wall (pressure dependendend darcy law) StefanW OpenFOAM Running, Solving & CFD 1 May 29, 2017 08:53
Wall Function Approach - log law Paul CFX 2 June 30, 2005 07:06


All times are GMT -4. The time now is 23:25.