Computer Simulation of Temperature Distribution during Cooling of the Thermally Insulated Room

This paper is devoted to modelling of temperature distribution and its time evolution in rooms with specific thermal insulation and heat transfer for different external conditions. The simulation results should help to design the room architecture and wall materials to reduce energy losses due to heating or cooling, and to increase the inside thermal comfort. For this purpose, a methodological procedure using real data processing in the COMSOL Multiphysics modelling environment and spatial visualization of temperature evolution is proposed. This paper describes a mathematical model for simulation of the temperature evolution inside a space with thermally insulated walls under selected outside conditions. Computer simulations are then used to assess the temperature distribution inside the room and the heat flow through the room walls. Results of the simulations are used for subsequent determination of the time needed for the desired decrease of air temperature inside the tested room during its cooling due to the low ambient temperature, which is related to the thermal stability of the building, specific heat capacity, and thickness of the thermal insulation. Under the studied conditions, the time to reach the temperature drops by 20 percent in a room with windows was from 1.4 to 1.8 times lower than that in the room without windows. The proposed methodology shows the flexibility of computer modelling in the design of insulated building systems. The mesh density testing was performed by comparing the air temperature evolution in the model of the selected mesh density and the model with its maximum value enabled by the size of computer memory. The maximum temperature deviation calculated for the mesh of the presented model was 0.57%.


