Open access peer-reviewed chapter

On Risk and Reliability Studies of Climate-Related Building Performance

By Krystyna Pietrzyk and Ireneusz Czmoch

Submitted: June 14th 2017Reviewed: October 13th 2017Published: February 28th 2018

DOI: 10.5772/intechopen.71684

Downloaded: 234

Abstract

A design strategy based on integration of the building form and structure with its external environment in order to take advantage of natural forces (wind and buoyancy effects) has been evaluated in terms of risk and reliability measures. Tools for the probabilistic analysis (First-Order Reliability Method (FORM), Monte Carlo) have been presented and applied in the probabilistic modelling and sensitivity analysis of the response function of the studied building physics problem. Sensitivity analysis of the influence of basic random variables on the probability distribution of a response function is straightforward in FORM methodology. The case-based studies of probabilistic modelling of uncertainties coupled to wind speed and temperature difference through the specified building/environment system have been presented (i.e., the distribution models of the air change rate ACH and the dynamic U value characterising thermal performance of dynamic insulation). Sensitivities of the probability model of ACH to the parameters of wind speed and temperature distributions have been estimated for the consecutive values of the air change rate using FORM methodology. Reliability of ACH turned out to be most sensitive to the shape parameter of the wind speed distribution (in two-parameter Weibull model). The probabilistic risk analysis along with the effective tools for sensitivity analysis can be used to support design decisions and also to develop better models for evaluation of building performance.

Keywords

  • building performance
  • environment
  • risk
  • reliability
  • probabilistic approximation
  • FORM
  • sensitivity
  • climate
  • climate mitigation
  • wind
  • air infiltration
  • ACH
  • dynamic U value

1. Introduction

Developing tools to support decision-making to ensure comfort and safety in built environment while taking into account climate change challenges becomes important. ‘Energy is the dominant contributor to climate change, accounting for around 60% of total global greenhouse gas emissions’ [1]. Reducing the carbon intensity of energy is a key objective in long-term climate goals. Hence, choosing a strategy based on integrating the building form and structure with its external environment in order to take advantage of natural forces (for natural ventilation, solar heating, etc.) is an example of design decisions leading towards mitigation of climate change.

Decisions concerning the choice of the design solutions for the particular project have to be taken under uncertainties related to the unknown and variable conditions, i.e. random climatic conditions, uncertain material performance, uncertain user behaviour, etc. Confronted with significant uncertainty, deterministic modelling supporting design process has been proved to be insufficient for decision-making. However, as it is said in [2] ‘Existing engineering-based models are unable to propagate uncertainties through the model, and are therefore limited in their ability to display the impact of uncertainties to decision makers’. Facing that challenge, the chapter includes the discussion of the models and the tools applied by the authors for the probabilistic transformation of uncertainties of climatic parameters through a building/environment system for the predictive modelling of building performance.

The method for the quantification of building performance in terms of probability of poor performance (failure) and satisfactory performance (safe behaviour, in general meaning) is presented. Next, the tools for the probabilistic analysis are described (FORM, Monte Carlo) in relation to probabilistic modelling and possible applications of sensitivity analysis. One of the important results of analysis is the probability distribution functions of different performances as the responses of building/environment systems to the environmental loads. Such analysis requires estimation of some climatic parameters in terms of frequency of occurrence and appropriate statistics.

The chapter includes the case-based studies of probabilistic transformation of uncertainties coupled to wind and temperature through the specified building/environment system to show the effect on the distribution model of the air change rate and further on the distribution model of the dynamic thermal transmittance (dynamic Uvalue) of the building envelope. Furthermore, the estimated distribution models could be included in risk/reliability calculations, carried out with FORM tools. The analysis of the sensitivity of the distribution of ACHwith respect to the randomness of wind speed and outdoor temperature exemplifies the potential of the FORM tools, which can be effectively used to find out the probabilistic characteristics typical for the combination of the important variables influencing climate-structure interaction.

2. Risk perspective on design for sustainable development

Design for sustainable development can be approached using risk analysis tools. To minimalise the risk of undesired consequences while increasing the chance to enhance the quality of life becomes the basic design objective. The design goal can be expressed in other terms—how to secure reliability of design under risk constraints [3]. To be clear about the terminology used further, some definitions are given below.

  • Risk—a state of uncertainty where some possible outcomes have an undesired effect or significant loss [4]. It can be expressed in terms of adverse consequences scaled by the probabilities of undesired outcomes.

  • Qualitative risk—Relative measure of risk based on ranking or separation into descriptive categories such as low, medium and high or on a scale, e.g., from 1 to 10. The example is the risk matrix of the failure mode [5] where failure severity is described by categories: minor, major, critical and catastrophic, whereas failure frequency is described in terms of very unlikely, remote, occasional, probably and frequent. A combination of both gives a qualitative assessment of risk as low, moderate or high.

  • Quantitative risk—The most common quantification of risk is the product of likelihood of occurrence and the effect of the hazardous event [6]. Risk is treated as convolution of hazard and vulnerability, and it refers to mutual conditioning of two phenomena [7].

  • Risk assessment is the determination of quantitative or qualitative estimate of risk related to the well-defined situation and the recognised hazard.

  • Probabilistic risk assessment denotes the methods dealing with computation both the likelihood of undesired event and severity of possible negative consequences due to occurrence of such event.

  • Risk management process—‘the systematic process of identification, assessment, monitoring and control of risk’ [8].

  • Reliability—Ability of a system to satisfactorily perform under the specified conditions of use over an intended period of time. It is quantified by the probability of satisfactory (or safe) performance [9]. Unreliability is measured by the probability of failure (undesired performance).

