Entropy of entanglement

CFD stands for computational fluid dynamics (and heat transfer). As per this technique, the governing differential equations of a flow system or thermal system are known in the form of Navier–Stokes equations, thermal energy equation and species equation with an appropriate equation of state. In the past few years, CFD has been playing an increasingly important role in building design, following its continuing development for over a quarter of a century. The information provided by CFD can be used to analyse the impact of building exhausts to the environment, to predict smoke and fire risks in buildings, to quantify indoor environment quality, and to design natural ventilation systems, etc.

Applications

Recently CFD finds very wide application in different areas of science and engineering; some examples are:

• Aerodynamics of aircraft and vehicles : lift and drag
• Hydrodynamics of ships
• Power plant : combustion in internal combustion engines and gas turbines
• Turbo machinery: Flows inside rotating passages, diffusers etc.
• Electrical and electronics engineering: cooling of equipment including microcircuits.
• Chemical process engineering: mixing and separation and polymer moulding.
• Marine engineering: loads on off-shore structure.
• Environmental engineering: distribution of pollutant and effluents.
• Hydrology and oceanography: flows in rivers, estuaries and oceans.
• Meteorology: weather prediction.
• Biomedical engineering: blood flows through arteries and veins.
• External and internal environment of buildings: wind loading, ventilation analysis and heating/cooling load calculations.

In early age of construction, the most of building related issues such as ventilation analysis, wind loading, wind environment etc. were conducted by the wind tunnel tests, but today all these test can be done effectively with CFD technique. CFD technique can resolve all above mentioned issues in very short time period and it is very economical as well as strong approach than the older one (experimental):. Recently Computational fluid dynamics is used as a sophisticated airflow modelling method and can be used to predict airflow, heat transfer and contaminant transportation in and around buildings. CFD plays an important role in building design, designing a thermally-conformable, healthy and energy-efficient building. CFD can examine the effectiveness and efficiency of various Heating ventilation and air conditioning (HVAC) systems by easily changing the different types and location of different components of diffuser types and locations, supply air conditions and system control schedules. Furthermore, CFD helps in developing passive heating/cooling/ventilation strategies (e.g. natural ventilation) by modelling and optimizing building site-plans and indoor layouts. Globally building sector shares approximately 40% of total energy consumption . In present era, there is a huge gap in energy consumption and energy production. As building sector share a huge amount of the total consumption, hence it becomes essential to investigate the optimum configuration for building to reduce the building's share of energy. In order to achieve this, CFD can play an important role. Energy simulation and CFD programs are important building design tools which are used for evaluation of building performance, including thermal comfort, indoor air quality mechanical system efficiency and energy consumption . . CFD in buildings mainly used for one or more followings purposes:

1. Thermal analysis: through walls, roof and floor of buildings
2. Ventilation analysis.
3. Orientation, site and location selection of buildings based on local geographical and environmental conditions.

Thermal analysis

In buildings, heat transfer takes place in its all modes i.e. conduction, convection and radiation. In order to reduce heat losses from buildings, CFD analysis can be done for the optimum configuration of composite walls, roof and floor. The differential form of the general transport equation is as follows .

${\frac {\partial {(\rho \phi )}}{\partial t}}+{div\,(\rho u\phi )}={div\,(k\,grad\,\phi )}+{S_{\phi }}\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,(1)$

The numerical solution of above equation can be obtained by finite difference method (FDM), finite volume method (FVM) and finite element method (FEM). In buildings, for heat transfer analysis, the scalar function ф in equation (1) is replaced by Temperature (T), diffusion coefficient Γ is replaced by thermal conductivity k and the source term $S_{\phi }$ is replaced by heat generation term e or by any heat radiation source $Q_{i}$ or by both (depending upon the nature of source available) and we have different forms of equation for different cases. For simplicity and easy understanding, only 1-Dimensional cases have been discussed.
In buildings the heat transfer analysis can be done for all parts of buildings (walls, roof and floor) in following two ways

1. Steady State Thermal Analysis
2. Transient Thermal Analysis

Steady state thermal analysis

The steady state thermal analysis consist the following type of governing differential equations.
Case-1: General steady state heat conduction equation.
For this case the governing differential equation (GDE) (1) becomes as follows:

${div\,(\rho uT)}={div\,(k\,grad\,T)}+{S_{T}}\,$

