# Permeability Measurement of Glass-Fiber Textiles Used in Composites Industry Using Radial Flow Experimental Setup and Comparison with Image-Based Numerical Methods

Mouadh Boubaker, Willsen Wijaya, Arthur Cantarel, Gerald Debenest, Simon Bickerton

► To cite this version:

Mouadh Boubaker, Willsen Wijaya, Arthur Cantarel, Gerald Debenest, Simon Bickerton. Permeability Measurement of Glass-Fiber Textiles Used in Composites Industry Using Radial Flow Experimental Setup and Comparison with Image-Based Numerical Methods. Sci, 2024, 6 (3), pp.49. 10.3390/sci6030049. [hal-05048472]

HAL Id: hal-05048472

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

Submitted on 29 Apr 2025

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 Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Permeability Measurement of Glass-Fiber Textiles Used in Composites Industry Using Radial Flow Experimental Setup and Comparison with Image-Based Numerical Methods

Mouadh Boubaker<sup>1,2</sup>, Willsen Wijaya<sup>3</sup>, Arthur Cantarel<sup>1,4</sup>, Gérald Debenest<sup>5,6,*</sup> and Simon Bickerton<sup>3</sup>

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

check for
updates

**Citation:** Boubaker, M.; Wijaya, W.; Cantarel, A.; Debenest, G.; Bickerton, S. Permeability Measurement of Glass-Fiber Textiles Used in Composites Industry Using Radial Flow Experimental Setup and Comparison with Image-Based Numerical Methods. *Sci 2024*, *6*, 49.
https://doi.org/10.3390/sci6030049

Academic Editor: Cristiano Fragassa

Received: 2 June 2024

Revised: 18 August 2024

Accepted: 19 August 2024

Published: 26 August 2024

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

**Copyright:** © 2024 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/).

**Abstract:** Permeability measurement of engineering textiles is a key step in preparing composite manufacturing processes. A radial flow experimental setup was used in this work to measure the unsaturated and saturated in-plane permeabilities of five different types of E-glass textiles and their ratios. In parallel, delayed tow saturation during the oil injection stage was visually observed to identify fabrics that exhibited a significant dual-scale effect. A numerical approach to determine the saturated permeability of a given fabric geometry at the mesoscale was tested and validated against analytical models found in the literature. It was then applied to a realistic geometry acquired from an E-glass plain weave textile using an X-ray microtomography scanner (µCT). Two numerical methods were adopted: the single-scale method, where the tows are considered impermeable, and the dual-scale method, where the permeability of the tows is taken into account. The numerical results from both methods were then compared with the experimental values and showed good agreement, especially with the second method.

**Keywords:** composites; permeability; saturated-unsaturated flow; radial flow; fibrous textiles; image-based numerical methods; numerical modeling

## 1. Introduction

Composite materials are increasingly and widely used in automobile, shipbuilding, aeronautics, and aerospace applications. These materials are manufactured from fibrous reinforcements that can take the form of woven fabrics, non-crimp fabrics, knits, braids, or random mats [1]. The fabrics are transformed into composite parts using different manufacturing processes. One of the most used and developed methods is liquid composite molding (LCM) due to its various advantages, such as being low-cost, clean, and more efficient than other methods. The LCM process is based on the impregnation of dry reinforcement with a low-viscosity resin [2]. During the impregnation step, the liquid resin flows through the reinforcement, which can be considered a porous medium. This flow is a key phenomenon that has a significant effect on the quality of the final part since dry zones (not fully covered by resin) can appear if flow parameters such as the injection pressure and the positions of the injection and vent gates are not well defined [3]. An important parameter to be considered in this step is the permeability of the preform, which is defined

<footnote n="1">¹ Institut Clément Ader (ICA), CNRS UMR 5312, University of Toulouse, 65000 Tarbes, France; mouadh.boubaker@iut-tarbes.fr (M.B.); arthur.cantarel@iut-tarbes.fr (A.C.)</footnote>

<footnote n="2">² Institut Universitaire Technologique of Tarbes, Université Paul Sabatier UPS, 65000 Tarbes, France</footnote>

<footnote n="3">³ Centre for Advanced Materials Manufacturing and Design, Department of Mechanical Engineering, The University of Auckland, Auckland 1010, New Zealand; willsen.wijaya@auckland.ac.nz (W.W.); s.bickerton@auckland.ac.nz (S.B.)</footnote>

<footnote n="4">⁴ Institut Universitaire Technologique of Tarbes, Université de Technologie Tarbes Occitanie Pyrénées UTTOP, 65000 Tarbes, France</footnote>

<footnote n="5">5 Toulouse INP, Université Paul Sabatier, 31400 Toulouse, France</footnote>

<footnote n="6">6 IMFT (Institut de Mécanique des Fluides de Toulouse), Université de Toulouse, Allée Camille Soula, 31400 Toulouse, France</footnote>

<footnote>* Correspondence: gerald.debenest@toulouse-inp.fr</footnote>
as the ability of the porous domain to let fluid flow through it. Permeability influences the resin flow front advancement through the fibrous reinforcement and ultimately controls the mold-filling time. This material property is crucial for composite fabrication optimization and is an important input for mold-filling simulations [4]. These simulations are used to accurately predict the flow front position, fill time, and formation of voids and dry spots. A good characterization of the material properties is important to run them properly.

In order to meet the growing demand for composites, different types of reinforcements are being developed with different properties and architectures. Therefore, an accurate and efficient characterization method, such as permeability measurement, is needed. There exists a large variability in permeability measurement methods, as no standard method has been achieved [4]. Errors in permeability prediction can lead to improper process design and manufacturing of defective parts.

The flow of resin through a fibrous network is usually considered a creeping flow associated with low Reynolds number (Re) values. Under these conditions, Darcy's law (Equation (1)) can be used to model the resin flow by relating the volume-averaged flow velocity u to the volume-averaged pressure P through the reinforcement permeability tensor K and the fluid viscosity µ. This law, based on experimental observations, is a phenomenological model [5] used to characterize fluid flow in porous media at the macro scale.

$$u = -\frac{K}{\mu}\nabla P \qquad (1)$$

Darcy's law was applied in our case to model fluid flow in fibrous textiles. It was originally introduced to model the flow in sands and is a theoretical foundation of groundwater hydrology. It is a phenomenological law that has been used by various researchers in fluid mechanics and reservoir engineering in order to solve complex problems of fluid flow in porous media such as rocks [6]. Several studies have been conducted to review the use of Darcy's law to model the flow of natural reservoirs [7] and underground fluids [8].

