In this study, an algebraic model and its experimental verification was carried out to investigate the effect of moisture content on the heat loss that takes place due to conduction of sock fabrics. The results show that increasing moisture content in the studied socks caused a significant increase in their conductive heat loss. Plain knitted socks with different fiber composition were wetted to a saturated level, and then their moisture content was reduced stepwise. When achieving the required moisture content, the socks samples were characterized by the Alambeta testing instrument for heat transfer. Three different existing modified mathematical models for the thermal conductivity of wet fabrics were used for predicting thermal resistance of socks under wet conditions. The results from both ways are in very good agreement for all the socks at a 95% confidence level. In the above-mentioned models, the prediction of thermal resistance presents newly a combined effect of the real filling coefficient and thermal conductivity of the so-called “wet” polymers instead of dry polymers. With these modifications, the used models predicted the thermal resistance at different moisture levels. Predicted thermal resistance is converted into heat transfer (due to conduction) with a significantly high coefficient of correlation.

#### Keywords

- Heat transfer
- mathematical modeling
- plain socks
- moisture content
- filling coefficient
- volume porosity

Heat transfer by conduction mostly depends on fiber's thermal conductivity, volume, orientation to the heat flow direction, and fabric construction [1, 2]. Thermal conductivity of the fabrics can be evaluated in different ways: experimental method, analytical solution method, and numerical method [3]. Most of the studies on thermal conductivity/thermal resistance in the wet state to date are experimental and reported that there was an increase in thermal conduction consequent to increasing the moisture content [4,5,6,7,8]. This study will provide a quantitative prediction of the insulation loss with the addition of water in socks. Sock insulation is linked to the material thickness and trapped air (porosity) in between the fibers. The increase in the insulation of footwear may not be obvious for the shoes only. It depends on the sock's insulation as well because low insulation shoes can get more insulation comparatively from socks than shoes. Footwear insulation is a vital feature for feet warmth; on the other hand, the sweat generation due to the motion of subjects can strongly affect the foot temperatures. Mostly the cold feet sensation is associated with low skin temperatures due to sweating [9]. Even well-insulated footwear will start feeling cold on wetting. Socks are made of fabrics where the absorbed moisture can strongly influence their thermal comfort properties, since a human foot could generate up to 30–50 g of sweat per hour in a hot environment [10, 11]. At a high physical activity, it could be 30 g · h^{−1} even in the cold environment [12]. For a single foot, the most recent study reports the range (10.3 ± 3.6 g · h^{−1}) with shoes compared to an uncovered foot (12.6 ± 3.7 g · h^{−1}) [13]. Owing to these high sweat rates, the heat loss may substantially increase. Prolonged damp and cold conditions can cause injuries similar to a trench foot. The trench foot, however, does not require a freezing temperature; it can occur at hot temperatures as well [14]. Many studies on thermal resistance prediction through empirical models are available in the literature; these models are specifically volume fractions and based on their respective thermal conductivities. Numerical approaches can deal with uneven profiles, solid/liquid/gas phases, different forms of heat transfer, number of boundary conditions, and uneven material properties. Numerical methods also have the potential to attain the utmost precision [15]. However, these methods are intrinsically more complex and awkward; further, under certain conditions, plain methods demonstrated to be more precise for much less stab [16]. Some researchers employed artificial neural networks (ANN) models for thermal resistance [17, 18] and thermal conductivity [19] predictions. In most of the studies, thermal resistance is predicted by statistical models [20]. Some researchers have predicted the thermal conductivity of wet fabrics with mathematical approaches. Schuhmeister [21] developed a relationship to calculate the thermal conductivity of the mixture of air and fiber with the following assumptions:

Fibers are distributed homogeneously in all directions;

One-third of fibers were placed in parallel; and

Two-thirds of fibers were placed in series or perpendicular to the heat flow.

