Thermal Conductivity and Microstructure of Novel
Flaxseed-Gum-Filled Epoxy Resin Biocomposite:
Analytical Models and X-ray Computed Tomography

Mohammed Zaidi, Dominique Baillis, Naim Naouar, François Delattre,
Michaël Depriester

► To cite this version:

Mohammed Zaidi, Dominique Baillis, Naim Naouar, François Delattre, Michaël Depriester. Thermal Conductivity and Microstructure of Novel Flaxseed-Gum-Filled Epoxy Resin Biocomposite: Analytical Models and X-ray Computed Tomography. Materials, 2023, 16 (18), pp.6318. 10.3390/ma16186318. hal-04213221

HAL Id: hal-04213221

https://hal.science/hal-04213221v1

Submitted on 21 Sep 2023

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
Article

# Thermal Conductivity and Microstructure of Novel Flaxseed-Gum-Filled Epoxy Resin Biocomposite: Analytical Models and X-ray Computed Tomography

Mohammed Zaidi <fnref n="1" />,<fnref n="2" />*, Dominique Baillis <fnref n="2" />, Naim Naouar <fnref n="2" /> ID, Michael Depriester <fnref n="3" /> and François Delattre <fnref n="1" />

**Abstract:** The growing awareness of the environment and sustainable development has prompted the search for solutions involving the development of bio-based composite materials for insulating applications, offering an alternative to traditional synthetic materials such as glass- and carbon-reinforced composites. In this study, we investigate the thermal and microstructural properties of new biocomposite insulating materials derived from flaxseed-gum-filled epoxy, with and without the inclusion of reinforced flax fibers. A theoretical approach is proposed to estimate the thermal conductivity, while the composite's microstructure is characterized using X-ray Computed Tomography and image analysis. The local thermal conductivity of the flax fibers and the flaxseed gum matrix is identified by using effective thermal conductivity measurements and analytical models. This study provides valuable insight into the thermal behavior of these biocomposites with varying compositions of flaxseed gum and epoxy resin. The results obtained could not only contribute to a better understanding the thermal properties of these materials but are also of significant interest for advanced numerical modeling applications.

**Keywords:** flax fiber; flaxseed gum; biocomposite material; thermal conductivity; analytical models; X-ray computed tomography

## 1. Introduction

Today, increasing global consumption is accelerating the scarcity of resources, and developing alternatives poses a major challenge, particularly in the context of climate change and the replacement of fossil fuels [1]. Consequently, it is time to promote the use of renewable resources whose agro-sourced origin would enable the integration of a circular economy with a regional and local dimension. The use of agro-resources, in particular agricultural residues, holds great interest for the development of biocomposite materials, offering both environmental benefits and competitive production costs [2]. Natural biocomposites have great potential in a wide range of applications, including in the automotive industry, packaging, and household goods [3,4]. Indeed, they can serve as cost-effective materials while offering a wide range of structural properties [5]. Fiber-reinforced biocomposites, especially those derived from flax, are emerging as promising substitutes for synthetic fibers in polymeric composites. They have attracted interest due to their ecological attributes, cost-effective production, and mechanical characteristics. Several studies have taken advantage of them in various applications, particularly the automotive industry and building insulation [6–8].

Due to their interesting mechanical and thermal properties, the characteristics of flax-fiber-based biocomposites was investigated by various researchers [9–13]. In [14], a

<figure><img src="image_2.png" /></figure>

### check for updates

**Citation:** Zaidi, M.; Baillis, D.; Naouar, N.; Depriester, M.; Delattre, F. Thermal Conductivity and Microstructure of Novel Flaxseed-Gum-Filled Epoxy Resin Biocomposite: Analytical Models and X-ray Computed Tomography. *Materials* **2023**, *16*, 6318. https://doi.org/10.3390/ma16186318

Academic Editor: Swarup Roy

Received: 2 August 2023

Revised: 11 September 2023

Accepted: 14 September 2023

Published: 20 September 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Unité de Chimie Environnementale et Interactions sur le Vivant (UR 4492, UCEIV), SFR Condorcet FR CNRS 3417, Littoral Côte d'Opale University, 145 Avenue Maurice Schumann, 59140 Dunkerque, France</footnote>

<footnote n="2">² LaMCoS, INSA-Lyon, CNRS UMR 5259, Université de Lyon, 69621 Villeurbanne, France</footnote>

<footnote n="3">³ Unité de Dynamique et Structure des Matériaux Moléculaires (UR 4476, UDSMM), Littoral Côte d'Opale University, 145 Avenue Maurice Schumann, 59140 Dunkerque, France</footnote>

<footnote>* Correspondence: medzaidi91@gmail.com</footnote>
biocomposite was developed using low-cost raw materials resulting from the flax industry. Non-woven flax fibers were chosen as a reinforcement, while mucilage polysaccharides, extracted from flaxseeds, served as a matrix. [15] investigated the possibility of using flax stems as reinforcement in a polylactic acid (PLA) matrix to produce a lightweight, fully bio-based composite with improved mechanical properties. Additionally, [16] provided an overview of flax fibers, the methods used to process these fibers, and the composites developed using different types of matrices. This overview aims to provide a succinct and fundamental resource for future research into flax-based composites.

Recently, [17,18] developed a new biocomposite-based flax fiber and flaxseed gum extracted from flax seeds as a matrix. This combination has exhibited interesting thermal and mechanical properties. An extensive thermal characterization of these biocomposite materials was conducted by considering various ratios of flaxseed gum, epoxy resin, and flax fibers. However, measuring the intrinsic thermal conductivity of individual components such as short flax fibers and flaxseed gum within the composite has proven to be challenging. Obtaining these thermal properties is essential, especially for modeling the thermal conductivity of the composites while varying their composition.

In the present study, a combined experimental and theoretical approach was pro-
posed to estimate the unknown thermal properties of these insulating biocomposites. This
approach is based on using analytical models that require the microstructural character-
istics for each component, including porosity, fiber volume fraction, and matrix volume
fraction [19,20].