2.1. Risk perspective on climate change challenges

Climate change threatens life on our planet. In view of high uncertainty, qualitative or semi-quantitative risk analysis based on the different scenarios is often applied. Following the quantitative definition of risk, one can write

Risk=PhazardConsequencesE1

Phazardis the probability of occurrence of undesired events leading to possible Consequenceslike loss, injury, or discomfort.

Risk reduction could be accomplished by decreasing the probability of undesired event as well as diminishing the scale adverse consequences. Risk reduction of climate change and its consequences can be accomplished by climate change mitigation (decrease of the probability of occurrence of adverse events) or climate change adaptation (decrease of the adverse consequences) described as follows:

Climate change mitigation—‘it consists of actions to limit the magnitude and/or rate of long-term climate change’ [10]. ‘It generally involves reductions in human (anthropogenic) emissions of greenhouse gases´ [11].

Climate change adaptation—‘anticipating the adverse effects of climate change and taking appropriate action to prevent or minimize the damage they can cause, or taking advantage of opportunities that may arise’ [12].

2.2. Risk assessment as a tool supporting design of buildings

Designing for the integration of the building form and structure with its external environment in order to use natural forces to secure comfort (passive strategies) is an example of activities towards mitigation of climate change. If it is supported by probabilistic prediction of a local climate changes, it can be viewed from the climate adaptation perspective too.

The needs for risk reduction related to the hazards induced by climate change become an important boundary condition in the modelling of building/environment system to support building design. Following the definition of risk (Eq. (1)), adverse consequences are indicated by the set {yes = 1, no = 0}, and as a result, the probability Phazardbecomes the discriminating factor for comparison of different design solutions. It means that certain design could be chosen on the basis of comparison of probability of unsatisfactory performance or reliability, evaluated by a set of alternative design proposals. In that way risk assessment becomes a tool supporting design of buildings.

3. Transformation of uncertainties: probabilistic approach

Probability is a measure of uncertainty about future events. Probability of a performance of a building/environment system depends on the theoretical model used and the randomness of the influencing parameters. The epistemic uncertainty about the theoretical models applied together with the aleatory uncertainty coupled to the randomness of important phenomena contributes to the final uncertain outcome, as the result of transformation of uncertainties throughout the model. Methods and tools for probabilistic reliability analysis can be used to estimate the probabilistic response of the structure to the random climatic load. They could be an important part of risk-based design process built upon the framework of risk management methodology as proposed in [13].

3.1. Probabilistic analysis with FORM

Development of reliability methods resulted in variety of powerful algorithms to estimate probability of failure for complicated physical and mathematical models of building systems incorporating random variables (i.e. properties or actions). FORM denotes ‘First-Order Reliability Method’, which has been developed by many researchers in about 40 years ago. Short description of the FORM basics as well as sensitivity tools is presented below. For details, check [14], and for application in building physics, look in [15]. First-order reliability method (FORM) is the most popular approach applied in practice.

Once the response of a system characterised by a set of basic random variables and a mathematical model describing the relationship among them has been established, the probability density function of the response can be estimated with help of FORM tools.

In general, the performance of a system is analysed in the space ΩX=XRnof basic random variables X. For a given failure mode or serviceability requirement, represented by the limit state surface gX=0, the space ΩXis divided into the safe subset, i.e. satisfactory performance, ΩS=XRngX>0, and the failure subset, ΩF=XRngX0, i.e., unsatisfactory performance. If all random variables are continuous with the multivariate joint probability density function (PDF) fXx, the failure probability is given by the integral

Pf=ΩFfxxdxE2

The integral (Eq. (2)) can be evaluated exactly for a few cases with the most important one: the linear limit state surface and multidimensional normal (Gaussian) distribution function of variables X.

FORM algorithm starts with the non-linear transformation. Non-normal random vector Xis transformed into a standard normal (Gaussian) vector Ywith zero mean and unit covariance matrix CYY=I. The limit state surface gx=0is mapped into a limit state surface Gy=0. Next, the design point y, i.e. the point on the limit state surface with the minimum distance to the origin of the Yspace, is determined by solving the non-linear optimisation problem with a non-linear constrain Gy=0:

β=minyTyforyonGy=0E3