This paper is aimed at characterizing the permeability K of a given textile using experimental and numerical methods. Different permeability determination methods can be found in the literature and can be classified into three main groups: analytical, experimental, and numerical methods. Different analytical models relate the permeability to the fiber volume fraction Vf and are based on simplifying assumptions where the fibers are modeled as parallel cylinders arranged in regularly ordered networks (e.g., quadratic and hexagonal) [9,10]. These models are idealized, valid only for a certain range of fiber volume fractions, and do not accurately represent the complex geometry of real preforms. Alternatively, many researchers have used experimental procedures to measure the permeability of fibrous reinforcements. These experiments can be classified according to the following parameters [11–13]:

* Type of fluid: gas or liquid.

* Direction of measurement: in-plane, out-of-plane, or 3D.

* Type of flow: radial or unidirectional.

* Injection boundary conditions: constant pressure or constant flow rate.

* Saturation state of the textile: saturated or unsaturated.

The experiments allow us to take into account the variabilities existing in the textiles, but they have the disadvantage of being labor- and, time-intensive and resource-consuming since they involve many sensors, transducers, and digital cameras [14]. Another drawback is the absence of a standardized approach. Several international benchmark exercises [15,16] were launched in order to study the variation in the experimental measurements on different geometries (non-crimp and woven fabrics) and found a coefficient of variation of more than 20%.

A promising alternative is the use of numerical methods that combine textile modeling and computational fluid dynamics CFD techniques. The potential of numerical methods is increasing with the growth in computing capabilities. Numerical permeability prediction
has many advantages, such as the ability to efficiently explore a large range of parameters without the need to resort to experimental measurements. There are, however, modeling challenges associated with the complex structure of reinforcements that include geometrical variability, anisotropic structure, and dual-scale porosity. The latter, in particular, necessitates the use of multiscale approaches in order to predict the virtual permeability accurately. The choice of the modeling scale (Figure 1) is important to optimize the numerical simulation since at the macroscale (i.e., component scale, m-range), the preformed geometrical details are not well represented. In contrast, at the micro-scale (fiber scale, µm range), the computational cost is prohibitively high [17]. The mesoscale models (tow scale, mm-range) provide a good compromise between accuracy and efficiency.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Modeling scales.</figcaption></figure>

In the first section of this paper, a 2D radial flow experimental procedure to determine the in-plane permeability tensor of different textiles is presented. The setup is designed to measure the unsaturated (transient) and saturated (steady-state) permeability in a single experimental run. The tests were carried out on different types of E-glass fabrics with various stitching styles and areal weights. The ratio of the unsaturated to saturated permeability for the different textiles was quantified at three different fiber volume fractions, and the flow front advancement was observed in order to identify the fabrics that exhibited a significant dual-scale effect (the difference in fluid velocity between the inter-tow and intra-tow space).

In the second section, two numerical methods for predicting the saturated permeability of a given textile geometry at the mesoscale are described. The first method is based on solving the Navier–Stokes equations in the inter-tow space while considering the tows as impermeable (i.e., single-scale flow). In the second method, the Brinkman equation is solved over the entire domain, and the flow inside the tows is considered (i.e., dual-scale flow). The numerical approach was first validated against analytical models found in the literature and then used to determine the permeability of a realistic geometry acquired from an E-glass plain weave fabric at three fiber volume fractions using an X-ray µCT scanner. The measured saturated permeabilities of the scanned textiles were compared with their corresponding predicted values. Good agreement was found, especially regarding the dual-scale case.

## 2. Radial Flow Experimental Setup for the Measurement of the Permeability of Fibrous Textiles

A 2D radial flow in-plane permeability experimental setup was used to measure the unsaturated and saturated permeability of five different E-glass fabrics at different fiber volume fractions: plain weave fabric, 0/90 biaxial stitched fabric, quad axial fabric, unidirectional reinforcement, and twill weave fabric. Examples are given in Figure 2.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** List of the used textiles.</figcaption></figure>

The test rig has been developed over the years in the Centre for Advanced Materials Manufacturing and Design (CAMMD) at The University of Auckland in order to efficiently characterize the in-plane permeability of fibrous textiles [18]. Figure 3 shows a schematic representation of this rig that was used by F. Zhang et al. [19] and S. Comas-Cardona [20] to investigate the transient permeability of Chopped Strand Mat (CSM) and other fabrics.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Schematic representation of the experimental facility [14].</figcaption></figure>

In the first part of this section, the experimental facility, including the different tools used and the steps of the experimental procedure, are presented. In the second part, the results are discussed, and the key findings of these experiments are presented.

In the scope of this paper, the test rig has been used to measure two different permeability values: the unsaturated or transient permeability through the following of the flow front, as will be explained later, and the saturated or steady-state permeability through the
injection of an additional quantity of oil after the saturation of sample which enables to compute the steady-state mass flow rate.

The radial flow technique was particularly chosen because of the relative simplicity of the setup compared to other techniques, in addition to its ability to compute the complete in-plane permeability tensor within a single test, which minimizes effort and cost.

## 2.1. *Experimental Facility*

The 2D in-plane testing facility was installed on an Instron 1186 Universal Testing Machine (UTM) to provide accurate control of the cavity thickness and measurement of the total compaction force applied. The whole facility was situated in an environmentally controlled room to obtain consistent measurements (a constant temperature of 22 °C).

Each sample was set between two plates: a transparent lower plate, Figure 4a, made of 55 mm thick glass with a side length of 350 mm that allows flow front detection through optical monitoring and a circular upper plate, Figure 4b, of a diameter of 250 mm and 30 mm thick that was parallelly aligned to the lower plate with a spherical alignment unit and a manufactured strain-gauged compressive alignment specimen. Through these alignment tools, a uniform cavity thickness with a tolerance of 5 µm could be achieved. These plates were designed to minimize deflection during the compaction step, allowing a maximum deflection of 0.07 mm under an average compaction pressure of 2 MPa [20]. The range of target thickness was around 3 mm which adheres to the ISO 4410:2023 [21] recommendations (Test methods for the experimental characterization of in-plane permeability of fibrous reinforcements for liquid composite molding), with target thickness precision being around 2%.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Upper and lower plates of the experimental facility.</figcaption></figure>

