Evaluating Fire Performance of Glass–Polyurethane Composite for Sustainable Cladding via Numerical and Empirical Simulation

The increased demand for cladding in high-rise buildings has prompted engineers to explore alternative products utilizing recycled materials. However, ensuring fire compliance in these alternative claddings, which are predominantly composed of low-volume polymer-based composites, poses a critical challenge. Traditional experimental methods for fire evaluation are costly, time consuming, and environmentally impactful. Considering this, a numerical approach was proposed for evaluating the fire performance of glass-polymer composite materials, which contain a high proportion of recycled glass and a lower percentage of rigid polyurethane. A cone calorimeter test was simulated using Computational Fluid Dynamics (CFD) software to investigate the flammability of the novel glass–polymer composite material. This validated numerical model was employed to assess the combustibility of the glass–polyurethane composite materials and identify influential parameters using the Design of Experiments (DoE) method. Statistical analysis revealed that three material properties, namely, the heat of combustion, the absorption coefficient, and the heat of reaction, significantly influenced the peak heat release rate (pHRR) of the glass–polyurethane composite materials compared to other properties. Based on these findings, an empirical equation was proposed that demonstrates a reasonable correlation with the pHRR of low-polymer recycled glass composite materials. The outcomes of this study hold considerable importance for understanding and predicting the combustibility behaviour of low-polymer–glass composites. By providing a validated numerical model and identifying critical material properties, this research contributes to the development of sustainable fire safety solutions for buildings, enabling the use of recycled materials and reducing reliance on conventional claddings.


