Next Article in Journal
Land Degradation Monitoring in the Ordos Plateau of China Using an Expert Knowledge and BP-ANN-Based Approach
Next Article in Special Issue
Experimental Analysis and Process Modeling of Carbon Dioxide Removal Using Tuff
Previous Article in Journal
On the Acceptance and Sustainability of Renewable Energy Projects—A Systems Thinking Perspective
Previous Article in Special Issue
Management Optimization of the Luminous Flux Regulation of a Lighting System in Road Tunnels. A First Approach to the Exertion of Predictive Control Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Fluid Dynamic Modelling of Thermal Periodic Stabilized Regime in Passive Buildings

Department of Astronautical, Electrical and Energy Engineering (DIAEE), Sapienza University of Rome, Via Eudossiana 18, 00184 Rome, Italy
*
Author to whom correspondence should be addressed.
Sustainability 2016, 8(11), 1172; https://doi.org/10.3390/su8111172
Submission received: 30 September 2016 / Revised: 27 October 2016 / Accepted: 10 November 2016 / Published: 12 November 2016

Abstract

:
The periodic stabilized regime is the condition where the temperature of each point of a certain environment varies following a periodic law. This phenomenon occurs in many practical applications, such as passive or ancient buildings not equipped with Heating, Ventilating and Air Conditioning HVAC systems and located in latitudes where the temperature greatly varies with Earth’s daily cycles. Despite that, the study of transient phenomena is often simplified, i.e., considering negligible the thermal response of the indoor microclimate. An exact solution to enclosures whose microclimate is free to evolve under a periodic stabilized regime does not exist nowadays, also from an analytical point of view. The aim of this study is to parametrically analyze the thermal variations inside a room when a transient periodic temperature is applied on one side. The phenomenon has been numerically studied through Computational Fluid Dynamics (CFD) and analytically validated using a function that reproduces the daily variation of the outdoor temperature. The results of this research would lay the groundwork to develop analytical correlations to solve and predict the thermal behavior of environments subject to a periodic stabilized regime.

1. Introduction

The indoor microclimate of buildings is greatly influenced by outdoor climate conditions and by the ability of the architectural envelope to dampen them. These physical phenomena involve different heat transfer processes (coupling convection/conduction) which can affect several aspects related to sustainable architecture and design of buildings and cities, such as the reduction of energy demand [1,2,3,4,5,6], the decrease of the environmental impact [7,8] and the analysis of the urban canyon effect [9,10]. Differentially heated enclosures have been widely studied in the past [11,12,13,14,15]. However, the complexity of the problem has led the researchers to apply several simplifications in the analysis of these phenomena, such as setting the outdoor conditions as steady or considering indoor microclimate as a controlled environment (where the interior conditions are set, and, if required, dynamically controlled by an HVAC system). Furthermore, the ISO 13786:2007(E) [16] supplies the procedure for the calculation of the dynamic thermal behavior of building components imposing the temperature on one side of the element as constant.
Time-periodic thermal conditions should be taken into account especially when the outdoor temperature greatly varies during the day [17,18,19,20,21]. One application is the so-called periodic stabilized regime, which is the condition where the average temperature varies with a periodic law in each point of a certain environment. This phenomenon finds practical applications in buildings not characterized by HVAC systems, as passive, ancient and/or historical buildings, where the temperature is free to evolve without thermal constraints.
The problem has been analytically studied taking into account the case of the semi-infinite solid [22,23,24]. An analytical solution for the transient temperature distribution in a semi-infinite solid is known only for three surface conditions: constant surface temperature, constant surface heat flux, and surface convection. To the best of our knowledge, there is not an exact analytical law to predict the indoor thermal behavior in enclosures where a varying temperature is applied at least on one external surface and the indoor temperature is free to evolve.
Time-periodic natural convection in 2D enclosures has been also numerically examined in the past literature [25,26,27,28,29,30,31,32]. However, these studies are more focused on the analysis of the fluid field inside the enclosure rather than supplying new analytical relations to solve the physical phenomenon of the periodic stabilized regime. This problem could be studied by coupling experimental measures and numerical simulations and comparing the results. Despite that, an experimental approach is not always applicable, for example when the case study is too complex to be experimentally duplicated or when no more existing buildings or sites are analyzed [33,34].
The aim of this study is to parametrically analyze the thermal variations inside a room when a transient periodic temperature is applied on one side. The results of this research would lay the groundwork to develop analytical correlations to solve and predict the thermal behavior of environments subject to time-dependent thermal forcing and where the indoor temperature is free to evolve, matching together all the approaches used up to now.

2. Materials and Methods

Nowadays, Computational Fluid Dynamic (CFD) is one of the main tools used in several studies to analyze heat transfer processes involving convection and conduction. However, these studies did not supply analytical correlations to predict the thermal behavior of enclosures where the temperature is free to evolve under a periodic stabilized regime. In this work, the indoor temperature trends inside the model of a simple room has been numerically investigated through CFD. Parametric analyses have been carried out by changing the thermo-physical and boundary conditions of the case study. Several sets of simulations have been, indeed, taken into account in order to understand what the thermal behavior of the indoor environment is when geometry and thermal inputs vary. Four sensitivity tests have been carried out changing one parameter at a time. All the case studies analyzed are summarized in Table 1.
In the first two cases (Case 1 and Case 2), the thermo-physical properties of the envelope have been varied, while in the last two cases (Case 3 and Case 4), different thermal inputs have been applied on the sidewalls of the model. The varying profile of the external daily temperature has been reproduced through a sinusoidal trend on the external surface of the enclosure.

2.1. Analytical Statement