Later, many researchers followed the footprints of Schuhmeister by changing the ratio of fibers placed in series to those placed in parallel [22, 23]. Recently, Militky considered placing 50% fibers in series and 50% in parallel to the heat flow [24]. R. S. Hollies and Herman Bogaty have suggested a parallel combination for measuring the effective thermal conductivity of moistening fabric by combining the volume fraction and thermal conductivity of water and polymer [25]. Mangat presented a number of mathematical models for thermal resistance (wet state) in the series and in parallel combinations of air, fiber, and water resistance. His predictions are in good correlation with the experiments [26, 27]. The problem with Mangat's models is the assumption made concerning the filling coefficient or porosity as constant components. However, varying the moisture levels can result in changes because water has a different density. Although his second assumption that the air is replaced by water is theoretically correct, but he did not quantify it. Hes, R. S. Hollies, and Herman Bogaty have ignored the series arrangement in their suggested models. It will predict the lower thermal resistance as heat will conduct along with the thickness of the fabric. S. Naka et al. suggested a theoretical approach for thermal conductivity prediction but they did not use it for calculations. They also involved the warp and weft fabric thickness in their suggested model. The model used by Dias and Delkumburewatte is a three parameters series model, which is a very simple approach; yet, they ignored the parallel conduction part so it can predict higher thermal resistance. Presently, Mansoor et al. [28, 29] have modified Maxwell–Eucken- (ME-)2, Schuhmeister, and Militky models by combining the water and fiber filling coefficient for predicting the thermal resistance of wet socks. We have used the same approach for the prediction of thermal resistance and heat transfer through conduction in this study. The general aims of this study are as follows:

To predict the heat transfer of wet socks through conduction by using the Mansoor et al. modified models.

To study conductive heat transfer (predicted/experimental) in the extended state (controlled moisture content in %) for simulating a real extension and minimizing the effect of the dimensional changes.

To ascertain the volume porosity of socks with and without extension in the model.

To ascertain the volume porosity and pore size distribution of socks by X-ray micro-tomography (MCT) scanning without extension.

By combining the fiber and water filling coefficients approach, below three models have predicted the reasonable thermal resistance for socks that are in agreement with the experimental results. These models are discussed in the forthcoming sections.

Maxwell introduced the two-phase concept for the determination of electrical conductivity [30]. Later on, Eucken used the same analogy for the thermal conductivity evaluation [31]. Brailsford and Major (Eq. [1]) have modified the ME models for thermal conductivity of a three-phase mixture by assuming first phase as continuous while other two as dispersed [32].

Later on Eq. (1) was generalized by Wang et al. [33] as shown by Eq. (2).

ME's Eq. (3) is obtained by assuming air and wet polymer as disperse and continuous phases, respectively, for Eq. (2). ME model (Eq. 3) can be used to describe an effective thermal conductivity of a two-component material with simple physical structures. Equation (3) represents a two-component system for effective thermal conductivity, based on volume fraction and respective. Many effective thermal conductivity models require the naming of continuous and dispersed phases. In materials with exterior porosity, individual solid particles are surrounded by a gaseous matrix, and hence the gaseous component forms the continuous phase and the solid component forms the dispersed phase [34]. For external porosity, they are considered as continuous and dispersed phases, respectively, as follows.

F_{wet polymer} and λ_{wet polymer} are calculated as per Eqs (13)–(15).

Schuhmeister's equation, given by Eq. (4), summarized the relationship between the thermal conductivity of fabric and the fabric structural parameters by an empirical equation [21].
_{fab} is the thermal conductivity of fabric, F_{wet polymer} is the conductivity of wet fibers, λ_{a} is the conductivity of air, F_{wet polymer} is the filling coefficient of the solid fiber, and F_{a} is the filling coefficient of air in the insulation.

Militky's equation, given by Eq. (7), summarized the relationship between the thermal conductivity of fabric by an empirical equation [24].
_{s} and λ_{p} are calculated using Eqs (5) and (6), respectively, λ_{fab} is the thermal conductivity of fabric, λ_{wet polymer} is the conductivity of wet fibers, λ_{a} is the conductivity of air, F_{wet polymer} is the filling coefficient of the solid fiber, and F_{a} is the filling coefficient of air in the insulation.

All the plain (single jersey) socks samples as shown in (Table 1 and Figures 1 and 2) were knitted on the same machine (Lonati Goal GL544S, 144Needles, Diameter 4”, 4Feed) settings by varying the main yarns to get the homogeneous samples with respect to specs and stretches for contrast comparison. “The yarn running at the surface of the sock is called the main yarn and the plaiting yarn (generally spandex covered polyamide or polyester filament yarn) runs inside the fabric providing stretch, elasticity, comfort and shape to the sock” [35, 36]. After knitting, all the samples were processed for washing in the same machine bath followed by tumble drying and boarding.

Sock samples specifications.

^{−2}) | |||||
---|---|---|---|---|---|

29.525/1 tex 100% Cotton spun yarn | 2.22/8.33/36/1 tex Polyester air covered Elastane (91:9)% | Cotton 80%, Polyester 18.20%, Elastane 1.8% | 129.88 | 0.95 | P1 |

29.525/1 tex 100% Viscose Spun yarn | Viscose 81.08%, Polyester 17.22%, Elastane 1.70% | 130.44 | 0.90 | P2 | |

