PDF

Thermal Performances of Roller Shutters
Introduction
During the design of a building, environmental issues have gained considerable influence in the entire project. One of the first concerns is to improve thermal performances. In this process, simulation software provide key tools for modeling thermal losses and performances in the building.
The international standard ISO 10077-2:2012 (Ref. 1) deals with thermal performances of windows, doors, and shutters. It provides computed values of the thermal characteristics of frame profiles in order to validate a simulation software.
COMSOL Multiphysics successfully passes the entire benchmark. This document describes two test cases of ISO 10077-2:2012 related to roller shutters only. Other test cases from this standard are available in the following applications:
Figure 1: 3D representation of the roller shutter box with shutters inside.
Model Definition
On each test case, a shutter section separates a hot internal side from a cold external side. After solving a model, two quantities are calculated and compared to the normative values:
air cavities
The roller shutter structure contains many cavities. The purpose is to ensure thermal insulation. According to the ISO 10077-2:2012 standard, cavities are modeled in different ways depending on their shapes.
The heat flow rate in cavities is represented by an equivalent thermal conductivity, keq, which includes the heat flow by conduction, convection, and radiation. It also depends on the geometry of the cavity and on the adjacent materials. The definition of keq is detailed in the next paragraphs.
Cavities are divided into three types:
unventilated cavities, completely closed or connected either to the exterior or to the interior by a slit with a width not exceeding 2 mm;
slightly ventilated cavities, connected either to the exterior or to the interior by a slit greater than 2 mm but not exceeding 10 mm;
well-ventilated cavities: corresponding to a configuration not covered by one of the two preceding types, it is assumed that the whole surface is exposed to the environment so that boundary conditions are applied to (see the Boundary conditions section below for more information).
For the main cavity within a roller shutter box, these rules are slightly different (see Figure 2):
Figure 2: Opening of a roller shutter box.
If e1 + e3 ≤ 2 mm, the cavity is considered as unventilated.
If etot ≤ 35 mm, the cavity is considered as slightly ventilated.
If etot > 35 mm, the cavity is considered as well-ventilated.
Unventilated Rectangular Cavity
For an unventilated rectangular cavity, the equivalent thermal conductivity is defined by:
where d is the cavity dimension in the heat flow rate direction, and R is the cavity thermal resistance given by:
Here, ha is the convective heat transfer coefficient, and hr is the radiative heat transfer coefficient. These coefficients are defined by:
where:
C1 = 0.025 W/(m·K)
C2 = 0.73 W/(m2·K4/3)
ΔT is the maximum surface temperature difference in the cavity
σ = 5.67·10-8 W/(m2·K4) is the Stefan-Boltzmann constant
Tm is the average temperature on the boundaries of the cavity
E is the intersurface emittance, defined by:
ε1 and ε2 are the surface emissivities (both are equal to 0.90 in this model)
F is the view factor of the rectangular section, defined by:
d is the cavity dimension in the heat flow rate direction
b is the cavity dimension perpendicular to the heat flow rate direction
Slightly Ventilated Rectangular Cavities
For a slightly ventilated cavity, the equivalent thermal conductivity is twice that of an unventilated cavity of the same size.
Nonrectangular Cavities
Nonrectangular cavities are transformed into rectangular cavities of same area and aspect ratio according to defined rules in ISO 10077-2:2012 presented below. Then, keq is evaluated following one of the two previous rectangular cases.
Figure 3: Nonrectangular cavity transformation.
Figure 3 shows a nonrectangular cavity of area A. Then, d and b are the depth and the width (in accordance with the direction of the heat flow) of the smallest rectangle than can contain of the nonrectangular cavity. The equivalent rectangular cavity, of size b × d and area A must satisfy:
Hence, b and d are given by:
Boundary conditions
The heat flux conditions for internal and external sides are given by the Newton’s law of cooling:
where Text is the exterior temperature (Text = Ti = 20°C for the internal side and Text = Te = 0°C for the external side). The standard defines thermal surface resistance, Rs, which is related to the heat transfer coefficient, h, by:
Internal and external thermal surface resistances are not equal.
Description of the two Applications
Figure 4 and Figure 5 depict the geometry of each model. Unventilated cavities are red-numbered while slightly ventilated cavities are green-numbered. Adiabatic boundaries are represented with striped rectangles.
Application 1: Roller Shutter Box
The first application studies the heat conduction in a roller shutter box. The main structure is made of polyvinyl chloride (PVC) which has a low thermal conductivity k of 0.17 W/(m·K). Inside the box, there is an insulation panel which has a very low thermal conductivity of 0.035 W/(m·K).
In this application, there are thirty-eight cavities. Thirty-seven of them are not connected to the exterior so they are considered as unventilated cavities. The main cavity is considered as slightly ventilated because of the large opening in the box (15 mm).
Figure 4: Geometry of the roller shutter box.
Application 2: PVC Shutter Profile
This application studies the heat conduction in a PVC shutter profile. The shutter is made of two PVC blocks which have a thermal conductivity of 0.17 W/(m·K).
In this application there are five cavities. They are not connected to the exterior so they are considered as unventilated cavities.
Figure 5: Geometry of the PVC shutter profile.
Results and Discussion
Temperature profiles
The temperature profiles for each model are shown in Figure 6 and Figure 7.
Figure 6: Temperature distribution in the roller shutter box.
Figure 7: Temperature distribution in the PVC shutter profile.
Quantities of interest
The quantities of interest are the following:
where is the heat flow rate through the shutter (in W/m), Te = 0°C is the external temperature and Ti = 20°C is the internal temperature.
where l is the projected length of the internal section perpendicularly to the heat flow direction (expressed in meters).
Table 1 and Table 2 compare the numerical results of COMSOL Multiphysics with the expected values provided by ISO 10077-2:2012.
Table 1: comparison between expected values and Computed values of Quantities in Application 1.
L2D (W/(m·K))
Uf (W/(m2·K))
Table 2: comparison between expected values and Computed values of Quantities in Application 2.
L2D (W/(m·K))
Uf (W/(m2·K))
The maximum permissible differences to pass this test case are 3% for the thermal conductance and 5% for the thermal transmittance. The measured values are completely coherent and meet the validation criteria.
Reference
1. European Committee for Standardization, ISO 10077-2:2012, Thermal performance of windows, doors and shutters – Calculation of thermal transmittance – Part 2: Numerical method for frames, 2012.
Application Library path: Heat_Transfer_Module/Buildings_and_Constructions/roller_shutter_thermal_performances
Modeling Instructions
Root
Start by opening the following prepared file. It already contains global definitions, geometries, local variables, selections, operators and material properties.
1
From the File menu, choose Open.
2
Roller Shutter Box
Roller Shutter Box (comp1)
In the Model Builder window, expand the Roller Shutter Box (comp1) node.
Definitions (comp1)
Variables 1
Define the thermal conductance of the section for the postprocessing part as follows.
1
In the Model Builder window, expand the Roller Shutter Box (comp1)>Definitions node, then click Variables 1.
2
In the Settings window for Variables, locate the Variables section.
3
Note that the heat flow rate through the internal and through the external boundaries are equal (in absolute value) because other boundaries are considered adiabatic.
4
In the Model Builder window, collapse the Roller Shutter Box (comp1)>Definitions node.
Heat Transfer in Solids and Fluids (ht)
Fluid 1
1
In the Model Builder window, expand the Heat Transfer in Solids and Fluids (ht) node, then click Fluid 1.
2
As there is no convection, a second order discretization of the temperature is set for better accuracy.
3
In the Model Builder window, click Heat Transfer in Solids and Fluids (ht).
4
In the Settings window for Heat Transfer in Solids and Fluids, click to expand the Discretization section.
5
From the Temperature list, choose Quadratic Lagrange.
Heat Flux 1
1
In the Physics toolbar, click  Boundaries and choose Heat Flux.
2
In the Settings window for Heat Flux, locate the Boundary Selection section.
3
From the Selection list, choose Exterior Side.
4
Locate the Heat Flux section. Click the Convective heat flux button.
5
In the h text field, type 1/Rse.
6
In the Text text field, type Te.
Heat Flux 2
1
In the Physics toolbar, click  Boundaries and choose Heat Flux.
2
In the Settings window for Heat Flux, locate the Boundary Selection section.
3
From the Selection list, choose Interior Side.
4
Locate the Heat Flux section. Click the Convective heat flux button.
5
In the h text field, type 1/Rsi.
6
In the Text text field, type Ti.
7
In the Model Builder window, collapse the Heat Transfer in Solids and Fluids (ht) node.
Study 1
The heat flow rate through the interior (or exterior) side of the section needs to be determined to calculate the thermal conductance of the section. In order to have enough precision in this value, the default relative tolerance of the solver has already been modified to 106. To access to this value, expand the Solver 1 node and click on the Stationary Solver 1 node. In the Stationary Solver settings window, locate the General section.
1
In the Home toolbar, click  Compute.
Results
Temperature (ht)
A Global Evaluation node is added in order to calculate the thermal conductance of the section and the thermal transmittance of the frame.
Thermal Properties, Roller Shutter Box
1
In the Model Builder window, expand the Results>Derived Values node.
2
Right-click Derived Values and choose Global Evaluation.
3
In the Settings window for Global Evaluation, type Thermal Properties, Roller Shutter Box in the Label text field.
4
Locate the Expressions section. In the table, enter the following settings:
5
Click  Evaluate.
Table
1
Go to the Table window.
The results should be close to the expected values in Table 1.
Results
Surface
1
In the Model Builder window, expand the Results>Temperature (ht) node, then click Surface.
2
In the Settings window for Surface, locate the Expression section.
3
From the Unit list, choose degC.
4
In the Temperature (ht) toolbar, click  Plot.
5
Click the  Zoom Extents button in the Graphics toolbar.
The current plot group shows the temperature distribution; compare with Figure 6.
The same simulation method is applied to the other benchmark. The instructions below describe the steps to achieve the calculations.
PVC Shutter Profile
Roller Shutter Box (comp1)
In the Model Builder window, collapse the Roller Shutter Box (comp1) node.
PVC Shutter Profile (comp2)
In the Model Builder window, expand the PVC Shutter Profile (comp2) node.
Definitions (comp2)
Variables 2
1
In the Model Builder window, expand the PVC Shutter Profile (comp2)>Definitions node, then click Variables 2.
2
In the Settings window for Variables, locate the Variables section.
3
4
In the Model Builder window, collapse the PVC Shutter Profile (comp2)>Definitions node.
Heat Transfer in Solids and Fluids 2 (ht2)
Fluid 1
1
In the Model Builder window, expand the Heat Transfer in Solids and Fluids 2 (ht2) node, then click Fluid 1.
2
As there is no convection, a second order discretization of the temperature is set for better accuracy.
3
In the Model Builder window, click Heat Transfer in Solids and Fluids 2 (ht2).
4
In the Settings window for Heat Transfer in Solids and Fluids, click to expand the Discretization section.
5
From the Temperature list, choose Quadratic Lagrange.
Heat Flux 1
1
In the Physics toolbar, click  Boundaries and choose Heat Flux.
2
In the Settings window for Heat Flux, locate the Boundary Selection section.
3
From the Selection list, choose Exterior Side.
4
Locate the Heat Flux section. Click the Convective heat flux button.
5
In the h text field, type 1/Rse.
6
In the Text text field, type Te.
Heat Flux 2
1
In the Physics toolbar, click  Boundaries and choose Heat Flux.
2
In the Settings window for Heat Flux, locate the Boundary Selection section.
3
From the Selection list, choose Interior Side.
4
Locate the Heat Flux section. Click the Convective heat flux button.
5
In the h text field, type 1/Rsi.
6
In the Text text field, type Ti.
7
In the Model Builder window, collapse the Heat Transfer in Solids and Fluids 2 (ht2) node.
Study 2
In the Home toolbar, click  Compute.
Results
A Global Evaluation node is added in order to calculate the thermal conductance of the section and the thermal transmittance of the frame.
Thermal Properties, PVC Shutter Profile
1
In the Results toolbar, click  Global Evaluation.
2
In the Settings window for Global Evaluation, type Thermal Properties, PVC Shutter Profile in the Label text field.
3
Locate the Data section. From the Dataset list, choose Study 2/Solution 2 (4) (sol2).
4
Locate the Expressions section. In the table, enter the following settings:
5
Click  Evaluate.
Table
1
Go to the Table window.
The results should be close to the expected values in Table 2.
Results
Surface
1
In the Model Builder window, expand the Results>Temperature (ht2) node, then click Surface.
2
In the Settings window for Surface, locate the Expression section.
3
From the Unit list, choose degC.
4
In the Temperature (ht2) toolbar, click  Plot.
5
Click the  Zoom Extents button in the Graphics toolbar.
The current plot group shows the temperature distribution; compare with Figure 7.