Theoretical Assessment of the Environmental Impact of the Preheating Stage in Thermoplastic Composite Processing: A Step toward Sustainable Manufacturing

: Manufacturing processes have always played a pivotal role in the life cycle assessment of products, necessitating focused efforts to minimize their impact on the environment. Thermoplastic composite manufacturing is no exception to this concern. Within thermoplastic composite manufacturing, the preheating process stands out as one of the most energy-intensive stages, significantly affecting the environment. In this study, a theoretical analysis is conducted to compare three modes of preheating: conductive, radiative, and convective modes, considering their energy consumption and environmental impact. The analysis reveals the potential for substantial energy savings and emissions reduction through the selection of a proper preheating mode. Since the analysis used in this study is theoretical, it facilitates a parametric study of different modes of preheating to assess how process parameters impact the environment. Moreover, this study includes a comparison between emissions from material production and the preheating process, highlighting the substantial contribution of the preheating process to the overall product life cycle assessment.


Introduction
Even though many countries around the globe have been striving to pass legislation to move toward net-zero emissions to control global warming, there is no indication that this ambitious objective has been achieved [1].To elaborate, emissions from the industry and energy production sectors have increased by 60% since the establishment of the United Nations Framework Convention on Climate Change in 1992 [2].This trend has concerned researchers and industrial leaders to take action to reduce emissions.Manufacturing-related industries are among the largest contributors to harmful gas emissions, accounting for 30% of greenhouse gas emissions [3].As such, developing sustainable manufacturing processes with reduced environmental impacts is extremely beneficial to address the concerning trend in industrial emissions [4].Per the US Department of Commerce [5], sustainable manufacturing is "the creation of manufactured products that use processes that are nonpolluting, conserve energy and natural resources, and are economically sound and safe for employees, communities, and consumers".One of the most effective ways in this direction is to integrate environmental impact indicators into the design process of products and their associated manufacturing processes.This strategy is referred to as Design For Environment (DfE), which can be implemented through life cycle assessment methods [6].
Life Cycle Assessment (LCA) is a technique to quantitatively evaluate the consumption of all resources and their corresponding emissions in a defined system [6].It consists of four stages, namely goal and scope, inventory analysis, impact analysis, and interpretation of the results [7][8][9].Two main techniques are employed for conducting LCAs: processlevel analysis and economic input-output analysis.Although the economic input-output analysis method provides a great insight to the environmental impacts of products and processes, most LCAs conducted in the literature are based on the process-level analysis techniques, in which the resource uses, environmental releases from the production process, and some important contributions from suppliers are assessed in detail [7].
The composite manufacturing industry has been using life cycle analysis techniques for many years [10][11][12][13][14]. Al-Lami et al. [15] developed an eco-efficiency assessment model for the manufacturing of Carbon Fiber-Reinforced Polymer (CFRP) composite parts based on life cycle assessment and life cycle cost analysis.In their work, they introduced an activitybased, bottom-up decision-making tool, through which the energy consumption and the cost associated with the manufacturing of carbon fiber composite parts were optimized.Khalil [16] performed a comparative impact assessment for the production of conventional aluminum alloy (AlMg3), used in aircraft fuselages, and its lighter counterpart, i.e., carbon fiber-reinforced polymer composites.Moreover, they demonstrated that the type of polymer used in the CFRP composite materials affects mid-point impact categories.To elucidate, it is shown that replacing the epoxy resin used in the CFRP polymeric phase with polypropylene polymer could result in a notable reduction in the energy intensity and mid-point impact categories associated with CFRP production.They concluded that the production of CFRP material with a polypropylene matrix leads to lower impact categories compared to those associated with AlMg3 production.Song et al. [7] analyzed the environmental impact of fiber-reinforced polymer composites made by the pultrusion process using a hybrid life cycle assessment technique, in which both process-level and economic input-output analysis techniques were employed.The intention behind this hybrid method was to make the life cycle assessment more efficient and accurate by minimizing the disadvantages of each method.
Despite the existence of many research works in the field of composite material LCA, research focused on the parametric study of the environmental impact assessment of composite manufacturing processes is modest.Hybrid analytical-numerical methods have proven to be useful tools that allow researchers to study manufacturing processes efficiently [17].The analysis of energy consumption, as an output of the manufacturing processes, is not an exception to this matter.In other words, hybrid analytical-numerical methods can be used to parametrically calculate the energy consumption of each manufacturing process.Doing so will pave the way to optimize manufacturing processes by controlling the process parameters to minimize their environmental impact.It is worth mentioning that most of the LCA research works reported in the literature heavily rely on published data or LCA databases, where the manufacturing processes are modeled as an aggregated unit [18].This limits researchers' ability to analyze how manufacturing process parameters can impact the LCA of a product.
In view of the extent of the literature, this research is presented to showcase the implementation of the hybrid analytical-numerical methods to parametrically study the effect of manufacturing processes and their process parameters on environmental impacts.For this study, the preheating stage of the thermoforming process for fiber-reinforced thermoplastic polymer composites is analyzed because of the high energy consumption associated with this stage of thermoplastic composite processing.The heating stage of thermoforming processes has been extensively investigated in the literature [19][20][21].Längauer et al. [19] developed an enhanced model for the infrared heating of thermoplastic composites, highlighting that the accuracy of simulation results can be improved by incorporating anisotropic thermal conductivity and temperature and pressure dependencies of specific heat capacity.Other related studies have also addressed the modeling of anisotropic thermal conductivity in thermoplastic composites and its impact on material processing.For instance, a novel model was developed in [22] to analytically determine the anisotropic thermal conductivity of thermoplastic composites.
The current research is structured as follows: first, the goal and scope of this study are discussed.Subsequently, inventory analysis will be conducted, in which the hybrid analytical-numerical method is put into practice to determine the energy consumption of the preheating stage of the thermoforming process.Three different modes of preheating are studied and compared from an energy consumption perspective.Afterward, the environmental impact of each mode of the preheating stage is calculated.The results are discussed subsequently, and some important conclusions are drawn.

Goal and Scope
The intended application of this study is to support composite manufacturing teams.The purpose for conducting this study is to help guide the selection of preheating processes towards reducing production environmental impact through energy consumption savings and emissions reduction.
The function of this stage in the thermoforming process is to heat the composite part to make it deformable.The significance of this particular stage of the process surpasses that of other stages in thermoplastic composite manufacturing.Hence, it must be analyzed to shed light on its environmental impacts.
In this study, the functional unit is defined as the thermal processing of a flat thermoplastic composite with dimensions of 1 × 1 × 0.001 m within a 10 s timeframe.This choice is made as a reasonable reference point to which the inventory data are related to ensure all modes of preheating are compared on a common basis.
The reference flow considered in this research is the amount of energy required to fulfill the functional unit.This study utilizes the engineering analysis method to determine the life cycle inventory for each mode of preheating.Therefore, the reference flow in this study is established through theoretical calculations.As part of the scope determination in the LCA, the identification of the initial system boundaries is another important step.For this study, a partial gate-to-gate LCA is adopted, as illustrated in Figure 1.The current research is structured as follows: first, the goal and scope of this study are discussed.Subsequently, inventory analysis will be conducted, in which the hybrid analytical-numerical method is put into practice to determine the energy consumption of the preheating stage of the thermoforming process.Three different modes of preheating are studied and compared from an energy consumption perspective.Afterward, the environmental impact of each mode of the preheating stage is calculated.The results are discussed subsequently, and some important conclusions are drawn.