In the literature, microstructure characterization of the biocomposites can be con-
ducted using X-ray Computed Tomography (X-ray CT) imaging and image analysis. This
technique enables the assessment of both quantitative and qualitative information in
three dimensions [21,22]. It was used to quantify the volume fraction, dimensions, and
orientation distribution of jute [23] and carbon fiber in composites [24]. To analyze poros-
ity, [25] investigated the size distribution and morphology of pore networks in carbon
fiber/epoxy composites using X-ray CT, while [26] analyzed the porosity and microstruc-
ture of natural coir fibers for their potential as reinforcements in composites. In the context
of improving braided composites, [12] studied the effect of voids in bio-based composites
formed from natural fibers and bio-resin. Additionally, image processing techniques were
used to calculate the size, pore size distribution, and shape of polycaprolactone-based
biocomposites [13].

The current study aimed to predict the effective thermal conductivity of a novel
flax-based biocomposite insulating material. The porosity and volume fraction of each com-
ponent were thoroughly obtained from image analysis. Such data are required for thermal
conductivity calculations. The originality of this study lies in the combined experimental
and analytical approach used to estimate the intrinsic conductivity of each component
of the material. Indeed, it is particularly difficult to measure them directly using typical
experimental devices. In addition, this study also aims to explore the potential applications
of these biocomposite as insulating materials in various applications.

## 2. Materials and Methods

Flaxseed gum polysaccharide (FG) used as a matrix for the composites was extracted according to the procedure described in [27]. The extraction process involved mixing 10% flax seeds with tap water and soaking for 1 h at 40 °C for 6 h with a stirring speed of 400 rpm. After filtering the seeds, the flax gum was precipitated with ethanol [8]. The epoxy resin (Epolam 2020) and its hardener were purchased from Sigma Aldrich, St. Louis, MO, USA. The flax fibers used for reinforcing composites were purchased from Van Robaeys Frère (Guillem, France) and were 1 mm in length. Biocomposite samples were prepared as follows: flaxseed gum (mucilage) was mixed with 160 mL of water and left to stand for 12 h. Then, epoxy resin and fibers (if used) were added to the polysaccharide solution to reach a mixture containing 20% organic matter. This mixture was poured into a silicone mold (44 mm in diameter and 2 mm in height) and frozen at -80 °C. Subsequently, the
samples were subjected to freeze drying for 72 h at 2 mbar (primary drying), followed by an additional 24 h at 0.002 mbar (secondary drying). Once the freeze-drying process was completed, the samples were baked in an oven at 80 °C for 5 h. Finally, the samples were polished to create cylindrical samples that were 20 mm in height with parallel top and bottom surfaces. Figure 1 shows the pictures of three samples and the images of their cross-section structure obtained with scanning electron microscopy (SEM).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Pictures and SEM images of the cross-section morphology of the three composites: FG20, FG80, and FFG.</figcaption></figure>

## 2.1. *Description of Materials*

Table 1 presents the composition, density, and measured thermal conductivity of the fabricated samples. A total of six samples were analyzed, with different mass ratio of flaxseed gum, epoxy resin, and flax fibers. The FG100, FG80, and FG20 samples do not contain flax fibers and are composed of 100%, 80%, and 20% by mass ratio of flaxseed gum, respectively, with the remaining proportion being epoxy resin. The fibrous flaxseed gum composites (FFG) contain a mixture of flaxseed gum, epoxy resin, and flax fibers. A dense sample of epoxy resin was also analyzed. The sample, referred to as “chopped fibers”, consists of compressed panels entirely made of random flax fibers measuring 1 mm in length.

<table><caption>**Table 1.** Mass fraction and thermo-physical properties of composite samples.</caption><thead><tr><th>Sample</th><th>m<sub>FG</sub><br>(%)</th><th>m<sub>Fibers</sub><br>(%)</th><th>m<sub>Epoxy</sub><br>(%)</th><th>Bulk Density<br>(g/cm<sup>3</sup>)</th><th>Thermal<br>Conductivity<br>(W·m<sup>-1</sup>·K<sup>-1</sup>)</th></tr></thead><tbody><tr><td>FG100</td><td>100.0</td><td>0.0</td><td>0.0</td><td>228.9 ± 9.3</td><td>0.054 ± 0.001</td></tr><tr><td>FG80</td><td>80.0</td><td>0.0</td><td>20.0</td><td>0.231 ± 8.2</td><td>0.065 ± 0.001</td></tr><tr><td>FG20</td><td>20.0</td><td>0.0</td><td>80.0</td><td>0.219 ± 3.0</td><td>0.057 ± 0.001</td></tr><tr><td>FFG</td><td>12.0</td><td>48.0</td><td>40.0</td><td>0.194 ± 5.3</td><td>0.064 ± 0.001</td></tr><tr><td>Epoxy resin (dense)</td><td>0.0</td><td>0.0</td><td>100</td><td>1.1</td><td>0.782 ± 0.001</td></tr><tr><td>Chopped fibers (1 mm<br>in length)</td><td>-</td><td>100.00</td><td>-</td><td>0.108</td><td>0.048 ± 0.001</td></tr></tbody></table>
Bulk density measurements (Table 1) were performed using an Electronic Digital LCD Gauge Stainless Nonius Caliper Micrometer 150 mm 6-inch MIS and a VWR LP 3102 balance. Composite samples were polished with a Struers Tegrapol 21 polisher and 80 g/cm² sandpaper to obtain parallel surfaces and a height of 20 ± 0.5 mm. For this, the volume of the samples was measured with calipers (4 samples per test) and weighed to an accuracy of 0.01 g.

Thermal conductivity measurements of materials were conducted using a C-Therm Technologies TCi Analyzer. The system is based on a modified transient plane source (MTPS) technique. It is important to note that all measurements were performed at ambient temperature. The experimental procedure conforms to ASTM D7984 [28].

## 2.2. Microstructural Characterization of Samples

