Thermal Properties of Hemp Shives Used as Insulation Material in Construction Industry

: The article presents the results of studies concerning raw hemp shives obtained from the Polish crop of industrial hemp as a loose-ﬁll thermal insulation material. The study focuses mainly on the measurements of the pore size distribution, thermal conductivity and air permeability of material. An increase in the value of the thermal conductivity coefﬁcient (0.049–0.052 W/(m · K)) was demonstrated with an increase in the bulk density. The porosity of the individual pieces of shives is 78.7% and the predominant number of pores is in the diameter range of 0.9–3 µ m. The paper also presents an example of the use of the tested material as thermal insulation of the wooden frame wall. The heat ﬂow analysis was performed in various wall variants (insulation thickness: 100, 200 and 300 mm and pressure difference 0, 5, 10 and 15 Pa). A clear inﬂuence of the variables on the temperature distribution was observed.


Introduction
Hemp shives are the wooden parts of the industrial hemp stem.These are fiber plants grown for the needs of various industries, including textile, food, pharmaceutical and construction.In Poland, these plants are cultivated for about 4 months.During this time, depending on the variety and sowing density, they can grow to a height of up to 4 m.From one hectare of cultivation, it is possible to obtain 6-9 tons of dry matter of hemp straw [1,2], of which about 75% is the wooden core of the stem yielding the shives [3,4].Hemp needs large amounts of carbon dioxide to grow.One ton of dry matter in hemp yield can store 1000-2900 kg of CO 2 [1,[5][6][7].The straw as raw material has therefore a highly negative carbon footprint.
Construction is an energy-consuming branch of the economy.Energy consumption for construction processes and the use of buildings correspond to approximately 40% of total world energy usage [8].Insulating materials used today, such as rock wool and polystyrene, have a high carbon footprint.It is reasonable to search for alternative insulation materials with minimal embodied energy.The use of semi-finished products obtained from straw (fiber and shives) in the production of building materials can significantly reduce the level of their embodied energy [4].
In the past century, they were usually treated as waste when obtaining the fiber from the stem.As waste, they were of poor quality.They contained a lot of dust, residues on fibers, whereas fractions were very irregular (e.g., presence of pieces several centimeters Energies 2022, 15, 2461 2 of 18 long).Such shives are mainly used as bedding for animals, as they are highly hygroscopic and thus improve hygiene in animal husbandry rooms [1].
Hemp shives are currently used (since the 1980s [9]) in construction as a lightweight aggregate the in composites based on lime, clay or magnesia binder.For this purpose, high quality of shives is required.They should be cleaned of dust (from plant and soil) as this increases the need for water and binder; moreover, they should be cleaned of fiber clumps (a small number of single, short fibers are recommended because they act as microreinforcement).Both finer and coarser fractions are used, but the shives mix should have an appropriate graining curve.
The shives in a dry form, without strict observance of the specified fraction and fiber content, can also be used as loose insulation in partitions in a wooden frame structure (insulation of walls, floors, ceilings and roofs).The effectiveness of the insulation will depend on the degree of shives density, which was confirmed in the tests on insulation made of flax shives [10].The efficiency of loose insulation is both affected by thermal conductivity and prevention against wind washing or air filtration through the material.The air filtration can significantly reduce the thermal resistance of the thermal insulation layer [11].Typically, this type of construction requires a wind barrier.
High porosity, in the range of 57-92% according to the literature [12][13][14], is one of the most characteristic features of shives.According to Collet et al. [15], the structure of hemp shives consists of a network of closed pores with a diameter predominantly from 5 to 50 µm.The pore walls of the shives are cellular.The cells in these walls are approximately 2.5 µm in diameter [15].Various methods are used to determine the porosity of the shives.Jiang et al. [16] examined the porosity using the mercury porosimetry intrusion method, obtaining a result of 76.67 ± 2.03%.It has also been found that testing with this method does not change the pore structure.The same shives were also tested for porosity by CT scanning (computed tomography), and a lower porosity of 50% was obtained.The reason was, among others, that this method only takes into account the pores larger than 2.54 µm.However, it gives more information about the shape of the pores and their interconnections in three dimensions.The pycnometric method was also employed to study the internal porosity, using toulene as a filling substance [17].Depending on the method, the porosity of shives can be tested on selected single shives particles or on a compacted mixture of shives.In the latter case, the pore volume between the shives is also taken into account.However, the literature does not always provide the information on the preparation of a sample for porosity testing.The total porosity and pore size distribution depend on the variety of hemp from which the shives were obtained, as the stalks differ in terms of parameters.The structure may also be influenced by the method of retting and the treatment of the hemp stem to obtain shives.Therefore, it is important to measure the parameters of the shives that are going to be used in practice.
High porosity influences of the low thermal conductivity, and thus makes the material applicable in building constructions.According to the literature [12,[18][19][20][21][22][23], the thermal conductivity coefficient of shives is in the range of 0.048-0.058W/(m•K).This parameter is usually tested on a sample of shives thickened to a different degree using a plate apparatus [12,18,19].When examining the material in this way, the thermal conductivity of the space between the shives is also taken into account.Tests are usually carried out on the shives dried to a constant weight.Balciunas et al. [12] proved that the shives seasoned in air-dry conditions, under the influence of air humidity, have a thermal conductivity 8.2-8.8% higher than the shives dried to constant weight.The direction of the capillary pores in relation to the heat flow stream also affects the thermal conductivity.The example of a hemp-lime composite showed that when the shives were arranged with fibers (and, also, capillaries) perpendicular to the heat flow, the composite was characterized by a lower thermal conductivity [24].
High porosity is responsible for the ability to absorb water.Arnaud et al. [25] proved that the shives are able to absorb significant amounts of water (2-3 times more in relation to the dry weight).After 10 min from immersion, their saturation reached 95%.Jiang et al. [16] showed that the capillaries extend along the hemp stem (along the direction of the fibers), but the walls between the capillaries contain a perforation that allows moisture to pass between the capillaries.On the example of the study of the hemp-lime composite, it was found that the shives in which the capillaries will be located perpendicular to the water table show greater water uptake capacity.Placing the shives with capillaries parallel to the water table allowed for the reduction of water absorption by about 18%, despite the presence of the aforementioned perforation in the capillary walls [24].
The study investigated the physical properties of hemp shives from the crops of industrial hemp grown in Poland.They were tested for use as a loose thermal insulation material.The thermal conductivity, porosity and air permeability of shives with varying degrees of density were investigated.There are publications on the thermal behavior of hemp shives-based building materials [26][27][28][29][30][31][32], but it is difficult to find the analyses of shives properties such as correlation between density and thermal conductivity; in addition, there are no results on air permeability.Another novelty in the work will correspond to characterization of the pore size distribution of shives.The investigated thermal parameters were used to simulate a construction (various variants) filled with shives in the Delphin software (http://bauklimatik-dresden.deaccessed on 11 February 2022).

Possible Applications of Shives in Construction Industry
In practice, hemp-lime composites are the main application of shives as lightweight aggregates.This composite may be used as the filling of the timber frame structure in the walls, floors and roof.It has a thermal conductivity coefficient in the range of 0.074-0.138W/(m•K) [26][27][28] depending on the composition used, e.g., shives to lime ratio.The composite also has a high thermal capacity of the order of about 1000-1600 J/(kg•K) [29][30][31][32] which positively influences the preservation of heat in the walls, affecting the thermal comfort in the rooms.
Hemp aggregate is also used in the composites based on the magnesia binder.These composites, due to their greater strength than those based on lime, are used in the production of prefabricated elements.Their thermal conductivity also depends on the composition and can be, for example, 0.074-0.082W/(m•K) [18,28].
Bricks consisting of clay and shives are also produced and researched scientifically.Depending on the density, their thermal conductivity can be in the range of about 0.09-0.18W/(m•K) [33].
Shives with finer fractions, often in combination with hemp fibers, are used to make thermal-insulating mortars and plasters [34][35][36].These additives also improve the resistance to shrinkage.
Hemp shives, such as flax shives [10], are used in their raw form as loose thermal insulation of partitions in a timber frame structure.In the past, shives dry mixed with lime or clay was used as insulation for ceilings.

Hemp Shives
The hemp shives of the Białobrzeskie variety of the Polish producer Podlaskie Konopie (Białystok, Poland) were used for the research.They were obtained with the help of a small decorative line, but adapted to the production of high-quality hemp shives for construction purposes.They are characterized by a low content of dust and fibers.The shives used in the tests are shown in Figure 1a, while the shives length distribution is shown in Figure 1b.To determine the mass share of different ranges of shives lengths, 180 pieces of shives randomly selected from the mixture provided by the manufacturer were measured.The literature describes the studies in which shives of similar size were used [35].The hemp shives of fractions 0-20 mm in length and a 4.4% proportion of length >20 mm in relation to the weight were tested.The largest share was constituted by the shives with a length of 5-10 mm.

Pore Size Distribution
The mercury porosimetry method was used to study the microstructure of hemp shives.Many authors indicate that the method of measuring porosity using mercury porosimetry is the right solution for the study of plant fibrous materials [16,[37][38][39].The porosity parameters were determined using the equation proposed by Washburn, which shows the relationship between the diameter of the pores and the pressure of mercury filling them (1): where R is radius of cylindrical pores (μm), F is shape factor (for cylindrical pore calculations equal to 2), γ is mercury surface tension (N/m), cosθ is the limit value of the contact angle and p is the mercury pressure (MPa).The test is applicable to the materials with open porosity, incompressible with cylindrical pore characteristics [40].The parameters determined during the study were: total porosity, total pore area, median pore diameter, bulk density and apparent density.An Autopore Mercury Porosimeter (Atlanta, GA, USA) was used to investigate the pores in the range from 300 μm to 6 nm during two measuring cycles-in a low pressure station (measuring range 0.014 to 0.21 MPa) and in a high pressure station (from 0.21 to 227.53 MPa).The samples of the woody parts of hemp shives were prepared for the research.The initial selection of shives type was made in order to the eliminate structures that could be damaged during the test.Before testing, the tested material was dried to a constant weight at a temperature of 60°C.The measurements were performed using a measuring vessel-designed for solid materials, with a capacity of 3 cm 3 .The weight of the samples was approximately 0.1872 g ± 0.016 g.Three pieces of shives of different sizes were randomly selected for the porosity test.They are designated as HS1, HS2 and HS3.

Thermal Conductivity
The study was carried out in order to determine the value of thermal conductivity of loose hemp shives as a function of density.The thermal conductivity study was implemented on the basis of international standard EN 12667 [41].The study was set up in the Laser Comp Fox 602 plate apparatus (TA Instruments, New Castle, DE, USA), the The literature describes the studies in which shives of similar size were used [35].The hemp shives of fractions 0-20 mm in length and a 4.4% proportion of length >20 mm in relation to the weight were tested.The largest share was constituted by the shives with a length of 5-10 mm.

Pore Size Distribution
The mercury porosimetry method was used to study the microstructure of hemp shives.Many authors indicate that the method of measuring porosity using mercury porosimetry is the right solution for the study of plant fibrous materials [16,[37][38][39].The porosity parameters were determined using the equation proposed by Washburn, which shows the relationship between the diameter of the pores and the pressure of mercury filling them (1): where R is radius of cylindrical pores (µm), F is shape factor (for cylindrical pore calculations equal to 2), γ is mercury surface tension (N/m), cosθ is the limit value of the contact angle and p is the mercury pressure (MPa).The test is applicable to the materials with open porosity, incompressible with cylindrical pore characteristics [40].The parameters determined during the study were: total porosity, total pore area, median pore diameter, bulk density and apparent density.An Autopore Mercury Porosimeter (Atlanta, GA, USA) was used to investigate the pores in the range from 300 µm to 6 nm during two measuring cycles-in a low pressure station (measuring range 0.014 to 0.21 MPa) and in a high pressure station (from 0.21 to 227.53 MPa).The samples of the woody parts of hemp shives were prepared for the research.The initial selection of shives type was made in order to the eliminate structures that could be damaged during the test.Before testing, the tested material was dried to a constant weight at a temperature of 60 • C. The measurements were performed using a measuring vessel-designed for solid materials, with a capacity of 3 cm 3 .The weight of the samples was approximately 0.1872 g ± 0.016 g.Three pieces of shives of different sizes were randomly selected for the porosity test.They are designated as HS1, HS2 and HS3.

Thermal Conductivity
The study was carried out in order to determine the value of thermal conductivity of loose hemp shives as a function of density.The thermal conductivity study was implemented on the basis of international standard EN 12667 [41].The study was set up in the Laser Comp Fox 602 plate apparatus (TA Instruments, New Castle, DE, USA), the heat flow meter method.
The measuring space of the apparatus is defined by dimensions: 600 mm × 600 mm × 190 mm.The lower plate is stationary, while the upper plate is movable up and down.The highly sensitive digital sensor converters (heat meters) are mounted on the surfaces of both plates.The Fox 602 instrument operates on the basis of one-dimensional Fourier law (2): where q is the heat flux density flowing through the sample (W/m 2 ), λ is the thermal conductivity coefficient of the sample (W/(m•K)) and dT/dx is the temperature gradient on the isothermal flat surface of the sample (K/m).The Fox 602 apparatus was calibrated before the measurements.The calibration sample (mineral wool) was certified by the Joint Research Centre, Institute for Reference Materials and Measurements in Geel, Belgium.According to the producer, the measuring uncertainty equals 1%.
Testing the samples involved placing them in a stabilizing frame.Every prepared sample (10 in total) was created by random selection of conditioned loose hemp shives.It was intentional to prepare the samples differing with density (from 109 to 124 kg/m 3 ).Every sample was tested 13 times.The loose material was placed evenly in the frame and subjected to initial manual compaction.The dimensions of the frame made of extruded polystyrene were as follows: 600 mm × 600 mm (external), and 520 mm × 520 mm (internal) as well as the height of 45 mm (Figure 2).Hemp shives are rigid enough to maintain the shape which was formed at the beginning of the measurement.The measuring space of the apparatus is defined by dimensions: 600 mm × 600 mm × 190 mm.The lower plate is stationary, while the upper plate is movable up and down.The highly sensitive digital sensor converters (heat meters) are mounted on the surfaces of both plates.The Fox 602 instrument operates on the basis of one-dimensional Fourier law (2): where q is the heat flux density flowing through the sample (W/m 2 ), λ is the thermal conductivity coefficient of the sample (W/(m•K)) and dT/dx is the temperature gradient on the isothermal flat surface of the sample (K/m).
The Fox 602 apparatus was calibrated before the measurements.The calibration sample (mineral wool) was certified by the Joint Research Centre, Institute for Reference Materials and Measurements in Geel, Belgium.According to the producer, the measuring uncertainty equals 1%.
Testing the samples involved placing them in a stabilizing frame.Every prepared sample (10 in total) was created by random selection of conditioned loose hemp shives.It was intentional to prepare the samples differing with density (from 109 to 124 kg/m 3 ).Every sample was tested 13 times.The loose material was placed evenly in the frame and subjected to initial manual compaction.The dimensions of the frame made of extruded polystyrene were as follows: 600 mm × 600 mm (external), and 520 mm × 520 mm (internal) as well as the height of 45 mm (Figure 2).Hemp shives are rigid enough to maintain the shape which was formed at the beginning of the measurement.
In order to avoid condensation in the tested material, the shives samples were separated from the lower apparatus plate with a thin foil.Then, the samples were compacted by the FOX 602.Due to the hardness of dry shives, it was only possible to achieve the density of material in the range 109-124 kg/m 3 .The measurements of the thermal conductivity coefficient were made at the average temperature of 10 °C, with 0 °C on the bottom plate and 20 °C on the top plate.The measurements were conducted automatically using the WinTherm32 software (http://lasercomp.comaccessed on 11 February 2022).For each of the 10 prepared samples, 13 readouts of thermal conductivity coefficient were obtained.

Air Permeability
The study was carried out in order to determine the value of air permeability κ (m 2 ) of loose hemp shives in the function of density.The values of κ were determined on the basis of the Forchheimer equation, which is an extension of Darcy's law, describing the effect of a square resistance on the pressure drop across the porous medium.For flows characterized with Reynolds number greater than one, derogation from the Darcy law is In order to avoid condensation in the tested material, the shives samples were separated from the lower apparatus plate with a thin foil.
Then, the samples were compacted by the FOX 602.Due to the hardness of dry shives, it was only possible to achieve the density of material in the range 109-124 kg/m 3 .The measurements of the thermal conductivity coefficient were made at the average temperature of 10 • C, with 0 • C on the bottom plate and 20 • C on the top plate.The measurements were conducted automatically using the WinTherm32 software (http://lasercomp.comaccessed on 11 February 2022).For each of the 10 prepared samples, 13 readouts of thermal conductivity coefficient were obtained.

Air Permeability
The study was carried out in order to determine the value of air permeability κ (m 2 ) of loose hemp shives in the function of density.The values of κ were determined on the basis of the Forchheimer equation, which is an extension of Darcy's law, describing the effect of a square resistance on the pressure drop across the porous medium.For flows characterized with Reynolds number greater than one, derogation from the Darcy law is experienced due where dp is the pressure difference between pressure taps (Pa), dx is the distance between pressure taps (m), µ is the dynamic viscosity (Pa•s), κ is the air permeability coefficient (m 2 ), β is the Forchheimer coefficient (1/m), ρ is the fluid density (kg/m 3 ) and v is the velocity of fluid flow (m/s).
The measuring system (Figure 3) was constructed based on previous investigations of many researchers, which is described in detail in the work [42].
compressor, equipped with the filter and dryer, and delivered to the apparatus by a pressure reducer.
Testing the samples involved placing them in a central part of the system.Every prepared sample was created by random selection of conditioned loose hemp shives.The authors prepared the samples differing with density, from 90 to 120 kg/m 3 with 5 kg/m 3 resolution.During the measurements, the pressure difference between two taps, located in central part of the apparatus, was registered.Simultaneously, air velocity on the exhaust was registered.In order to avoid discrepancy of the density and viscosity of the air during investigation, a constant temperature +10 °C of compressed air was maintained.Every single calculated value of permeability was determined from the approximation of (3) on the basis of at least 800 pairs of air velocities and air pressures.
Due to the technical declaration of producers, the standard measurement uncertainty of the pressure difference is ±2 Pa in the measuring range, while the uncertainty of air velocity is ±0.04 m/s.Pressure uncertainty can be practically neglected in the achieved range (up to 6000 Pa) of the pressure difference, but the uncertainty of air velocity is very important.The authors estimated this influence: underestimation of air velocity can result in up to 29% overestimation of air permeability, while overestimation of air velocity can lead up to 55% of underestimation of air permeability.

Pore Size Distribution
After conducting the tests of open porosity on three samples of hemp material, the curves of the cumulative intrusion volume were obtained, presented in Figure 4.The characteristics of the pore distribution curves in all tested samples were similar.The Testing the samples involved placing them in a central part of the system.Every prepared sample was created by random selection of conditioned loose hemp shives.The authors prepared the samples differing with density, from 90 to 120 kg/m 3 with 5 kg/m 3 resolution.During the measurements, the pressure difference between two taps, located in central part of the apparatus, was registered.Simultaneously, air velocity on the exhaust was registered.In order to avoid discrepancy of the density and viscosity of the air during investigation, a constant temperature +10 • C of compressed air was maintained.Every single calculated value of permeability was determined from the approximation of (3) on the basis of at least 800 pairs of air velocities and air pressures.
Due to the technical declaration of producers, the standard measurement uncertainty of the pressure difference is ±2 Pa in the measuring range, while the uncertainty of air velocity is ±0.04 m/s.Pressure uncertainty can be practically neglected in the achieved range (up to 6000 Pa) of the pressure difference, but the uncertainty of air velocity is very important.The authors estimated this influence: underestimation of air velocity can result in up to 29% overestimation of air permeability, while overestimation of air velocity can lead up to 55% of underestimation of air permeability.

Pore Size Distribution
After conducting the tests of open porosity on three samples of hemp material, the curves of the cumulative intrusion volume were obtained, presented in Figure 4.The characteristics of the pore distribution curves in all tested samples were similar.The highest total porosity was found in the HS1 sample, while the smallest in the HS2 sample.From the cumulative intrusion vs. pore size figure, it can be clearly seen that the material pore dimensions range between 50 and 0.08 µm.
Energies 2022, 15, 2461 7 of 19 highest total porosity was found in the HS1 sample, while the smallest in the HS2 sample.
From the cumulative intrusion vs. pore size figure, it can be clearly seen that the material pore dimensions range between 50 and 0.08 μm.On the basis of the analysis of the obtained results (Figure 5.) the domination of pores in the range of 3-0.9 μm is visible, their share is about 40% of the total porosity.An additional peak is seen in the range of 50-10 μm.The pores in the range of 0.9-0.1 μm account for about 30% of the total porosity.The smallest pores in the range from 0.08 μm constitute a small share.The additional parameters determined by the measuring device are summarized in Table 1.The total porosity is 78.70% ± 1.65%.The bulk density is about 0.25 g/cm 3 .The median pore diameter is approx.1754.2 nm ± 280 nm.Bourdot et al. [13] showed higher values of porosity of hemp shives and the influence of shives length on porosity.The 0-5 mm shives were characterized by the porosity of 89.3%, while the shives of the 5-20 mm length were characterized by the porosity equal to 91.3%.Rahim et al. [14] also showed a higher porosity of shives, equal to 90.1%.However, this study does not provide a  On the basis of the analysis of the obtained results (Figure 5.) the domination of pores in the range of 3-0.9 µm is visible, their share is about 40% of the total porosity.An additional peak is seen in the range of 50-10 µm.The pores in the range of 0.9-0.1 µm account for about 30% of the total porosity.The smallest pores in the range from 0.08 µm constitute a small share.
highest total porosity was found in the HS1 sample, while the smallest in the HS2 sample.From the cumulative intrusion vs. pore size figure, it can be clearly seen that the material pore dimensions range between 50 and 0.08 μm.On the basis of the analysis of the obtained results (Figure 5.) the domination of pores in the range of 3-0.9 μm is visible, their share is about 40% of the total porosity.An additional peak is seen in the range of 50-10 μm.The pores in the range of 0.9-0.1 μm account for about 30% of the total porosity.The smallest pores in the range from 0.08 μm constitute a small share.The additional parameters determined by the measuring device are summarized in Table 1.The total porosity is 78.70% ± 1.65%.The bulk density is about 0.25 g/cm 3 .The median pore diameter is approx.1754.2 nm ± 280 nm.Bourdot et al. [13] showed higher values of porosity of hemp shives and the influence of shives length on porosity.The 0-5 mm shives were characterized by the porosity of 89.3%, while the shives of the 5-20 mm length were characterized by the porosity equal to 91.3%.Rahim et al. [14] also showed a higher porosity of shives, equal to 90.1%.However, this study does not provide a The additional parameters determined by the measuring device are summarized in Table 1.The total porosity is 78.70% ± 1.65%.The bulk density is about 0.25 g/cm 3 .The median pore diameter is approx.1754.2 nm ± 280 nm.Bourdot et al. [13] showed higher values of porosity of hemp shives and the influence of shives length on porosity.The 0-5 mm shives were characterized by the porosity of 89.3%, while the shives of the 5-20 mm length were characterized by the porosity equal to 91.3%.Rahim et al. [14] also showed a higher porosity of shives, equal to 90.1%.However, this study does not provide a methodology for testing porosity.It is possible that these results relate to the porosity of the compacted shives sample and that the volume of the pores between the shives is also taken into account.In the authors' own research, the results concern only shives (three randomly selected shives: HS1, HS2 and HS3).

Thermal Conductivity
Figure 6 presents the results of measurements of the mean thermal conductivity coefficient of 10 samples with different bulk densities.The range of achieved density was: 109-124 kg/m 3 .The obtained mean thermal conductivity coefficients of tested samples were in the range 0.049-0.052W/(M•K) with a little tendency to rise with density.This tendency may be due to reducing the pore size between separate shives.Along with reducing the pore size between the shives, the effect of heat transport by convection is reduced [43].
Energies 2022, 15, 2461 8 of 19 methodology for testing porosity.It is possible that these results relate to the porosity of the compacted shives sample and that the volume of the pores between the shives is also taken into account.In the authors' own research, the results concern only shives (three randomly selected shives: HS1, HS2 and HS3).

Thermal Conductivity
Figure 6 presents the results of measurements of the mean thermal conductivity coefficient of 10 samples with different bulk densities.The range of achieved density was: 109-124 kg/m 3 .The obtained mean thermal conductivity coefficients of tested samples were in the range 0.049-0.052W/(M•K) with a little tendency to rise with density.This tendency may be due to reducing the pore size between separate shives.Along with reducing the pore size between the shives, the effect of heat transport by convection is reduced [43].Delhomme et al. [44] tested the thermal conductivity of the shives with a bulk density of 106 kg/m 3 and 138 kg/m 3 , obtaining the results of 0.079 W/(m•K) and 0.1 W/(m•K), respectively.However, these were the shives with a smaller fraction (over 80% of the shives mixture mass passed through a sieve with a mesh of 4.75 mm).Balcianus et al. [12], on the other hand, tested the mixtures of shives with various ranges of 2.5-5 mm, 5-10 mm and 10-20 mm.The thermal conductivity of the shives in a dry state was 0.053 W/(m•K), 0.055 W/(m•K) and 0.059 W/(m•K), respectively.Other papers also confirm that the results of thermal conductivity are quite varied.Different types of hemp, different structure, pore distribution and fraction make it possible to obtain significantly different results of thermal conductivity by examining two different mixtures of shives with the same bulk density (Table 2).The obtained own results of the thermal conductivity Delhomme et al. [44] tested the thermal conductivity of the shives with a bulk density of 106 kg/m 3 and 138 kg/m 3 , obtaining the results of 0.079 W/(m•K) and 0.1 W/(m•K), respectively.However, these were the shives with a smaller fraction (over 80% of the shives mixture mass passed through a sieve with a mesh of 4.75 mm).Balcianus et al. [12], on the other hand, tested the mixtures of shives with various ranges of 2.5-5 mm, 5-10 mm and 10-20 mm.The thermal conductivity of the shives in a dry state was 0.053 W/(m•K), 0.055 W/(m•K) and 0.059 W/(m•K), respectively.Other papers also confirm that the results of thermal conductivity are quite varied.Different types of hemp, different structure, pore distribution and fraction make it possible to obtain significantly different results of thermal conductivity by examining two different mixtures of shives with the same bulk density (Table 2).The obtained own results of the thermal conductivity coefficient are similar to those presented in Table 2; however, they provide a broader view of the variability of this parameter depending on the bulk density of the one type of shives.Unfortunately, the cited publications did not investigate the pore size distribution.Brzyski et al. [10] studied the dependence of the thermal conductivity of flax shives on the degree of compaction.The shives with a density of 131-159 kg/m 3 were characterized by a thermal conductivity coefficient in the range of 0.047-0.053W/(m•K).No clear relationship was found.
The discrepancy in the results presented in this section is related to different varieties and fractions of shives, but also to their anisotropic properties.It is a material with directed fibers (the direction of the fibers is in the direction of the growth of the hemp stalk (Figure 7)).A sample of shives subjected to the thermal conductivity test contains random shanks.In an article by Brzyski et al. [24] it was proven that the composite based on limestone binder and hemp shives, compacted parallel to the heat flux, was characterized by a thermal conductivity coefficient lower by 12-16% (depending on the fraction of shives) compared to the composite compacted perpendicular to the direction of heat flow.Unfortunately, the cited publications did not investigate the pore size dist Brzyski et al. [10] studied the dependence of the thermal conductivity of flax shiv degree of compaction.The shives with a density of 131-159 kg/m 3 were characte a thermal conductivity coefficient in the range of 0.047-0.053W/(m•K).N relationship was found.
The discrepancy in the results presented in this section is related to different and fractions of shives, but also to their anisotropic properties.It is a mater directed fibers (the direction of the fibers is in the direction of the growth of th stalk (Figure 7)).A sample of shives subjected to the thermal conductivity test random shanks.In an article by Brzyski et al. [24] it was proven that the composi on limestone binder and hemp shives, compacted parallel to the heat fl characterized by a thermal conductivity coefficient lower by 12-16% (dependin fraction of shives) compared to the composite compacted perpendicular to the d of heat flow.

Air Permeability
Registered pressure drops and air velocity corresponding to them were use approximation of the Forchheimer equation (Equation (3)) to determine permeability coefficient κ of hemp shives.Different densities in the range of 90 ÷ 1 were investigated.The samples with higher density showed better resistanc filtration.It is well seen by increasing the pressure difference with increasing d the material.The results of the air permeability for investigated densities of hem are shown in Table 3.

Air Permeability
Registered pressure drops and air velocity corresponding to them were used in the approximation of the Forchheimer equation (Equation (3)) to determine the air permeability coefficient κ of hemp shives.Different densities in the range of 90 ÷ 120 kg/m 3 were investigated.The samples with higher density showed better resistance to air filtration.It is well seen by increasing the pressure difference with increasing density of the material.The results of the air permeability for investigated densities of hemp fibers are shown in Table 3.Air permeability coefficient for the hemp shives sample with the lowest density is 4.259744 × 10 −8 m 2 and for the sample with the highest density is 1.832226 × 10 −9 m 2 .The air permeability coefficient values decrease with increasing density of material.The same relationship, except for selected single measurements, was observed in other studies [45][46][47].Brzyski et al. investigated the air permeability of other insulation materials in the form of cellulose in the 30-70 kg/m 3 density range.The air permeability coefficient was from 2.209767 × 10 −7 m 2 to 2.896716 × 10 −10 m 2 , respectively [45].In another work [46], insulation from loose, unprocessed hemp fibers in the density range of 25-70 kg/m 3 was investigated.The obtained results are closer to those of hemp shives, between 9.775609 × 10 −8 m 2 and 1.053169 × 10 −8 m 2 .

An Impact of Convection on Heat Transfer through Hemp Shives Insulation
Convection plays a significant role in heat transfer in the porous structures.Mostly, this effect is considered in thermal insulations, as they have the highest thermal resistance within the whole building element.The convection effect is smaller if the pores between the fibers are closed due to material compacting.While it is possible to easily compact straw materials [48] it is not in the case of hard shives.In the presented laboratory measurements, the authors did not manage to compact hemp shives more into the 124 kg/m 3 density.This is why the convection analysis is important in this work.
The effect of convection on heat transfer in porous structures is described with the Rayleigh (Ra) number.The Rayleigh number relating to fibrous materials can be described as (4): where β is the thermal expansion coefficient (1/K), g is the gravity acceleration (m/s 2 ), δ is the insulation thickness (m), κ is the air permeability coefficient (m 2 ), ρ is the dry air density (kg/m 3 ), c p is the specific heat of dry air (J/(kg•K)), ν is the kinematic viscosity (m 2 /s), λ is the thermal conductivity coefficient of dry air (W/(m•K)) and ∆T is the temperature difference on both sides of insulation (K).
It is assumed that the heat transfer through the structures is based mainly on conduction if the Ra number is below a certain critical value.If the Rayleigh number is above the critical value, then the heat flow is mainly by convection.
According to ISO 10456 [49], it is recommended to include convection in heat transfer calculations for Ra number higher than 2.5 in the horizontal direction, in the vertical direction downwards Ra > 30, and in the vertical direction upwards Ra > 15.
On the basis of Formula (4), it could be concluded that convection takes part in heat exchange only in porous insulations with high air permeability, in thick layers and additionally with a high temperature difference.In fact, the contribution of forced convection in the heat transfer process in fibrous thermal insulation depends solely on the air permeability of the material and the exposure of the material to the air stream.
On the other hand, due to the increasing thermal requirements for building elements, thus, decreasing permissible heat transfer coefficients, the thermal insulation will be sig-nificantly thicker.It will additionally increase the share of possible convection in the heat transfer.For this reason, it is worth analyzing the Ra number in the designed buildings insulated with air-permeable materials.
Figure 8 presents the calculated values of Ra number for different densities of thermal insulation made of loose hemp shives for 3 different thickness of the insulation: 100, 200 and 300 mm.The analysis of the values shows that as the density increases, the Rayleigh number and thus the share of convection in heat transfer in the insulation decreases.The maximal calculated Ra numbers for 90 kg/m 3 density of hemp shives are: 8.94 for 100 mm, 17.88 for 200 mm and 26.82 for 300 mm, but minimal calculated for 120 kg/m 3 are in the range 0.385-1.154.It shows how large is the impact of density on the possible air filtration in loose material.It should be noted that even if the Ra number is smaller than the critical values pointed in the ISO 10456 [49] standard, air filtration may occur in the insulation with any pressure difference.The detailed analysis is presented in the next section.
convection in the heat transfer process in fibrous thermal insulation depends solely on the air permeability of the material and the exposure of the material to the air stream.
On the other hand, due to the increasing thermal requirements for building elements, thus, decreasing permissible heat transfer coefficients, the thermal insulation will be significantly thicker.It will additionally increase the share of possible convection in the heat transfer.For this reason, it is worth analyzing the Ra number in the designed buildings insulated with air-permeable materials.
Figure 8 presents the calculated values of Ra number for different densities of thermal insulation made of loose hemp shives for 3 different thickness of the insulation: 100, 200 and 300 mm.The analysis of the values shows that as the density increases, the Rayleigh number and thus the share of convection in heat transfer in the insulation decreases.The maximal calculated Ra numbers for 90 kg/m 3 density of hemp shives are: 8.94 for 100 mm, 17.88 for 200 mm and 26.82 for 300 mm, but minimal calculated for 120 kg/m 3 are in the range 0.385-1.154.It shows how large is the impact of density on the possible air filtration in loose material.It should be noted that even if the Ra number is smaller than the critical values pointed in the ISO 10456 [49] standard, air filtration may occur in the insulation with any pressure difference.The detailed analysis is presented in the next section.

Preprocessing
The thermal simulations were performed to present the hemp shives properties as thermal insulation in wooden frame construction.The simulations were performed using Control Volume Method in the Delphin 6.1 software, under steady state conditions in two cases: no air filtration and air filtration through the entire element.
Physical data for the simulation were adopted on the basis of laboratory measurements, presented in the article.The physical properties of materials, laboratory investigated are presented in Table 4.

Thermal Simulation 5.1. Preprocessing
The thermal simulations were performed to present the hemp shives properties as thermal insulation in wooden frame construction.The simulations were performed using Control Volume Method in the Delphin 6.1 software, under steady state conditions in two cases: no air filtration and air filtration through the entire element.
Physical data for the simulation were adopted on the basis of laboratory measurements, presented in the article.The physical properties of materials, laboratory investigated are presented in Table 4.The simulation model (Figure 9) corresponds to a wood frame wall filled with loose hemp fibers.Two examples differing in shives density were investigated: 110 and 120 kg/m 3 .
The simulation model (Figure 9) corresponds to a wood frame wall filled with loose hemp fibers.Two examples differing in shives density were investigated: 110 and 120 kg/m 3 .
For each density, three cases were modelled: wall with insulation thickness of 0.1, 0.2 and 0.3 m.For each thickness 4 conditions were applied: steady state conditions without air filtration and three cases of air filtration through the entire wall with pressure differences: 5, 10 and 15 Pa.The model consists of three layers: OSB from outside (0.02 m), loose insulation made of hemp shives (0.1, 0.2 and 0.3 m) and gypsum board from inside (0.012 m).The models height is 1.0 m.Automatic discretization of models with stretch factor 1.3 divided them into 3232 elements for the model with a thermal insulation thickness of 0.1 m as well as 0.2 m, and 3542 for the model with a thickness of 0.3 m.
Only temperatures and pressure difference were set in the model.The external temperature corresponds to the computational temperature of Warsaw: te = −20.00°C, internal as for living room ti = 20.00°C, while pressure difference Δp = 0.00, 5.00, 10.00 and 15.00 Pa correspond to possible conditions during wintertime.The pressure difference is subjected to 2 mm leakages situated 100 mm from the bottom and the top of the model.

Model Description
The energy balance is given with the Equation ( 4) is described in detail in [46].

( ) (
) where ρU is the internal energy density (J/m 3 ), λ is the heat conduction coefficient (W/(m•K), ∇T is the temperature gradient (K/m), κ is the intrinsic permeability of a porous medium (m 2 ), cp is the specific heat of dry air (J/(kg•K)), ΔT is the temperature difference (K), ν is the kinematic viscosity of air (m 2 /s), ∇p is the pressure gradient (Pa/m), ρg is the dry air density (kg/m 3 ), g is the gravity acceleration (m/s 2 ) and ρr is the internal heat source density(W/m 3 ).For each density, three cases were modelled: wall with insulation thickness of 0.1, 0.2 and 0.3 m.For each thickness 4 conditions were applied: steady state conditions without air filtration and three cases of air filtration through the entire wall with pressure differences: 5, 10 and 15 Pa.The model consists of three layers: OSB from outside (0.02 m), loose insulation made of hemp shives (0.1, 0.2 and 0.3 m) and gypsum board from inside (0.012 m).The models height is 1.0 m.Automatic discretization of models with stretch factor 1.3 divided them into 3232 elements for the model with a thermal insulation thickness of 0.1 m as well as 0.2 m, and 3542 for the model with a thickness of 0.3 m.
Only temperatures and pressure difference were set in the model.The external temperature corresponds to the computational temperature of Warsaw: te = −20.00• C, internal as for living room ti = 20.00• C, while pressure difference ∆p = 0.00, 5.00, 10.00 and 15.00 Pa correspond to possible conditions during wintertime.The pressure difference is subjected to 2 mm leakages situated 100 mm from the bottom and the top of the model.

Model Description
The energy balance is given with the Equation ( 4) is described in detail in [46].
where ρU is the internal energy density (J/m 3 ), λ is the heat conduction coefficient (W/(m•K), ∇T is the temperature gradient (K/m), κ is the intrinsic permeability of a porous medium (m 2 ), c p is the specific heat of dry air (J/(kg•K)), ∆T is the temperature difference (K), ν is the kinematic viscosity of air (m 2 /s), ∇p is the pressure gradient (Pa/m), ρ g is the dry air density (kg/m 3 ), g is the gravity acceleration (m/s 2 ) and ρ r is the internal heat source density(W/m 3 ).

Simulation Results and Discussion
Figures 10-12 present the heat flux profiles through the wall model depending on the thickness of the insulation (100 mm, 200 mm and 300 mm), density of hemp shives insulation (110 kg/m 3 and 120 kg/m 3 ) and air pressure (0 Pa, 5 Pa, 10 Pa and 15 Pa).The averaged heat flux was calculated for the internal surfaces of the model.The graphs show the period of 42 h, because this is the maximum time of stabilization of the heat flux flow (in the case of 300 mm thickness of the insulation-see Figure 12).In most cases, the stabilization was much faster, but in all cases the same period was used for comparative purposes.

Simulation Results and Discussion
Figures 10-12 present the heat flux profiles through the wall model depending on thickness of the insulation (100 mm, 200 mm and 300 mm), density of hemp shives i lation (110 kg/m 3 and 120 kg/m 3 ) and air pressure (0 Pa, 5 Pa, 10 Pa and 15 Pa).The a aged heat flux was calculated for the internal surfaces of the model.The graphs show period of 42 h, because this is the maximum time of stabilization of the heat flux flow the case of 300 mm thickness of the insulation-see Figure 12).In most cases, the sta zation was much faster, but in all cases the same period was used for comparative poses.Energies 2022, 15, 13

Simulation Results and Discussion
Figures 10-12 present the heat flux profiles through the wall model depending on thickness of the insulation (100 mm, 200 mm and 300 mm), density of hemp shives i lation (110 kg/m 3 and 120 kg/m 3 ) and air pressure (0 Pa, 5 Pa, 10 Pa and 15 Pa).The a aged heat flux was calculated for the internal surfaces of the model.The graphs show period of 42 h, because this is the maximum time of stabilization of the heat flux flow the case of 300 mm thickness of the insulation-see Figure 12).In most cases, the sta zation was much faster, but in all cases the same period was used for comparative poses.Energies 2022, 15, 13

Simulation Results and Discussion
Figures 10-12 present the heat flux profiles through the wall model depending on thickness of the insulation (100 mm, 200 mm and 300 mm), density of hemp shives i lation (110 kg/m 3 and 120 kg/m 3 ) and air pressure (0 Pa, 5 Pa, 10 Pa and 15 Pa).The a aged heat flux was calculated for the internal surfaces of the model.The graphs show period of 42 h, because this is the maximum time of stabilization of the heat flux flow the case of 300 mm thickness of the insulation-see Figure 12).In most cases, the sta zation was much faster, but in all cases the same period was used for comparative poses.In the presented figures (Figures 10-12) it can be clearly seen that air filtration causes the rise of heat fluxes for both analyzed densities of hemp shives, in the case of all thicknesses.Comparing the two tested insulation densities, it is visible that for a lower density, the increase in heat flux occurs more rapidly after air filtration begins.The higher the pressure difference, the faster the changes occur.The heat fluxes depending on the pressure difference also reach higher values in this case.This is due to greater air permeability.For both densities, there is a tendency that as the thickness of the insulating material increases, the increase in the heat flux value takes longer to reach the maximum.This is characteristic of loose materials subjected to air filtration.Similar observations in the case of laboratory tests on blowing loose mineral wool were made by Kosi ński et al. [50].It is worth noticing that in the no-filtration cases, the difference between the results is no higher than 3% in the corresponding insulation thicknesses.The difference between the results is up to 62% in the air filtration cases.
The dependences of the heat flux on the pressure difference was analyzed.They were made for two thicknesses (100 and 300 mm) of insulation and two material densities of 110 and 120 kg/m 3 (Figures 13 and 14).The abscissa axis shows the values obtained for the reference variant, without the pressure difference on both sides of the wall, and on the ordinate axis the values of the heat flux at the following pressure differences on both sides of the wall: 5, 10 and 15 Pa.The data are presented for a time interval of 70 h with an interval of one hour.In the presented figures (Figures 10-12) it can be clearly seen that air filtration causes the rise of heat fluxes for both analyzed densities of hemp shives, in the case of all thicknesses.Comparing the two tested insulation densities, it is visible that for a lower density, the increase in heat flux occurs more rapidly after air filtration begins.The higher the pressure difference, the faster the changes occur.The heat fluxes depending on the pressure difference also reach higher values in this case.This is due to greater air permeability.For both densities, there is a tendency that as the thickness of the insulating material increases, the increase in the heat flux value takes longer to reach the maximum.This is characteristic of loose materials subjected to air filtration.Similar observations in the case of laboratory tests on blowing loose mineral wool were made by Kosiński et al. [50].It is worth noticing that in the no-filtration cases, the difference between the results is no higher than 3% in the corresponding insulation thicknesses.The difference between the results is up to 62% in the air filtration cases.
The dependences of the heat flux on the pressure difference was analyzed.They were made for two thicknesses (100 and 300 mm) of insulation and two material densities of 110 and 120 kg/m 3 (Figures 13 and 14).The abscissa axis shows the values obtained for the reference variant, without the pressure difference on both sides of the wall, and on the ordinate axis the values of the heat flux at the following pressure differences on both sides of the wall: 5, 10 and 15 Pa.The data are presented for a time interval of 70 h with an interval of one hour.It was found that the increase of the pressure gradient on both sides of the simulated wall has a great influence on the heat exchange flux.In a model with a small thickness (100 mm) of insulation, this difference is smaller than in the case of a thicker insulation (300 mm).In the case of the insulation with a density of 120 kg/m 3 , the increase in flux over time was linearly proportional to the reference flux (at 0 Pa), and the maximum values obtained after the longest exposure time did not exceed 3 W/m 2 (18%) compared to the initial state.In the case of the insulation with a density of 110 kg/m 3 , the differences were more significant and for the 15 Pa pressure gradient reached 49%, it should be noted that after a period of about 8 h the differences are constant in relation to each other.In the case of a thicker wall (Figure 14), the influence of the pressure difference is more pronounced.For higher density insulation, a linear dependence of the heat flux was found, different to the reference situation (0 Pa pressure gradient).It was found that for the 15 Pa pressure gradient, the difference between the heat fluxes were 88%, and for the density 110 kg/m 3 it was 309%.It was found that the increase of the pressure gradient on both sides of the simulated wall has a great influence on the heat exchange flux.In a model with a small thickness (100 mm) of insulation, this difference is smaller than in the case of a thicker insulation (300 mm).In the case of the insulation with a density of 120 kg/m 3 , the increase in flux over time was linearly proportional to the reference flux (at 0 Pa), and the maximum values obtained after the longest exposure time did not exceed 3 W/m 2 (18%) compared to the initial state.In the case of the insulation with a density of 110 kg/m 3 , the differences were more significant and for the 15 Pa pressure gradient reached 49%, it should be noted that after a period of about 8 h the differences are constant in relation to each other.
In the case of a thicker wall (Figure 14), the influence of the pressure difference is more pronounced.For higher density insulation, a linear dependence of the heat flux was found, different to the reference situation (0 Pa pressure gradient).It was found that for the 15 Pa pressure gradient, the difference between the heat fluxes were 88%, and for the density 110 kg/m 3 it was 309%.
Figure 15   In the case of a thicker wall (Figure 14), the influence of the pressure difference is more pronounced.For higher density insulation, a linear dependence of the heat flux was found, different to the reference situation (0 Pa pressure gradient).It was found that for the 15 Pa pressure gradient, the difference between the heat fluxes were 88%, and for the density 110 kg/m 3 it was 309%.The temperature distribution presented in Figure 15 indicates the range of changes in heat distribution in the frame walls insulated with loose hemp shives subjected to air filtration.It is clearly visible that with the rising air pressure, the larger area of the model The temperature distribution presented in Figure 15 indicates the range of changes in heat distribution in the frame walls insulated with loose hemp shives subjected to air filtration.It is clearly visible that with the rising air pressure, the larger area of the model is affected by the cold outside air.Interestingly, the changes in the internal temperature distribution for the thinnest insulation model appear to be the smallest.
The analysis of the temperature changes on the internal surface in the contact point between the leakages and the gypsum board show a clear influence of the wind pressure on the temperature drop on the inner surface of the wall.Larger temperature drops were observed in the wall with the shives of lower density.This is understandable, because the looser structure of the insulation results in larger spaces between the shives, which allows the wind to flow more easily.The temperature drop in all cases occurs in direct proportion to the increase in wind pressure.At the strongest wind action, the temperature dropped, depending on the thickness of the glass, to 9.6-11.2• C for a wall comprising the shives with density of 110 kg/m 3 and to 15.3-15.6 • C for a wall insulated with the shives characterized by density of 120 kg/m 3 .The temperature without filtration was 17.9-19.2• C for a wall comprising the shives with density of 110 kg/m 3 and 17.1-18.3• C for a wall with the shives characterized by density of 120 kg/m 3 .Interestingly, in the case of a wall insulated with the shives of lower density, the temperature decreased as the wall insulation thickness increased.The temperature stabilization time is comparable (approx.10 h), regardless of the wind pressure.
In another study [10], the temperature on the inner surface of a flat roof insulated with flax shives (lambda value equals 0.049 W/(m•K) and density equals 145 kg/m 3 ) with a layer thickness of 15 cm, 20 cm, 25 cm and 30 cm were investigated.The temperature ranged from 18.6 • C to 19.3 • C with −20 • C outside and 20 • C inside.However, the influence of the wind was not taken into account in these analyses.

Conclusions
The research on the selected physical properties of hemp shives as the natural thermal insulation material enabled formulating the following conclusions: • Thermal conductivity 0.049-0.052W/(m•K), measured in the density range of 109-124 kg/m 3 , shows a little tendency to rise with density.The density of 109 kg/m 3 was the minimum value that could be achieved with this test.

•
Air permeability in the range of 4.259744 × 10 −8 -1.832226 × 10 −9 m 2 , measured in the density range of 90-120 kg/m 3 , shows a decreasing tendency with increasing density.The density of 90 kg/m 3 was the minimum value that could be achieved with this test.

•
The porosity tests showed that the individual elements of the hemp shives have a porosity of more than 78%.The largest share, i.e., about 40%, corresponds to the pores with a diameter in the range of 0.9-3 µm.

•
The results of the Rayleigh number calculation show that along with the rising density of shives, the share of convection will be limited.However, the results of the simulations indicate that in order to minimize the heat losses caused by convection, the thermal insulation made of shives should be covered with air-tight membranes.
The tested physical parameters, depending on the density of the shives mixture, made it possible to simulate the heat flux through a partition insulated with shives of different density.On the basis of the simulation results, the following conclusions were drawn:

•
In the case of insulation with lower density, the increase in heat flux occurs more rapidly after air filtration begins.An increase in the pressure difference enhances the rate of changes.The maximum value of heat flux is reached after a longer time, the thicker the insulation is.

•
Larger temperature drops on the inner surface of the wall, due to the air pressure were observed in the wall with shives of lower density.
The results of these studies indicate that hemp shives can be used well as a loose fill thermal insulation material in frame constructions.Currently, the authors continue the research on determining the moisture properties of hemp shives.The results, combined with simulation calculations will give broader view of the performance of loose hemp shives in constructions.

Figure 1 .
Figure 1.Hemp shives used in investigation (a) and their length (b).

Figure 1 .
Figure 1.Hemp shives used in investigation (a) and their length (b).

Figure 2 .
Figure 2. Heat flow meter-Fox 602 apparatus with hemp shives sample formed in the XPS frame.

Figure 2 .
Figure 2. Heat flow meter-Fox 602 apparatus with hemp shives sample formed in the XPS frame.
square resistance.This phenomenon is described by the Forchheimer equation presented by (3):

Figure 3 .
Figure 3.The apparatus for air permeability measurements filled with hemp shives.

Figure 3 .
Figure 3.The apparatus for air permeability measurements filled with hemp shives.Pressure difference was measured by a Testo 435-4 electronic differential pressure gauge (produced by Testo SE & Co. KGaA, Titisee-Neustadt, Germany), with measuring range 0-25.0 hPa and resolution 0.01 hPa.Air velocity was recorded by a TA7 thermoanemometer (produced by Airflow Developments Limited, High Wycombe, UK) with measuring range 0-2.0 m/s and resolution 0.01 m/s.Air was compressed in the screw compressor, equipped with the filter and dryer, and delivered to the apparatus by a pressure reducer.Testing the samples involved placing them in a central part of the system.Every prepared sample was created by random selection of conditioned loose hemp shives.The authors prepared the samples differing with density, from 90 to 120 kg/m 3 with 5 kg/m 3 resolution.During the measurements, the pressure difference between two taps, located in central part of the apparatus, was registered.Simultaneously, air velocity on the exhaust was registered.In order to avoid discrepancy of the density and viscosity of the air during investigation, a constant temperature +10 • C of compressed air was maintained.Every single calculated value of permeability was determined from the approximation of (3) on the basis of at least 800 pairs of air velocities and air pressures.Due to the technical declaration of producers, the standard measurement uncertainty of the pressure difference is ±2 Pa in the measuring range, while the uncertainty of air velocity is ±0.04 m/s.Pressure uncertainty can be practically neglected in the achieved range (up to 6000 Pa) of the pressure difference, but the uncertainty of air velocity is very important.The authors estimated this influence: underestimation of air velocity can result in up to 29% overestimation of air permeability, while overestimation of air velocity can lead up to 55% of underestimation of air permeability.

Figure 6 .
Figure 6.Thermal conductivity of hemp shives for different achieved bulk densities of samples (error bars mean standard deviation).

Figure 6 .
Figure 6.Thermal conductivity of hemp shives for different achieved bulk densities of samples (error bars mean standard deviation).

Figure 8 .
Figure 8. Dependence of the Rayleigh number on the density and thickness of thermal insulation made of hemp shives.

Figure 8 .
Figure 8. Dependence of the Rayleigh number on the density and thickness of thermal insulation made of hemp shives.

Figure 10 .
Figure 10.Heat conduction flux through the wall with 100 mm thickness of hemp fiber insula with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 11 .
Figure 11.Heat conduction flux through the wall with 200 mm thickness of hemp fiber insula with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 10 .
Figure 10.Heat conduction flux through the wall with 100 mm thickness of hemp fiber insulation with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 10 .
Figure 10.Heat conduction flux through the wall with 100 mm thickness of hemp fiber insul with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 11 .
Figure 11.Heat conduction flux through the wall with 200 mm thickness of hemp fiber insul with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 11 .
Figure 11.Heat conduction flux through the wall with 200 mm thickness of hemp fiber insulation with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 10 .
Figure 10.Heat conduction flux through the wall with 100 mm thickness of hemp fiber insula with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 11 .
Figure 11.Heat conduction flux through the wall with 200 mm thickness of hemp fiber insula with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 12 .
Figure 12.Heat conduction flux through the wall with 300 mm thickness of hemp fiber insulation with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 12 .
Figure 12.Heat conduction flux through the wall with 300 mm thickness of hemp fiber insulation with density of: 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 13 .
Figure 13.Comparison of the dependence of the heat flux in a wall with 100 mm insulation thickness on the pressure difference in the form of a comparison of the reference variant (0 Pa) with variants with a pressure difference of 5, 10 and 15 Pa, for a density of 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 13 . 19 Figure 14 .
Figure 13.Comparison of the dependence of the heat flux in a wall with 100 mm insulation thickness on the pressure difference in the form of a comparison of the reference variant (0 Pa) with variants with a pressure difference of 5, 10 and 15 Pa, for a density of 110 kg/m 3 (a) and 120 kg/m 3 (b).Energies 2022, 15, 2461 15 of 19

Figure 14 .
Figure 14.Comparison of the dependence of the heat flux in a wall with 300 mm insulation thickness on the pressure difference in the form of a comparison of the reference variant (0 Pa) with the variants characterized by a pressure difference of 5, 10 and 15 Pa, for a density of 110 kg/m 3 (a) and 120 kg/m 3 (b).
presents the temperature distribution in the 110 kg/m 3 model for three different insulation thicknesses: 100, 200 and 300 mm in the case of no air filtration and the filtration caused by different air pressure: 5 Pa, 10 Pa and 15 Pa.

Figure 14 .
Figure 14.Comparison of the dependence of the heat flux in a wall with 300 mm insulation thickness on the pressure difference in the form of a comparison of the reference variant (0 Pa) with the variants characterized by a pressure difference of 5, 10 and 15 Pa, for a density of 110 kg/m 3 (a) and 120 kg/m 3 (b).

Figure 15
presents the temperature distribution in the 110 kg/m 3 model for three different insulation thicknesses: 100, 200 and 300 mm in the case of no air filtration and the filtration caused by different air pressure: 5 Pa, 10 Pa and 15 Pa.

Figure 15 .
Figure 15.The temperature distribution in the 110 kg/m 3 models: no filtration (a), 5 Pa pressure difference (b), 10 Pa pressure difference (c) and 15 Pa pressure difference (d).

Figure 15 .
Figure 15.The temperature distribution in the 110 kg/m 3 models: no filtration (a), 5 Pa pressure difference (b), 10 Pa pressure difference (c) and 15 Pa pressure difference (d).

Table 1 .
Structure characteristics determined by mercury porosimetry.

Table 1 .
Structure characteristics determined by mercury porosimetry.

Table 2 .
Thermal conductivity of hemp shives, a literature review.

Table 3 .
Results of the air permeability for investigated densities of hemp shives.

Table 4 .
Physical properties of materials used in the simulation.