Goal and Scope
The intended application of this study is to support composite manufacturing teams.The purpose for conducting this study is to help guide the selection of preheating processes towards reducing production environmental impact through energy consumption savings and emissions reduction.
The function of this stage in the thermoforming process is to heat the composite part to make it deformable.The significance of this particular stage of the process surpasses that of other stages in thermoplastic composite manufacturing.Hence, it must be analyzed to shed light on its environmental impacts.
In this study, the functional unit is defined as the thermal processing of a flat thermoplastic composite with dimensions of 1 × 1 × 0.001 m within a 10 s timeframe.This choice is made as a reasonable reference point to which the inventory data are related to ensure all modes of preheating are compared on a common basis.
The reference flow considered in this research is the amount of energy required to fulfill the functional unit.This study utilizes the engineering analysis method to determine the life cycle inventory for each mode of preheating.Therefore, the reference flow in this study is established through theoretical calculations.As part of the scope determination in the LCA, the identification of the initial system boundaries is another important step.For this study, a partial gate-to-gate LCA is adopted, as illustrated in Figure 1.It is assumed that the preheating process is the dominant unit in the thermoplastic composite processing affecting the environmental impact, and thus, the contributions of other processes are not taken into account at this stage.As previously mentioned, the sole input considered in this analysis is the energy required to heat a single composite laminate with dimensions of 1 × 1 × 0.001 m, which represents the process output.Please note that the energy needed to produce the composite material is not factored into this analysis to isolate the environmental impact of the preheating mode.While sixteen impact categories associated with the preheating process are reported in Section 4, the primary focus is on climate change.This analysis employs a midpoint impact analysis.

Inventory Analysis
As explained in the Goal and Scope section of the paper, the life cycle inventory analysis in this study relies on the engineering analysis method, in which theoretical calculations are employed.To this end, the theoretical modeling of the preheating modes is first carried out.Subsequently, the theoretical results are validated using the published data.Finally, the energy consumption for each mode of preheating is calculated.

Thermal Simulation of Preheating Stage
The analysis presented in this section focuses on the thermal simulation of the preheating stage of the thermoforming process by employing the Finite Difference (FD) method with MATLAB implementation.The investigation aims to gain insights into the temperature distribution within the composite part during the preheating process.Three distinct preheating modes are studied, namely, the conductive, convective, and radiative modes.Figure 2 illustrates the schematic representation of the preheating process for each mode.While heat transfer within the composite part mainly occurs through conduction, the method by which heat is supplied to the composite part varies based on the selected preheating mode, necessitating a comprehensive analysis of each scenario.
It is assumed that the preheating process is the dominant unit in the thermoplastic composite processing affecting the environmental impact, and thus, the contributions of other processes are not taken into account at this stage.As previously mentioned, the sole input considered in this analysis is the energy required to heat a single composite laminate with dimensions of 1 × 1 × 0.001 m, which represents the process output.Please note that the energy needed to produce the composite material is not factored into this analysis to isolate the environmental impact of the preheating mode.While sixteen impact categories associated with the preheating process are reported in Section 4, the primary focus is on climate change.This analysis employs a midpoint impact analysis.

Inventory Analysis
As explained in the Goal and Scope section of the paper, the life cycle inventory analysis in this study relies on the engineering analysis method, in which theoretical calculations are employed.To this end, the theoretical modeling of the preheating modes is first carried out.Subsequently, the theoretical results are validated using the published data.Finally, the energy consumption for each mode of preheating is calculated.

Thermal Simulation of Preheating Stage
The analysis presented in this section focuses on the thermal simulation of the preheating stage of the thermoforming process by employing the Finite Difference (FD) method with MATLAB implementation.The investigation aims to gain insights into the temperature distribution within the composite part during the preheating process.Three distinct preheating modes are studied, namely, the conductive, convective, and radiative modes.Figure 2 illustrates the schematic representation of the preheating process for each mode.While heat transfer within the composite part mainly occurs through conduction, the method by which heat is supplied to the composite part varies based on the selected preheating mode, necessitating a comprehensive analysis of each scenario.

Analysis of Boundary Nodes
The preheating process can be conducted through different modes of heat transfer, each characterized by distinct formulations.This section explores the various preheating modes and provides a detailed explanation of their corresponding boundary conditions.Additionally, explicit finite difference (FD) formulations for these preheating modes are derived to aid in accurate and efficient numerical simulations.The explicit method is selected because it is easier to use, especially in computer coding [23,24].To facilitate the analyses presented in the subsequent sections, an essential assumption is made, suggesting a uniform temperature distribution across the surfaces (boundaries) of the composite laminate.This assumption is considered reasonable, particularly when the thermal elements of the preheating stage have been activated for a sufficient duration, leading to the attainment of a steady-state temperature throughout the boundaries.Therefore, no temperature gradients are induced in-plane within the composite laminate.As a result, the analysis simplifies to the calculation of a 1-D temperature distribution solely through the laminate thickness.While this assumption may impact the accuracy of the results, the comparative nature of this study, which focuses on evaluating the relative performance of different preheating modes, would further justify its use.

