The tendency of the 21st century in industrialised nations is to limit the impact of human activity on the environment. It is thus within this context that construction industries had to innovate its construction practices in order to improve the use of energy in buildings as well as proposing cutting-edge technologies on materials that meet new user demands and environmental impact laws. Greenhouse gases emission, which makes up one of the main cause of environmental pollution, is constantly on the rise. These gases touch a wide range of sectors such as: agriculture, waste treatment plant, industries, transportation, and building construction etc. The global emission rate of greenhouse gas (GHG) in the Building sector is 4.1 Gt eq CO2 (gigatonnes in CO2 equivalent) in 2016  . According to Meggers et al.  the building sector alone is responsible for 35% of world emission of carbon dioxide (CO2). This figure is steadily on the rise and can reach 50% by the year 2050. In Africa, particularly in Cameroon, the figures are alarming. Since 2012, energy production is not commensurate to its consumption. The highest energy consumption units have artificial air-conditioners to ensure better heat control. For a better control of energy use in a building, the use of a good local isolating material with high emitting heat ability and hygroscope is recommended. Heat emission will help curb the high energy consumption in the building currently with air-conditioner that doesn’t meet current international standards. This situation is very regular in the Sahelian regions of Cameroon with its hot tropical climate  . However, in order to reduce global warming and improve on the heat emission in buildings, laterite soil mixture is used to make composite material applied in the interior of the building to ensure its better heat insulation.
Our research is therefore, based on this concept of the used of laterite composite material as insulant in the interior of buildings. The sample composite material used in this study is laterite bricks mixed with thatch fibres all from the locality of Meiganga. In the Adamawa region of Cameroon, otherwise known as Watershed, 80% of houses are built with laterite brick not industrially treated nor having any geotechnical, mechanical, thermophysical or acoustic characterisation prepared initially. The bricks used do not meet civil engineering standards. A lot of research works have been carried out in the domain of characterisation of composite material made from laterite soil. Meukam et al.  have shown that cement stabilized clay bricks mixed with pozzolan have good mechanical strengths and low thermal conductivities. Lamharouet et al.  demonstrated how thermal characteristics of Larache locality (Morocco) change considerably when stabilized with a proportion of cement. Bal et al.  on their part studied the lateritic clay bricks of Senegal. They used the asymmetric hot plate method to estimate the thermal conductivities when used in a mixture with millet pod fibres. They went further to show that mixture with that fibre drastically reduces its thermal conductivity. But as limitation, his work didn’t show to which extend the fibre content of the composite material created respected Civil Engineering Standards. Meukam et al.  worked on the thermal-mechanical characteristics of laterite bricks in the Centre region (Yaounde) of Cameroon solidified with cement on one hand and on the other with the use of a substance like pouzzolan. Thermophysical characterisation with box method confirmed that using pozzolan, leads to better insulating material. Damfeu et al.  modelled the thermal properties of insulating plant fibres using the box asymmetric hot plate and radial flow methods. They demonstrated that the apparent thermal conductivity of kapok and coconut fibre could be used as substitute to synthetic insulating materials such as polyester or glass fibre. Probably, the use of an additional element in making a laterite brick would have led to good insulating composite material. Damfeu et al.  in another write up show the modelling and assessing of aggregate thermal properties of construction materials and went further to justify that red pouzzolan aggregates has a smaller thermal effusivity and conductivity than the black aggregate pouzzolan. Their complementary role in the making of laterite brick could therefore improve its thermal performance. Elhamdouni et al.  studied the thermophysical and mechanical characterisation of reinforced clay bricks with alfa or straw fibres. The results showed that the use of straw fibres help reduce the thermal conductivity energy to 2% - 7% as compared to alfa fibre. Nevertheless, they went further to demonstrate that alfa fibre can improve mechanical performance of clay material as compared to clay straw material (with an increase of 8% - 17% of the flexural strength and 6% - 18% of compressive strength). Although most of these studies seek to improve the thermal insulating capacity of the sample materials, the findings still did not corollate with that of the pilot studies, especially with that of the apparent thermal conductivity. These results do not always reflect the conditions of its use compared to Civil Engineering standards.
This work has as main objective the enhancement of thatch fibres as thermal insulation in buildings to reduce the quantity usually burnt during the long dry season. It was assessed when mixed with laterite bricks to get the thermophysical properties. This material can be used to curb energy consumption when used as filling for buildings.
2. Materials and Methods of Mechanical Characterization
Around the Meiganga vicinity, more than 80% of houses are built with laterite soil as opposed to the West Region of Cameroon where the laterite soil is rather dry. Although the objective is to get a composite material with insolation properties, it’s however imperative to begin with a mastery of the geotechnical properties of the sample laterite soil. Cameroon Civil Engineering standards stipulates (CS: 102 - 114) that, filling material for building such as compressed terracotta brick wall recommended for making laterite bricks should be used with plasticity index of Ip ≤ 30%   . The laterite (Figure 1) studied in this research project was taken from Meiganga. The water content of laterite soil sample as of escavation is 15% got from the Equation (1):
where: is the mass of soil sample in its natural state and is the mass of the soil sample after passing through heat for 105˚C for 24 hours.
In order to determine the physico-chemical and mineralogical characteristics of this laterite, a grading analysis (according to BS 1377: 1975  ) was carried out in a bit to check if the excavated laterite sample contained much clay and therefore, needed rectification before use. Following the size of particles obtained with the help of an electric sieve i.e. with diametre between 1 ≤ D (μm) ≤ 100 undergone a particle size grading study and that of diametre between 100 ≤ D (μm) ≤ 6500 underwent a sediment metric particle size grading study. Two equations were used to calculate the diametres of the particles used to represent the diagram of the texture of this laterite.
- The diametre of particles is given by the stockes Law
- For the percentage of the diametre particles inferior to D, we use the equation
where: is the bulk specific weight of the solid particles; is the bulk specific weight of water; g is accelerating gravity in cm・s−2; U is the falling speed of particles in cm・s−1, η is the fluid viscosity measured in poise; C3 is the conversion coefficient from 1/1000; C4 is the reading for water only and L is the vertical reading 1/1000 beyond one.
The results from Equations (2) and (3) enabled us to represent the texture diagram presented in Figure 2.
Distribution of particle size grading analysis of laterite soil as illustrated in Table 1.
In analysing the diagram in Figure 2, we noticed that the texture curve of the laterite soil did not entirely fall under the grading range, moreover, the analysis of the summary in (Table 1) show the clay nature of the soil with a (38.85%) which is well above the limit recommended by the Cameroonian standard CN: 102 - 114 for the design of BTC which is 30%    . This result is rather interesting to us because in the samples to be made, thatch fibres would be integrated without the use of stabilisers like cement or lime. It is needless to treat laterite before using because the more the use, the more clay it contains. We equally brought out the limitations of Atterberg’s work at the level of water and plastic limitation in laterite soil. Not having used the shock number (N) method directly equal to 25 shocks, we used the equation to calculate the liquidity and the plasticity limit of each test:
where: Mt tare mass; Mth is tare mass + wet sample; Mts is tare mass + dry sample
Placticity index IP of laterite came from the equation:
Table 2 shows limitation of Atterberg laterite
We can then conclude that though the laterite texture curve does not fall within the grading range, its plastic index meets BTC conception standards that stipulates that Atterberg’s limitation must be between 10 < Ip < 20. The laterite was not adulterated with in anyway.
Figure 1. Lateritic soil.
Figure 2. Diagram showing texture of Laterite from Meiganga.
Table 1. Summary analysis of particle size grading.
Late: Lateritic soil.
Table 2. Atterberg limitations.
2.1.2. Thatch Fibres
Several researchers have worked on fibres as natural insulators that help improve thermal performance of a composite material. Petkova et al.  worked on the thermal properties of insulation material such as straw and even went further to show that the use of straw (less than 0.5% in mass) improves heat accumulation and thermal insulating properties of the ecological composite material. Azahry et al.  carried out a study of the energy efficiency and thermal properties of the clay-straw of composite material thereby; high-lighting its important thermal properties for heat storage and temperature regulation during winter day and night. He also proved that this material can withstand extreme temperature (hot or cold) without any cooling or heating system. Kossi et al.  looked at the mechanical and thermal characterisation of compact soil blocks with straw or cement that have different percentages and even went further to show that the thermal conductivity decreases in relation to the proportion of fibres used. But our study dwelt on natural fibres from thatch as shown in Figure 3.
The process of extracting dry fibres uses as follow:
- The fibres were picked wet (Figure 3(a));
- Chopped at regular length and dried in the open (Figure 3(b));
- After the crushed stem mixed with lateritic soil (Figure 3(c));
- They were later dried in an oven at 45˚C (Figure 3(d)) for 48 hours;
- After that they were packed in plastic bags so as to keep their water content close to zero.
Measuring the apparent density of dry fibres (Table 3) shows the light fibres.
The dry fibres obtained are then mixed with the laterite at different proportions to get different samples of composite material.
2.1.3. Mixed Composition of Flexural and Compressive Strengths
The experimental setup used in making samples were made of three cuboid-shaped metallic mould of (width 1, length L and thickness h): 4 cm × 8 cm × 2 cm (French standards Crater standard NFP-31-212  ) for the Bending test (Figure 4(a)), another square shape of 4 cm × 4 cm × 4 cm (standard NFP-98-335  ) for the Compression test (Figure 4(b)) and of rectangular shape 10 cm × 10 cm × 3 cm for the thermophysical characterisation tests (Figure 4(c)). The hydraulic press used for compacting is suitable for sample design for three-point bending, tensile and simple compression tests. The specimens produced were maintained at 32 days of cure. Table 4 shows the different samples made.
(a) (b) (c) (d)
Figure 3. Natural thatch fibres: (a) stem of fibres; (b) dry stem; (c) crushed stem mixed with lateritic soil; (d) oven for drying sample.
(a) (b) (c)
Figure 4. Samples for: (a) bending test; (b) compression test; (c) thermophysical characterisation test.
Table 3. Apparent density of dry fibres.
Ms: dry mass of fibres; Vf: volume of fibre; ρapp(f): Apparent density of dry fibre.
Table 4. Samples composition.
mL (g): mass of laterite: mF (g): mass of fibres:; V (l): volume of water.
2.2. Mechanical Characterisation Approach
Mechanical characterisation is of prime importance in Civil Engineering. Depending on the nature of the material to be used, either as a load-bearing structure or as a filling of the building, this characterisation helps to verify whether the material meets Civil Engineering standards. In our case, filling materials were used in the interior of the building. We were therefore interested in two mechanical characteristics: Flexural and the Compressive strength.
2.2.1. Flexural Strength
Figure 5 shows the research tool used for measuring flexural strength
F (N) is flexural samples load (measured by CONTROLS.CAT machine No. T1004, Mat No. 3558), distance from cylindrical support, length of block and the average thickness of block. The flexural strength of MPa calculated using the equation:
2.2.2. Compressive Performance
Figure 6 shows the research tool (cubic block).
If F is the compressive break-down load (observed in reading the machine CONTROLS No. of rings 4117, Fmax = 30 kN), we then obtain RC constraint starting from the relation (h2 is the section):
3. Thermophysical Characterisation Approach
3.1. Testing Device
The testing tool used to show the thermophysical characterisation is the asymmetric hot plane method   . It’s a transient characterisation approach used to get the apparent thermal conductivity of a material (Equation (21)) by estimating the thermal effusivity E and the volumetric heat capacity ρCP in relation to the testing temperature Texp(t) and modelled temperature Tmodel(t). The testing instruments used for measuring the thermal properties are shown in Figure 7.
Figure 5. Research tool for flexural strength.
Figure 6. Compressive research tool.
Figure 7. A cross-section of test instruments.
Material with specifications of (10 cm × 10 cm × 3) is placed on a heat sensor between two polystyrene blocks of (10 cm × 10 cm × 5 cm). A generator is used to heat the resistor. The increase in temperature at the centre of the heat resistor was as a result of type K thermocouple that registers the test temperature on the hot surface of the material. The experimental temperatures Texp(t) were registered using the acquisition module TC 08-USB Picolog. The dropping of experimental and simulated temperature obtained after modelling the testing instrument helped us to estimate E and ρCP. It is therefore important to note that the temperature above the polystyrene blocks remains at its initial state, to arrive at this; we placed two aluminium blocks of 10 cm × 10 cm × 4 cm above and below the polystyrene blocks.
3.2. Quadrupoles Formalism of Heat Transfer
Given that all the blocks used for the pilot test were thick, we could therefore model heat transfer to other blocks using the heat equation.
where modelises the temperature thermogram on the heated face of the composite material while does same on the heat sensor. By assuming that the heat transfer remains unidirectional for a time t, solving this 1D equation can be done using the quadrupole formalism developed by Maillet et al.  . It is generally known that solving Equation (9) using the quadrupole formalism in the Laplace space leading us to a relation whose form is determined by Equation (10) below:
where: M is matrix or matrix product.
We can use this formalism according to the characterisation method used in evaluating thermophysical parameters E, ρCP and λ. In this study, the asymmetrical hot plate method enables us to estimate by using the complete 1D model, the parametres E and ρCP yielding to λ.
3.2.1. Simplified 1D Model
Considering the sample and polystyrene insulating block as a semi-infinite medium and hoping that the heating sensor is a thin device (Figure 8), the equations can be written following the thermal flow traces across the sample (Equation (11)) or across the heating block (Equation (12)):
Figure 8. Across of the heat flow.
where: E and Epo are respectively the thermal effusivity of samples and the polystyrene blocks while P is the Laplace parametre; S is the surface area of the sample which is supposed to be equal to Sh, area of heating element and Spo, area of the polystyrene blocks.
Combining Equations (11) and (12) yields to:
For a considerable period of time (p®0), the analytic inversion of the above equation, using Laplace table gives:
The pre-estimated thermal effusivity may be calculated with the simplified 1D model from numerical calculation of the slope δ(t) of the linear part of the curve when Texp(t) and Tsinf(t) are overlapped.
The pre-estimated value of volumetric heat capacity given by the equation (15), may be approximated from the slope η(t) of the linear part of the curve when Texp(t) and Tsinf(t) are overlapped.
These pre-estimated values (Equations (15) and (16)) are then used as initial values in the estimation code of E and ρCp for complete model.
3.2.2. Completed 1D Model
With the Complete model, the quadrupole method no longer uses the sample and insulating block as in semi-infinite medium. Consequently, matrix M representing these two blocks is shown in the Equation (17):
The solving of the matrix system equation is given by Equations (18) and (19):
The total heat flux density in the Laplace field is calculated from relation:
Combining relations (18); (19) and (20), the system leads to:
The theoretical curve was obtained by numerical inversion of the Equation (21) using the De Hoog algorithm  . The Levenberg Marquert algorithm   used with Matlab helps to calculate the E and ρCp parametres with the complete model. This assessment can only drop when converged (reducing the sum of the quadratic difference: Equation (22)) between the experimental curve and theoretical curve, taking as initial values the pre-estimated values (Equations (15) and (16)).
The calculated values of apparent thermal conductivity of composite material can also be inferred from the Equation (23):
4. Research Procedure for Measuring Thermal Capacity Cp
Establishing thermal capacity of composite material with different fibre proportions that was achieved had as objective to confirm the results of the pilot study carried out with Equation (23). This test achievement was made known through its water content, colometry sizing which helps to validate the results. To establish the water content, several methods do exist among which the dehydration method (DM) uses in this research. Three research instruments (Figure 9) were used to get the results of the pilot study: drying oven regulated at 105˚C for 24 h to dry the sample; Omega Engineering probe; Dewar Vapour calorimetre 12Vmax.
4.1. Method of Determining CP
This method was inspired by the manner in which objects transfer heat. The Cp was measured in cold water as well as normal temperature. The thermal heat emitted by the object is as follows:
In cases where heat loss is neglected within the exchange system, we can therefore conclude that the energy produced by the object is exactly equivalent to that absorbed by water:
where Qhw equals energy emitted by the object, Qcw energy absorbed by water and Qca energy absorbed by calorimeter.
For calorimeter sizing, the following equation was used:
If the mass of hot water is equal to that of cold water, we’ll then say that:
To determine the Cp sample; we have:
where: m1 mass of hot water, m2 mass of cold water, θc hot water temperature, θwo Cold water temperature, θeq balance point, Calorie capacity of calorimeter
4.2. Experimental Result of Cp
To validate the research method used, we tested the value of the thermal capacity of the calorimeter used. The specification sheet provided shows CP = 30 J・K−1. (Table 5).
We notice a correlation between the average value of the thermal capacity of calorimeter calculated with dehydration method and that shown on the researcher’s specification sheet. The relative difference obtained between the two values is less than 3%. This result enables us to validate the results of the thermal capacities of the composite materials studied. Table 6 presents the set of experimental results obtained for samples E0, E1, E2, E3 and E4.
5. Discussion of Findings
Measurements were carried out on five different samples E0, E1, E2, E3 and E4.
(a) (b) (c)
Figure 9. Testing material for establishing thermal heat capacity; (a) Adjustable oven; (b) Omega engineering probe; (c) Dewar vapor calorimetre.
Table 5. Equilibrium temperatures and thermal heat capacity Cp of the calorimeter obtained.
Mce: Mass of cold water; Mhe: Mass of hot water; Tce: cold water temperature; T: hot water temperature; D: deviation.
Table 6. Experimental thermal capacity of samples.
Three were performed on each sample as well as the average value of each of the samples. For the mechanical tests, Table 7 summarizes all the experimental results.
We notice that the apparent density ρapp and the porosity ε were subsequently obtained from the Equation (30):
where: Vapp is the apparent density of the composite material, Mh wet density of composite material, Ms is dry density of composite material and ρe is the mass density of water.
It should be noted that control protocol on building site recommended by MIPROMALO (Mission de Promotion des Matériaux Locaux)  BTC dry compressive strength and with a 14 day greater than 2 MPa. Generally, cases with filling material of 0.6 ≤ RC ≤ 1 MPa blocks should be used to construct houses with first floor and large borders. Figure 10 represents the compressive strength of bricks made according to the proportion of fibres. We noticed that as quantity of fibres increases, so is the solidity of the compressive strength of the composite material. The reverse was supposed to be true, but since the lateritic soil used is very rich in clay, these thatch fibres from the Adamawa region, known for its firm stems (especially when dried) would however, naturally increase the strength of these bricks. At the same time adding fibres increase the porosity (Figure 11) of the material and being too clayish, it therefore creates micropores that help increase the insulation power as well as its flexural strength. (Figure 12) Looking at the most general case, that is: standards of construction on the compressive samples RC ≥ 2 MPa of the observation results, we found out that the brick can be used even at a grade of 4% in fibres content.
Figure 10. Compressive strength of 28 days according to the fibrous content.
Figure 11. Compressive strength of 28 days according porosity of composites materials.
Figure 12. Flexural strength to 28 days according to the fibrous content.
Table 7. Result of mechanical experiment.
As concerns thermophysical characterisation, we initially verified if the new model could react to thermophysical properties when calculated with transient hot plate method. Since these thermophysical parametre are not yet known, we noticed a reduced sensitivity of the thermophysical parametres by using the result obtained with the simplified model (Equations (15) and (16)) Values Epres and (ρCp)pres were calculated within a period of time [t1, t2] as much as the thermogram Tasym(t) and Tsinf(t) superimposed (Figure 13(a)). For example, with the case of sample E3, the numerical calculations of slope δ and η (Figure 13(b)) between limits 100 - 150 s gave us: Epres = 878.34 J・m−2 K−1・s−1/2 and (ρCp)pres = 1.531 × 106 J・m−3・K−1.
The reduced sensitivities of the temperature of parameter E, ρCp and SRchs were calculated as presented in Figure 14.
Figure 13. Temperature curve for sample E3 (a) T = f(t); (b) T = f(rac(t)) with residues curves.
Figure 14. Reduced sensitivities curves: (a) for samples E0; (b) for samples E3.
- We noticed that the reduced sensitivities of T to the heating sensor becomes constant after 10 s, this shows that the inertia of the heating sensor could be neglected only 10 s after having measured T. This shows that, although the effect of the inertia becomes important as the proportion of fibres in the material increases, the developed model cannot provide a reliable estimate of the thermal contact resistance at the interface of the heat sensor/material.
- Sensitivities of T to ρCp and to E are decorrelated, but T is sensible to E only between 0 and 500 s and sensible to ρCp practically after 300 s. The link observed between ρCp and Ri is as a result of strong thermal initial observed 100 s after having started registering the temperature. This could be fully seen on the residual graph that shows the plotting in Figure 14(b). After having well focused the residual at the centre between 100 s and 220 s and the temperature being very sensitive to E within the gap limit (Figure 14(a) and Figure 14(b), a pre-estimate of E was carried out on all the samples within the set interval of 100 - 220 s. On the contrary ρCp pre-estimate was done from the moment there was no correlation with Ri, i.e. between 400 s and 450 s.
Contrary to the simplified model which does not reduce the sum of quadratic differences (Equation (22)) between the model and the experimental temperature, the completed model instead reduces inversely. We therefore observed a convergence between pre-estimated and experimental results as seen in Figure 15, explaining why the completed model portrays a reliable estimation of the thermophysical parameters ρCp and E.
In a bid to validate and confirm the experimental results with the asymmetric hot plate method, we compared the estimated value of (ρCp)est with that of ρapp × (Cp)DM where the value of (Cp)DM and papp are respectively shown in Table 6 and Table 7. We equally compared the value of λest (Equation (23)) with λexp as shown by the Equation (32):
All the results obtained are shown in Table 8.
When comparing results from both methods which aim at determining the volumetric heat capacity, we can clearly see the mean value of apparent thermal conductivity obtained by the asymmetric hot plate method on one hand and that obtained by the Equation (32), perfectly match (Figure 16). As the use of fibre increases the porosity of composite material, so too decreases the thermal conductivity with the addition of fibres. The results of the experiment showed that blocks E4 have a much lower thermal conductivity than the blocks E0. We came to the conclusion that the use of plant fibres increases the weightlessness of the material without reducing its resistance to compression. This can also mean a reduced thermal conductivity.
Table 8. Results of the thermophysical properties.
Simpl: simplified model; compl: completed model; est: estimated; exp: experimental; D: standard deviation.
Figure 15. Results mixture with the complete model for the E4 sample.
Figure 16. Thermal conductivity of materials according to the fibrous content (ASYM: obtained from Equation (23); DM: obtained with Equation (32)).
This research work sets out to investigate the energetic and mechanical reinforcement of sun-dried bricks made with laterite and different proportions of thatch fibres from the Meiganga locality of the Adamawa region in Cameroon. Mechanical characterisation studies have shown that with up to 4% fibre incorporated, the composite materials produced have a very good compressive and flexural strength that meets civil engineering standards. The asymmetric hot plate method used to calculate the thermophysical parametres of these different materials was validated with another method, i.e. the dehydration method, used for determining the heat capacity. With the use of apparent thermal conductivity study, we came to the conclusion that with fibres used, a better thermal insulating composite material was the bi-product. These materials can then contribute to ensuring a better thermal comfort of the interior of a building. The extremely hot climate of this region gives credence to the treatment and use of plant fibres to make sun-dried bricks. These could be used to reduce energy consumption thus, limiting the emission of greenhouse gas (especially CO2) when air-conditioners are used.
T. Temperature (˚C)
E. Thermal effusivity (J・m−2・C−1・s−1⁄2)
λ. Thermal conductivity (W・m−1・K−1)
h. Convective heat loss coefficient (W・m−2・˚C−1)
R. Thermal contact resistance (˚C−1・W−1)
Cp. Heat capacity (J・kg−1・K−1)
ρ. Density (kg・m−3)
e. Thickness (m)
θ. Laplace transform of temperature
P. Laplace parameter
Φ. Laplace transform of heat flux
φ0. Heat flux density (W・m−2)
ρCp. Volumetric heat capacity (J・m−3・K−1)
P0. Insulating blocks
comp. Completed model
simpl. Simplified model
h. Heation element
DM. Dehydration Method
Sinf. Semi-infinite medium
 Meggers, F., Leibundgut, H., Kennedy, S., Qind, M., Schlaiche, M., Sobekf, W. and Shukuyag, M. (2012) Reduce CO2 from Buildings With Technology to Zero Emissions. Sustainable Cities and Society, 2, 29-36.
 Osseni, S.O.G, Apovo, B.D., Ahouannou, C., Sany E.A. and Jannot, Y. (2016) Caractérisation thermique des mortiers de ciment dopés en fibres de coco par la méthode du plan chaud asymétrique à une mesure de température. Afrique Science, 12, 119-129.
 Meukam, P., Noumowe, A., Jannot Y. and Duval, R. (2003) Caractérisation thermophysique et mécanique des briques de terre stabilisées en vue de l’isolation thermique de batiment. Materials and Structures, 36, 453-446.
 Lamkharouet, N., Boussaid, S., Ezbakhe, H., EL Bakkouri, A., Ajzoul, T. and El Bouard, A. (2001) Etude thermique de la terre de Larache stabilisée au Ciment. Laboratoire de Thermique, Energie solaire & Environnement, Faculté des Sciences, Université Abdelmalek Essaadi, Tétouan, 541-546.
 Bal, H., Jannot, Y., Gaye, S. and Demeurie, F. (2013) Measurement and Modelisation of the Thermal Conductivity of a Wet Composite Porous Medium: Laterite Based Bricks with Millet Waste Additive. Construction and Building Materials, 41, 586-593.
 Meukam, P., Jannot, Y., Noumowe, A. and Kofane, T.C. (2004) Thermophysical Characteristics of Economical Building Materials. Construction and Building Materials, 18, 437-443.
 Damfeu, J.C., Meukam, P. and Jannot, Y. (2016) Modeling and Measuring of the Thermal Properties of Insulating Vegetable Fibers by the Asymmetrical Hot Plate Method and the Radial Flux Method: Kapok, Coconut, Groundnut Shell Fiber and Rattan. Thermochimica Acta, 630, 64-77.
 Damfeu, J.C., Meukam, P., Jannot, Y. and Wati, E.(2017) Modelling and Experimental Determination of Thermal Properties of Local Wet Building Materials. Energy and Buildings, 135, 109-118.
 Elhamdouni, Y., Khabbazi, A., Benayad, C., Mounir, S. and Dadi, A. (2017) Thermophysical and Mechanical Characterization of Clay Bricks Reinforced by Alfa or Straw Fibers. IOP Conference Series: Materials Science and Engineering, Vol. 186, Beijing, 12-35.
 Azhary El, K., Chihab, Y., Mansour, M., Laaroussi, N. and Garoum, M. (2017) Energy Efficiency and Thermal Properties of the Composite Material Clay-Straw. Energy Procedia, 141, 160-164.
 Kossi, I.B., Florent, P.K. and Ouedraogo, E. (2014) Thermal and Mechanical Study of the Adobe Stabilized with Straws and/or Cement at Different Rate. International Journal of Engineering and Innovative Technology, 4.
 De Hoog, F.R. and Stokes, A.N. (1982) A Improved Method for Numerical Inversion of Laplace Transforms. SIAM Journal on Scientific and Statistical Computing, 3, 357-366.