Linearization of Thermal Equivalent Temperature Calculation for Fast Thermal Comfort Prediction

: Virtual simulations and calculations are a key technology for future development methods. A variety of tools and methods for calculating thermal comfort have not gained sufﬁcient acceptance in practice due to their inherent complexity. This article investigates alternative means of determining thermal comfort, namely, the linearization of the equivalent temperature calculation. This enables a wide range of users to evaluate thermal comfort in a fast and easy manner, for example, for energy efﬁciency simulation. A ﬂow and thermal model were created according to the requirements of DIN EN ISO 14505 to determine heat transfer coefﬁcients under calibration conditions. The model to simulate the equivalent temperature in calibration conditions comprises a geometrically realistic 3D model of a human test person according to the standard. The inﬂuence of the turbulence model, as well as the inﬂuence of the equivalent temperature on the heat transfer coefﬁcient in calibration conditions, was investigated. It was found that the dependence of the equivalent temperature is mandatory. The dependence between the heat transfer and the equivalent temperature was taken into account with a continuous linearization approach. An equation-based implementation methodology is proposed, enabling a quick implementation of comfort evaluation in future simulation models. Two test cases show the capabilities of the new model and its application in future work.


Introduction
Virtual development tools will significantly change engineering methods in the next few years. Simulations will not only shorten development cycles but also reduce development costs. In many areas, a simulation-based approach is already replacing experimentation, and the proportion of simulations will continuously increase. However, the sub-discipline of technical calculation should not only be enhanced in the field of simulation techniques but also in the area of post-processing. Subjective sensations from the field of ergonomics, such as thermal comfort, still require physical testing for evaluation in many places today. Although efforts are being made in many areas to derive thermal comfort from measurable variables, this approach has not yet been fully established in the simulation-based development process of the industry. The approach of equivalent temperature (T eq ) for thermal comfort evaluation defined in the standard DIN EN ISO 14505 offers a comprehensible way of achieving this. However, recent publications such as Voelker [1] and Streblow [2] only show an evaluation for a few discrete ambient temperatures and a rather complex iterative calculation for the heat transfer coefficients (heq). In contrast, a continuous evaluation of the equivalent temperature is mandatory for real-world applications, e.g., in car interiors with many possible air temperatures. The approach described in this paper aims to simplify the heq calculation by linearization with a negligible loss in accuracy to allow for a fast thermal comfort evaluation using the concept of equivalent temperature, applicable for arbitrary real-world ambient conditions.

The DIN EN ISO 14505
The DIN EN ISO 14505 is based on the calculation and assessment of the equivalent temperature. The equivalent temperature is a physical temperature quantity based on the dry heat exchange of the human body. It combines the effects of convection and radiation on the heat exchange of the human body. The climate conditions are assessed in terms of local effects on specific parts of the body. The local equivalent temperatures determine the degree to which the relevant parts of the body parts are within acceptable levels of heat loss. The equivalent temperature allows a quantitative assessment of the heat exchange. Higher T eq values indicate lower heat losses, whereas lower T eq values indicate higher heat losses. The relationship between local T eq values for 16 body segments and the perceived thermal sensation was determined with measurements on subjects subjected to 1 h of a variety of thermal conditions. The necessary link between measurement, calculated T eq , and the subjective thermal sensation is then completed by thermal comfort diagrams. The underlying measurements in the DIN EN ISO 14505 are widely accepted in the scientific context; therefore, this work will focus on the calculation of the equivalent temperature and use the underlying data from the DIN EN ISO 14505 as validation. The calculation of the local equivalent temperature is mainly based on the work of Nilsson [6,7]. Nilsson proposed characterizing the thermal comfort by the heat transfer coefficients at the body. For this purpose, he compared test person studies with measurements and simulations to describe thermal comfort.

Related Work Regarding Equivalent Temperature Calculation and the Application of T eq Calculation
To improve the methodology further with respect to the thermal physiology of the human body, the Fiala model is often used in current research [8]. In addition, other authors and Fiala himself have also worked with nodal networks to calculate thermal physiology [9][10][11][12][13][14][15][16]. Several authors use the relation from Nilsson to calculate the thermal comfort for all kinds of applications, such as in buildings, aircraft, and transportation [17][18][19][20]. Wölki et al. [21] have even used the method as a control target value for thermal comfort. Other authors focused on another influence on thermal comfort-namely, the contact resistance between surfaces and the skin [22][23][24].

