Near-wall treatment for k-omega models

From CFD-Wiki

Revision as of 12:22, 31 October 2011 (view source) AlmostSurelyRob (Talk | contribs) ← Older edit		Latest revision as of 10:40, 31 March 2021 (view source) ChaoLW (Talk | contribs) (the formula for kp is wrong)
(19 intermediate revisions not shown)
Line 1:		Line 1:
-	==Standard wall function approach==	+	{{Turbulence modeling}}
-	There are two possible ways of implementing wall functions in a finite volume code:	+	As described in [[Two equation turbulence models]] low and high reynolds number treatments are possible.
-	* Additional source term in the momentum equations.	+	==Standard wall functions==
-	* Modification of turbulent viscosity in cells adjacent to solid walls.	+	Main page: [[Two equation models#Near-wall treatments| Two equation near-wall treatments]]
-	The source term in the first approach is simply the difference between logarithmic and linear interpolation of velocity gradient multiplied by viscosity (the difference between shear stresses). The second approach does not attempt to reproduce the correct velocity gradient. Instead, turbulent viscosity is modified in such a way as to guarantee the correct shear stress.	+	For <math>k</math> the boundary conditions imposed at the solid boundary are:
		+	:<math>
		+	\begin{matrix}
		+	\frac{\partial k}{\partial n} = 0 & & \frac{\partial \omega}{\partial n} = 0
		+	\end{matrix}
		+	</math>
		+	where <math>n</math> is the normal to the boundary.
-	Using the compact version of log-law	+	Moreover the centroid values in cells adjacent to solid wall are specified as
-	<table width="70%"><tr><td>	+
-	<math>	+
-	\frac{U}{u_\tau} = \frac{1}{\kappa} \ln E y^{+}	+
-	</math></td><td width="5%">(1)</td></tr></table>	+
-	where <math>E=9.8</math> is equivalent to additive constants, and and using <math>\tau_w = \rho u_\tau^2</math> we obtain:	+
	:<math>		:<math>
-	\tau_w = \frac{\rho u_\tau \kappa U}{\ln Ey^{+}},	+	\begin{matrix}
		+	k_p = \frac{u^2_\tau}{\sqrt{C_\mu}},
		+	&&
		+	\omega_p = \frac{u_\tau}{\sqrt{C_\mu}\kappa y_p} = \frac{\sqrt{k_p}}{{C_\mu^{1/4}}\kappa y_p}.
		+	\end{matrix}
	</math>		</math>
-	On the other hand, the linear interpolation for shear stress, rembembering that <math>U|_{y=0}=0</math>, is:	+	In the alternative approach <math>k</math> production terms is modified.
		+	==Automatic wall treatments==
		+	The purpose of automatic wall treatments is to make results insensitive with respect to wall mesh refinement. Many blending approaches have been proposed. The one by Menter takes advantage of the fact that the solution to <math>\omega</math> equations is known for both viscous and log layer
	:<math>		:<math>
-	\tau_w = (\nu_t + \nu)\frac{U_p}{y_p}.	+	\begin{matrix}
		+	\omega_\text{vis} = \frac{6\nu}{\beta y^2} &  \omega_\text{log} = \frac{u_\tau}{C_\mu^{1/4} \kappa y}
		+	\end{matrix}
	</math>		</math>
-		+	where <math>y</math> is the cell centroid distance from the wall. Using this a blending can take the following form:
-	Comparing the above equations we obtain an expression for turbulent viscosity can be obtained:	+
-	<table width="70%"><tr><td>	+
-	<math>	+
-	\nu_t = \nu\left( \frac{y^{+}\kappa}{\ln Ey^{+}} - 1\right).	+
-	</math> </td><td width="5%">(3)</td></tr></table>	+
-	Note that <math>u_\tau</math> has been been incorporated in <math>y^+</math>. The latter remains the only unknown in the equation and has to be estimated for the current velocity field. In the standard approach this cannot be done explicitly and instead an implicit way of obtaining <math>y^{+}</math> has to be employed.	+
-		+
-	After multiplying log law (1) by <math>y_p/\nu</math> and after reorganising some terms we get:	+
-	<table width="70%"><tr><td>	+
-	<math>	+
-	\frac{\kappa U_p y_p}{\nu} = y^{+}\ln{Ey^{+}}.	+
-	</math></td><td width="5%">(3)</td></tr></table>	+
-	This equation can be solved numerically with respect to <math>y^+</math> for example via root searching algorithms e.g. Newton method for specified <math>U_p</math>, <math>y_p</math> and <math>\nu</math>. One iteration in a Newton method for (3) is	+
	:<math>		:<math>
-	y^{+}_{n+1} = \frac{\frac{\kappa U_p y_p}{\nu} + y^{+}_{n}}{1 + \ln E y_n^{+}}.	+	\omega_p = \sqrt{\omega_{\text{vis}}^2 + \omega_{\text{log}}^2},
	</math>		</math>
-	Thus obtained <math>y^{+}</math> is then substituted to (2) Eventually the estimated <math>u_\tau</math> serves also to define the values of turbulent quantities in the cell adjacent to the wall:	+	Note that for low <math>y</math> values the <math>1/y^2</math> will dominate and therefore viscous value of <math>\omega</math> will be reproduced. Conversely, for larger values of <math>y</math>, <math>1/y</math> will be dominant and logarithmic value will be recovered.
		+
		+	Subsequently Menter proposes also blending for friction velocity. Friction velocity for viscous and logarithmic region are:
	:<math>		:<math>
-	k_p = \frac{u_\tau^2}{\sqrt{C_\mu}}	+	\begin{matrix}
		+	u^\text{vis}_\tau = \frac{U}{y^{+}}  & &  u_\tau^\text{log} = \frac{U}{\log E y^{+}}
		+	\end{matrix}
	</math>		</math>
		+	And the blending suggested:
	:<math>		:<math>
-	\omega_p = \frac{\sqrt{k_p}}{{C_\mu^{1/4}}\kappa y_p},	+	u_\tau = \sqrt[4]{(u_\tau^{\text{vis}})^4 + (u_\tau^{\text{log}})^4},
	</math>		</math>
-	which are the values for <math>k</math> and <math>\omega</math> according to Wilcox2006 asymptotic analysis of log layer. These wall functions for <math>k</math> and <math>\omega</math> are the results of the solution of model equation for the logarithmic layer.
-
-	The above methodology is known to produce spurious results in separated flows, where, by definition, <math>u_\tau = 0</math>. Many extension of this approach has been proposed.
-
-
	==FLUENT==		==FLUENT==
Line 68:		Line 64:
	<br>		<br>
-	Reference:	+	== References ==
-	FLUENT 6.2 Documentation, 2006 <br>	+
		+	* {{reference-paper|author=Menter, F., Esch, T.|year=2001|title=Elements of industrial heat transfer predictions|rest='COBEM 2001, 16th Brazilian Congress of Mechanical Engineering.'}}
		+	* {{reference-paper|author=ANSYS|year=2006|title=FLUENT Documentation|rest=''}}

---

Latest revision as of 10:40, 31 March 2021

Turbulence modeling

Turbulence

RANS-based turbulence models Linear eddy viscosity models Algebraic models Cebeci-Smith model Baldwin-Lomax model Johnson-King model A roughness-dependent model One equation models Prandtl's one-equation model Baldwin-Barth model Spalart-Allmaras model Two equation models k-epsilon models Standard k-epsilon model Realisable k-epsilon model RNG k-epsilon model Near-wall treatment k-omega models Wilcox's k-omega model Wilcox's modified k-omega model SST k-omega model Near-wall treatment Realisability issues Kato-Launder modification Durbin's realizability constraint Yap correction Realisability and Schwarz' inequality Nonlinear eddy viscosity models Explicit nonlinear constitutive relation Cubic k-epsilon EARSM v2-f models model model Reynolds stress model (RSM)

Large eddy simulation (LES) Smagorinsky-Lilly model Dynamic subgrid-scale model RNG-LES model Wall-adapting local eddy-viscosity (WALE) model Kinetic energy subgrid-scale model Near-wall treatment for LES models

Detached eddy simulation (DES)

Direct numerical simulation (DNS)

Turbulence near-wall modeling

Turbulence free-stream boundary conditions Turbulence intensity Turbulence length scale

As described in Two equation turbulence models low and high reynolds number treatments are possible.

Contents 1 Standard wall functions 2 Automatic wall treatments 3 FLUENT 4 References

Standard wall functions

Main page: Two equation near-wall treatments

For the boundary conditions imposed at the solid boundary are:

where is the normal to the boundary.

Moreover the centroid values in cells adjacent to solid wall are specified as

In the alternative approach production terms is modified.

Automatic wall treatments

The purpose of automatic wall treatments is to make results insensitive with respect to wall mesh refinement. Many blending approaches have been proposed. The one by Menter takes advantage of the fact that the solution to equations is known for both viscous and log layer

where is the cell centroid distance from the wall. Using this a blending can take the following form:

Note that for low values the will dominate and therefore viscous value of will be reproduced. Conversely, for larger values of , will be dominant and logarithmic value will be recovered.

Subsequently Menter proposes also blending for friction velocity. Friction velocity for viscous and logarithmic region are:

And the blending suggested:

FLUENT

Both k- omega models (std and sst) are available as low-Reynolds-number models as well as high-Reynolds-number models.

The wall boundary conditions for the k equation in the k- omega models are treated in the same way as the k equation is treated when enhanced wall treatments are used with the k- epsilon models.

This means that all boundary conditions for
- wall-function meshes will correspond to the wall function approach, while for the
- fine meshes, the appropriate low-Reynolds-number boundary conditions will be applied.

In Fluent, that means:

If the Transitional Flows option is enabled in the Viscous Model panel, low-Reynolds-number variants will be used, and, in that case, mesh guidelines should be the same as for the enhanced wall treatment
(y+ at the wall-adjacent cell should be on the order of y+ = 1. However, a higher y+ is acceptable as long as it is well inside the viscous sublayer (y+ < 4 to 5).)

If Transitional Flows option is not active, then the mesh guidelines should be the same as for the wall functions.
(For [...] wall functions, each wall-adjacent cell's centroid should be located within the log-law layer, 30 < y+ < 300. A y+ value close to the lower bound y+ = 30 is most desirable.)

References

• Menter, F., Esch, T. (2001), "Elements of industrial heat transfer predictions", 'COBEM 2001, 16th Brazilian Congress of Mechanical Engineering.'.
• ANSYS (2006), "FLUENT Documentation", .