Method
When using any CFD method, we first need to make a grid over the domain. The number of grid nodes in each direction depends on the type of flow, and what accuracy or reliability the user wants for the predictions. The gov erning equations have been discretized on this grid, and the finer the grid, the more accurate the results. If the user only wants an overall picture of the flow field, a 30 x 30 x 30 grid could be sufficient. The grid should be refined in regions where large gradients are expected. This normally means that the grid should be made finer near walls and in the region of the inlet.
It should also be remembered that the discretization scheme influences the accuracy of the results. In most CFD codes, different discretization schemes can be chosen for the convective terms. Usually, one can choose between firstorder schemes (e. g., the firstorder upwind scheme or the hybrid scheme) or secondorder schemes (e. g., a secondorder upwind scheme or some modified QUICK scheme). Secondorder schemes are, as the name implies, more accurate than firstorder schemes. However, it should also be remembered that the secondorder schemes are numerically more unstable than the firstorder schemes. Usually, it is a good idea to start the computations using a firstorder scheme. Then, when a converged solution has been obtained, the user can continue the calculations with a secondorder scheme.
When a coarse grid is used, w’all functions are used for imposing boundary conditions near the walls (Section 11.2.3.3). The nondimensional wall distance should be 30 < Y+ < 100, where Y+ = u*y/v. We cannot compute the friction velocity w. before doing the CFD simulation, because the friction velocity is dependent on the flow. However, we would like to have an estimation of Y* to be able to locate the first grid node near the wall at 30 < Y+ < 100. It we can estimate the maximum velocity in the boundary layer, the friction velocity can be estimated as 0.04f/max. After the computation has been carried out, we can verify that 30 < Y+ < 100 is satisfied for the grid nodes adjacent to the walls.
If the user wants higher accuracy in the predictions, and if he/she is willing to pay for the increase in accuracy in terms of increased computing cost, more grid nodes should be used. In general, more grid nodes should be located where the flow is complex. For an empty room without furniture or persons, this normally means that more grid nodes should be placed near walls and in




! FIGURE 11.4 Aspect ratio and biasing for control volumes. Aspect ratio is defined as Ax/Ay for cell /. Biasing in the X direction is defined as A^+j/Ax^,
The inlet region. If heat transfer at walls is to be predicted, it is usually very important to refine the grid near the wails. When using a fine grid near the walls, lowReynoldsnumber (LRN) models should be used (see Section 11.2.3.3). In this case the grid nodes should be located so that Y+ — 1 near the wall. As previously explained, we cannot compute the friction velocity U, in advance, but it can be estimated. After the computation has been carried out, the friction velocity can be computed by using the fact that the velocity in the viscous sublayer is linear (valid for y+ < 5), given by
A high aspect ratio (e. g., above 100) usually does not pose any problems (Fig. 11.4). Biasing causes larger problems, both in terms of convergence and accuracy. Preferably, biasing should be smaller than approximately 1.1, that is, the cell size between two adjacent cells in one direction should not increase or decrease by more than 10%.
I 1.2.2.2 Equations
The incompressible, timeaveraged continuity and the NavierStokes equations can be written as
(11.t) 
3^ = 0 OXj
,0U, DU V{ dxf + dxt 
(11.2) 
—! + — r(U U,) = +
Dt oXj 1 ‘ pdXj dXj
+ r7ret)
‘jlsLL —, Prdx, Ur 
(11.3) 
— + —(UT) = — Bt Bar, I ‘ dx:
Where
Pressed in °C then /3 = 1/(273 + Tref). The X2 coordinate increases vertically upwards. Turbulent stresses and turbulent heat fluxes Ujt appear on the righthand side of Eqs. (11.2) and (11.3). They act as additional diffusion terms due to correlations between fluctuating velocities and temperature, which are unknown. Transport equations can be derived for Utuj and N !. The higherorder unknown correlations in the resulting equations are then modeled; this type of model is called a Reynolds stress model.[3]
+ Tref) 
In ventilation problems, it is often sufficient to use simpler turbulence models, such as eddyviscosity models. and Ujt are then replaced by introducing a turbulent viscosity; Eqs. (11.2) and (11.3) can be rewritten as
(11.5)
(11.6)
Commonly used eddyviscosity turbulence models are the Ke model and the Ko) model. The eddy viscosities for these models have the form
(11.7) 
K2 , , , = *~e model
Il, = C..Q— Ka> model.
^^ <o
A derivation can be found in the literature.2^
If the flow is isothermal, there is no need to solve for the temperature equation (Eq. (11.6)). In this case the last term in Eq. (11.5) is also dropped. If, however, the thermal comfort is simulated, then the temperature equation must be solved. In ventilation the temperature variations are normally small, which means that it is sufficient to account for density variation only in the gravitation term (the last term in Eq. (11.5)). The gravitation term acts in the vertical direction, and in Eq. (11.5) it is assumed that the X2 coordinate is directed vertically upward. Tref denotes a reference temperature, which should be constant. It does not influence the predicted results, except that the pressure level is changed. It could, however, affect convergence rate (i. e., increase the number of required iterations required to reach a converged solution), and it should be chosen to a reasonable value, such as the inlet temperature.
The temperature equation is derived from the energy’ equation, in which the units of each term is joules per unit volume per second, J/(m3 s) = W/m5. The temperature equation above has been divided by the specific heat Cf) and density P (assumed to be constant), and thus the units of each term in Eq.
(11.6) is °C/s. If a heat source Q is to be added in a cell, it should be divided by
Posted in INDUSTRIAL VENTILATION DESIGN GUIDEBOOK