In the study of heat transfer by conduction, the Fourier’s equation admits analytical solutions only for well-defined problems, usually characterized by simple geometries, simple boundary conditions, materials with uniform properties and independent temperature [22,23,24].
Many heat transfer problems are time-dependent. It is possible to analyze transient conduction phenomena referring to a semi-infinite solid that extends its dimensions to infinity with the exception of one of them.
A typical case of transient conduction is described by a body of a wide thickness delimited by a flat surface where a transient periodic temperature is imposed on. The other surfaces of the solid are at a distance large enough not to consider their influence. The temperature in all the points of the body varies with the same periodic law of the imposed temperature after a certain time. This condition is called periodic stabilized regime.
Suppose to analyze a semi-infinite wall, whose material is homogeneous and isotropic, delimited by a plan surface. The thermo-physical characteristics of the medium (cp, λ, ρ) are independent from temperature and time. Imagine to impose on its surfaces a sinusoidal temperature in the normal x direction, which follows the analytical Equation (1):
T ( x , τ ) = T m + Θ sin ( ω τ )
After an initial transient period, the temperature reaches a periodic stabilized regime. Considering a time τ in the periodic stabilized regime, since the temperature does not depend on the spatial coordinates y and z, the solution to the Fourier’s equation is (2):
T ( x , τ ) = T m + Θ e β x sin ( ω τ β x )
Transient conduction in a semi-infinite solid can be solved also in three different conditions on the surface [24]: a constant temperature TsTi (Figure 1a); a constant heat flux q 0 (Figure 1b); an exposure to a fluid whose temperature is TTi and a convective coefficient h (Figure 1c).
These three cases can be solved respectively through the following equations [24]:
T ( x , τ ) T s T i T s = e r f ( x 2 α τ )
T ( x , τ ) T i = 2 q 0 ( α τ / π ) 1 / 2 λ exp ( x 2 4 α τ ) q 0 x λ e r f c ( x 2 α τ )
T ( x , τ ) T i T T i = e r f c ( x 2 α τ ) [ exp ( h x λ + h 2 α τ λ 2 ) ] [ e r f c ( x 2 α τ + h α τ λ ) ]
The mathematical calculations to obtain the Equations (3)–(5) can be found in [22,23,24].
The thermo-physical characteristics of the solid body, such as the diffusivity, the phase lag and the thermal damping factor, are noticeable in all the equations. It means that the initial sinusoidal temperature is attenuated and delayed proportionally to the thickness and the physical and thermal properties of the wall.
To the best of our knowledge, an analytical solution to enclosures subject to a periodic stabilized regime does not exist nowadays.

2.2. Numerical Approach

The commercial software Ansys Fluent v. 14.5 [35] has been used for the transient simulations. The governing equations used for the gas phase are the unsteady form of the continuity, compressible Navier–Stokes, and energy equations for a perfect gas with the Sutherland model for the viscosity and temperature-dependent properties. The equations that govern the phenomenon (momentum (6), continuity (7) and energy conservation (8)) can be expressed as follows:
u ¯ j u ¯ i x j = 1 ρ p ¯ x i + μ ρ 2 u ¯ i x i x j x j ( u ¯ i u ¯ j ) + g i
u ¯ i x i = 0
u ¯ i T ¯ x i + x i ( k T ¯ x i ) = 0
where i = 1, 2, 3 for the three-dimensional (3D) case and i = 1, 2 for the two-dimensional (2D) case. The Fourier law has been used for the analysis of the solid body.
Time-accurate solutions have been obtained using the implicit coupled scheme with second-order accuracy in time with the SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) segregated scheme, used for the gas phase until convergence was achieved. For time-accurate solutions, the fluxes for all equations at the cell faces have been interpolated by using the second-order upwind scheme. Pressure equation and the Laplace equation for the conduction have been also computed by using second-order accuracy. The Boussinesq approach has been applied to the model and RNG k-ε turbulence model [36,37] has been chosen to study the fluid turbulence after the calculation of the dimensionless groups for the natural convection (Ra = 2.54 × 1010). Basing on convergence and computational stability criteria, the RNG k-ε model has been considered the most suitable for this case in terms of indoor flow field characterization [38]. The simulations have run with the unsteady regime with an initial time step size between 10−7 and 10−5 s and it has been increased at 20 s after the residuals have reached the stability. This time step has been obtained after several tests based on numerical stability and numerical residual convergence. Convergence criteria have been reached when all the residuals were below 10−6.

3. Validation Cases

Heat conduction problems are often solved through many assumptions and simplifications, such as the steady state condition. Despite that, several applications involve time-dependent variation of the temperature, especially when it changes in a periodic manner. As previously said, there is not an exact analytical law to predict the indoor thermal behavior in enclosures where the indoor temperature is free to evolve subject to a periodic stabilized regime.
The numerical analyses carried out in this research have been validated through the simulation of simple cases where the analytical solution is known. Two case studies have been taken into account [34]:
  • Case A: heat transfer through a solid where a time-dependent sinusoidal temperature is applied on one side;
  • Case B: heat transfer through a solid and an adjacent fluid domain.
The solid body has been assumed isotropic and homogeneous. The material of the solid is the aluminum, which is characterized by: ρ = 2719 kg/m3, cp = 871 J/(kgK), λ = 202.4 W/(mK).
The geometry of Case A is represented by a rectangular domain having dimensions of L1 = 5.0 m and L2 = 3.0 m (Figure 2). A time-dependent sinusoidal temperature Tsin has been applied on one sidewall in the normal x direction, which follows the analytical Equation (1), where Tm = 298.15 K, Θ0 = 5 K and ω = 7.3 × 10−5 rad/s.
A fixed temperature Twall = 298.15 K has been applied on the opposite sidewall, while the top and the bottom walls have been considered as adiabatic.
Since the geometry is very simple, a mesh with a structured grid, composed by 10,500 cells (interval size 0.035 m), has been chosen for the discretization.
Figure 3 shows the temperature trends in the center of the solid domain (x = 2.5 m, y = 1.5 m) in both numerical and analytical analyses. The analytical data are derived (2).
The results show that the two temperatures have the same trends and almost the same values. The comparison between analytical and numerical data has been performed also through the root mean square error (RMSE) in a time period of 24 h after the stabilization of the simulation. The obtained RMSE value is 0.027 K. Therefore, CFD simulations can be considered accurate enough for the analysis of the transient heat conduction through a solid.
The geometry of Case B is represented by a solid domain with a rectangular shape and a thickness of s = 1 m, while the fluid domain dimensions are L1 = 5.0 m and L2 = 3.0 m (Figure 4). The heat exchange coefficient has been imposed equal to 10 W/(m2K). The temperature of the fluid domain has been kept constant and equal to Tfluid = 293.15 K, while the solid temperature Tin has been left free to evolve, starting from a temperature of 215.15 K. The top and the bottom walls have been considered adiabatic.
As the previous case, the geometry is very simple. For this reason, a mesh with a structured grid, composed by 12,700 cells (interval size 0.035 m), has been chosen for the discretization.
Figure 5 shows the temperature trends in the center of the solid domain (x = 0.5 m, y = 1.5 m) in both numerical and analytical analyses. In this case, it is possible to use the Equation (5), where λ, α, h, and T are assumed constant. The assumption of constant h and T is almost never encountered in real cases, as for the one-dimensional application.
Analytical and numerical temperature trends greatly differ only before the stabilization of the simulation (around 2000 h). The root mean square error (RMSE) has been evaluated in a time period of 400 h after 2000 h of simulated time and the obtained value is 0.0001 K. It is possible to assume that the numerical model approximates accurately the analytical one also in this case.
The CFD model has been considered validated, accepting that the obtained values are affected by a minimum error. Therefore, it has been applied to a similar real model where a solid envelope adjacent to the flow domains is subject to thermal sinusoidal boundary conditions varying with time.