A Navitar SV5C10 0238 digital camera (manufactured by NAVITAR, INC in Rochester, NY, USA) was attached to the setup via an L-bend with its lens pointing to a 45° angled mirror (with respect to the horizontal plate) to capture the flow front evolution during the test fluid infusion. Two Mitutoyo IDU25 touch sensors 1 µm resolution were used to measure the distance between the parallel plates; they were attached to the top plate with their probe touching the upper surface of the bottom plate. National Instruments LabVIEW saved the captured data (time, pressure, and oil mass) for post-processing.

The different measurement accuracies of the experimental facility are summarized in the following table (Table 1).
<table><caption>**Table 1.** Measurement accuracies.</caption><thead><tr><td>Measurement</td><td>Accuracy</td></tr></thead><tbody><tr><td>Sample mass</td><td>±0.01 g</td></tr><tr><td>Injected oil mass</td><td>±0.1 g</td></tr><tr><td>Cavity thickness</td><td>±5 µm</td></tr><tr><td>Fluid pressure</td><td>±1%</td></tr><tr><td>Temperature</td><td>±0.1 °C</td></tr></tbody></table>

A general overview of the experimental setup is shown in Figure 5.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** 2D radial flow experimental setup.</figcaption></figure>

## 2.2. *Experimental Procedure*

Utilizing a hydraulic press machine and a dedicated circular cutting knife, Figure 6a, circular stacks with a diameter of 200 mm of fibrous textiles were prepared.

A hole of 15 mm diameter was punched at the center of each sample to serve as an injection port for the test fluid, Figure 6b. The mass of each sample was measured before starting the experiment using a digital balance (which measures up to 0.01 g precision). The sample was then placed on the bottom glass plate, Figure 7, and must be accurately placed at the center to ensure a pure in-plane flow during the injection step. The layers were assembled in a 0/90/0/90 sequence to enforce isotropic in-plane flows and, thus, circular flow fronts.
<figure><img src="image_10.png" /><figcaption>Figure 6. Sample preparation.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Figure 7. A sample of four circular layers placed on the lower plate of the experimental setup.</figcaption></figure>

The position of the bottom plate was controlled by the Instron machine in order to achieve the target fiber volume fraction, $V_f$, as given by Equation (2):

$$V_f = \frac{m}{\rho Ah} \qquad (2)$$

where $m$ is the mass of the sample, $\rho$ is the fiber density (characteristic of the used fabric), $A$ is the sample's planar area, and $h$ is the target compaction thickness. In this equation, the term $\frac{m}{\rho}$ provides the actual volume of the sample, while the term $A \cdot h$ is the total volume occupied by the sample plus air/oil; thus, their ratio provides the actual fiber volume fraction. The following table (Table 2) summarizes the average sample mass and the average cavity thickness for each fiber volume fraction for the different textiles:

The test fluid was then injected radially from the oil pot into the sample initially filled with air via a circular port of 15 mm in diameter, centered on the top plate. The test fluid used was a Newtonian Mobil DTE series HEAVY mineral oil ISO VG 100, whose temperature-dependent viscosity, Figure 8, was determined using a Parr Physica UDS200 rheometer and is given by the following power-law:

$$\mu = -0.00000000519 \cdot T^5 + 0.00000111 \cdot T^4 - 0.0000969 \cdot T^3 + 0.00448 \cdot T^2 - 0.114 \cdot T + 1.38 \quad (3)$$
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Viscosity of Mobil heavy mineral oil as a function of temperature.</figcaption></figure>

The oil was supplied from an enclosed metallic pot driven by pressurized air regulated by a FESTO proportional pressure controller. During the injection, the fluid flowed through the specimen radially, and the flow front progression was captured and analyzed to determine the in-plane permeability of the sample.

The oil was injected from the sample center at a constant pressure of 50 kPa, and the outlet was placed at atmospheric pressure. The flow front progression was recorded as a series of images at a frequency of 1 Hz, Figure 9.

<table><caption>**Table 2.** Average sample mass and cavity thickness for different fiber volume fractions.</caption><thead><tr><th>Textile</th><th>V<sub>f</sub></th><th>Average Sample Mass [g]</th><th>Average Cavity Thickness [mm]</th></tr></thead><tbody><tr><td rowspan="3">Plain Weave fabric</td><td>0.5</td><td>106.8</td><td>2.7</td></tr><tr><td>0.55</td><td>107.1</td><td>2.43</td></tr><tr><td>0.6</td><td>107</td><td>2.25</td></tr><tr><td rowspan="3">0/90 Biaxial stitched fabric</td><td>0.5</td><td>114.275</td><td>2.86</td></tr><tr><td>0.55</td><td>114.33</td><td>2.6</td></tr><tr><td>0.6</td><td>113.9</td><td>2.38</td></tr><tr><td rowspan="3">Quad Axial stitched fabric</td><td>0.5</td><td>109.3</td><td>2.74</td></tr><tr><td>0.55</td><td>109.33</td><td>2.49</td></tr><tr><td>0.6</td><td>109.6</td><td>2.29</td></tr><tr><td rowspan="3">Unidirectional reinforcement fabric</td><td>0.5</td><td>101.45</td><td>2.544</td></tr><tr><td>0.55</td><td>101.56</td><td>2.315</td></tr><tr><td>0.6</td><td>101.375</td><td>2.118</td></tr><tr><td rowspan="2">Twill Weave fabric</td><td>0.5</td><td>74.16</td><td>1.86</td></tr><tr><td>0.55</td><td>74.6</td><td>1.69</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>(a) Initial image</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>(b) Intermediate image</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>(c) Final image</figcaption></figure>

<caption>**Figure 9.** Flow front progression.</caption>
In addition to the flow front progression, the static inlet pressure measured by a pressure sensor located near the injection gate and the oil temperature measured by a thermocouple in the pressure pot were also recorded at regular time intervals.

A MATLAB 2023b script developed in CAMMD was used to process the captured flow front images [14]. The algorithm comprises the following steps starting from the initial image captured by the camera, which was converted to grayscale, Figure 10a:

* The possible distortions due to the camera and mirror positioning are corrected using a geometric transformation routine applied to the input images, Figure 10b.

* A binary image is obtained from the subtraction of the transformed image from its corresponding initial image, Figure 10c.

* The boundary pixels of the binary image are fitted with a direct least-square ellipse, from which various geometrical parameters are extracted: major radius a, minor radius b, and major axis orientation θ, Figure 10d.