The samples underwent characterization using X-ray CT and image analysis to analyze their microstructure and quantify the porosity and volume fraction of each component. The X-ray CT scans were conducted at INSA Lyon, France, using a laboratory X-ray scanner. The data acquisition system recorded a total of 1120 projections, evenly distributed over a 360° rotation along the vertical axis of the sample. To obtain the 3D reconstructions, the recorded projections were processed using Xact software (http://xactsoftware.co.in/ accessed on 5 September 2023).

To optimize both the quality of the images and the number of scans, the samples were scanned in smaller representative elementary volumes (REVs) of approximately 1 cm³. For each sample, a total of four REVs scans were randomly performed, with 8 µm in voxel size. It is important to note that during the image analysis, pores, and flax fibers in contact with the edges of the REV were excluded from the calculation.

Segmentation of the images was conducted using the LABKIT plug-in [29], which is an integrated feature within the ImageJ software [30]. An example of the quality of segmentation is shown in Figure 2, where black voxels represent void spaces, and grey voxels correspond to the solid phase within the sample.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Example of FG20 sample segmentation using LABKIT plug-in. (A) Original X-ray CT image, and (B) binary image.</figcaption></figure>

To determine the volume fraction of each component within the composite, we divided the volume of the segmented component by the total volume of the REV. This process was repeated for each REV, and the average volume fraction was then calculated.

## 2.3. Analytical Approach for Estimating Thermal Conductivity

This section focuses on the proposed theoretical approach based on three steps. Firstly, analytical models allow the estimation of the effective thermal conductivity of fibrous composites, including chopped flax fibers (in air). Then, analytical models are used to
determine the effective thermal conductivity of the porous matrix of composites. Finally, mixing laws are used to predict the local thermal conductivity (non-porous) of the matrix, $\lambda_{FG(s)}$, with a given FG100/epoxy volume fraction ratio.

### 2.3.1. Thermal Conductivity of the Fibrous Media

The analytical model of Schuetz and Glickman [31] allows for the prediction of the thermal conductivity of fibrous media. This model was previously used in similar types of bio-based composites containing juncus maritimus fiber media [32]. The effective thermal conductivity is determined using the following expression:

$$ \lambda_{eff} = \lambda_{FG}(1 - \nu_{Fb}) + \nu_{Fb} \left( \frac{2 - f_s}{3} \right) \lambda_{Fb(s)} \quad (1) $$

with $1 - \nu_{Fb} = \varphi + \nu_{FG(s)}$,

where,

$\lambda_{eff}$: is the effective thermal conductivity of the fibrous composite.

$\lambda_{FG}$: is the effective thermal conductivity of the porous flaxseed gum (FG) matrix.

$\lambda_{Fb(s)}$: is the local thermal conductivity of the flax fibers (s = solid).

$\nu_{Fb}$: is the volume fraction of the flax fibers.

$\varphi$: is the porosity.

$f_s$: is a morphological parameter ($f_s = 1$ for fibers).

It can be noted that this model can be used in the case of fibrous composites in air [26].

In this case, $\lambda_{FG}$ is replaced by $\lambda_g$:

$$ \lambda_{eff} = \lambda_g(1 - \nu_{Fb}) + \nu_{Fb} \left( \frac{2 - f_s}{3} \right) \lambda_{Fb(s)} \quad (2) $$

with $1 - \nu_{Fb} = \varphi$.

where,

$\lambda_g$: is the thermal conductivity of the air.

### 2.3.2. Thermal Conductivity of the Porous Flaxseed Gum Matrix

As presented in the materials and methods section, the non-fibrous composites are made from flaxseed-gum-filled epoxy resin. To determine the conductivity of the non-porous phase in the matrix $\lambda_{FG(s)}$, three analytical models were selected: Russell, Maxwell, and Bruggeman [33,34]. These models were chosen based on their capability for estimating thermal conductivity for different shapes of pores within the composite material.

Russell proposed a mathematical model (3) to predict the thermal conductivity in composites composed of cubic cells arranged in a row, i.e., series and parallel arrangements:

$$ \lambda_{FG} = \frac{\lambda_{FG(s)} \left[ \lambda_{FG(s)} + \varphi^{\frac{2}{3}} (\lambda_g - \lambda_{FG(s)}) \right]}{\lambda_{FG(s)} + (\lambda_g - \lambda_{FG(s)}) \left( \varphi^{\frac{2}{3}} - \varphi \right)} \quad (3) $$

where the components are as follows:

$\lambda_{FG}$: is the effective thermal conductivity of the porous flaxseed gum matrix.

$\lambda_{FG(s)}$: is the local thermal conductivity of the non-porous flaxseed gum matrix.

$\varphi$: is the porosity of the flaxseed gum matrix.

The Maxwell model, also known as the Maxwell–Eucken model, is an effective medium approximation used to estimate the thermal conductivity of porous materials. This
model (4) assumes that the fluid phase is composed of randomly dispersed spheres [33]. The equation of the effective thermal conductivity of a porous matrix can be expressed as follows:

$$ \lambda_{FG} = \lambda_{FG(s)} \frac{\lambda_g + 2\lambda_{FG(s)} + 2\varphi(\lambda_g - \lambda_{FG(s)})}{\lambda_g + 2\lambda_{FG(s)} - \varphi(\lambda_g - \lambda_{FG(s)})} \quad (4) $$

The Bruggeman model, also referred to as the Bruggeman homogenization model, is used to predict the effective thermal conductivity of porous materials that contain spherical or cylindrical inclusions [35]. In this study, the pores can be assumed to be nearly spherical, and the expressions for the effective thermal conductivity can be expressed as follows (5):

$$ \lambda_{FG} = \frac{\lambda_{FG(s)} \left[ 1 - \left( 1 - \frac{\lambda_g}{\lambda_{FG(s)}} \right)^{\frac{2}{3}} \varphi \delta \right]}{\left[ 1 + (\delta - 1) \varphi \right]} \quad (5) $$

where $\delta$ is determined from the Equation (6) for spherical particles:

$$ \delta = \frac{3\lambda_{FG(s)}}{2\lambda_{FG(s)} + \lambda_g} \quad (6) $$

### 2.3.3. Local Thermal Conductivity of the Flaxseed Gum Matrix

To estimate the local thermal conductivity of the non-porous matrix, $\lambda_{FG(s)}$, with a given volume fraction ratio of FG and epoxy resin, the analytical models require determination of the intrinsic thermal conductivity of FG and epoxy resin in the matrix. The thermal conductivity of the solid phase can be calculated using the mixing law of series (7) or parallel (8) mixing laws:

$$ \frac{1}{\lambda_{FG(s)}} = \frac{X}{\lambda_{FG100(s)}} + \frac{(1-X)}{\lambda_{Epoxy(s)}} \quad (7) $$

$$ \lambda_{FG(s)} = X\lambda_{FG100(s)} + (1-X)\lambda_{Epoxy(s)} \quad (8) $$

where,

X: is the volume fraction of the pure flaxseed gum (FG100) in the solid matrix.

$\lambda_{FG100(s)}$: is the intrinsic thermal conductivity of the non-porous FG100.

$\lambda_{Epoxy(s)}$: is the intrinsic thermal conductivity of the non-porous epoxy resin.

## 3. Results and Discussion

### 3.1. Microstructure Analysis

Figure 3 shows the microstructure of the samples obtained from X-ray CT images with a voxel size of 8 µm. The porosity network within flaxseed-gum-filled epoxy composites, encompassing different ratios of flaxseed gum/epoxy, highlights the presence of three distinct types of pores: (i) inter-lamellar pores separating two adjacent lamellae, with an average distance ranging from 2 to 300 µm; (ii) spherical and sub-spherical globular pores resulting from air trapped in the matrix during the fabrication process, which have an average maximum diameter of 1500 µm; and (iii) microcracks that developed due to the thermal stresses during the freezing process, with an orientation identical to that of the inter-lamellar pores.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** 2D visualization of the sample's microstructure and pores shape.</figcaption></figure>

The parallel lamellar alignment structure showed in samples FG20, FG802 and FG100 is a result of the extraction of lyophilized ice crystals in the vertical direction of the specimens. This similar lamellar structure was observed in previous studies involving graphene/epoxy and ceramic composites [36,37]. The architecture of FG20 displays a well-organized arrangement of lamellae, while FG80 and FG100 show some defects and microcracks. Reducing the epoxy content from 80% to 20% leads to a trend towards a less structured microstructure. In the fibrous composite (FFG), the incorporation of flax fibers into the flaxseed gum matrix resulted in substantial alterations to its microstructure. The images also revealed more spherical pores compared to the FG20 and FG80 samples, indicating good compatibility between the fibers and the matrix. Additionally, the absence of cracks and deformations in the fibrous composite suggests good structural cohesion.

Table 2 summarizes the porosity and volume fraction of the flaxseed-gum-filled epoxy and fibers in the samples obtained from the segmentation method. The porosity values range from ~61% ± 3.0 for FG20 to ~68% ± 2.4 for FG80 and FG100. This increase in porosity was attributed to alterations in the microstructure, such as the opening of interlamellar pores (Figure 3), as well as the development of some microcracks in FG80 and FG100 due to the increase in the volume fraction of the flaxseed gum in the matrix. These microcracks create additional pathways for fluid flow within the matrix, leading to an overall higher porosity. These results also highlight the impact of the mucilage/epoxy ratio on the porosity within the matrix. In addition, the panel made from chopped fibers (in air) has a higher porosity (76% ± 1.0) due to the absence of a fiber-binding component.

<table><caption>**Table 2.** Porosity and volume fraction of the flaxseed-gum-filled epoxy and fibers in the samples.</caption><thead><tr><th>Sample</th><th>ν<sub>FG</sub> (%)</th><th>ν<sub>fibers</sub> (%)</th><th>Porosity (%)</th></tr></thead><tbody><tr><td>FG100</td><td>32.0</td><td>-</td><td>68.0 ± 2.4</td></tr><tr><td>FG80</td><td>32.0</td><td>-</td><td>68.0 ± 2.5</td></tr><tr><td>FG20</td><td>39.0</td><td>-</td><td>61.0 ± 3.0</td></tr><tr><td>FFG</td><td>19.0</td><td>16.0</td><td>65.0 ± 1.1</td></tr><tr><td>Chopped fibers (1 mm)</td><td>-</td><td>24.0</td><td>76.0 ± 1.0</td></tr></tbody></table>
## 3.2. *Thermal Conductivity Estimation*

To determine the effective thermal conductivity of heterogeneous materials using analytical models, it is necessary to know the thermal conductivity of each component constituting the composites. As these are porous materials, it is essential to determine the conductivity of the non-porous phase of the matrix, as well as that of the flax fibers. Furthermore, as the matrix itself is made up of pure flaxseed gum (FG100) and epoxy resin, it is crucial to identify the conductivity of each of these components separately. The following sections discuss in detail the method used to determine the thermal conductivity of each component of the composite, using analytical models and mixing laws.

### 3.2.1. Thermal Conductivity of the Non-Porous Flaxseed-Gum-Filled Epoxy Matrix

As shown in Figure 3, the segmented images of the non-fibrous composites revealed two clearly identifiable phases: a continuous solid phase composed of flaxseed-gum-filled epoxy and a dispersed phase corresponding to the pores. However, since flaxseed gum and epoxy resin have a close attenuation coefficient, the images showed only one component, making it difficult to distinguish between the two components. Thus, both flaxseed gum and epoxy were considered as a one solid component. To determine the local thermal conductivity of this non-porous phase, the three analytical models Russell (3), Maxwell (4), and Bruggeman (5) (Cf. Section 3.2), based on different pore shapes, were used. As shown in Figure 4, the histograms present the variation in the local thermal conductivity of the two heterogeneous samples FG20 and FG80. These results indicate minimal discrepancies among the conductivities obtained from the three models. This suggests that the shape of the pores has limited influence on the thermal conductivity variation. The porosity has a more significant impact on conductivity than pore shape. Thus, in the next section, the usual Maxwell model is retained for its compatibility with the spherical shape of the pores present in the matrix, compared with the other tested models. The conductivity of the non-porous samples FG80 and FG20 are 0.181 and 0.123 W·m⁻¹·K⁻¹, respectively.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Estimated local thermal conductivity of the non-porous flaxseed gum matrix.</figcaption></figure>

### 3.2.2. Local Thermal Conductivity of Chopped Flax Fibers

The thermal conductivity measurements obtained from panels of chopped fibers (1 mm in length) were used to estimate the intrinsic thermal conductivity of individual flax fibers (without porosity). The analytical model of Schuetz and Glickman (2) was applied to a panel containing randomly oriented fibers in air (Figure 5B). This model considers the thermal conductivity of the solid fibers $\lambda_{Fb(s)}$, the porosity, and the morphological parameter $f_s$ of the solid phase. Considering the experimental thermal conductivity
(0.048 W·m⁻¹·K⁻¹) and the porosity (76% ± 1.0) obtained from X-ray CT images, the local thermal conductivity of the fibers, λ<sub>Fb(s)</sub>, was determined using the following expression:

$$ \lambda_{Fb(s)} = \frac{3}{2 - \text{fs}} \frac{\lambda_{Fb} - \varphi \lambda_g}{1 - \varphi} = 0.35 (\text{W} \cdot \text{m}^{-1} \cdot \text{K}^{-1}) \quad (9) $$

where,

λ<sub>Fb</sub>: is the effective thermal conductivity of the flax fibers (in air).

λ<sub>g</sub>: is the thermal conductivity of the air (λ<sub>g</sub> = 0.026 W·m⁻¹·K⁻¹).

<figure><img src="image_8.png" /><figcaption>Figure 5. X-ray CT images of the segmented flax fibers from (A) the fibrous composite (FFG) and (B) the panel's fibers (in air).</figcaption></figure>

In the literature and to our knowledge, there are no previous studies available in the literature that specifically provide the local thermal conductivity of randomly oriented short flax fibers. However, for well-oriented flax fibers, a study conducted by [38] provided results of experimental measurements of the thermal conductivity of flax fibers. In the longitudinal direction, the thermal conductivity was 1.232 (W·m⁻¹·K⁻¹), while in the transverse direction, it was 0.17 (W·m⁻¹·K⁻¹). These measurements diverge from our results because we used shorter and randomly oriented fibers, compared to their well-oriented fibers.

Note that in the Section 3.2.4, a second method for determining fiber conductivity based on the experimental conductivity of the fibrous composite will be proposed. The results obtained from both methods will then be compared and discussed.

### 3.2.3. Intrinsic Thermal Conductivity of the Pure Flaxseed Gum (FG100) and Epoxy Resin

To estimate the matrix's thermal conductivity for a given flaxseed gum/epoxy ratio, it is necessary to identify the intrinsic conductivity of the pure flax gum (FG100) and the epoxy resin present in the matrix. For this purpose, the experimental conductivities of the three samples, FG100, FG80, and the epoxy resin, were used. The effective conductivity of their solid phase (non-porous) was deduced from Maxwell's model by knowing the porosity (Table 2). Then, knowing the volume fractions of the pure flaxseed gum (FG100) and epoxy in the matrix, their local conductivity (FG100 and epoxy) was deduced using the series and parallel mixing laws (Equations (7) and (8)).

The volume fraction of epoxy resin v<sub>Epoxy(s)</sub> in each sample was determined using Equation (10) [32], where m<sub>Epoxy(s)</sub> represents the mass fraction of the epoxy and ρ<sub>matrix(s)</sub>/ρ<sub>Epoxy(s)</sub> is the ratio of the density of the matrix to that of the epoxy (Table 1).

$$ v_{Epoxy(s)} = m_{Epoxy(s)} \frac{\rho_{matrix(s)}}{\rho_{Epoxy(s)}} \quad (10) $$
where,

$$ \rho_{matrix(s)} = \frac{\rho_{matrix}}{1 - \varphi} $$

with

$$ \nu_{FG(s)} = 1 - \nu_{Epoxy(s)} $$

By knowing the volume fractions (Table 3) and the conductivity of the non-porous fraction of the three samples, FG100, FG80, and the epoxy resin (Table 1), the thermal conductivity of FG100 and epoxy resin can be obtained from mixing laws. From the linear trend curves of Equations (11) and (12), the slope and y-intercept can be deduced. This yields $\lambda_{FG100}$ and $\lambda_{Epoxy}$ (Figure 6).

$$ \frac{1}{\lambda_{FG(s)}} = X \left( \frac{1}{\lambda_{FG100(s)}} - \frac{1}{\lambda_{Epoxy(s)}} \right) + \frac{1}{\lambda_{Epoxy(s)}} \quad (11) $$

$$ \lambda_{FG(s)} = X(\lambda_{FG100(s)} - \lambda_{Epoxy(s)}) + \lambda_{FG100(s)} \quad (12) $$

<table><caption>Table 3. Volume fraction of FG100 and epoxy resin in the matrix.</caption><thead><tr><th>Sample</th><th>ν<sub>FG100</sub> (%)</th><th>ν<sub>Epoxy</sub> (%)</th></tr></thead><tbody><tr><td>FG100</td><td>100.00</td><td>0.00</td></tr><tr><td>FG80</td><td>80.30</td><td>19.70</td></tr><tr><td>FG20</td><td>59.20</td><td>40.80</td></tr><tr><td>FFG</td><td>35.00</td><td>65.00</td></tr><tr><td>Epoxy resin</td><td>0.00</td><td>100.00</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>Figure 6. The linear trend curves obtained from (A) parallel (Equation (11)) and (B) series mixing laws (Equation (12)).</figcaption></figure>

For FG100 sample, the thermal conductivities ($\lambda_{FG100}$) were around 0.056 and 0.057 (W·m⁻¹·K⁻¹) for both parallel and series models, respectively. Thus, an average value of 0.0565 W·m⁻¹·K⁻¹ can be considered for FG100. The thermal conductivity of epoxy resin is around 0.28 and 0.288 (W·m⁻¹·K⁻¹) for parallel and series models, respectively. Thus, a mean value of 0.284 W·m⁻¹·K⁻¹ can be considered for epoxy. It can be noted that in the literature, the thermal conductivity of epoxy varies between 0.18 and 0.26 (W·m⁻¹·K⁻¹) [39,40]. These values are of the same order of magnitude as that obtained and tend to show the validity of the approach.

### 3.2.4. Local Thermal Conductivity of Flax Fibers Based on the Effective Conductivity of the Fibrous Composite

In Section 3.2.2, the local thermal conductivity of the flax fibers was estimated from the experimental conductivity of the chopped fibers in air. An alternative method is also proposed for determining the local conductivity of fibers from the fibrous composite (FFG) (Figure 5A).
Given the known thermal conductivities of FG100 and epoxy resin, and their respective volume fractions (Table 3), the thermal conductivity of the non-porous FFG matrix was deduced: $\lambda_{FG(s)}$ was around 0.180 W·m⁻¹·K⁻¹. Considering that the porosity of FFG is 65% (Table 2) and by applying Maxwell’s model, the estimated thermal conductivity of the porous matrix, $\lambda_{FG}$, was 0.069 W·m⁻¹·K⁻¹.

Using the fiber volume fraction obtained through image analysis, the conductivity of the porous matrix (0.069 W·m⁻¹·K⁻¹) and the experimental effective conductivity of the FFG ($\lambda_{eff}$), the local conductivity of the fibers, $\lambda_{Fb(s)}$, was deduced from Equation (13). The obtained value was approximately 0.123 (W·m⁻¹·K⁻¹).

$$ \lambda_{Fb(s)} = \frac{3}{2-fs} \frac{\lambda_{eff} - (1 - \nu_{Fb})\lambda_{FG}}{\nu_{Fb}} = 0.123(W \cdot m^{-1} \cdot K^{-1}) \quad (13) $$

Note that this value is lower than the $\lambda_{Fb(s)} = 0.35 (W \cdot m^{-1} \cdot K^{-1})$ obtained in Section 3.2.2. In [32,41], this overestimation of the thermal conductivity of fibers in air was explained by the potential contribution of radiation inside highly porous Juncus maritimus fibers. Accordingly, the obtained value (0.123 W·m⁻¹·K⁻¹) is close to the average thermal conductivity of treated flax fibers 0.1187 (W·m⁻¹·K⁻¹), obtained by [19].

### 3.2.5. Estimating the Thermal Conductivity of the Porous Matrix for Different Volume Fractions of Flaxseed Gum

Figure 7 illustrates the evolution of flaxseed gum conductivity as a function of FG100 volume ratio. The conductivities of FG100 and epoxy were obtained by fitting linear trend curves using the two mixing laws in series and parallel, as discussed in Section 3.2.3. This curve represents all experimental and theoretical data points concerning matrix thermal conductivity in both fibrous and non-fibrous composites. A good agreement between experimental and theoretical results was observed. Since the parallel and series mixing laws yield values close to each other, the average value between the two models could be considered as a reliable estimation of the thermal conductivity of the composite matrix.

<figure><img src="image_10.png" /><figcaption>Figure 7. Comparison of flaxseed gum matrix thermal conductivity to experimental data.</figcaption></figure>

## 4. Conclusions

The present study aims to estimate the thermal conductivity of newly developed insulating composites made from flax fibers and a flaxseed-gum-filled epoxy matrix. An analytical approach was applied to estimate the local thermal conductivity of the matrix
and fibers using analytical models. In addition, the microstructure of the composites was characterized using X-ray CT and image analysis.

The composite's microstructure shows a lamellar structure with three types of pores: inter-lamellar, spherical and microcracks. The parallel alignment of the lamellae is due to the freeze-drying process. The lamellar structure of the flaxseed gum matrix, containing 80% epoxy (FG20), appears more organized and uniform than FG20 and FG100. The incorporation of reinforced fibers in the matrix (FFG sample) results in a structural modification that reinforces the composite's cohesion.

To overcome the lack of experimental data on the local thermal conductivity of the solid (non-porous) phase in composites, analytical models based on the two-phase continuous/discontinuous principle were used in this study. The results revealed that the shape of the pores has a negligible effect on the conductivity variation, while the porosity significantly influences the thermal conductivity of the composite material.

Since the matrix is made up of two homogeneous components, FG100 and epoxy resin, determining the conductivity of each component is crucial to assessing the overall conductivity of the matrix. Consequently, the average conductivity of pure flaxseed gum FG100 (0.145 W·m⁻¹·K⁻¹) and pure epoxy resin (0.284 W·m⁻¹·K⁻¹) was estimated using both parallel and series mixing laws.

The local thermal conductivity of flax fibers was determined from chopped flax fibers (in air) 1 mm long. Glickman's analytical model was used for this calculation, revealing a thermal conductivity of 0.35 W·m⁻¹·K⁻¹. However, this value may be overestimated due to the radiation contribution. An alternative estimation was therefore proposed, based on the fibrous composite (FFG). As expected, the fiber's thermal conductivity obtained was around 0.123 W·m⁻¹·K⁻¹. This value is close to the average thermal conductivity of the pure flaxseed gum, FG100 (0.145 W·m⁻¹·K⁻¹).

Obtaining the thermal conductivity of this composite material can be challenging, as it is difficult to directly measure the intrinsic thermal conductivity of each (non-porous) component using conventional experimental devices. However, determining the intrinsic conductivity of each component of the composite material could provide the basis for more advanced numerical modeling of the composite material in the future.

Finally, the thermal properties of these materials highlight their ability to provide effective thermal insulation, allowing the development of lightweight biocomposites capable of competing with synthetic foams commonly used in the automotive and building insulation fields. The renewable nature, natural origin, and biodegradability of these materials make them attractive for the development of competitive products in these sectors.

**Author Contributions:** Conceptualization, M.Z. and D.B.; methodology, M.Z. and D.B.; software, M.Z.; validation, M.Z., D.B. and F.D.; formal analysis, M.Z. and D.B.; investigation, M.Z. and D.B.; resources, F.D. and M.D.; data curation, M.Z.; writing—original draft preparation, M.Z. and D.B.; writing—review and editing, F.D and N.N.; visualization, N.N. and M.D.; supervision, F.D.; project administration, F.D. and D.B.; funding acquisition, F.D. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Region Hauts-de-France and Van Robaeys Frères.

**Institutional Review Board Statement:** Not applicable.

**Data Availability Statement:** No data available.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Philp, J.C.; Bartsev, A.; Ritchie, R.J.; Baucher, M.-A.; Guy, K. Bioplastics science from a policy vantage point. *New Biotechnol.* **2013**, *30*, 635-646. [CrossRef]

2. Phiri, R.; Mavinkere Rangappa, S.; Siengchin, S.; Oladijo, O.P.; Dhakal, H.N. Development of sustainable biopolymer-based composites for lightweight applications from agricultural waste biomass: A review. *Adv. Ind. Eng. Polym. Res.* **2023**, *in press*. [CrossRef]
3. Akter, M.; Uddin, M.H.; Tania, I.S. Biocomposites based on natural fibers and polymers: A review on properties and potential applications. *J. Reinf. Plast. Compos.* **2022**, *41*, 705–742. [CrossRef]

4. Andrew, J.J.; Dhakal, H.N. Sustainable biobased composites for advanced applications: Recent trends and future opportunities—A critical review. *Compos. Part C Open Access* **2022**, *7*, 100220. [CrossRef]

5. Yan, L.; Chouw, N.; Jayaraman, K. Flax fibre and its composites—A review. *Compos. Part B Eng.* **2014**, *56*, 296–317. [CrossRef]

6. Raja, R.A.A.; Sunil, J. Estimation of Thermal Conductivity of Nanofluids Using Theoretical Correlations. *Int. J. Appl. Eng. Res.* **2018**, *13*, 7932–7936.

7. Audenaert, A.; De Cleyn, S.H.; Buyle, M. LCA of low-energy flats using the Eco-indicator 99 method: Impact of insulation materials. *Energy Build.* **2012**, *47*, 68–73. [CrossRef]

8. Braulio-Gonzalo, M.; Bovea, M.D. Environmental and cost performance of building's envelope insulation materials to reduce energy demand: Thickness optimisation. *Energy Build.* **2017**, *150*, 527–545. [CrossRef]

9. Moudood, A.; Rahman, A.; Öchsner, A.; Islam, M.; Francucci, G. Flax fiber and its composites: An overview of water and moisture absorption impact on their performance. *J. Reinf. Plast. Compos.* **2019**, *38*, 323–339. [CrossRef]

10. Baley, C.; Bourmaud, A.; Davies, P. Eighty years of composites reinforced by flax fibres: A historical review. *Compos. Part A Appl. Sci. Manuf.* **2021**, *144*, 106333. [CrossRef]

11. Baley, C.; Gomina, M.; Breard, J.; Bourmaud, A.; Drapier, S.; Ferreira, M.; Le Duigou, A.; Liotier, P.J.; Ouagne, P.; Soulat, D.; et al. Specific features of flax fibres used to manufacture composite materials. *Int. J. Mater. Form.* **2019**, *12*, 1023–1052. [CrossRef]

12. Melenka, G.W.; Bruni-Bossio, B.M.; Ayranci, C.; Carey, J.P. Examination of voids and geometry of bio-based braided composite structures. *IOP Conf. Ser. Mater. Sci. Eng.* **2018**, *406*, 012012. [CrossRef]

13. Wang, Y.; Chang, H.-I.; Wertheim, D.F.; Jones, A.S.; Jackson, C.; Coombes, A.G.A. Characterisation of the macroporosity of polycaprolactone-based biocomposites and release kinetics for drug delivery. *Biomaterials* **2007**, *28*, 4619–4627. [CrossRef]

14. Alix, S.; Marais, S.; Morvan, C.; Lebrun, L. Biocomposite materials from flax plants: Preparation and properties. *Compos. Part A Appl. Sci. Manuf.* **2008**, *39*, 1793–1801. [CrossRef]

15. Verstraete, S.; Buffel, B.; Madhav, D.; Debruyne, S.; Desplentere, F. Short Flax Fibres and Shives as Reinforcements in Bio Composites: A Numerical and Experimental Study on the Mechanical Properties. *Polymers* **2023**, *15*, 2239. [CrossRef] [PubMed]

16. Zhu, J.; Zhu, H.; Njuguna, J.; Abhyankar, H. Recent development of flax fibres and their reinforced composites based on different polymeric matrices. *Materials* **2013**, *6*, 5171. [CrossRef]

17. Musa, C. Élaboration et Caractérisation de Matériaux Composites Biosourcés à Base de Mucilage et de Fibres de lin. Ph.D. Thesis, Université du Littoral Côte d'Opale, Côte d'Opale, France, 2019.

18. Musa, C.; Kervoëlen, A.; Danjou, P.-E.; Bourmaud, A.; Delattre, F. Bio-based unidirectional composite made of flax fibre and isosorbide-based epoxy resin. *Mater. Lett.* **2020**, *258*, 126818. [CrossRef]

19. Li, X.; Tabil, L.G.; Oguocha, I.N.; Panigrahi, S. Thermal diffusivity, thermal conductivity, and specific heat of flax fiber–HDPE biocomposites at processing temperatures. *Compos. Sci. Technol.* **2008**, *68*, 1753–1758. [CrossRef]

20. Nait-Ali, B.; Haberko, K.; Vesteghem, H.; Absi, J.; Smith, D.S. Thermal conductivity of highly porous zirconia. *J. Eur. Ceram. Soc.* **2006**, *26*, 3567–3574. [CrossRef]

21. Hanhan, I.; Agyei, R.; Xiao, X.; Sangid, M.D. Comparing non-destructive 3D X-ray computed tomography with destructive optical microscopy for microstructural characterization of fiber reinforced composites. *Compos. Sci. Technol.* **2019**, *184*, 107843. [CrossRef]

22. Panerai, F.; Ferguson, J.C.; Lachaud, J.; Martin, A.; Gasch, M.J.; Mansour, N.N. Micro-tomography based analysis of thermal conductivity, diffusivity and oxidation behavior of rigid and flexible fibrous insulators. *Int. J. Heat Mass Transf.* **2017**, *108*, 801–811. [CrossRef]

23. Jiang, N.; Li, Y.; Li, D.; Yu, T.; Li, Y.; Xu, J.; Li, N.; Marrow, T.J. 3D finite element modeling of water diffusion behavior of jute/PLA composite based on X-ray computed tomography. *Compos. Sci. Technol.* **2020**, *199*, 108313. [CrossRef]

24. Emerson, M.J.; Jespersen, K.M.; Dahl, A.B.; Conradsen, K.; Mikkelsen, L.P. Individual fibre segmentation from 3D X-ray computed tomography for characterising the fibre orientation in unidirectional composite materials. *Compos. Part A Appl. Sci. Manuf.* **2017**, *97*, 83–92. [CrossRef]

25. Yang, P.; Elhajjar, R. Porosity Content Evaluation in Carbon-Fiber/Epoxy Composites Using X-ray Computed Tomography. Available online: https://www.ingentaconnect.com/content/tandf/lpte/2014/00000053/0000003/art00001 (accessed on 15 July 2023).

26. Tran, L.Q.N.; Minh, T.N.; Fuentes, C.A.; Chi, T.T.; Van Vuure, A.W.; Verpoest, I. Investigation of microstructure and tensile properties of porous natural coir fibre for use in composite materials. *Ind. Crops Prod.* **2015**, *65*, 437–445. [CrossRef]

27. Dubois, F.; Musa, C.; Duponchel, B.; Tidahy, L.; Sécordel, X.; Mallard, I.; Delattre, F. Nuclear Magnetic Resonance and Calorimetric Investigations of Extraction Mode on Flaxseed Gum Composition. *Polymers* **2020**, *12*, 2654. [CrossRef] [PubMed]

28. D7984 Standard Test Method for Measurement of Thermal Effusivity of Fabrics Using a Modified Transient Plane Source (MTPS) Instrument. Available online: https://www.astm.org/d7984-21.html (accessed on 15 July 2023).

29. Arzt, M.; Deschamps, J.; Schmied, C.; Pietzsch, T.; Schmidt, D.; Tomancak, P.; Haase, R.; Jug, F. LABKIT: Labeling and Segmentation Toolkit for Big Image Data. *Front. Comput. Sci.* **2022**, *4*, 10. [CrossRef]

30. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. *Nat. Methods* **2012**, *9*, 671–675. [CrossRef]
31. Glicksman, L.R.; Torpey, M.; Marge, A. Means to Improve the Thermal Conductivity of Foam Insulation. *J. Cell. Plast.* **1992**, *28*, 571–583. [CrossRef]

32. Saghrouni, Z.; Baillis, D.; Naouar, N.; Blal, N.; Jemni, A. Thermal Properties of New Insulating Juncus Maritimus Fibrous Mortar Composites/Experimental Results and Analytical Laws. *Appl. Sci.* **2019**, *9*, 981. [CrossRef]

33. Maxwell, J.C. *A Treatise on Electricity and Magnetism*; Clarendon Press: Oxford, UK, 1873.

34. Pietrak, K.; Wiśniewski, T. A review of models for effective thermal conductivity of composite materials. *J. Power Technol.* **2015**, *95*, 14–24.

35. Lakhtakia, A.; Michel, B.; Weiglhofer, W.S. Bruggeman formalism for two models of uniaxial composite media: Dielectric properties. *Compos. Sci. Technol.* **1997**, *57*, 185–196. [CrossRef]

36. Liu, P.; Li, X.; Min, P.; Chang, X.; Shu, C.; Ding, Y.; Yu, Z.-Z. 3D Lamellar-Structured Graphene Aerogels for Thermal Interface Composites with High through-Plane Thermal Conductivity and Fracture Toughness. *Nano-Micro Lett.* **2020**, *13*, 22. [CrossRef] [PubMed]

37. Zhang, G.; Chen, H.; Yang, S.; Guo, Y.; Li, N.; Zhou, H.; Cao, Y. Frozen slurry-based laminated object manufacturing to fabricate porous ceramic with oriented lamellar structure. *J. Eur. Ceram. Soc.* **2018**, *38*, 4014–4019. [CrossRef]

38. Bodaghi, M.; Delfrari, D.; Lucas, M.; Senoussaoui, N.-L.; Koutsawa, Y.; Uğural, B.K.; Perrin, H. On the relationship of morphology evolution and thermal conductivity of flax reinforced polypropylene laminates. *Front. Mater.* **2023**, *10*, 1150180. [CrossRef]

39. Xu, S.; Liu, H.; Li, Q.; Mu, Q.; Wen, H. Influence of magnetic alignment and layered structure of BN&Fe/EP on thermal conducting performance. *J. Mater. Chem. C* **2016**, *4*, 872–878. [CrossRef]

40. Wan, X.; Demir, B.; An, M.; Walsh, T.R.; Yang, N. Thermal conductivities and mechanical properties of epoxy resin as a function of the degree of cross-linking. *Int. J. Heat Mass Transf.* **2021**, *180*, 121821. [CrossRef]

41. Borges dos Santos, M.; Moura, L.M.; Baillis, D. Identification of the Radiative Parameters-Albedo and Optical Thickness—Of the Juncus maritimus Fiber. *Materials* **2023**, *16*, 1891. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.