The Study of Footwear Thermal Insulation Using Thermography and the Finite Volume Method

The paper concerned the issue of ergonomics and thermal comfort of footwear. The presented work was an attempt to apply selected modeling CAD software and the finite volume method for analysis of heat transfer through footwear between the shoes’ user and his surroundings. Based on the actual two selected sport shoes, three-dimensional models have been designed, taking into account their different geometry and raw material composition. The models used simplifications related to the internal construction (in microscale) of the specific parts of the footwear (sole, insole, shank, lining and tongue). The main purpose of the work was to determine the influence of the simplifications used on the accuracy of simulation as a tool to predict thermal insulation of real shoes determined by means of thermal imaging camera. The simulations results were in correlation with the experiment and showed that applied software can be an effective tool in studying the thermophysiological properties of footwear.


Introduction
Thermal comfort of footwear is one of relevant properties that are taken into account during qualitative evaluation of its functionality and suitability. It results from many various factors related to the activity of the human body. Thermal comfort depends on many factors such as: temperature and humidity inside the footwear, its geometry and thermal properties of the raw materials from which it is made. The thermal properties of the raw material determine the efficiency of heat and water exchange between the user and his environment, thus affecting the microclimate of the user's foot, which determines the overall comfort of use of the footwear.
Over the decades, increase in production of plastic and synthetic footwear has been observed. The upper parts are mostly made of leather-like materials, while the soles are made from rubber and synthetic materials. Footwear made of synthetic materials does not provide the right comfort during use. Application of artificial leather instead of natural leather (providing a proper thermal balance between the user and its environment), despite enormous progress in the chemical industry, is still unsatisfactory. Usually, the shank is made of genuine leather or synthetic material leather. Soles are usually made from thermoplastic synthetic polymers and from rubber [1]. Currently, the footwear design process is based on advanced computer applications. They enable precise construction of three-dimensional models of shoes as well as conducting virtual researches on footwear functional properties such as mechanics [2][3][4], thermal comfort [5][6][7][8][9][10], and sensory comfort [11][12][13][14][15][16][17]. Covill et al. [5] developed two-dimensional and three-dimensional thermal models of in-shoe climate were using finite element analysis. The thermal properties of the upper shoe, sole, and air were considered. Dry heat flux from the foot was calculated based on typical blood flow in the arteries on the foot. By means of the thermal models developed, in-shoe temperatures were predicted to cover various locations for various controlled ambient temperatures. Kuklane et al. [6] validated predictions with the model proposed by Lotens (including factors such as blood flow to feet, thermal insulation of footwear and environmental climatic conditions) using actual measurements on subjects exposed to cold environments. The aim of the paper was to check how well the insulation values that were measured on the thermal foot model fit into the prediction model. In other work [7], influence of moisture inside the footwear on the thermal insulation was analyzed. Five types of footwear with various insulation levels were tested. Covill et al. [9] used finite element analysis to describe the heat transfer in footwear. Studies were conducted to determine the temperature distribution in footwear with a variety of environmental temperature and footwear properties considered. Finite element models describing the heat transfer between the foot, sock and shoe were presented with the conductivity, specific heat and density properties of each material presented taken from the literature. On the basis of foot geometry, these models predicted temperatures within footwear with conduction, convection and radiation. One of the most commonly used research methods for evaluation of footwear thermal insulation is the thermal model of the human foot. Endrusick et al. [17] investigated the effects of a simulated environmental exposure on the thermal insulation of an improved version of a military boot designed for use in cold and wet environments. A thermal foot model measured total boot and selected part (toes, heels) thermal resistance values of the standard (non-removable insulation) and the new version of the US Army Intermediate Cold Wet Boot. Kuklane [18] presented work focused on local thermal comfort of feet with special attention to the effects of physical changes on footwear thermal properties. Presented studies have shown that insulation values measured on a thermal foot model correlated with foot skin temperatures measured on subjects. Bergquist and Holmér [19] tested 5 different types of cold-protective footwear with regard to their resistance to dry heat loss (i.e., the insulation) with a new electrically heated foot model. The model was able to simulate 'walking' movements in order to provide a more realistic simulation of wear conditions.
The purpose of the work is related to the important problem of thermal comfort [20], and it is a continuation of work on modeling physical phenomena related to the ergonomics of clothing [21][22][23][24][25]. In this work, studies over problems of footwear thermal insulation were presented. The work presents the results of the attempt to use the selected CAD software to simulate thermal insulation of footwear. The main intention of the work was to analyze of heat transfer occurring in footwear during its use and check the usability of the tool in the case of simple models taking into account the shoe geometry and the basic physical parameters of raw materials affecting thermal insulation of analyzed shoes like density, thickness, thermal conductivity coefficient, specific heat. Two selected sport shoes with a different geometry and composition were tested for their thermal insulation. The studies were divided into two parts: experiment and simulation. In experimental part, thermal insulation of actual shoes was measured using a thermal imaging camera. In the simulation part, 3-D models of both shoes were designed and their thermal insulation in selected specific constant ambient conditions was calculated by the finite volume method.

