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

Blog Entries

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

Old

Closing on wall functions - part 8: coupled/thin wall boundary conditions

Posted April 27, 2022 at 09:37 by sbaffini (NuTBox)
Updated May 16, 2022 at 12:51 by sbaffini

It might happen that wall functions for temperature or scalars are needed when the assigned boundary condition is not simply the assigned flux or temperature/scalar, but rather a more complex one. One example is when the wall is a coupled one, with either a fluid or a solid on the other side, or it has some thickness, possibly with source terms in it and, say, radiation or convective boundary conditions assigned on the other side, or maybe some other combination of the above.

In all...
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 229 Comments 0 sbaffini is offline Edit Tags
Old

Closing on wall functions - part 7: starting from a profile

Posted April 26, 2022 at 13:12 by sbaffini (NuTBox)
Updated April 27, 2022 at 19:57 by sbaffini

It might happens that one doesn't have a turbulent viscosity profile but actually has just an equilibrium profile for velocity or temperature. More specifically:

T^+ = Pr \left(\frac{{s_T^{-1}}^+}{y^+}\right)  y^+

with its obvious extension to the velocity case. In order to go back to the framework presented here one should notice that:

\frac{d}{dy^+}\left({s_T^{-1}}^+\right) = \frac{1}{1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}}...
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 174 Comments 0 sbaffini is offline Edit Tags
Old

Closing on wall functions - part 6: viscous dissipation

Posted April 24, 2022 at 11:20 by sbaffini (NuTBox)
Updated June 1, 2022 at 13:28 by sbaffini

One thing which is missing in the previous derivations is the viscous dissipation term in the temperature equation. Let's reconsider the initial temperature equation when it is present:

\frac{d}{dy}\left[C_p\left(\frac{\mu}{Pr}+\frac{\mu_t}{Pr_t}\right)\frac{dT}{dy}\right]=F_T - \frac{d}{dy}\left[\left(\mu+\mu_t\right)U\frac{dU}{dy}\right]

A first integration leads to:

\left(\frac{C_p \mu}{Pr}\right)\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)\frac{dT}{dy}=q_w+...
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 231 Comments 0 sbaffini is offline Edit Tags
Old

Closing on wall functions - part 5: testing scripts

Posted April 23, 2022 at 22:26 by sbaffini (NuTBox)
Updated June 3, 2022 at 09:02 by sbaffini

I provide here a set of MATLAB scripts to test all the claims made in the first 4 parts.

The first group of scripts is actually made of functions, that you are not supposed to directly call or modify:
  • muskersp.m: returns \left(\frac{{s_{U,T}^i}^+}{{y^+}^{i+2}}\right), \left(\frac{{p^i}^+}{{y^+}^{i+2}}\right) and {q^i}^+ as shown here. It only works for N up to 0 (constant non equilibrium terms) EDIT: There is an apparently innocuous mistake in the limiting behavior of s,
...
Attached Files
File Type: zip wallfunction.zip (14.2 KB, 181 views)
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 195 Comments 0 sbaffini is offline Edit Tags
Old

Closing on wall functions - part 4: a Musker-Monkewitz wall function

Posted April 23, 2022 at 18:53 by sbaffini (NuTBox)
Updated April 26, 2022 at 17:49 by sbaffini

In the Musker-Monkewitz wall function the following assumption is made on the turbulent viscosity profile:

\frac{\mu_t}{\mu} = \frac{\left(\kappa y^+\right)^3}{\left( \kappa y^+\right)^2+\left(\kappa a \right)^3-\left(\kappa a \right)^2}

where \kappa is the von Karman constant and a is a constant that specifies the y^+ for which \frac{\mu_t}{\mu} =1 but, in practical terms has the same role of y_v^+ in the standard wall function of the previous post....
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 304 Comments 0 sbaffini is offline Edit Tags

All times are GMT -4. The time now is 09:21.