## Abstract

The infrared optical properties of textiles are of great importance in numerous applications, including infrared therapy and body thermoregulation. Tuning the spectral response of fabrics by the engineering of composite textile materials can produce fabrics targeted for use in these applications. We present spectroscopic data for engineered polyester fabric containing varying amounts of ceramic microparticles within the fiber core and report a spectrally-dependent shift in infrared reflectance, transmittance and absorptance. A thermal transport model is subsequently implemented to study the effect of these modified properties on the spectral distribution of infrared radiation incident upon the wearer of a garment constructed of this fabric.

© 2017 Optical Society of America

## 1. Introduction

The engineering of textiles to exhibit desired optical properties has long been explored, with success demonstrated in various methods including: co-spinning of different materials, introduction of gradients in fiber cross-sectional shape, inclusion of inorganic content, or use of dyes and other additives [1,2]. Recently, increased attention has been paid to the interactions between textile products and infrared radiation, due to the potential impact in applications ranging from body cooling to industrial drying processes to infrared therapy [3–6]. Mid-infrared optical properties are of particular importance when considering interactions with the human body, as at nominal skin temperature much of the body’s emissive radiative power is centered in the mid-infrared between 7 and 14 µm, and infrared heat losses account for approximately 50% of body cooling in typical indoor conditions [7,8].

In the work presented herein, we utilized Fourier transform infrared (FTIR) spectroscopy to assess the spectral optical properties – namely the near-normal reflectance, transmittance, and emittance – of textile fabrics knitted with varying percentages of ceramic-bearing polymeric fibers. The biological effects of garments manufactured from these fabrics are presently under investigation [6]. The fabrics, which we also studied in our previous work [3], were nominally identical in thickness, basis weight, knit structure, and color. All fabrics consisted of 8% elastane fibers and 92% polyethylene terephthalate (PET, i.e., polyester) based fibers; the only difference between the samples was the fraction of these fibers that contained a core of polyester blended with ceramic particles (e.g., titanium dioxide) for optical property modification.

The measured optical properties were subsequently utilized in a first-principle-based heat transfer model, accounting for emitted, reflected and transmitted radiant energy as well as thermal transport via conduction and convection, to determine the magnitude and spectral distribution of infrared energy received by a wearer of each fabric in the 0.8 – 16.7 µm spectral region. It is demonstrated that with increasing added ceramic content, the modified textile fabrics reflect and transmit less infrared energy and thus absorb more. This is particularly evident at wavelengths shorter than 6 µm. As a result, the modified fabrics with added ceramic content can absorb incoming radiation from the sun in the near-IR region, and emit increased levels of infrared energy at longer wavelengths, as compared to the same fabric without ceramic content. Thermal modeling indicates that this effect persists over a broad range of environmental and fit factors; such a fabric could be utilized when the physical properties of a polyester blend warrant its selection over alternative fabric materials, yet modified near- and mid-infrared optical properties are desired.

## 2. Experimental methods

A Nicolet iS50 FTIR spectrometer (Thermo-Fisher Scientific) was used to measure the spectral reflectance (ρ_{λ}) and transmittance (τ_{λ}) of the fabric samples with differing percentages of ceramic-bearing fibers. Because of the diffuse nature of reflection and transmission from composite surfaces [9], a mid-IR IntegratIR^{TM} integrating sphere (Pike Technologies) with a reflective gold inner coating and a deuterated triglycine sulfate detector was used to measure the reflected or transmitted infrared radiation from the fabric samples. Mid-IR measurements were conducted across a wavenumber range of 4000 – 600 cm^{−1} (corresponding to wavelengths, λ, ranging from 2.5 – 16.7 µm), using a potassium bromide (KBr) beam splitter. Near-IR measurements, ranging from wavenumbers of 12000 – 4000 cm^{−1} (0.8 < λ < 2.5 µm) were conducted using a near-IR Nicolet iS50 NIR integrating sphere module (Thermo-Fisher Scientific) with an indium gallium arsenide detector.

A schematic of the FTIR setup, indicating the sample placement relative to the IR source and detector for the various reflectance and transmittance measurements is shown in Fig. 1.

Because the near-IR integrating sphere did not allow placement of the sample directly between the IR source and the sphere inlet, the transmittance was instead assessed indirectly by comparing the reflectance where transmitted IR was permitted to escape through the sample (ρ_{λ}) to the reflectance with a reflective backing behind the sample (ρ′_{λ}). The spectral transmittance can be shown to be related to the ρ_{λ} and ρ′_{λ} reflective measurements illustrated in Fig. 1(a) and Fig. 1(c), respectively, by:

Details of this relationship are included in the Appendix. The methodology of calculating transmittance from the reflectance measurements with and without a 100% reflective backing was validated through comparison of indirect to direct transmittance measurements with the mid-IR integrating sphere setup where such direct measurement was possible.

To minimize the effect of absorption by water vapor and other trace gas compounds, the integrating sphere was surrounded with desiccant and the interior was swept with 5 scfh (standard cubic feet per hour) of inert nitrogen gas during all measurements.

Four different fabric samples were examined, each containing 0, 0.82, 0.91 or 1.18 weight percent inorganic content from modified ceramic-bearing fibers. For each fabric sample, absorption and transmission spectra were calculated by averaging the respective spectra measured for 5 unique specimens. For each specimen, the absorption and transmission spectra were obtained from 256 scans at a resolution of 16 cm^{−1}.

## 3. Experimental results

The results of the FTIR measurements are reported in this section. Although reflectance and transmittance are often used interchangeably with reflectivity and transmissivity, respectively, the terminology ending in ‘-ance’ is used in this paper to note that textiles are a mixture of fibers and air, and not pure substances. Because radiation properties such as reflectance and transmittance of non-conducting materials such as textiles do not vary for incidence angles less than 70° from normal [10], the values measured according to the schematic shown in Fig. 1 are reported as near-normal properties.

While not directly measured, the spectral absorptance (α_{λ}) of each fabric can be determined because the reflectance, transmittance and absorptance must sum to unity (ρ_{λ} + τ_{λ} + α_{λ} = 1). Also, according to Kirchoff’s law, the spectral near-normal emittance (ε_{λ}) must equal the near-normal absorptance [10].

#### 3.1 Mid-IR reflectance, transmittance and emittance

The spectral mid-IR properties of the textile samples are shown in Fig. 2. The most significant difference between the various fabrics appears in wavenumbers greater than approximately 1700 cm^{−1} (wavelengths below 6 µm). Each of the fabric samples exhibit similar overall spectral patterns, with a series of narrow absorption peaks between 3000 – 2000 cm^{−1} wavenumbers. Increasing amounts of ceramic content, however, generally decreases the reflectance and transmittance, and increases absorptance, in the 4000 – 2000 cm^{−1} region. Similar behavior occurs at smaller wavenumbers; however, the effect is less pronounced because the overall absorptance is already high at these longer wavelengths.

To better illustrate how the optical properties align with infrared radiation emitted by the human body, Fig. 3 plots these properties as a function of wavelength, as opposed to wavenumber, for the fabric with the maximum amount of added ceramic (1.18 wt. %) and with no added ceramic. The figure also illustrates the spectral distribution of intensity of radiation emitted by a blackbody at a nominal skin temperature of 35°C. This emission is described by Planck’s distribution:

where the first and second radiation constants are C_{1}= 3.742 × 10

^{8}W·µm

^{4}/m

^{2}and C

_{2}= 1.439 × 10

^{4}µm

^{4}·K, respectively [10].

As can readily be seen in Fig. 3, the spectral region where the impact of added ceramic is greatest (λ < 6 µm) corresponds to an area with less significant magnitude of radiation at nominal body temperature; less than 5% of the total radiation of a blackbody at 35°C is contained in this region [10]. To provide a representative value for an optical property, it must be spectrally weighted by its source of radiation. For example, the average reflectance over a wavelength range of λ_{1}-λ_{2} is:

When the broader 2.5 – 16.7 µm mid-IR spectral range is considered, the impact of added ceramic content on emittance is even greater, as evidenced by the steeper gradient for this wavelength range in Fig. 4, due to the more significant difference in properties where λ < 6 µm. Table 1 summarizes the spectral average mid-IR properties of the four fabric samples with varying wt. % added ceramic content. It should be noted that the wt. % values reported in our prior study [3] did not account for the inorganic additive present in standard polyester; hence the values reported in Table 1 vary slightly from what was previously reported for the same samples. However, these slight variations do not affect the conclusions from that study, which remain unchanged.

#### 3.2 Near-IR reflectance, transmittance and emittance

At its relatively modest temperature, the human body does not emit significant radiation at near-IR wavelengths; however, optical properties in this region are of interest due to interactions with radiation from hotter sources (most notably, the sun). The reflectance and transmittance of the fabric samples in the near-IR region are shown in Fig. 5.