The hyperplane tangential to the limit state surface at the point yis given by formula

βαTy=0E4

where αis a unit outward normal vector to the hyperplane and βis the distance between the hyperplane and the origin. Since the random vector Y=YXhas standard normal distribution, the first-order approximation of the failure probability is given by

PfPβαTY0=ΦβE5

where Φ(…) is the Laplace function.

The non-linear constrained optimisation problem (Eq. (3)) can be solved with many standard procedures as well as algorithms developed especially for this purpose, e.g. algorithm for the case of independent, non-normal random variables [16] and algorithm for problems with incomplete probability information [17].

All such solvers are iterative: for the assumed value of design point xk, the values of limit state function gxkand its gradient gxkare determined. Next, a new position of design point xk+1is derived, and the process continues until the convergence criteria are fulfilled. If the state of the analysed system is described by the performance function defined by analytical formula, then the gradient can be evaluated easily, and one of algorithms solving the optimisation problem (Eq. (3)) can be applied directly. Otherwise the stochastic finite element method should be applied in order to calculate the value of the limit state function and its gradient vector at following values of design points.

The first-order reliability index βand the failure probability PfΦβdepend on:

  • Parameters p=p1pdof the probability distributions of basic random variables, e.g. mean value, standard deviation, skewness or location, scale and shape parameters

  • Any deterministic parameters Θ=θ1θgdefining the form of the performance function gxθ1θg

Practical experience shows that the failure probability is usually a strongly non-linear function of the parameter θ, whereas the reliability index βis a rather linear function of the parameter θ. Thus the change in the failure probability due to the change of the parameter θcan be approximated as follows:

Pfθ+ΔθΦβΔβΦβΔθE6

The sensitivity measures of the first-order reliability index do not depend on the curvature of the limit state surface gx=0at the design point. Therefore, the application of sensitivity measures is limited to small changes of the values of the parameters.

The sensitivity of the first-order approximation of the failure probability PfΦβis directly related to the sensitivity of the reliability index β, since

dPf=φβE7

If θis a parameter of the limit state function gxθ, then derivative of the reliability index with respect to the parameter θis equal to

=1GyθθGyθE8

where vector Ycontains independent standard normal variables related to the vector of basic random variables by transformation Y=YX, and the limit state surface gxθ=0defined in the space Xhas been mapped into the surface Gyθ=0. Since the FORM index βis equal to the minimum distance between the origin of the Yspace and the limit state surface Gyθ=0, thus the design point yis laying on the limit state surface; see Figure 1:

β=GyθGyθyE9

Figure 1.

Illustration of sensitivity indices αi (modified from [14]).

The limit state surface in the Xspace of basic random variables gxθ=0does not depend on any parameter pikof a random variable Xiwith the distribution function Fixipik. However, the limit state surface Gyθ=0depends on parameter pikdue to the transformation Y=YX.

The derivative of the reliability index with respect to the parameter pikis given by relation

βpik=1βyTpikyE10

The derivative of the vector ywith respect to parameter pikhave to be evaluated for each specific transformation Y=YX.For details, see [14].

Sensitivity analysis shows how the uncertainty in the output response function of a system can be allocated to the different uncertainties in the basic variables. Sensitivity analysis is straightforward for FORM methodology. The influence of the basic random variable yi on the statistics of the response can be quantified by the sensitivity indices αi[18]:

αi=dyiforyi=yiE11

where yis the design point in the space of normalised reduced random variables.

For uncorrelated random variables, the sensitivity vector α coincides with the direction cosines vector of the random variables [18]. Illustration of sensitivity indices is given in Figure 1.

3.2. FORM versus Monte Carlo simulation

An alternative technique applied for probability estimation of risk or reliability is Monte Carlo simulation (MCS). For the purpose of the zero–one indicator-based MCS, Eq. (2) defining the failure probability is given as follows:

Pf=ΩFfX(x)dx=kRnkhK(k)dkE12

where the random vector K has the non-negative sampling density function hK(k) and is defined by the transformation

k=I(g(k)0)fX(k)hK(k)E13

and Iuis an indicator function:

Iu=1ifu00ifu>0E14

In this way the failure probability is equal to the expectation of random vector with the non-negative sampling density function hKk:

Pf=EKE15

The average of N simulated values of the random vector K is the estimator of the failure probability, which variance is equal to

VarPf̂=1NN1i=1NkiEKi2E16

Monte Carlo simulation technique is a powerful tool to calculate the probability of failure for the system described by non-continuous performance function as well as discrete random variables. However, the basic drawback of the MCS is long CPU time calculation, if the failure probability is of the orders 102106, since the sample size must be very large in order to obtain estimation of failure probability with low variance and narrow confidence interval. Various variance reduction techniques have been suggested to increase the efficiency of MCS. The basic idea is to assume a sampling density function hYythat reduces the variance of the estimator Pf̂. In the case of highly complicated systems, when time-consuming method must be applied to evaluate a single value of the limit state function, the MCS with the variance reduction technique is still an approach demanding a lot of computer time. Another drawback of the MCS, especially important, in the context of the chapter, is lack of the sensitivity analysis tools. It is simply impossible to run billions of simulation in order to study sensitivity of the system with respect to specific parameters.