4. 3D-2D Simplification

The periodic stabilized regime has been numerically applied to a model of a simple closed room with a rectangular plan of 5 × 5 m (W × D) and H = 3 m. The room model is influenced by an outdoor thermal forcing applied on a solid wall, while the opposite wall is adjacent to an environment at a fixed temperature. The other walls have been considered adjacent to environments having the same temperature of the room; therefore, it has been possible to assume that there is not heat transfer through those walls (adiabatic process).
The 3D geometry of the room has been simplified in a 2D computational domain representing the longitudinal section on the central xy plane (Figure 6).
In this work, this simplification has been possible since the considered room has a preferential symmetry plan along the z axis. The choice of a 2D model allows lower computational costs without compromising the accuracy of the results.
In order to demonstrate the effectiveness of this choice, both the 3D and 2D models have been simulated with the same boundary conditions and simulation settings. The results have been shown in terms of temperature trends at the end of the solid domain (probe_wall) and in the center of the fluid domain (probe_fluid) and in terms of contours of the velocity magnitude in the central plane xy (Figure 6).
The two computational models have been discretized using an orthogonal grid (Figure 7).
A grid independence test has been performed for the 2D model in order to demonstrate the independence of the results from the mesh spacing. A fixed temperature has been applied on the external wall of the solid domain and the static temperature in the middle of the fluid domain (probe_fluid) has been chosen as a parameter for the test. Assuming the finest mesh considered, Mesh D (mesh spacing equal to 0.0175 m), as the pivot case, it is possible to observe that the differences percentage Δ% with the other three meshes are always very low (Table 2). The mesh with the lower error, Mesh A (mesh spacing equal to 0.035 m), has been used for these CFD analyses and the 3D model has been discretized considering the same interval size of the mesh, with an amount of 155,344 cells.
Both models are composed by two computational domains, one for the solid wall and one for the fluid inside the room. The same boundary conditions have been taken into account (Figure 8). A time-dependent sinusoidal temperature Tsin has been applied on the external surface of the solid domain following the analytical Equation (1), while a fixed temperature Twall has been applied on the right sidewall of the fluid domain. The top and the bottom walls have been considered adiabatic assuming those surfaces adjacent to environments at the same temperature, as usually happens in a room of a building. Also, the other two sides of the 3D model have been considered adiabatic. Neither a fixed temperature nor a convection coefficient have been applied to the indoor environment. Therefore, the temperature is free to evolve under the influence of the external temperature.
Figure 9 shows the temperature trends in probe_wall and probe_fluid for the two models.
The results show that the temperatures have almost the same values and trends. The root mean square error (RMSE) is 0.2 K for probe_fluid and 0.1 K for probe_wall. Therefore, it is possible to assume that the values obtained from the 2D numerical model describe the thermal phenomena with the same accuracy of the 3D model.
This assumption has been supported also by the analysis of the fluid field and the convective motions (Figure 10 and Figure 11).
The air motions have been analyzed in two different time steps which represent the highest and the lowest peaks of the sinusoidal temperature Tsin. The results are almost the same for the two models. The fluid fields show in both time steps two areas characterized by an acceleration of the air. The most noticeable difference between the two models is the air velocity of the acceleration zones, which reaches highest values in the 2D model. This phenomenon is due to an increase of kinetic energy in the 2D model, where the air motions are distributed only in two directions. However, these differences have been considered small enough to be ignored. Therefore, the two models show comparable results also from a fluid dynamic point of view.

5. Analytical Evaluation of the Solid Envelope

The analysis of building materials and building elements is a crucial aspect in the evaluation of the effects of outdoor microclimatic conditions on indoor environment. For this reason, the temperature trends at the end of the different analyzed walls have been analytically calculated considering a semi-infinite solid and applying Equation (2). This calculation has been performed for the walls in Case 1 and Case 2 in order to compare the variation between an ideal geometry and a more realistic one. In both cases, a sinusoidal temperature Tsin given by the Equation (1), where Tm = 298.15 K, Θ0 = 5 K and ω = 7.3 × 10−5 rad/s, has been applied on the sidewall of the solid domain.
Three different materials have been taken into account: aluminum, masonry and wood. Thermo-physical properties of the materials (density, specific heat capacity, conductivity) have been derived from UNI 10351:1994 [39] and UNI EN ISO 10456:2007 [40] (Table 3).
These materials have been chosen because of their different capability in transmitting, damping and delaying a thermal input. Considering a fixed thickness of the solid wall s = 0.5 m, the properties of the wall vary as shown in Table 4.
The temperature trends at the end of the wall have been shown in Figure 12 for the three considered materials.
Another important aspect that influences the thermal behavior of the building envelope is the thickness of the building elements. Three different thicknesses s of the solid domain have been taken into account (0.1 m, 0.5 m and 1.0 m), while the material has been set as aluminum (Table 3). The thermo-physical properties of the solid wall vary as shown in Table 5.
The temperature trends at the end of the wall for the three considered thicknesses s are shown in Figure 13.
The reported analytical data underline how and to what extent the external temperature is perceived on the internal surface of the wall in the semi-infinite solid case. Knowing Θ0 and the thermal damping factor σ of each building element, it has been possible to calculate the semi-amplitude of the sinusoidal wave applied on the external surface once it is perceived on the internal surface.
The solid wall of aluminum delays and damps a little the temperature Tsin and it is possible to expect that indoor temperature largely follows the sinusoidal trend of the outdoor thermal forcing. On the contrary, the solid walls composed by masonry and wood greatly impact on the outdoor temperature. For this reason, it is perceived almost as a constant temperature. In the same way, the higher the thickness s is, the less the indoor temperature varies following the external sinusoidal trend.