The fabric with no added ceramic exhibits significantly higher reflectance and transmittance, and therefore lower absorptance, than fabrics that do contain fibers impregnated with ceramic particles. The magnitude of the impact on ceramic content in the near-IR region is even more significant than the impact in the mid-IR region shown in Fig. 2; this will have significant impact on the solar radiation absorbed by and transmitted through the fabric, as approximately 50% of the solar spectrum falls in the near-IR region [11].

## 4. Impact of optical properties on infrared irradiation received by body

#### 4.1 Heat transfer model

To understand the impact of the change in radiative optical properties on the wearer of a fabric, an energy-conservation-based heat transfer model was developed. This model accounts for all possible modes of radiative interaction – reflection, absorption, and transmission – between the body, fabric and environment, as well as non-radiative heat transfer via conduction and convection. The fabric is modeled as being in instantaneous thermal equilibrium with its surroundings. The ambient surroundings are approximated as being sufficiently large relative to the fabric surface such that radiation is exchanged with the fabric as if the surroundings behave like a blackbody at *T _{amb}*. The fabric is assumed to be sufficiently thin such that its temperature is uniform and cross-fabric thermal gradients can be neglected. Evaporative cooling is neglected in the model, as the goal of the model is to compare the relative radiative heat transfer rates for fabrics with different ceramic contents, and evaporative cooling is expected to be relatively independent of this effect [12]. Under these assumptions, the skin-fabric-environment system and pertinent mechanisms of heat exchange are shown in Fig. 6.

Each radiative heat flux component, shown in red, includes reflected and transmitted radiation in addition to emitted radiation (i.e., it represents a radiosity [10]). Radiative energy is transferred from the skin to the fabric (*J _{sf}*), the fabric to the skin (

*J*), the fabric to the ambient surroundings (

_{fs}*J*), and from the ambient surroundings to the fabric (

_{fa}*J*). The radiative heat flux from the surroundings can include both solar irradiation (

_{af}*J*) and irradiation from the surroundings (

_{solar}*J*) at ambient temperature (

_{env}*T*). The non-radiative heat fluxes due to convection and/or conduction (

_{amb}*Q*and

_{i}*Q*) are shown in blue. The fabric of thickness

_{o}*t*is located between the wearer (at temperature

*T*) and the surroundings (at T

_{skin}*) and is separated from the wearer’s skin by a gap of*

_{amb}*b*thickness.

At steady state, the energy conservation on the fabric can be represented, on a per unit area basis [W/m^{2}] by:

The general form of the non-radiative heat fluxes is *Q = h _{eff} ·*(

*T*) where

_{surface}– T_{surroundings}*h*is the average convective heat transfer coefficient and is a function of the geometry and flow characteristics. Conduction across the air gap for

_{eff}*Q*is a special case of

_{i}*h*, where

_{eff}= k_{air}/b*k*is the thermal conductivity of air, that occurs when the gap

_{air}*b*is sufficiently small such that viscous forces damp out buoyancy effects [10,11]. It is important to note that in the case of

*T*>

_{amb}*T*, the direction of heat transfer will change; however, this is accounted for appropriately in the model (i.e.,

_{fab}*Q*and

_{o}*Q*will be negative, indicating heat flow in the opposite direction of the arrow shown in Fig. 6). The total radiative fluxes (represented by uppercase

_{i}*J’*s in Fig. 6) are determined by integrating the spectral flux, or the radiative heat flux leaving surface

*x*in the direction of surface

*y*is

*J*= ∫

_{xy}*j*(λ)·dλ. Similarly, the solar irradiation is determined by integrating the spectral solar intensity normal to the fabric surface, or

_{xy}*J*= ∫

_{solar}*j*(λ)·cos(

_{solar}*θ*)·dλ. The spectral solar intensity, accounting for atmospheric transmission losses, is taken from ASTM G-173 [13].

_{s}The spectral radiative heat flux leaving each surface has a component of emission, reflection, and transmission and can be written as:

Equations (8)-(10) are circularly defined with dependence upon one another. They are made independent by substituting Eq. (10) into Eq. (9) and solving explicitly for *j _{sf}* (λ):

An iterative solution procedure can be used to solve for the required fabric temperature to satisfy the energy balance given in Eq. (6). The specific heat transfer correlations used to determine the effective heat transfer coefficients in determining the non-radiative heat transfer components *Q _{i}* and