4. Case-based risk/reliability studies of climate-related building performance

4.1. Building/Environment system performance

In the context of the ventilation design, air infiltration constitutes an important complement to air exchange. Furthermore, air infiltration can influence on the properties (thermal and structural) of building components. For some building technologies, e.g. lightweight timber frame with mineral wool filling, and loose mineral wool layers for roof insulation, the dependence of the thermal properties of building components on air infiltration can be observed; thus the interaction between, e.g. thermal transmittance and air infiltration should be taken into account. Therefore, it is important to apply a systemic approach to building/environment system performance taking into account different aspects of building physics. Due to the random natural driving forces governing the rate of air infiltration, the approach based on probabilistic methodology seems to be very well suited to handle these phenomena.

A building can be seen as a system transforming as well as resisting different loads (static and dynamic loads—caused by flow of air, heat and moisture) which is designed to ensure safe and comfortable living conditions inside the enclosure. The structure has to be designed in such a way that the possibilities of adverse consequences of this transformation, for example, loss of stability of the structure, inadequate ventilation or mould growth inside a building, have been minimalised. This systemic approach provides a proper theoretical tool for the analysis of the interrelations between the structure, its environment and its performance. An example of systemic model of a building, applicable in building physics studies, is shown in Figure 2 [19].

Figure 2.

Building/Environment system applied in a traditional building physics analysis [19].

The local environmental conditions interact with building structure to form a microclimate around a building. Sources of heat, air and moisture, including the products of HVAC systems as well as user behaviour, build up the internal load. Physical boundary conditions define the level of integration of the structure with the environment.

The output of the system can be described by the performance of the building (structure and enclosure). The performance can be considered in terms of safety, comfort and energy consumption and described by various parameters depending on physical conditions of the building structure and inside air. Those parameters should fulfil the performance requirements in order to prevent undesired performance (failure state) occurrence.

4.2. Case description

4.2.1. Description of the test house

The object of the study is a timber-framed low-rise naturally ventilated building with aspect ratio 2 and slope of the roof of 45° [20]. The building site in the district of Gothenburg has been considered and can be described as a semi-urban area with the surface roughness equal to 0.3 m. Example has been worked out for wind blowing from the west. It is assumed that the building is surrounded by other obstructions (other buildings, topography, vegetation, trees etc.) equivalent to half of its height. The following input data are used: volume of the house V=486m3, area of the building envelope A=336m2and internal temperature Tint=20°C.

The house was constructed in 1979 with the intention of using it for experimental studies in building physics with focus on ventilation and energy saving. The garage with doors facing south is located in the extended south part of the concrete cellar as shown in Figure 3.

Figure 3.

Object of the study ? the Building/Environment system [20].

4.2.2. Measurement programme

The following parameters have been measured, as is shown in Figure 4: (1) leakage characteristics of the house using blower door tests, (2) mean value of pressure difference across the six building components with Validyne pressure transducers, (3) wind speed and wind direction with the anemometer located on a small hill about 25mfrom the house, (4) internal and external temperatures and (5) limited number of tracer gas measurements of ACH. The measurement programme has been carried out during 8 months. As a result, hourly mean data have been registered.

The results of the pressure drop measurements have been used to validate the air infiltration through the envelope. An opening under the garage door has been treated separately in the calculation model for air change rate [20, 21].

4.3. Modelling of air change rate

The applied infiltration model takes into account the contribution of wind and stack effect to the total air change rate (ACH) in the following form [22]:

ACH=ACHs2+ACHw2E17

where ACHsis the air change rate caused by stack effect and ACHwis the air change rate caused by wind.

The model refers to low-rise building with light-weight construction, single ventilation zone, single temperature zone and steady-state conditions of air flow.

The infiltration model developed by Pietrzyk [20] indicates the air change rate ACHas a random function of three basic random variables: temperature difference, wind speed and wind direction. Wind direction is divided into eight sectors and is treated as a uniformly distributed within each sector. Finally, the air change rate conditioned by the wind direction sector is given by the following expression:

ACHd=s1T2+s2T+s3T1,5+wd,1vd4+wd,2vd2+wd,3vd3E18

where dis a wind direction sector; s1,s2,s3,wd,1,wd,2andwd,3are the deterministic coefficients related to the house dimensions, position of neutral pressure layer, level of external and internal pressure coefficients; Tis an ext.-int. temperature difference; and vdenotes wind speed.

Distributions of air change rate averaged over one hour (1-h) periods at a randomly chosen time in the year have been estimated with the help of the model described by Eq. (18). One-hour mean data ensure steady-state conditions of the airflow through the building envelope. Wind is the most important source of variations in the process of air exchange. However, according to wind energy spectrum presented in [23] for the frequency range 0.000140.0033cycles/hourrelated to time interval from 5 min to 2 h, the wind speed varies slightly. This range is called spectral gap. Measurements carried out for periods of that duration can be regarded as representing the steady-state conditions [24].