6. CFD Results: Sensitivity Analysis for Different Geometrical Parameters and Boundary Conditions

6.1. Materials Sensitivity Test

The materials described in Table 3 have been used for the solid wall of Case 1. The sinusoidal temperature Tsin described in the previous chapter has been applied on the external side of the solid wall; the temperature of the right wall of the fluid environment is Twall = 293.15 K, while the top and the bottom walls have been set adiabatic.
Figure 14 shows the temperature trends at the end of the wall and in the center of the fluid domain for different solid wall materials compared with the external sinusoidal temperature.
The numerical results confirm the analytical data shown in Table 3 and Table 4, underlining the relationship between thermo-physical properties of building materials and architectural elements.
The characteristics of the sinusoidal wave at the end of the wall and in the center of the fluid domain have been derived from the numerical results and analyzed in-depth in Table 6.
As expected, the semi-amplitude ϴ of the thermal wave decreases between x = 0.0 m and x = 2.5 m. This decrement is easily shown in the aluminum and the masonry cases, while it is very slight in the wood wall.
After passing through the masonry and wood walls, the sinusoidal wave is almost totally damped. On the contrary, in the case of the aluminum, the fluid environment damps the sinusoidal wave more than the wall. In this case, the thermal wave Tsin quickly passes through the solid domain and it is subject to a small damping, because of the high thermal conductivity of the material. Therefore, the resulting thermal wave dissipates a great part of its energy in the fluid environment, which is characterized by a lower thermal conductivity. For the other two materials, Tsin arrives into the fluid environment with almost completely lowered oscillations.
This behavior is noticeable also in the variation of Tm for the three materials, which decreases with the decrease of the thermal conductivity. Higher energy levels correspond to higher Tm, as it happens for aluminum.
Finally, regarding the phase lag, masonry and wood have similar behavior. The phase lag, indeed, remains almost the same between the end of the wall and the center of the fluid domain in both cases. On the contrary, the phase lag increases for aluminum. This phenomenon is probably due to the influence of Tsin on the indoor environment. In the cases where masonry and wood have been used, the entering temperature is almost constant and has almost the same values. Therefore, it is easy to suppose that the fluid field is almost the same. In the aluminum case, indeed, the great oscillation of the entering temperature causes different fluid motions, which take to different thermo-fluid dynamic behavior.
The comparison between the numerical results in probe_wall and the analytical data (Table 4 and Figure 12) shows how the fluid environment is able to limit the energy losses from the solid wall. In the comparison with the analytical results, the phase lag is lower and the thermal damping factor and the semi-amplitude are higher in the numerical simulations. This phenomenon is emphasized by the decrease of the thermal conductivity of the material, which is related to a lower energy level dissipated by the wall and, consequently, to lower heat levels transferred to the fluid environment. This suggests that the analytical solution for the semi-infinite solid is too approximated for the case of a wall adjacent to a fluid environment.

6.2. Thickness Sensitivity Test

Three different wall thicknesses s have been modeled in Case 2. The applied thermal boundary conditions are the same used in Case 1, while the material has been set as aluminum (Table 3); this choice has been made since the other two materials greatly damp and delay the external sinusoidal wave, approximating it to a constant temperature. For this reason, the indoor thermal variations are too small to be perceived. The obtained temperature trends are shown in Figure 15.
The characteristics of the sinusoidal wave at the end of the wall and in the center of the fluid domain have been derived from the numerical results (Table 7).
Comparing the numerical results and the analytical data (Table 5 and Figure 13), the trend of thermal damping factor at thickness variation is the most evident phenomenon. Based on the analytical data, the sinusoidal wave is expected to be highly damped with the increase of thickness. Despite that, the thermal damping factor is almost the same in the numerical results in both probe_wall (90%–100%) and probe_fluid (45%–50%). Furthermore, unlike Case 1, the numerical values of phase lag in the same point of the wall are higher than the analytical ones, and this difference is emphasized consequently to the increase of the thicknesses s, passing from an increase of thermal lag of 6% (s = 0.1 m) to 80% (s = 1.0 m). On the contrary, thermal damping factor and semi-amplitude are lower compared to the analytical ones. Also in this case, the difference increases with the increase of thickness s. These differences are caused by the convective phenomena in the fluid environment, which are not taken into account in the analytical calculation.
The numerical results also suggest that the higher the thickness is, the lower the energy levels of the entering wave are. This phenomenon causes lower heat exchanges among the fluid particles of the enclosure, damping less the temperature. On the contrary, basing on the same assumption, a lower thickness allows higher energy levels of the entering wave and a consequent higher heat exchange and a higher damping of the temperature.
Also in this case, it is possible to assume that the analytical solution for the semi-infinite solid is too approximated for the case of a wall adjacent to a fluid environment.

6.3. External Sinusoidal Temperature Tsin Sensitivity Test

The third and fourth sensitivity analyses concern the variation of the thermal boundary conditions, while the thermo-physical properties of the solid wall have been left unaltered. In the third case the sinusoidal wave Tsin applied on the external surface of the solid wall has been changed varying the initial Tm,0, while Θ0 = 5 K and ω = 7.3 × 10−5 rad/s have been left unaltered in all the simulations. Three Tm,0 have been considered: 295.15 K, 298.15 K and 301.15 K. The temperature on the right wall of the fluid environment is Twall = 293.15 K, while the top and the bottom walls have been set adiabatic. The material has been set as aluminum (Table 3) and the thickness of the wall s is equal to 1.0 m.
The obtained temperature trends are shown in Figure 16.
The characteristics of the sinusoidal wave at the end of the wall and in the center of the fluid domain have been derived from the numerical results (Table 8).
Varying Tm,0 of the sinusoidal wave Tsin, the temperature trends in probe_wall and probe_fluid are damped and delayed in the same way, since the thermo-physical properties of the building materials and elements are the same. The only noticeable difference is the variation of Tm which increases with the Tm,0 of the external temperature. Furthermore, Tm does not increase proportionally to Tm,0 since it is influenced also by Twall. The results underline an interesting relation applicable in the analyzed case: the temperature Tm of the environment can be calculated, as the average between the external Tm,0 and the fixed Twall basing on the distance from the considered walls. This phenomenon is not computable when the analytical solution is applied.

