Implementation of Numerical Model for Prediction of Temperature Distribution for Metallic-Coated Firefighter Protective Clothing

: The aim of this study is to predict the distribution of temperature at various positions on silver-coated firefighter protective clothing when subjected to external radiant heat flux. This will be helpful in the determination of thermal protective performance. Firefighter clothing consists of three layers, i.e., the outer shell, moisture barrier and thermal liner. The outer shell is the exposed surface, which was coated with silver particles through a physical vapor deposition process called magnetron sputtering. Afterwards, these uncoated and silver-coated samples were exposed to radiant heat transmission equipment at 10 kW/m 2 as per the ISO 6942 standard. Silver-coated samples displayed better thermal protective performance as the rate of temperature rise in silver-coated samples slowed. Later, a numerical approach was employed, contemplating the impact of metallic coating on the exterior shell. The finite difference method was utilized for solving partial differential equations and the implicit method was employed to discretize the partial differential equations. The numerical model displayed a good prediction of the distribution of temperature at different nodes with respect to time. The comparison of time vs. temperature graphs at different nodes for uncoated and silver-coated samples acquired from numerical solutions showed similar patterns, as witnessed in the experimental results.


Introduction
Clothing is one of the necessities of human life.Clothing creates a barrier between the human body and the surrounding environment and plays an integral role in maintaining the thermal balance of the human body [1].In hot climates, clothing assists in maintaining the thermoregulatory process of the wearer's body by creating a microclimate between the user's body and the clothing assembly.This is a very important factor for protective clothing when used in operational activities by rescuers, firefighters, and steel industry workers.Firefighters endure several challenges like extreme radiant heat flux densities, flames, chemical spillage, contact with hot objects and flash fires during their duties [2,3].Firefighter clothing consists of three layers, i.e., the outer shell, moisture barrier and thermal barrier.Each layer has its own unique functionalities [4].
The exterior shell consists of those materials which, when in contact with flame and heat, do not burn or degenerate, i.e., they inhibit ignition when they have contact with flames and must have characteristics of water repellence and decent thermal insulation.
Generally, fibers like combinations of Nomex and Kevlar (Nomex III A), Zylon, polybenzimidazole (PBI), kermel and some flame-retardant finishes like Pyrovatex and Proban are engaged to improve the thermal protective performance [5].The moisture barrier comprises a microporous membrane fixed between the outer shell and thermal barrier.This layer is responsible for permitting water vapor to pass through it.However, it is impermeable to liquid water.The fundamental objective of the moisture barrier is to secure the body of firefighters from liquefied chemicals.Moisture barriers are obtainable in the market as Action, Proline, Goretex, Cross tech and Neo Guard [5].The thermal barrier guards the bodies of firefighters by obstructing the environmental heat, and this barrier has flame-retardant fibers and their blends.It can be in the form of lining, knitted, non-woven, quilted batting, laminated woven and spun-laced fabric [5][6][7].For firefighter protective clothing, thermal protective performance (TPP) is a key issue.TPP means the ability of the firefighter clothing ensemble to protect the body of the firefighter against second-degree burns.Better thermal protective performance means more time can be spent by firefighters performing their duties without enduring severe burn injuries [8,9].There are several ways to improve the thermal protective performance of firefighter protective clothing.The thermal protective performance of firefighter protective clothing can be enhanced by increasing the thickness of firefighter protective clothing [10,11].However, if an increase in thickness results in a significant enhancement of the corresponding weight of the textile substrate, it might result in a decline in the thermal protective performance [12,13].With an increase in the number of air gaps between layers of protective clothing, the thermal protective performance of firefighter protective clothing is largely enhanced.However, the size limit and placement of air gaps in several layers are very critical for improvement in the thermal protective performance of firefighter protective clothing [14,15].By employing phase-change materials in firefighter protective clothing ensembles, the thermal protective performance can be improved.However, the influence of phase-change materials (PCMs) lasted a very short duration of time [15,16].Furthermore, by laminating the outer shell of firefighter clothing with aluminum foil, the thermal protective performance was increased.However, the breathability of firefighter protective clothing was compromised [16].Therefore, researchers are trying metallic particles for the impregnation of textile substrates, which might improve the thermal protective performance without significant differences in air permeability and water vapor permeability.Metallic particles of silver, aluminum oxide and titanium dioxide have elevated melting points and excellent reflective properties.A physical vapor deposition (PVD) methodology called magnetron sputtering can be used for the coating of these particles on textile substrates.In metallic particles, silver metal has a high melting point, i.e., approximately 962 • C. Furthermore, silver metal has outstanding reflective properties [17,18].Due to this reason, silver particles might be an appropriate prospect as a deposition layer on the exterior side of the outer shell of firefighter clothing.The thermal protective property of clothing materials against radiant heat flux can be measured according to ISO 6942/EN 366 [19,20].In the past two decades, many numerical models have been employed to study heat transfer in protective clothing.In 1997, Torvi established a heat transmission model for single-fabric ensembles.This model reported calculation methods for radiant heat transmission in textile substrates along with boundary conditions [21].Afterwards, several scientists tried to develop more accurate numerical models.A transient heat model was established by Kukuck and Prasad to study the impact of the rate of convective heat/rate of radiant heat flux on the thermal protective performance of textile substrates [22].Zhu, Zhang, and Song utilized a novel numerical model to investigate heat transmission in a cylinder sheathed by a flame-resistant textile substrate [23].Mercer and Sidhu reported the transfer of heat through protective clothing with deposited phase-change materials (PCMs) [24].A clothing numerical model proposed by Song et al. was utilized for the simulation of heat transmission in firefighter protective clothing [25].Ghazy and Bergstrom reported a modified heat transmission model for air space gaps in clothing ensembles based on a radiation heat transmission equation [26].Su, He and Li reported a heat storage model that was useful in the analysis of heat storage when subjected to heat and discharge of heat without pressure during cooling [27].Su, He and Li also developed a heat transmission model for protective clothing to determine the impact of compression and surface temperature on skin burns [28].In 2018, Su et al. established another numerical model for thermal protective fabric when exposed to heated steam.This model was based on the assumption of the effect of the flow jet stream between firefighter fabric and steam nozzles due to pressure differences and constant equilibrium between the three phases.The authors also employed a bio skin model and the Henriques burn integral in a steam model for the prediction of skin burns [29].
In 2012, Ackerman et al. reported that a horizontal bench-top tester could investigate the steam protective efficiency of fabric when the steam contacted the fabric with a range of pressure between 69 and 620 kPa.The results revealed that density, thickness, and air permeability played pivotal roles in protection against pressurized steam [30].Furthermore, a vertical system was employed to analyze these aforementioned factors by Yoo et al. and Havenith [31].In a study conducted by Desruelle and Schmid in 2004 [32], a vertical test system with controlled steam pressure and spraying distance was used to simulate different exposure situations.The purpose of that study was to examine the effect of warm steam on skin burns.Moreover, Sati et al. [33] established a simulation model for the steam protection of clothing ensembles under elevated steam pressure.The transmission of heat in one layer of firefighter protective clothing upon exposure to radiant heat and flame tests was determined by Zhu and Zhang [34].Yang et al. established a new approach to predict the thermal stress and skin damage of firefighter protective clothing when exposed to low-level radiation.The thermal flux generated under the fabric was studied by employing bench-top test equipment and using a thermoregulation simulation to model the temperature of skin in several parts of the human body.The acquired skin temperature data were employed to calculate skin burns and core body temperature for the detection of symptoms of thermal stress to determine the maximum exposure time [35,36].In 2018, Mandal et al. established separate models for the determination of thermal protective and thermophysiological comfort performances of fire-retardant fabrics [37].
Puszkarz and Machnowski used 3D geometry and morphology models of multilayer firefighter clothing ensembles by selected Computer-Aided Design (CAD) software 24.2.Finite volume methodology was employed for the simulation of heat transmission through multilayer ensembles when subjected to flames [38].Malaquais et al. studied the impact of dispersal of water on firefighter protective clothing in several environmental conditions, i.e., steady and unsteady states.Numerical methodology was utilized to model the transmission of heat and moisture in firefighter protective clothing.The clothing ensembles were exposed to heat fluxes ranging from 0 to 80 kW/m 2 .It was noted that the existence of water in the exterior shell enhanced the time for second-degree burns to occur.However, the presence of water in the inner shell had adverse impacts for high levels of heat flux [39].Kabir and Aghdam developed a Bézier-based multi-step method to investigate the nonlinear vibration and post-buckling structures of Euler-Bernoulli composite beams reinforced with graphene nano-platelets (GnPs) [40].Bert and Malik recommended the differential quadrature method as a substitute for numerical solution methodology for the initial and boundary value problems of applied sciences and engineering [41].
From the previous studies, it was noted that no research was performed on the prediction of the thermal protective performance of metallic-coated firefighter protective assemblies through numerical models.The present study deals with a numerical model for the prediction of temperature distributions at several positions along with the thickness of uncoated and silver-coated samples at different intervals of time.The numerical model implemented by Su et al. [27,42] was used with slight modifications for the estimation of the temperature distribution in uncoated and silver-coated firefighter protective clothing samples.In the end, a comparison was made between uncoated and silver-coated samples for the rate of the rise in temperature at different positions with respect to time intervals with the help of numerical solutions.The radiant heat transmission equations mentioned by Su et al. [27,42] and Torvi et al. [21] were employed to illustrate the radiant heat flux density transmitted towards the firefighter clothing assembly from a heat source.