29.525/1 tex 100% Spun Polyester | Polyester 98.38%, Elastane 1.62% | 125.70 | 0.95 | P3 |

Images for all the samples are taken using a camera (ProgRes-CT3) attached to a macro-scope (NAVITAR) with the lights and stand (HAISER Germany) from front and backside of the socks under extension (Figure 2). Main and plaited yarns could be observed from front side and back side of the socks easily. Main yarn is more prominent from front side, except that some of the plaited yarn loops could be seen over the main yarn for samples P1, P2, and P3 due to contrast of color. Plaiting yarn along with the main yarn could be easily distinguished from the back side of these samples.

Sock's structure is important due to several advantages. Physically, it presents properties of comfort such as high elasticity, conformity with the shape of the body, softer hands feel, and others. In general, heat and mass transmission rate is dependent mainly on the fabric's geometrical parameters, namely, thickness and porosity [37]. Porosity (ɛ) is the volumetric ratio of the pores accessible by total volume [38]. The porosity of the fabrics can be calculated by air permeability, image processing, and geometrical modeling approaches [39]. Volume porosity of the socks was determined according to Eq. (8) [40, 41].
_{0} is fiber density (kg · m^{−3}) and ρ is fabric density (kg · m^{−3}).

3D porosity of the socks was investigated by using an x-ray computed MCT SKYSCAN 1272 system. In this system, radiation is converted into an electrical signal between the x-ray source and the detector when the specimen revolves on a vertical axis. 2D images are taken in several steps during this rotary motion. Reconstruction software generates a 3D model of the actual specimen from these images [42]. The following are the common settings for all the tested samples: image pixel size, −3.0 μm; lower grey threshold, −33; upper grey threshold, −255; rotation step, −0.2°; rotation degrees, −180°; frame averaging, −3; exposure, −672 ms; voltage source, −50 kV; source current, −200 uA.

For the extension simulation, the socks were loaded on a dummy leg (Salzmann MST Switzerland) [43] of medium size (24 cm) as per specification of the standard method (RAL-GZ-387/1). Then worn socks are marked as per the testing template. After unloading, the socks were extended to the marked circle with the help of an embroidery hoop as shown in Figure 3. Sock samples were wet to the saturated level (100% moisture content) by BS EN ISO 105-X12 standard test method. The established technique for preparing a wet fabric of the known oven-dry fabric weight is used and thus the fabric is thoroughly wetted in distilled water. The wet pick-up is brought to 100 ± 0.5% by putting wet testing fabric on a blotting paper. The evaporation of the moisture content below the specified level was avoided by using polyethylene bags. Further, we tested again for the above-mentioned tests under extension at different moisture levels.

A testing method was selected for this research as per the condition of worn socks and limitations of the manikins. Socks worn inside shoes show 1st-order boundary conditions and constant, different temperatures on both surfaces of the fabric (such as Alambeta). Socks were worn (calf area) partly under 3rd-order boundary conditions; conduction inside = convection outside (thermal foot model, permetest). The condition is more clearly illustrated in Figure 4. Furthermore, short testing time (almost keeping the specific moisture content) distinguish the Alambeta from other skin models and manikins. So, Alambeta was selected especially for wet testing.

The thermal resistance (Rct) of the developed samples was measured by Alambeta tester [44], which rapidly measures both steady-state and transient-state thermal properties. This instrument simulates the heat flow q (W · m^{−2}) from the human skin to the fabric during a short initial contact in the absence of body movement and external wind flow. When the two bodies are brought into ideal contact, the time course of temperatures is determined by solving the known one-dimensional partial differential equation Eq. (9) according to the relation subsisting between them.

The measuring head drops down, touches the fabrics, the heat flow levels are processed, and the thermo-physical properties of the measured specimens are evaluated [5]. The measurement lasts for several minutes only. Thus, reliable measurements on wet fabrics are possible, since the sample moisture during the measurement is almost constant. As mentioned earlier, socks are worn inside the shoes under 1st-order boundary conditions, and Alambeta-testing corresponded well with the use of socks inside a shoe (boundary conditions of first-order). Thermal resistance assessed using the Alambeta tester, which enables fast measurement of both steady-state and transient-state thermal properties, are shown in (Figure 5). This diagram demonstrated the maximum, q_{max}; dynamic (transient), q_{dyn}; and heat flow, q_{steady}.

Steady-state heat flow has been shown by Eq. (11).

Thermal resistance (m^{2} · K · W^{−1}) is used to express the heat insulation properties of the fabric. Thermal resistance of textiles is affected by fiber conductivity, fabric porosity, and fabric structure. It is also a function of fabric thickness, as shown by Eq. (12).