6.4. External Fixed Temperature Twall Sensitivity Test

Finally, the influence of the temperature Twall on the interior environment has been investigated. Thermo-physical properties of materials and building elements are the same as Case 3, while the sinusoidal temperature Tsin follows Equation (2), where Tm = 295.15 K, Θ0 = 5 K and ω = 7.3 × 10−5 rad/s. Three Twall have been considered: 293.15 K, 303.15 K and 323.15 K.
The obtained temperature trends are shown in Figure 17.
The characteristics of the sinusoidal wave at the end of the wall and in the center of the fluid domain have been derived from the numerical results (Table 9).
The results show the same physical phenomena observed in Case 3. The indoor temperatures, indeed, are characterized by the same thermal damping factors and phase lags, and their Tm is function of the average temperature of Tsin and Twall.
This consideration is helpful when materials with low thermal damping factors and high phase lags are taken into account (such as the previously analyzed wood and masonry). In that case, indeed, the external sinusoidal temperature is perceived inside as a constant temperature equal to Tm and the indoor temperature is equal to the average of the all the external temperatures considered. As for the previous cases, this phenomenon is not computable when the analytical solution is applied.

7. Conclusions

The indoor microclimate of buildings is influenced by several parameters. Especially when the indoor conditions are subject to an outdoor transient thermal forcing and they are free to evolve in a periodic stabilized regime, the heat transfer phenomena are very complex. For this reason, several studies use some simplifications in the investigation of these physical problems. To the best of our knowledge, an exact solution to solve this phenomenon does not exist also from the analytical point of view. The aim of this research is to propose a first approach to solving the thermal phenomena inside environments where the temperature is free to evolve, such as buildings not equipped with HVAC systems.
The coupling of conduction and convection inside a simplified room has been investigated through Computational Fluid Dynamics. Different preliminary analyses have been taken into account in order to validate the numerical results and to evaluate possible discrepancies between the 3D geometry of the room and its simplification in a 2D numerical model. The numerical model has been studied taking into account four configurations and performing parametric analysis in order to understand the response of the indoor environment when thermal inputs and thermo-physical properties of the envelope change.
In the first and the second configurations (Case 1 and Case 2), the thermo-physical properties of the solid domain have been varied considering three different construction materials (aluminum, masonry and wood) and three different thicknesses of the wall (0.1 m, 0.5 m and 1 m.) The results of these analyses underline a strong dependence of the indoor temperature on thermo-physical properties of the solid envelope. Furthermore, the entering thermal wave is also greatly influenced by the energy exchanges with the fluid particles in the indoor environment. These heat exchanges are shown in terms of decrease of the thermal damping factor related to the variation of thermal conductivity.
In the third and fourth sensitivity analyses (Cases 3 and Case 4), the thermal boundary conditions have been varied considering three Tm,0 of the sinusoidal temperature (295.15 K, 298.15 K and 301.15 K) and three different Twall (293.15 K, 303.15 K and 323.15 K). In these cases, the results underline that the thermal damping factor and the phase lag do not change with these variations, while the Tm in the middle of the fluid environment (probe_fluid) changes proportionally to the average of Tm,0 and Twall.
The comparison between the obtained numerical data and the analytical ones for the four cases have shown that the case of the semi-infinite solid is too approximate to correctly describe the temperature trend inside a fluid environment adjacent to a wall in a periodic stabilized regime. More in detail, the analytical calculation does not take into account the convective phenomena in the fluid environment resulting from the heat levels exchanged with the wall. This situation is clear when the thermo-physical characteristics of the solid vary, especially in Case 2, where the phase lag and thermal damping factor do not change as expected from the analytical calculations. For this reason, further developments of the research are needed in order to formulate more accurate predictions of the periodic stabilized regime.
This work is a first step to better understanding these physical phenomena and it would lay the groundwork to develop analytical correlations to predict the trend of the temperature inside an environment free from thermal control and subject only to the varying external thermal conditions.

Author Contributions

The study has been conceived and designed by Fabio Nardecchia. Fabio Nardecchia, Francesca Pagliaro and Benedetta Mattoni have elaborated the case studies, performed the numerical simulations and analyzed analytical and numerical data. Lucia Cellucci has done the graphical works. Fabio Nardecchia, Francesca Pagliaro and Benedetta Mattoni have written the paper. English corrections have been undertaken by Francesca Pagliaro and Benedetta Mattoni. Franco Gugliermetti and Fabio Bisegna have revised the research work and the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

cpspecific heat at constant pressure (J/kgK)
erfGaussian error function
erfccomplementary error function, defined as e r f c 1 e r f
Ddepth (m)
ggravitational acceleration (m/s2)
hconvection heat transfer coefficient (W/m2K)
Hheight (m)
kturbulent kinetic energy (m2/s2), defined as k = u i u j ¯ 2
ppressure (Pa)
RaRayleigh number
sthickness (m)
Ttemperature (K)
ufluid velocity component (m/s)
Uthermal transmittance (W/m2K), defined as U = 1 1 h 1 + s λ + 1 h 2
Wwidth (m)
x,y,zspatial coordinates
Greek letters
αthermal diffusivity (m2/s), defined as α = λ ρ c p
βxphase shift (rad), defined as β = π α τ 0
Θsemi-amplitude of a sinusoidal temperature
λthermal conductivity (W/mK)
μdynamic viscosity (kg/s m)
ρmass density (kg/m3)
σthermal damping factor, defined as σ = e β x
τtime (s)
τrphase lag (h), defined as τ r = β x ω = x 2 τ 0 α π
ωangular velocity (rad/s)
Subscript
0initial condition
Eexternal
iindex; internal
ssurface
mmean value
free stream condition
Overbar
¯average conditions