<figure><img src="image_16.png" /><figcaption>Figure 10. Image processing algorithm steps: (a) Original grayscale image, (b) Image after geometric transformation, (c) Binary image after subtraction, (d) Image with ellipse fitting.</figcaption></figure>

The unsaturated permeability was computed using the formula given by Weitzenbock et al. [22], derived from Darcy's law for radial flow in equations (Equations (4) and (5)).

$$K_{11} = (a^2[2\ln(\frac{a}{r_i}) - 1] + r_i^2) \frac{\mu(1 - V_f)}{4tP_{inj}} \quad (4)$$

$$K_{22} = (b^2[2\ln(\frac{b}{r_i}) - 1] + r_i^2) \frac{\mu(1 - V_f)}{4tP_{inj}} \quad (5)$$

where $K_{11}$ and $K_{22}$ are the major and minor in-plane principal permeability components, a and b are the major and minor radii of the fitted flow front ellipse, $r_i$ is the radius of the injection hole, $\mu$ is the test oil viscosity deduced from the measured oil temperature using the viscosity-temperature curve (Figure 8), t is the time from the start of the injection, and $P_{inj}$ is the inlet injection pressure. In the Weitzenbock formula [22], the pressure used corresponds to the pressure drop, but it is approximated to the injection pressure $P_{inj}$ in our case, we assume that the outlet is at atmospheric pressure (i.e., zero gauge pressure).

For each test, the permeability tensor is computed for each flow front image and plotted against the time steps. The reported permeability values are given by the average of the final one-third of the computed values to eliminate any flow stabilization effects at the start of the experiment. As can be observed in Figure 11, the measured permeability values become stable after a short initial transient effect.

This setup also allows the measurement of isotropic saturated permeability. After the sample is fully saturated, the oil injection is continued in order to compute the steady-state mass flow rate, which is found by fitting a linear regression to the mass-time curve. The saturated permeability is deduced using Darcy's law through equation (Equation (6)) [23].
$$K_{\text{saturated}} = \frac{Q \cdot \mu \cdot \ln \left( \frac{r_0}{r_i} \right)}{2\pi \cdot \rho \cdot h \cdot P_{\text{inj}}} \quad (6)$$

where Q is the steady-state mass flow rate, μ is the oil viscosity, r₀ is the sample's outer radius, rᵢ is the internal radius taken as the radius of the punched hole, ρ is the oil density, h is the cavity thickness, and Pᵢₗᵢ is the injection pressure.

<figure><img src="image_17.png" /><figcaption>Figure 11. Output of the post-processing MATLAB algorithm.</figcaption></figure>

## 2.3. Results

These experiments were carried out in order to quantify the difference between the transient and steady-state permeabilities of the tested textiles and to observe the flow front advancement during the injection step. The flow front shape provides an idea of fabrics where the dual-scale effect is significant. Three different fiber volume fractions Vf (0.5, 0.55 and 0.6) were considered, except for the twill weave fabric, for which 16 layers were used for each sample because of its low linear density; as a result, the quantity of fabric was not sufficient to carry out the tests at three fiber volume fractions. The isotropic unsaturated permeability used for the comparison is the average (geometric mean) of the values found in the principal directions $K_{\text{unsat}} = \sqrt{K_{11} \cdot K_{22}}$.

The ratios $K_{\text{unsat}}/K_{\text{sat}}$ are different from one textile to another and were found to vary between 0.9 and 1.42, which is the range found in the literature. Researchers found ratios between 0.25 and 4 [2,24], and for most of them [25,26], the ratio was around 1. The difference in this ratio from one textile to another can be explained by the difference in the size of the inter-tow space, which affects the importance of the capillary effects and their effect on the transient flow and, as a result, on the value of the unsaturated permeability. Table 3 summarizes the experimental findings. Each test was repeated three to four times, and the relative standard deviation was between 3% and 16%, which is comparable to the standard deviations (around 15%) found in the experimental benchmark [16] and provides confidence in the repeatability of our results. The discrepancy between different repeats can be attributed to the manual preparation and handling of samples, uncertainties in sensor readings, and the effects of ambient conditions.

Specific to the E-glass plain weave textile (800 g/m²), Wijaya [14] performed detailed unsaturated and saturated permeability characterizations of the same fabric in both the warp and weft directions. In the study, no relative change in orientation within a given sample stack was permitted (i.e., the warp and weft tows of all layers were nominally aligned in their respective directions). The results of the said work served as a basis for comparison with our values since the isotropic value that we are measuring should lie
more or less between the weft and warp permeability given by his tests. The results plotted in the following Figure 12 confirm the good agreement between both measurements and provide more confidence in the accuracy and repeatability of our results.

<table><caption>**Table 3.** Summary of the experimental results.</caption><thead><tr><th>Textile</th><th>V<sub>f</sub></th><th>K<sub>sat</sub><br>[10<sup>-10</sup>m<sup>2</sup>]</th><th>K<sub>unsat</sub><br>[10<sup>-10</sup>m<sup>2</sup>]</th><th>K<sub>unsat</sub>/K<sub>sat</sub></th></tr></thead><tbody><tr><td rowspan="3">E-glass Plain weave (800 g/m<sup>2</sup>)</td><td>0.5</td><td>2.916 ± 0.19</td><td>3.16 ± 0.24</td><td>1.08</td></tr><tr><td>0.55</td><td>1.12 ± 0.11</td><td>1.16 ± 0.14</td><td>1.04</td></tr><tr><td>0.6</td><td>0.482 ± 0.06</td><td>0.458 ± 0.054</td><td>0.95</td></tr><tr><td rowspan="3">0/90 Biaxial stitched fabric (825 g/m<sup>2</sup>)</td><td>0.5</td><td>1.39 ± 0.2</td><td>1.6 ± 0.3</td><td>1.15</td></tr><tr><td>0.55</td><td>0.672 ± 0.097</td><td>0.729 ± 0.11</td><td>1.08</td></tr><tr><td>0.6</td><td>0.326 ± 0.044</td><td>0.34 ± 0.011</td><td>1.04</td></tr><tr><td rowspan="3">Quad Axial Stitched Fabric (851 g/m<sup>2</sup>)</td><td>0.5</td><td>0.62 ± 0.026</td><td>0.649 ± 0.038</td><td>1.045</td></tr><tr><td>0.55</td><td>0.27 ± 0.012</td><td>0.286 ± 0.007</td><td>1.05</td></tr><tr><td>0.6</td><td>0.171 ± 0.016</td><td>0.177 ± 0.0178</td><td>1.03</td></tr><tr><td rowspan="3">Unidirectional reinforcement fabrics (850 g/m<sup>2</sup>)</td><td>0.5</td><td>3.5 ± 0.22</td><td>4.97 ± 0.8</td><td>1.42</td></tr><tr><td>0.55</td><td>1.91 ± 0.057</td><td>2.47 ± 0.097</td><td>1.29</td></tr><tr><td>0.6</td><td>1.12 ± 0.16</td><td>1.33 ± 0.218</td><td>1.19</td></tr><tr><td rowspan="2">Twill weave (295 g/m<sup>2</sup>)</td><td>0.5</td><td>0.398 ± 0.004</td><td>0.38 ± 0.045</td><td>0.96</td></tr><tr><td>0.55</td><td>0.2 ± 0.022</td><td>0.183 ± 0.023</td><td>0.91</td></tr></tbody></table>