Materials and Methods
The firefighter clothing samples used in this study consisted of an outer shell, moisture barrier and thermal barrier.These samples were provided by Vochoc Company, Plzen, Czech Republic.The specifications of these samples are given in Table 1 below.The outer shell of these samples was coated with a 1 µm coating of silver particles through magnetron sputtering.Two arrangements of fabric assemblies as given in Table 2 were made from these samples.A schematic diagram of the samples, along with their codes, is presented in Figure 1.ing samples.In the end, a comparison was made between uncoated and silver-coated samples for the rate of the rise in temperature at different positions with respect to time intervals with the help of numerical solutions.The radiant heat transmission equations mentioned by Su et al. [27,42] and Torvi et al. [21] were employed to illustrate the radiant heat flux density transmitted towards the firefighter clothing assembly from a heat source.

Materials and Methods
The firefighter clothing samples used in this study consisted of an outer shell, moisture barrier and thermal barrier.These samples were provided by Vochoc Company, Plzen, Czech Republic.The specifications of these samples are given in Table 1 below.The outer shell of these samples was coated with a 1 μm coating of silver particles through magnetron sputtering.Two arrangements of fabric assemblies as given in Table 2 were made from these samples.A schematic diagram of the samples, along with their codes, is presented in Figure 1.

Coating of Samples through Magnetron Sputtering
Magnetron sputtering is a physical vapor deposition process, where magnetically confined plasma is generated near the surface of the target metal (silver).Highly energetic ions with a positive charge collide with the surface of the negatively charged target metal.As a result, atoms are ejected from the target material, which then deposit on the surface of the substrate.Magnetron sputtering employs a strong magnetic field to confine electrons near the surface of the target material (metal).In consequence, the efficiency of Micro 2024, 4 372 the ionization process and deposition rate is enhanced, and plasma is created at a lower pressure [44,45].A schematic diagram of magnetron sputtering is shown in Figure 2. Coating of silver particles on the outer shell of firefighter clothing samples was conducted by the DC-reactive magnetron sputtering equipment (Prevac, Rogów, Poland).At first, the deposition chamber was evacuated at 4 Pascals.The distance between the sample holder and the target (silver metal) was maintained at 15 cm.Later, these samples were placed in a physical vapor deposition (PVD) chamber.The silver particle layer was acquired by sputtering of pure 99.99 percent silver (Ag) in the presence of 100 percent Argon gas floating at a speed of 10 sccm (standard cubic centimeter per minute).The thickness of the layer was 1 um, and it was almost uniform throughout the surface of the fabric.At the time of sputtering, the working pressure was maintained at a pressure of 1.3 Pascals.The power of the magnetron supply was maintained at 615 watts, and the negative substrate bias voltage was maintained at 300 volts.After sputtering, the color of the coated samples changed to silver on the surface of the fabric [43].
getic ions with a positive charge collide with the surface of the negatively charged target metal.As a result, atoms are ejected from the target material, which then deposit on the surface of the substrate.Magnetron sputtering employs a strong magnetic field to confine electrons near the surface of the target material (metal).In consequence, the efficiency of the ionization process and deposition rate is enhanced, and plasma is created at a lower pressure [44,45].A schematic diagram of magnetron sputtering is shown in Figure 2. Coating of silver particles on the outer shell of firefighter clothing samples was conducted by the DC-reactive magnetron sputtering equipment Prevace.At first, the deposition chamber was evacuated at 4 Pascals.The distance between the sample holder and the target (silver metal) was maintained at 15 cm.Later, these samples were placed in a physical vapor deposition (PVD) chamber.The silver particle layer was acquired by sputtering of pure 99.99 percent silver (Ag) in the presence of 100 percent Argon gas floating at a speed of 10 sccm (standard cubic centimeter per minute).The thickness of the layer was 1 um, and it was almost uniform throughout the surface of the fabric.At the time of sputtering, the working pressure was maintained at a pressure of 1.3 Pascals.The power of the magnetron supply was maintained at 615 watts, and the negative substrate bias voltage was maintained at 300 volts.After sputtering, the color of the coated samples changed to silver on the surface of the fabric [43].