References

  1. Fantozzi, F.; Galbiati, P.; Leccese, F.; Salvadori, G.; Rocca, M. Thermal analysis of the building envelope of lightweight temporary housing. J. Phys. Conf. Ser. 2014, 547, 491–502. [Google Scholar] [CrossRef]
  2. Carbonara, E.; Tiberi, M.; Astiaso Garcia, D. Analysis of energy performance improvements in Italian residential buildings. Energy Procedia 2015, 82, 855–862. [Google Scholar] [CrossRef]
  3. Evangelisti, L.; Guattari, C.; Gori, P.; de Lieto Vollaro, R. In situ thermal transmittance measurements for investigating differences between wall models and actual building performance. Sustainability 2015, 7, 10388–10398. [Google Scholar] [CrossRef]
  4. Burattini, C.; Nardecchia, F.; Bisegna, F.; Cellucci, L.; Gugliermetti, F.; de Lieto Vollaro, A.; Salata, F.; Golasi, I. Methodological approach to the energy analysis of unconstrained historical buildings. Sustainability 2015, 7, 10428–10444. [Google Scholar] [CrossRef]
  5. Galli, G.; Vallati, V.; Recchiuti, C.; de Lieto Vollaro, R.; Botta, F. Passive cooling design options to improve thermal comfort in an Urban District of Rome, under hot summer conditions. Int. J. Eng. Technol. 2013, 5, 4495–4500. [Google Scholar]
  6. De Lieto Vollaro, R.; Evangelisti, L.; Carnielo, E.; Battista, G.; Gori, P.; Guattari, C.; Fanchiotti, A. An integrated approach for an historical buildings energy analysis in a smart cities perspective. Energy Procedia 2014, 45, 372–378. [Google Scholar] [CrossRef]
  7. Cumo, F.; Astiaso Garcia, D.; Stefanini, V.; Tiberi, M. Technologies and strategies to design sustainable tourist accommodations in areas of high environmental value not connected to the electricity grid. Int. J. Sustain. Dev. Plan. 2015, 10, 20–28. [Google Scholar] [CrossRef]
  8. Pagliaro, F.; Cellucci, L.; Burattini, C.; Bisegna, F.; Gugliermetti, F.; de Lieto Vollaro, A.; Salata, F.; Golasi, I. A methodological comparison between energy and environmental performance evaluation. Sustainability 2015, 7, 10324–10342. [Google Scholar] [CrossRef]
  9. De Lieto Vollaro, A.; Galli, G.; Vallati, A.; Romagnoli, R. Analysis of thermal field within an urban canyon with variable thermophysical characteristics of the building’s walls. J. Phys. Conf. Ser. 2015. [Google Scholar] [CrossRef]
  10. De Lieto Vollaro, A.; Galli, G.; Vallati, A. CFD analysis of convective heat transfer coefficient on external surfaces of buildings. Sustainability 2015, 7, 9088–9099. [Google Scholar] [CrossRef]
  11. Acharya, S.; Tsang, C.H. Influence of wall conduction on natural convection in an inclined square enclosure. Wärme Stoffübertragung 1987, 21, 19–30. [Google Scholar] [CrossRef]
  12. Hanjalić, K.; Kenjereš, S.; Durst, F. Natural convection in partitioned two-dimensional enclosures at high Rayleigh numbers. Int. J. Heat Mass Transf. 1996, 39, 1407–1427. [Google Scholar] [CrossRef]
  13. Bilgen, E. Natural convection in enclosures with partial partitions. Renew. Energy 2002, 26, 257–270. [Google Scholar] [CrossRef]
  14. Mobedi, M. Conjugate natural convection in a square cavity with finite thickness horizontal walls. Int. Commun. Heat Mass 2008, 35, 503–513. [Google Scholar] [CrossRef]
  15. Varol, Y.; Oztop, H.F.; Koca, A. Effects of inclination angle on natural convection in composite walled enclosures. Heat Transf. Eng. 2011, 32, 57–68. [Google Scholar] [CrossRef]
  16. International Organization for Standardization (ISO). Thermal Performance of Building Components—Dynamic Thermal Characteristics—Calculation Methods; BS EN ISO 13786:2007; International Organization for Standardization: Geneva, Switzerland, 2007. [Google Scholar]
  17. Nardecchia, F.; Bisegna, F.; Gugliermetti, F. A novel approach to CFD analysis of the urban environment. J. Phys. Conf. Ser. 2015, 655, 1–11. [Google Scholar] [CrossRef]
  18. Salata, F.; Golasi, I.; de Lieto Vollaro, E.; Bisegna, F.; Nardecchia, F.; Coppi, M.; Gugliermetti, F.; de Vollaro, A.L. Evaluation of different urban microclimate mitigation strategies through a PMV analysis. Sustainability 2015, 7, 9012–9030. [Google Scholar] [CrossRef]
  19. Nardecchia, F.; Bisegna, F.; Gugliermetti, F. How temperature affects the airflow around a single-block isolated building. Energ Build. 2016, 118, 142–151. [Google Scholar] [CrossRef]
  20. Evangelisti, L.; Guattari, C.; Gori, P. Energy retrofit strategies for residential building envelopes: An Italian case study of an early-50s building. Sustainability 2015, 7, 10445–10460. [Google Scholar] [CrossRef]
  21. Ciampi, M.; Fantozzi, F.; Leccese, F.; Tuoni, G. On the optimization of building envelope thermal performance. Multi-layered walls design to minimize heating and cooling plant intervention in the case of time varying external temperature fields. Civ. Eng. Environ. Syst. 2003, 20, 231–254. [Google Scholar] [CrossRef]
  22. CarsLaw, H.S.; Jaeger, J.C. Conduction of Heat in Solids, 2nd ed.; Oxford University Press: London, UK, 1959. [Google Scholar]
  23. Gröber, H.; Erk, S. Fundamentals of Heat Transfer, 3rd ed.; McGraw-Hill: New York, NY, USA, 1961. [Google Scholar]
  24. Incropera, F.P.; de Witt, D.P. Fundamentals of Heat and Mass Transfer, 5th ed.; John Wiley and Sons: New York, NY, USA, 1985. [Google Scholar]
  25. Kazmierczak, M.; Chinoda, Z. Buoyancy-driven flow in an enclosure with time periodic boundary conditions. Int. J. Heat Mass Transf. 1992, 35, 1507–1518. [Google Scholar] [CrossRef]
  26. Lage, J.L.; Bejan, A. The resonance of natural convection in an enclosure heated periodically from the side. Int. J. Heat Mass Transf. 1993, 36, 2027–2038. [Google Scholar] [CrossRef]
  27. Antohe, B.V.; Lage, J.L. A dynamic thermal insulator: Inducing resonance within a fluid saturated porous medium enclosure heated periodically from the side. Int. J. Heat Mass Transf. 1994, 37, 771–782. [Google Scholar] [CrossRef]
  28. Lakhal, E.K.; Hasnaoui, M.; Vasseur, P.; Bilgen, E. Natural convection in a square enclosure heated periodically from part of the bottom Wall. Numer. Heat Transf. A Appl. 1995, 27, 319–333. [Google Scholar] [CrossRef]
  29. Abourida, B.; Hasnaoui, M.; Douamna, S. Transient natural convection in a square enclosure with horizontal walls submitted to periodic temperatures. Numer. Heat Transf. A Appl. 1999, 36, 737–750. [Google Scholar]
  30. El Ayachi, R.; Raji, A.; Hasnaoui, M.; Bahlaoui, A. Combined effect of radiation and natural convection in a square cavity differentially heated with a periodic temperature. Numer. Heat Transf. A Appl. 2008, 53, 1339–1356. [Google Scholar] [CrossRef]
  31. Wang, G.; Zeng, M.; Ren, Y.C.; Ozoe, H.; Wnag, Q.W. Transient heat flux measurement of natural convection in an inclined enclosure with time-periodically-varying wall temperature. Exp. Therm. Fluid Sci. 2011, 35, 105–111. [Google Scholar] [CrossRef]
  32. Raji, A.; Hasnaoui, M.; Firdaouss, M.; Ouardi, C. Natural convection heat transfer enhancement in a square cavity periodically cooled from above. Numer. Heat Transf. A Appl. 2013, 63, 511–533. [Google Scholar] [CrossRef]
  33. Pagliaro, F.; Bukowiecki, E.; Bisegna, F.; Gugliermetti, F. The architecture of warehouses: A multidisciplinary study on thermal performances of Portus’ roman store buildings. J. Cult. Herit. 2014, 16, 560–566. [Google Scholar] [CrossRef]
  34. Pagliaro, F.; Nardecchia, F.; Bisegna, F.; Gugliermetti, F. CFD analysis for the validation of archaeological hypotheses—The indoor microclimate of ancient storage-rooms. J. Archaeol. Sci. 2016, 73, 107–119. [Google Scholar] [CrossRef]
  35. ANSYS Inc. ANSYS Fluent Software Package, User’s Manual, version 14.5; ANSYS Inc.: Canonsburg, PA, USA, 2012. [Google Scholar]
  36. Choudhury, D. Introduction to the Renormalization Group Method and Turbulence Modeling; Technical Memorandum TM-107; Fluent Inc.: Canonsburg, PA, USA, 1993. [Google Scholar]
  37. Yakhot, V.; Orszag, S.A. Renormalization group analysis of turbulence: Basic theory. J. Sci. Comput. 1986, 1, 1–51. [Google Scholar] [CrossRef]
  38. Roy, J.C.; Boulard, T. CFD prediction of the natural ventilation in a tunnel-type greenhouse: Influence of wind direction and sensibility to turbulence models. Acta Hortic. 2010, 691, 457–464. [Google Scholar] [CrossRef]
  39. Ente Nazionale Italiano di Unificazione (UNI). Building Materials. Thermal Conductivity and Vapour Permeability; UNI 10351:1994; Ente Nazionale Italiano di Unificazione: Milan, Italy, 1994. [Google Scholar]
  40. International Organization for Standardization (ISO). Building Materials and Products—Hygrothermal Properties—Tabulated Design Values and Procedures for determining Declared and Design Thermal Values; UNI EN ISO 10456:2007; International Organization for Standardization: Geneva, Switzerland, 2007. [Google Scholar]
