CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > Turbulence boundary conditions

Turbulence boundary conditions

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
 
(7 intermediate revisions not shown)
Line 5: Line 5:
==Fully developed turbulent pipe-flow inlet==
==Fully developed turbulent pipe-flow inlet==
-
For fully developed turbulent pipe flow the turbulence inlet properties can be estimated using the model presented by Basse in Table 1 [10].
+
For fully developed turbulent pipe flow the turbulence inlet properties can be estimated using the model presented by Basse in Table 1 of [10].
:<math>l_{AA} = 0.14 \; \kappa_g \times \delta</math>
:<math>l_{AA} = 0.14 \; \kappa_g \times \delta</math>
Line 12: Line 12:
:<math>\epsilon_{AA} = C_{\mu,AA}^\frac{3}{4} \times \frac{k_{AA}^\frac{3}{2}}{l_{AA}}</math>
:<math>\epsilon_{AA} = C_{\mu,AA}^\frac{3}{4} \times \frac{k_{AA}^\frac{3}{2}}{l_{AA}}</math>
-
Where the von Kármán number <math>\kappa_g</math> is a function the friction Reynolds number <math>Re_\tau</math> following [11]:
+
The model parameters <math>\kappa_g</math>, <math>A_g</math>, <math>B_g</math> and <math>C_g</math> can be computed using the following function:
 +
 
 +
:<math>Q(Re_\tau) = a + b \cdot tanh(c \cdot [Re_\tau - d])</math>
 +
 
 +
Where the a, b, c and d constants have been fitted using Princeton Superpipe measurements [2]:
 +
 
 +
{| class="wikitable" style="margin-left:20px; text-align: center;"
 +
|-
 +
! style="width:80px;"|Parameter !! style="width:60px;"|a !! style="width:60px;"|b !! style="width:60px;"|c !! style="width:60px;"|d
 +
|-
 +
| <math>\kappa_g</math> || style="text-align:right;"|−1.18 || style="text-align:right;"|1.52 || style="text-align:right;"|2.15e-4 || style="text-align:right;"|-8786
 +
|-
 +
| <math>A_g</math> || style="text-align:right;"|2.21 || style="text-align:right;"|-0.60 || style="text-align:right;"|3.97e-5 || style="text-align:right;"|11186
 +
|-
 +
| <math>B_g</math> || style="text-align:right;"|1.28 || style="text-align:right;"|-0.32 || style="text-align:right;"|5.85e-5 || style="text-align:right;"|4609
 +
|-
 +
| <math>C_g / \sqrt{Re_\tau}</math> || style="text-align:right;"|1.03 || style="text-align:right;"|-0.91||  style="text-align:right;"|3.30e-5 || style="text-align:right;"|-11755
 +
|}
 +
 
 +
:<math>\delta</math> is the boundary layer thickness, which in fully developed pipe flow is the radius, or half the [[hydraulic diameter]].
 +
 
 +
:<math>\lambda</math> is the friction factor.
-
:<math>\kappa_g(Re_\tau) = −1.18 + 1.52 \cdot tanh(2.15 \cdot 10^{-4} \cdot (Re_\tau + 8785.94))</math>
 
-
:<math></math>
 
-
:<math></math>
 