Performance criteria in terms of ACHshould take into account the minimum threshold evaluated with respect to unhygienic conditions. Then, probability of unsatisfactory performance is equal to PACH<threshold.

Figure 4 presents how the building response such as ACHdepends on the uncertain environmental conditions. The wind speed is traced from the meteorological station to the site and eventually to the building envelope which in turn influences the microclimatic conditions near to structure. The zone of wind-structure interaction is included in the model of designed system (see boundary conditions of the system presented by the solid lines). Serviceability performance due to wind action can be evaluated in terms of probability of undesired performance (failure). It is worth noticing that measurement data have been used to model the building performance as well as to validate the results of analysis carried out with the help of the established model.

Figure 4.

Transformation of uncertainty within the modelling of building performance.

The probability density function for air change rate as a function of basic random variables 1-h mean wind speed and 1-h mean temperature difference at time points chosen randomly during the year has been estimated with the help of the FORM sensitivity analysis for the performance function

gx1x2a=ACHx1x2aE19

where ACHx1x2is given by Eq. (18), x1=Tand x2=vd.

The parametric sensitivity analysis applied to FORM measures (reliability index or failure probability) is used in order to determine the probability density function for the random response ACH=ACHx1x2. The cumulative distribution function of random function ACH=ACHx1x2is actually equivalent to the probability of failure defined for the performance function Eq. (19):

FACHa=PACHx1x2a=Pgx1x2a0E20

Thus, the cumulative probability function can be estimated with the help of the FORM analysis:

FACHaΦβaE21

where Φuis the Laplace function and the reliability index βahas been determined for the limit state surface gx1x2=ACHx1x2a=0defined for a given value of parameter a.

Following the sensitivity measures presented earlier in the chapter, the probability density function of the random response ACHcan be estimated with the help of formula:

fACHaφβadβada=φβaGyafory=yE22

where φuis the probability density function of the standard Gaussian distribution, Gyais the limit state function in the space Y=YXof normalised random variables and yis the design point, i.e. the point on the surface Gya=at the shortest distance to the origin of the coordinate system. The value of probability density function of random response function fACH(ACH) can be obtained by means of FORM sensitivity analysis for consecutive values of parameter a; for details, see [14].

4.3.1. Wind transformation ? climate/local climate/microclimate

The input basic random variable for the infiltration model is wind speed in the vicinity of the building envelope. Wind speed and direction are usually measured at the meteorological stations. The mean value of 1-h mean wind speed can be evaluated from the mean value of 10-min mean wind speed obtained from meteorological station using the principle that the mean velocity increases by 5% when the averaging period is reduced from 1 h to 10 min. The transformation of the abovementioned data to the site of the building is often needed especially for wind that changes drastically due to the roughness of the ground surface.

The hourly mean wind speed vis assumed to follow the two-parameter Weibull distribution with probability density function as follows:

fvcλ=λcvcλ1expvcλE23

where λis a shape parameter and cis a scale parameter.

In general, parameters λ,cfor wind speed averaged over 1-h should be estimated for different wind direction sectors that shall result in quite different Weibull distributions due to two reasons: directional variability in the terrain surrounding the house and the predominance of certain wind directions. In particular, most building sites are subjected to sheltering effects from topography, trees and buildings. The roughness of the ground surface changes the mean wind speed and its turbulent characteristics and is described by the surface roughness height (aerodynamic roughness length) denoted z0. Roughness height depends on the mean element height of the roughness field. The results of laboratory measurements show that the value of z0is approximately equal to 1/30 of the height of the roughness elements. Table 1 presents the classification of roughness height for different types of surfaces with reference to the categories of terrain roughness used in Swedish Code [24].

Types of surface roughnessHeight z0 (m)Category
Calm open sea, water area0.0001I
Sand surface (smooth)0.001I
Snow surface0.003I
Bare soil0.005I
Airport runway area, mown grass0.01I
Farmland with very few buildings, trees, etc.0.03I
Farmland with open appearance0.05R
Farmland with closed appearance0.1II
Many trees and bushes0.2II
Shelter belts0.3II
Suburbs0.5II
City, forest1.0III

Table 1.

Roughness height for different types and categories of surfaces, acc. to Swedish Code [24].

Transformation of the wind speed between terrains of different surface roughness is possible due to the similarity theory [25], based on the equilibrium boundary layer height, which is according to [26] equal to 1200 m. The wind flows with the gradient velocity vgalong the isobars:

vg=uκlnu fz0AuE24

where z0is the surface roughness height (m); κis Karman’s constant, κ=0,4; Auis const.; assumed −1; fis Coriolis parameter (1/s), f=1,12104for latitude of order of 50° [25]; and vgis gradient velocity (m/s).