Figure 1. Transient conduction in a semi-infinite solid depending on different surface conditions: constant temperature TsTi (a); constant heat flux q 0 (b); exposure to a fluid whose temperature is TTi and a convective coefficient h (c). Original figure from [24].
Figure 1. Transient conduction in a semi-infinite solid depending on different surface conditions: constant temperature TsTi (a); constant heat flux q 0 (b); exposure to a fluid whose temperature is TTi and a convective coefficient h (c). Original figure from [24].
Sustainability 08 01172 g001
Figure 2. Case A—Geometry, mesh and boundary conditions.
Figure 2. Case A—Geometry, mesh and boundary conditions.
Sustainability 08 01172 g002
Figure 3. Case A: Temperature trends comparison between CFD and analytic analyses.
Figure 3. Case A: Temperature trends comparison between CFD and analytic analyses.
Sustainability 08 01172 g003
Figure 4. Case B: Geometry, mesh and boundary conditions.
Figure 4. Case B: Geometry, mesh and boundary conditions.
Sustainability 08 01172 g004
Figure 5. Case B: Temperature trends comparison between CFD and analytic analysis.
Figure 5. Case B: Temperature trends comparison between CFD and analytic analysis.
Sustainability 08 01172 g005
Figure 6. Geometry: 3D model (a); 2D model (b).
Figure 6. Geometry: 3D model (a); 2D model (b).
Sustainability 08 01172 g006
Figure 7. Mesh: 3D model (a); 2D model (b).
Figure 7. Mesh: 3D model (a); 2D model (b).
Sustainability 08 01172 g007
Figure 8. Geometry, mesh and boundary conditions of the 2D model.
Figure 8. Geometry, mesh and boundary conditions of the 2D model.
Sustainability 08 01172 g008
Figure 9. Temperature trends for the 3D model (symbols) and the 2D model (lines).
Figure 9. Temperature trends for the 3D model (symbols) and the 2D model (lines).
Sustainability 08 01172 g009
Figure 10. Contours of the velocity magnitude expressed in m/s in the highest peak of the sinusoidal temperature Tsin: 3D model (a); 2D model (b).
Figure 10. Contours of the velocity magnitude expressed in m/s in the highest peak of the sinusoidal temperature Tsin: 3D model (a); 2D model (b).
Sustainability 08 01172 g010
Figure 11. Contours of the velocity magnitude expressed in m/s in the lowest peak of the sinusoidal temperature Tsin: 3D model (a); 2D model (b).
Figure 11. Contours of the velocity magnitude expressed in m/s in the lowest peak of the sinusoidal temperature Tsin: 3D model (a); 2D model (b).
Sustainability 08 01172 g011
Figure 12. Analytical temperature trends at the end of the wall (s = 0.5 m) for three different materials.
Figure 12. Analytical temperature trends at the end of the wall (s = 0.5 m) for three different materials.
Sustainability 08 01172 g012
Figure 13. Analytical temperature trends at the end of the wall for three different thicknesses s.
Figure 13. Analytical temperature trends at the end of the wall for three different thicknesses s.
Sustainability 08 01172 g013
Figure 14. Case 1: Temperature trends in probe_wall (a) and probe_fluid (b).
Figure 14. Case 1: Temperature trends in probe_wall (a) and probe_fluid (b).
Sustainability 08 01172 g014
Figure 15. Case 2: Temperature trends in probe_wall (a) and probe_fluid (b).
Figure 15. Case 2: Temperature trends in probe_wall (a) and probe_fluid (b).
Sustainability 08 01172 g015
Figure 16. Case 3: Temperature trends in probe_wall (a) and probe_fluid (b).
Figure 16. Case 3: Temperature trends in probe_wall (a) and probe_fluid (b).
Sustainability 08 01172 g016
Figure 17. Case 4: Temperature trends in probe_wall (a) and probe_fluid (b).
Figure 17. Case 4: Temperature trends in probe_wall (a) and probe_fluid (b).
Sustainability 08 01172 g017
Table 1. Case studies summary.
Table 1. Case studies summary.
Case NameSensitivity Test
Case 1Material of the solid domain
Case 2Thickness of the solid domain
Case 3External sinusoidal temperature Tsin
Case 4External fixed temperature Twall
Table 2. Mesh sensitivity test for the 2D model.
Table 2. Mesh sensitivity test for the 2D model.
Mesh AMesh BMesh CMesh D
Cells of the mesh116,80292073046,720
Interval size0.0350.070.140.0175
Δ%0.02%0.17%0.30%-
Table 3. Thermo-physical properties of the chosen materials.
Table 3. Thermo-physical properties of the chosen materials.
AluminumMasonryWood
ρ27191700700
λ202.40.580.17
cp8718502310
Table 4. Thermo-physical properties of the solid wall as the material changes.
Table 4. Thermo-physical properties of the solid wall as the material changes.
AluminumMasonryWood
U5.800.970.33
σ72.180.860.01
τr1.2518.1835.22
ϴ3.610.040.00
Table 5. Thermo-physical properties of the solid wall as the thicknesses changes.
Table 5. Thermo-physical properties of the solid wall as the thicknesses changes.
s = 0.1 ms = 0.5 ms = 1.0 m
U5.865.805.71
σ93.6972.1852.09
τr0.251.252.49
ϴ4.683.612.60
Table 6. Case 1: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
Table 6. Case 1: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
AluminumMasonryWood
Tsinprobe_wallprobe_fluidprobe_wallprobe_fluidprobe_wallprobe_fluid
σ-99.1549.141.520.760.010.01
τr-0.551.2318.0518.6031.1431.26
ϴ5.04.962.460.080.040.000.00
Tm298.10298.13295.67296.39294.76294.96294.04
Table 7. Case 2: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
Table 7. Case 2: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
s = 0.1 ms = 0.5 ms = 1.0 m
Tsinprobe_wallprobe_fluidprobe_wallprobe_fluidprobe_wallprobe_fluid
σ-99.6649.4199.1549.1493.8446.45
τr-0.040.740.551.231.662.43
ϴ54.982.474.962.464.692.32
Tm298.15298.13295.67298.13295.67298.12295.66
Table 8. Case 3: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
Table 8. Case 3: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
Tm,0 = 295.15Tm,0 = 298.15Tm,0 = 301.15
Tsinprobe_wallprobe_fluidprobe_wallprobe_fluidprobe_wallprobe_fluid
σ-93.8646.2993.8446.4593.8446.48
τr-1.592.351.662.431.592.36
ϴ54.692.324.692.324.692.32
Tm298.15295.14294.17298.12295.66301.10297.17
Table 9. Case 4: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
Table 9. Case 4: Characteristics of the sinusoidal wave in probe_wall and probe_fluid.
Twall = 293.15 KTwall = 303.15 KTwall = 323.15 K
Tsinprobe_wallprobe_fluidprobe_wallprobe_fluidprobe_wallprobe_fluid
σ-93.8646.2993.8445.9993.8351.02
τr-1.592.351.602.401.612.24
ϴ54.692.314.692.304.692.55
Tm298.15295.14294.17295.20299.19295.32309.20

