Turbulence Intensity Calculator

Definition and Physical Meaning

Turbulence intensity is defined as:

I = \fracu'_rmsU_mean = \frac√(\frac1)3(u'^2 + v'^2 + w'^2)U

For isotropic turbulence: I = \frac√(\frac2k)3U

Where: • u'_rms = root-mean-square of velocity fluctuations (m/s) • U_mean = mean velocity (m/s) • k = turbulent kinetic energy (m²/s²) • u', v', w' = velocity fluctuation components

**Physical interpretation:** • Low intensity (< 1%): Very low turbulence (wind tunnels, smooth flows) • Medium intensity (1-5%): Typical external flows • High intensity (5-20%): Internal flows, downstream of obstacles • Very high (> 20%): Behind grids, in mixers, combustors

I = \fracu'_rmsU = (√(2k/3))/(U)

Turbulent Kinetic Energy

Turbulent kinetic energy per unit mass:

k = (1)/(2)\overlineu_i' u_i' = (1)/(2)(\overlineu'^2 + \overlinev'^2 + \overlinew'^2)

For isotropic turbulence (\overlineu'^2 = \overlinev'^2 = \overlinew'^2): k = (3)/(2)\overlineu'^2

**Relationship to turbulence intensity:** k = (3)/(2)(U · I)^2

**Typical values:** • Free stream, far from disturbances: k ≈ 0.001 m²/s² • Boundary layer, attached flow: k ≈ 0.1 m²/s² • Separated regions, wakes: k ≈ 1-10 m²/s² • High-intensity mixers: k ≈ 10-100 m²/s²

k = (3)/(2)(U · I)^2

Turbulent Dissipation Rate

Dissipation rate (ε) represents energy transfer from large to small scales:

\varepsilon = \frack^3/2\ell

Where \ell is turbulence length scale.

**Length scale estimation:** • External flow: \ell ≈ 0.07 L (7% of characteristic dimension) • Internal flow (pipes): \ell ≈ 0.07 D_h (7% of hydraulic diameter) • Boundary layer: \ell ≈ 0.4 \delta (40% of BL thickness) • Mixing layer: \ell ≈ layer thickness

**Alternative formulation:** \varepsilon = C_\mu^3/4 \frack^3/2\ell

Where C_\mu = 0.09 (standard k-ε model constant).

**Practical estimate:** For turbulence intensity I and length scale \ell: \varepsilon = C_\mu^3/4 ((\frac3)/(2)U^2I^2)^3/2\ell

\varepsilon = \frack^3/2\ell

Specific Dissipation Rate (ω)

For k-ω models, specific dissipation rate:

\omega = (\varepsilon)/(C_\mu k) = \frack^1/2C_\mu^1/2 \ell

**Alternative formulation:** \omega = (√(k))/(\ell)

**Relationship to other variables:** \varepsilon = C_\mu k \omega

**Typical values:** • Free stream: \omega ≈ 1-10 s⁻¹ • Boundary layer: \omega ≈ 10-100 s⁻¹ • Near wall: \omega → ∞ (approaches \omega = 6\nu/(β y^2)) • High turbulence: \omega ≈ 100-1000 s⁻¹

**For inlet boundaries:** \omega = \frack^0.5C_\mu^0.25 \ell

\omega = (√(k))/(C_\mu^1/4) \ell

Practical Turbulence Intensity Values

**Low turbulence intensity (< 1%):** • High-quality wind tunnels • Free stream far from disturbances • Smooth external flows • Example: Aircraft cruise conditions ≈ 0.1%

**Medium turbulence intensity (1-5%): ** • Atmospheric boundary layer • Smooth pipes and ducts • Most external aerodynamic flows • Example: Car on highway ≈ 1-3%

**High turbulence intensity (5-20%):** • Flow downstream of screens/grids • Rough pipes and ducts • Flow around buildings • Example: Urban environment ≈ 10-15%

**Very high turbulence intensity (> 20%):** • Mixing tanks and agitators • Combustion chambers • Immediate wake behind bluff bodies • Example: Stirred reactor ≈ 30-50%

**Rule of thumb for pipes:** I ≈ 0.16 Re_D^-1/8

For Re_D = 50,000: I ≈ 4.4\%

I ≈ 0.16 Re_D^-1/8

Turbulence Length Scale Selection

**External flows:** • Boundary layer: \ell = 0.4\delta where \delta is BL thickness • Free stream: \ell = 0.07L where L is body length • Wake: \ell ≈ wake width

**Internal flows:** • Fully developed pipe: \ell = 0.07D_h • Entrance region: \ell = 0.07D_h to 0.4D_h • Channel: \ell = 0.07H where H is channel height

**Complex geometries:** • Use 7% rule: \ell = 0.07L_char • For multiple scales: use smallest relevant dimension • If uncertain: \ell = 0.07L_char is reasonable default

**Effect on results:** • Larger \ell: Slower turbulence decay • Smaller \ell: Faster decay • Impact diminishes as flow develops • Most significant near inlet (first 5-10 characteristic lengths)

\ell = 0.07 L_characteristic

CFD Boundary Condition Specification

**k-ε models (specify k and ε):** k = (3)/(2)(UI)^2 \varepsilon = C_\mu^0.75 \frack^1.5\ell

**k-ω models (specify k and ω):** k = (3)/(2)(UI)^2 \omega = \frack^0.5C_\mu^0.25\ell

**Direct specification:** • Turbulence intensity: I (as percentage) • Turbulence length scale: \ell (in meters) • Most solvers calculate k, ε, ω automatically

**Sensitivity:** • Results near inlet sensitive to inlet turbulence • Far from inlet: boundary layer or mixing determines local turbulence • Conservative approach: test range of values (e.g., I = 1%, 5%, 10%) • For critical regions: measure or use validated literature values

Worked Example

Inputs

• velocity: 10
• turbulenceIntensity: 5
• lengthScale: 0.1

Result: k = 0.225 m²/s², ε = 0.234 m²/s³, ω = 2.88 s⁻¹

Explanation

**Example 1: External Flow - Flow Over Vehicle**

Air flow over vehicle: • Velocity: U = 10 m/s • Turbulence intensity: I = 5\% (typical road conditions) • Characteristic length: L = 5 m (vehicle length) • Turbulence length scale: \ell = 0.07L = 0.35 m

**Step 1: Calculate turbulent kinetic energy** k = (3)/(2)(UI)^2 = (3)/(2)(10 × 0.05)^2 = (3)/(2)(0.5)^2 = (3)/(2) × 0.25 = 0.375\text m^2\text/s^2

**Step 2: Calculate dissipation rate** \varepsilon = C_\mu^0.75 \frack^1.5\ell = 0.09^0.75 \frac(0.375)^1.50.35 = 0.164 × (0.229)/(0.35) = 0.107\text m^2\text/s^3

**Step 3: Calculate specific dissipation (for k-ω)** \omega = (√(k))/(C_\mu^0.25)\ell = (√(0.375))/(0.09^0.25) × 0.35 = (0.612)/(0.183) = 3.34\text s^-1

**CFD inlet conditions:** • k-ε model: k = 0.375 m²/s², \varepsilon = 0.107 m²/s³ • k-ω model: k = 0.375 m²/s², \omega = 3.34 s⁻¹ • Turbulence intensity: I = 5\%, length scale: \ell = 0.35 m

---

**Example 2: Internal Flow - Pipe Flow**

Water flow in pipe: • Velocity: U = 2 m/s • Pipe diameter: D = 0.1 m • Reynolds number: Re_D = 200,000

**Step 1: Estimate turbulence intensity** I = 0.16 Re_D^-1/8 = 0.16 × (200,000)^-0.125 = 0.16 × 0.281 = 0.045 = 4.5\%

**Step 2: Determine length scale** \ell = 0.07D_h = 0.07 × 0.1 = 0.007\text m = 7\text mm

**Step 3: Calculate k and ε** k = (3)/(2)(UI)^2 = (3)/(2)(2 × 0.045)^2 = (3)/(2) × 0.0081 = 0.0122\text m^2\text/s^2

\varepsilon = 0.164 \frack^1.5\ell = 0.164 × \frac(0.0122)^1.50.007 = 0.164 × (0.00135)/(0.007) = 0.0316\text m^2\text/s^3

**Step 4: Calculate ω** \omega = (√(0.0122))/(0.09^0.25) × 0.007 = (0.1105)/(0.00365) = 30.3\text s^-1

---

**Example 3: Wind Tunnel Test Section**

Low-turbulence wind tunnel: • Velocity: U = 50 m/s • Turbulence intensity: I = 0.1\% (high-quality tunnel) • Test section size: 1 × 1 m • Length scale: \ell = 0.07 × 1 = 0.07 m

**Calculate turbulence parameters:** k = (3)/(2)(50 × 0.001)^2 = (3)/(2) × 0.0025 = 0.00375\text m^2\text/s^2

\varepsilon = 0.164 × \frac(0.00375)^1.50.07 = 0.164 × (0.000229)/(0.07) = 0.000536\text m^2\text/s^3

\omega = (√(0.00375))/(0.183 × 0.07) = (0.0612)/(0.0128) = 4.78\text s^-1

**Interpretation:** • Very low k value reflects minimal turbulence • Low \varepsilon means slow energy cascade • These conditions preserve laminar flow or delay transition • Critical for accurate drag measurements on models

Second Scenario

Inputs

• velocity: 7.5
• turbulenceIntensity: 5
• lengthScale: 0.1

Result: k = 0.225 m²/s², ε = 0.234 m²/s³, ω = 2.88 s⁻¹

Explanation

This scenario uses different inputs (velocity = 7.5, turbulenceIntensity = 5, lengthScale = 0.1) to show how changing one variable affects the turbulence intensity result. Run the calculator above with these values to get the exact updated output with step-by-step work.