Further Development Regarding Heat Transfer and Influence on Heat Transfer
However, the equations for the thermal comfort currently used by many authors are based on the geometry of Nilsson's thermal manikin. The manikin consisted of rectangular blocks, so there is a major geometrical difference to real human body parts and, therefore, their heat transfer coefficients. Taghinia [25] studied the effect of simplification in relation to heat transfer but not in relation to thermal comfort. A study by Yousaf investigated the influence of the k-omega SST and Standard k-epsilon turbulence model on the flow field around a female manikin [26]. According to Lee et al. [27], turbulence intensity plays a major role in regard to the heat transfer at an airspeed above 0.3 m/s, but they did not put their results into perspective with thermal comfort. A similar study from Voelker et al. [1] with a coupling approach was completed on the full body and not for single body parts. Gao et al. [28] made a study on the heat transfer regression for single body parts at sitting and standing positions for a female body without a regression for the heat transfer in calibration conditions itself. Other influences such as the difference with respect to the clothing factor, but setting the convective and radiative heat transfer coefficient as a constant was investigated by several authors [29][30][31]. Streblow and Voelker accounted for the dependence of heat transfer coefficient in regard to the equivalent temperature in calibration conditions for discrete T eq , but not for a continuous T eq [1,2]. The author could not find investigations about the influence of the turbulence model in regard to the thermal comfort of a manikin. With the thermal resistance network between the skin, clothing, and surrounding air in mind, the convective and radiative thermal resistance will be important at low air velocities and the clothing factor for higher velocities. The aim of the present work is to linearize the equivalent temperature calculation to resolve the dependency between the heat transfer coefficient and the equivalent temperature to make the approach closer to the real-world behavior for low velocities. Furthermore, all of the previous studies share the necessity to implement the calculation of the thermal comfort in a rather complex manner to predict the local thermal comfort [2].

Further Development Regarding Calculation Effort
Similar work to reduce the calculation effort for the heat transfer coefficients in real conditions has been completed by several authors in different applications but not in the context of equivalent heat transfer calculation [32,33]. Table 1 gives an overview of the works that are most closely related. The remainder of the paper is structured as follows: Section 3 highlights the simulation models and gives an overview of the topic before concluding in Section 4, describing the linearization procedure. Section 5 gives some application examples as well as validation and performance benchmarks of the model.

Simulation Models
The calculation of the equivalent temperature is based on the basic assumption that the heat flow between body parts and the environment in real conditions corresponds to the heat flow in homogenous calibration conditions of a uniform enclosure. The wall temperature of this respective enclosure under these uniform (calibrated) conditions is called equivalent temperature. Comfort diagrams for each segment of the body map the equivalent temperatures to a subjective thermal comfort vote. The correlation between subjective votes and equivalent temperature is evaluated by subject testing under welldefined personal (clothing insulation, metabolic activity, etc.) and climatic test conditions.
The projection between the real conditions and the calibration conditions is realized by evaluating the heat transfer coefficient h eq , defined in calibration conditions as with respect to the heat flow . Q eq , temperature in calibration conditions T eq , and the skin temperature T s .
For a constant skin temperature, emission factors, and inflow velocity, the emitted heat flow is only dependent on the temperature of the calibration conditions. Thus, the heat transfer coefficient in calibration conditions is only dependent on the equivalent temperature. In order to calculate the equivalent temperature in real conditions, the heat flow in working conditions . Q real = . Q eq = .

Q
(4) must correspond to that in equivalent conditions. Equation (1) can be solved using Equations (3) and (4) according to the equivalent temperature to The calculated equivalent temperature is inserted into the DIN EN ISO 14505 comfort diagrams, and the thermal comfort can be assessed from the resulting diagram.

Numerical Fluid Dynamics Model
For the determination of the heat transfer resistances in calibration conditions, thermal fluid flow simulations are used according to the standard DIN EN ISO 14505 [5]. The continuity equation results from observations on an infinitesimal control volume dV with the velocity v and the density ρ The momentum conservation in integral form for the fluid is given with the identity tensor I, the pressure p, the viscous stress tensor T, and body forces f b as [34] Energies 2021, 14, 5922

of 25
The energy conservation of the infinitesimal volume element for the total energy per unit mass E can be written under consideration of the Fourier heat conduction and a source term, as [34] Modeling the medium as an ideal gas allows for the consideration of the influence of pressure and temperature on the air density [34] to consider the buoyancy effects in the flow field. The mass, momentum, and energy conservation equations result in a non-linear equation system, the Navier-Stokes equations. The equations are solved numerically with second-order upwind discretization on a grid with polyhedron elements. Due to the strong coupling between velocity and temperature field, the equations have to be solved by an implicitly coupled approach.
The surface-to-surface radiation is calculated by ray tracing using the enclosure theory with view factors. The radiation power P 1−2 from one surface patch dS 1 to another dS 2 is where β is the angle between the surface normal and a line between the two surfaces, and the length L of this line. The total intensity i 1 is defined as the radiative energy passing through an area per unit solid angle, per unit of the area projected normal to the direction of passage, and per unit of time.
The view factor F is defined as the ratio of the total radiation emitted by patch 1 to the radiation received by patch 2 The calculation of the view factor F 1−2 between the surfaces is based on the topological conditions and is calculated by The approximation of the integral is completed by ray tracing, where each patch sends out a specified number of beams at a discretized hemisphere over the patch.
The view factors are calculated once in the initialization phase. During the calculation, the heat exchange is iteratively calculated based on the view factor matrix, assuming a radiative equilibrium. The heat flow exchanged by radiation is used as a boundary condition on the faces. For more details on the solving approach, see [35] or [36].
The equivalent heat transfer coefficient h eq , defined in Equation (13), is based on the total exchanged heat flow, incorporating the sum of convection, conduction in the thermal boundary layer, and radiation. The reference temperature for the heat transfer coefficient is the equivalent temperature, i.e., the wall temperature of the room.