Introduction
The building and construction industry has witnessed an increased demand for cladding materials that must also meet fire compliance standards.However, the resources for current cladding products are finite, necessitating the exploration of innovative approaches to recover waste products and upcycle them into high-volume cladding materials.Among the abundant waste materials, glass fines show great potential for utilization in new cladding products, wherein they will replace traditional aggregate components [1].Nevertheless, a significant challenge arises when incorporating glass fines into cladding composites, as they must be bound by combustible polymer matrices to meet structural design requirements [2][3][4][5][6][7][8].The resulting glass-polymer composite cladding must also pass rigorous fire compliance tests specific to the building and construction industry [9][10][11][12][13][14][15].However, the experimental evaluation of each potential formulation is both costly and time consuming.To overcome this challenge, it is imperative for engineers and researchers to develop computational models that can accurately predict the combustion behaviour of these glass-polymer composite claddings, enabling cost-effective and efficient engineering practices for assessing fire compliance.
The development of a reliable computational model for predicting the combustion behaviour of cladding materials offers numerous benefits to the building and construction industry.Firstly, accurate predictions facilitate improved fire safety by enabling designers and architects to make informed decisions regarding the selection of cladding materials, ultimately enhancing overall fire safety and potentially saving lives.Additionally, computational models help identify the most effective materials for preventing or slowing down the spread of flames, allowing designers and builders to select materials that meet fire safety requirements while minimising costs.Furthermore, the environmental impact of fires can be reduced by choosing cladding materials that are less likely to contribute to fire propagation, thus limiting the release of pollutants and greenhouse gases.By rapidly testing various cladding materials using computational models, the development cycle can be expedited, leading to the creation of more efficient and effective materials.Lastly, accurate predictions of cladding behaviour in fires provide building owners and managers with a better understanding of the associated risks, allowing them to take appropriate measures for risk mitigation.
In response to the demand for modelling tools that can accurately predict the combustion behaviour of engineering materials in real-life fire scenarios, various computational modelling approaches have been developed.A notable approach involves the use of numerical analysis and Computational Fluid Dynamics (CFD) models to assess the fire performance of fibre-reinforced polymer (FRP) composites in prefabricated modular units and various other engineering products prone to fire [4,[16][17][18].Nguyen et al. successfully demonstrated the validity of their CFD models by showing the similarities between the experimentally measured and numerically predicted fire reaction properties of organoclay/glass FRP composites [16].Alfakhry employed computer modelling and fire simulation tools, such as Pyrosim, FDS and Smokeview, to evaluate the fire behaviour of different materials used for external cladding [19].The objective of their study was to identify the most effective cladding material in terms of its resistance to fire propagation and ability to prevent the rapid spread of fire to upper stories.By simulating and analysing various materials, including aluminium composite panels, cement plastering, limestone, and perforated yellow bricks, researchers have aimed to provide information for the selection of cladding materials that can enhance fire safety in buildings.
Numerical simulations utilizing the CFD code Fire Dynamics Simulator (FDS) have also been used to investigate the behaviour of engineering materials in intermediate and large-scale fire tests.Incorporating Large Eddy Simulation (LES) techniques to accurately represent turbulent flows, CFD simulations focus on the performance of aluminium composite panels and polymer composite cladding in building fire scenarios [20][21][22][23].By utilizing FDS, large-scale tests were successfully simulated, providing insights into potential improvements to materials' fire performance without relying solely on physical experiments.These advancements in numerical modelling and simulation techniques offer valuable tools for understanding and enhancing the fire performance of engineering materials in practical applications.
This study contributes to the research field of fire engineering by reporting the accurate modelling of the combustion behaviour of the newly developed glass-polyurethane composite cladding.The developed models, both numerical and empirical, can be utilized to predict the flammability properties of other formulations containing the same constituent materials, thus reducing the need for extensive experimental testing.To achieve this, it was also crucial to determine critical parameters governing the combustion response, which can minimize the number of experiments required for input parameter collection.This study utilised the LES modelling technique in FDS to simulate the combustion behaviour of newly developed glass-polyurethane composite cladding, comprising 95% glass fines Polymers 2023, 15, 3635 3 of 20 and 5% polyurethane by weight.Statistical analysis based on the Design of Experiment (DoE) approach was used to determine the sensitivity of the peak heat release rate (e.g., maximum flame intensity) to the input parameters, including the material properties.Based on the findings from the statistical analysis, an empirical equation was proposed for screening newly developed cladding materials against fire safety requirements.This empirical predictive model enables the cost-effective evaluation of the combustion behaviour of glass-polyurethane composite claddings with varying polymer and glass fines content.Ultimately, the empirical model will support the rapid screening of new cladding materials for fire-threatened infrastructure, thereby contributing to the sustainable use of recycled materials in building claddings.
The high costs associated with fire compliance tests pose significant challenges for engineers aiming to introduce recycled cladding composites into the market.Traditional experimental methods for evaluating the fire performance of cladding materials are not only costly but also time consuming and environmentally unfriendly.This study sought to address the need for a more efficient and sustainable approach by developing a comprehensive numerical and empirical modelling framework.By accurately capturing the complex interactions of heat, mass, and momentum within the material, the CFD model provides valuable insights into the fire performance of claddings.By identifying the input parameters that significantly influence flame intensity through the empirical model, we can effectively streamline and speed up the screening process for new and potential cladding materials.This approach not only reduces the cost associated with fire performance testing but also saves valuable time in the evaluation of newly developed cladding materials.The outcomes of this study carry significant implications for the industry, academia, and regulatory bodies involved in building cladding design and fire safety.By understanding the combustibility behaviour of low-polymer-volume glass composites, this research paves the way for the sustainable utilization of recycled materials in cladding applications.Our findings offer engineers and architects a more informed basis for selecting and designing fire-resistant cladding systems, promoting the adoption of environmentally friendly practices and reducing reliance on conventional and finite materials.

Materials and Manufacture of Glass-Polyurethane Composite Cladding
Glass particles (or fines) used in manufacturing glass-polyurethane composite cladding were sourced from RepurposeIT (Melbourne, Australia).Glass particles with granulometric sizes of ~300 µm were obtained via the mechanical grinding of kerbside glass collections.As-received glass particles were washed under running water and oven-dried at 100 • C for 12 h.Isocyanate and polyol components from which the binder polyurethane resin was developed were supplied by Nuplex (Melbourne, Australia).Dried glass particles (95 wt.%) were compounded with isocyanate/polyol resin (5 wt.%) using a mechanical mixer operated at 100 rpm.The glass/polyurethane slurry was transferred to a 100 × 180 × 12 mm steel mould.A steel caul plate was placed above the glass/polyurethane slurry, and the mixture was subjected to an elevated temperature (160 • C) and compression pressure (6.5 MPa) moulding for 1 min [1].The resultant glass-polyurethane composite cladding was demoulded and left to cure at room temperature (23 • C) for 7 days.Cone calorimeter specimens (100 × 100 × 9 mm in thickness) were cut from the master panel using a diamond saw.Diamond-saw-cut glass-polyurethane composite cladding specimens were dried and conditioned at room temperature over 24 h prior to the cone calorimeter and other physical and thermal tests.

Cone Calorimeter Evaluation of Glass-Polyurethane Composite Cladding
A cone calorimeter was used to evaluate the fire reaction properties of glass-polyurethane composite cladding according to AS 3837 [24].Cone calorimeter specimens were wrapped in aluminium foil and placed on a steel-frame sample holder that was placed ~25 mm below the cone heater.The horizontally mounted specimen was subjected to an incident radiant heat flux of 50 kW/m 2 representative of a moderate-intensity fire.The heat release rate (HRR) and mass loss (ML) data were recorded as a function of exposure time.The HRR data were derived from oxygen consumption rates [25].The cone calorimeter test was stopped immediately after the burning specimen flamed out.

Thermal Stability and Thermal Properties of Glass-Polyurethane Composite Cladding
The thermal stability of glass-polyurethane composite cladding was evaluated using Netzsch TGA/DSC Jupiter STA 449 F5 thermogravimetric analysis (TGA) equipment (Bayern, Germany).The glass-polyurethane composite cladding was mechanically ground before ~15 mg of the powdered sample was transferred to a crucible.The crucible containing the test sample was inserted into the TGA furnace, which was continuously purged using nitrogen gas flowing at 50 mL/min.The remaining mass was measured as a function of temperature between 50 and 850 • C at a constant heating rate of 20 K/min.The reference temperature, T ref , and pyrolysis temperature range, T pyro , which are input parameters for the numerical fire simulation model, were derived from the TGA data.The reference temperature is defined as the temperature at which the peak mass loss rate is achieved during a TGA experiment.The pyrolysis range is defined by the onset and end point of a specific thermal degradation stage, as discussed later in Section 4.
The heat flow in the glass-polyurethane composite cladding was measured using the DSC131 Evo differential scanning calorimeter (DSC) according to the EN ISO 11357-1 standard [26].For the DSC experiment, ~5 mg of the test material was placed in a sample holder, which was then inserted into a DSC chamber.The chamber was heated from ambient temperature (23 • C) to 300 • C at a constant heating rate of 10 K/min.The heating segment was followed by a cooling step in which the temperature was reduced from 300 • C to ~50 • C at a constant cooling rate of 10 K/min.The heat of reaction, h r , of the glasspolyurethane composite cladding material was derived from DSC heat flow measurement.The specific heat capacity, C p , was calculated via the Rule of Mixtures (RoM) utilising input parameters acquired from the literature for glass particles [27] and polyurethane [28].
The thermal conductivity, k, of the glass-polyurethane composite cladding was measured via the transient line source method using a Thermtest TLS 100 transient thermal conductivity meter equipped with a 50 mm long rock needle sensor probe employed according to ASTM D5334-14 [29].Glass-polyurethane composite cladding specimens measuring 5 × 5 × 1 cm were used in this experiment.The sensor probe consisted of a thin heating wire and a temperature sensor, both of which were sealed in a steel tube.The heating wire applied thermal energy, while the sensor measured the sample temperature at a single point.The thermal conductivity of the glass-polyurethane composite cladding was calculated based on the heating power of the needle and the slope of the temperature-log (time) profile.

Absorption Coefficient of Glass-Polyurethane Composite Cladding
The absorption coefficient of the glass-polyurethane composite cladding was measured using a Fourier Transform infrared (FTIR) spectrometer equipped with an attenuated total reflectance (ATR) lens, as specified in ASTM E1252 [30].The recorded FTIR spectrum was an average of 8 scans collected between 4000 and 400 cm −1 at a resolution of 2 cm −1 .The emissivity, ε, was estimated using the RoM approach based on literature-derived values for polyurethane [31] and glass [27].

Numerical Simulation of Fire Reaction Properties of Glass-Polyurethane Composite Cladding
The use of numerical models to simulate the flammability properties of engineering materials offers a practical and more affordable alternative to the often expensive and labour-intensive process of physical measurement.The CFD simulation tool is valuable for fire analysis as it can provide accurate insights into the fire behaviour of materials with different compositions.However, not all CFD software products can simulate real fires, and some are better suited for fire analysis than others.Among the various categories of available CFD software, including Simula XFLOW, Ansys CFD, COMSOL Multiphysics, OpenFOAM, FireFOAM, CFAST, and FDS, FireFOAM, CFAST, and FDS are specifically designed for fire simulation [32][33][34][35][36].These programs use numerical methods to solve Navier-Stokes equations that are suitable for low-speed and thermally driven flow.The CFD software focuses on smoke and heat transport and is designed to describe the evolution of fire.By developing CFD simulations, researchers can study the performance of multiscale dimensional materials or systems under different fire loads.
In this study, a PyroSim numerical model was developed for simulating the fire reaction behaviour of glass-polyurethane composite cladding.PyroSim is an FDS-based graphical user interface software that offers advanced features for simulating complex fire scenarios.There are three different numerical approaches commonly used in CFD simulation: Direct Numerical Simulation (DNS), Reynolds-Averaged Navier-Stokes equations (RANS), and LES [37].DNS resolves the governing equations of fluid flow, from the smallest Kolmogorov scale up to the integral scale, without using turbulent models.While DNS is useful for fundamental turbulence studies, its high computational cost limits its feasibility at larger scales.RANS, another commonly used approach, averages out the turbulence effects and solves for mean flow properties.While this approach is less computationally expensive than DNS or LES, it relies on the use of turbulence models to solve the equations, which can introduce additional uncertainty.Also, RANS models do not resolve the turbulent eddies responsible for temporal fluctuations, and they are not included in Pyrosim [38].On the other hand, LES captures important features of the turbulent flow field, such as eddy structures, without explicitly resolving all the turbulence scales.LES reduces the range of length scales that need to be computed by filtering out small scales and replacing them with sub-grid models.This makes LES computationally less expensive than DNS while still allowing it to provide sufficient accuracy, making it a more practical tool for studying turbulent flows for wide-ranging engineering problems [39][40][41][42][43].In summary, the choice of a specific numerical approach depends on the fire simulation needs and available computational resources.In this study, an LES-based model was used for the numerical simulation of the flammability of glass-polyurethane composite cladding.
A numerical analysis of the fire behaviour of glass-polyurethane composite cladding subjected to an incident radiant heat flux of 50 kW/m 2 under cone calorimeter conditions was simulated in PyroSim 2022, as illustrated in Figure 1a.The first step in setting up the numerical simulation model involved defining the domain size and grid size.In CFD simulations, both the domain and the grid sizes can impact the accuracy of the corresponding analysis.The first step involved the creation of an appropriate domain size representing the cone calorimeter specimen and thus measuring 100 × 100 × 9 mm.The X and Y dimensions were fixed at 160 mm, while the Z dimension was incrementally changed from 80 to 180 mm.The HRR-time profiles for each step change in the Z dimension are presented in Figure 1b.Convergence of the HRR-time profiles was observed for all simulated cases in which the Z dimension was at least 140 mm.As such, a domain size of 160 × 160 × 140 mm was adopted.After an appropriate domain size was established, the next step involved the determination of a suitable grid size.Grid sizes can be determined via the integral length (IL) method or according to the turbulent resolution (TR) concept.For fire simulations, the grid size is typically selected based on the ratio of the characteristic fire diameter, D * , to the grid cell size, δx, which is typically between 4 and 16 [21,44,45].The characteristic fire diameter, D * , is described in Equation (1): where (Q ) is the fire heat release rate (kW), ρ is the ambient temperature density of air (kg/m 3 ), c is the specific heat capacity of air (kJ/kgK), T is the ambient temperature (K), and g is the acceleration due to gravity (m/s 2 ).However, this approach is not suitable for scenarios in which the fire heat release rate is unknown.In such cases, the appropriate grid size can be determined by measuring the amount of resolved kinetic turbulent energy After an appropriate domain size was established, the next step involved the determination of a suitable grid size.Grid sizes can be determined via the integral length (IL) method or according to the turbulent resolution (TR) concept.For fire simulations, the grid size is typically selected based on the ratio of the characteristic fire diameter, D * , to the grid cell size, δx, which is typically between 4 and 16 [21,44,45].The characteristic fire diameter, D * , is described in Equation (1): Polymers 2023, 15, 3635 7 of 20 where ( .Q) is the fire heat release rate (kW), ρ ∞ is the ambient temperature density of air (kg/m 3 ), c p is the specific heat capacity of air (kJ/kgK), T ∞ is the ambient temperature (K), and g is the acceleration due to gravity (m/s 2 ).However, this approach is not suitable for scenarios in which the fire heat release rate is unknown.In such cases, the appropriate grid size can be determined by measuring the amount of resolved kinetic turbulent energy in the LES.The resolution quality is affected by the size of the fire and the grid size, which can be evaluated based on the resolved energy of the eddies for a specific grid size.
To determine the turbulent kinetic energy (TKE) in CFD simulations, one must carry out post-processing calculations involving directional velocities as in Equation ( 2): where ∼ u, ∼ v, and ∼ w are turbulent velocities in the X, Y, and Z directions, respectively.The turbulent resolution is calculated according to Equation (3): where k sgs is the sub-grid kinetic energy computed using Equation ( 4): where µ t is the turbulent viscosity, ρ is the density of air, C v is Deardorff's eddy viscosity constant, and ∆ is grid size.The accuracy of the simulation results is influenced by the resolution of the turbulent energy in the flow field, which is typically expressed as a percentage of the resolved turbulent kinetic energy with respect to the total energy.In LES, it is generally recommended to achieve a resolution of around 80% of the total kinetic energy for the accurate simulation of turbulent flows [46].Therefore, it was crucial to select appropriate simulation parameters and grid size to achieve sufficient turbulent resolution in our LES simulations.
A biased grid size approach may be implemented in CFD simulation to minimise simulation time.In this approach, the domain space is divided into sections representing fine and coarse domain spaces.It is imperative to maintain smaller grid sizes within the domain space closest to the fire initiation point (e.g., closest to the heat-exposed surface).The domain space between 0 and 60 mm along the Z dimension was modelled using a finer 4 mm grid size.The remaining domain space (e.g., 60 mm < Z < 140 mm) was modelled using a relatively coarse 8 mm grid size, as shown in Figure 2a.The HRR-time profile derived from the biased grid size approach matched the data derived from a fire simulation model in which a 4 mm grid size was used throughout the domain space, as shown in Figure 2b.The fire development stages at different times during a typical CFD analysis are shown in Figure 3.The developing fire was contained in the bottom domain space for exposure times of up to 265 s.Even after 315 s of thermal exposure, only a very small portion (<5% vol.%) of the developing fire encroached into the upper domain space.Moderateintensity fires generally achieve peak heat release rates at exposure times less than 315 s.As such, the biased grid size approach was implemented in all numerical simulations performed in this study due to its acceptable accuracy and reduced simulation time.
exposure times of up to 265 s.Even after 315 s of thermal exposure, only a very small portion (<5% vol.%) of the developing fire encroached into the upper domain space.Moderate-intensity fires generally achieve peak heat release rates at exposure times less than 315 s.As such, the biased grid size approach was implemented in all numerical simulations performed in this study due to its acceptable accuracy and reduced simulation time.Following the definition of the domain conditions, a 3D obstruction model was developed.This step was followed by the assignment of the fuel and reaction type as well as the physical and thermal properties of the test material.An incident radiant heat flux of 50 kW/m 2 was applied, followed by the setting of the simulation time, type, and step time.All domain faces were exposed to an air atmosphere except for the bottom domain.The back surface of the simulated cone calorimeter specimen was insulated.The radiant heat flux was continuously applied to the top surface of the specimen, simulating the cone calorimeter experimental test without a specimen holder.

Results and Discussion
The glass-polyurethane composite cladding density at ambient conditions was found to be ~1552 kg/m 3 .The TGA mass-temperature and corresponding derivatized thermogravimetric (dTG) profiles are shown in Figure 4.The reference temperature, T , used as an input in the numerical model is the temperature at the peak mass loss rate and was recorded as 320 °C (573 K).The pyrolysis temperature range, T , calculated as the difference between temperatures T and T (in Figure 4), was 250 °C (523 K).The heat of the reaction, h , was derived from the DSC heat flow data and corresponded to 240 kJ/kg.The thermal conductivity, k, of the glass-polyurethane composite cladding material at Following the definition of the domain conditions, a 3D obstruction model was developed.This step was followed by the assignment of the fuel and reaction type as well as the physical and thermal properties of the test material.An incident radiant heat flux of 50 kW/m 2 was applied, followed by the setting of the simulation time, type, and step time.All domain faces were exposed to an air atmosphere except for the bottom domain.The back surface of the simulated cone calorimeter specimen was insulated.The radiant heat flux was continuously applied to the top surface of the specimen, simulating the cone calorimeter experimental test without a specimen holder.

Results and Discussion
The glass-polyurethane composite cladding density at ambient conditions was found to be ~1552 kg/m 3 .The TGA mass-temperature and corresponding derivatized thermogravimetric (dTG) profiles are shown in Figure 4.The reference temperature, T ref , used as an input in the numerical model is the temperature at the peak mass loss rate and was recorded as 320 • C (573 K).The pyrolysis temperature range, T pyro , calculated as the difference between temperatures T 2 and T 1 (in Figure 4), was 250 • C (523 K).The heat of the reaction, h r , was derived from the DSC heat flow data and corresponded to 240 kJ/kg.The thermal conductivity, k, of the glass-polyurethane composite cladding material at ambient temperature was measured as 0.23 W/(mK), while the specific heat capacity, C p , was estimated to be 0.85 kJ/(kgK).A calculated emissivity (ε) value of 0.85 and an experimentally measured absorption coefficient, α, of 10 m −1 were used as input parameters.The heat of combustion, h c , measured via cone calorimetry was ~16,338 kJ/kg.The physical and thermal material properties used in the numerical fire simulation model are given in Table 1.The input parameters representing the lower (−12.5%)and upper (+12.5%)bounds relative to the nominal input values are also presented in Table 1.The lower and upper bound values were used for parametric studies, as discussed in Section 4.1.was estimated to be 0.85 kJ/(kgK).A calculated emissivity (ε) value of 0.85 and an experimentally measured absorption coefficient, α, of 10 m −1 were used as input parameters.The heat of combustion, h , measured via cone calorimetry was ~16,338 kJ/kg.The physical and thermal material properties used in the numerical fire simulation model are given in Table 1.The input parameters representing the lower (−12.5%)and upper (+12.5%)bounds relative to the nominal input values are also presented in Table 1.The lower and upper bound values were used for parametric studies, as discussed in Section 4.1.

Model Validation
The experimentally measured HRR-time profile for the glass-polyurethane composite cladding is shown in Figure 5 together with the numerically simulated data.When exposed to an incident radiant heat flux of 50 kW/m 2 under cone calorimeter conditions, autoignition occurred after 50 s of thermal exposure.The ignition of the glass-polyurethane composite cladding material generated additional thermal energy due to the combustion of volatiles released from the pyrolysis of the polyurethane binder.The combustion process led to a rapid increase in the heat release rate with a peak value (pHRR) of 45 kW/m 2 measured after 245 s of radiant heat exposure.Due to the depletion of the polyurethane component, the HRR rapidly declined to values just above 5 kW/m 2 a short period after the peak HRR event was achieved.The peak HRR and the time taken to reach this event are useful parameters for predicting the fire growth rate (FIGRA), which is the propensity of fire to increase in size or intensity.In this study, the fire growth rate was calculated by dividing the pHRR by the time taken to arrive at this event after ignition.FIGRA is an important parameter in fire safety engineering, as it can be used to assess the potential hazards of a fire and for the design of effective fire protection and suppression systems.For example, a high fire growth rate indicates that a fire may spread rapidly and, therefore, requires a more aggressive suppression strategy to control it.A FIGRA value of 0.23 kJ/m 2 s 2 was determined from the experimentally measured data.Another critical parameter in the assessment of cladding material fire hazards is the total heat release (THR).The total heat release was computed by integrating the area under the HRR-time graph in Figure 5, yielding the data presented in Figure 6 and a value of 11.6 MJ/m 2 .The total heat release represents the total energy contributed by the material to a fire and is essential in assessing the fire performance of building materials such as claddings.High THR values indicate that a material is likely to release a large amount of heat energy during combustion, which can contribute to the rapid spread of fire and increase the risk of structural damage or collapse.

Assessing the Sensitivity of the Numerical Model Results to Changes in Input Parameters
The Minitab ® software was used, employing the Design of Experiment (DoE) method, to perform statistical analysis of the pHRR data estimated from the numerical fire simulations of the glass-polyurethane composite cladding.The first step in the statistical analysis process involved conducting a screening analysis to identify and eliminate parameters with an insignificant effect on the pHRR.The initial screening was then followed by fractional factorial analysis, wherein parameters that had a significant influence on the pHRR values were identified.A 95% confidence level was used for the screening analysis, while

Assessing the Sensitivity of the Numerical Model Results to Changes in Input Parameters
The Minitab ® software was used, employing the Design of Experiment (DoE) method, to perform statistical analysis of the pHRR data estimated from the numerical fire simulations of the glass-polyurethane composite cladding.The first step in the statistical analysis process involved conducting a screening analysis to identify and eliminate parameters with an insignificant effect on the pHRR.The initial screening was then followed by fractional factorial analysis, wherein parameters that had a significant influence on the pHRR There was good agreement between the numerical fire simulation response and the experimentally measured data for the glass-polyurethane composite cladding, as shown in Figures 5 and 6.While the numerically simulated HRR-time profile has shifted to higher temperatures, the shape of the curve resembles that of the experimentally measured data.The numerical model revealed delays in the fire development and flame out stages, indicating a potential limitation of the model in simulating these phases of fire propagation.The numerically predicted pHRR value of 46 kW/m 2 matches the experimentally measured value of 45 kW/m 2 .The numerically derived FIGRA (0.21 kJ/m 2 s 2 ) and THR (10.2 MJ/m 2 ) values were within the experimental error range of the data derived from the physical experiment.The fair resemblance between the HRR-time profiles and the derived fire reaction properties indicates the high accuracy of the FDS-based fire simulation model developed in this study.Following the satisfactory validation of the numerical simulation model, we explored its sensitivity to variations in the input parameters using Minitab ® 2022 statistical analysis tools as discussed below.

Assessing the Sensitivity of the Numerical Model Results to Changes in Input Parameters
The Minitab ® software was used, employing the Design of Experiment (DoE) method, to perform statistical analysis of the pHRR data estimated from the numerical fire simulations of the glass-polyurethane composite cladding.The first step in the statistical analysis process involved conducting a screening analysis to identify and eliminate parameters with an insignificant effect on the pHRR.The initial screening was then followed by fractional factorial analysis, wherein parameters that had a significant influence on the pHRR values were identified.A 95% confidence level was used for the screening analysis, while factorial analysis was performed with a higher confidence level of 99%, accepting a lower error of 1%.To avoid undertaking the large number of simulations required for a full factorial analysis, a 1/8 fractional factorial approach was used, covering the maximum variation with a high confidence level.All input parameters provided to the numerical fire model were assumed to vary linearly between the lower and upper bounds and to have a significant impact on the pHRR.The input parameters and corresponding lower and upper bound levels used for parametric analysis are given in Table 1.
During the screening step, twenty-one simulations were performed, with the pHRR serving as the response variable.The design matrix, normalised parameters relative to the nominal input values, and the pHRR response for the twenty-one simulations are given in Table 2.The entries −1, 0, and 1 represent the lower bound, nominal, and upper bound values for the input parameters, respectively.An analysis of variance (ANOVA) test was performed on the design matrix, and the corresponding pHRR responses are given in Table 2. ANOVA data derived from the screening stage, which involved twenty-one simulations, are presented in Table 3.The data indicate that the specific heat capacity (p = 0.001), emissivity (p = 0.022), absorption coefficient (p = 0.000), heat of combustion (p = 0.000), reference temperature (p = 0.001), pyrolysis range (p = 0.000), and the heat of reaction (p = 0.000) all had a significant influence on the pHRR of the glass-polyurethane composite cladding material when the significance level was set at 5%.Density and thermal conductivity were found to have insignificant effects on the pHRR response.The screening analysis demonstrated a high R 2 value of 98% and an Adjusted R 2 value of 97%, indicating a significant amount of variation in the response function.This is an indication that the model is a good fit for the data and that the selected independent variables are important in predicting the dependent variable.The residual plot, which is shown with an even balance of pHRR response values on either side of the zero residual line in Figure 7, is a good indication of the reliability and validity of the statistical analysis conducted.Two input parameters that had a negligible impact on the pHRR response were eliminated before a detailed analysis was performed using the fractional factorial method.The fractional factorial analysis assessed the significance of the seven parameters on the pHRR of the glass-polyurethane composite cladding based on sixteen runs.Table 4 presents the design matrix and the pHRR response derived from the fractional factorial analysis.Two input parameters that had a negligible impact on the pHRR response were eliminated before a detailed analysis was performed using the fractional factorial method.The fractional factorial analysis assessed the significance of the seven parameters on the pHRR of the glass-polyurethane composite cladding based on sixteen runs.Table 4 presents the design matrix and the pHRR response derived from the fractional factorial analysis.Table 4. Design matrix and pHRR responses for sixteen fractional factorial analysis runs.

Run
To maintain brevity, ANOVA data derived from the fractional factorial assessment are not presented here.Instead, only those parameters found to have significant effects on the pHRR value are discussed.The heat of combustion (p = 0.000), absorption coefficient (p = 0.010), and heat of reaction (p = 0.004) were found to have significant effects on the pHRR of the glass-polyurethane composite cladding material at a 1% significance level.The R 2 (82%) and adjusted R 2 (77%) values were within an acceptable range, indicating a strong correlation between the independent variables and the dependent variable and the development of a highly effective model for explaining the variation in the response function.The residual plot of each factor on the pHRR is shown in Figure 8 based on the fractional factorial analyses, and, again, an even balance of pHRR response values on either side of the zero residual line is a good indication of the reliability and validity of the statistical analysis conducted.
Polymers 2023, 15, 3635 14 of 20 14 To maintain brevity, ANOVA data derived from the fractional factorial assessment are not presented here.Instead, only those parameters found to have significant effects on the pHRR value are discussed.The heat of combustion (p = 0.000), absorption coefficient (p = 0.010), and heat of reaction (p = 0.004) were found to have significant effects on the pHRR of the glass-polyurethane composite cladding material at a 1% significance level.The R 2 (82%) and adjusted R 2 (77%) values were within an acceptable range, indicating a strong correlation between the independent variables and the dependent variable and the development of a highly effective model for explaining the variation in the response function.The residual plot of each factor on the pHRR is shown in Figure 8 based on the fractional factorial analyses, and, again, an even balance of pHRR response values on either side of the zero residual line is a good indication of the reliability and validity of the statistical analysis conducted.It was crucial for us to identify the key parameters with significant effects on the pHRR in order to reduce the complexity of further analysis.In this case, the fractional factorial analysis led to the identification of three critical parameters that had a significant influence on the pHRR of the glass-polyurethane composite material: the heat of reaction, the heat of combustion, and the absorption coefficient.These parameters should be carefully considered during material selection to ensure the resulting cladding material satis- It was crucial for us to identify the key parameters with significant effects on the pHRR in order to reduce the complexity of further analysis.In this case, the fractional factorial analysis led to the identification of three critical parameters that had a significant influence on the pHRR of the glass-polyurethane composite material: the heat of reaction, the heat of combustion, and the absorption coefficient.These parameters should be carefully considered during material selection to ensure the resulting cladding material satisfies the desired level of fire performance.The effects of the three critical input parameters on the pHRR value were assessed by varying their magnitudes between the lower and upper bounds.The pHRR value was selected for this analysis since it provides valuable information on flammability and potential fire hazards in the event of a fire.The pHRR value provides an indication of flame intensity and the extent to which a developing fire can spread.As shown in Figure 9, when a linear variation in the input parameters was assumed, the pHRR monotonically increased with the increase in the heat of combustion and the absorption coefficient.In contrast, the increase in the heat of the reaction led to lower pHRR values.Since the heat of combustion is the amount of heat released when a substance undergoes complete combustion with oxygen, it is conceivable that the HRR values should increase with the increase in the heat of combustion.Materials with high absorption coefficients are more effective at absorbing radiant energy, which, in turn, accelerates pyrolysis and increases pHRR values.When the heat required to break down the polymer binder during pyrolysis is high, reduced amounts of combustion volatiles are produced, resulting in lower pHRR values.That is, when creating new glass-polymer composite cladding formulations, attention should be paid to the heat of combustion of polymer binders.Polymers with relatively high heat of combustion values should be avoided.On the other hand, thermally stable polymers (e.g., polymers with an elevated heat of reaction) are more suitable since the design objective for cladding material is to minimise the fire growth rate and flame intensity.Further, the amount of polymer and the resultant absorption coefficient of the glass-polymer composite are crucial parameters that should be carefully considered to meet stringent fire compliance requirements for cladding applications.Since the heat of combustion is the amount of heat released when a substance undergoes complete combustion with oxygen, it is conceivable that the HRR values should increase with the increase in the heat of combustion.Materials with high absorption coefficients are more effective at absorbing radiant energy, which, in turn, accelerates pyrolysis and increases pHRR values.When the heat required to break down the polymer binder during pyrolysis is high, reduced amounts of combustion volatiles are produced, resulting in lower pHRR values.That is, when creating new glass-polymer composite cladding formulations, attention should be paid to the heat of combustion of polymer binders.Polymers with relatively high heat of combustion values should be avoided.On the other hand, thermally stable polymers (e.g., polymers with an elevated heat of reaction) are more suitable since the design objective for cladding material is to minimise the fire growth rate and flame intensity.Further, the amount of polymer and the resultant absorption coefficient of the glass-polymer composite are crucial parameters that should be carefully considered to meet stringent fire compliance requirements for cladding applications.

Estimating Maximum Flame Intensity Using an Empirical Model
Numerical fire simulation models offer a cost-effective, accurate, and quicker alternative to performing physical experiments in relation to assessing the fire performance of new cladding materials.Further, numerical models can simulate complex and detailed physical processes and are easily adaptable to different scenarios and conditions.Despite a suite of advantages of numerical models over experimental data acquisition, there is room to further simplify predictive fire simulation models to reduce computational costs and time.Empirical models are generally simpler and easier to understand than numerical models, as they are based on direct relationships between input and output variables.Empirical models are often computationally quicker than their numerical counterparts.If the data used to develop the model are representative of the conditions of interest and the empirical model has been validated, this approach can provide quick, approximate predictions for a range of conditions.A linear empirical equation (Equation ( 5)) describing the relationship between the three significant parameters and the pHRR was developed using a multivariable data-fitting method based on the linear regression model.
The numerical fire model was used to simulate six different scenarios with three significant input parameters varied simultaneously.The same parameters were also used as an input for the empirical model.The pHRR values estimated using the numerical and empirical models were compared.The results showed strong agreement between the data estimated using the two models, with the data points closely following the x = y line, as shown in Figure 10.Although numerical models are more suitable for detailed fire scenario analysis and investigating underlying mechanisms compared to empirical models, this study demonstrates the potential of empirical models for the initial fire performance screening of new cladding materials.The empirical model presented herein can be used for the rapid screening of different cladding formulations such as with respect to changes in the type and weight percent of the polymer binder.This study validates the feasibility of using both numerical and empirical models to predict the maximum flame intensity (pHRR) in the event of glass-polyurethane composite cladding catching fire.the data used to develop the model are representative of the conditions of interest and the empirical model has been validated, this approach can provide quick, approximate predictions for a range of conditions.A linear empirical equation (Equation ( 5)) describing the relationship between the three significant parameters and the pHRR was developed using a multivariable data-fitting method based on the linear regression model.The numerical fire model was used to simulate six different scenarios with three significant input parameters varied simultaneously.The same parameters were also used as an input for the empirical model.The pHRR values estimated using the numerical and empirical models were compared.The results showed strong agreement between the data estimated using the two models, with the data points closely following the x = y line, as shown in Figure 10.Although numerical models are more suitable for detailed fire scenario analysis and investigating underlying mechanisms compared to empirical models, this study demonstrates the potential of empirical models for the initial fire performance screening of new cladding materials.The empirical model presented herein can be used for the rapid screening of different cladding formulations such as with respect to changes in the type and weight percent of the polymer binder.This study validates the feasibility of using both numerical and empirical models to predict the maximum flame intensity (pHRR) in the event of glass-polyurethane composite cladding catching fire.

Conclusions
In this study, the fire performance of a new cladding composite material incorporating recycled glass particles and polyurethane as a binder was assessed using a combination of experimental, numerical, and empirical methods.The findings of this study have several important implications and provide insights into the combustion behaviour of glass-polymer composite materials.The significance of this study lies in the use of numer-

Conclusions
In this study, the fire performance of a new cladding composite material incorporating recycled glass particles and polyurethane as a binder was assessed using a combination of experimental, numerical, and empirical methods.The findings of this study have several important implications and provide insights into the combustion behaviour of glasspolymer composite materials.The significance of this study lies in the use of numerical simulations as a cost-effective and environmentally friendly alternative to traditional experimental tests.The validated numerical model allowed us to conduct a parametric study, which identified the key material parameters that significantly influence the peak heat release rate (pHRR), a critical indicator of thermal performance.By incorporating these influential parameters, including the heat of combustion, absorption coefficient, and heat of reaction, an empirical equation was proposed to evaluate the combustion behaviour of glass-polyurethane composite materials in terms of pHRR.
The application of numerical simulations and the development of the empirical equation offer valuable benefits for researchers and practitioners in the field.They enable a more efficient assessment of material performance, saving time and resources while providing insights into the thermal behaviour of newly developed materials during a fire.This knowledge allows individuals to make informed decisions regarding the suitability of glass-polymer composite materials for specific applications, thereby enhancing fire safety in the building and construction industry.However, it is important to acknowledge the limitations of this study.The accuracy of the proposed empirical equation is contingent upon the range of parameters investigated.Therefore, future research should explore a wider range of glass and polymer proportions to enhance the precision and applicability of this equation across different compositions.
In conclusion, this study contributes to the understanding of the thermal behaviour of glass-polymer composite materials and their suitability for building claddings.The combination of numerical simulations, parametric studies, and the development of an empirical equation provides a valuable tool for evaluating combustion behaviour and improving fire safety in construction.Future research should focus on expanding the parameter range to further refine the proposed equation and advance the use of recycled glass-polymer composites in building materials.

Figure 1 .
Figure 1.Numerical simulation model for cone calorimeter conditions in the PyroSim software.(a) the adopted domain size and (b) calculated HRR-time profiles for different domain sizes.

Figure 1 .
Figure 1.Numerical simulation model for cone calorimeter conditions in the PyroSim software.(a) the adopted domain size and (b) calculated HRR-time profiles for different domain sizes.

Figure 2 .
Figure 2. (a) Biased grid size modelling approach and (b) HRR-time profiles derived from numerical fire simulations in which uniform and biased grid size approaches were considered.

Figure 2 .Figure 2 .
Figure 2. (a) Biased grid size modelling approach and (b) HRR-time profiles derived from numerical fire simulations in which uniform and biased grid size approaches were considered.

Figure 3 .
Figure 3. Fire development stages for a biased grid size numerical fire model simulating fire reaction behaviour under cone calorimeter conditions at (a) t = 265 s, (b) t = 315 s, and (c) t = 365 s.

Figure 3 .
Figure 3. Fire development stages for a biased grid size numerical fire model simulating fire reaction behaviour under cone calorimeter conditions at (a) t = 265 s, (b) t = 315 s, and (c) t = 365 s.

Figure 4 .
Figure 4. (a) Mass-temperature and (b) dTG profiles of the glass-polyurethane composite cladding.

Figure 4 .
Figure 4. (a) Mass-temperature and (b) dTG profiles of the glass-polyurethane composite cladding.

Figure 6 .
Figure 6.Experimentally and numerically simulated THR-time profiles for the glass-polyurethane composite cladding.

Figure 6 .
Figure 6.Experimentally and numerically simulated THR-time profiles for the glass-polyurethane composite cladding.

Figure 6 .
Figure 6.Experimentally and numerically simulated THR-time profiles for the glass-polyurethane composite cladding.

Figure 7 .
Figure 7. Residual plot for pHRR from screening analysis of the glass-polyurethane composite cladding.

Figure 7 .
Figure 7. Residual plot for pHRR from screening analysis of the glass-polyurethane composite cladding.

Polymers 2023, 15 , 3635 15 of 20 Figure 9 .
Figure 9. Variation in the pHRR response with changes in the heat of combustion, absorption coefficient, and heat of reaction of the glass-polyurethane composite cladding.

Figure 9 .
Figure 9. Variation in the pHRR response with changes in the heat of combustion, absorption coefficient, and heat of reaction of the glass-polyurethane composite cladding.

Figure 10 .
Figure 10.Comparison between the pHRR values predicted using the numerical fire model and those estimated from the empirical equation for the glass-polyurethane composite cladding.

Figure 10 .
Figure 10.Comparison between the pHRR values predicted using the numerical fire model and those estimated from the empirical equation for the glass-polyurethane composite cladding.

Table 1 .
Physical and thermal input parameters for the numerical fire simulation model.Lower-and upper-bound values are 12.5% lower and greater, respectively, than nominal input values.

Table 1 .
Physical and thermal input parameters for the numerical fire simulation model.Lowerand upper-bound values are 12.5% lower and greater, respectively, than nominal input values.

Table 3 .
ANOVA data from the screening analysis of the effect of different parameters on the pHRR response value.