Radiant Heat Transmission Equipment and Measurement Protocol
The X637 B radiant heat transmission machine uses the ISO 6942 standard to measure the transmission of heat through a material or material assembly.The temperature of the room was maintained between 15 °C and 35 °C.The radiant heat testing equipment was employed to investigate the radiant heat flux density through the material or material assembly according to the ISO 6942 standard.The apparatus consisted of six carbide rods serving as radiation heat sources, a small, curved copper plate calorimeter, and a movable test frame with a cooling device and sample holders as shown in Figure 3.At first, calibration was performed when the moveable screen was withdrawn and reverted to a point when a rise in temperature of 30 °C was reached, and the incident heat flux density Qₒ was measured.Afterwards, the sample was affixed to one side of the plate of the sample holder and held in contact with the face of the calorimeter, applying a mass of 200 g.The movable screen was withdrawn, and the starting point of the radiation was recorded.The movable screen was returned to its closed position after a temperature rise of about 30 °C was reached.The size of the sample was 230 mm × 80 mm.The time for

Radiant Heat Transmission Equipment and Measurement Protocol
The X637 B radiant heat transmission machine uses the ISO 6942 standard to measure the transmission of heat through a material or material assembly.The temperature of the room was maintained between 15 • C and 35 • C. The radiant heat testing equipment was employed to investigate the radiant heat flux density through the material or material assembly according to the ISO 6942 standard.The apparatus consisted of six carbide rods serving as radiation heat sources, a small, curved copper plate calorimeter, and a movable test frame with a cooling device and sample holders as shown in Figure 3.At first, calibration was performed when the moveable screen was withdrawn and reverted to a point when a rise in temperature of 30 • C was reached, and the incident heat flux density Q o was measured.Afterwards, the sample was affixed to one side of the plate of the sample holder and held in contact with the face of the calorimeter, applying a mass of 200 g.The movable screen was withdrawn, and the starting point of the radiation was recorded.The movable screen was returned to its closed position after a temperature rise of about 30 • C was reached.The size of the sample was 230 mm × 80 mm.The time for temperature accelerations of 12 • C and 24 • C in the calorimeter was determined, and conclusions are mentioned in the form of the radiant heat transmission index (RHTI 12 and RHTI 24), the percentage heat transmission factor (percentage TF Q o ) and the transmitted heat flux density Q c [48].Before experimentation, all samples were pre-conditioned for 24 h at a temperature of 20 • C and had a relative humidity of 65 percent [48].Three samples were required for testing at each level of heat flux density.Incident heat flux density is evaluated using Equation (1).
The percentage age heat transmission factor [Percentage TF Qo] for the incident heat flux density level is explained by Equation ( 4) [48]: Figure 3. Radiation heat testing equipment [43].
These samples were exposed to a 10 kW/m 2 heat flux.To determine the temperature distribution between multilayer clothing assemblies, two K-type thermocouples were placed between the outer shell and moisture barrier and between the moisture barrier  The percentage age heat transmission factor [Percentage TF Q o ] for the incident heat flux density level is explained by Equation (4) [48]: These samples were exposed to a 10 kW/m 2 heat flux.To determine the temperature distribution between multilayer clothing assemblies, two K-type thermocouples were placed between the outer shell and moisture barrier and between the moisture barrier and thermal barrier as shown in Figure 4. Thermocouple 1 measures the temperature of the outer shell and thermocouple 2 measures the temperature of the moisture barrier.
To simplify the explanation with theoretical equations, the following assumptions have been made: i.
Transmission of heat takes place in one dimension only.ii.
The transfer of mass is negligible.iii.
Radiation only penetrates through the exterior shell of the multilayer assembly as almost 95% of incident energy is in the form of radiation that is absorbed after covering a distance equivalent to the outer shell thickness.iv.
Optical characteristics like transmissivity, reflectivity and absorptivity are assumed to be constant [21,27,42,43].The temperature of uncoated and silver coated samples against time after exposure to 10 kW/m 2 are shown in Figures 5 and 6, respectively.To simplify the explanation with theoretical equations, the following assumptions have been made: i.
Transmission of heat takes place in one dimension only. ii.
The transfer of mass is negligible. iii.
Radiation only penetrates through the exterior shell of the multilayer assembly as almost 95% of incident energy is in the form of radiation that is absorbed after covering a distance equivalent to the outer shell thickness.
The temperature of uncoated and silver coated samples against time after exposure to 10 kW/m 2 are shown in Figure 5 and Figure 6 respectively.To simplify the explanation with theoretical equations, the following assumptions have been made: i.
Transmission of heat takes place in one dimension only. ii.
The transfer of mass is negligible. iii.
Radiation only penetrates through the exterior shell of the multilayer assembly as almost 95% of incident energy is in the form of radiation that is absorbed after covering a distance equivalent to the outer shell thickness.
The temperature of uncoated and silver coated samples against time after exposure to 10 kW/m 2 are shown in Figure 5 and Figure 6 respectively.A quantitative assessment of the rate of the rise in temperature in uncoated samples and silver-coated samples is shown in Table 3.It was noted from Table 3 that the RHTI 24 [sec] value of silver-coated samples was greater compared to uncoated samples A. This indicated that the time for a rise of 24 °C was greater in silver-coated samples as  A quantitative assessment of the rate of the rise in temperature in uncoated samples and silver-coated samples is shown in Table 3.It was noted from Table 3 that the RHTI 24 [sec] value of silver-coated samples was greater compared to uncoated samples A. This indicated that the time for a rise of 24 • C was greater in silver-coated samples as compared to uncoated samples.The transmitted heat flux density Q c [kW/m 2 ] was lower in the case of silver-coated samples as compared to uncoated samples.The lower value of transmitted heat flux density indicated that a lower amount of radiant heat flux density was passed through silver-coated samples for a given value of time.Similarly, the lower value of the percentage heat transmission factor for silver-coated samples indicated that a lower amount of heat was transmitted in silver-coated samples.This might be due to the good reflective properties of silver particles, which dissipated the external radiant heat flux.The reflective property of silver-coated samples was measured with a Mid-Infrared Integrator with an FTIR spectrometer as mentioned in Table 3. Figure 7 shows the time vs. temperature values of uncoated and silver-coated samples when subjected to 10 kW/m 2 .It was noted that a flatter curve was obtained for silvercoated samples A1.On the other hand, a relatively steep curve was obtained for uncoated samples A. This indicated that the rate of the rise in temperature of silver-coated samples A1 took place at a slower rate as compared to uncoated samples A [43].