Materials
The subject of this research was two selected sport shoes ( Fig. 1): Shoe A made in China and Shoe B made in Vietnam. The shoes differ in both geometry and raw material composition presented in with a different geometric structure and raw composition ( Table 1).

Experiment
To assess the thermal insulation, shoes were placed in a room with a constant air temperature: 20°C. During measurement, a rubber container with water at a temperature of 35°C (similar to human foot temperature) was placed inside the shoe. The applica- tion of a flexible container ensured equal direct contact with all the inside walls of the footwear, and it provided a conduction heat transfer through shoe walls to the outside. The temperature of the water in the container was controlled by means of a thermometer. After reaching the steady state, the average temperature of outer surface of three selected parts of shoe (sole, shank, tongue) was measured using a thermal imaging camera (FLIR SC5000 model made in USA) and Altair-Thermal Image Analysis Software. Additionally, temperature distribution on the outer surface of the shoe was determined. To determine the temperature distribution, the entire surface of the shoe, the tested shoe (after reaching the steady state), was placed in several orientations relative to the thermal imaging camera. Because different parts of both shoes were made of different raw materials, it would be impossible to assume different values in the settings of the thermal imaging camera during a given measurement. Because the emissivity values did not differ significantly, the averaged emissivity value from all raw materials of tested shoe was applied. The errors of the temperature measurement resulting from thermal imaging camera specifications were ± 1°C. The presented results were the average of three independent experiments. Scheme of the experiment is presented in Fig. 2.

Model Design
Two 3-D models of both shoes (presented in Fig. 3) were developed in SolidWorks 2014 software based on geometry determined on the basis of images of actual shoes. Both models had the following geometry parameters of the actual shoes: (1) shape of the sole (thickness, profiles, length and width), (2) shape of shank (thickness, profiles), (3) shape of the insole (thickness, profiles length and width), (4) shape of the tongue (thickness, profiles), and (5) shape of the lining (thickness, profiles). The thickness of the individual parts of the shoes has been determined on the basis of their images and application of the ImageJ software (average values are shown in Table 1). The initial phase of the design was to create two-dimensional sketches mapping of complicated geometry of the mentioned parts of shoes (for example, the shape of the soles in the orthogonal projection) using Non-Uniform Rational B-Spline (NURBS) curves. As a result of the applied operations, two-dimensional sketches were transformed into zero-thickness surface objects. The final stage of the designing was to transform these surfaces into three-dimensional objects. Selected stages of designing footwear models are presented in Fig. 4. In both models, simplifications regarding the geometry (no holes for laces, no laces, no decorating elements on outer surface of the shank and no treads on the lower surface of the sole) and concerning internal structure (shank, tongue, lining and insole were mapped as a continuous, homogenous objects) of real shoes were applied. In particular, the simplification concerning internal construction obviously had to affect the thermal insulation properties of footwear. All the parts of footwear mapped in the model (sole, insole, lining, shank, tongue) have, in fact, a complex, layered internal structure with inhomogeneous porosity. The air (characterized by a much lower coefficient of thermal conductivity compared to typical raw materials for the production of footwear) filling space between layers and pores increases the thermal insulation of the footwear.