Introduction
Urban areas are developing rapidly from a technological viewpoint [1].Nowadays, over half of the world's population live in cities.Furthermore, the ongoing trend toward urbanisation will create a staggering increase in city populations of 2.8 billion by 2050.The constant increase in population living in cities has led to increased energy demands from households.Apart from industry and transport, households are currently major energy consumers.Heating accounts for the biggest energy and fuel consumption in European households.Its consumption ranged between 60 and 70% in most European Union counties in 2014, which is a leading cause of global warming.
Reducing the carbon dioxide emissions produced by households is another reason for creating the concept of smart cities, which use smart, meaningful and gentle modern technology and approaches to improve life, to manage energy more effectively, to conserve natural resources and ensure energy sustainability [2].An integral part of smart cities will be smart buildings, which will utilise monitoring devices that will track usage, and empower users and service providers to better control and reduce energy demand [3][4][5].
Minimising the energy demands of buildings needs to be focused on when designing new buildings and during house renovations [6,7].The thermal stability of buildings depends on both the thermal stability of the building's shell and the heat capacity of the internal structures and equipment.The heat resistance of buildings can also be significantly influenced by the appropriate thermal insulation of their walls by a suitable heat-insulating material that can reduce heat losses in winter and undesirable heating in summer [8,9].Phase change materials can also be used as passive energy storage for building systems [10].The impact of coupled heat and moisture transfer affects a building's energy consumption [11].Testing the optimal thermal performance of buildings requires us to model the energy transfer between a building and its surroundings.However, solving the problems of heat transfer in buildings is very challenging.Various heat exchange processes are possible between the indoor building environment and the external environment.The heat flow is conducted by the walls, roof, ceiling, floor, and so on.Heat transfer also takes place on different surfaces by convection and radiation.Heat is also added to space due to the presence of human occupants, and due to the use of lights and equipment [12].
Many building energy simulation programs (e.g., BLAST, DeST, DOE-2, ECOTECT, IES, EnergyPlus, TAS, TRNSYS, FLUENT, ESP-r [13], etc.) have been developed over the past few decades [14][15][16][17].However, the code for energy simulation in buildings often uses heat transfer models based on a one-dimensional heat flow, which can give unreliable results or can calculate only stationary processes [18].Therefore, the TRNSYS software (18, University of Wisconsin-Madison, Madison, WI, USA) is designed mainly to assess the performance of thermal and electrical energy systems under transient conditions.In addition, SolidWorks Flow Simulation software was tested as a possible tool for modelling the energy balance in building structures [19].Compared with this software, the current version of COMSOL Multiphysics (5.3, COMSOL, Inc., Stockholm, Sweden) includes modules which allow simultaneous assessment of the heat and mass transfer for both stationary and time-dependent processes.In this regard, the paper by [20] deals with a two-dimensional transient modelling of energy and mass transfer in the porous building components using COMSOL Multiphysics.The authors in [21] deal with the study of moisture hysteresis influence on mass transfer through bio-based building materials in dynamic state.In [22], COMSOL Multiphysics was used for testing of an impact of paper and wooden collections on humidity stability and energy consumption in museums and libraries humidity stability and energy consumption.Similarly, in [23], the multiphysics modelling of heat and moisture induced stress and strain of historic building materials and artifacts was tested.The paper by [24] presents a renewed approach to the development of two layers super-insulated precast concrete structural sandwich panels using the thermal path method.The heat transfer was experimentally measured and the results were verified using a 3D component of COMSOL Multiphysics software.The research by [25] aims at simulating the process of spatial heat transfer in a multilayer non-uniform structure of an external cast-in-place framed wall produced from polystyrene concrete with a stay-in-place formwork.Based on the physico-mathematical model developed with the use of ANSYS and COMSOL software complexes, parametric analysis of the impact of various factors on thermal behavior of the external wall was performed with the account of heat-stressed frame elements.Meanwhile, the authors in [26] present results of the investigation of evaporation cooling on drying capillary active building materials by numerical simulation and non-destructive measurements.The paper by [27] presents an investigation regarding the effect of indoor thermal to building's occupant satisfaction in a building that is mechanically ventilated.The study by [28] deals with the transient heat transfer in a multi-layered building wall through the facades of the buildings located in the city of Bechar in southwest Algeria.This physical model was presented to find the variation of the transient temperature in these structures and the heat flux through these elements, which depends on the air temperature of the inner surface and the instantaneous climatic conditions of the air outside.A combined finite element analysis and experimental validation approach to estimating effective edge conductivities of vacuum insulation panels embedded in foam-VIP composites is presented in [29].An integral component of numerical simulation using finite elements is the boundary conditions, which, in the case of simulating the thermal performance of a building element, are usually expressed in terms of the external surface temperature as a function of time.The purpose of this study is to examine the effect of the accuracy of the boundary conditions on the thermal performance simulation of building elements.
In study [30], the possibility of using COMSOL Multiphysics in simulations of the heat transfer in a building environment was verified by comparative Building Energy Simulation Test ("BESTEST") [31], which is a methodology that was developed by the International Energy Agency and that was used for systematically testing whole building energy simulation programs and diagnosing sources of predictive disagreement.An accuracy of the room model presented in this paper was verified by experimental testing on a real room with the same construction as the construction of a model formulated in the COMSOL Multiphysics interface.Comparison of the data obtained by computer simulation for this model with data obtained by experimental testing is presented in [32].
In addition to the computational kernel that enables complex multiphysical tasks, the current COMSOL Multiphysics version offers the Application Builder interface, which enables us to quickly create easy-to-use applications based on working models [33].From this aspect, we intend to use COMSOL Multiphysics to create a complex software tool for the assessment of energy performance of selected models of building parts.Simultaneously, we would like to ensure the possibility of universal use of this tool, not only for specialists who are trained to control the software but also for external users via a web interface [32].For this purpose, we are preparing a set of room modules in which users can specify both geometric and physical properties.
In this paper, we follow the paper [34], which presents the interim results of the research on the developing methodological procedure which could be used for assessment of the thermal stability of buildings with regards to their thermal accumulative parameters concerning compliance with valid European technical standards [35] used in building industry and architecture.We describe computer application programmed in the Application Builder of the COMSOL Multiphysics software, which can simulate cooling of a room as a 3D model under the selected winter condition.The first part of this paper includes the theoretical background of the studied problem.The second part describes a model of the studied room.The following parts demonstrate the use of the proposed computer application by testing the thermal stability of the room in case of its cooling concerning the outside wall's thermal insulation properties.
The aim of the simulations is to assess the temperature distribution of air inside the room and the heat losses in the winter season dependent on thermophysical properties of a thermal insulation of the outside wall.The Heat Transfer Module of COMSOL Multiphysics 5.3 [36] with the laminar flow was used to build a model of the tested room.This provides several tools for modelling coupled heat transport in building materials and adjacent air domains by taking into account heat storage and latent heat effects.Its computing core is based on numerical solving of the studied models by using the finite element method.