An average thermal conductivity for different fibers (within socks) at different moisture levels will be calculated as per Eq. (13).
_{w} = water filling coefficient; F_{fib1} = 1st-fiber filling coefficient; F_{fib2} = 2nd-fiber filling coefficient; l_{w} = water thermal conductivity; I_{fib1} = 1st-fiber thermal conductivity; and I_{fib2} = 2nd-fiber thermal conductivity. Filling coefficients for water, fiber, wet polymer, and the air are calculated as per steps given in Table 2.

Filling coefficients.

_{w} = water filling coefficient | _{fib} = fiber filling coefficient | |
---|---|---|

Moisture content | % | % |

Mass | gram | gram |

Area | m^{2} | m^{2} |

Areal density | ||

Volumetric density | ||

Filling coefficient |

Air filling coefficient (F_{a}) is calculated as per Eq. (14).

Filling coefficient for wet polymer will be calculated as per Eq. (15). This value will be used as input in all above models for measurement of thermal resistance in wet states.

Thermal conductivity of water and air are taken as 0.600 and 0.026 (W · m^{−1} · K^{−1}), respectively, while the density of water is 1,000 (kg · m^{−3}). Different values were found for the thermal conductivity of textile fibers. However, the following values of density [45] and thermal conductivity have been taken for different fibers in this study; these are given in Table 3.

Different fibers’ properties.

^{−3}) | ^{−1} · K^{−1}) | |
---|---|---|

Cotton | 1,540 | 0.50 [46] |

Viscose | 1,530 | 0.50 [22, 45, 46] |

Polyester | 1,360 | 0.40 [45] |