Transmission of Heat from Heating Source to Firefighter Clothing Assembly
As described earlier, uncoated and silver-coated fabric assemblies were subjected to a heat flux of 10 kW/m 2 with the help of a heating source.This process of heat transfer in firefighter clothing assemblies due to heat fluxes induced by heating rods can be represented with the help of an energy equation as mentioned by Su et al. [27,42]: where ρ fb , Cp fb and λ fb are the density, specific heat capacity and thermal conductivity of firefighter fabric, respectively.∂q rad−abs is the absorbed portion of radiant heat flux from the heating source to the outer shell.Equation ( 6) can also be written as follows:

Transmission of Heat from Heating Source to Firefighter Clothing Assembly
As described earlier, uncoated and silver-coated fabric assemblies were subjected to a heat flux of 10 kW/m 2 with the help of a heating source.This process of heat transfer in firefighter clothing assemblies due to heat fluxes induced by heating rods can be represented with the help of an energy equation as mentioned by Su et al. [27,42]: ∂T ∂t where ρ fb , Cp fb and λ fb are the density, specific heat capacity and thermal conductivity of firefighter fabric, respectively.∂q rad−abs is the absorbed portion of radiant heat flux from the heating source to the outer shell.Equation ( 6) can also be written as follows: where T(x, t) is the temperature field under investigation and t is the total time during which we investigate the temperature field.This equation is a second-order transient parabolic differential equation.It was assumed that textile material was homogeneous and isotropic.Therefore, the values of λ fb , C pfb and ρ fb are assumed to be constant.∇ designates the Vector Hamilton operator in a Cartesian coordinate system.Q rad−abs (x, t) is the part of radiation absorbed by the outer shell from the radiant heat source.It was assumed that the textile material was homogeneous and isotropic.Therefore, the values of λ fb , C pfb and ρ fb are assumed to be constant [49][50][51][52].
It is also possible to rewrite Equation ( 7) in the following way: where the symbol ∆ designates the Laplace operator ∆T = ∂ 2 T(x,t) The thermal energy from the heating source to the surface of the cloth is transported due to thermal radiation and convection of air.In this scenario, the mode of radiation is dominant as compared to convection.The radiant heat transmission between the heating source and the outer shell is dependent on the differentiation of temperature and radiation view factor, as shown in Figure 8 [27,42,43].