Model Description
The basic parameterised 3D model of the studied room was first formulated and tested by the standard Model Builder of COMSOL Multiphysics.Its geometric layout is shown in Figure 1.This pilot model includes only the basic elements that influence heat transfer between inside and outside environments.Therefore, it consists of side walls, floor, windows, doors and thermally insulated ceiling and outside wall with thermal insulation under the window.It also includes two heaters, which are needed to test the heat transfer in the case of heating the room.The model supposes that the tested room is surrounded by two neighboring rooms and by another one below.There is a corridor behind the wall with the door.The dimensions of all elements can be changed as required.Because of assumed time-dependent temperature distributions of all of the elements of the model, and also variable temperatures of air inside and outside the room, the studied problem has to be solved as non-stationary heat transfer using relation [37] to describe these processes: where symbol stands for the material density, c p is the specific heat capacity, T denotes the temperature, t is the time, v is the fluid velocity, λ is the thermal conductivity, and Φ denotes the inner heat-generation rate per unit volume.
The first part of the Equation (1) describes heat transfer by conduction in solids.The second part describes heat convection in fluids, and the third part describes heat accumulation in the mass of a specific domain.The sum of these processes is equal to the domain heat source.
It is assumed that the boundary conditions can be of the following types [37,38]: 1. Specified temperature (T s ) on a surface: T s = T (x, y, z, t) .2.
Heat flow density (q s ) on a surface: q s = q(x, y, z, t).

3.
Convective boundary condition: q x n x + q y n y + q z n z = h (T s − T e ) + q r .Symbol h stands for the heat transfer coefficient, T s denotes the surface temperature, T e is the convective exchange temperature, q r denotes the incident radiant heat flow per unit surface area, and q x , q y , q z are components of the heat flow density.
Finally, for the initial temperature field of a body at the time t = 0, it holds [37] where T 0 stands for the initial temperature of the body.

Simulation Conditions
Size and physical properties of the geometrical elements of the tested room are listed in Table 1. Figure 2 presents the temperature of air in the adjoining rooms and the temperature of the air outside the building.
Further conditions, geometry and material properties of the studied model include: • Room dimensions: length 8.7 m, width 5 m, height 3 m.

•
The initial temperature of the air inside the tested room 21 Table 1.Size and physical properties of geometrical elements of the studied room.

Geometrical Element Size [m]
Thermal Density Specific Emissivity Conductivity Heat Capacity Rear

Domain Setting
In the presented room model, the heat energy is transferred by conduction, convection and radiation.The conduction describes heat passage through the walls, thermal insulation, windows, door, and heaters of the room.The convection describes heat transfer in the internal room air.The sum of these processes is described by simplified Equation ( 3)

Boundary Setting
Convective heat flux on external boundaries and boundaries related to internal air is represented by Neumann boundary condition where T stands for the boundary temperature and T in f is the external temperature.
A condition that describes heat transfer by radiation is applied on boundaries with internal air where ε stands for the emissivity, σ is Stephan-Boltzmann constant, T is the boundary temperature, and T amb is the ambient temperature.A continuity boundary condition is used for the rest of the boundaries.The boundary conditions that can be defined by user for geometrical elements of the studied room model are described in Table 2.In the other parts of the model, it is assumed that the condition is not applicable or it can be neglected, or this condition in COMSOL Multiphysics is automatically inserted independently of the user based on other defined input parameters and conditions.Similarly, the continuity conditions are automatically defined.

Computer Model Construction
The basic model of the room was built in the standard Model Builder interface of COMSOL Multiphysics.It was then transformed into a comfortable computer application in the Application Builder interface, so that a wide range of users can perform the simulations, regardless of their knowledge of programming the models in COMSOL Multiphysics.The user interface of the programmed application is "interactive" and it contains the cards in which it is possible to change the input parameters of the model, view results of the simulations and export the output data (as shown in video presentation attached as the Supplement S1).
The source code of the application is available as the Supplement S5.Sequential sorting of cards allows the sequencing of modifying the model according to the user's requirements, its subsequent solution and the assessment of results of the simulations.This process includes the following steps: 1.
Defining the dimensions and location of all model elements.2.
Defining the physical properties (specific heat capacity, thermal conductivity, density) of all model elements.