Validation of the theoretical models is done by comparing results (x) with results obtained by experiments (y) for a set of parallel determinations. If both methods (theoretical and experimental) lead to same results, then the dependence of y on x is linear ((y = β_{1}x + β_{2}) with zero intercept β_{2} = 0 and unit slope β_{1} = 1. This validation is done by the joint confidence region for intercept and slope because estimators are correlated. Assumptions for this composite inference will be :

Null hypothesis H_{0}: β_{2} = 0 and β_{1} = 1

Alternative hypothesis H_{1}: β_{2} ≠ 0 and β_{1} ≠ 1.

Level of significance: ∝ = 0.05

Test statistics:

Critical region:

Check the value from table for Fisher–Snedecor F-distribution F_{0.95} (m, n–m).

Conclusion:

If the calculated value as per Eq. (16) is less than the critical value, then accept the null hypothesis H_{0}: β_{2} = 0 and β_{1} = 1. It means intercept and slope are not significantly different from 0 and 1, respectively, at a 95% confidence level. A simultaneous test of the composite hypothesis confirmed that a new laboratory method (by theoretical model) is in agreement with the results of a standard one (experimental). Further, if the calculated value is higher than the critical, then alternative hypothesis H_{1}: b_{2} ¹ 0 and b_{1} ¹ 1 will be accepted with the conclusion that theoretical model results are not in agreement with the experimental results [47].

Volume porosity (with and without extension) was determined by a semi-empirical approach as per Eq. (8) and fabric areal density and thickness need to be determined experimentally. Figure 6 shows the final (from the top, wales, and course views) images of all the tested socks scanned by MCT scanner (SkyScan 1272) as 2D and converted into 3D using NRecon. A sample size of 5 × 5 mm has been used for scanning these images.

For porosity quantification, distribution of the pores, and pore thickness, the above images were analyzed by using software recommended by the manufacturer (BRUKER), which is CTAn. The color coded images (Figure 7) were generated by CTVox by using the data provided by CTAn. The measurement of the 3D pore thickness is referred to as “sphere-fitting” and this thickness is considered as the diameter of the largest enclosed sphere [48].

The results of the volume porosity demonstrated that extended socks have higher porosity (Figure 8). This increase in the porosity was also reported by Abdolmaleki et al. at different extension levels for loose knitted fabrics [49]. Porosity falls in the range of 78% and 90% without and with extension, respectively. Guidoin et al. stated that knitted fabrics’ porosity lies between 67% and 84% and even 90% is not uncommon [40]. Extension causes an increase in the pore size (space between loops) of the fabric and decreases the fabric thickness. It leads to a decrease in the volume of the fiber (solid part) and increases the volume of air corresponding to porosity. Porosity measured by MCT is in agreement with theoretical porosity (without extension) at a 95% confidence level for all the socks. Thus, the thermal resistance model's prediction in the next sections is based on this porosity model. This comparison is logical and it further validated that the used model for the calculation of porosity is correct. The difference is in the range of (0.1400–4.3715%) for all the socks except P1. The value of 7.4256% lower porosity is measured by MCT with respect to the predicted value, which is close to the difference observed by Doczyova et al. – i.e., 6% during porosity comparison of knitted structures [50].

Figures 9, 11, and 13 clearly demonstrate that as the moisture (%) increases, heat transfer also increases, irrespective of sock fiber composition. For all the models, the input thermal conductivity and filling coefficients were measured in wet polymer at different moisture levels. The correlation between experimental and predicted models was checked by coefficient of determination (R^{2}). The values of coefficient of determination (Figures 10, 12, and 14) for all the three modified models (ME-2, Schuhmeister, and Militky) showed that these models could make reasonable predictions of heat transfer in the dry as well as the wet condition at different moisture levels for all the major fiber blends used for socks. Coefficient of determination falls between 0.78 and 0.90 for all the samples.

Following are the assumptions made for the development of theoretical models for the prediction of heat transfer in the wet state:

Fabric thickness assumed as constant.

No Free convection (as Rayleigh Number < 1,000).

The difference between the temperatures on both surfaces of the fabric 1st-order boundary conditions is considered constant.

To simplify the model, fiber filling coefficient is assumed as constant.

Air and water filling coefficients are variable.

Fiber (polymer) and water filling coefficients are combined as wet polymer filling coefficient.

Values for thermal conductivity of wet polymer in relation to water and fibers are combined as per their volume.

No dimensional changes occurred at different moisture levels as tested in extended state.

Fabric areal density and thickness are measured in the extended state.

Alambeta's thickness is considered.

The predicted and experimental heat transfer by conduction of P1 (cotton 80%, polyester 18.20%, and elastane 1.8%) at various moisture levels is given in Figure 9. All the three, namely, Maxwell-modified, Militky-modified, and Schuhmeister models have the best prediction at different moisture levels for the P1 sample. Militky-modified, Schuhmeister-modified, and ME-2–modified models have R^{2} values of 0.8776, 0.8766, and 0.8586, respectively, as shown in Figure 10. The heat loss is increasing with the increase of moisture level (Figure 9). About 50% increase in the heat transfer is observed at 30% moisture content. This increase is uniform till 20% moisture content, and then a rapid rise is observed between 20% and 30% moisture content. This trend is evident in Figure 9 (green square legends). Overall, Schuhmeister has the highest prediction due to 67% consideration of thermal conduction in series followed by the Militky-modified model. This means that as the portion of series in consideration decreases, thermal conductivity increases and further the result is in line with previous investigations of fiber alignment in series that has 2–3 times lower thermal conductivity than parallel [23, 51]. The findings are in accordance with the findings of Wang et al. [33]. They have predicted the thermal conductivity with respect to porosity by using different combinations and models, i.e., ME-1, ME-2, series, parallel, EMT, series + parallel, ME-1 + ME-2, etc. According to their findings, the series model has the lowest thermal conductivity, followed by ME-2, while the parallel model has the highest prediction. Reddy and Karthikeyan [52] also have the same findings during their study for predicting the thermal conductivity of frozen and unfrozen food materials.

Validation of the theoretical models is done by comparing the results (x) with results obtained by experiments (y) for a set of parallel determinations. This validation is done by joint confidence region. On substitution into Eq. (16):

In the case of ME-2–modified model, (F_{1} = 0.5174) is lesser than the quantile of the Fisher–Snedecor F-distribution F_{0.95} (2, 3) = 9.5521, so the null hypothesis H_{0} cannot be rejected. It means the predicted heat transfer with the ME-2–modified model is not significantly different from the experimental results. Similarly, on substitution into Eq. (16) for Militky-modified model, we obtain:

In Militky-modified model, (F_{1}= 5.8238) is lesser than the quantile of the Fisher–Snedecor F-distribution F_{0.95} (2, 3) = 9.5521, so the null hypothesis H_{0} cannot be rejected. Again, it is verified that the predicted heat transfer with the Militky-modified model is not significantly different from the experimental results. On substitution into Eq. (16), for composite confidence region validation with Schuhmeister-modified model, we obtain:

In Schuhmeister-modified model, (F_{1} = 18.5481) is greater than the quantile of the Fisher–Snedecor F-distribution F_{0.95} (2, 3) = 9.5521, so the null hypothesis H_{0} cannot be accepted. It means that the predicted heat transfer with the Schuhmeistermodified model is significantly different from the experimental results.

In the case of P2 sock (viscose 81.08%, polyester 17.22%, and elastane 1.77%), Militky-modified model has the best prediction at 11.45% and 19.50% moisture levels, as shown in (Figure 11). ME-2 modified has a better heat transfer prediction at 30.30, 40.17%, and 49.80% moisture levels. All three models have a reasonable prediction with R^{2} > 0.93, as shown in (Figure 12). Similar to the P1 sample, a rapid increase in the heat transfer with the increased moisture content is also observed, between 20% and 30% moisture content. This increase is in agreement with Naka and Kamata's study and close to the value reported by Mangat, i.e., 70% [27, 53]. Schuhmeister-modified model has the highest prediction followed by Militky-modified and ME-2–modified models at all the moisture levels. Over and over again lowest to the highest prediction of heat transfer order by different models has verified the findings of Finck [51], Bogaty et al. [23], Wang et al. [33], and Reddy [52]. From these studies, it has been established that series alignment has predicted the lowest thermal conductivity followed by ME-2; combinations of ME-2, ME-1, EMT, series, and parallel; EMT; ME-1; and parallel.

The values of composite confidence region for slope and intercept at 95% confidence level validated all the theoretical models except Schuhmeister-modified model, as calculated value F_{1} = 35.5029 is higher than the critical value F_{0.95} (2, 3) = 9.5521. This means that the heat transfer prediction with Schuhmeister-modified model is not significantly correct with respect to experimental results. As per acceptance of the null hypothesis H_{0}, ME-2–modified model is validated as having lower F_{1}, i.e., 0.7598 than the critical value 9.5521. In case of Militky-modified model, (F_{1} = 7.3972) is lower than the quantile of the Fisher–Snedecor F-distribution F_{0.95} (2, 3) = 9.5521, so the null hypothesis H_{0} cannot be rejected. On substitution in to Eq. (16):

F_{1} for Militky-modified model:

F_{1} for Schuhmeister-modified model:

Figure 13 depicts theoretical and experimental thermal resistances of P3 socks (polyester 98.38% and elastane 1.62%) at various moisture levels. Militky-modified, Schuhmeister-modified, and ME-2–modified models have R^{2} values 0.8425, 0.841, and 0.7671 respectively (Figure 14). The increase in the heat transfer is slower and uniform between 5–10% and 20–50% moisture content levels. However, this rise (42%) is fast between 10% and 20% moisture content, as evident from experimental green square legends (Figure 13). Unlike P1 and P2, 50% of the heat transfer increase in P3 is observed at 50% moisture content due to the hydrophobic nature of polyester. Once more Schuhmeister-modified model has a higher prediction at all the moisture levels except 5% and 10% moisture content. It has predicted 0.5–2 times lower thermal conductivity. This is in accord with Mao and Russel's study [54]. They have observed 0.5–3 times lower thermal conductivity prediction for 100% polyester spacer fabric with Schuhmeister's model. They have not incorporated moisture content. Even then their predictions are very high with respect to experiments.

As per the values of the constructed confidence region for slope and intercept at 95% confidence level, all the theoretical models except Schuhmeister-modified model have been validated. It means the intercepts (β_{2}) and slopes (β_{1}) are not significantly different from 0 and 1, respectively. So, the heat transfer prediction with ME-2 and Militky-modified models is not significantly different with respect to experimental results for sample P3. Calculated values of F_{1} also justify that the ME-2–modified model has top prediction among all others followed by Militky-modified and Schuhmeister-modified models. To summarize, ME-2–modified model has the better forecast for sample P3 than both other models. Test statistics (calculated F_{1}) values are 0.6071, 0.8773, and 10.9260 for ME-2–, Militky-, and Schuhmeister-modified models, respectively, against the critical value of the Fisher–Snedecor F-distribution F_{0.95} (2, 4) = 6.9443. F_{1} is calculated for ME-2 on substitution in to Eq. (16):

F_{1} for Militky modified:

F_{1} for Schuhmeister modified:

A semi-empirical approach was used to model the heat transfer prediction of plain knitted socks in the wet state. This work focuses on the heat transfer prediction of socks in the wet state. The effects of fiber dimensions, orientation, and yarn hairiness were not taken into account. Both theoretical porosity and experimental for socks were calculated, but the heat transfer prediction is based on theoretical results. Image processing and analysis were used to obtain data on the porous structure of fabrics. Validation of the models has been done through the coefficient of determination (R^{2}) and inference statistics, i.e., hypothesizing slope = 1 and intercept = 0 at 95% confidence interval.

By adopting this new approach of feeding the wet polymer filling coefficient and the thermal conductivity instead of dry polymers, different models can provide a justified prediction of heat transfer under wet conditions as well. All the models have a reasonable coefficient of determination, i.e., R^{2} range is in between 0.76 and ~0.95 for all the sock samples at different moisture levels. However, as far as the validation through hypothesis, i.e., slope = 1 and intercept = 0 is concerned; Schuhmeister-modified model could not qualify for all the socks. The polymer filling coefficient remains constant while water and air filling coefficients are changing with the variation of moisture that leads to changing the thermal conductivity. A higher value of moisture causes an increase in the heat transfer. A 50% increase in heat transfer was observed at 30% moisture content in all the samples, except P3 (polyester) socks. We hope that this study has uncovered many worthy avenues for future investigations. The inquisitive readers will no doubt have ideas of their own, but the following are some suggestions for research of possible interest:

This study was conducted by assuming thickness and GSM as constant. A separate study could be planned to identify the effect of swelling on the thickness, especially in hydrophilic fabrics.

Future studies could be planned for examining other types of fabrics and mathematical models by adopting this approach.

Shoes could be included with the addition of more boundary conditions.

#### Different fibers’ properties.

^{−3}) | ^{−1} · K^{−1}) | |
---|---|---|

