Objective Review Reynolds Navier Stokes Equations RANS Learn

Objective • Review Reynolds Navier Stokes Equations (RANS) • Learn about General Transport equation • Start with Numerics

From the previous class Reynolds Averaged Navier Stokes equations (incompressible flow) Reynolds stresses total 9 - 6 are unknown same Total 4 equations and 4 + 6 = 10 unknowns We need to model the Reynolds stresses !

From the previous class Modeling of Reynolds stresses Eddy viscosity models Average velocity Boussinesq eddy-viscosity approximation Is proportional to deformation Coefficient of proportionality k = kinetic energy of turbulence Substitute into Reynolds Averaged equations

From the previous class Reynolds Averaged Navier Stokes equations Continuity: 1) Momentum: 2) 3) 4) Similar is for STy and STx 4 equations 5 unknowns → We need to model

Modeling of Turbulent Viscosity Fluid property – often called laminar viscosity Flow property – turbulent viscosity MVM: Mean velocity models TKEM: Turbulent kinetic energy equation models Additional models: LES: RSM: Large Eddy simulation models Reynolds stress models

Kinetic energy and dissipation of energy Kolmogorov scale Eddy breakup and decay to smaller length scales where dissipation appear

One equation models: Prandtl Mixing-Length Model (1926) Vx y x l Characteristic length (in practical applications: distance to the closest surface) -Two dimensional model -Mathematically simple -Computationally stable -Do not work for many flow types There are many modifications of Mixing-Length Model: - Indoor zero equation model: t = 0. 03874 V l Distance to the closest surface Air velocity

Two equation turbulent model Kinetic energy Energy dissipation From dimensional analysis constant We need to model Two additional equations: kinetic energy dissipation

Reynolds Averaged Navier Stokes equations Continuity: 1) Momentum: 2) 3) 4) General format:

General CFD Equation Values of , , eff and S , eff Continuity 1 0 0 x-momentum V 1 + t - P/ x+Sx y-momentum V 2 + t - P/ y- g (T∞-Twall)+Sy z-momentum V 3 + t - P/ z+Sz T-equation T / l + t/ t ST k-equation k ( + t)/ k G- +GB -equation ( + t)/ [ (C 1 G-C 2 )/k] +C 3 GB( /k) Species C ( + t)/ c SC Age of air t + t Equation t = C k 2/ , G= t ( Ui/ xj + Uj/ xi) Ui/ xj S , GB=-g( /CP)( t/ T, t) T/ xi C 1=1. 44, C 2=1. 92, C 3=1. 44, C =0. 09 , t=0. 9, k =1. 0, =1. 3, C=1. 0

Numerics

Finite Volume Method - Conservation of f for the finite volume Divide the whole computation domain into sub-domains h w One dimension: n e W P dx s E dx Dx l e w - Finite volume is a fixed space in the flow domain with imaginary boundaries that allow the fluid to flow in and out. - Integral conservation of the quantities such as mass, momentum and energy. f

General Transport Equation -3 D problem steady-state H N W P E S L Equation for node P in the algebraic format:

1 -D example of discretization of general transport equation Steady state 1 dimension (x): dxw P W w Dx dxe E e Point W and E represent the cell center of the west and east neighbors of cell P and w, e the neighboring surfaces. Integrating with Gaussian theorem on this control volume gives: To obtain the equations for the value at point P, assumptions are used to convert the surface values to the center values.

Convection term dxw P W dxe E Dx w – Central difference scheme: - Upwind-scheme: If Vx>0 and If Vx<0 and e

Diffusion term dxw P W dxe E Dx w e

Summary: Steady–state 1 D I) X direction Convection term - Upwind-scheme: a) If Vx > 0, and If Vx < 0, and W dxw P dxe E Dx w e Diffusion term: b) When mesh is uniform: DX = dxe = dxw c) Source term: Assumption: Source is constant over the control volume

1 D example - uniform mesh After substitution a), b) and c) into I): We started with partial differential equation: same and developed algebraic equation: We can write this equation in general format: Unknowns Equation coefficients

1 D example multiple (N) volumes N unknowns 1 2 3 i Equation for volume 1 N equations Equation for volume 2 ……………… Equation matrix: For 1 D problem 3 -diagonal matrix N-1 N

3 D problem Equation in the general format: H N W P E S L Wright this equation for each discretization volume of your discretization domain A F 60, 000 cells (nodes) N=60, 000 elements x = 7 -diagonal matrix This is the system for only one variable ( ) When we need to solve p, u, v, w, T, k, e, C system of equation is larger

Convection term dxw P W dxe E Dx w – Central difference scheme: - Upwind-scheme: and e