*Q*, as well as the procedure for adjusting the epidermis skin temperature based on the external heat transfer rate, can be found in the Appendix.

_{o}#### 4.2 Spectral shift in incident infrared radiation to the body

Of interest for therapeutic applications is the magnitude and spectral distribution of radiation received by the body as a result of wearing a fabric, namely *j _{fs}* (λ) given by Eq. (10). A common set of baseline input parameters, provided in Table 2, are used in the thermal model. The fabric gap, ambient temperature, air velocity and solar incidence angle can all be adjusted in the model to evaluate sensitivity to these environmental and fit factors.

Under the baseline conditions of Table 2, Fig. 7 shows the resulting spectral distribution of infrared radiation incident on the skin *j _{fs}* (λ) for the fabrics with the maximum (1.18 wt. %) and minimum (0 wt. %) added ceramic content.

In the near-IR spectrum, the body will receive more infrared when wearing a garment composed of fabric with no added ceramic, due to the higher transmittance in this spectral region relative to the fabric with added ceramic. In contrast, the fabric with 1.18 wt. % added ceramic is able to absorb more of this solar near-IR radiation which is then re-emitted by the fabric at longer wavelengths. This shift in spectral incident radiation can be seen in the mid-IR portion of Fig. 7, where there is an increase in *j _{fs}* for wavelengths greater than 4 µm.

To determine if the increase in incident mid-IR radiation persists over a variety of environmental and fit factors, the model was run at different skin-fabric gap distances, ambient temperatures, wind speeds, and incident solar angles. The results of these simulations, used to calculate the total incident radiation in the mid-IR spectrum from 2.5 – 16.7 µm by integrating over this wavelength region for both the 0 wt. % and 1.18 wt. % added ceramic fabrics, are reported in Table 3. For the baseline case, all model input parameters match those listed in Table 2. For the other cases, only the parameter noted is changed from the baseline; the other parameters are held at the baseline state.

Table 3 illustrates that the spectral shift in incident radiation to the wearer of a ceramic-embedded garment to the mid-IR region is present under a variety of scenarios. The incidence angle of solar irradiation appears to have the most significant effect on the magnitude of the shift; a similar effect to a less-normal incidence angle would be expected for cloudy days where the intensity of solar irradiation is reduced.

Increased solar absorption and subsequent re-emission of energy at longer wavelengths may be beneficial in numerous textile applications. The impacts of infrared radiation consisting of wavelengths ranging from 3 – 12 µm on tissue oxygenation and cell stimulation are under investigation [6]. Additionally, performance thermal outerwear designed for cold-weather applications could benefit from additional absorption from the solar spectrum to warm the outer surface of the garment and provide this heat to the wearer of the garment. Experimental studies of such ceramic-modified garments incorporating the use of thermal manikins or similar techniques should be conducted to assess these potential applications.

In summary, this study employed spectrophotometric measurements of textile fabrics modified with ceramic particles, finding favorable agreement with prior measurements using a scientific grade infrared camera on such fabrics. The increase in absorptance and emittance was mainly found to be focused in wavelengths below 6 µm, below the typical wavelength of human body radiation. However, in applications where the fabric receives radiation in this near-IR spectral region (e.g., via sunlight), thermal modeling indicates that the wearer will receive increased infrared radiation at wavelengths in the mid-IR (2.5 – 16.7 µm) portion of the electromagnetic spectrum, encompassing the wavelengths of human body radiation.

## Appendix

#### A.1 Mathematical relationship to determine reflectance through the indirect method

Equation (1) shows the mathematical relationship between the spectral reflectivity with (ρ′_{λ}) and without (ρ_{λ}) a reflective backing (cf. Fig. 1). This relationship can be derived by considering the additional radiation reflected by the sample when a highly reflective backing is present, as shown in Fig. 8.

From Fig. 8, it can readily be seen that:

Because the reflectivity must be less than 1, each term inside the brackets moving rightward becomes increasingly smaller, and Eq. (12) can be rewritten as:

#### A.2 Heat transfer correlations for non-radiative thermal transport coefficients

The general form of the convection heat flux (units, W/m^{2}) from a surface to its surrounding medium is given by Newton’s law of cooling:

*Re*) and/or Grashof number (

_{L}*Gr*) are first calculated to understand the characteristics of the heat and momentum transfer, and along with the fluid’s Prandtl number (

_{L}*Pr*) these values are used to calculate an average dimensionless heat transfer coefficient (

*Nu*=

*f*(

*Re*)). The dimensional average heat transfer coefficient $\overline{h}$ is then the Nusselt number multiplied by the characteristic length (

_{L},.Gr_{L}, Pr*L*) divided by the fluid thermal conductivity$\left(\overline{h}=\overline{Nu}\cdot {L}_{c}/{k}_{fluid}\right)$ [10,11].

_{c}In certain geometries, such as the gap between two parallel plates at different temperatures (i.e., between the skin and fabric), Eq. (14) can be employed regardless of whether conduction or convection is the dominant mode of heat transfer. If dimensionless analysis indicates that viscous dissipation forces outweigh the buoyancy forces generated by local differences in fluid density, conduction will prevail over convection and the Nusselt number will become unity (when $Nu=1=\overline{h}\cdot {L}_{c}/{k}_{fluid}$, it follows that $\overline{h}={k}_{fluid}/{L}_{c}$, where the characteristic length *L _{c}* is the width of the gap between the plates). It can readily be seen that in this scenario, Eq. (14) becomes $Q=\overline{h}\cdot \Delta T={k}_{fluid}\cdot \Delta T/{L}_{c}$ which is simply Fourier’s law for steady-state one-dimensional conduction through a solid of uniform thermal conductivity.

### A.2.1 Outer convection heat transfer coefficient

Convective heat transfer from the fabric to the ambient surroundings can occur by forced convection (resulting from external flow, such as wind, with a “free stream” velocity *u _{∞}*), free convection (driven by local buoyancy differences induced by a temperature gradient), or a combination of both. To determine which of these (or both) is significant; the Reynolds and Grashof numbers must be calculated [10,11].

*L*is the characteristic length (in this case, the width of the plate

_{c}*w*in Table 2), μ is the dynamic viscosity of the fluid, g is acceleration due to gravity, β is the coefficient of thermal expansion (reciprocal of absolute temperature for an ideal gas), and ν is the kinematic viscosity (ν = μ/ρ). All fluid properties are evaluated at the “film temperature,” or the average of the fabric and ambient temperature. Three scenarios are possible and must be considered in estimating the average convective heat transfer coefficient:

- 1. The wind velocity
*u*is sufficiently high for forced convection to dominate natural (free) convection. This is the case when G_{∞}*r*⁄_{L}*Re*_{L}^{2}≪ 1. If G*r*⁄_{L}*Re*_{L}^{2}< 0.1, the average Nusselt number should be calculated based on forced convection correlations. If the Reynolds number is less than the critical Reynolds number of*Re*= 5.5 × 10_{crit}^{5}, the flow is laminar and Eq. (17) applies. If not, the flow is turbulent and Eq. (18) is used [11]. - 2. In cases where the wind velocity is low, natural convection may dominate forced convection (G
*r*⁄_{L}*Re*_{L}^{2}≫1). If G*r*⁄_{L}*Re*_{L}^{2}> 10, the dimensionless Rayleigh number (*Ra*) must be calculated and used to determine the average Nusselt number from natural convection correlations. The Rayleigh number is first used to calculate both the laminar and turbulent Nusselt numbers [11]._{L}= Gr_{L}·Prwhere

*C*and_{lam}*C*are dimensionless parameters given by:_{turb,V}The average Nusselt number can then be determined by asymptotically averaging the laminar and turbulent Nusselt numbers through the following empirical formula:

- 3. If the Grashof number is of the same order of magnitude as the square of the Reynolds number (0.1< G
*r*⁄_{L}*Re*_{L}^{2}< 10), the average Nusselt number can be determined by asymptotically averaging the natural and forced convection Nusselt numbers from Eqs. (17)–(23) through the following empirical formula [11]:Once the average Nusselt number is determined, the outer convection heat transfer coefficient can be determined [10,11].

### A.2.2 Convection heat transfer coefficient between skin and fabric

Similar dimensional analysis can be used to determine the inner convection heat transfer coefficient. It is assumed that the fabric is impermeable and, as such, no bulk flow exists in the space between the fabric and the skin. For the case of two vertically-oriented parallel plates, the convection coefficient will again depend on the Rayleigh number. In this case, the characteristic length is the gap between the parallel plates (*b*), as opposed to the plate width [10,11]:

The average Nusselt number and heat transfer coefficient will depend on the magnitude of the *Ra _{b}* compared to a critical Rayleigh number (

*Ra*

_{cri}_{t}= 1708). If

