National University of Singapore · FACULTY OF ENGINEERING

BPS5234A Chap.2 Turbulence Models and Near-Wall Treatment

- one subject, every graph, every model, every mark
6 Chapters4-page Bible
Our own words - no uploaded lecturer files
Updated for this semester
Chapter 2 of 12 · BPS5234A

Turbulence Models and Near-Wall Treatment

Why a model is needed at all

Turbulent airflow in a room or a street contains eddies from the size of the space down to millimetres. Resolving all of them directly is far beyond any practical computer for building-scale problems, so the flow is split into a mean and a fluctuating part and the equations are averaged.

Averaging leaves one new set of unknowns, the Reynolds stresses, and every turbulence model is a way of supplying them.

The families and what they cost

Eddy viscosity models assume turbulent stress behaves like viscous stress with a much larger, flow-dependent viscosity.

The k–&epsilon family computes that viscosity from turbulent kinetic energy and its dissipation and is the workhorse for design and compliance. The shear stress transport model blends a near-wall formulation with k–&epsilon in the free stream and copes better with separating flow and with heat transfer at walls.

Reynolds stress models drop the eddy viscosity idea and solve for each stress, which helps in swirling or strongly curved flow.

Large eddy simulation resolves the large eddies in time and models only the small ones, at a cost one to three orders of magnitude higher.

The evidence the course offers

Benchmark comparisons show where the cheap model fails: a ceiling jet predicted too thick, stagnation-point heat transfer over-predicted in an impinging jet, and vortex shedding averaged away.

They also show that for many mean-flow design questions the cheaper model is adequate.

The wall region

Near a surface the velocity changes steeply. Wall functions bridge that layer if the first cell sits in the right zone, and outdoor studies must also make the ground roughness consistent with the incoming wind profile.

In this chapter

What this chapter covers

  • 01

    Reynolds decomposition and the closure problem

  • 02

    Eddy viscosity and the k-epsilon variants

  • 03

    SST, Reynolds stress models and large eddy simulation

  • 04

    Benchmark evidence on jets and stratified rooms

  • 05

    Wall functions and the y-plus zones

  • 06

    Sand-grain and aerodynamic roughness in wind studies

Worked example · free

How large is the turbulent viscosity?

Q [5 marks]. In a ventilated room the k–&epsilon model returns k = 0.5 m²/s² and ε = 0.1 m²/s³ at a point. With Cμ = 0.09, air density 1.2 kg/m³ and molecular viscosity 1.8 × 10−5 Pa s, find the turbulent viscosity and compare it with the molecular value. The step weighting shown here is an AskSia study allocation and is not a University marking scheme.
Figure 2.1.Figure 2.1. A velocity record at one point: RANS solves for the dashed mean and models the effect of the fluctuations about it. mean value instantaneous velocity fluctuation u′ about the mean time u at one point
Figure 2.1. A velocity record at one point: RANS solves for the dashed mean and models the effect of the fluctuations about it.
  • 1Write the eddy viscosity relation of the k-epsilon model.
  • 2Substitute and evaluate the turbulent viscosity.
  • 2Form the ratio to molecular viscosity and interpret it.
The model sets turbulent viscosity as Cμ ρ k² / ε. Substituting, 0.09 × 1.2 × 0.5² / 0.1 = 0.09 × 1.2 × 0.25 / 0.1 = 0.27 Pa s. Dividing by the molecular viscosity gives 0.27 / 1.8 × 10−5 = 15,000. Turbulent mixing therefore transports momentum about fifteen thousand times more effectively than molecular viscosity at this point, which is why the closure, and not the air’s own viscosity, controls the spreading of jets and plumes in a room. It also shows why the turbulent value must be computed everywhere: it is a property of the flow, not of the fluid.
Sia tip — Square k before multiplying; leaving it unsquared under-predicts the turbulent viscosity by a factor equal to k.
Glossary

Key terms

Reynolds Decomposition
Splitting an instantaneous flow variable into a mean value and a fluctuation about it before averaging the equations.
Eddy Viscosity
A flow-dependent viscosity used to model turbulent stress in proportion to the mean rate of strain.
Turbulent Kinetic Energy
The kinetic energy per unit mass contained in velocity fluctuations, written k.
SST Model
A two-equation closure that uses a k-omega form near walls and a k-epsilon form away from them.
Large Eddy Simulation
A transient approach that resolves the large turbulent structures and models only the sub-grid scales.
Aerodynamic Roughness Length
The height parameter describing terrain roughness in a logarithmic wind profile.
Log-Law Region
The part of the near-wall layer where velocity varies with the logarithm of wall distance, where standard wall functions apply.
FAQ

Turbulence Models and Near-Wall Treatment FAQ

Which turbulence model should a building study use by default?

For mean-flow design and compliance work, a realizable or RNG k-epsilon model is the usual starting point and is what most of the course labs use. Move to SST when separating and reattaching flow or wall heat transfer decide the answer, and reserve Reynolds stress models and LES for anisotropic flows, transient phenomena or research.

Why does k-epsilon over-predict heat transfer at the stagnation point of an impinging jet?

Heat transfer there is governed by a thin viscous sub-layer, and the standard model assumes fully turbulent flow, so it produces too much turbulence near the wall. The lecture shows the effect against measured Nusselt numbers for a circular jet.

What does y-plus tell you about a mesh?

It is the non-dimensional distance of the first cell centre from the wall. Values between about 40 and 130 place the cell in the log-law region where standard wall functions work; enhanced wall treatment and low-Reynolds models need much smaller values.

Why can an inflow wind profile change before it reaches the buildings?

If the wall function on the ground uses a sand-grain roughness that does not correspond to the aerodynamic roughness implied by the inlet profile, the flow adjusts as it travels across the empty upstream region, so the buildings receive a different profile from the one specified.

Is LES always more accurate than RANS for a building?

It resolves unsteady structures and peaks that RANS averages away, and benchmarks rate it highly for shedding and stratified flow, but it costs far more and needs finer meshes and time steps. For comparing design options on mean quantities, a well-verified RANS study is often the better use of time.

Study strategy

Assessment move

Build a one-page table of your own with four rows, one per model family, and three columns: what it solves, a flow it handles well, and a flow it gets wrong with the reason. Fill the reasons from the benchmark examples rather than from memory, then test yourself on the reasons only.

A+Everything unlocked
Unlocks this Bible + all 25 of your National University of Singapore subjects - and 1,000+ Bibles across every Australian university.
Sia - your BPS5234A tutor, unlimited, worked the way the exam marks it
The full 4-page Bible + practice bank with worked solutions
Chrome extension - sync your LMS so Sia knows your deadlines
Bilingual EN / Chinese on every Bible and every Sia answer
$0.99 Trial
30-day money-back · cancel in one tap · how it works