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

Closing on wall functions - part 6: viscous dissipation

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

Rate this Entry

Closing on wall functions - part 6: viscous dissipation

Posted April 24, 2022 at 10:20 by sbaffini
Updated June 1, 2022 at 12: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+ yG_T - \mu\left(1+\frac{\mu_t}{\mu}\right)U\frac{dU}{dy}

where yG_T is the integral of F_T and G_T has the same form we assumed for F_T, just with its coefficients divided by i+1. Then a further integration leads to:

T\left(y\right) = T_w + \left(\frac{Pr}{\mu C_p}\right) \int_{0}^{y}\frac{q_w+ zG_T}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}dz - \left(\frac{Pr}{C_p}\right) \int_{0}^{y}\frac{\left(1+\frac{\mu_t}{\mu}\right)}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}U\frac{dU}{dz}dz

The first integral in the equation above has been the subject of the previous posts, while the neglected viscous dissipation term is the last term in the equation. Now, there are two ways which could be used to proceed. A first one involves invoking the velocity equation and transforming the viscous dissipation term as follows:

D\left(y\right)=-\left(\frac{Pr}{C_p}\right) \int_{0}^{y}\frac{\left(1+\frac{\mu_t}{\mu}\right)}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}U\frac{dU}{dz}dz = -\left(\frac{Pr}{\mu C_p}\right)\int_{0}^{y}\frac{\left(\tau_w+zG_U\right)U}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}dz

where G_U has the same form as G_T, but with the velocity coefficients. At this point, then, one can note that the integral has the exact same form of the first one, but with the additional velocity U in it. I have not attempted to solve it with none of the models presented but, even if doable, I expect high levels of cumbersomeness.

One could approximate the integral by assuming a constant, representative, value for U, say, \widetilde{U} and take it out of the integral. In this case, while one should still be careful in mixing the velocity numerator with the temperature denominator in the integrand, the resulting integral would be formally identical to the first one, which we have solved in previous posts.

A second, preferred way to proceed would instead first recognize that the integral above can be transformed as:

D\left(y\right)=-\left(\frac{Pr}{C_p}\right) \int_{0}^{y}\frac{\left(1+\frac{\mu_t}{\mu}\right)}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}U\frac{dU}{dz}dz = -\left(\frac{Pr}{C_p}\right) \int_{0}^{y}\frac{\left(1+\frac{\mu_t}{\mu}\right)}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}\frac{d}{dz}\left(\frac{U^2}{2}\right)dz

Hence integration by parts leads to the following result:

D\left(y\right)=-\left(\frac{Pr}{2C_p}\right)\left\{U^2\frac{\left(1+\frac{\mu_t}{\mu}\right)}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}-U_w^2-\int_{0}^{y}U^2\frac{d}{dz}\left[\frac{\left(1+\frac{\mu_t}{\mu}\right)}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}\right]dz\right\}

At this point we are still left with an even more cumbersome integral to evaluate, if possible at all. However, this integral is exactly 0 in two special cases: fully laminar ones (\mu_t/\mu=0) and when Pr/Pr_t=1. Another advanatge of this form is that if we now assume a constant, representative, value for U, say, \widetilde{U} and take it out of the integral, we are left with no more work to do and just obtain:

D\left(y\right)=-\left(\frac{Pr}{2C_p}\right)\left[\left(U^2-\widetilde{U}^2\right)\frac{\left(1+\frac{\mu_t}{\mu}\right)}{\left(1+\frac{Pr}{Pr_t}\frac{\mu_t}{\mu}\right)}-\left(U_w^2-\widetilde{U}^2\right)\right]

which also always works in the cases where the integral is exactly 0. Of course, we are now left with determining \widetilde{U}^2, but this seems a more reasonable task, also because it must be representative just for the region where derivative in the integrand is not 0 and not the full y length. A value for \widetilde{U} commonly used with standard wall functions is the velocity at {y_T}^+. By extension, for the Musker-Monkewitz wall function one could use the velocity at y^+=\widetilde{a}, the profile constant modified by the Pr/Pr_t ratio (see here). In more general cases, it seems that a good generalization for \widetilde{U} is the velocity at the point where \mu_t/\mu=Pr_t/Pr, which is where the term under the derivative reaches half its maximum excursion.

However, this whole presentation is an approximate solution to the original problem for the viscous dissipation, that's why I left it out of the general discussion of the previous posts.
Posted in Uncategorized
Views 736 Comments 0 Edit Tags Email Blog Entry
« Prev     Main     Next »
Total Comments 0

Comments

 

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