:<math></math>
:<math></math>
:<math></math>
:<math></math>
Line 60: Line 78:
{{reference-paper|author=[11] Basse, N.T.|year=2023|title=Supplementary Information: An algebraic non-equilibrium turbulence model of the high Reynolds number transition region|rest=https://www.researchgate.net/publication/373108195_Supplementary_Information_An_algebraic_non-equilibrium_turbulence_model_of_the_high_Reynolds_number_transition_region}}
{{reference-paper|author=[11] Basse, N.T.|year=2023|title=Supplementary Information: An algebraic non-equilibrium turbulence model of the high Reynolds number transition region|rest=https://www.researchgate.net/publication/373108195_Supplementary_Information_An_algebraic_non-equilibrium_turbulence_model_of_the_high_Reynolds_number_transition_region}}
 +
 +
{{reference-paper|author=[12] Hultmark M, Vallikivi M, Bailey SCC and Smits AJ.|year=2013|title=Logarithmic scaling
 +
of turbulence in smooth- and rough-wall pipe flow|rest=J. Fluid Mech. 728, 376-395}}

Latest revision as of 23:45, 1 December 2024

This section is under construction, please do not trust the information available here yet

Introduction

Fully developed turbulent pipe-flow inlet

For fully developed turbulent pipe flow the turbulence inlet properties can be estimated using the model presented by Basse in Table 1 of [10].

l_{AA} = 0.14 \; \kappa_g \times \delta
I_{AA} = \sqrt{\left[ B_g + \frac{3}{2} A_g - \frac{8 C_g}{\sqrt[3]{Re_\tau}}\right] \times \frac{\lambda}{8}}
k_{AA} = U_m^2 \; I_{AA}^2
\epsilon_{AA} = C_{\mu,AA}^\frac{3}{4} \times \frac{k_{AA}^\frac{3}{2}}{l_{AA}}

The model parameters \kappa_g, A_g, B_g and C_g can be computed using the following function:

Q(Re_\tau) = a + b \cdot tanh(c \cdot [Re_\tau - d])

Where the a, b, c and d constants have been fitted using Princeton Superpipe measurements [2]:

Parameter a b c d
\kappa_g −1.18 1.52 2.15e-4 -8786
A_g 2.21 -0.60 3.97e-5 11186
B_g 1.28 -0.32 5.85e-5 4609
C_g / \sqrt{Re_\tau} 1.03 -0.913.30e-5 -11755
\delta is the boundary layer thickness, which in fully developed pipe flow is the radius, or half the hydraulic diameter.
\lambda is the friction factor.


Information below is old and will be reformulated

For fully developed duct flow the turbulence intensity at the core can be estimated as [1]:

I = 0.16 \; Re_{d_h}^{-\frac{1}{8}},

where Re_{d_h} is the Reynolds number based on the pipe hydraulic diameter d_h. Additional details on the derivation can be found in [2].

Russo and Basse published a paper [3] where they derive turbulence intensity scaling laws based on CFD simulations and Princeton Superpipe measurements. The turbulence intensity over the pipe area is defined as an arithmetic mean (AM). The measurement-based scaling laws are:

I_{\rm Smooth~pipe~axis} = 0.0550 \; Re^{-0.0407}
I_{\rm Smooth~pipe~area,~AM} = 0.227 \; Re^{-0.100}

Scaling using other turbulence intensity definitions is investigated in [4,5]. Here, it is also found that turbulence intensity scales with the friction factor, both for smooth- and rough-wall pipe flow. Code for an example in [5] can be found in [6]. A high Reynolds number transition in the scaling has been characterized in [7,8]. Turbulence intensity scaling extrapolated to extreme Reynolds numbers is studied in [9].

References

[1] ANSYS, Inc. (2022), "ANSYS Fluent User's Guide, Release R1", Equation (7.71).

[2] Basse, N.T. (2022), "Mind the Gap: Boundary Conditions for Turbulence Modelling", https://www.researchgate.net/publication/359218404_Mind_the_Gap_Boundary_Conditions_for_Turbulence_Modelling.

[3] Russo, F. and Basse, N.T. (2016), "Scaling of turbulence intensity for low-speed flow in smooth pipes", Flow Meas. Instrum., vol. 52, pp. 101–114.

[4] Basse, N.T. (2017), "Turbulence intensity and the friction factor for smooth- and rough-wall pipe flow", Fluids, vol. 2, 30.

[5] Basse, N.T. (2019), "Turbulence intensity scaling: A fugue", Fluids, vol. 4, 180.

[6] Basse, N.T. (2019), "Python code to calculate turbulence intensity based on Reynolds number and surface roughness.", https://www.researchgate.net/publication/336374461_Python_code_to_calculate_turbulence_intensity_based_on_Reynolds_number_and_surface_roughness.

[7] Basse, N.T. (2021), "Scaling of global properties of fluctuating and mean streamwise velocities in pipe flow: Characterization of a high Reynolds number transition region", Physics of Fluids, vol. 33, 065127.

[8] Basse, N.T. (2021), "Scaling of global properties of fluctuating streamwise velocities in pipe flow: Impact of the viscous term", Physics of Fluids, vol. 33, 125109.

[9] Basse, N.T. (2022), "Extrapolation of turbulence intensity scaling to Re_tau>>10^5", Physics of Fluids, vol. 34, 075128.

[10] Basse, N.T. (2023), "An Algebraic Non-Equilibrium Turbulence Model of the High Reynolds Number Transition Region", Water 2023, 15, 3234. https://doi.org/10.3390/w15183234.

[11] Basse, N.T. (2023), "Supplementary Information: An algebraic non-equilibrium turbulence model of the high Reynolds number transition region", https://www.researchgate.net/publication/373108195_Supplementary_Information_An_algebraic_non-equilibrium_turbulence_model_of_the_high_Reynolds_number_transition_region.

[12] Hultmark M, Vallikivi M, Bailey SCC and Smits AJ. (2013), "Logarithmic scaling of turbulence in smooth- and rough-wall pipe flow", J. Fluid Mech. 728, 376-395.

Retrieved from "https://mail.cfd-online.com/Wiki/Turbulence_boundary_conditions"
My wiki