uis friction velocity depending on the surface shear stress τ0as given in Eq. (25):

u=τ0/ρE25

where τ0is surface shear stress (kg/ms2) and ρis air density (kg/m3).

The mean velocity profile uznear the ground, where zis the height above the ground (z<100m), can be expressed by the log-law model described by Eq. (24) assuming ideal conditions, i.e. the uniform height of roughness field and the neutrally stable atmosphere when thermal gradient is weak or absent.

uz=uκlnzz0E26

Eq. (26) is used for wind speeds greater than 10 m/s. For low wind speeds, the influence of thermal gradient both for unstable and stable atmosphere should be taken into account. The modified logarithmic formula can be found in [27].

The 10-min mean wind velocity measured at a meteorological station (usually at the level of z=10mabove the ground) for upwind surface roughness z0mcan be transformed to any other location described by upwind surface roughness height z0sthrough similarity of the wind speed at the gradient height for all terrain types [27]. Hence, the gradient velocity takes the same value for both locations and can be expressed by Eq. (27):

vg=umκlnum fz0mAu=usκlnus fz0sAuE27

where usis friction velocity at the building site (m/s) and umis friction velocity at the meteorological station (m/s).

The friction velocity at the meteorological station umis computed from Eq. (28), which has been derived on the basis of Eq. (27), by substituting the friction velocity at the site with friction velocity at the meteorological station.

um=umzκlnzz0mE28

where umzis 10-min mean wind speed measured at the meteorological station at the height z. The mean wind velocity uszat the site and at the height zcharacterised by upwind surface roughness z0scan be estimated from Eq. (26):

usz=usκlnzz0sE29

The ratio between wind velocity at the site and the wind velocity measured at the meteorological station denoted as ηis a function of the surface roughness z0mand z0s:

η=uszz0sumzz0mE30

It can be shown that non-linear relationship ηumcan be approximated with errors of order of 7% or less by a constant factor ηfor specified surface roughness at the building site. As the surface roughness appears in an implicit form in the expression for wind velocity (Eq. (29)), an analytical expression is not available. Instead, values of the factor ηhave been computed for different combinations of the surface roughness at the site and at the meteorological station (Table 2).

z0m = 0.01z0m = 0.05z0m = 0.3
z0s = 0.011.001.161.53
z0s = 0.050.861.001.31
z0s = 0.30.660.761.00

Table 2.

Values of the ratio η corresponding to different roughness conditions.

Simple wind transformation between categories of roughness is possible for z<20z0[28]. Thus, for z=10m, the transformation is valid for z0<0,5m. In the case of non-homogeneous upwind terrain, implementation for multiple roughness changes is required [26].

Values of wind speed measured at the meteorological station can be transformed using Eqs. (27)(29). The statistical parameters of wind speed evaluated for the site, i.e. the mean value μusand the standard deviation σus,can be easily evaluated on the basis of the mean wind speed measured at the meteorological station μumand the standard deviation of the wind speed measured at the meteorological station σum, since the wind speed at site is related the wind speed measured at the meteorological station by simple Eq. (30). Hence:

μus=ημumandσus=ησumE31

Concluding, the shape parameter λcharacterising the distribution of wind speed at the site of the building remains the same as the shape parameter λ for the meteorological station. The scale parameter for the distribution of wind speed at the site is equal to ηc(Eq. (23)).

Hence, probability density function of wind speed transformed to building location is given by 2p–Weibull probability model:

fvηcλ=ληcvηcλ1expvηcλE32

where cis a scale parameter and λis a shape parameter of the PDF of wind speed measured at the nearest meteorological station.

Modelling of microclimate around the structure takes into account the influence of structure form, orientation and the quality of the surrounding. Usually, the effect of wind pressure on the façade is estimated with the help of the tabulated values of wind pressure coefficients. In the analysed case pressure differences across the six building components on the structure were measured.

4.3.2. Air flow through the building envelope (influence of wind and temperature)

Some building performance aspects are dependent on the wind-structure interaction. Wind together with temperature difference causes airflow through building envelope.

The probability distribution model of external temperature depends on the specific geographical region. For temperate regions characterised by four seasons evenly distributed over the year, the normal (Gaussian) model with probability density function φTμTσT, given by Eq. (33), can be used for 1-h mean external temperature at “a random time” [29, 30]:

fTμTσT=1σT2πexp12TμTσT2E33

Also the full-scale measurements carried out near Gothenburg indicate [20] that the outdoor temperature can be approximated by the normal distribution.

Climatic data consist of 40-year record of observations made on meteorological stations at the airport in Säve, near Gothenburg. External temperature at the building site has been assumed to be equal to the temperature measured at the meteorological station, and its randomness is modelled by the normal distribution with the mean value of 11.1and the standard deviation of 6.1as shown in Figure 5.

Figure 5.

Normal PDF of ext. temperature T (*C) (left) and PDF of wind speed (m/s) for data coming from the Säve meteorological station (dashed line) and for local wind (solid line).