*Ra*≤

_{b}*Ra*, the viscous dissipation forces of the fluid in the gap will outweigh buoyancy-driven forces. As a result, convection cannot be established and conduction is the dominant heat transfer mode across the gap (

_{crit}*Nu*=1). For

*Ra*>

_{b}*Ra*, the Nusselt number will be greater than unity due to the presence of convective flows, but the appropriate correlation depends on the magnitude of

_{crit}*Ra*. The overall relations for various levels of

_{b}*Ra*can be summarized as [11]:

_{b}Once the average Nusselt number is determined from Eqs. (27)–(29), the average inner convection coefficient between the skin and the fabric can be determined [10,11]:

### A.2.3 Determination of skin temperature

The human skin consists of the inner dermis and the outer epidermis, which are roughly 2 mm and 0.15 mm thick, respectively. The epidermis thermal conductivity (*k _{epidermis}*) is 0.21 W⁄m∙K [14]. It is assumed that the dermis is maintained at the consistent temperature of 35°C due to vasoconstriction and vasodilation. However, it is necessary to account for the temperature drop across the epidermis where there is no blood vessel or capillary thermoregulation (Fig. 9).

The epidermis is modeled as a solid material experiencing one-dimensional, steady-state conduction. As such, the heat flow from the dermis must match the net heat leaving the skin towards the fabric (*Q _{loss,net}* =

*Q*+

_{i}*J*–

_{sf}*J*). Fourier’s law of conduction then can be used to determine the temperature drop across the epidermis:

_{fs}## Acknowledgments

Partial funding for this work was provided by Hologenix, LLC (Santa Monica, CA). We thank Prof. Andrei Fedorov (Georgia Institute of Technology, Atlanta, GA) for insightful discussions regarding the heat transfer modeling of radiative energy exchange.

## References and links

**1. **S. B. Warner, *Fiber Science* (Prentice Hall, 1995) pp. 213–229.

**2. **US Patent, No: 5344297 (1994).

**3. **M. A. Pooley, D. M. Anderson, H. W. Beckham, and J. F. Brennan III, “Engineered emissivity of textile fabrics by the inclusion of ceramic particles,” Opt. Express **24**(10), 10556–10564 (2016). [CrossRef] [PubMed]

**4. **W. W. Carr, D. S. Sarma, M. R. Johnson, B. T. Do, V. A. Williamson, and W. A. Perkins, “Infrared absorption studies of fabrics,” Text. Res. J. **67**(10), 725–738 (1997). [CrossRef]

**5. **P. C. Hsu, A. Y. Song, P. B. Catrysse, C. Liu, Y. Peng, J. Xie, S. Fan, and Y. Cui, “Radiative human body cooling by nanoporous polyethylene textile,” Science **353**(6303), 1019–1023 (2016). [CrossRef] [PubMed]

**6. **F. Vatansever and M. R. Hamblin, “Far infrared radiation (FIR): its biological effects and medical applications,” Photonics Lasers Med. **4**(4), 255–266 (2012). [PubMed]

**7. **J. D. Hardy and E. F. Dubois, “Regulation of heat loss from the human body,” Proc. Natl. Acad. Sci. U.S.A. **23**(12), 624–631 (1937). [CrossRef] [PubMed]

**8. **C.-E. A. Winslow, A. P. Gagge, and L. P. Harrison, “The influence of air movement upon heat losses from the clothed human body,” Am. J. Physiol. **127**(3), 508–515 (1939).

**9. **ASTM E408–13, *Standard Test Methods for Total Normal Emittance of Surfaces Using Inspection-Meter Techniques* (ASTM International, 2013).

**10. **F. P. Incropera and D. P. Dewitt, *Fundamentals of Heat and Mass Transfer, 5th Ed.* (John Wiley & Sons, 2002).

**11. **G. Nellis and S. Klein, *Heat Transfer* (Cambridge University Press, 1997).

**12. **E. A. Arens and H. Zhang, “The skin’s role in human thermoregulation and comfort” in *Thermal and Moisture Transport in Fibrous Materials*, N. Pan and P. Gibson, Eds. (Woodhead Publishing Limited, 2006).

**13. **ASTM G173–03, *Spectral Solar Irradiance* (ASTM International, 2003).

**14. **A. C. Guyton and J. E. Hall, *Textbook of Medical Physiology* (W.B. Saunders Company, 2000).