Boundary Nodes of the Conductive Mode of Preheating
To develop the formulation of the interface boundary condition for this mode of preheating, a control element, as illustrated in Figure 3, is introduced around the upper interface of the laminate and the hot plate.The relationship of the heat transfer rates entering the element and the rate of energy generated and stored within the element is achieved by following the framework outlined in Equation (1).
where q Top , q Bot , .
q, ρ c , C p c , ρ HP−Top , C p HP−Top , ∆t, T, i, and ∆V(= A∆z) arethe rate of heat transfer from the top hot plate to the control element, the rate of heat transfer from the composite laminate to the control element, the rate of heat generated in the control element, the density of the composite laminate, the specific heat capacity of the composite laminate, the density of the hot plate, the specific heat capacity of the hot plate, the time interval, the nodal temperature, the time step counter, and the volume of the control element (i.e., the surface area of the interface multiplied by the space interval), respectively.Expanding Equation (1) further using an explicit finite difference formulation and considering the thermoplastic nature of the polymer matrix in the composite laminate (resulting in  = 0) leads to a relationship described in Equation (2): An assumption is introduced, where the temperatures of nodes M + 1 and M − 1 are higher than that of node M, facilitating heat transfer solely to the control element.
Expanding Equation (1) further using an explicit finite difference formulation and considering the thermoplastic nature of the polymer matrix in the composite laminate (resulting in .q = 0) leads to a relationship described in Equation (2): here, k HP−Top and k c denote thermal conductivity of the hot plate and composite laminate, respectively.Let us assume ρ c C p c + ρ HP−Top C p HP−Top = β, ά = k c β , and τ = ά ∆t ∆z 2 .By isolating T i+1 M , we obtain an expression elucidating the temperature at the interface node: Similarly, an analogous equation is derived for the lower interface: The aforementioned equations were derived based on the assumption of having a perfect interface, which does not hold true in the conductive preheating stage.To take this imperfection into account, a thermal contact resistance is introduced, defined as follows: where R Contact and h Contact are thermal contact resistance and thermal contact conductance, respectively.To incorporate this thermal contact resistance into our analysis, an approach involving an additional equivalent composite laminate thickness is adopted [25].This involves the introduction of an additional composite laminate thickness (l Add−Com−Thick ) at each interface, the magnitude of which is determined by Equation (6): This consideration accounts for the imperfect nature of the interface between the hot plates and the composite laminate [26].Figure 4 shows the added layer of composite which accounts for the contact resistance at the hot plate and composite laminate interface.
While the explicit method offers simplicity in comparison to the implicit method, its applicability is constrained by inherent time stability issues.Consequently, an upper threshold for the time step (∆t) must be established to ensure stability.To determine this upper limit, the stability criterion is put into practice, requiring that all coefficients of T i m in the T i+1 m equation exceed or equal zero [25].This criterion for the conductive mode of preheating leads to: While the explicit method offers simplicity in comparison to the implicit method, its applicability is constrained by inherent time stability issues.Consequently, an upper threshold for the time step (∆) must be established to ensure stability.To determine this upper limit, the stability criterion is put into practice, requiring that all coefficients of  in the  equation exceed or equal zero [25].This criterion for the conductive mode of preheating leads to:

Boundary Nodes of the Convective Mode of Preheating
To model this boundary condition, a control element, as illustrated in Figure 5, is introduced at the upper boundary of the laminate.The formulation of the heat transfer rates entering the element and the rate of energy stored within the element is achieved by following the framework outlined in Equation (8).
In this context,  ,  , and ∆(=  ∆ 2 ) correspond to the rate of heat transfer from the composite laminate to the control element, the rate of heat transfer from the free stream to the control element, and the volume of the control element (i.e., the laminate s surface area multiplied by half of the space interval), respectively.