3.
Defining the initial and boundary conditions.4.
Solving the model under modified conditions. 5.
Displaying the results as both 3D and 2D plots and export the output data into the TXT files.
This application allows the user to enter the dimensions of all of the model's elements, the material properties of all of the solids (thermal conductivity, density, and specific heat capacity).The properties of the air are automatically specified from Material Library of COMSOL Multiphysics.The user can also enter initial temperatures of all elements and heat transfer coefficients inside and outside the building.The temperature of the air outside the building was obtained from real data and inserted into the model in the Model Builder interface.It can be changed as needed.
The application can compute and display the time evolution of the air temperature for a selected place (point) inside the room.It also provides information about the temperature heat flow distribution in the room at the selected time of its cooling or heating.The output data can then be exported to a TXT file and further processed and evaluated by MATLAB (R2018a, The MathWorks, Inc., Natick, MA, USA) or another computer algebraic software program.

Simulation Setup
A free tetrahedral mesh specified by a matrix of about 4,200,000 elements (degrees of freedom) was used to solve the finite element problem in the model environment.This matrix includes information about the physics, material property and boundary conditions based on the number of elements and the discretisation order.In this sense, the size of this matrix is closely related to the complexity of the whole problem, demands for the size of computer memory, and time necessary for numerical modelling of the whole spatial model.The meshed model is shown in Figure 3a,b.A blue color indicates an air domain.The used mesh quality was influenced by computer capability.A simple test was performed to determine the air temperature deviation in the study room, dependent on the mesh quality, and to find the optimal mesh density subsequently.The results are shown in Figure 3c.The temperature deviation at the given time of cooling of the room was calculated with regards to the temperature of the model with a minimum element size of 1.49 mm and the maximum element size of 149 mm.These values correspond to the defined size of the model elements and the quality of the finest mesh that could be set with regards to the memory of the computer used.Figure 3c shows an increase of temperature deviation with the increasing maximum and minimum element size (i.e., with decreasing mesh density).For the most coarse meshed model, the maximum temperature deviation is 0.93%, and for the model with the finest mesh, the maximum temperature deviation is 0.57%.With regards to the accuracy of the calculation, the simulation time and the available memory of the used computer, the maximum element size was set to the size of 186 mm.Further choices included the selection of the minimum element size of 1.86 mm, the curvature factor of 0.2, and the maximum element growth rate of 1.3.
The average simulation time took about 110 min for the selected sampling period of 4 h.The model was calculated using COMSOL Multiphysics 5.3 on a computer with the processor Intel(R) Core(TM) i7-8550U CPU 1.80 GHz, 4 cores.The available memory was 16.22 GB.The whole simulation was performed by the PARADISO solver which is a part of the COMSOL Multiphysics software.
An absolute tolerance for control of the largest allowable absolute error at any step of simulation was 0.001.
The quadratic Lagrange discretisation method was used for solving of the studied model.Further discretisation methods tested during simulations included linear and quadratic serendipity, Lagrange cubic discretisation, cubic serendipity, Lagrange quartic, quartic serendipity, and Lagrange quintic discretisation.Best results were obtained by the Lagrange quadratic and quadratic serendipity discretisation methods.The size of computer memory and time necessary for numerical modelling were similar for both of these methods.Errors were recorded for linear discretisation under the conditions considered.Too extensive computing time or lack of computer memory were the most restrictive factors for other discretisation methods.