Temperature difference across the building envelope is also described by the normal PDF but shifted towards positive values by the average value of internal temperature.

The Weibull probability density function for a local wind speed has been evaluated on the basis of 10-min mean values of wind speed measured at meteorological station (see Figure 5). The meteorological station is located at the airport with assumed surface roughness 0.01. The ratio between wind velocity at the site and the velocity measured at the meteorological station has been calculated and is equal to 0.66 (Table 2). Probabilistic models of local wind speed together with wind speed measured at the meteorological station are given in Table 3.

Mean valueStandard dev.Scale parameterShape parameter
Meteo5.653.306.351.77
Local3.732.184.191.77

Table 3.

Stochastic parameters of the wind speed.

The probability density function of the random function ACH(Figure 6) has been evaluated using FORM approach (Eq. (22)). Probabilistic inference leads to the conclusion that the randomness of ACHis best described by the log-normal distribution with the mean value of 0.73 and the standard deviation of 0.38. Mean value and standard deviation are denoted, respectively, by μand σ. The PDF of the air change rate due to stack effect ACHsand the PDF of air change rate due to wind ACHware also shown in Figure 6. Randomness of air change rate due to stack effect can be described by the normal distribution whereas due to wind by the Weibull distribution skewed to the right.

Figure 6.

The probability density function for ACHsleft,ACHwmiddleandACH (right) established with the help of FORM analysis.

4.3.3. Sensitivity analysis of the probabilistic variability of air change rate with respect to the variability of wind and temperature

Dependence of ACHon the mean values of input variables follows the trends showed by sensitivity indices for individual variables (see Figures 79) [31]. For the values of ACHabove 1.0, where αΔTapproaches 0 and αvis equal to 1, the changes of reliability indices are dependent almost only on the changes of wind speed. Concluding, the wind velocity and temperature difference contribute significantly to the variability of the air change rate with sensitivity indices up to 0.8 for ΔT (for lower ACH) and up to near to 1 for wind speed (for higher ACH(Table 4)).

Figure 7.

Course of sensitivity index α,for variables∆Tandv.

Figure 8.

ACHsensitivity to theμ∆Tsolidand theσ∆TdashedleftandACHsensitivity to theμvsolidand theσv (dashed) (right).

Figure 9.

ACH sensitivity to scale (dashed) and shape (solid) parameter of Weibull distribution of wind speed.

ACH = 0.32ACH = 0.64ACH = 3.0
αΔT0.80.40.0
αv0.60.81.0

Table 4.

Some results from Figure 7.

Sensitivity of ACHdistribution with respect to mean values and standard deviations of input variables leads to the following conclusions: (1) strong dependence on wind variation, (2) temperature difference variations affect only low values of ACH(up to 0.4), (3) variations of ΔTaffect the lowest values of the ACHdistribution, and (4) variations of the wind speed are significant for performance studies of ACHwithin the whole range of wind speed values.

Figure 9 shows the measure of sensitivity of the PDF of ACHwith respect to scale or shape parameter of the Weibull distribution of wind speed for consecutive values of air change rate. The changes of shape parameter are the most important for the distribution of air change rate, especially for the threshold values close to the tail of distribution.

4.4. Probabilistic modelling of airflow-dependent thermal transmittance

For lightweight timber frame with mineral wool filling, the dependence of the thermal properties of building components on air infiltration is well acknowledged. An example is so-called dynamic wall [32], specially designed to save energy. In such a wall, the ventilation air passes through the insulation exchanging heat with a porous material reducing its conduction heat loss. The air entering the building is preheated by the conduction heat of the insulation (infiltration case), or the air leaving the building heats up the insulating material (exfiltration case) [33]. In the case of dynamic wall, the thermal transmittance becomes the most interesting parameter that can vary with the climatic data. Dynamic wall as a natural heat exchanger is a future of high-performance housing. The interaction between thermal transmittance and airflow through the components should be considered while calculating heat loss through a building envelope. A modelling approach based on probabilistic methods is proposed in [34].

Probabilistic model of dynamic Uvalue takes into account only some of the uncertainties related to the properties of the thermal insulation described by the thermal transmittance U0, the climatic load and the internal load coming from the building installations and occupants’ behaviour (ventilation strategy). The model described by Eq. (34) can be used to estimate a probability distribution of the dynamic Uvalue of the building envelope consisting of i‐thelements with total area of Atot:

U=14inNuiUi0AiE34

The Nusselt number Nuiis equal to 1 for the element without convection flow. In general, the value of the Nusselt number depends on the velocity of the airflow through the insulation, the direction of the flow and the thickness and the density of the insulation.

The example of approximation of the probability density functions of a dynamic Uvalue has been carried out with the help of FORM techniques. PDF of dynamic Uvalue has been evaluated using FORM sensitivity analysis (see Section 3.1.1). It depends on statistical parameters of the joint distribution of two random variables: thermal transmittance U0(varying with the temperature) and wind as well as buoyancy-driven airflow in terms of air change rate ACH(see Figure 10). It has been assumed that stochastic information is limited to the parameters of marginal probability density functions of those variables and the correlation coefficient between them.