Another objective of the experiments was to observe the flow front advancement in the different fabrics in order to identify those where the dual-scale effect, i.e., the difference in flow front progress between the inter-tow and intra-tow space, is stronger. This effect is clearly visible in the unidirectional reinforcement (Figure 13a) and biaxial stitched fabrics (Figure 13b) and less so in the plain weave textile (Figure 13c), whereas the flow front was nearly compact and uniform for the twill weave and quad axial stitched fabrics (Figure 14).

<figure><img src="image_18.png" /><figcaption>**Figure 12.** *Cont.*</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Figure 12.** Comparison with the results obtained by W. Wijaya [14].</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Figure 13.** Oil flow fronts with visible dual-scale effect.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Figure 14.** Compact oil flow fronts.</figcaption></figure>

These observations are consistent with the architecture of the textiles since the free space or inter-tow space is larger for the unidirectional, biaxial, and plain weave textiles
than for the quad axial and twill weave textiles. The oil flowed faster in the inter-tow spaces than in the intra-tow spaces; such a difference in velocity was evident in these textiles.

### 3. Permeability Determination Using Numerical Simulation

Numerical simulation is largely used at the component scale to simulate the LCM process; however, its adoption for textile permeability characterization is not well developed despite its increasing potential. It appears as an attractive proposition for theoretical and experimental alternatives because it neither generates large amounts of material waste (test fluid, textiles, etc.) nor uses complex equipment, and it allows it to carry out many parametric studies efficiently [4].

In this section, two numerical methods for saturated permeability determination are presented. The first method, the single-scale approach, is based on solving the Navier-Stokes equations in the inter-tow free space and neglecting the flow inside the fiber tows (i.e., solid impermeable tows). This assumption is justified by the low permeability of the tows (around $10^{-13} \text{m}^2$), which is around two orders of magnitude lower than the overall textile permeability [27]. The second method, the dual-scale approach, is based on solving the Brinkman equation on the whole domain (including the intra-tow domain) using a specific multiscale OpenFoam solver that allows simulating the inter- and intra-tow flows by considering the computational domain as a porous medium.

Both methods were applied to a scanned 3D geometry acquired from a 4-layer stack E-glass plain weave specimen with an areal weight of $800 \text{ g.m}^{-2}$ and a tow linear density of $2400 \text{ g.km}^{-1}$. The numerical values are compared to the experimental values obtained with the experimental setup described in the previous section that was used to measure the weft and warp permeabilities of a sample of four layers of the same plain weave textile.

#### 3.1. Numerical Approach

The general procedure for numerically determining the permeability of a given flow geometry is as follows. A pressure gradient $\Delta P$ is imposed between the inlet and the outlet of a given domain across a flow distance $L_{flow}$. The steady-state volumetric flow rate $Q$ through a cross-sectional area normal to the flow direction $A_{flow}$ is numerically computed by solving the corresponding fluid equations assuming creeping flow (i.e., incompressible and non-inertial characterized by a low Reynolds number).

The flow domain is considered saturated by a fluid of viscosity $\mu$ so that the capillary effects are neglected [28], and the permeability is deduced from Darcy's law (Equation (1)) for linear flow according to the following relationship:

$$K = \frac{Q \cdot \mu \cdot L_{flow}}{A_{flow} \cdot \Delta P} \quad (7)$$

The flow through the structure is solved using the computational fluid dynamics CFD open-source software OpenFOAM v7 based on the Finite Volume Method FVM.

#### 3.2. Comparison with Analytical Models

The first test to validate the numerical approach was performed on a square array of impermeable cylinders, Figure 15a, in order to compare the transverse permeability of this configuration with the analytical models found in the literature. Using the symmetry properties of the configuration, a minimal representative geometry, black square in Figure 15a, was used in the simulations in order to minimize the computing effort.

The simulations were carried out at different fiber volume fractions ranging from 0.2 to 0.7, and the computed non-dimensionalized permeability values, $K/r^2$, were plotted in Figure 16 as a function of fiber volume fraction. In the same graph, the analytical predictions from models such as Gebart, Gutowski, Kozeny-Carman, etc. [9] have also been plotted for comparison. It can be seen that a good agreement was found between the numerical and analytical values.
<figure><img src="image_22.png" /><figcaption>(a) Square arrangement of cylinders</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>(b) Minimal representative geometry</figcaption></figure>

<figure><img src="image_24.png" /><figcaption>**Figure 16.** Analytical models vs numerical simulation using OpenFoam for transversal flow in a square array of fibers.</figcaption></figure>
A similar numerical versus analytical investigation was also carried out on the flow around a hexagonal array of fibers. A minimal representative geometric domain (black square in Figure 17a) is also considered in this case, Figure 17b. The results of the simulation, Figure 18, are compared to those of the Gebart model and show very good agreement.

<figure><img src="image_25.png" /><figcaption>Figure 17. Simulation of flow around a hexagonal array of cylinders.</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>Figure 18. Analytical models vs numerical simulation using OpenFoam for transversal flow in a hexagonal array of fibers.</figcaption></figure>
The objective of this section is to validate the numerical approach to simple geome-
tries, solving the Navier–Stokes equations in the free space between the cylinders. In
the next section, this method, in addition to the dual-scale method based on the Darcy–
Brinkman equations, is described, and the flow simulation was performed in a realistic
3D textile geometry.

3.3. *Image-Based Permeability Determination*