Structure of the Simulation Model
The manikin geometry resembles a 95-percentile male RAMSIS model geometry and is geometrically processed as a solid body [37]. The subdivision of the body areas was carried out as suggested by Nilsson, corresponding to the comfort diagrams in DIN EN ISO 14505 [4,5], see Table 2. The grid independence of the mesh was investigated with a study of four different base sizes, using the total heat flux of the body as an indicator for independence; see Figure 1.

Structure of the Simulation Model
The manikin geometry resembles a 95-percentile male RAMSIS model geometry and is geometrically processed as a solid body [37].
The subdivision of the body areas was carried out as suggested by Nilsson, corresponding to the comfort diagrams in DIN EN ISO 14505 [4,5], see Table 2. x Whole Body The grid independence of the mesh was investigated with a study of four different base sizes, using the total heat flux of the body as an indicator for independence; see   The surface of the body is meshed with a base size of 5 mm polyhedron elements. The volume mesh uses a growth factor of 1.2 until reaching 30 mm in 2 m distance to the RAMSIS. To account for realistic flow behavior around the RAMSIS model, a prism layer mesh with 12 layers is included. The total mesh results in about 2 million cells.

Simulation and Linearization of the Calibration Conditions
The flow diagram of the linearization procedure that the authors used can be seen in Figure 2.

Simulation and Linearization of the Calibration Conditions
The flow diagram of the linearization procedure that the authors used can be seen in Figure 2. For the determination of the heat transfer coefficients, the following conditions are specified by the standard DIN EN ISO 14505-2 [5]. For a constant skin temperature of 34 °C, different equivalent temperatures were specified to determine the dependence of the heat transfer coefficient. Because the heat flux depends on the temperature difference between the skin and the surroundings where the equivalent temperature is evaluated, it is possible to account for inhomogeneous temperature distributions when using a specific skin temperature in Equation (24). In this work, a homogenous temperature is used for simplification. Figure 3 shows the heat flux on the RAMSIS geometry and the velocity fields, which result at an equivalent temperature of 25 °C. Three characteristic rays are formed as wake areas, in which the temperature of the flow in +z direction has a longer length scale until mixing with the direct surroundings. The jets (white arrows in Figure 3) are significantly longer than the distance the flow travels along the body. The flow heats up faster per unit of travel as it passes the body and then releases the heat through mixing in the wake zone. The wake formation is induced by the head and the knee area. In addition, a small trailing zone is created in the vicinity of the feet. The trailing zones appear due to the slow flow speeds, as the air follows the body shape and then leaves at geometrically characteristic edges. The representation of the heat flows can be used to derive initial model-like impressions of the temperature sensation. For example, a test person in this configuration gives off more heat via the stomach than via the sides.
Calculate the heatflux in calbiration conditions according to the norm for various equivalent temperatures.
Calculate h eq as a function of T eq for the results of step 1. Use a polynomal regression h eq = f (T eq )= x 1 * T eq +x 0 for the dataset Calculate the heatflux for the conditions you want to predict the thermal comfort. İ nsert the h eq equation, and the known skin temperature into the equation for thermal comfort according to DIN EN ISO 14505-2 , = +ḣ For a constant skin temperature of 34 • C, different equivalent temperatures were specified to determine the dependence of the heat transfer coefficient. Because the heat flux depends on the temperature difference between the skin and the surroundings where the equivalent temperature is evaluated, it is possible to account for inhomogeneous temperature distributions when using a specific skin temperature in Equation (24). In this work, a homogenous temperature is used for simplification. Figure 3 shows the heat flux on the RAMSIS geometry and the velocity fields, which result at an equivalent temperature of 25 • C. Three characteristic rays are formed as wake areas, in which the temperature of the flow in +z direction has a longer length scale until mixing with the direct surroundings. The jets (white arrows in Figure 3) are significantly longer than the distance the flow travels along the body. The flow heats up faster per unit of travel as it passes the body and then releases the heat through mixing in the wake zone. The wake formation is induced by the head and the knee area. In addition, a small trailing zone is created in the vicinity of the feet. The trailing zones appear due to the slow flow speeds, as the air follows the body shape and then leaves at geometrically characteristic edges. The representation of the heat flows can be used to derive initial model-like impressions of the temperature sensation. For example, a test person in this configuration gives off more heat via the stomach than via the sides.  As a validation of the flow field, the qualitative results of PIV measurements by Mitterhofer et al. [38] were used. In their experiment, the same phenomena could be observed as in this study.

Turbulence Models
In fluid flow simulations, different methods are used to solve the non-linear Navier-Stokes equation system. A huge challenge is the correct calculation of turbulence effects using semi-empirical models and assumptions. In the industry, the RANS method (Reynolds-averaged Navier-Stokes) for turbulence modeling has become generally accepted, in which the time-dependent quantities velocity, pressure, and temperature are divided into an average and a fluctuation term and are then averaged over time in order to eliminate the fluctuation terms. A good overview of the topic can be found in the book of Ferziger and Perić [34].
Due to the averaging, Reynolds stress terms are introduced in the equation system, which has to be modeled to solve the system of equations. The equation system contains more unknown quantities than equations; hence, the Reynolds stress needs to be modeled. A common model is the Boussinesq approximation that treats the Reynolds stresses as viscous stresses [34], which defines the Reynolds stress tensor to Different methods for turbulence modeling have been established, where k-epsilon and k-omega models in different variations are most common. Kappa or k is the turbulent kinetic energy, epsilon is the turbulent dissipation rate, and omega is the specific dissipation rate. As a validation of the flow field, the qualitative results of PIV measurements by Mitterhofer et al. [38] were used. In their experiment, the same phenomena could be observed as in this study.