Results
The main results of simulations are shown in Figures 4-6. Figure 4 shows a decrease of the air temperature in the middle the room as a result of heat transfer to the unheated adjoining rooms and to the environment outside the building.Figure 5 depicts the temperature and heat flux curves through the outside wall after two days of cooling the room.These curves can also be obtained for any other time of cooling by entering the required value of the appropriate box of the computer application.The time evolution of temperature and heat flow provides important information about the wall's thermal accumulation properties and the time needed to achieve steady heat transfer.
Similarly, Figure 6a shows sections with the distribution of the air temperature inside the room as well as temperatures in the walls.Figure 6a shows initial temperature distribution and Figure 6b shows temperature distribution after one day of cooling.This data can also be obtained for any other time of cooling.
The output data from the simulations using COMSOL Multiphysics were processed with the use of MATLAB software, as shown below.In the presented case, the width of the room is 5 m and the thermal insulation thickness is 400 mm.The supposed temperature everywhere outside the room is 0 • C. The other parameters are the same as in the previous case.The output data from the simulations before processing in MATLAB software is included as Supplements S3, S4.   Figure 7 shows the temperature distributions in the selected section of the room processed in MATLAB as color slices with the desired color resolution.A video sequence of these slices calculated for the cooled room with time steps of four hours forms the temperature evolution, which is shown in Supplement S2.Figures 8 and 9 depict the temperature of the air inside the cooled room with respect to the specific heat capacity of the outside wall thermal insulation.Figure 8a depicts the time evolution of the air temperature in the middle of the room as a 3D plot.The 2D temperature distribution is shown in Figure 8b.The velocity of cooling the room decreases with increasing specific heat capacity of the thermal insulation of the outside wall, which is reflected by increasing temperature of the air inside the room.Figure 9 presents deviations of the air temperature while cooling the thermally insulated room from the air temperature in the room with non-insulated outside wall.It is evident that the temperature deviation increases as the thermal insulation capacity of the thermal insulation increases.Under these conditions, the maximum temperature deviation is 1.3 • C in the room thermally insulated by the material of the specific heat capacity of 6000 J The degree of utilisation of heat gains or thermal heat losses is directly related to the thermal inertia of the building.Based on the inner heat capacity of the building structure, a time to achieve desired temperature decrease can be determined.In the presented example, the desired temperature was computed from the polynomial, which can approximate temperature evolution of the air inside the tested room during its cooling.For this purpose, the output data from the simulations were approximated by polynomial of the 6th degree using the MATLAB software.The applied procedure for computing the described model and conditions is shown in Figure 10.Tables 3 and 4 show the impact of the specific heat capacity of the thermal insulation on the desired time of the temperature decrease.For testing, the specific heat capacity of the insulation was 500-6000 J • kg −1 • K −1 .The thermal insulation thickness was 200 mm, 300 mm, and 400 mm.The desired time computed for tested room with windows in the external wall (see Figure 1) was compared with a model of a room with an outside wall without windows, which is completely covered by thermal insulation.In the model of room of windows, an area of windows and the window frames represents 48 percent of the outside wall area.The area of the thermal insulation under the windows represents 35 percent of the outside wall area.The results for desired time of the 20 percent temperature decrease are shown in Table 3. Table 4 shows the results for desired time of the 50 percent temperature decrease.The results for room with windows are presented in the left-hand columns of Tables 3 and 4, the results for room without windows are shown in the right-hand columns of the tables.In both tables, the results confirm increasing time of the room temperature drop with increasing specific heat capacity and the thickness of the thermal insulation of the outside wall.It is also evident that the windows have a considerable influence on the heat loss of the room.The time to reach the temperature drops by 20 percent was in the room with windows from 1.4 to 1.8 times lower than in the room without windows.The time to reach the temperature drops by 50 percent was in the room with windows from 1.7 to 2.9 times lower than in the room without windows under the studied conditions.The mesh statistic for determination of desired temperature decrease with respect to the geometry of the model and the thickness of the thermal insulation are given in Table 5.

Discussion
Heat accumulation properties of building structures' materials significantly affect the cost of heating or cooling of buildings.In practice, the energy performance of buildings is assessed, in particular by theoretical calculations based on valid technical standards for construction and architecture.
However, these standards are often very demanding and, moreover, they consider many simplifications compared to real conditions, which can significantly affect the accuracy of the results obtained.In contrast, the use of suitable software tools for modelling physical processes is an important advantage, the possibility to perform simulations for marginal conditions inserted in the form of data obtained by real measurement on the object under consideration.However, for objective assessment of the thermal properties of buildings, the output data from computer simulations should be subsequently processed by the uniform methodical procedures of the technical calculations generally used in engineering practice.Therefore, research is focused on testing the factors that influence the thermal stability of a building.We have used a combination of computer simulations in the COMSOL Multiphysics user interface and theoretical calculations using the European technical standards for the building industry.
This paper presents a study of the thermal stability of the model with various geometry sizes and material properties of all elements of the room.Compared to the previously published models of a heated (cooled) room [32,34], the thermal insulation of the outer wall is now incorporated into the model, which allows its much wider use by testing of thermal distribution inside the studied room and searching for the optimal thermal insulation material properties.Additionally, the graphical user interface created with Application Builder provides engineers and scientists the opportunity to use the model without necessarily having knowledge of COMSOL Multiphysics control.Accuracy of the basic model that was used for transformation into the computer application was verified by experimental testing [32].
The current results are presented in Figures 7-10 and in Supplement S2.The mesh density analysis was performed by comparing the air temperature evolution in the model of the selected mesh quality and in the model with the maximum mesh density that could be set with regards to available memory of the used computer.The maximum temperature deviation calculated for the mesh of the presented model was 0.57%.The mean and median error values were 0.08% and 0.04%, respectively.