The 3D representation of the textile meshes was acquired from a 4-layer stack E-glass plain weave specimen. During the stacking process, no relative interlayer change in orientation was allowed. Note that the same textile was experimentally characterized. The plain weave has a nominal areal weight of 800 g.m⁻² and a tow linear density of 2400 g.km⁻¹. The acquisition was realized using an X-ray microtomography (µCT) scanner, which provided volumetric data of the textile at grayscale, e.g., Figure 19, at three different compaction levels corresponding to fiber volume fractions of 0.5, 0.55, and 0.6. An image processing algorithm was then used to separate the inter-tow space from the tow data and to obtain the final 3D geometry. The physical dimensions of the scanned specimen were 25.14 mm × 26.04 mm with a variable height for each target fiber volume fraction (between 2.18 mm and 3.28 mm). The voxel dimensions were 2012 × 2084 pixels in the weft and warp directions and a variable number of pixels in the through-thickness direction (175 for the highest volume fraction 0.6 and 211 for the smallest one 0.5). The acquisition procedure is described in detail in the paper of Wijaya et al. [29].

For the given 3D geometry data, three different stereolithography (STL) files suitable for numerical simulation were generated using a MATLAB algorithm. These STL files corresponded to the volume of the warp tows, weft tows, and inter-tow space (Figure 20). The algorithm is based on a 3D image processing methodology developed by Wijaya [29], which allows the segmentation of the reconstructed grayscale data of a woven textile into three groups: weft and warp tows and inter-tow space. Segmentation is a crucial step in the processing of grayscale images and affects the accuracy of numerical simulations. The methodology used separates the different zones with satisfactory quality, but more sophisticated methods using AI, for example, are being developed. Guibert et al. [30] discussed the effects of image processing techniques on phase identification in grayscale scanned geometries in their paper.

<figure><img src="image_27.png" /><figcaption>**Figure 19.** 3D grayscale tomography reconstruction of the scanned sample.</figcaption></figure>

3.3.1. Single-Scale Flow

Here, the tows were considered solid impermeable zones (the fluid flows only in the inter-tow free space). In order to optimize the simulation time and allocated computational resources, only the inter-tow space was meshed (Figure 21), and the Navier–Stokes equations were solved in this domain. This approximation allows to reduce the computation resources used and is justified by the low permeability of the tows (10<sup>-13</sup>–10<sup>-14</sup> m<sup>2</sup>) [27].

The numerical approach described was used to compute the saturated permeability of the flow domain using the boundary conditions described in Figure 22. Note that the permeability values in the weft and warp directions were computed by specifying the axis of the pressure gradient, and, hence, the flow direction.
<figure><img src="image_28.png" /><figcaption>**Figure 20.** 3D geometry of the textile used for the numerical simulation.</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>**Figure 21.** Inter-tow space.</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>**Figure 22.** Boundary conditions for numerical permeability determination.</figcaption></figure>
The viscosity of the simulated test fluid was set to $\mu = 1$ mPa.s (close to the typical viscosity of industrially used resins [31]), and a pressure gradient $\Delta P = P_{in} - P_{out} = 990$ Pa is imposed between the inlet and the outlet.

The permeability values in the weft and warp directions were computed using Equation (6). The flow in the meshed domain was simulated using the steady-state solver *simpleFoam* in OpenFoam v7 [32]. The flow rate was computed by integrating the flow velocity at the outlet surface. The results of the numerical permeability in the weft and warp directions, compared to the corresponding experimental values for different fiber volume fractions, are shown in Figure 23a and Figure 23b, respectively. The numerical method underestimates the experimental values, as expected, since the flow inside the tows was neglected in the simulation. A large scatter between the numerical simulation and the experimental results for the warp direction is observed (up to 56%) compared to the weft direction (maximum of 25%), which can be attributed to the insufficient scan accuracy and the small size of the flow channels in the warp direction. Therefore, these narrow geometries are more difficult to mesh and can eventually be artificially removed if they are badly shaped, and the permeability in the warp direction is largely underestimated [31]. This problem can be solved if the flow inside the tows is also simulated, as explained in the following section.

<figure><img src="image_31.png" /><figcaption>**Figure 23.** Comparison between numerical and experimental values.</figcaption></figure>
### 3.3.2. Dual-Scale Flow

Here, the tows were considered permeable, and the simulated fluid was allowed to flow between and inside the tow spaces. The Brinkman equation was adapted to solve such a multiscale flow problem. The same boundary conditions used in the single-scale simulations, Figure 22, were again used. The complete domain (i.e., inter-tow space, weft, and warp tows), as depicted in Figure 24, was meshed, and the interface elements between different spaces were chosen depending on the surface orientation of the tows’ geometry using the specific command `useSurfaceOrientation` in OpenFOAM. This command allows one to select all cells inside a closed surface, the tows in this case, in order to assign them specific properties. The permeability of the tows was computed using the Gebart analytical model, assuming a hexagonal arrangement of fibers. The fiber radius of the used glass fabric, needed as an input for the Gebart model, was measured using microscopy images, and an average of 12.12 µm was found [14]. The local values of fiber volume fraction and permeability are given by the scan, whereas in the simulation, average and constant values are computed and assigned to the weft and warp tows, assuming that permeability tensors are diagonal for simplicity reasons since the off-diagonal permeability values do not affect the fluid flow in this case [33]. The values of permeability assigned to the weft and warp tows are given in equations 8 and 9. A geometric average was used to compute these values because of the non-periodic nature of the voxel’s permeability values. This approach was adapted from the averaging method used in random permeability distribution media in reservoir engineering [34].

$$ K_{\text{tows,warp}} = \begin{array}{c} 3.69 \cdot 10^{-12} \\ 0 \\ 0 \end{array} \begin{array}{c} 0 \\ 7.68 \cdot 10^{-13} \\ 0 \end{array} \begin{array}{c} 0 \\ 0 \\ 7.68 \cdot 10^{-13} \end{array} \text{m}^2 \quad (8) $$

$$ K_{\text{tows,weft}} = \begin{array}{c} 1.78 \cdot 10^{-12} \\ 0 \\ 0 \end{array} \begin{array}{c} 0 \\ 3.59 \cdot 10^{-13} \\ 0 \end{array} \begin{array}{c} 0 \\ 0 \\ 3.59 \cdot 10^{-13} \end{array} \text{m}^2 \quad (9) $$