Numerical Solution
To solve Equation ( 6), the finite difference method was utilized.The implicit method was employed for the discretization of the second-order parabolic differential equation.The values of temperature at discrete points (xi and tj) are indicated by T i j .
Due to the scattering and absorption of radiation, Beer's law was utilized to explain the extinction of the incident thermal radiation ( Q rad ) in multilayer fabric systems [42,53].Beer's law is mostly used for liquid solutions; however, it can also be utilized for porous mediums like textiles [43].

Numerical Solution
To solve Equation ( 6), the finite difference method was utilized.The implicit method was employed for the discretization of the second-order parabolic differential equation.The values of temperature at discrete points (x i and t j ) are indicated by T j i .Due to the scattering and absorption of radiation, Beer's law was utilized to explain the extinction of the incident thermal radiation (Q rad ) in multilayer fabric systems [42,53].
Beer's law is mostly used for liquid solutions; however, it can also be utilized for porous mediums like textiles [43].
∂Q rad−abs (x, t) ∂x = [q rad−rhs − q rad−surr ](e −γ fb i∆x )γ fb ( 9) σ is the Stefan-Boltzmann constant (5.678 × 10 −8 W/m 2 K 4 ).T rhs is the temperature of the radiant heat source in K; T surr is the temperature of the surrounding atmosphere in K; A fb is the area of firefighter fabric in m 2 ; A rhs is the area of the radiant heat source in m 2 ; ε rhs is the emissivity of the radiant heat source; ε o.shell is the emissivity of the outer shell; V rhs-o.shell is the view factor from the radiant heat source towards the outer shell [27,42]; V o.shell-surr is the view factor from the outer shell towards the surrounding atmosphere.
The view factor is a fraction of the radiation leaving one surface which is received by another surface.It was assumed that the heating source and exterior shell were parallel coaxial disks of unequal radii [42,54].
L is the distance between the heating rods and the multilayer clothing assembly in meters, r o.shell is the radius of the outer shell and r rhs is the radius of the heating source.
γ fb is the extinction coefficient of the outer shell, which is illustrated by the following equation: where τ is the transmissivity of the outer shell and h shell is the thickness of the outer shell [49].Some important values are mentioned in Table 4 [43].Equation ( 7) can also be written as Partial derivatives in Equation (15) were substituted by an implicit finite divided scheme, generating a discrete description related to position x i and time t j : For further solutions, the Thomas algorithm method was employed to solve the abovementioned equation.The Thomas algorithm is a tridiagonal matrix algorithm, utilized for the solution of second-order partial differential equations.It provides an effective way to resolve tridiagonal systems of equations during the discretization of second-order partial differential equations.The Gaussian elimination method is streamlined by this algorithm [56] as follows: Thickness of fabric Number of nodes The fabric assembly is divided into 5 nodes as seen in Figure 9.
The fabric assembly is divided into 5 nodes as seen in Figure 9.For solving partial differential equations, a node means a point in the discretized domain where the computation of solutions took place.For implicit methodology, the solution at a particular node is dependent on the solution at the adjacent neighboring nodes, and equations associated with solutions must be resolved in a simultaneous manner.Due to this reason, a system of equations was generated that must be resolved in an iterative manner, and the solution at every node was updated till convergence was achieved [43,57].For solving partial differential equations, a node means a point in the discretized domain where the computation of solutions took place.For implicit methodology, the solution at a particular node is dependent on the solution at the adjacent neighboring nodes, and equations associated with solutions must be resolved in a simultaneous manner.Due to this reason, a system of equations was generated that must be resolved in an iterative manner, and the solution at every node was updated till convergence was achieved [43,57].

Boundary Conditions
The initial temperature was 301.45 K and its schematic diagram is shown in Figures 10 and 11, respectively.λ o.shell .
∂T(x, t) ∂x = −α(T o.shell − T surr ) + q rad−abs (17) where α represents the coefficient of heat transfer, T surr is the temperature of the surrounding atmosphere and q rad−abs is the part of radiation absorbed by the outer shell [27,42].
ε T.B is the emissivity of the thermal barrier and ε cal is the emissivity of the calorimeter.T T.B is the temperature of the thermal barrier in Kelvin and T cal is the temperature of the calorimeter in Kelvin [27].
The distribution of temperature at different nodes in different intervals of time for uncoated and silver-coated samples is depicted in Figures 10 and 11  where  represents the coefficient of heat transfer, T surr is the temperature of the s rounding atmosphere and q rad−abs is the part of radiation absorbed by the outer sh [27,42].

Results and Discussion
Distributions of temperatures for uncoated and silver-coated samples were determined at different intervals of time with the help of numerical solutions, and later, these results were compared with each other.
It can be witnessed from Figures 12 and 13 that there was a sequential increment in the rise in temperature for different curves with an increase in time.However, this increment was not in a linear form.It was also noticed that as the number of nodes increased, there was a decline in the values of temperature for the same time period.This indicates a decline in temperature values with an increase in the thickness of firefighter clothing assemblies for both uncoated samples A and silver-coated samples A1.
It can be noticed from Figure 14 that in the first sixteen seconds, the curves of uncoated samples A and silver-coated samples A1 overlap with each other.Afterwards, a gap appears between the curves of silver-coated and uncoated samples.This gap continued increasing with the increase in time.
It can be witnessed from Figures 15-18 that in the first 15 s, both curves of uncoated and silver-coated samples superimposed each other.Afterwards, the curve of uncoated samples started to get steeper as compared to the curves of silver-coated samples.It was also observed that as the number of nodes increased, the curve of both uncoated and silver-coated samples became flatter as compared to previous nodes.The slackness in the curve indicates that the rate of the rise in temperature takes place slowly.Consequently, the thermal protective performance will be better because less heat is transmitted towards the fabric assembly.
results were compared with each other.
It can be witnessed from Figures 12 and 13 that there was a sequential increment in the rise in temperature for different curves with an increase in time.However, this increment was not in a linear form.It was also noticed that as the number of nodes increased, there was a decline in the values of temperature for the same time period.This indicates a decline in temperature values with an increase in the thickness of firefighter clothing assemblies for both uncoated samples A and silver-coated samples A1.It can be noticed from Figure 14 that in the first sixteen seconds, the curves of uncoated samples A and silver-coated samples A1 overlap with each other.Afterwards, a gap appears between the curves of silver-coated and uncoated samples.This gap continued increasing with the increase in time.It can be noticed from Figure 14 that in the first sixteen seconds, the curves of uncoated samples A and silver-coated samples A1 overlap with each other.Afterwards, a gap appears between the curves of silver-coated and uncoated samples.This gap continued increasing with the increase in time.ver-coated samples became flatter as compared to previous nodes.The slackness in the curve indicates that the rate of the rise in temperature takes place slowly.Consequently the thermal protective performance will be better because less heat is transmitted towards the fabric assembly.ver-coated samples became flatter as compared to previous nodes.The slackness in the curve indicates that the rate of the rise in temperature takes place slowly.Consequently, the thermal protective performance will be better because less heat is transmitted towards the fabric assembly.It is evident from Figure 19 that the pattern of the curves for the uncoated and silver-coated samples is almost the same after the initial 5 s.Later, the curve of uncoated samples starts to rise by creating a gap between the curves of uncoated A and silver coated samples A1.The wideness of this gap between the two curves increases till 100 s This indicates the better thermal protective performance of silver-coated samples A1.It is evident from Figure 19 that the pattern of the curves for the uncoated and silvercoated samples is almost the same after the initial 5 s.Later, the curve of uncoated samples starts to rise by creating a gap between the curves of uncoated A and silver-coated samples A1.The wideness of this gap between the two curves increases till 100 s.This indicates the better thermal protective performance of silver-coated samples A1.It can be inferred that the numerical model solution predicts the temperature distribution at different nodes for uncoated samples A and silver-coated samples A1 at different intervals of time.This model also shows the silver-coated samples A1 incur a less steep curve of temperature values as compared to uncoated samples A at different nodes.This flatness of the curve was also witnessed in the case of silver-coated samples A1 for experimental work.With the help of the numerical model, it was possible to estimate the temperature distribution at different places along with the thickness of the fabric, which was not possible in experimental work.

Conclusions
Firefighters are always putting their lives in danger when performing their duties.For firefighter protective clothing, thermal protective performance is the most important property.It is the capability of firefighter clothing to protect firefighters from seconddegree burns.The thermal protective performance can be improved by metallic silver It can be inferred that the numerical model solution predicts the temperature distribution at different nodes for uncoated samples A and silver-coated samples A1 at different intervals of time.This model also shows the silver-coated samples A1 incur a less steep curve of temperature values as compared to uncoated samples A at different nodes.This flatness of the curve was also witnessed in the case of silver-coated samples A1 for experimental work.With the help of the numerical model, it was possible to estimate the temperature distribution at different places along with the thickness of the fabric, which was not possible in experimental work.

Conclusions
Firefighters are always putting their lives in danger when performing their duties.For firefighter protective clothing, thermal protective performance is the most important property.It is the capability of firefighter clothing to protect firefighters from second-degree burns.The thermal protective performance can be improved by metallic silver coatings through magnetron sputtering.It was noted from a radiant heat transmission machine that the rate of the rise in temperature was low in silver-coated samples because of excellent reflective properties.A numerical model and its solutions were established using appropriate equations for the transmission of heat.The distribution of temperature through the multilayer sandwich structure was determined with the help of the numerical model and its solutions.Previously, no attempt was made to predict the distribution of temperature for firefighter clothing coated with low-emissivity silver particles.The structure of the sandwich assembly was divided into five nodes and finite difference methodology was selected to solve partial differential equations.Discretization of the partial differential was carried out by implicit methodology.The Thomas algorithm method was used for further solutions of the equations.Thus, the temperature distribution of multilayer firefighter ensembles was studied through these techniques.The outcomes of this study might be beneficial for researchers and industrial partners to predict the thermal protection of firefighter protective clothing with lower-emissivity coatings.

Figure 1 .
Figure 1.Schematic diagram of arrangement of fabric assemblies along with their codes.

Figure 1 .
Figure 1.Schematic diagram of arrangement of fabric assemblies along with their codes.

Figure 3 .
Figure 3. Radiation heat testing equipment [43].A = the area of the copper plate in m 2 , a = the absorption coefficient of the painted surface of the calorimeter, M = the mass of the copper plate in kg, C p = the specific heat of copper (0.385 kJ/kg • C), R = the rate of the rise in the calorimeter temperature in the linear region in • C/s.The transmitted flux density Q c , in kW/m 2 , is evaluated by Equation (2):

K=
Mean rate of escalation of the calorimeter temperature in • C/s in the region between a 12 • C and 24 • C rise.RHTI 12 = Threshold time in sec when the temperature of the calorimeter increases by 12 • C. RHTI 24 = Threshold time in sec when the temperature of the calorimeter increases by 24 • C.

Figure 5 .
Figure 5. Temperature of uncoated samples with respect to time when samples are exposed to 10 kW/m 2 [43].

Figure 5 .
Figure 5. Temperature of uncoated samples with respect to time when samples are exposed to 10 kW/m 2 [43].

Figure 6 .
Figure 6.Temperature of silver-coated samples with respect to time when samples are exposed to 10 kW/m 2 [43].

Figure 7 .
Figure 7. Temperature of firefighter samples with respect to time when samples are exposed to 10 kW/m 2 [43].

Figure 8 .
Figure 8. Schematic diagram showing equations involved in transmission of heat from heat source towards outer shell [43].

Figure 8 .
Figure 8. Schematic diagram showing equations involved in transmission of heat from heat source towards outer shell [43].

Figure 9 .
Figure 9. Schematic diagram of uncoated and silver-coated samples with representation of nodes.Figure 9. Schematic diagram of uncoated and silver-coated samples with representation of nodes.

Figure 9 .
Figure 9. Schematic diagram of uncoated and silver-coated samples with representation of nodes.Figure 9. Schematic diagram of uncoated and silver-coated samples with representation of nodes.
Number of time steps = Time f inal −Time initial ∆t = 100 sec−0 sec 1 = 100.Interval of time ∆t = 1 sec.T j i corresponds to the temperature at node i. x = i.∆x and t = j.∆t.

Figure 11 .
Figure 11.Schematic diagram showing equation for boundary condition of thermal barrier [43].

Figure 12 .
Figure 12.Temperature distribution in uncoated samples with respect to time [43].

Figure 14 .Figure 13 . 15 Figure 13 .
Figure 14.Comparison of temperature distribution of uncoated and silver-coated samples at node [0] [43].It can be witnessed from Figures 15-18 that in the first 15 s, both curves of uncoated

Figure 14 .Figure 14 .
Figure 14.Comparison of temperature distribution of uncoated and silver-coated samples at node [0] [43].It can be witnessed from Figures 15-18 that in the first 15 s, both curves of uncoated and silver-coated samples superimposed each other.Afterwards, the curve of uncoated samples started to get steeper as compared to the curves of silver-coated samples.It was also observed that as the number of nodes increased, the curve of both uncoated and sil-

Figure 17 .
Comparison of temperature distribution of uncoated and silver-coated samples at node[3] [43].

Table 3 .
Calorimetric properties of samples.