Case-2: Steady state heat conduction equation (no heat generation)
For this case the governing differential equation (GDE) (1) becomes as follows:

${div\,(\rho uT)}={div\,(k\,grad\,T)}\,$

Case-3: Steady state heat conduction equation (no heat generation and no convection)
For this case the governing differential equation (GDE) (1) becomes as follows:

${div\,(k\,grad\,T)}=0\,$

Transient thermal analysis

The transient thermal analysis consist the following type of governing differential equations.
Case-1: Transient heat conduction
For this case the governing differential equation (GDE) (1) becomes as follows:

${\frac {\partial {(\rho T)}}{\partial t}}+{div\,(\rho uT)}={div\,(k\,grad\,T)}+{S_{T}}\,$

Case- 2: Transient heat conduction (no heat generation) For this case the governing differential equation (GDE) (1) becomes as follows:

${\frac {\partial {(\rho T)}}{\partial t}}+{div\,(\rho uT)}={div\,(k\,grad\,T)}$

Case-3: Transient heat conduction (no heat generation and no convection)
For this case the governing differential equation (GDE) (1) becomes as follows:

${\frac {\partial {(\rho T)}}{\partial t}}={div\,(k\,grad\,T)}\,$

We can solve these above mentioned governing differential equation (GDE) equations using CFD technique.

Ventilation analysis

The ventilation study in buildings is done to find the thermally comfortable environment with acceptable indoor air quality by regulating indoor air parameters (air temperature, relative humidity, air speed, and chemical species concentrations in the air). CFD finds an important role in regulating the indoor air parameters to predict the ventilation performance in buildings. The ventilation performance prediction provides the information regarding indoor air parameters in a room or a building even before the construction of buildings .
These air parameters are crucial for designing a comfortable indoor as well as outdoor environment. This is because the design of appropriate ventilation systems and the development of control strategies need detailed information regarding the following parameters;

• Airflow
• Contaminant dispersion
• Temperature distribution

The aforesaid information’s are also useful for an architect to design the building configuration. From the last three decade, the CFD technique is widely used with considerable success in building.
Recently ventilation and its related fields has becomes a great part of wind engineering. A ventilation study can be done using wind tunnel investigation (experimentally) or by CFD modeling (theoretically). Natural ventilation system is always preferred over the forced ventilation system, as it causes to save the burning of fuel, which is economical as well as nature friendly. In present era, due to development of a lot of CFD software and other building's energy simulation software, it becomes quite easy to assess the possibility of natural/forced ventilation system in a building. CFD analysis is quite useful than the experimental approach because here we can find other related relations among the variables in post-processing. The data obtained either experimental or numerically is useful in two ways .
1) better comfort of user.
2) It provides the data which is used as input to the heat balance calculation of the buildings .

Orientation, site and location selection

Figure-1 (a): Flow around a building (collection of air at height and delivering it at ground level)

Figure-1 (b): Flow around a building (center of the front face )

Earlier, the choice of dwelling location was made mindful of the need for water, so most of earlier development started in valley area. In present era, due to advancement in science and technology, it becomes easy to select a proper orientation, site and location of buildings based on local geographical and environmental conditions. In selection of building site and location, wind loading plays and important role. In case of two buildings at a location exists side by side having some gap, when volume of wind blows round the ends of building through the gap is, in first instant the sum of flow around each building separately, then its velocity must increase above that around the end of a single building at the expense of pressure loss. So, there will be a built of pressure, entering the gap, which will lead to higher wind loads on the sides of buildings. When wind blows over the face of a high rise building, a vortex is created by the downward flow on the front face (as shown in figure-1). The wind speed in the reverse direction near the ground level may have 140 percent of the reference wind speed. So, if any building exist in such region, then that may be subjected to damage (especially the roof of building may get severe damage). Such damage to buildings can be prohibited successfully, if the effects of wind loading are considered in the early stage of construction of a building. In early age of construction all these wind loading effects were determined by the wind tunnel test, but today all these test can be successfully through CFD analysis. The importance of providing pleasant environment to buildings is increasing and architect and wind engineer are often asked to look over the design (orientation, site, location and gaps between the surrounding buildings) in the formative stage of buildings & planning stage of construction. So by using CFD analysis, we can find the suitable information (local wind velocity, convective coefficients, and solar radiation intensity) for orientation, site and location selection of buildings.