Physical Basis of Heat Flow Simulation
The simulations were performed using Solidworks Flow Simulation 2014 software applying computational fluid dynamics (CFD) analysis [26] to enable simulation of heat transfer by the finite volume method.
The software allows predicting simultaneous heat transfer in solid and fluid media with energy exchange between them. Solidworks Flow Simulation solves energy conservation equations and Navier-Stokes formulas describing fluid flow. The software calculates temperature fields created by: (1) heat transfer in solids (conduction), (2) free, forced, and mixed convection, and (3) radiation in both the steady state and transient state [27]. Solidworks Flow Simulation uses the Navier-Stokes equations, which are formulations using momentum and energy conservation and mass laws for fluid flows. Additionally, the software applies fluid state equations defining the nature of the fluid by empirical dependencies of fluid viscosity, thermal conductivity, and density on temperature [27]. The Navier-Stokes equations are supplemented by definitions of thermophysical properties and state equations for the fluids. Solidworks Flow Simulation allows to gas and liquid flow modeling with density, viscosity, thermal conductivity, specific heats, and species diffusivities as functions of pressure, temperature, and species concentrations in fluid mixtures.
Materials used in the production of footwear often have a complicated internal structure. For example, the porous structure makes the thermal properties of footwear depend not only on the properties of applied raw material but also on the air inside the pores. Then heat is transported through both (solid body) and fluid (gas), with simultaneous exchange between these environments. Heat transfer in fluids is described by means of the following conservation equation: where S i S where The constant C μ is determined according to [27] as equal to C μ 0.09, whereas σ c 0.9. The equations describe both laminar and turbulent flows. Moreover, transitions from one case to another and back are possible, k is turbulent kinetic energy, while ζ , J·kg −1 ·s −1 is turbulence dissipation (rate at which turbulence kinetic energy is converted into thermal internal energy). The parameters k and μ t are zero for pure laminar flows. The anisotropic heat conductivity in solid media is described by the following equation: where e cT . It is assumed that the heat conductivity tensor is diagonal to the considered coordinate system and that the heat transport within analyzed material is direction independent, i.e., we introduce an isotropic medium and can denote λ 1 λ 2 λ 3 λ. The energy exchange between the fluid and solid media is computed via the heat flux in the direction normal to the solid/fluid interface, taking into account the solid surface temperature and the fluid boundary layer characteristics, if necessary [27].

Conditions of Simulations
The aim of simulations was to determine the temperature of the outer surface of three selected parts of the shoes (sole, shank, tongue) using calculation on designed models.  The initial conditions of simulation (presented in Fig. 5) imitated real conditions during use of the footwear in four selected ambient temperature -10°C, 0°C, 10°C, and 20°C, pressure 1013.25 hPa, while assumed temperature inside of the shoes was equal to human foot temperature: 35°C. The raw materials with physical properties presented in Table 2, have been assigned to the tested models.
During calculations, the analyzed shoe model was placed in the center of computational domain (presented in Fig. 6

) divided into cells of three types (fluid cells, solid cells, and partial cells).
Fluid cells were filled with air, solid cells filled with shoe, while the partial cells contained both air and shoe. The number of cells (presented in Table 3) was dependent on the number and size of the elements making up the model of particular shoe.