Boundary Nodes of the Convective Mode of Preheating
To model this boundary condition, a control element, as illustrated in Figure 5, is introduced at the upper boundary of the laminate.The formulation of the heat transfer rates entering the element and the rate of energy stored within the element is achieved by following the framework outlined in Equation (8).Employing an explicit finite difference formulation to expand Equation (8) leads t the derivation of a resulting relationship, as shown in Equation ( 9): where ℎ and  denote convective heat transfer coefficient and the temperature of th In this context, q Bot , q Conv , and ∆V(= A ∆z 2 ) correspond to the rate of heat transfer from the composite laminate to the control element, the rate of heat transfer from the free stream to the control element, and the volume of the control element (i.e., the laminate's surface area multiplied by half of the space interval), respectively.
Employing an explicit finite difference formulation to expand Equation (8) leads to the derivation of a resulting relationship, as shown in Equation ( 9): where h and T ∞ denote convective heat transfer coefficient and the temperature of the free stream, respectively.Let us assume and τ = α ∆t ∆z 2 .Solving for T i+1 M yields the following expression for the temperature at the uppermost node: Similarly, an analogous equation is derived for the lowermost node: Note that in this scenario, we assume the laminate's periphery is secured within a frame-style blank holder system, ensuring that both surfaces of the laminate remain exposed to hot air.To establish the upper threshold of the time step suitable for this particular setup, we again employ the stability criterion, the result of which is summarized in Equation (12).

Boundary Nodes of the Radiative Mode of Preheating
To determine the temperature profile within the composite laminate during the radiative mode, a control element, depicted in Figure 6, is introduced at the upper boundary of the laminate, allowing us to formulate a heat transfer rate balance equation, as detailed in Equation (13).
here, q Bot , q Rad−Abs−ex , and ∆V(= A ∆z 2 ) represent the rate of heat transfer from the composite laminate to the control element, the rate of absorbed heat transfer from the radiative source to the boundary control element, and the volume of the control element (i.e., the laminate's surface area multiplied by half of the space interval), respectively.
Utilizing an explicit finite difference formulation to expand Equation ( 13) yields a resultant relationship, as described in Equation ( 14): where ε e , σ, F, and B ex denote, in turn, effective emissivity, Stefan-Boltzmann constant, view factor, and the absorbed fraction of the radiation striking the boundary layer of the composite laminate.Effective emissivity in Equation ( 14) is calculated as follows: where ε IR Heater and ε c are the emissivity of the infrared heater and the composite laminate, respectively [27].The absorbed fraction of the radiation striking the boundary layer of the composite laminate is calculated using the Beer-Lambert law [28]: where Á is the absorption coefficient of the material.By isolating T i+1 M , we obtain an expression elucidating the temperature at the uppermost node:  Utilizing an explicit finite difference formulation to expand Equation ( 13) yield resultant relationship, as described in Equation ( 14): where  ,  ,  , and  denote, in turn, effective emissivity, Stefan-Boltzmann c stant, view factor, and the absorbed fraction of the radiation striking the boundary la of the composite laminate.Effective emissivity in Equation ( 14) is calculated as follow where  and  are the emissivity of the infrared heater and the composite la inate, respectively [27].The absorbed fraction of the radiation striking the boundary la of the composite laminate is calculated using the Beer-Lambert law [28]: where  is the absorption coefficient of the material.By isolating  , we obtain an pression elucidating the temperature at the uppermost node: Similarly, an analogous equation is derived for the lowermost node: Similarly, an analogous equation is derived for the lowermost node: IR Heater (18) Applying the Beer-Lambert law to fiber-reinforced polymer composites is more complex than applying it to pure polymers.This complexity arises from the heterogenous nature of composites, which consist of fibers and polymer matrix.To use the Beer-Lambert law in composites, special considerations or assumptions shall be made.One such assumption is the homogeneity of the composite material, which allows for the definition of an effective absorption coefficient.This coefficient represents an average value that accounts for the combined effects of both the polymer and fiber.Due to intricate structure of composites, this effective absorption coefficient is typically determined empirically.
The determination of the upper limit for the time step (∆t) within this specific boundary condition mode relies on an iterative trial-and-error process to ensure the stability of the finite difference model.While in other modes of preheating, the interior nodes may not be directly affected by the boundary conditions, the dynamics of radiative heat transfer make it distinct.Radiative heat transfer stands apart from other mechanisms due to its unique ability to transmit through boundary layers and penetrate into lower levels, thus influencing the nodal temperatures of the interior nodes directly.Consequently, when computing the temperature profile of these interior nodes, it is important to account for the radiation energy transferred via the boundary layers.This concept will be elaborated upon in subsequent sections.

Analysis of Interior Nodes
While heat transfer from the environment to the composite laminate could occur through one or a combination of the mechanisms elucidated in Section 3.1, the primary mode of heat transfer within the laminate is through conduction.This statement holds true for both conductive and convective boundary conditions.However, in the case of radiation, an additional mechanism of heat transfer emerges within the composite laminate that needs to be considered in our analysis.Subsequent sections delve into a comprehensive explanation of heat transfer mechanisms pertaining to the interior nodes.

Interior Nodes of the Conductive and Convective Modes of Preheating
The primary mechanism of heat transfer within the composite laminate is conduction, which is governed by the partial differential equation, described in Equation (19).
Similarly to the boundary cases, we use an explicit finite difference formulation to solve this particular partial differential equation.Equation (20) outlines the explicit finite difference formulation applicable to the interior nodes of the composite laminate during the conductive and convective modes of preheating: Interior Nodes of the Radiative Modes of Preheating When analyzing the temperature profile of the interior nodes within the composite laminate under the radiative preheating mode, an additional energy transfer term comes into play.This term accounts for the rate of radiation transmitted from the boundary layers to the interior layers.In the radiative preheating mode, as the heating element emits energy towards the composite surface, the majority of this energy is absorbed by the laminate (typically around 95% for a standard composite sheet [28]), with a smaller fraction being reflected.Among this non-reflected energy, a portion is absorbed by the laminate itself, while the remainder transmits through the material.
As previously detailed in Section 3.1.3,a fraction of the radiative energy gets absorbed by the boundary layer, characterized by its absorption rate q Rad−Abs−ex .The rest of the energy is conveyed to the interior layers, where a fraction is absorbed by these layers.Additionally, a portion of the radiative energy, unabsorbed by either the boundary or interior layers, exits the laminate.By considering this energy transfer narrative, the heat transfer rate balance equation for the interior nodes can be expressed as follows (see Figure 7): In the equation, q Top , q Bot , q Rad−Abs−Top−in , q Rad−Abs−Bot−in denote the respective rates of conductive heat transfer from node m + 1 to the control element and from node m − 1 to the control element, and the respective radiative heat transfer rates from the top boundary layer to the control element and from the bottom boundary layer to the control element.
The expansion of Equation ( 21) using an explicit finite difference formulation leads to the derivation of a resultant relationship, as elaborated upon in Equation ( 22).
q Rad−Abs−Top−in m and q Rad−Abs−Bot−in m are calculated as follows: Additionally, a portion of the radiative energy, unabsorbed by either the boundary or inte rior layers, exits the laminate.By considering this energy transfer narrative, the heat transfe rate balance equation for the interior nodes can be expressed as follows (see Figure 7): In the equation,  ,  ,  ,  denote the respective rates of conductive heat transfer from node  + 1 to the control element and from node  − 1 to the control element, and the respective radiative heat transfer rates from the top boundary layer to the control element and from the bottom boundary layer to the contro element.Using the Beer-Lambert law [28], the nodal absorbed fraction of the radiation striking the composite laminate by the interior layers is calculated as follows: where Á∆z) .Let us assume that B m−Total = B m−Top + B m−Bot .By isolating T i+1 m in Equation ( 22), we obtain an expression elucidating the temperature at the interior nodes of the composite laminate exposed to the radiative boundary condition:

Finite Difference Model Validation
In the process of validating the finite difference formulations developed for the preheating modes, the support of experimental data from existing literature was leveraged [24,29].MATLAB software (R2023a) was employed to simulate thermal phenomena during the preheating stage by using the finite difference formulations derived in this study.The validation process was first carried out for both the conductive and convective preheating modes through a comparison of results with those previously presented in [24].The similarity in configurations between our study and the referenced results facilitated this validation.
Table 1 summarizes the key parameters and their corresponding values (extracted from the reference data) used for modeling the conductive heat transfer mode of APC-2/S4 (PEEK-Carbon Fiber Reinforced Prepreg).Figure 8 presents a comparison between the simulation results and the experimental data as presented in [24].The figure demonstrates a reasonable alignment between the experimental data and the numerical results derived from the finite difference model developed in MATLAB for the midpoint of the composite laminate.An identical approach was undertaken to validate the MATLAB results pertaining to the convective preheating mode.Table 2 provides an overview of the parameters and corresponding values adopted for modeling the convective heat transfer mode.Displayed in Figure 9 is a juxtaposition of simulation results and experimental data presented in [24] for the convective preheating mode.This figure highlights a favorable agreement between the experimental data and the numerical results derived through the MATLAB finite difference model for the midpoint of the composite laminate.An identical approach was undertaken to validate the MATLAB results pertaining to the convective preheating mode.Table 2 provides an overview of the parameters and corresponding values adopted for modeling the convective heat transfer mode.Displayed in Figure 9 is a juxtaposition of simulation results and experimental data presented in [24] for the convective preheating mode.This figure highlights a favorable agreement between the experimental data and the numerical results derived through the MATLAB finite difference model for the midpoint of the composite laminate.
While [24] provided a reasonable basis for the experimental investigations of the conductive and convective preheating modes, the research presented in [29] is notable for its depth and comprehensiveness in studying the radiative preheating mode.Thus, the experimental investigation presented in [29] is selected to validate the MATLAB finite difference code.Table 3 outlines the parameters and their corresponding values (extracted from the reference data) employed for modeling the radiative heat transfer mode.As depicted in Figure 10, the numerical model s temperature predictions for the surface point show a satisfactory alignment with the measured experimental data reported in [29], for the time interval after 25 s.A possible rationale for the discrepancy between numerical and experimental results during the first 25 s could be attributed to an assumption made in the MATLAB simulation-namely, the assumption of a steady radiation source right from the beginning of the process.In practice, this may not hold true, as the radiative heater might require a certain duration to attain a stable energy emission towards the composite laminate.While [24] provided a reasonable basis for the experimental investigations of the conductive and convective preheating modes, the research presented in [29] is notable for its depth and comprehensiveness in studying the radiative preheating mode.Thus, the experimental investigation presented in [29] is selected to validate the MATLAB finite difference code.Table 3 outlines the parameters and their corresponding values (extracted from the reference data) employed for modeling the radiative heat transfer mode.As depicted in Figure 10, the numerical model's temperature predictions for the surface point show a satisfactory alignment with the measured experimental data reported in [29], for the time interval after 25 s.A possible rationale for the discrepancy between numerical and experimental results during the first 25 s could be attributed to an assumption made in the MATLAB simulation-namely, the assumption of a steady radiation source right from the beginning of the process.In practice, this may not hold true, as the radiative heater might require a certain duration to attain a stable energy emission towards the composite laminate.

Material
Analyzing the environmental impact of each preheating mode necessitates the calculations of the energy consumption.This study considers a flat 8-harness satin woven glassfiber-reinforced polyetherimide (GF/PEI) composite produced by Toray with dimensions of 1 × 1 × 0.001 m.The properties and the processing temperature of this material were reported in [29].This specific composite laminate is chosen to serve as a representative for the blanks employed in composite thermoforming processes.The properties of the GF/PEI composite laminate used to calculate the energy consumption during the preheating phase are summarized in Table 4.

Material
Analyzing the environmental impact of each preheating mode necessitates the calculations of the energy consumption.This study considers a flat 8-harness satin woven glass-fiber-reinforced polyetherimide (GF/PEI) composite produced by Toray with dimensions of 1 × 1 × 0.001 m.The properties and the processing temperature of this material were reported in [29].This specific composite laminate is chosen to serve as a representative for the blanks employed in composite thermoforming processes.The properties of the GF/PEI composite laminate used to calculate the energy consumption during the preheating phase are summarized in Table 4.

Analysis of Energy Consumption of Modes of Preheating
The analysis of energy consumption for the preheating modes requires us to adjust their process parameters to ensure all modes of preheating are compared on the same grounds.To elucidate, the processing parameters for each preheating mode were adjusted to align their production rates.With a processing temperature of 350 • C for the GF/PEI composite laminate, MATLAB model parameters were tuned to achieve a midpoint tem-perature of 350 • C within 10 s.This 10 s interval solely accounts for the net heating time, excluding process setup.

Energy Consumption of the Conductive Mode of Preheating
The simulation parameters to obtain a midpoint temperature of 350 • C in 10 s were tailored in the MATLAB finite difference code as detailed in Table 5. Figure 11 represents the temperature profile of the composite laminate corresponding to the conductive preheating mode with the tailored parameters.

Analysis of Energy Consumption of Modes of Preheating
The analysis of energy consumption for the preheating modes requires us to adjust their process parameters to ensure all modes of preheating are compared on the same grounds.To elucidate, the processing parameters for each preheating mode were adjusted to align their production rates.With a processing temperature of 350 °C for the GF/PEI composite laminate, MATLAB model parameters were tuned to achieve a midpoint temperature of 350 °C within 10 s.This 10 s interval solely accounts for the net heating time, excluding process setup.

Energy Consumption of the Conductive Mode of Preheating
The simulation parameters to obtain a midpoint temperature of 350 °C in 10 s were tailored in the MATLAB finite difference code as detailed in Table 5. Figure 11 represents the temperature profile of the composite laminate corresponding to the conductive preheating mode with the tailored parameters.These parameters are used to calculate the energy required to fulfill one functional unit.To determine this energy, the following relationship is used based on the heat transfer mechanisms shown in Figure 12: q in−HP = q out−Conv + q out−Rad + q out−Laminate (28) where q in−HP , q out−Conv , q out−Rad , and q out−Laminate are, in turn, the rate of heat transfer entering the hot plate, the rate of heat transfer existing the hot plate due to natural convection, the rate of heat transfer existing the hot plate due to radiation, and the rate of heat transfer exiting the hot plate to heat up the composite laminate.
These parameters are used to calculate the energy required to fulfill one functional unit.To determine this energy, the following relationship is used based on the heat transfer mechanisms shown in Figure 12: where  ,  ,  , and  are, in turn, the rate of heat transfer entering the hot plate, the rate of heat transfer existing the hot plate due to natural convection, the rate of heat transfer existing the hot plate due to radiation, and the rate of heat transfer exiting the hot plate to heat up the composite laminate.This analysis assumes that the hot plates temperature remains constant during preheating.Equation ( 28) can be further expanded to: here, ℎ ,  ,  ,  , and  denote the coefficient of natural convection heat transfer, the ambient temperature, the coefficient of emissivity of the aluminum hot plate, the final temperature of the laminate, and the initial temperature of the laminate, respectively.Note that in this analysis, the startup energy required to heat up the plates has been disregarded, as this energy, when divided by the considerable number of parts processed daily, becomes negligible.
For calculating the energy required per functional unit,  is multiplied by ∆ (10 s) and then by 2 to account for two hot plates in the conductive preheating mode.Using the summarized values in Table 6, the energy consumption per functional unit is 0.17 kWh.This analysis assumes that the hot plates' temperature remains constant during preheating.Equation ( 28) can be further expanded to: here, h n , T ∞−n , ε HP , T c−Final , and T c−Initial denote the coefficient of natural convection heat transfer, the ambient temperature, the coefficient of emissivity of the aluminum hot plate, the final temperature of the laminate, and the initial temperature of the laminate, respectively.Note that in this analysis, the startup energy required to heat up the plates has been disregarded, as this energy, when divided by the considerable number of parts processed daily, becomes negligible.For calculating the energy required per functional unit, q in−HP is multiplied by ∆t (10 s) and then by 2 to account for two hot plates in the conductive preheating mode.Using the summarized values in Table 6, the energy consumption per functional unit is 0.17 kWh.

Energy Consumption of the Convective Mode of Preheating
The approach for determining the energy consumption of the convective preheating mode closely mirrors that of conductive preheating.Parameters derived from MATLAB are summarized in Table 7, achieving a 350 • C midpoint temperature within 40 s.For analysis, it is assumed that four parts are processed simultaneously in the oven, rendering an average processing time of 10 s for each laminate, consistent with other preheating modes.Figure 13 demonstrates the temperature profile of the composite laminate in convective preheating mode using the parameters summarized in Table 7. Units N.A.

Energy Consumption of the Convective Mode of Preheating
The approach for determining the energy consumption of the convective preheating mode closely mirrors that of conductive preheating.Parameters derived from MATLAB are summarized in Table 7, achieving a 350 °C midpoint temperature within 40 s.For analysis, it is assumed that four parts are processed simultaneously in the oven, rendering an average processing time of 10 s for each laminate, consistent with other preheating modes.Figure 13 demonstrates the temperature profile of the composite laminate in convective preheating mode using the parameters summarized in Table 7.To determine the energy required to fulfill one functional unit in a convection oven (Q in−Oven ), the following terms of energy are considered in the analysis: where Q out−Laminate , Q out−Loss , and Q out−Fan represent the energy needed to heat up the laminate, the heat (energy) loss occurring in the oven during the process, and the energy required by the fans to generate the necessary air velocity for the hot air (T ∞ = 459 °C) to obtain h = 30 W m 2 • • C , respectively.Generally, the average air velocity governs the heat transfer capability of the convective oven [30].Heat transfer correlations have been formulated and can be employed to estimate the range of convective heat transfer coefficient as a function of the average air velocity and pressure [31,32].For turbulent flow over a flat plate, the subsequent correlation has been proposed: here, l c is a characteristic flow length, k Air is the thermal conductivity of hot air, and Re and Pr are the Reynolds number and Prandtl number, respectively.The latter two dimensionless numbers are defined as follows: where ρ Air , v Air , and µ Air denote the density of hot air, the hot air velocity, and the viscosity of hot air, respectively.By employing Equations ( 31)-( 33), the average velocity of the hot air (v Air ) required to obtain a specific convective heat transfer coefficient (h) is determined using the following equation: On the other hand, the airflow velocity generated by an axial fan, v Fan , can be calculated using the equation below [30]: where P Fan and A Fan are, in turn, the power of the fan's motor and the cross-sectional flow area of the fan.Using the continuity equation, where v Air A Oven = v Fan A Fan , the following equation is derived: here, A Oven is the cross-sectional area of the oven.Now, we can expand Equation (30) further to calculate the energy required to fulfill one functional unit in the convective mode of preheating: Similarly, in these calculations, the start-up energy required to raise the oven temperature to the processing level is neglected, as this energy, when divided by the large number of parts processed per day, becomes negligible.The energy needed to fulfill one functional unit is calculated using the values outlined in Table 8, resulting in a total of 1.79 kWh.To ascertain the energy consumption of the radiative energy source, the heat transfer balance for a unit of the radiation source fulfilling one functional unit is determined using the following equation: where A IR Heater represents the surface area of the infrared heater and ε IR Heater is the coefficient of emissivity of the infrared heater.The energy required to fulfill a single functional unit can be calculated by multiplying q in−Rad by the processing time (∆t), which is 10 s.Then, this product is multiplied by 2 due to the presence of two infrared heaters in the radiative preheating mode.Utilizing the values summarized in Table 10, the energy needed to fulfill one functional unit amounts to 0.46 kWh.

Impact Assessment of Modes of Preheating
Life Cycle Impact Assessment (LCIA) is a critical step in evaluating the environmental impact of the preheating modes based on the theoretical energy consumption results obtained in Section 3.2.OpenLCA 2.0.2 software, integrated with the Product Environmental Footprints (PEF) database, is employed in this study to determine how the amount of energy required to fulfill one functional unit using different modes of preheating impacts the environment.Environmental footprints originated from the European Commission's Single Market for Green Products initiative, which aims to develop a common methodology for the quantitative assessment of environmental impacts of products and processes to support their assessment and labeling [33].
The impact assessment method employed in this database relies on the PEF Mid-Point Indicator methodology [34].As outlined in Table 11, PEF includes sixteen impact categories at the midpoint level, which are determined for the different modes of preheating based on their corresponding energy consumptions.Note that for this study, the impact assessment of each preheating mode is conducted by considering two basic processes to generate the reference flow: the medium voltage grid-mix electricity generation and delivery process (1-60 kV), as well as the conservation from medium voltage to low voltage process (240-600 V) that is commonly used by the majority of the industry sectors employing the thermoforming process.

Results and Discussion
The global market size of thermoplastic composites was estimated at USD 32.21 billion in 2021 and is expected to increase to USD 62.62 billion in 2030 [35].The demand for thermoplastic composites is on the rise on account of their high impact resistance, unique processing possibilities, lightweight properties, strength, and environmental advantages [36].When converting the estimated dollar value to the estimated volume of thermoplastic composites used in the market, it is projected that over 1 million cubic meters of thermoplastic composites will be used annually.This volume can be hypothetically equated to 1 billion units of the relatively thin 1 × 1 × 0.001 m laminates investigated in this study.
Since thermal processing is the dominant method used for thermoplastic composites before shaping, a straightforward upper-lower bound calculation suggests that a potential reduction of 1.45 billion kg of CO 2 per year is achievable if conductive preheating is chosen over the convective mode of preheating for relatively thin laminates in the global thermoplastic composite industry.This type of theoretical analysis underscores the importance of the decisions made by manufacturing leaders in incorporating environmental impacts into the design of products and their manufacturing processes.
While the radiative mode of preheating is the most-commonly employed technique for heating thermoplastic composites [37], a comparison of energy consumption between the radiative and conductive modes reveals room for improvement, with a well-considered choice between these two preheating methods.According to Hertwich [38], "Sustainable solutions, even if they do not achieve universal market penetration, are important because they indicate a development alternative, influence the development of conventional products, and may be adopted in situations of crisis".Therefore, it is evident that even small improvements in the manufacturing processes of thermoplastic composites can be significant steps toward creating a more sustainable world.This is further emphasized when considering Figure 15, which schematically illustrates the relationship between effects and emissions.
they indicate a development alternative, influence the development of conventional products, and may be adopted in situations of crisis".Therefore, it is evident that even small improvements in the manufacturing processes of thermoplastic composites can be significant steps toward creating a more sustainable world.This is further emphasized when considering Figure 15, which schematically illustrates the relationship between effects and emissions.As depicted in the graph, the impact on the environment varies depending on the amount of emissions.When emissions are minimal, the effect on the environment may be negligible.However, as the accumulated emissions increase over several decades or centuries, the adverse impact becomes substantially larger [39].Hence, regardless of our As depicted in the graph, the impact on the environment varies depending on the amount of emissions.When emissions are minimal, the effect on the environment may be negligible.However, as the accumulated emissions increase over several decades or centuries, the adverse impact becomes substantially larger [39].Hence, regardless of our current level of emissions, we should strive to minimize the emissions from our processes to delay the significant adverse environmental consequences of our actions.This study demonstrates the potential for notable reductions in environmental impacts through the thermal processing of thermoplastic composites.
In light of the preceding discussion, it is of utmost importance to develop practical techniques that harness the benefit of conductive preheating in thermoplastic composite processing, particularly for relatively thin to medium-thickness laminates.Although there are limitations associated with thermal processing using the conductive preheating method [40], such as the risk of composite material sticking to hot plates, the substantial difference in environmental impact between these methods justifies further research to address these challenges.Companies engaged in thermal processing of thermoplastic composites are strongly encouraged to proactively explore and address the technical challenges associated with adopting the conductive preheating mode, particularly when processing relatively thin to medium-thickness composite laminates.
On another note, when comparing the environmental impact of material production for thermoplastic composites, specifically PEI and fiberglass in this case, with the various preheating methods, it becomes evident that the preheating process can consume energy equivalent to as much as 13 percent of the material production (i.e., 14.17 kg of CO 2 emissions for material production [33] compared to 1.79 kg of CO 2 emissions when using convective preheating for the relatively thin laminate cases).This comparison also highlights the significant role of the preheating process throughout the entire life cycle of thermoplastic composite parts, which must be taken into serious consideration.
Viewing this study from another perspective helps us recognize how process parameters can be adjusted to reduce the environmental impacts of each process.Historically, most optimization efforts have been primarily focused on enhancing process productivity, reducing production cost, and improving final product quality [41].The environmental impact perspective received comparatively less attention or was often overlooked.
An analysis of the conductive preheating process reveals the interrelation between energy consumption and, consequently, the environmental impact with various process parameters.Figure 16a presents a spider diagram that illustrates the influence of process parameters on energy consumption of the conductive preheating mode.The graph shows the sensitivity of energy consumption, and consequently, environmental impact, to changes in process parameters.
In this study, One-Factor-At-A-Time (OFAT) sensitivity analysis is employed, in which a Sensitivity Index (SI) is defined as follows: where Γ, ∆Γ, ϱ, and ∆ϱ are the output parameter at the base value of the input, the change in the output parameter, the base value of the input parameter, and the change in the input parameter from the base value, respectively.This approach is a straightforward method for assessing the impact of individual input parameters on process output.As demonstrated in Figure 16a, the most influential process parameters in the conductive mode of preheating are the temperature of the hot plate and the interface surface property.This highlights the importance of optimizing these process parameters to minimize energy consumption and environmental impacts in this preheating mode.
The same analysis was conducted for the convective mode of preheating, as illustrated in Figure 16b.The sensitivity analysis reveals that the velocity of the hot fluid and its viscosity have the most significant influence on energy consumption.While other parameters do affect energy consumption, these two parameters have a notably greater impact.It is worth noting that these process parameters are interrelated with the temperature of the oven; therefore, an integrated process parameter optimization approach should be considered for realistic and practical results.
In the sensitivity analysis of the radiative mode of preheating as shown in Figure 16c, it is evident that the view factor has the most pronounced impact on energy consumption in this mode.Hence, to reduce energy consumption in this process, the radiation set-up should be designed to maximize the view factor, bringing it close to unity.This theoretical study lays the foundation for a parametric optimization process with the primary goal of minimizing energy consumption at this stage of the manufacturing process and its associated environmental impacts.
While this theoretical investigation highlights the importance of selecting the conductive mode of preheating for relatively thin composite laminates, it is crucial to note that the choice of the most efficient preheating method also depends on the thickness of the laminate.As shown in Figure 17a, the trend of energy consumption for different preheating modes would vary with changes in laminate thickness.The analysis of this trend reveals that the conductive mode is the most efficient method for preheating composite laminates within a relatively thin to medium-thickness range (i.e., 1-10 mm thick laminates).However, as the laminate thickness increases, the convective and radiative modes of preheating become more competitive and even more efficient, as depicted in Figure 17b, where the normalized energy consumption for laminates with different thicknesses referenced to the conductive mode is shown.Therefore, it is imperative to determine the most efficient mode of preheating for thermoplastic composite laminates based on their material and geometrical properties to ensure minimal environmental impacts associated with the manufacturing process selection.where the normalized energy consumption for laminates with different thicknesses referenced to the conductive mode is shown.Therefore, it is imperative to determine the most efficient mode of preheating for thermoplastic composite laminates based on their material and geometrical properties to ensure minimal environmental impacts associated with the manufacturing process selection.

Conclusions
In the current study, the preheating stage of the thermoforming process for thermoplastic composites was theoretically modeled to investigate its environmental impact.The analysis revealed that the preheating process significantly contributes to the life cycle assessment of thermoplastic composite parts.To elucidate, the environmental impact of the preheating process using the convection method for relatively thin laminates amounts to as much as 13 percent of the composite material production.Furthermore, the study demonstrated that conductive preheating has the lowest energy consumption (i.e., 0.17 kWh) and environmental impact for relatively thin laminates, while the convective preheating mode exhibits the highest (i.e., 1.79 kWh).This substantial difference underscores the importance of considering the environmental impact of these manufacturing methods when designing manufacturing processes for the thermoplastic composites.The radiative mode of preheating, which is the most common technique in thermal processing of thermoplastic composites, falls in the middle position in terms of energy consumption and environmental impact for relatively thin laminates.However, the difference with the conductive preheating method remains significant (i.e., 2.7 times larger).Additionally, a sensitivity analysis was conducted, revealing how sensitive energy consumption and environmental impact are to the manufacturing process parameters.For conductive preheating, the temperature of the hot plates and their interface properties with the composite laminate were found to have the highest impact.In the case of convection preheating, fluid velocity and viscosity had the most significant impact.In radiative preheating, the view factor emerged as one of the most influential factors.It was also discussed that laminate thickness impacts the energy consumption trend in different modes of preheating, making the convective and radiative modes of preheating more competitive or even more efficient as laminate thickness increases.

Funding:
The author declares that no funds, grants, or other support were received during the preparation of this manuscript.Data Availability Statement: Data are contained within the article.

Conflicts of Interest:
The author declares no conflicts of interest.

Nomenclature q Top
Rate of heat transfer from top to the control element q Bot Rate of heat transfer from bottom to the control element

Figure 1 .
Figure 1.Identification of the initial system boundaries, representing a partial gate-to-gate LCA.Figure 1. Identification of the initial system boundaries, representing a partial gate-to-gate LCA.

Figure 1 .
Figure 1.Identification of the initial system boundaries, representing a partial gate-to-gate LCA.Figure 1. Identification of the initial system boundaries, representing a partial gate-to-gate LCA.

Figure 2 .
Figure 2. Schematic representation of the preheating process, illustrating different modes of heat transfer (conductive, radiative, and convective) in a composite laminate.The space domain of the laminate has been discretized, and a space interval of ∆z is considered.

Figure 3 .
Figure 3. Schematic representation of the interface between the hot plate and the composite laminate.An assumption is introduced, where the temperatures of nodes M + 1 and M − 1 are higher than that of node M, facilitating heat transfer solely to the control element.

Figure 3 .
Figure 3. Schematic representation of the interface between the hot plate and the composite laminate.An assumption is introduced, where the temperatures of nodes M + 1 and M − 1 are higher than that of node M, facilitating heat transfer solely to the control element.

Figure 4 .
Figure 4. Schematic representation of the equivalent additional contact resistance layer.

Figure 4 .
Figure 4. Schematic representation of the equivalent additional contact resistance layer.

3 Figure 5 .
Figure 5. Schematic representation of the boundary layer in the convective mode of heat transfe An assumption is introduced, where the temperatures of nodes M + 1 and M − 1 are higher than tha of node M, facilitating heat transfer solely to the control element.

Figure 5 .
Figure 5. Schematic representation of the boundary layer in the convective mode of heat transfer.An assumption is introduced, where the temperatures of nodes M + 1 and M − 1 are higher than that of node M, facilitating heat transfer solely to the control element.

2ε e σB ex ∆z τ k c T 4
IR Heater (17) J. Manuf.Mater.Process.2024, 8, x FOR PEER REVIEW 9 oradiative source to the boundary control element, and the volume of the control elem (i.e., the laminate s surface area multiplied by half of the space interval), respectively.

Figure 6 .
Figure 6.Schematic representation of the infrared heater and composite laminate configuration

Figure 6 .
Figure 6.Schematic representation of the infrared heater and composite laminate configuration.

Figure 7 .
Figure 7. Schematic representation of the infrared heater and composite laminate configuration fo the calculation of interior nodal temperature profile.An assumption is made in which heat only flows toward the control element.The expansion of Equation (21) using an explicit finite difference formulation lead to the derivation of a resultant relationship, as elaborated upon in Equation (22).   −  ∆ +    −  ∆ +  +  =   ∆  −  ∆ (22  and  are calculated as follows:

Figure 7 .
Figure 7. Schematic representation of the infrared heater and composite laminate configuration for the calculation of interior nodal temperature profile.An assumption is made in which heat only flows toward the control element.

Figure 10 .
Figure 10.MATLAB finite difference simulation results vs. experimental results for the radiative mode of preheating.Experimental results reprinted from Composites Part C: Open Access, 5, Nardi, D.; Sinke, J., Design analysis for thermoforming of thermoplastic composites: prediction and machine learning-based optimization, 100126, 2021.

Figure 10 .
Figure 10.MATLAB finite difference simulation results vs. experimental results for the radiative mode of preheating.Experimental results reprinted from Composites Part C: Open Access, 5, Nardi, D.; Sinke, J., Design analysis for thermoforming of thermoplastic composites: prediction and machine learning-based optimization, 100126, 2021.

Figure 11 .
Figure 11.MATLAB temperature profile prediction of the composite laminate using the parameters listed in Table 5 (Conductive mode of preheating).

Table 5 .Figure 11 .
Figure 11.MATLAB temperature profile prediction of the composite laminate using the parameters listed in Table 5 (Conductive mode of preheating).

Figure 12 .
Figure 12.The mechanism of heat transfer in the hot plate mold.

Figure 12 .
Figure 12.The mechanism of heat transfer in the hot plate mold.

Table 7 .
Simulation parameters to obtain a midpoint temperature of 350 • C in 10 s through the convective mode of preheating.

Figure 13 .
Figure 13.MATLAB temperature profile prediction of the composite laminate using the parameters listed in Table 7 (Convective mode of preheating).

Table 7 .
Simulation parameters to obtain a midpoint temperature of 350 ℃ in 10 s through the convective mode of preheating.

Figure 13 .
Figure 13.MATLAB temperature profile prediction of the composite laminate using the parameters listed in Table 7 (Convective mode of preheating).

Figure 16 .
Figure 16.Process parameter sensitivity analysis: (a) Conductive, (b) Convective, and (c) Radiative mode of preheating.Please note that to better illustrate the sensitivity of the preheating process to each parameter, the graphs are drawn based on a logarithmic modified sensitivity index approach.That is, the shown SI for each parameter is equal to Log (SI + A small positive constant offset).However, the SI values reported for each parameter on the axes of the graphs are the actual sensitivity indices.

Figure 16 .
Figure 16.Process parameter sensitivity analysis: (a) Conductive, (b) Convective, and (c) Radiative mode of preheating.Please note that to better illustrate the sensitivity of the preheating process to each parameter, the graphs are drawn based on a logarithmic modified sensitivity index approach.That is, the shown SI for each parameter is equal to Log (SI + A small positive constant offset).However, the SI values reported for each parameter on the axes of the graphs are the actual sensitivity indices.

Figure 17 .
Figure 17.Laminate thickness analysis: (a) The trend of energy consumption for different preheating modes versus laminate thickness, (b) Normalized energy consumption for laminates with different thicknesses referenced to the conductive mode.

Table 1 .
Key parameters and corresponding values utilized for the modeling of the conductive heat transfer mode.

Table 2 .
Parameters and corresponding values adopted for modeling the convective heat transfer mode.

Table 2 .
Parameters and corresponding values adopted for modeling the convective heat transfer mode.

Table 3 .
Parameters and corresponding values employed for the modeling of the radiative heat transfer mode.

Table 3 .
Parameters and corresponding values employed for the modeling of the radiative heat transfer mode.

Table 4 .
The properties of the GF/PEI composite laminate used to calculate the energy consumption during preheating.

Table 4 .
The properties of the GF/PEI composite laminate used to calculate the energy consumption during preheating.

Table 5 .
Simulation parameters to obtain a midpoint temperature of 350 • C in 10 s through the conductive mode of preheating.

Table 6 .
Parameters used in the conductive mode of preheating to calculate energy consumption per functional unit.

Table 6 .
Parameters used in the conductive mode of preheating to calculate energy consumption per functional unit.

Table 9 .
Simulation parameters to obtain a midpoint temperature of 350 • C in 10 s through the radiative mode of preheating.

Table 10 .
Parameters used in the radiative mode of preheating to calculate energy consumption per functional unit.

Table 11 .
Impact categories of different modes of preheating at the midpoint level.
Rate of heat generated in the control element Add−Com−Thick Additional composite laminate thickness at the interface q Conv Rate of heat transfer from the free stream to the control element h Convective heat transfer coefficient T ∞ Temperature of the free stream q Rad−Abs−ex Radiative absorption rate of boundary layer Respective radiative heat transfer rates from bottom boundary layer to control element q in−HP Rate of heat transfer entering the hot plate q out−Conv Rate of heat transfer existing the hot plate due to natural convection q out−Rad Rate of heat transfer existing the hot plate due to radiation q out−Laminate Rate of heat transfer exiting the hot plate to heat up the composite laminate h n Coefficient of natural convection heat transfer T ∞−n Ambient temperature ε HP Coefficient of emissivity of the aluminum hot plate T c−Final Final temperature of the laminate T c−Initial Initial temperature of the laminate Q in−Oven Energy required to fulfill one functional unit in a convection oven Q out−Laminate Energy needed to heat up the laminate Q out−Loss Heat loss occurring in the oven during the process Air Viscosity of hot air v Fan Airflow velocity generated by an axial fan P Fan Power of the fan's motor A Fan Cross-sectional flow area of the fan A Oven Cross-sectional area of the oven A IR Heater Surface area of the infrared heater Γ Output parameter at the base value of the input ∆Γ Change in the output parameter ϱ Base value of the input parameter ∆ϱ Change in the input parameter from the base value . q