1. Introduction
Liquid transfer moulding techniques, including resin transfer moulding (RTM), are increasingly used for the production of fibre-reinforced polymer composites. Their advantages over hand lamination (high productivity at a lower cost, increased fibre-to-resin ratio, controlled dimensional tolerances providing outstanding reproducibility) allow their application in the production of structural components for the automotive and aircraft industries [
1,
2]. These transport industries have strict safety requirements, which can only be met if the flame retardant (FR) properties of the polymer composite matrix materials are improved. For this purpose, a wide range of FR additives are available, offering the possibility to consider the effect of the FR on the mechanical and thermal properties, environmental issues creating risks for human life and the environment, waste treatment, and recycling [
3,
4,
5]. When the RTM method is applied for the processing of an FR composite, besides these properties, the phase of the additive and its effect on the viscosity of the resin also have to be taken into account. During the RTM process, a liquid resin is injected into a preform consisting of fibre reinforcement, and therefore, solid phase additives may be filtered by the reinforcement layers, leading to the non-uniform distribution of additives, and uneven fire performance as a result. For this reason, liquid-phase flame retardants can be more easily integrated into the RTM process, as only their effect on resin viscosity has to be taken into account. Hand lamination is less sensitive to increased viscosities, and as the reinforcement layers are impregnated one by one, no filtration of the solid additives occurs. On the other hand, hand lamination is less productive and the achievable maximum fibre content is lower. Therefore, RTM is preferred in many industries. In many applications, different additives have to be incorporated into the polymer matrix to enhance existing or provide new properties [
6,
7,
8,
9,
10]. The most common additive types applied in thermoset polymer composites are:
formulation additives:
air release agents eliminating or preventing the formation of air bubbles (e.g., silicon-containing or silicon-free polymers)
dispersants reducing the viscosity and improving the flow behaviour (e.g., various low to high molecular mass polymers)
wetting agents and surface modifiers preventing surface defects (e.g., polysiloxanes, polyacrylates)
performance additives:
light stabilizers (e.g., hindered amines)
thermal stabilizers (e.g., phenol and phosphite type antioxidants)
flame retardants (e.g., metal hydroxides, various halogen- and phosphorus-containing additives)
antistatic additives (e.g., quaterner ammonium salts of fatty acids, alkyl sulfonates, alkyl phosphates, alkyl dithiocarbamates, alkyl carboxylates)
conductive additives (e.g., quaterner ammonium salts of fatty acids, alkyl sulfonates, alkyl phosphates, alkyl dithiocarbamates, alkyl carboxylates)
colourants:
The proper dispersion of the additives is a particular challenge when solid-phase additives are added in relatively large amounts to the polymer matrix, as in the case of many common FRs as aluminium trihydroxide, magnesium hydroxide, or ammonium polyphosphate. Accordingly, it is important to explore the possibilities to improve their distribution in composites produced by RTM. Therefore, the present article reviews:
the parameters influencing particle distribution during RTM,
the models established for the description of particle distribution in composites prepared by RTM,
the analytical methods for determining particle distribution in composites manufactured by RTM,
the possible solutions to improve the particle distribution of solid-phase additives.
As the literature on the distribution of FRs is quite limited, publications related to the distribution of other solid additives are also discussed; with the primary aim of facilitating the production of flame retarded structural composites by RTM in the case of solid FRs.
2. Parameters Influencing Particle Distribution during RTM
Liquid moulding techniques, such as resin transfer moulding (RTM), are commonly used to manufacture fibre-reinforced composites. In the RTM process, a fibrous preform (lay-up sequence) is inserted into a rigid half-mould, and the mould is then closed. Following this, resin is injected into the mould with a positive gradient pressure through the gate points, replacing the air entrapped within the preform. At this point, the curing phase begins. After curing, the mould is opened and the composite part is removed. In some cases, a post-curing phase is necessary in order to guarantee a high degree of cross-linking. If solid phase additives are applied in resin formulations, a homogenous and stable suspension of the particles is of paramount importance in these processes [
11,
12]. To avoid the filtration of the solid particles by the reinforcement, it is important to understand the mechanism of filtering and the process parameters determining its extent. Fundamentally, there are three possibilities in particle-filled resin systems [
13]:
No retention: Particles can flow freely through fibre reinforcement. Uniform particle distribution is possible.
Deep filtration: Distribution is not uniform. The particles are gathered at the inlet, and the concentration decreases with a higher filtration length.
Cake filtration: Particles cannot enter the fibre reinforcement. A cake of particles is formed outside the preform entry.
The most relevant process parameters influencing distribution are [
1]:
the pressure applied during RTM,
the temperature of the resin/suspension,
the viscosity of the resin/suspension,
the particle size of the solid additive,
the pore size of the fibrous media,
the permeability of the fibre reinforcement,
the fibre content of the composite,
the orientation of the fibres,
the location of the inlet,
the filling time.
However, the process parameters mentioned above are closely related to each other, e.g., when the temperature increases, the viscosity of the resin decreases, the resin can flow in the cavity more easily, the filling time decreases, and the injection pressure is lower [
1]. In the reverse case, when the temperature is lower, the viscosity increases, and so do the filling time and the pressure [
1]. There are also connections between the other parameters, which should be taken into consideration.
The ratio between the mean particle diameter and the pore size of the fibrous media is another important aspect to be considered [
14,
15]. If the resin is filled with large particles (d ≥ 30 µm), volume phenomena prevail over surface phenomena, while for small particles (d~1 µm), surface interactions are dominant. Within these two particle diameters, both volume and surface phenomena can occur [
1,
14]. Generally, cake filtration occurs when the particle size is larger than the pore size, in which case particles cannot enter the fibre layer and accumulate, forming a cake. When the resin is filled with small particles, the suspension can flow through the fibrous media and only little filtration takes place. For an intermediate ratio of the particle diameter to pore size, deep filtration is the main mechanism [
16]. When deep filtration occurs, larger particles are gradually deposited. This can narrow and clog the flow channels within the fibre bundles, which ultimately leads to cake filtration [
16].
3. Modelling Particle Distribution in RTM
The filtration of solid particles is often modelled with analytical and/or numerical simulations in the literature. The main model for one-dimensional flow and filtration is Darcy’s law (1):
where
Q (m
3/s) is the volumetric flow rate,
A (m
2) is the cross-sectional area of the porous channel,
η (Pas) is the viscosity of the fluid, and Δ
p (Pa) is the pressure drop measured over a flow distance
L (m).
K (m
2) represents the permeability of the porous media. Lefevre et al. efficiently coupled Darcy’s law (flow model) to a filtration model to describe particle distribution in the cross direction of a composite manufactured by RTM [
17,
18]. Reia da Costa et al. presented analytical and numerical models of liquid moulding of nanoparticle filled thermosets [
19]. These models were also based on Darcy’s law. The authors found a non-linear model incorporating variations in permeability, viscosity, and porosity as a function of local nanoparticle loading. It is obvious that the filtration of the particles by the reinforcement leads to variations in material properties. According to Reia da Costa et al. [
19], the narrowing of the reinforcement flow channels caused by the accumulation of nanoparticles results in a reduction of porosity as the resin flow front progresses. They described the relationship between porosity and permeability with the Kozeny-Carman relation, although it should be noted that under the conditions used in their work, the influence on porosity and permeability is negligible due to the small particle retention concentration compared to the volume fraction of the liquid resin. The authors found that the main parameter in this non-linear model was the viscosity of the suspension. Variation in the suspended particle concentration lead to variations in the viscosity. They found a non-linear dependence between the particle concentration and the viscosity. They compared this model to the linear analytical solution in order to validate it and found that these two models are suitable for the simulation of flow and filtration in the liquid moulding of nanoparticle-filled resins. Analytical approximation can be used for short filling lengths, while the non-linear model can be used for infiltration lengths in the meter range. The models developed can serve as a basis for process design or optimization for the production of fibre-reinforced nanoparticle-filled composites where an even distribution of particles is required.
Hwang et al. developed a direct numerical simulation technique to describe particle filtration and flow [
20]. They determined hydrodynamic interactions between filler particles and fluid using Stokes-Brinkman coupling to describe the flow in a dual-scale porous media (macro-pores in the intertow space; micro-pores in the intratow space). They modelled fluid flow with Stokes flow, and used the Brinkmann equation as a momentum equation, considering the porous media as a continuum body that can be characterized with its permeability. They found that particle deposition was enhanced in the porous media with a higher permeability. Two mechanisms for this enhancement were reported: one is the penetrating flow into the porous tow and the other is the formation of a downward flow because of the presence of the particles near the porous media. Downward flow occurs in the gap between the particle and porous wall in the direction of particle movement. According to Hwang et al., this flow is responsible for the early deposition of the particle, because this kind of flow is not present in the absence of the particle and yields additional downward penetrating flow through the porous media. In another paper, the fluid-flow-driven motion of particles was studied with Voronoi discretization and minimization of the dissipation rate of energy; they described two-dimensional particle flow and filtration through fibre bundles [
21,
22]. The results of the simulation were compared to experimental observations. The simulation suggested that clogging occurs if a space is narrower than twice the diameter of the largest particles. Three particle size regions were considered by Frischfelds et al. [
21,
22]: a medium one with the largest particle size, with approximately half the average spacing between the fibres in the bundles; a second one with twice as large particles as the first set; and the third with half as large particles as the medium set of particles. When they doubled the size of the particles during flow simulation, the overall clogging increased. The particles formed bridges in the channels between fibre bundles. This built up a pressure gradient that pumped particles into the larger openings inside the bundle. Particles enter spaces of the fibre bundle until they get trapped and the further motion of particles at that location is not possible. Experimental results showed a higher accumulation of particles at the front of the fibre bundle, as predicted by the model. Elgafy et al. used a simulation model based on the Eulerian multiphase flow approach to investigate the flow characteristics of carbon nanoparticle filled fluids around a carbon microfibre matrix, and also investigated the interactions between microfibre walls and nanoparticles [
23]. They found that the interactions between the microfibre sidewalls and the interfacial fluid layers have a tendency to reduce the flow velocity, causing an attraction of the flow to the microfibre sidewalls. After some time, the flow passages may be blocked. Chohra et al. suggested a so called “sieve mechanism” to describe the macroscopic filtration of the particles: large particles accumulate in front of small pores, causing a concentration distribution within a dual-scale fibrous media [
24]; the ratio of particle size to pore size obviously matters. The authors also performed an experimental study on the behaviour of microparticles in dual-scale fibrous media to validate the proposed model. They injected particle-filled resin suspensions along the perpendicular direction of the fibre mats and measured the amount of filtered particles within each layer after rinsing the layers and collecting the particles. Nordlund et al. applied microscopic imaging and microparticle image velocimetry measurements to investigate filtration, using a colourless glycerol/water mixture and fluorescent particles [
25]. In the two studies mentioned above, liquid suspensions were investigated instead of fully cured composites. Erdal et al. presented another macroscopic model to investigate particle filtration during RTM, using Darcy’s law and particle retention kinetics [
26]. This model was improved by Lefevre et al. to take liquid retention into consideration [
17]. The authors found another type of behaviour: the filler concentration along the composite part was merely varying compared to the model presented by Erdal et al. [
26]. They found a decrease in filler concentration near the inlet, which was attributed to retention, whereas the accretion near the outlet was assumed to depend on a liquid depletion mechanism in the suspension in the flow front area. These phenomena were taken into account during the improvement of the filtration model. In another paper, Lefevre et al. proposed an enhanced model, which takes into account the changes in viscosity and permeability [
18]. According to them, the flow results from a competition between a low permeability and viscosity at the mould inlet and a high permeability and viscosity at the outlet. As the fibrous preform retains fillers, the permeability of the new filter system consisting of the fibrous preform and the retained particles had to be evaluated. Their model also provided clogging detection by monitoring the evolution of the porosity with respect to time: when the porosity was zero, there was no more available volume for the suspension flow.
4. Analytical Methods for Determining Particle Distribution in Composites Manufactured by RTM
Several techniques exist for the investigation of particle distribution. However, most of them are based on optical methods such as scanning electron microscopy (SEM) or transmission electron microscopy (TEM). In many research papers, SEM is combined with other techniques, e.g., wavelength dispersive spectrometry (WDS) and electron probe microanalysis (EPMA). A combination of SEM/EPMA/WDS can also be used in fully cured composite parts.
Yum et al. developed a new methodology for measuring the distribution of particles in fully cured composites via elemental analysis [
16]. The matrix resin used was an epoxy system consisting of a bisphenol-A resin with an amine hardener. The fillers were spherical titanium dioxide or carbon nanotube (CNT)-silver particles, and dual-scale fibre mats consisting of E-glass and carbon fibre were used. The width of the flow channels in the intertow region may be up to several hundreds of microns, while in the intratow region, it is a few microns. For qualitative analysis, SEM images were combined with EPMA. This combination, called mapping analysis, provides a visual representation of fibre geometry and the distribution of the particles, as shown in
Figure 1. Quantitative analysis is also possible; the filler concentration (particle size range from nanometers to micrometers) on the fully cured composite parts can be measured both in the intertow and intratow region with the use of WDS provided by EPMA. Yum et al. also stated that the degree of dispersion of the filler particles in the resin influences the mechanism of filtration. The resin suspension can be classified into two categories according to the degree of dispersion: well-dispersed particles and clusters of particles [
16,
27]. It is obvious that the ratio of particle size to pore size determines the filtration mechanism. Clusters are mostly captured by the fibrous preform, and because of their size, cake filtration occurs. Well-dispersed particles—as expected—can flow through the pores of the fibrous media; the main filtration mechanism is deep filtration. With the combination of SEM, EPMA, and WDS, the spatial distribution of the filler particles can be investigated in a fully cured composite part both inside and outside the fibrous media [
16].
Fernberg et al. also investigated the mechanisms of particle distribution in infusion moulded composites [
28]. They used aluminatrihydroxide (ATH) as a flame retardant; the resin suspension was processed by infusion moulding techniques, including resin transfer moulding. They took into consideration the viscosity of the suspension. They found several models that rely on empirical calibrations [
29,
30,
31], but the following formula had the best correlation with the experimental data (Equation (2)) [
30]:
where
η (Pas) is the viscosity of the suspension,
ηs (Pas) is the viscosity of the resin,
Φ (-) is the phase volume of the filler,
Φm (-) is the maximum packing fraction of particles, and [
η] (Pas) is the intrinsic viscosity. The latter is rather a limit of viscosity than an accurate value; it has a physical meaning: it is the viscosity of extremely dilute solutions of a resin in a solvent and is related to the molecular weight of the resin [
32]. Equation (2) shows that if
Φ =
Φm, then the viscosity of suspension reaches a very high level, mathematically
η → ∞. If there is no filler in the resin (
Φ = 0), we cannot talk about resin suspension, only about neat resin; Equation (2) interprets this too: in this case,
η =
ηs. Between these two limits, Equation (2) describes the evolution of suspension viscosity depending on the filler content. Intrinsic viscosity can be determined by the extrapolation of experimental values. Fernberg et al. used an orthophtalic resin designed for vacuum infusion with methyl-ethyl-ketone peroxide hardener, ATH-particles with different sizes, and plain glass fibre weave with an average weight of 816 g/m
2. Rheological measurements of the resin suspensions with different filler types at various filler mass fractions were carried out. Although there were some discrepancies from the theory expected, Fernberg et al. found that viscosity increases with the decrease of particle size at the same filler volume fraction. At lower concentrations, the suspensions act like Newtonian fluids, but at high concentrations of smaller particles (0.5–5 µm), the suspensions showed a Bingham-type behaviour [
18,
33,
34,
35]. This means that no viscous flow occurs under a certain shear stress limit. This Bingham-type behaviour can have a big influence on the mould filling process in liquid infusion moulding processes. The microscopy analysis of the final microstructure of the infusion moulded composites showed that particle-rich areas are always in contact with a fibre bundle; ATH particles can also be observed within the bundles, preferably in resin-rich areas. The authors found that the permeability of the fibre mats does not decrease with flow length, contrary to what they expected. Fernberg et al. expected a decrease in permeability with an increasing flow front position in the case of particle filled resin suspensions because of particle filtration. Interestingly, the permeability did not decrease with flow length; however, after an initial lower permeability, it increased with flow length, and then reached a constant value. According to Fernberg et al. [
28], this indicates that particle filtering over time or length only has a minor influence on the overall macroscopic flow behaviour. However, they did not find any sign of systematic differences in the particle distribution along the flow direction. Nevertheless, if the strain rate in a section of the flow channel becomes too low (for example: inadequate flow channel geometry), no flow will occur in that section. In this case, the changes in flow channel geometry can be considered part of the porous structure, decreasing the permeability and providing resistance to flow, thus allowing filtration.
Filtration can be investigated in fully cured composites with a so called “burn-off” technique. In this test, a certain composite part is cut into samples layer by layer and both the matrix and the fibre material of each sample are burned off, so that the mass of the filtered microparticles can be measured in each layer. For this method, both the matrix and the fibrous media should be completely removed by calcination, while filler particles should be non-flammable. Lefevre et al. validated their models with this experimental method [
17,
18]. In their experiments, the matrix was polyester and the fibre reinforcement was from synthetic polyethylene terephthalate (PET), because PET fibres can be completely removed by burning; the non-flammable filler particles were glass beads. Louis et al. manufactured resin transfer moulded epoxy composites filled with aluminium oxide (Al
2O
3) nanoparticles and used aramid fabric as reinforcement [
13]. After curing, they examined the local particle distribution by the above-mentioned burn-off method. They measured the resin flow time and particle concentration as process parameters to determine their impact on the final particle distribution in fully-cured composite parts. The burn-off technique is an adaptation of the conventional matrix burn-off (e.g., ASTM D2584 [
36]). Because of the aramid reinforcement, they raised the temperature to 1000 °C during the burn-off test, allowing full thermal decomposition. As alumina was used as the filler, only its non-combustible particles remained after burning. They found that 3 h are sufficient at 1000 °C for the burn-off of nominal sample sizes of 0.3 cm × 1.0 cm × 2.0 cm. They also stated that this temperature and method are suitable for composites made with carbon fibres. They observed a non-uniform particle distribution: the highest particle concentration occurred near the inlet and it decreased with an increasing distance from the inlet. They found that with an increasing resin infusion time, the local particle concentration in the laminate increased. Although particle distribution was non-uniform, a desired minimum amount of particles was found across the entire length of the laminate. Louis et al. [
13] also investigated the effect of particle concentration in the resin. Their aim was to see whether a lower particle concentration is less sensitive to filtration. Process conditions were controlled, e.g., the same volume of resin and the same overflow conditions (Louis et al. [
13] also investigated the effect of fixed resin overflows (5%, 30%, 50%) on the mass percentage of particles retained). However, with a reduced particle concentration, particle filtration still occurred. They concluded that a lower particle concentration does not reduce particle clogging and filtration, and that these are independent of the particle concentration in the resin suspension.
Garay et al. investigated the influence of calcium carbonate on RTM processing and the properties of the moulded composites [
37]. Their aim was to investigate resin characteristics, reinforcement permeability, and the mechanical properties of the final composite, while varying the amount of calcium carbonate. Orthophtalic polyester resin was used as the matrix, a combination mat (E-glass fibre mats with polypropylene (PP) flow media core) as reinforcement, and calcium carbonate as the filler, with a mean diameter of 18.22 µm. Particle distribution in the cured composite parts was investigated with the “burn-off” technique according to ASTM D2734 and ASTM D5630. They found that with an increasing calcium carbonate content in the resin, the permeability decreased and a longer time was necessary to fill the mould. The presence of the filler particles between the fibres reduces the porosity of the reinforcement and hinders the flow of resin [
25]. Flow in a fibrous media can be classified into two types:
The reinforcement used in the experiment was a glass mat with non-continuous, randomly oriented fibre bundles. The presence of the filler was expected to slow down both macro- and micro-flows, although particles were also expected to be found between and within fibre bundles [
28]. Garay et al. did not find particle filtration, regardless of the amount of calcium carbonate in the resin suspension, which agrees with the results that Fernberg et al. presented [
28]. However, the mass fraction of the glass fibre was not uniform along the mould because of the randomly oriented non-continuous glass fibres and also due to the fact that Garay et al. only used one layer. The results of their study showed that the presence of calcium carbonate impacted the RTM process and the properties of the composite. With the addition of the inorganic filler, the gel temperature decreased and the gel time increased. By increasing the CaCO
3 content in the resin, the maximum curing temperature decreased and the reactivity of the system decreased. CaCO
3 has a higher thermal conductivity coefficient (0.24–0.30 W/Km) in comparison to polyester (0.15–0.24 W/Km) and it acts as a kind of a physical reaction inhibitor. The presence of the filler increased the viscosity of the resin suspension according to Equation (2), exceeding the desired value for RTM at higher filler loadings. Furthermore, the addition of CaCO
3 considerably decreased the permeability of the reinforcement and increased the mould filling time. Although no particle filtration was observed, due to the fibre mat used, these variations in the properties might lead to difficulties in proper RTM processing and manufacturing.
In addition to the different measurement methods mentioned above (SEM, TEM, WDS, EPMA etc.), there are some alternative methods to investigate particle distribution in fibre-reinforced particle-filled composites. Moreover, these alternative techniques have some advantages over methods involving electron microscopy, such as less complicated sample preparation. A very promising non-destructive technique for microstructure characterization is micro X-ray computed tomography (X-ray CT), although it is only suitable for small specimens [
39,
40]. Confocal fluorescence microscopy (CFM) provides microstructural characterization similar to electron microscopic methods (TEM, SEM), but as its resolution is lower, larger areas can be investigated. For relevant characterization, the material imaged must be autofluorescent or contain fluorescent particles of known properties (size, interactions with the matrix, and the fibres) [
41]. Unfortunately, there is a risk that the filtration of these fluorescent particles occurs during RTM processing, which makes confocal fluorescence microscopy unsuitable for structure characterization. Another possible alternative is confocal Raman imaging (CRI). The distribution of the components can be identified with this method at a microscopic scale based on their chemical composition [
42,
43]. CRI is a combination of confocal microscopy and Raman spectrometry, which analyses the inelastic scattering of coherent monochromatic light produced by a laser. In CRI, Raman spectra are collected with high throughput. The entire Raman spectrum is collected for each pixel, providing a chemical map of the sample [
44,
45]. Schmidt et al. reported that confocal Raman imaging can be used effectively for the spatial imaging of local phases in a polymeric matrix [
46]. Gallos et al. applied CRI to characterize the microstructure of polycaprolactone (PCL) composites reinforced with natural lignocellulosic hemp fibres [
47]. Although the samples investigated were prepared with conventional injection moulding for thermoplastics, the paper provides a possible use of the method for the microstructural investigation of polymer composites. They found that CRI is capable of providing information for the entire multicomponent structure. Their results showed that the high sensitivity of the method allows a qualitative evaluation of particle distribution in the composite materials. They claim that CRI can be used for tracking non-fluorescent additives spread in a polymeric matrix over a large surface area. They also investigated the noise that autofluorescence can cause in the Raman signal. For this purpose, they compared CRI with a reference technique based on confocal laser scanning microscopy (CLSM). Their results demonstrated two advantages of CRI against CLSM:
the ability to characterize small particles that are not visible with CLSM,
less sample preparation before the measurements.
Figure 2 shows the comparison between CRI and CLSM analysis of the same sample. It is conspicuous that confocal Raman imaging highlights the presence of small objects, particles, while CLSM does not show these.
Shojaee et al. studied the dispersion and distribution of graphene in an epoxy matrix [
48]. To produce phase images, they used confocal Raman imaging. Based on the images taken by CRI, they found not only different distributions, but also different agglomerate sizes in the case of different graphene concentrations, as shown in
Figure 3. They stated that the average size of graphene agglomerates increases with an increasing graphene concentration (
Figure 4.). Additionally, Shojaee et al. found a jump in the average agglomerate size above a graphene loading of 0.1 wt%. According to them, this jump indicates the optimum graphene concentration for homogenous dispersions. Their findings were in good agreement with the results of another study by Rafiee et al. [
49]. Their research focused on graphene/epoxy systems, but obviously confocal Raman imaging can be used to study dispersion and distribution in other additive/polymer systems with appropriate Raman signals. They claimed that CRI provides an accurate, reproducible, and non-destructive way to quantify graphene dispersion in polymer composites. With this method, homogenous distribution can be achieved in composites. Their study suggests that the use of confocal Raman imaging should be taken into consideration as a possible and accurate test method for the investigation of particle distribution in thermosets and their composites.
The main areas of use, along with the advantages and disadvantages of the above mentioned analytical methods for determining particle distribution in composites manufactured by RTM, are summarized in
Table 1.