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

thermal wall function calculation?

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

Like Tree1Likes
  • 1 Post By LuckyTran

Reply
 
LinkBack Thread Tools Display Modes
Old   November 11, 2012, 10:56
Default thermal wall function calculation?
  #1
Senior Member
 
François Grégoire
Join Date: Jan 2010
Location: Laval University, Canada
Posts: 325
Rep Power: 7
macfly is on a distinguished road
Hi,

I'm trying to understand how exactly Fluent calculates T^*. From the procedure explained in the attached pdf (above Eqn. 4-285), here is what I understand:

1. Pr is determined from fluid conditions.

2. from the procedure, "given the molecular Prandtl number, the thermal sublayer thickness, y_T^*, is computed from the intersection of the linear and logarithmic profiles, and stored." => does it mean that egality between linear and logarithmic profiles of (4-283) is posed, and this egality is solved for y^*? i.e.:
Pr y^*=Pr_t\left[\frac{1}{\kappa}ln(Ey^*)+P\right] solved for y^* (those are the simplified profiles for incompressible flow).

3. the y^* determined from the above egality is the thermal sublayer thickness: y^*=y^*_T. Then, this y_T^* is compared against the actual y^* value (from turbulence model) and either the linear or logarithmic thermal profile of (4-283) is choosen for the calculation of T^*.

Am I right, is it how T^* is determined? I'm not sure if my step 2. above is right, please correct me if I'm wrong.

Regards,
François
Attached Files
File Type: zip twf.zip (79.1 KB, 14 views)

Last edited by macfly; November 11, 2012 at 16:26.
macfly is offline   Reply With Quote

Old   November 13, 2012, 20:16
Default
  #2
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 497
Rep Power: 10
LuckyTran is on a distinguished road
You are pretty much correct. Except that you may still be slightly confused on step 2.

Once the thermal properties are known, the egality of
y_T^* is computed and then stored. This means equation the linear and log laws, and solving for the value of y^* at which they intersect. Note that this special intersection is then called y_T^*. This value is only computed once. The thermal wall functions only needs the y* value, but it needs a way of determining which law to use.

During the computation, the local
y^* is compared to the pre-calculated y_T^* value to determine which law to use.

The thermal wall functions are identical to the velocity wall functions. However, for the velocity functions we know beforehand the intersection of the linear and log regions since these are independent of fluid properties. The thermal law of the wall instead can only be computed after all thermal boundary conditions have been applied.

Last edited by LuckyTran; November 13, 2012 at 21:10.
LuckyTran is offline   Reply With Quote

Old   November 13, 2012, 20:59
Default
  #3
Senior Member
 
François Grégoire
Join Date: Jan 2010
Location: Laval University, Canada
Posts: 325
Rep Power: 7
macfly is on a distinguished road
Quote:
y^*T
You meant y^*_T.

Quote:
This value is only computed once.
Computed once at the beginning of the iteration process, right?


About step 2, I'm still under the impression that in order to find where the 2 profiles intersect, a clever way is to equate their equations. What am I missing? Fluent equates the profiles and performs some iteration process (like narrowing the right y^* value starting with a low and high value...)?
macfly is offline   Reply With Quote

Old   November 13, 2012, 21:08
Default
  #4
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 497
Rep Power: 10
LuckyTran is on a distinguished road
Quote:
Originally Posted by macfly View Post
You meant y^*_T.


Computed once at the beginning of the iteration process, right?


About step 2, I'm still under the impression that in order to find where the 2 profiles intersect, a clever way is to equate their equations. What am I missing? Fluent equates the profiles and performs some iteration process (like narrowing the right y^* value starting with a low and high value...)?
Thanks, y_T^*, I'm not used to typing equations.

Correct, computed only once. In other words, it is not recomputed on subsequent iterations. This computation is probably done during the initialization step.

You must equate the equations,that is the only way to do it, I would not describe that as a clever equating. I imagine Fluent calculates it by some iterative process like you described, but ignoring the details of how it is calculated, it is a simple find the intersection of two curves problem.

Pretty much any algorithm will work. Bisection works. Fixed-point iteration and Newton-Raphson are also guaranteed to work since one of the equation is a line, and that would converge very quickly (Richardson extrapolation can be used to speed up convergence also). I am sure that Fluent has these algorithms pre-packaged in the background somewhere since they are also needed when solving the large systems of equations during the iteration process.
camoesas likes this.

Last edited by LuckyTran; November 13, 2012 at 23:51.
LuckyTran is offline   Reply With Quote

Old   November 13, 2012, 23:36
Default
  #5
Senior Member
 
François Grégoire
Join Date: Jan 2010
Location: Laval University, Canada
Posts: 325
Rep Power: 7
macfly is on a distinguished road
All right, thanks for the clarifications. Case closed.
macfly is offline   Reply With Quote

Reply

Tags
thermal wall function, turbulence

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 10:31
UDF for wall slipping HFLUENT Fluent UDF and Scheme Programming 0 April 27, 2011 13:03
OpenFOAM static build on Cray XT5 asaijo OpenFOAM Installation 9 April 6, 2011 13:21
ParaView for OF-1.6-ext Chrisi1984 OpenFOAM Installation 0 December 31, 2010 07:42
Quick Question - Wall Function D.Tandra Main CFD Forum 2 March 16, 2004 05:29


All times are GMT -4. The time now is 03:51.