<figure><img src="image_32.png" /><figcaption>Figure 24. Geometry used in the dual-scale case.</figcaption></figure>

An open-sourced multiscale Darcy–Brinkman solver called `hybridPorousInterFoam` was used for this simulation. This solver allows the simulation of two-phase flows in hybrid systems containing both porous matrices and solid-free regions [35]. A unique set of partial differential equations is used to represent the flow in both regions using the same mesh without the need to define complex interfacial conditions.

As can be observed in Figure 25, the permeability values obtained by this numerical method are closer to the experimental values than those obtained by the single-scale tests. The deviation from the experimental values is reduced to 11% for $V_f = 0.5$ and 8% for $V_f = 0.6$ in the weft direction, and it is around 20% in the warp direction; this is a significant improvement in permeability prediction accuracy compared to the single-scale simulation results. Nonetheless, there is still a deviation that can be reduced by refining
the mesh of the computational domain further, as demonstrated by Guibert et al. [36]. The experimental/numerical deviations in this case can be explained by the assumption of constant permeability and porosity in the weft and warp directions, neglecting local variations in these parameters.

<figure><img src="image_33.png" /><figcaption>**Figure 25.** Comparison between numerical simulation results and experimental values considering dual-scale flow.</figcaption></figure>

**4. Conclusions**

A 2D radial flow experimental setup was used to investigate the discrepancy between the unsaturated and saturated permeabilities and to observe the flow front during the fluid injection of different types of textiles. The experimental setup detailed in this article provides robust and repeatable unsaturated and saturated permeability data with reduced errors. The ratio between both permeabilities was between 0.9 and 1.42, which lay in the range found in the literature and was coherent with the flow front observations; the difference between both permeabilities was higher for the Unidirectional and Biaxial textiles corresponding to higher permeability ratios, for example. It was also confirmed that the flow front shape and the ratio $\frac{K_{unsat}}{K_{sat}}$ depend on the properties of the used textiles, such as
the stitching style, areal weight, and especially the size of the inter-tow space, which shows
the importance of the capillary effects on the fluid flow in porous media.

Two numerical methods for saturated permeability prediction were proposed. The first method is based on solving the Navier–Stokes equations in the inter-tow space and neglecting the permeability of the tows. In the second method, the flow inside the tows was considered using a specific OpenFOAM solver based on solving the Brinkman equations in multiscale domains. Both methods were used to measure the permeability of the scanned textile, and the results were compared with the experimental values given by the previously described setup. The dual-scale method gave results closer to those of the experiments, which shows the importance of using a multiscale approach to model the flow of fibrous reinforcements. The main challenge with these methods is their dependence on the quality of the scanned geometry, segmentation algorithm, and mesh refinement. Permeability determination using numerical simulations is a topic that is under development, as no standard method currently exists. The second part of this article contributes to the standardization efforts of virtual permeability determination, as international benchmark exercises are being launched to identify the different numerical methods that can be used for permeability measurement based on real images. The solver used in the dual-scale numerical method can also be used to simulate the resin flow front advancement during impregnation in order to visualize the filling scenario and to identify the locations of the potential dry spots, which can be the subject of upcoming work.

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

**Funding:** This research received no external funding.

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

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** No new data were created.

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

## References

1. Michaud, V. A Review of Non-saturated Resin Flow in Liquid Composite Moulding processes. *Transp. Porous. Med.* **2016**, *115*, 581–601. [CrossRef]

2. Pillai, K.M. Modeling the Unsaturated Flow in Liquid Composite Molding Processes: A Review and Some Thoughts. *J. Compos. Mater.* **2004**, *38*, 2097–2118. [CrossRef]

3. Vernet, N.; Ruiz, E.; Advani, S.; Alms, J.B.; Aubert, M.; Barburski, M.; Barari, B.; Beraud, J.M.; Berg, D.C.; Correia, N.; et al. Experimental determination of the permeability of engineering textiles: Benchmark II. *Compos. Part A Appl. Sci. Manuf.* **2014**, *61*, 172–184. [CrossRef]

4. Syerko, E.; Schmidt, T.; May, D.; Binetruy, C.; Advani, S.G.; Lomov, S.; Silva, L.; Abaimov, S.; Aissa, N.; Akhatov, I.; et al. Benchmark exercise on image-based permeability determination of engineering textiles: Microscale predictions. *Compos. Part A Appl. Sci. Manuf.* **2023**, *167*, 107397. [CrossRef]

5. Henry Darcy et les Ecoulements de Fluides en Milieu Poreux. Available online: https://www.researchgate.net/publication/266868182_Henry_Darcy_et_les_ecoulements_de_fluides_en_milieu_poreux (accessed on 20 August 2024).

6. Ketata, C.; Satish, M.G.; Islam, M.R. Revising Darcy's law: A necessary step toward progress in fluid mechanics and reservoir engineering. *WIT Trans. Eng. Sci.* **2006**, *52*, 1–8. Available online: https://www.witpress.com/elibrary/wit-transactions-on-engineering-sciences/52/16016 (accessed on 18 May 2024).

7. Chang, Z.; Sun, H.G.; Zhang, Y.; Zheng, C.; Min, F. Spatial Fractional Darcy's Law in Natural Reservoirs. *Phys. A Stat. Mech. Its Appl.* **2019**, *519*, 119–126. [CrossRef]

8. Hubbert, M.K. Darcy's law and the field equations of the flow of underground fluids. *Hydrol. Sci. J.* **1957**, *2*, 23–59. [CrossRef]

9. Karaki, M.; Younes, R.; Trochu, F.; Lafon, P. Progress in Experimental and Theoretical Evaluation Methods for Textile Permeability. *J. Compos. Sci.* **2019**, *3*, 73. [CrossRef]
10. Gebart, B.R. Permeability of Unidirectional Reinforcements for RTM. *J. Compos. Mater.* **1992**, *26*, 1100–1133. [CrossRef]

11. Sharma, S.; Siginer, D.A. Permeability Measurement Methods in Porous Media of Fiber Reinforced Composites. *Appl. Mech. Rev.* **2010**, *63*, 020802. [CrossRef]

12. Diallo, M.L.; Gauvin, R.; Trochu, F. Key factors affecting the permeability measurement in continuous fiber reinforcements. In *Proceedings of the ICCM 1997, Gold Coast, Australia, 14–18 July 1997*; pp. 441–445.