Turbulence Models
In fluid flow simulations, different methods are used to solve the non-linear Navier-Stokes equation system. A huge challenge is the correct calculation of turbulence effects using semi-empirical models and assumptions. In the industry, the RANS method (Reynolds-averaged Navier-Stokes) for turbulence modeling has become generally accepted, in which the time-dependent quantities velocity, pressure, and temperature are divided into an average and a fluctuation term and are then averaged over time in order to eliminate the fluctuation terms. A good overview of the topic can be found in the book of Ferziger and Perić [34].
Due to the averaging, Reynolds stress terms are introduced in the equation system, which has to be modeled to solve the system of equations. The equation system contains more unknown quantities than equations; hence, the Reynolds stress needs to be modeled. A common model is the Boussinesq approximation that treats the Reynolds stresses as viscous stresses [34], which defines the Reynolds stress tensor τ to Different methods for turbulence modeling have been established, where k-epsilon and k-omega models in different variations are most common. Kappa or k is the turbulent kinetic energy, epsilon is the turbulent dissipation rate, and omega is the specific dissipation rate. The k-epsilon model is considered in both standard and Abe-Kondoh-Nagano formulation with a modified damping term near the wall. According to the literature, this model is particularly suitable for the simulation of heat transfer phenomena, as the model is derived for low Reynolds numbers [35,39].

k-Epsilon Model
The k-omega model was originally derived from Wilcox and then further developed by Menter [40]. In this work, the industry standard for turbulent flows, the k-omega model, in the form of the well-known Menter SST model, is used [41].

Laminar Model
In the laminar model, the turbulent friction terms in the Navier-Stokes equations are not calculated. Therefore, these flows are similar to potential flows.

Turbulent Viscosity
Depending on the model, the turbulent viscosity is calculated with the transported quantities.

Model Comparison
The models are compared with approximate empirical equations from the work of Nilsson. In Nilsson's work, Nusselt correlations for simple bodies, such as cuboids and cylinders, are used to determine the heat transfer within a flow model [6].
If the heat transfer from the flow simulation is compared to the heat transfer from Nilsson [6], it is obvious that there are differences for the extremities such as arms and legs. The scatter within the turbulence models is smaller than the deviation to the empirical equations of Nilsson but still around 15%; it is assumed that the difference to Nilsson can be explained because of the more simple geometry Nilsson used-see Figure 4. This deviation could be attributed to the different positions of the arms and legs. In the Nilsson investigation, horizontal legs and arms are assumed, whereas in this work, the arms and legs are assumed to be in an ergonomic seating position. Therefore, an exact geometric mapping of the conditions is necessary to determine the heat transfer in calibration conditions. For all turbulence models, the qualitative course corresponds to that of the empirical equations. Therefore, it can be concluded that a realistic geometry surface mesh is more important than the turbulence model used. For the calculations in this work, the kω-SST all y+ model was used.
The heat transfer coefficients are averaged over the corresponding body part surface. Figure 5 shows the different heat transfers across the individual body parts for different equivalent temperatures.  The heat transfer coefficients are averaged over the corresponding body part surface. Figure 5 shows the different heat transfers across the individual body parts for different equivalent temperatures.   The heat transfer coefficients are averaged over the corresponding body part surface. Figure 5 shows the different heat transfers across the individual body parts for different equivalent temperatures. If the heat transfer coefficients are plotted against the absolute magnitude of the temperature difference between the equivalent and the skin temperature, it can be observed that the heat transfer in an equivalent temperature range from −25 to 25 • C can be well approximated with the square of the temperature, as Figure 6 shows. For a better representation, different body parts are grouped together; the simulated values are shown as points and the approximation as a line in the same color. The coefficient of determination R 2 is above 0.99. The difference in the heat transfer coefficient in regard to the equivalent temperature is around 10%, between 20 and 25 • C. For example, when a calculated comfortable equivalent temperature for a 34 • C body part is 25 • C, not incorporating this behavior would lead to an error that is almost 50% of the width of the comfort area. If the heat transfer coefficients are plotted against the absolute magnitude of the temperature difference between the equivalent and the skin temperature, it can be observed that the heat transfer in an equivalent temperature range from −25 to 25 °C can be well approximated with the square of the temperature, as Figure 6 shows. For a better representation, different body parts are grouped together; the simulated values are shown as points and the approximation as a line in the same color. The coefficient of determination R² is above 0.99. The difference in the heat transfer coefficient in regard to the equivalent temperature is around 10%, between 20 and 25 °C. For example, when a calculated comfortable equivalent temperature for a 34 °C body part is 25 °C, not incorporating this behavior would lead to an error that is almost 50% of the width of the comfort area. The parameters a, b, and c are shown in Table 3 and Figure 7 as a result of the quadratic approximation.   A coefficient matrix K contains the parameters for the quadratic regression for each body segment via the index j and is constructed as follows.
With the coefficient matrix K, the equivalent heat transfer coefficient can be calculated as A coefficient matrix K contains the parameters for the quadratic regression for each body segment via the index j and is constructed as follows.
With the coefficient matrix K, the equivalent heat transfer coefficient can be calculated as h eq,j = T T eq,j · K j · T eq,j If the approach is inserted into Equation (5), the following equation for the equivalent temperature of the different parts of the body is obtained: Since the goal of the present work is to reduce the effort for the T eq calculation, the approximation degree of h eq is reduced to a linear approach. h eq = x j,1 · T eq,j + x j,0 The following coefficients x j,1 and x j,0 result from linear approximation with R 2 > 0.95, see Table 4. In the diagrams (Figures 8 and 9), the linear approximation of the h eq is shown for the two body groups. The coefficient of determination is R 2 > 0.96. ℎ = ,1 ⋅ , + ,0 (2 The following coefficients xj,1 and xj,0 result from linear approximation with R² >0. see Table 4. In the diagrams (Figures 8 and 9), the linear approximation of the heq is shown for t two body groups. The coefficient of determination is R² > 0.96.  The linear approximation is, however, less accurate than the quadratic approximation. Nevertheless, the errors introduced by the linear approximation are smaller than the variations due to measurement errors in the calibration chamber or the variance within the turbulence models. As well, variations in the assessment of thermal comfort votes by subject testing in a climatic chamber are large compared with the linear approximation error. Hence, the accuracy of the presented linear approximation is deemed sufficient for the presented context.
The formula for the equivalent temperature using linear approximation results in , = + ,1 ⋅ , + ,0 After rearranging the equation according to the unknown temperature , and solving the resulting new equation, the equivalent temperature for the individual body parts can be written as The implementation in common flow simulations as user code is now possible as soon as only the physically reasonable values of the positive root are considered. Additionally, for implementation purposes, a limiter can be introduced to limit the denominator of the terms in the free flow range.
In order to implement the equivalent temperature approach, the two parameters are assigned to the corresponding body parts, and then a new scalar is created using Equation (26). The heat flow of the respective face element can be used. The linear approximation is, however, less accurate than the quadratic approximation. Nevertheless, the errors introduced by the linear approximation are smaller than the variations due to measurement errors in the calibration chamber or the variance within the turbulence models. As well, variations in the assessment of thermal comfort votes by subject testing in a climatic chamber are large compared with the linear approximation error. Hence, the accuracy of the presented linear approximation is deemed sufficient for the presented context.
The formula for the equivalent temperature using linear approximation results in T eq,j = T s + Q j x j,1 · T eq,j + x j,0 After rearranging the equation according to the unknown temperature T eq,j and solving the resulting new equation, the equivalent temperature for the individual body parts can be written as The implementation in common flow simulations as user code is now possible as soon as only the physically reasonable values of the positive root are considered. Additionally, for implementation purposes, a limiter can be introduced to limit the denominator of the terms in the free flow range.
T eq,j = Q + T s · x j,0 min x j,1 , 10 −5 + x j,0 − T s · x j,1 2 · min x j,1 , 10 −5 2 − x j,0 − T s · x j,1 2 · min x j,1 , 10 −5 (24) In order to implement the equivalent temperature approach, the two parameters are assigned to the corresponding body parts, and then a new scalar is created using Equation (26). The heat flow Q of the respective face element can be used.

Validation
In the following section, the application and validation of the linearized approach from Section 3 are demonstrated for thermal comfort assessment of real case scenarios in two exemplary applications. It will be demonstrated that both simulation experiments can be improved with respect to local thermal comfort prediction without adding additional simulation complexity by applying the linearized equivalent temperature approach. Each scenario contains numerical manikins for thermal comfort evaluation. The results from these experiments are validated and are therefore used to validate the approach presented in this work as well.

Model Setup
Atish and Upender [42] published a study to predict thermal comfort with the Simplified Thermoregulatory Bio-heat (STB) equation for a cubic room with non-symmetric thermal boundaries. The authors varied the inlet temperature from a cubic room in 10 K steps between 10 and 50 • C while keeping the right-hand side (RHS) of the room at 40 • C and the left-hand side (LHS) of the room at 10 • C, both sides from the perspective of the manikin. The body temperature is changing based on the computation of the STB equation according to the original work. The thermal comfort in the original study was evaluated using Fanger's model for the steady-state solutions.
The local thermal comfort for this test case is calculated using the linearized equivalent temperature calculation. Based on these results, the global thermal comfort has been calculated for validation purposes. Since no radiation was implemented in the original model and Fanger's model is used for a global thermal comfort calculation, the prediction of local thermal comfort is firstly introduced with the linearized equivalent temperature calculation in the following section. The computational mesh consists of polyhedral cells with a base size of 10 mm in the computational domain. A prism-based mesh was used on all elements near the surfaces except at the inlet and outlet. The manikin surfaces have a custom control, as shown in Figure 10, with a finer size of 5-50 mm. In advance of the final results, a grid study has been completed to confirm grid independence for the heat fluxes. The kω-SST (Menter) turbulence model is used for the cubic room in combination with segregated fluid temperature and surface-to-surface radiation because it is commonly used in similar applications. The air is modeled as an ideal gas with a segregated fluid temperature, gravity, and surface-to-surface radiation. For the thermal comfort eval-  The kω-SST (Menter) turbulence model is used for the cubic room in combination with segregated fluid temperature and surface-to-surface radiation because it is commonly used in similar applications. The air is modeled as an ideal gas with a segregated fluid temperature, gravity, and surface-to-surface radiation. For the thermal comfort evaluation, the proposed method is used.
The boundary conditions are according to [42]: • A constant velocity inlet with 0.5 m/s in front of the manikin; • An outlet behind the manikin; • Left-hand side wall at 10 • C, right-hand side wall at 40 • C; the remaining walls are adiabatic.

Global Results
After calculating the equivalent temperature, the PMV for each body segment was calculated based on the formulas from Tomschke [43] and averaged over the body. Then, the PPD was calculated based on the formula according to ASHRAE [3].
Comparing the curve for the PPD in relation to the inlet temperature, both slopes look similar. The position of the minimum PPD is comparable with the linearized T eq approach when simulating with thermal radiation to the STB simulation without radiation. The biggest difference is at about 30 • C, where the PPD with STB is 43 %, and the PPD for the linearized temperature approach is 27 %, as shown in Figure 11. The PPD for the linearized approach looks symmetrical around the minimum PPD, which is at 30 • C inlet temperature. The steeper increase in the PPD around that temperature in the STB model could be related to the missing radiation in the setup of the original test case or physiology aspects. The linearized equivalent temperature approach shows, for the global thermal sensation, almost the same results as the STB model with a coupled PMV and PPD calculation. The following chapter will show the error of the STB model in relation to local thermal comfort due to the missing radiation.

Local Thermal Comfort
Because of the neglected radiation, the original work by Atish and Upender [42] did not study the difference between the left-hand side and right-hand side of the body, even if the boundary conditions were asymmetrical. Figure 12 shows the calculated equivalent temperature for five different inlet temperatures for each body segment. For a better in-

Local Thermal Comfort
Because of the neglected radiation, the original work by Atish and Upender [42] did not study the difference between the left-hand side and right-hand side of the body, even if the boundary conditions were asymmetrical. Figure 12 shows the calculated equivalent temperature for five different inlet temperatures for each body segment. For a better interpretation of the plot, dotted lines connect the body parts for every inlet temperature. According to Figure 11, the lowest PPD for the test persons is at an inlet temperature of 25 • C. In the local thermal comfort diagram, a difference between the left-hand side and the right-hand side for each extremity can be seen. The difference between larger extremities such as the LHS and RHS thigh and upper arm is higher than for smaller extremities such as the hands, feet, or calves. Keeping this observation in mind, a PMV of -1 for the left side and a PMV of +1 for the right side would result in a PMV of 0, which would suggest false results with respect to thermal sensation.   Figure 12 shows the difference in the equivalent temperature for the LHS and RHS of the manikin at the comfort diagram. The equivalent temperature at the manikin surfaces that faces the colder LHS is between 5 and 10 degrees lower than at the side facing the hotter RHS. Even at a specific extremity, a difference can be seen in the surface plot in Figure 13, which would probably be neglected when simply averaging over an extremity. The difference is based on a dot product between the surface normal of a location on the manikin surface and the incoming radiation in the Y-direction.  The heat flux for the manikin is defined as positive for receiving incoming radiation and negative for outgoing radiation. Figure 14 shows these radiation heat fluxes at the manikin's surface, which can radiate about −55 W/m² to the colder LHS, while it receives 5 W/m² from the warmer side. The heat flux for the manikin is defined as positive for receiving incoming radiation and negative for outgoing radiation. Figure 14 shows these radiation heat fluxes at the manikin's surface, which can radiate about −55 W/m 2 to the colder LHS, while it receives 5 W/m 2 from the warmer side. Comparing the radiation heat fluxes to the total heat flux at the boundary in Figure  15, the amount of heat the manikin received by radiation can be concluded to be a major driver for the thermal comfort in this specific case. The results for a global thermal sensation showed nearly no differences between the linearized, simplified approach and the detailed, complex model. Looking at the local thermal comfort, however, the results are diverging. With the linearized thermal comfort prediction approach for a manikin in the cubic room, a difference in thermal comfort could be examined between the LHS and RHS setup. The differences between the two sides of the body are related to the difference in radiative heat flux based on the asymmetrical thermal boundary conditions at the walls. Comparing the radiation heat fluxes to the total heat flux at the boundary in Figure 15, the amount of heat the manikin received by radiation can be concluded to be a major driver for the thermal comfort in this specific case. The results for a global thermal sensation showed nearly no differences between the linearized, simplified approach and the detailed, complex model. Looking at the local thermal comfort, however, the results are diverging. With the linearized thermal comfort prediction approach for a manikin in the cubic room, a difference in thermal comfort could be examined between the LHS and RHS setup. The differences between the two sides of the body are related to the difference in radiative heat flux based on the asymmetrical thermal boundary conditions at the walls. driver for the thermal comfort in this specific case. The results for a global thermal sensation showed nearly no differences between the linearized, simplified approach and the detailed, complex model. Looking at the local thermal comfort, however, the results are diverging. With the linearized thermal comfort prediction approach for a manikin in the cubic room, a difference in thermal comfort could be examined between the LHS and RHS setup. The differences between the two sides of the body are related to the difference in radiative heat flux based on the asymmetrical thermal boundary conditions at the walls.

Setup
Cook et al. [44] describe a realistic scenario of a rectangular environmental chamber with a fan and a coupled simulation model to test and validate human comfort. The setup is adopted, and unknown settings are estimated as given in Table 5. The k-omega SST (Menter) turbulence model is used with surface-to-surface radiation, ideal gas for air,

Environmental Chamber Setup
Cook et al. [44] describe a realistic scenario of a rectangular environmental chamber with a fan and a coupled simulation model to test and validate human comfort. The setup is adopted, and unknown settings are estimated as given in Table 5. The k-omega SST (Menter) turbulence model is used with surface-to-surface radiation, ideal gas for air, gravity, and a coupled solver to run a steady simulation. The fan is modeled as a simplified cylindrical body with an intake and outlet on the frontal surfaces connected by a fan interface with a generic fan curve. Transient effects such as the fan swirl are neglected. The mesh consists of polyhedral cells with prism layer cells towards all surfaces except for the ceiling opening, the wall inlet, and the fan front and back faces. The number of prism layers is eight on the manikin and four for every other surface. The mesh base size is 50 mm with a target/minimum size of 15/10 mm on the manikin and fan surfaces. A cylindrical volumetric mesh refinement with 20 mm is used for the fan jet. These sizes result in a mesh containing around 700,000 volume cells, depicted in a plane section in Figure 16a. On the manikin surface, there are eight prism layers, as shown in Figure 16b. The final mesh settings are the results of a mesh study to find mesh independence. cept for the ceiling opening, the wall inlet, and the fan front and back faces. The number of prism layers is eight on the manikin and four for every other surface. The mesh base size is 50 mm with a target/minimum size of 15/10 mm on the manikin and fan surfaces. A cylindrical volumetric mesh refinement with 20 mm is used for the fan jet. These sizes result in a mesh containing around 700,000 volume cells, depicted in a plane section in Figure 16a. On the manikin surface, there are eight prism layers, as shown in Figure 16b. The final mesh settings are the results of a mesh study to find mesh independence.

Global Results
The vertical temperature distributions in the chamber in Figure 17 result from two steady simulations with the fan on or off. Without the fan, a distinct plume of hot air is rising up from the manikin, and the temperature stratification is stronger than with a running fan, similar to the results from Cook et al. [44].

Global Results
The vertical temperature distributions in the chamber in Figure 17 result from two steady simulations with the fan on or off. Without the fan, a distinct plume of hot air is rising up from the manikin, and the temperature stratification is stronger than with a running fan, similar to the results from Cook et al. [44]. With an activated fan, the air temperature is slightly higher and more uniform throughout the chamber. The fan jet is oriented towards the manikin and partially increases the air velocity ( Figure 18) around some of the manikin's body segments.

Thermal Comfort
The calculation of the equivalent temperature Teq yields an overall thermally neutral comfort level for the manikin with a deactivated fan around the neutral line ( Figure 19). The hands tend to have higher Teq values because of their proximity to both legs ( Figure  20). When the fan is activated, the predicted comfort level is lower. Especially, the righthand side calf and left-hand side knee show low Teq values (Figure 19), which translate to a too-cold and cold comfort level for the right-hand side calf and left-hand side thigh, respectively. This is mainly due to the increased heat transfer caused by the higher veloc- With an activated fan, the air temperature is slightly higher and more uniform throughout the chamber. The fan jet is oriented towards the manikin and partially increases the air velocity ( Figure 18) around some of the manikin's body segments. With an activated fan, the air temperature is slightly higher and more uniform throughout the chamber. The fan jet is oriented towards the manikin and partially increases the air velocity ( Figure 18) around some of the manikin's body segments.

Thermal Comfort
The calculation of the equivalent temperature Teq yields an overall thermally neutral comfort level for the manikin with a deactivated fan around the neutral line ( Figure 19). The hands tend to have higher Teq values because of their proximity to both legs ( Figure  20). When the fan is activated, the predicted comfort level is lower. Especially, the righthand side calf and left-hand side knee show low Teq values (Figure 19), which translate to a too-cold and cold comfort level for the right-hand side calf and left-hand side thigh, respectively. This is mainly due to the increased heat transfer caused by the higher veloc-

Thermal Comfort
The calculation of the equivalent temperature T eq yields an overall thermally neutral comfort level for the manikin with a deactivated fan around the neutral line ( Figure 19). The hands tend to have higher T eq values because of their proximity to both legs ( Figure 20).
When the fan is activated, the predicted comfort level is lower. Especially, the right-hand side calf and left-hand side knee show low T eq values ( Figure 19), which translate to a too-cold and cold comfort level for the right-hand side calf and left-hand side thigh, respectively. This is mainly due to the increased heat transfer caused by the higher velocity fan jet directly hitting these body segments. The arms, hands, and legs are cold but still comfortable, with the feet, right upper arm, upper back, and head being thermally neutral. From these Teq values, a global PMV and PPD index can be calculated, such as the previous generic cubic room study. The PMV value yields a positive, neutral comfort level with a deactivated fan and a tendency to a negative slightly cold comfort level with an activated fan (Table 6).  From these Teq values, a global PMV and PPD index can be calculated, such as the previous generic cubic room study. The PMV value yields a positive, neutral comfort level with a deactivated fan and a tendency to a negative slightly cold comfort level with an activated fan (Table 6).  From these T eq values, a global PMV and PPD index can be calculated, such as the previous generic cubic room study. The PMV value yields a positive, neutral comfort level with a deactivated fan and a tendency to a negative slightly cold comfort level with an activated fan (Table 6). To compare the Dynamic Thermal Sensation (DTS) values from Cook et al. [44], their transient values were extrapolated with an exponential function (Equation (27)) to estimate a final value that can be compared to a steady-state simulation. This results in the final DTS values of −0.19 without a fan and −0.47 with a fan. Despite the transient formulation of the DTS index, it uses the same 7-point ASHRAE scale as the PMV index ranging from −3 to +3 and can therefore be compared with each other in terms of a thermal comfort prediction using the final value over a long period of time. This leads to a different comfort prediction when focusing on the global comfort level, where both settings (with/without fan) show a negative DTS value equivalent to a neutral to slightly cold comfort level, and the PMV values from T eq are positive or negative. However, the lower comfort level with the activated fan is clearly visible in this T eq and Cook et al.'s [44] approach.

Computational Runtimes
The total runtime of the steady-state simulations with 1000 iterations is around 50 min on an Intel ® Xeon ® Gold Linux workstation with 16 cores. Using the High Performance Computing Center Stuttgart (HLRS) HAWK cluster, the total runtime can be further reduced to, e.g., 7 min on 256 cores. Cook et al. [44] had a longer runtime, especially with the fan running about 1.5 days on 12 cores. This is a significant reduction for a fast comfort prediction with the T eq approach compared to the DTS method (Table 7), with respect to achieving similar results in terms of the PPD comfort index.

Summary
In the DIN EN ISO 14505, the equivalent temperature method is proposed for the assessment of thermal comfort. Comfort diagrams, based on the work of Nilsson, indicate comfort in real conditions by means of the equivalent temperature [8]. This paper presented an approach to simplify the calculation of the equivalent temperature significantly by linearizing the heat transfer coefficient calculation based on the calculated heat transfer resistances in calibration conditions according to the standard DIN EN ISO 14505 for a temperature window of −25 to +25 • C. Subsequently, a quadratic and linear approach is chosen for the approximation of the heat transfer coefficients. The linear approach shows a deviation of less than 5% in the heat transfer coefficients, which is small compared to the influences of different turbulence models or inaccuracies in the comfort assessment in general. Two test cases show the application of the implemented linearization for a generic cubic room and an environmental chamber, respectively. The first test case is extended by a local thermal comfort prediction with the linearized T eq approach. The second test case is used to show the significant reduction in computation time for a faster thermal comfort prediction.

Conclusions
The linearization approach helps to reduce the effort to predict the thermal comfort considerably for the user, which will allow a broader application in the future. This work shows the procedure to linearize the thermal comfort prediction as well as the implementation in a state-of-the-art CFD software. The simulation of thermal comfort will, in the future, tend to have more complex models with the influence of moisture and thermophysiology. All of these effects will be easy to implement in the presented approach using the linearization according to the presented method. Comparing the results with the results of an STB thermophysiology model shows that the contribution in the calibration conditions is comparably small. Effects such as sweating, shivering, and higher dimensional effects that do not depend on the equivalent temperature directly can be neglected in the linearization approach as long as they represent only a minor contribution to the calibration conditions. To implement these effects of the human body for the working conditions, they can just be incorporated in the application simulations while keeping the linearization approach for the equivalent conditions with the linearization coefficients presented in this work. Nevertheless, it should be noted that there may be environmental conditions, such as high humidity, where the effects play a major role. Comparing the results for the environmental chamber with respect to the calculation time, the linearized approach shows its advantages above coupled approaches like the DTS method. In summary, the presented linearization is the first step to make thermal comfort predictions more attainable for a larger audience. Because no partial differential equation has to be solved, the linearization method improves the application of the equivalent temperature method further because of its easy implementation. Thanks to the increased performance, the linearized T eq approach enables more complex thermal comfort simulations in the future, such as the interaction of many occupants in a room, for example. The performance increase is not a disadvantage in relation to the accuracy. It can be shown that the influence of the turbulence model on the heat transfer is higher than the influence of the linearization. Rather, the linearization includes the dependence of the heat transfer on the equivalent temperature, which has often been neglected so far. This approach enables future work to consider the simulative thermal comfort evaluation with an easy-to-implement calculation approach. Since this work is only focused on the simulation practice of a thermal comfort prediction, there should be additional experimental validation to demonstrate the real-world accuracy. As the described approach is in accordance with the DIN EN ISO 14505, any possible difference can only arise from the minor loss in accuracy due to the linearization of the heat transfer coefficients h eq or the standard itself. Possible real-world differences could arise from environmental specific effects on the aerodynamics, time-dependent and varying skin temperatures, and the consideration of human-related differences in clothing, stature, and body shape, which are not fully considered in the standard.