Figure 10.

Probability density functions of dynamic U value for the cases of infiltration (left) and exfiltration (right) approximated for the building located near to Gothenburg for western winds.

Probability density functions of thermal transmittance depend on the direction of the airflow through the envelope as well as on the probability model of the air change rate. Respectively, to the contribution of the natural forces (wind, temperature) and mechanical forces, different probability distributions (normal, log-normal, Weibull and gamma) can be fitted to model randomness of the air change rate [35]. In general, the probability density functions of the dynamic Uvalue are skewed to the left—in the case of infiltration—and are skewed to the right, in case of exfiltration. The specific character of the relationship between Nusselt number and air change rate may explain these results. For the case of infiltration, the best fit according to the Kolmogorov-Smirnov test has been obtained for the Weibull distribution, while for the exfiltration case, the three-parameter gamma (or alternatively Gumbel) distribution has been obtained (see Figure 10).

The model could be further developed to include uncertainties due to other mechanisms and factors, e.g. influence of wind or radiation on external heat transfer coefficient or the influence of non-homogeneity of the material characteristics.

The probabilistic model for estimation of heat loss accounting for interactions between ventilation and transmission heat losses has been presented in [20, 33]. The model predicts the probability density function of the heat loss distribution over a specified period of time (e.g. a heating season) on the basis of the design parameters of the house, temperature characteristics of the site as well as the air change rate due to mechanical or natural ventilation. The probability of heat loss exceeding certain number of kW can be compared for different design options concerning various ventilation strategies (natural or/and mechanical ventilation) and various transmittance properties (tight envelope contra dynamic wall) of the building envelope. Hence, rational engineering decisions promoting low-energy solution contributing to climate change mitigation can be taken into account in the design process.

5. Conclusions

Risk analysis together with appropriate tools can support building design strategies concerning climate change mitigation and adaptation. Lower levels of uncertainty can be handled by means of risk analysis based on system’s risk or reliability estimations. In the case of higher order of uncertainties [36], other strategies could be developed based on the concept of resilience.

Risk analysis of building performance enables the selection of the best design based on comparison of probabilities of undesired performance estimated for alternative design solutions. Systemic approach gives opportunity to identify important relationships between variables. For example, air infiltration as a result of climate/structure interaction may be a significant variable in the thermal performance of building envelope. However, in order to handle the whole complexity of the real system, multivariable decision models for different design solutions should be further developed.

The examples of dynamic Uvalues resulting in the different characters of distribution models for the cases of infiltration and exfiltration show that the probabilistic methods and tools can be effectively used to establish the probabilistic characteristics typical for the combination of the important variables influencing climate-structure interaction.

The sensitivity measures are important in the case of risk or reliability-based design. Sensitivity analysis of the distribution of a response variable (random function) with respect to the basic random variables and its parameters is straightforward for FORM methodology, whereas it is not easy in the case of the Monte Carlo simulation. The results of case studies show that the air change rate distribution depends on the temperature difference Tand the wind velocity significantly. Dependence of PDF model of ACH on the mean value of input variables is similar in the trends for both studied variables: temperature and wind speed. Sensitivity analysis of ACHprobability distribution model to standard deviations of input variables shows the high contribution of wind speed and limited to low values of ACH(up to 0.4) influence of temperature difference.

Approximate transformation of wind speed data from the meteorological station to a specific location, where analysed building is situated, can be carried out by multiplying wind speed by a constant factor η(with 7% error or less), established for a specific ranges of roughness conditions. The transformation of the probabilistic model of 10-min mean wind speed from meteorological station to the probabilistic model of hourly mean speed for the site of the building results in change of the scale parameter, while the shape parameter remains the same. The form of PDF for ACHas well as reliability index is sensitive to the value of the shape parameter of the Weibull distribution of wind speed and much less sensitive to the scale parameter. Hence, the transformation of probabilistic model of wind speed to the local site seems to be robust for the analysed case.

As it was shown, the sensitivity analysis has helped to understand the relationships between model inputs. It can also help to test the model outcome in terms of its robustness in the presence of uncertainty. The probabilistic risk analysis along with the effective tools for sensitivity analysis can be used to support design decisions and also to develop better models for evaluation of building performance.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Krystyna Pietrzyk and Ireneusz Czmoch (February 28th 2018). On Risk and Reliability Studies of Climate-Related Building Performance, Risk Assessment, Valentina Svalova, IntechOpen, DOI: 10.5772/intechopen.71684. Available from:

chapter statistics

234total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Machinery Safety Requirements as an Effective Tools for Operational Safety Management

By Hana Pacaiova

Related Book

First chapter

Introductory Chapter: Earthquakes, Life at Risk

By Valentina Svalova

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us