Share and Cite

MDPI and ACS Style

Nardecchia, F.; Mattoni, B.; Pagliaro, F.; Cellucci, L.; Bisegna, F.; Gugliermetti, F. Computational Fluid Dynamic Modelling of Thermal Periodic Stabilized Regime in Passive Buildings. Sustainability 2016, 8, 1172. https://doi.org/10.3390/su8111172

AMA Style

Nardecchia F, Mattoni B, Pagliaro F, Cellucci L, Bisegna F, Gugliermetti F. Computational Fluid Dynamic Modelling of Thermal Periodic Stabilized Regime in Passive Buildings. Sustainability. 2016; 8(11):1172. https://doi.org/10.3390/su8111172

Chicago/Turabian Style

Nardecchia, Fabio, Benedetta Mattoni, Francesca Pagliaro, Lucia Cellucci, Fabio Bisegna, and Franco Gugliermetti. 2016. "Computational Fluid Dynamic Modelling of Thermal Periodic Stabilized Regime in Passive Buildings" Sustainability 8, no. 11: 1172. https://doi.org/10.3390/su8111172

APA Style

Nardecchia, F., Mattoni, B., Pagliaro, F., Cellucci, L., Bisegna, F., & Gugliermetti, F. (2016). Computational Fluid Dynamic Modelling of Thermal Periodic Stabilized Regime in Passive Buildings. Sustainability, 8(11), 1172. https://doi.org/10.3390/su8111172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop