Analytical Solution for Adiabatic Surface Temperature (AST)

In this contribution an analytical solution of the equation of an ideal surface of a perfect insulator is introduced. The new solution bases on the solution of the heat balance equation, which is fourth degree polynomial. The resulting expression for adiabatic surface temperature (AST) is discussed and verified. This solution can be easily incorporated, with negligible computational cost, inside the computational fluid dynamics solvers for fire simulations. Finally, since AST is close to the temperature measured by plate thermometers, AST is easy to measure in any arbitrary point in a computational domain. This makes it possible to validate numerical results with the temperatures measured by plate thermometers in experiments. Furthermore, an expression, which is given in the end of the article, describes the physical quantity with a closed-form solution. This may be potentially used for analyses performed for better understanding of physical process, e.g. during the experiments involving localized fires.


Introduction
Since Professor Ulf Wickstro¨m introduced the concept of adiabatic surface temperature (AST) to the fire science community in 2007 [1], it appears to be a very efficient way for expressing heat exposure of the solid surfaces both in real experiments, using plate thermometers [2][3][4] which measure the temperature that is close to the AST [5], and in numerical analyses [6][7][8][9]. Adiabatic surface temperature can be also used as a single thermal boundary conditions when calculating temperature of structures exposed to fire [10,11], which is one of the biggest advantages of this concept. Thus, such computational fluid dynamics (CFD) numerical codes as fire dynamics simulator incorporated AST as the variable that can be obtained at solid surfaces [12]. Nonetheless, simulating fires is computationally very expensive and the size of numerical grid usually does not allow to accurately model structural elements. In case of steel sections they usually even do not exist in CFD models, and then it is not possible to explicitly describe the thermal conditions at their section surfaces. In this case, many researches extract necessary quantities from CFD analyses, such as: incident radiations, gas velocities, gas temperatures, which finally ends up with a big amount of data that have to be posteriorly processed [13]. This problem has been already discussed by Joakim Sandstro¨m in his thesis [14]. He developed a method for numerical evaluation of AST at arbitrary points and directions in the computational domain, not necessarily connected to any actual surface. That method requires additional computation time, in order to numerically calculate the AST by solving the heat balance equation. In this paper, an analytical solution for the same problem is introduced, scientifically discussed and verified. This analytical approach, resulting in closed-form solution for AST, is advantageous not only from the numerical point of view, but may be also beneficial for future scientists examining the influence of natural fires on the dynamics of heat distribution in space.

Heat Balance Equation
There are two contributions of the total (net) heat flux _ q 00 net to a surface: convective heat flux _ q 00 conv and radiative heat flux _ q 00 rad . The sum _ q 00 net ¼ _ q 00 rad þ _ q 00 conv can be obtained by the following formula: where e is the emissivity of the surface, _ q 00 inc the incident radiation, r the Stefan-Boltzmann constant, h c the convective heat transfer coefficient. Finally T g is the gas temperature and T s is the surface temperature.
The total net heat flux to an ideal surface of a perfect insulator is by definition zero, thus introducing concept of AST T AST , following expression has to be fulfilled (according to Wickstro¨m [1]): After elementary algebraical operations, Eq. (2) can be written as: what is the fourth order polynomial equation with T AST as variable, that may be written in the form: where subsequent coefficients of the polynomial are: Since a is always positive for gray and black bodies (e > 0), Eq. (4) has always four roots, but some of the roots may be complex numbers.

Solution
In order to solve an Eq. (4), the WolframAlpha on-line system for symbolic and numeric mathematics has been used [15]. Finally the output from the computational engine is simplified and the solutions of Eq. (4) can be written in the following forms: Using coefficients a, b and c as previously introduced, a, b and c are obtained using following formulas: Analytical Solution for Adiabatic Surface Temperature 415

Discussion
From above consideration, it can be seen that there are four equivalent solutions resulting with four potential values for adiabatic surface temperature. Those four solutions, without the scientific deliberation, are only unprofitable mathematical formulas. In order to give the physical sense, it is necessary to come back to the physical meaning of the coefficients within those formulas. Since a ¼ er and, from the definition, ðe; rÞ > 0, the coefficient a is always positive. Similarly, since b ¼ h c , to have physical meaning, it must be positive (for h c ¼ 0 we have trivial solution explained in next section). On the other hand c ¼ ÀðerT 4 r þ h c T g Þ is valid for temperatures in Kelvin and always gives negative value.
Considering now relations (10)(11)(12) it is clearly seen, that for positive values of a and b, and negative value of c, the coefficients a and c are clearly positive, and coefficient b is negative. Let us now look at the coefficient M. If M takes the real value, then the expression under square root must be non-negative. At the same time, expressions (5-8) may be evaluated only if M is positive. Thus, after standard algebraic operations M can be rearranged to: Since ac > 0, the value of M > 0 () bc þ a 2 > 0. Substituting expressions (10)(11)(12) into that right-hand side inequality we obtain: Let us now conduct basic mathematical operations on the inequality (14) as following: Finally, for physical values of a, b and c, inequality (15) is always true, so that coefficient M is always positive. Coming back now to expressions (5)(6)(7)(8) it is clearly visible that: expressions (5) and (6) have no real evaluation; expression (7) gives negative (unphysical) value; the only expression that may give real and physical result is expression (8). Thus, let us now recall expression (8) as the exact solution for evaluation of adiabatic surface temperature: Eq. (16) is finally proven by tests over a range of e=h c ratios. In Figure 1, it is shown, that the dependence between AST and parameters of Eq. (2) are the shape of those referred by Sandstro¨m [14]. Additionally, quantitative comparison between AST obtained with solution (16) and approximate solution obtained using Newton-Raphson method is considered in next section.

Exact Versus Approximate Solution
The biggest advantage of using a closed-form exact solution is no need for performing time consuming, iterative computations, in order to obtain the result. Currently, AST is calculated mostly using the Newton-Raphson method. According to author's tests, AST can be approximated within no more than three iterations, with the error less than 0:5 C. This is valid for different configurations of gas and incident heat flux temperatures, when the incident heat flux temperature is taken as the first guess and the derivative of the function is calculated directly (not approximated). That means, using a closed-form solution can speed-up the process of computations about three times; the solution is not sensitive to ''the first guess'' and the approximation of derivative; finally solution is pure physical, not affected by approximation error. Comparison between exact and approximate solution is shown in Figure 2 (graphs overlap). Analytical Solution for Adiabatic Surface Temperature 417

Fire Safety Engineering Application
The usefulness of AST concept in fire safety engineering (FSE) is related to the heat transfer calculation. During the experiments, the crucial issue is the determination of heat boundary conditions on the specimen surface. Wickstro¨m, in [16], summarizes how important the determination of heat flux is and how many problems are related with measurements of the physical quantities in experiments. Thus, the plate thermometers are used as the devices that can measure the exposure of a surface both to convection and radiation. A good and simple example of an application of AST concept is given by Bystro¨m et al. [17]. Authors compare a temperature data measured using plate thermometer and two types of thermocouples in a test carried out in cone calorimeter with burning specimen under the cone shape radiation panel. Bystro¨m summarizes, that plate thermometer is sensitive both to convection and radiation in a similar way as real specimen and the AST can be measured using plate thermometers even under harsh fire conditions. In the same way, one can imagine the test, where the AST is to be specified at a specimen surface with the conditions given by incident radiation and gas temperature. Then AST may be calculated directly from Eq. 16 and can be used as a single quantity for heat transfer calculations.

Final Remarks
An analytical solution of the polynomial expressing the heat balance between the total (net) heat flux approaching the ideal surface of perfect insulator is introduced. The resulting expressions for roots of heat balance equation are discussed and the only expression having physical meaning is chosen as the one, that describes adiabatic surface temperature. Results obtained by proposed solution are quantitatively and qualitatively checked with respect to previous researchers Comparison between exact and approximate solutions for T g ¼ 20 C and T g ¼ 500 C (graphs overlap).

418
Fire Technology 2017 findings. It is shown, that since the approximate iterative solution is sensitive to ''the first guess'' and the approximation of derivative, proposed analytical solution gives exact results for the whole range of applications. This solution can be easily incorporated, without any computational cost, inside the computational fluid dynamics solvers for fire simulations. Finally, enabling measurements of AST in any arbitrary points and directions of computational domain, makes it possible to validate numerical results with the results given by plate thermometers in real experiments with satisfactory accuracy. Because the given expression describes the physical quantity in an exact way, it may be potentially used for analyses performed for better understanding of physical process, e.g during the real experiments.