13. Khan, M.A.A. In-Plane Permeability Measurement of Biaxial Woven Fabrics by 2D-Radial Flow Method. *Sci. Eng. Compos. Mater.* **2021**, *28*, 153–159. [CrossRef]

14. Wijaya, W. Permeability of 2D Woven Composite Textile Reinforcements: Textile Geometry and Compaction, and Flow Modelling. Ph.D. Thesis, The University of Auckland, Auckland, New Zealand, 2020.

15. May, D.; Aktas, A.; Advani, S.G.; Berg, D.C.; Endruweit, A.; Fauster, E.; Lomov, S.V.; Long, A.; Mitschang, P.; Abaimov, S.; et al. In-plane permeability characterization of engineering textiles based on radial flow experiments: A benchmark exercise. *Compos. Part A Appl. Sci. Manuf.* **2019**, *121*, 100–114. [CrossRef]

16. Arbter, R.; Béraud, J.M.; Binetruy, C.; Bizet, L.; Bréard, J.; Comas-Cardona, S.; Demaria, C.; Endruweit, A.; Ermanni, P.; Gommer, F.; et al. Experimental determination of the permeability of textiles: A benchmark exercise. *Compos. Part A Appl. Sci. Manuf.* **2011**, *42*, 1157–1168. [CrossRef]

17. Li, C.; Cantarel, A.; Gong, X. A study on resin infusion and effects of reinforcement structure at dual scales by a quasi-realistic numerical simulation method. *J. Compos. Mater.* **2020**, *54*, 4157–4171. [CrossRef]

18. Tournier, L. Method for Optical Measure of Local Reinforcement Permeability and Infusion Behavior Prediction in Liquid Molding Processes. Master's Thesis, Centre for Advanced Composite Materials, The University of Auckland, Auckland, New Zealand, 2010.

19. Zhang, F.; Comas-Cardona, S.; Binetruy, C. Statistical modeling of in-plane permeability of non-woven random fibrous reinforcement. *Compos. Sci. Technol.* **2012**, *72*, 1368–1379. [CrossRef]

20. Comas-Cardona, S.; Cosson, B.; Bickerton, S.; Binetruy, C. An optically-based inverse method to measure in-plane permeability fields of fibrous reinforcements. *Compos. Part A Appl. Sci. Manuf.* **2014**, *7*, 1–48. [CrossRef]

21. ISO 4410:2023; Test Methods for the Experimental Characterization of In-Plane Permeability of Fibrous Reinforcements for Liquid Composite Moulding. International Organization for Standardization: Geneva, Switzerland, 2023.

22. Weitzenbock, J.R.; Shenoi, R.A.; Wilson, P.A. Radial flow permeability measurement. Part A: Theory. *Compos. Part A Appl. Sci. Manuf.* **1999**, *30*, 781–796. [CrossRef]

23. Pomeroy, R.; Grove, S.; Summerscales, J.; Wang, Y.; Harper, A. Measurement of permeability of continuous filament mat glass-fiber reinforcements by saturated radial airflow. *Compos. Part A Appl. Sci. Manuf.* **2007**, *38*, 1439–1443. [CrossRef]

24. Foley, M.F.; Gutowski, T. The effect of process variables on permeability in the flexible resin transfer molding (FRTM) process. In *Proceedings of the 36th International SAMPE Symposium and Exhibition, San Diego, CA, USA, 15–18 April 1991*; pp. 326–340.

25. Gebart, B.R.; Lidström, P. Measurement of in-plane permeability of anisotropic fiber reinforcements. *Polym. Compos.* **1996**, *17*, 43–51. [CrossRef]

26. Pollard, M. Permeabilities of Fiber Mats Used in Resin Transfer Molding. In *Proceedings of the 24th International SAMPE Technical Conference, Toronto, ON, Canada, 20–22 October 1992*; pp. T408–T420.

27. Geoffre, A.; Wielhorski, Y.; Moulin, N.; Bruchon, J.; Drapier, S.; Liotier, P.J. Influence of intra-yarn flows on whole 3D woven fabric numerical permeability: From Stokes to Stokes-Darcy simulations. *Int. J. Multiph. Flow* **2020**, *129*, 103349. [CrossRef]

28. Schmidt, T.; May, D.; Duhovic, M.; Widera, A.; Hümbert, M.; Mitschang, P. A combined experimental-numerical approach for permeability characterization of engineering textiles. *Polym. Compos.* **2021**, *42*, 3363–3379. [CrossRef]

29. Wijaya, W.; Kelly, P.A.; Bickerton, S. A novel methodology to construct periodic multi-layer 2D woven unit cells with random nesting configurations directly from µCT-scans. *Compos. Sci. Technol.* **2020**, *193*, 108125. [CrossRef]

30. Guibert, R.; Nazarova, M.; Voltolini, M.; Beretta, T.; Debenest, G.; Creux, P. Influence of Standard Image Processing of 3D X-ray Microscopy on Morphology, Topology and Effective Properties. *Energies* **2022**, *15*, 7796. [CrossRef]

31. Han, W. Composites Manufacturing by Injection of Reactive Thermoplastic Resin in a Fibrous Preform. Ph.D. Thesis, Université Paul Sabatier—Toulouse III, Toulouse, France, 2022.

32. Peralta, C.; Nugusse, H.; Kokilavani, S.P.; Schmidt, J.; Stoevesandt, B. Validation of the simplefoam (rans) solver for the atmospheric boundary layer in complex terrain. In *ITM Web of Conferences*; EDP Sciences: London, UK, 2014; Volume 2.

33. Guibert, R.; Horgue, P.; Debenest, G.; Quintard, M. A comparison of various methods for the numerical evaluation of porous media permeability tensor from pore-scale geometry. *Math. Geosci.* **2016**, *48*, 329–347. [CrossRef]

34. Ahmed, T. *Reservoir Engineering Handbook*; Gulf Professional Publishing: Houston, TX, USA, 2018.

35. Carrillo, F.J.; Bourg, I.C.; Soulaine, C. Multiphase flow modeling in multiscale porous media: An open-source micro-continuum approach. *J. Comput. Phys.* **X** **2020**, *8*, 100073. [CrossRef]

36. Guibert, R.; Nazarova, M.; Horgue, P.; Hamon, G.; Creux, P.; Debenest, G. Computational permeability determination from pore-scale imaging: Sample size, mesh and method sensitivities. *Transp. Porous Media* **2015**, *107*, 641–656. [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.