Conclusions
The present results form the basis for more complex simulations of thermal properties inside newly designed buildings.The heat transfer incorporated into the assembled model will be enriched by information related to moisture transfer and various types of heating as well.The future work will be devoted to the subsequent processing of output data obtained by computer simulations to determine the additional criteria and parameters that are used to assess the heat losses of buildings.This information will enable a more detailed balance of energy to optimize its requirements for building heating (cooling) or architectural changes that can eliminate thermal losses of buildings under given conditions.In this respect, the MATLAB computational tools will be used to analyse the temperature evolution in the room.The influence of permanent or cyclic heating of the room by modifying the boundary conditions will be studied as well.
It is assumed that the future work will include studies related to visualization in engineering [39], modelling of temperature distribution both inside the room and in the walls [30,32,34], and the study of comfortable indoor environment to keep emotional and physical states of individuals healthy [40].Simulations or thermal imaging [39] for breathing analysis using thermal and depth imaging camera video records will be included in the future research too.

Supplementary Materials:
The following are available online at https://www.mdpi.com/1996-1073/11/11/3205/s1, Video S1: Presentation of the computer application for simulation of cooling the room with thermally insulated outside wall.Video S2: Temperature evolution in the section of the cooled room.Data S3: Data for Figure 7-Distribution of the air temperature in selected times of cooling the room.Data S4: Data for Figures 8-10-Temperature of the air in the middle of the room concerning the specific heat capacity of the outside wall thermal insulation.The source code of the programmed software application S5 is available at https://drive.google.com/file/d/1riwSI3AmQ1jZhPGjT0af5l_nrhR5-IwS/view?usp=sharing.

Figure 1 .
Figure 1.Geometric layout of the studied model.

Figure 2 .
Figure 2. Temperatures of air in the adjoining rooms and the temperature of the air outside the building-data used as the boundary conditions for the tested model.

Figure 3 .
Figure 3. Meshed model: (a) sketch of the model; (b) close-up view of a section near the outer wall; and (c) percentage temperature deviations related to the mesh densities.

Figure 4 .
Figure 4. Time evolution of the air temperature in the center of the room: (a) sketch of the monitored location (the red point) inside the room; (b) dependence of the air temperature on time of cooling the room.

Figure 5 .
Figure 5. Temperature and heat flow in the outside wall after two days of cooling the room: (a) position of the monitored place (the red line); (b) temperature and heat flow distribution.

Figure 6 .
Figure 6.Slices of temperature inside the cooled room: (a) initial temperature distribution; (b) temperature distribution after one day of cooling.

Figure 7 .
Figure 7. Distribution of the air temperature in selected times of cooling the room: (a) the spatial location of the plane in which the temperature was determined; (b) slices with the temperature distributions.

Figure 8 .
Figure 8.Time evolution of the air temperature in the middle of the room concerning the specific heat capacity of the outside wall thermal insulation as (a) 3D plot and (b) 2D temperature distribution.

Figure 9 .Figure 10 .
Figure 9. Deviation of the air temperature while cooling the thermally insulated room from the room with non-insulated outside wall as (a) 3D plot and (b) 2D temperature distribution.
• C. • Temperature of air in the adjoining rooms of the building 0 • C. • Velocity of the airflow inside the tested room 0.03 m • s −1 .• Heat transfer coefficient inside the building 8 W

Table 2 .
Boundary conditions defined by user for geometrical elements of the studied room model.

Table 3 .
Computed time of the 20 percent temperature decrease.t D200 is the cooling time for a 200 mm thick wall insulation, t D300 for a 300 mm thick wall insulation, and t D400 for a 400 mm thick wall insulation.

Table 4 .
Computed time of the 50 percent temperature decrease.t D200 is the cooling time for a 200 mm thick wall insulation, t D300 for a 300 mm thick wall insulation, and t D400 for a 400 mm thick wall insulation.

Table 5 .
The mesh statistic for determination of desired temperature decrease for a 200 mm (I 200 ) thick wall insulation, for a 300 mm (I 300 ) thick wall insulation, and for a 200 mm (I 400 ) thick wall insulation.The numbers given indicate the number of elements.