The polar bear and several other Arctic mammals use fur composed of hollow-core fibers to survive in extremely cold environments. Here, we use finite element analysis to elucidate the role that the hollow core plays in regulating thermal transport. Specifically, we establish a three-dimensional model of a textile based on fibers with various core diameters and study transverse heat transport. First, these simulations revealed that textiles based on hollow-core fibers conduct significantly less heat than their solid-core counterparts with fibers with a core-to-fiber diameter ratio of 0.95, reducing thermal transport by 33%. In addition to this decrease in thermal transport, the mass per area of textiles is substantially reduced by making them hollow core. This led us to consider the performance of multi-layer textiles and to find that four-layer hollow-core textiles can exhibit a four-fold decrease in heat flux relative to single-layer solid-core textiles with the same mass per area. Taken together, these simulations show that hollow-core fibers are well suited for thermal insulation applications in which gravimetric thermal insulation is a priority.
I. INTRODUCTION
The demand for weight-efficient thermally insulating materials is particularly acute due to changing climate conditions and persistent energy challenges.1,2 Textiles, in particular, are widely employed for their thermally insulating properties in fields such as performance clothing,3,4 shelters,5–7 medical dressings,8,9 and cold-weather outerwear.10,11 Such textiles are generally evaluated based on their weight, cost, appearance, and thermal insulation.2,12–15 In the drive to improve these metrics, researchers have turned to biomimicry, a field where motifs found in nature are adapted to become synthetic structures.16–30 In the context of thermal insulation, the polar bear is a particularly inspiring natural example, as these animals maintain a body temperature of 36.9 °C with minute fluctuations of ∼0.5 °C even in environments as cold as −40 °C.31–38 A number of factors have been implicated in this feat, including thick layers of subcutaneous fat,32 the large thickness of their fur,33 the ability of this fur to be transparent to thermal radiation, allowing such radiation to be absorbed by the skin,34–36 and the hollow structure of the fibers in their fur.37,38 Among these adaptations, the hollow fibers of their fur stand out as particularly intriguing for textile design. Further evidence of this adaptation’s utility is that it is also seen in other cold-climate mammals such as reindeer, caribou, and moose.39–42 However, the specific degree to which the hollow nature of the fur contributes to its thermal insulation has not been reported, nor have the thermal properties of textiles based on hollow fibers been analyzed in detail. Melt spinning has been used to make hollow polyester fibers,43 and these have been experimentally shown to have improved thermal resistance relative to solid-core fibers in non-woven fabrics.44 However, the main commercial avenue for hollow-core fibers has been filtration due to their high surface area-to-volume ratio.45 Given recent innovations in advanced fiber design and extrusion techniques,46–52 fundamental knowledge and design heuristics are needed to guide the development of advanced textiles.
Here, we analyze the thermal insulation of textiles based on fibers with hollow and solid cores to determine the role of the hollow core in their thermal properties. Initially, a parametric model is constructed to represent textiles with varying core diameters and layers. Next, a series of multiphysics simulations are performed to determine the effective medium thermal conductivity of each textile geometry. A notable result of this study is that textiles based on hollow-core fibers conduct significantly less heat than their solid-core counterparts, with fibers with a core-to-fiber diameter ratio of 0.95 exhibiting a 33% reduction in thermal transport. A more striking comparison emerged from considering multi-layered textiles. For example, multi-layer hollow-core textiles exhibited a four-fold decrease in heat flux relative to single-layer solid-core textiles with the same mass per area. These simulations highlight that hollow-core fibers are well suited for applications in which gravimetric, rather than volumetric, thermal insulation is a priority, as is the case for many clothing applications.
II. METHODS
The unit cell had a cross-sectional area of 4D × 4D, which is as closed-packed as possible, preserving the same fiber density in both lateral dimensions. This unit cell was selected to preserve the oscillatory nature of the weft’s path while maintaining contact between the warp and weft, thus providing a realistic depiction of fibers in a textile.
The modeled textile was used as the basis for a simulation of heat transport [Figs. 1(b) and 1(c)]. The top and bottom boundaries of this region were set to constant temperatures to model thermal diffusion in a cold environment such as the natural habitat of the polar bear. In particular, the upper boundary was held at a constant temperature of 253.15 K to represent free air. The bottom boundary was held at 310.15 K to represent the physiological temperature of mammal skin. Periodic boundary conditions were applied to all other sides of the simulation region to simulate the behavior of a continuous textile. Air was used to fill the interstitial spaces around and inside the fibers. The thermal conductivities of air and the fiber were set to 0.025 and 0.24 mW/mm⋅K (value taken for cotton fiber),53 respectively. We assumed that in this thin boundary layer near the textile, free convection was negligible,54 which was confirmed using multiphysics simulations. The overall study volume was 0.512 mm3, with a length of 0.8 mm, a width of 0.8 mm, and a height of 0.8 mm. A constant fiber diameter of 0.2 mm was selected. Finite element analysis (FEA) simulations of heat transfer were carried out using COMSOL Multiphysics using the steady-state thermal transport module. The unstructured tetrahedral mesh was constructed using the built-in physics controller with a minimum element size of 8 μm and a maximum element size of 64 μm. The properties of the model and the simulation are detailed in Table I.
Parameter . | Value . | Rationale . |
---|---|---|
Textile geometry | Standard weave | Architecture is based on a standard woven pattern attainable using a loom.55 Fiber arrangement emulates the interwoven paths of the warp and weft fibers |
Fiber diameter (D) | 0.2 mm | Chosen in accordance with the ∼0.175 mm diameter of polar bear fur and standard fiber diameters used in woven textiles, ensuring relevance to both biological mimicry and textile manufacturing56 |
Core diameter (d) | 0–0.19 mm | The core diameters were selected to encompass d/D ∼ 0.4 exhibited by polar bear fur while enabling a thorough examination of core geometry57 |
Number of layers | 1–4 | Typical multi-layer woven fabrics are in this range58 |
Thermal conductivities | Air: 0.025 mW/mm⋅K, fiber (cotton): 0.24 mW/mm⋅K | Values were selected to reflect potential applications in thermal insulation, with cotton chosen for its prevalent use in clothing |
Simulation mesh type | Unstructured tetrahedral mesh | Chosen for its flexibility in handling complex geometries |
Grid size (minimum element size) | 8 µm | This size was chosen to be smaller than the smallest geometric feature, which was the 10 μm wall of the hollow fibers with d/D = 0.95 |
Grid size (maximum element size) | 64 µm | Chosen to be smaller than the largest geometric feature, which was the 200 μm diameter of solid-core fibers |
Height of simulation domain | 0.2 mm larger than the textile thickness | The bottom of the textile was assumed to rest against the hot boundary (e.g., skin), while the top of the simulation domain was separated from the top of the textile by 0.2 mm to account for a hydrodynamic boundary layer in air |
Temperature boundary conditions | Top (253.15 K) and bottom (310.15 K): Fixed, sides: Periodic | The top and bottom boundaries were set to constant temperatures to model thermal gradients in Arctic environments, while periodic boundary conditions on the sides simulate the continuous nature of the textile structure |
Parameter . | Value . | Rationale . |
---|---|---|
Textile geometry | Standard weave | Architecture is based on a standard woven pattern attainable using a loom.55 Fiber arrangement emulates the interwoven paths of the warp and weft fibers |
Fiber diameter (D) | 0.2 mm | Chosen in accordance with the ∼0.175 mm diameter of polar bear fur and standard fiber diameters used in woven textiles, ensuring relevance to both biological mimicry and textile manufacturing56 |
Core diameter (d) | 0–0.19 mm | The core diameters were selected to encompass d/D ∼ 0.4 exhibited by polar bear fur while enabling a thorough examination of core geometry57 |
Number of layers | 1–4 | Typical multi-layer woven fabrics are in this range58 |
Thermal conductivities | Air: 0.025 mW/mm⋅K, fiber (cotton): 0.24 mW/mm⋅K | Values were selected to reflect potential applications in thermal insulation, with cotton chosen for its prevalent use in clothing |
Simulation mesh type | Unstructured tetrahedral mesh | Chosen for its flexibility in handling complex geometries |
Grid size (minimum element size) | 8 µm | This size was chosen to be smaller than the smallest geometric feature, which was the 10 μm wall of the hollow fibers with d/D = 0.95 |
Grid size (maximum element size) | 64 µm | Chosen to be smaller than the largest geometric feature, which was the 200 μm diameter of solid-core fibers |
Height of simulation domain | 0.2 mm larger than the textile thickness | The bottom of the textile was assumed to rest against the hot boundary (e.g., skin), while the top of the simulation domain was separated from the top of the textile by 0.2 mm to account for a hydrodynamic boundary layer in air |
Temperature boundary conditions | Top (253.15 K) and bottom (310.15 K): Fixed, sides: Periodic | The top and bottom boundaries were set to constant temperatures to model thermal gradients in Arctic environments, while periodic boundary conditions on the sides simulate the continuous nature of the textile structure |
III. SIMULATION RESULTS AND DISCUSSION
In order to estimate the effect of porosity on the magnitude of thermal conductivity, we simulated the heat transport across textiles composed of fibers with solid cores and textiles composed of fibers with hollow cores [Fig. 2(a)]. Specifically, we found that solid-core fibers exhibited a heat flux Qc = 1.46 mW/mm2 under the simulation conditions, while fibers with d/D = 95% exhibited Qc = 0.98 mW/mm2, showing a ∼33% reduction in heat transport. To justify the use of steady-state thermal simulations, we perform transient thermal transport analysis and find that all structures reach steady-state conditions in <16 ms [Fig. 2(b)]. To explore this process in more depth, we repeated this simulation while sweeping d from 0.02 to 0.19 mm in steps of 0.01 mm. The transition from large heat transport with solid cores to smaller heat transport with hollow cores was found to be monotonic and smooth, with the most dramatic change at d/D → 1 [Fig. 2(c)]. This result makes it clear that hollow fibers result in lower heat transfer, motivating their study for use in harsh weather environments. The reduction in heat transport with increasing porosity can be attributed to the hollow cores interrupting the flow of heat through the fibers.
Effective medium theory can be used to understand and quantify the degree to which being hollow improves fiber performance relative to other possible configurations of materials. Specifically, we may compute an effective-medium thermal conductivity ke that is based on Qc = keΔT/h, where h is the height of the simulation domain and ΔT is the temperature difference across the domain. For this analysis, we employ h = 0.6 mm, so that we only consider heat flow across the region with fibers. Based on this, we find that the solid core fiber exhibited ke = 0.077 mW/mm⋅K, while the textile with d/D = 95% exhibited ke = 0.030 mW/mm⋅K. Therefore, the total thermal transport can be reduced by more than a factor of two by making the fibers hollow cores.
All the observed structures exhibited Ef ∼ 50%, showing that being hollow does not make the architecture intrinsically more insulating per unit mass. In fact, lower Φf, which corresponds to larger hollow cores, is slightly less efficient. Instead, the main reason why these are better insulators is that more of their total volume is made up of air.
To test the hypothesis that textiles consisting of layers of hollow fibers will further reduce thermal transport, we developed computational models of textiles with two, three, and four layers of warp and weft.59 The practicality of multi-layer woven structures is supported by prior work realizing such structures using weaving.60,61 These computations revealed that Qc decreased with increasing textile thickness. For instance, at d/D = 50%, we found Qc = 1.367 mW/mm2 for a single layer and Qc = 0.859 mW/mm2 for a two-layer textile. The 37% percent decrease in Qc while the thickness doubled is a smaller decrease than one would expect if the textile was behaving as a simple effective medium. This is because the air boundary layer above the textile has a constant thickness. This result emphasizes the common principle of conductive thermal insulation, which is that air has the lowest thermal conductivity in the system, so engineering the system to have more entrapped air is a major benefit. To explore how this trend depends on d/D and the number of layers, we repeated these computations with up to four layers with d varied from 0% to 95% of D [Fig. 3(a)]. As expected, larger d or more layers decreased thermal transport.
While it is clear that hollow fibers or multiple layers improve thermal insulation, this is but one consideration when designing textiles. Another key consideration is the mass per unit area. To quantify this, we define the textile mass per unit skin area Mt. It is not a priori clear that trends in Mt will have a simple connection to trends in Qc, as adding layers or decreasing d will affect Mt. Therefore, we compute Mt for each textile and plot Qc vs Mt [Fig. 3(b)]. When visualized together, it is clear that adding layers vastly reduces thermal transport. Perhaps more strikingly, the hollow core fibers are sufficiently light that multiple layers can obtain lower heat transport than a single layer of solid core textile without increasing the mass. In particular, a four-layer fabric with d/D = 0.85 has the same Mt as a single layer textile based on a single fiber while being four times less thermally conductive. These results highlight the opportunity of hollow-core fibers to provide a path to thermal insulation that is highly efficient per unit mass.
In order to explore the role of fiber material, we repeated this analysis by varying fiber thermal conductivity kf. Fiber kf was set to tabulated values for common textile materials.62 The resulting thermal transport increased monotonically with kf [Fig. 3(c)]. For fibers with very small or very large kf, being hollow did not affect thermal transport, as in both limits, the thermal conductivity of air limited heat flow. For intermediate values of kf, hollow fibers exhibited markedly reduced thermal transport. This finding suggests that for fibers with kf ∼ 0.75 mW/mm⋅K, the impact of fibers being hollow is especially high in terms of managing thermal transport.
To objectively compare the thermal insulation performance of all geometries that were considered, we use a figure of merit (FOM) defined using FOM = ΔT/(QC⋅MT). This reflects the thermal insulation provided per unit mass of a given textile and monotonically increases with increasing d/D [Fig. 3(d)]. Interestingly, the number of layers considered has a weak effect on this FOM, as the added layer increases thermal insulation and mass in nearly equal measure. Comparing the FOM computed for textiles based on hollow cotton fibers with tabulated experimental values for a variety of textiles,58,63 it is clear that being hollow can increase this FOM by >10×. This is a particularly exciting prospect when considering performance textiles such as those based on solid-core Kevlar-PBI/Nomex fibers as a starting part for improvement.
IV. CONCLUSIONS
In summary, this work utilized extensive finite element analysis to determine the role that hollow cores have on the thermal transport of woven textiles. First, we found that making the fibers hollow offers a modest reduction in thermal transport while slightly reducing the effective medium efficiency. In other words, the textile did not insulate better than would be expected given that more of the more thermally conductive fiber had been replaced by less thermally conductive air. This observation points to hollow cores being primarily useful as a tool for light-weighting insulating materials. The main quantitative result of this study is that transitioning to hollow-core fibers can reduce thermal transport by 33% and decrease mass density by >90% relative to solid-core fibers. We anticipate that the dual benefits of better thermal insulation and lighter mass will motivate further study into how to efficiently realize hollow-core structures for thermal insulation applications. To emphasize this point, we compute that four-layer textiles based on hollow-core fibers can reduce the overall thermal conductivity by more than a factor of four without adding to the mass per unit area. As the textile architecture studied here can be readily woven using looms and associated techniques, there are no fundamental barriers to experimentally realizing the structures discussed here on a large scale. Given the high importance of insulation normalized by mass in, for example, performance clothing or thermal insulation in aerospace, these results motivate the study of advanced fiber geometries to incorporate structures such as hollow cores.
ACKNOWLEDGMENTS
The authors acknowledge support from the Toyota Research Institute and the Department of Mechanical Engineering at Boston University. The authors acknowledge Dr. Timothy J. Lawton for helpful discussions.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
A.D.A.: Investigation, software, analysis, and writing. K.A.B.: Conceptualization, writing, and supervision.
Adedire D. Adesiji: Formal analysis (equal); Investigation (equal); Software (equal); Writing – original draft (equal); Writing – review & editing (equal). Keith A. Brown: Conceptualization (equal); Supervision (equal); Writing – review & editing (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon request.