CFD approach for heat transfer analysis in buildings

CFD technique can be used for the analysis of heat transfer in each part of a building.CFD technique finds the solution by following ways:

1. Discretization of the governing differential equation using numerical methods (Finite difference method has been discussed).
2. Solve the discretized version of equation with high performance computers.

Discretization of the governing differential equations for the steady state heat transfer analysis

Consider a building having a plane wall with thickness L, heat generation e and constant thermal conductivity k. The wall is subdivided into M equal regions of thickness $\Delta x$ = X/T in x-direction, and the divisions between the regions are selected as nodes as shown in figure-2.

Figure-2:the nodal points and volume elements for the finite difference formulation of 1-D conduction in a plane wall

The whole domain of wall in x-direction is divided in elements as shown in figure and the size of all interior elements is same while for exterior elements it is half.

Now to obtain the FDM solution for the interior nodes, consider the element represented by the node m which is surrounded by neighboring nodes m-1 and m+1. The FDM technique presumes that temperature varies linearly in walls (shown in figure-3). FDM solution is:( for all interior nodes except to 0 and last node)

${\frac {(T_{m-1}^{i}-2T_{m}^{n}+T_{m}^{i})}{\Delta {x}^{2}}}+{\frac {e}{k}}=0$

Figure-3:Linear temperature variation in finite difference formulation

Boundary conditions

Above equation is valid only to interior nodes only. To obtain the solution for exterior nodes we have to apply the boundary conditions (as applicable), which are as follows. 1.Specified heat flux boundary condition

$q_{0}A+kA{\frac {(T_{1}-T_{0})}{\Delta {x}}}+{\frac {e_{0}}{2}}A\Delta {x}=0$

When boundary is insulated (q=0)

$kA{\frac {(T_{1}-T_{0})}{\Delta {x}}}+{\frac {e_{0}}{2}}A\Delta {x}=0$

2. Convective boundary condition

$hA{(T_{\infty }-T_{0})}+kA{\frac {(T_{1}-T_{0})}{\Delta {x}}}+{\frac {e_{0}}{2}}A\Delta {x}=0$

3. Radiation boundary condition

$\epsilon \sigma A{(T_{sur}^{4}-T_{0}^{4})}+kA{\frac {(T_{1}-T_{0})}{\Delta {x}}}+{\frac {e_{0}}{2}}A\Delta {x}=0$

4. Combined convective and radiation boundary condition (shown in figure-4).

$hA{(T_{\infty }-T_{0})}+\epsilon \sigma A{(T_{sur}^{4}-T_{0}^{4})}+kA{\frac {(T_{1}-T_{0})}{\Delta {x}}}+{\frac {e_{0}}{2}}A\Delta {x}=0$

or when radiation and convection heat transfer coefficient are combined, above equation becomes as follows;

$hA_{combined}{(T_{\infty }-T_{0})}+kA{\frac {(T_{1}-T_{0})}{\Delta {x}}}+{\frac {e_{0}}{2}}A\Delta {x}=0\,$

Figure-4: Schematic for the FDM formulation of combined convective and radiation on the left boundary of a plane wall

5. Combined convective, radiation and heat flux boundary condition

$q_{0}A+hA{(T_{\infty }-T_{0})}+\epsilon \sigma A{(T_{sur}^{4}-T_{0}^{4})}+kA{\frac {(T_{1}-T_{0})}{\Delta {x}}}+{\frac {e_{0}}{2}}A\Delta {x}=0$

6.Interface boundary condition : when there is an interface (in composite walls) of different walls having different thermo-physical properties, the two different solid media A and B are assumed to be perfect contact and thus have same temperature at interface at node m (as shown in figure-5).

$k_{A}A{\frac {(T_{m-1}-T_{m})}{\Delta {x}}}+k_{B}A{\frac {(T_{m+1}-T_{m})}{\Delta {x}}}+{\frac {e_{A,m}}{2}}A\Delta {x}+{\frac {e_{B,m}}{2}}A\Delta {x}=0$

Figure-5: Schematic for the FDM of the interface boundary condition for two mediums A and B having perfect thermal contact