Cotton | 1,540 | 0.50 [46] |

Viscose | 1,530 | 0.50 [22, 45, 46] |

Polyester | 1,360 | 0.40 [45] |

#### Sock samples specifications.

^{−2}) | |||||
---|---|---|---|---|---|

29.525/1 tex 100% Cotton spun yarn | 2.22/8.33/36/1 tex Polyester air covered Elastane (91:9)% | Cotton 80%, Polyester 18.20%, Elastane 1.8% | 129.88 | 0.95 | P1 |

29.525/1 tex 100% Viscose Spun yarn | Viscose 81.08%, Polyester 17.22%, Elastane 1.70% | 130.44 | 0.90 | P2 | |

29.525/1 tex 100% Spun Polyester | Polyester 98.38%, Elastane 1.62% | 125.70 | 0.95 | P3 |

#### Filling coefficients.

_{w} = water filling coefficient | _{fib} = fiber filling coefficient | |
---|---|---|

Moisture content | % | % |

Mass | gram | gram |

Area | m^{2} | m^{2} |

Areal density | ||

Volumetric density | ||

Filling coefficient |

Apparel Industry in the EU–China Exports and Circular Economy Automatic Identification Of Wrist Position In A Virtual Environment For Garment Design Pressure Evaluation Of Seamless Yoga Leggings Designed With Partition Structure Experimental and Modelling Studies on Thermal Insulation and Sound Absorption Properties of Cross-Laid Nonwoven Fabrics Tensile Properties Analysis Of 3D Flat-Knitted Inlay Fabric Reinforced Composites Using Acoustic Emission Optimization of Sodium Lignosulfonate Treatment on Nylon Fabric Using Box–Behnken Response Surface Design for UV Protection A Study on the Woven Construction of Fabric Dyed With Natural Indigo Dye and Finishing for Applying to Product Design for Home Textile Products A Calculation Method for the Deformation Behavior of Warp-Knitted Fabric Nondestructive Test Technology Research for Yarn Linear Density Unevenness Numerical Simulation and Analysis of Airflow in the Condensing Zone of Compact Spinning with Lattice Apron Blend Electrospinning of Poly(Ɛ-Caprolactone) and Poly(Ethylene Glycol-400) Nanofibers Loaded with Ibuprofen as a Potential Drug Delivery System for Wound Dressings Application of Plasticized Cellulose Triacetate Membranes for Recovery and Separation of Cerium(III) and Lanthanum(III) Analysing Service Quality and its Relation to Customer Satisfaction and Loyalty in Sportswear Retail Market A Review on the Performance and Comfort of Stab Protection Armor A Fabric-Based Integrated Sensor Glove System Recognizing Hand Gesture Developing Real Avatars for the Apparel Industry and Analysing Fabric Draping in the Virtual Domain Simulations of Heat Transfer through Multilayer Protective Clothing Exposed to Flame Determination of Sewing Thread Consumption for 602, 605, and 607 Cover Stitches Using Geometrical and Multi-Linear Regression Models Evaluation of Functional Insoles for Protective Footwear Under Simulated Use Conditions Designing a Three-Dimensional Woven Fabric Structure as an Element of a Baby Stroller Computer-Assisted Modeling and Design of Compression Garments with Graded Unit Compression Application of Physical Vapor Deposition in Textile Industry Modeling Lean and Six Sigma Integration using Deep Learning: Applied to a Clothing Company Comparative Analysis of Structure and Properties of Stereoscopic Cocoon and Flat Cocoon Effect of Water pH on Domestic Machine Washing Performance of Delicate Textiles Effect of Different Yarn Combinations on Auxetic Properties of Plied Yarns Analysis of Heat Transfer through a Protective Clothing Package Smart Textile for Building and Living Investigation of Twist Waves Distribution along Structurally Nonuniform Yarn 3D Body Scan as Anthropometric Tool for Individualized Prosthetic Socks Preliminary Experimental Investigation of Cut-Resistant Materials: A Biomimetic Perspective Durable Wash-Resistant Antimicrobial Treatment of Knitted Fabrics Modeling Supply Chain Sustainability-Related Risks and Vulnerability: Insights from the Textile Sector of Pakistan Numerical Simulation of Fiber Motion in the Condensing Zone of Lateral Compact Spinning with Pneumatic Groove Study on the Thermal and Impact Resistance Properties of Micro PA66/PU Synergistically Reinforced Multi-Layered Biaxial Weft Knitted Fabric Composites Improvement of Physical Properties of Viscose Using Nano GeO _{2}as Doping MaterialFea-Based Structural Heat Transfer Characteristic of 3-D Orthogonal Woven Composite Subjected to the Non-Uniform Heat Load Bending Failure Behavior of the Glass Fiber Reinforced Composite I-Beams Formed by a Novel Bending Pultrusion Processing Technique Comfort-Related Properies of Cotton Seersucker Fabrics Economical and Social Dimensions of Unionization in Turkish Textile and Clothing Sector Conductive Heat Transfer Prediction of Plain Socks in Wet State A Novel Foam Coating Approach to Produce Abrasive Structures on Textiles Textronic Solutions Used for Premature Babies: A Review Effect of Lycra Weight Percent and Loop Length on Thermo-physiological Properties of Elastic Single Jersey Knitted Fabric Texture Representation and Application of Colored Spun Fabric Using Uniform Three-Structure Descriptor Analysis of Mechanical Behavior of Different Needle Tip Shapes During Puncture of Carbon Fiber Fabric Approach to Performance Rating of Retroreflective Textile Material Considering Production Technology and Reflector Size Influence of Multilayer Interlocked Fabrics Structure on their Thermal Performance Prediction of Standard Time of the Sewing Process using a Support Vector Machine with Particle Swarm Optimization A Novel Theoretical Modeling for Predicting the Sound Absorption of Woven Fabrics Using Modification of Sound Wave Equation and Genetic Algorithm Ag Coated Pa-Based Electro-Conductive Knitted Fabrics for Heat Generation in Compression Supports Design Method of Circular Weft-Knitted Jacquard Fabric Based on Jacquard Module Image Analysis as a Method of the Assessment of Yarn for Making Flat Textile Fabrics Investigation of Heat Transfer in Seersucker Woven Fabrics using Thermographic Method Research into the Textile-Based Signal Lines Made Using Ultrasonic Welding Technology Transformable Warning Clothing for Children with Active Light Sources Regenerated Cellulose/Graphene Composite Fibers with Electroconductive Properties High-Performance Workwear for Coal Miners in Northern China: Design and Performance Evaluation Comfort-Related Properties of Double-Layered Woven Car Seat Fabrics Experimental Investigation of the Wettability of Protective Glove Materials: A Biomimetic Perspective An Integrated Lean Six Sigma Approach to Modeling and Simulation: A Case Study from Clothing SME Mechanical Properties of Composites Reinforced with Technical Embroidery Made of Flax Fibers Consumer Adoption of Fast-Fashion, Differences of Perceptions, and the Role of Motivations Across the Adoption Groups A New Consumer Profile Definition Method Based on Fuzzy Technology and Fuzzy AHP Optimal Design of a Novel Magnetic Twisting Device Based on NSGA-II Algorithm Microscopic Analysis of Activated Sludge in Industrial Textile Wastewater Treatment Plant Evaluation of Physical and Mechanical Properties of Cotton Warps Under a Cyclic Load of Stretch-Abrasion Theoretical and Experimental Evaluation of Thermal Resistance for Compression Bandages Effects of Flocks Doping on the Dynamic Mechanical Properties of Shear Thickening Gel Estimation of Seams in Paraglider Wing Sensitivity of Aerodynamic Characteristics of Paraglider Wing to Properties of Covering Material Numerical Investigation of Heat Transfer in Garment Air Gap Determination of State Variables in Textile Composite with Membrane During Complex Heat and Moisture Transport Design and Performance Evaluation of Protective Clothing for Emergency Rescue Biological Properties of Knitted Fabrics Used in Post-Burn Scar Rehabilitation Fabrication and Characterization of Fibrous Polycaprolactone Blended with Natural Green Tea Extracts Using Dual Solvent Systems