Inlet Turbulence Parameters

CFD simulations require specification of turbulence variables at inlet boundaries.

When measured turbulence data is available you can explicitly specify turbulence variables. For example, eddy viscosity for the Spalart-Allmaras (SA) model, turbulent kinetic energy and eddy frequency/dissipation rate for k-ε and k-ω based models. When measured data is not available, there are estimations for the turbulence values that are based on the turbulence intensity and turbulence length scale or eddy viscosity ratio.

Turbulence Intensity

The turbulence intensity (I) is defined as the ratio of the root-mean-square of the turbulent velocity fluctuations ( u ' ) and the mean velocity ( u ¯ ): (1)
I u ' u ¯

where u = ( u x ' 2 ¯ + u y ' 2 ¯ + u z ' 2 ¯ ) 3 = 2 k 3 and u ¯ = ( u ¯ x 2 + u ¯ y 2 + u ¯ z 2 ) .

For fully developed internal flows the turbulence intensity can be estimated as (2)
I = 0.16   R e h 1 / 8

where R e h = ρ u   ¯ D h μ is the Reynolds number and D h is a hydraulic diameter.

The following turbulence intensity values can be assumed
  • 5 percent < I < 20 percent for rotating machineries, for example, turbines and compressors.
  • 1 percent < I < 5 percent for internal flows.
  • I ~0.05 percent for external flows.

Turbulence Length Scale

The turbulence length scale represents the characteristic size of the turbulent eddies within a flow field. This parameter is often used to characterize the nature of the turbulence and appears in some form in nearly every turbulence model. RANS turbulence models have different definitions of the turbulence length scale based on the turbulence model and the type of application.

Estimated turbulence length scales based on the flow conditions:
  • l = 0.038 D h for fully developed internal flows.
  • l = 0.22 δ for developing flows.

where δ 0.382 x R e x 1 / 5 is the turbulent boundary layer thickness over a flat plate and R e x = ρ u ¯ x μ is the Reynolds number and x is the distance from the start of the boundary layer.

Eddy Viscosity Ratio

The eddy viscosity ratio (   ν r ) is the ratio between the eddy viscosity ( ν t ) and fluid kinematic viscosity ( ν ). The eddy viscosity ratio is set between 1 and 10 for internal flows, while it should be between 0.2 and 1.3 for external flows.

Inlet Turbulence Specification for the SA Model

The SA model has three options to estimate the turbulence variable at the inlet.
  • Option 1: Eddy viscosity ν t
  • Option 2: Turbulence intensity I and length scale l . For the eddy viscosity, ν t = 3 2 u ¯   I   l .
  • Option 3: Viscosity ratio   ν r . For the eddy viscosity, ν t = ν   ν r

Inlet Turbulence Specification for the k-ε Models

  • Option 1: Kinetic energy k and dissipation rate ε
  • Option 2: Turbulence intensity I and length scale l . For the kinetic energy, k = 3 2 ( u ¯ I ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8 qacaWGRbGaeyypa0ZaaSaaaeaacaaIZaaabaGaaGOmaaaacaGGOaGa bmyDa8aagaqea8qacaWGjbGaaiyka8aadaahaaWcbeqaa8qacaaIYa aaaaaa@3DF5@ . For the dissipation rate, ε = C μ k 3 / 2 l , where C μ = 0.09.
  • Option 3: Turbulence intensity I and viscosity ratio   ν r . For the kinetic energy, k = 3 2 ( u ¯ I ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8 qacaWGRbGaeyypa0ZaaSaaaeaacaaIZaaabaGaaGOmaaaacaGGOaGa bmyDa8aagaqea8qacaWGjbGaaiyka8aadaahaaWcbeqaa8qacaaIYa aaaaaa@3DF5@ . For the dissipation rate, ε = C μ   k 2 ν   ν r .

Inlet Turbulence Specification for the k-ω Models

  • Option 1: Kinetic energy k and eddy frequency ω
  • Option 2: Turbulence intensity I and length scale l . For the kinetic energy, k = 3 2 ( u ¯ I ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8 qacaWGRbGaeyypa0ZaaSaaaeaacaaIZaaabaGaaGOmaaaacaGGOaGa bmyDa8aagaqea8qacaWGjbGaaiyka8aadaahaaWcbeqaa8qacaaIYa aaaaaa@3DF5@ . For the eddy frequency, ω = k l     C μ 1 / 4 .
  • Option 3: Turbulence intensity I and viscosity ratio   ν r . For the kinetic energy, k = 3 2 ( u ¯ I ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8 qacaWGRbGaeyypa0ZaaSaaaeaacaaIZaaabaGaaGOmaaaacaGGOaGa bmyDa8aagaqea8qacaWGjbGaaiyka8aadaahaaWcbeqaa8qacaaIYa aaaaaa@3DF5@ . For the dissipation rate, ω = k ν   ν r .