Results and Comments
The images obtained by means of thermal imaging camera (presented in Fig. 7) show temperature distribution on outer surface of both shoes in ambient temperature: 20°C after a period of about 15 min when the steady state was reached. The area marked in pink color indicates a temperature of 35°C of container filled with a water placed inside the shoes. The obtained images showed that in case of both shoes the highest temperature (red color) was observed in the middle part of the shank and the lower part of the tongue, where temperatures ranged from 33°C to 28°C. The simulation and experiment results obtained for ambient normal conditions show the heterogeneity of the temperature distribution on all three analyzed parts of the footwear. In addition, areas with the same temperature on the corresponding parts of the shoe differ in shape and temperature value. As can be seen in the experimental and simulated temperature distributions, the area with the highest temperature was the middle part of the shank, while the coolest parts were: the sole, toe cap, quarter, collar and upper part of tongue. Differences in shape of the corresponding temperature fields between experiment and simulation resulted probably from the heterogeneous internal structure of parts of real shoes (sole, insole, shank, tongue), while the both shoe models establish their uniform internal structure. In Table 4, temperature values (minimum, average, maximum) obtained using simulated and experimental distributions on analyzed parts of both shoes determined in normal ambient conditions were presented. Experimental average temperature T av of the given part of shoe (sole, shank, tongue) was calculated as a weighted average temperature of individual pixels of thermogram included in the analyzed part. Simulated average temperature was calculated on the basis of simulated distributions in the same way like experimental. The extreme values, on the other hand, were points with highest and lowest temperatures for the analyzed part of the shoe, which were determined on the basis of the experimental and simulated distributions. In addition, the percentage differences between simulated and experimental temperature values are shown in Table 4.
In Fig. 10, average temperature on outer surface of three selected parts of tested shoes in four ambient temperatures computed by simulations and determined in the experiment was presented.  Simulation results revealed that the average temperatures of three selected parts of shoes increased in higher ambient temperatures. For both tested shoes, the largest temperature differences determined for different ambient conditions occurred for soles, which probably resulted from their much bigger thickness in comparison with thickness of shanks and tongues. However, these differences were clearly bigger for sole of Shoe B, which could be the result of their bigger thickness and lower heat conduction coefficient of raw material (polyurethane), and the sole of Shoe A was made of synthetic rubber. In addition, for shank and tongue in Shoe B temperature differences were larger than for Shoe A. Although these parts in both shoes were made of the same raw materials, they clearly differed in thickness.
Outcomes of thermal imaging conducted for the ambient temperature 20°C confirmed that in the two analyzed shoes a sole was the most thermo-insulating part, while a tongue was the least one. The results of the experiment showed that average temperature of three analyzed parts was always lower than the calculated one. The percentage differences between simulated and experimental values of average temperature of the given part of footwear (both obtained in normal climate) were as follows: T sole 3.6 %, T shank 4.7 %, T tongue 2.9 % for Shoe A and T sole 3.5 %, T shank 3.3 %, T tongue 2.2 % for Shoe B. However, the differences in average values are not sufficient to fully assess the usefulness of the designed models. The comparison of extreme temperature values (Table 4) showed that the differences between the results of the experiment and modeling are greater. The biggest differences between the results of the simulation and the experiment were observed in the case of minimum temperatures (16.5 %-21.9 % in the case of Shoe A and 17.3 %-23.6 % in the case of Shoe B). Significantly smaller differences were observed in the case of maximum temperatures (3.7 %-9.6 % in the case of shoe A and 3.9 %-5.9 % in the case of shoe B). Small differences in the case of simulation and experimental values of average temperatures and much clearer differences in the case of extreme temperatures can be a decisive influence of the complex internal structure of the analyzed footwear which should be taken into account when analyzing the evaluation of its thermal insulation. For both analyzed shoes, the greatest differences in extreme temperatures were found in the shank and the tongue, i.e., parts with a much higher porosity than the sole. The main reason for this could be the fact that both parts are characterized by higher porosity, more complicated geometry, and a more heterogeneous internal structure than the sole.
Despite applied simplifications, simulation results obtained at normal climate correlated with the experiment and showed influence of geometry and raw material composition of the three analyzed parts of both shoes on their thermal insulation.

Conclusions
In the work, the issues related to the thermal comfort of footwear were studied. An attempt to use the selected CAD software to design three-dimensional models of two selected actual sports shoes (differing in geometry and raw material composition) in order to determine their thermal insulation properties using the finite volume method was presented. Designed models mapped both the geometry and the raw material composition of the real footwear. In first applied approach for both shoe models, significant simplifications in mapping of internal structure of analyzed parts (shank, tongue, lining and insole were modeled as a continuous, homogenous, non-porous objects) were applied. The main goal of this approach was to attempt to assess the impact of the simplifications used on the accuracy of simulation as a tool to predict thermal insulation of real shoes determined by means of thermal imaging camera. Obtained simulation results led to the following conclusions: 1. The simulations outcomes were in correlation with the experiment; however, in the case of both shoes, the thermal insulation calculated by means of simulation was slightly lower than that determined in the experiment. 2. Not taking into account in both models the complex internal structure of real shoes, in particular heterogeneity, porosity and thus the free spaces filled by air decreased simulated thermal insulation. 3. Obtained differences between experimental and the simulation temperature distributions for footwear models indicate the need to correct these models, probably in the direction of increasing its similarity to the structure parameters of actual shoes as: a. heterogeneous internal structure b. heterogeneous porosity 4. Obtained differences between experimental and the simulation results may also be caused by the method of conducting experimental research (application of a common average value of emissivity for all raw materials forming shoes) 5. Despite the above-mentioned simplifications of applied software and developed geometrical footwear models allowed to predict the average temperatures of analyzed parts of real shoes with the error on the level of 2.9 % − 4.7 % (Shoe A) and 2.2 % − 3.5 % (Shoe B) in relation to the verification experiment. 6. The presented models could be useful in predicting the thermal comfort of the designed and manufactured footwear containing raw materials like natural leather, synthetic rubber, polyurethane and polyamide.