In above equations q_0 = denotes specified heat flux is in $(W/m^{2})$, h =convective coefficient, $h_{combined}$ = combined convective and radiation heat transfer coefficient,$T_{s}ur$ = Temperature of surrounding surface,$T_{(}\infty )$ =Ambient Temperature, $T_{0}$ = Temperature of at initial node. Note: For interior side of wall we can apply the suitable boundary condition from above (as applicable), in that case $T_{(}\infty )$ will be replaced by $T_{r}$ (Room Temperature), $T_{0}$= will be replaced by $T_{l}$ (Temperature of last node).

Discretization of the governing differential equations for the transient heat transfer analysis

Transient thermal analysis is more important than the steady thermal analysis, as this analysis include the variable ambient condition with time. In transient heat conduction, the temperature changes with time as well as position. The finite difference solution of transient heat conduction requires discretization in time in addition to space, as shown in figure-6.

Figure-6: FDM fotirmulation of time dependent problem involves discrete points in time as well as in space

The nodal points and volume elements for the transient FDM formulation of 1-D conduction in a plane wall exist as shown in the figure-7.

Figure-7:The nodal points and volume elements for the transient FDM formulation of 1-D conduction in a plane wall

For this case the FDM explicit solution for equation (1) will be as follows,

$kA{\frac {(T_{m-1}^{i}-T_{m}^{i})}{\Delta {x}}}+kA{\frac {(T_{m+1}^{i}-T_{m}^{i})}{\Delta {x}}}+{e_{m}}A\Delta {x}=(\rho c_{p}\Delta xA){\frac {(T_{m}^{i+1}-T_{m}^{i})}{\Delta x}}$

The above equation can be solved explicitly for the temperature $(T_{m}^{i+1})$ to give

${T_{m}^{i+1}}=\tau {(T_{m+1}^{i}-T_{m}^{i})}+{(1-2\tau )}T_{m}^{i}+\tau {\frac {(e_{m}\Delta {x}^{2})}{k}}$

where,

$\tau ={\frac {(\alpha \Delta t)}{\Delta x^{2}}}\,$

and

$\alpha ={\frac {k}{\rho c_{p}}}\,$

here,$\tau$ represents the cell Fourier no,$\alpha$ represents thermal diffusivity, $c_{p}$ represents specific heat at constant pressure, $\Delta t$ represents time step,$\Delta x$ represents space step.

Above equation is valid for all interior nodes and to find the relation for first and last node, apply boundary conditions (as applicable) as discussed in steady state heat transfer. For a convective & radiation boundary if solar radiation data $q_{solar}$\, in ($(W/m^{2})$) is available and absorptivity-transmissivity constant K is known, the relation for temperature is obtained as follows;

$hA{(T_{\infty }^{i}-T_{0}^{i})}+\kappa Aq_{sol}=(\rho c_{p}\Delta xA){\frac {(T_{1}^{i}-T_{0}^{i})}{\Delta x}}$

Note: the thermal analysis for the roof and floor of a building can be done in same way, as discussed for walls.\\

See also

• Computational fluid dynamics
• CFDTemplate:Dn
• Natural ventilation
• JPMorgan Chase Tower (Houston)
• Environmental Systems Design, Inc.
• Dynamic insulation
• Thermal management of high-power LEDs
• Vented balance safety enclosure
• Different Types of Boundary Conditions in Fluid Dynamics
• Wind tunnel
• Greenhouse

External links

• [1]-Analysis of Building Performance using Computational Fluid Dynamics
• [2] - Images for cfd analysis in buildings
• [3]-DesignBuilder CFD
• [4]-CFD ANALYSIS OF FLOW AND CONCENTRATION FIELDS AROUND A BUILDING WITH A ROOF STACK
• [5]-CFD AS A BUILDING SERVICES ENGINEERING TOOL
• [6]-Performance of Coupled Building Energy and CFD Simulations
• [7]-APPLICATION OF CFD IN BUILDING PERFORMANCE SIMULATION FOR THE OUTDOOR ENVIRONMENT
• [8]-Integrating CFD and building simulation
• [9]-Building Design and Analysis
• [10]-Building Simulation
• [11]-Capability in Computational Fluid Dynamics (CFD) for Building Services
• [12]- CFD in HVAC and green building
• [13]-CFD Prediction of Airflow in Buildings for Natural Ventilation
• [14]-CFD Advantages and Practical Applications
• [15]- Numerical coupling of Modelica and CFD for building energy supply systems

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.