1 Introduction
Nowadays, as the increasing fossil fuel continuous depletion [
1-
3] and global severe environmental problems [
4-
6], exploring novel energy-saving and environmentally-friendly materials become essential for the sustainable development of the economy and society [
7-
9]. The burgeoning progress of microscale and nanoscale porous polymer materials in recent years has demonstrated that this kind of materials with its low density [
10-
14], high surface area [
15-
19] and perfect thermal insulation properties [
20-
24] for possible challenging social and economic issues. Generally speaking, different materials, processing methods, and environmental conditions also influence the morphologies of porous polymers. Figure 1 represents the SEM images of various micro/nano-porous polymers’ morphological characteristics [
24-
32]. Meanwhile, the morphology of the porous polymers could be divided into open-cell and closed-cell morphologies. Open-cell porous polymers have highly interconnected [Figs. 1(a)–(g)], while closed-cell porous polymers have isolated pores in which the gas is enclosed [Figs. 1(h) and (i)] [
33,
34].
Actually, porous polymers are widely used in various burgeoning fields, such as gas capture [
35-
38], catalysis [
39-
41], water treatment [
42-
44], sensors [
45-
47], molecular separation [
48-
52], food industry [
53-
55], proton conduction [
56-
58], pharmaceutical industry [
59-
61], drug delivery [
62-
64], supercapacitors [
65-
67], synthetic chemistry [
68-
70], petrochemical engineering [
71-
73], Li-on batteries [
74-
78], environmental protection [
79-
82], and some others [
83-
91].
In recent years, in order to solve the increasingly severe energy shortage problem, the potential of polymer foam in the field of energy utilization and management system has been dramatically highlighted, such as outer layer of buildings [
92-
97], photoenergy conversion in solar cells [
98-
102], electrochemical energy storage [
103-
108], personal thermal management devices [
109-
114], and some others [
24,
115-
119]. Therefore, it is necessary to build an accurate and appropriate thermal conductivity prediction model of micro/nano-porous polymers for better designing and applying its heat transfer characteristics [
120-
123] in energy conversion and storage systems [
124-
128], as shown in Fig. 2 [
79,
122,
123,
129-
133].
This review summarizes the thermal conductivity prediction models of porous polymers, to help subsequent scholars to conveniently carry out their research. Firstly, it reviews the process of establishing the overall equivalent thermal conductivity calculation model for micro/nano porous polymers. Then, the thermal conductivity prediction calculation models are introduced according to thermal conduction and thermal radiation separately. In addition, the thermal conduction part is divided into the gaseous thermal conductivity calculation model, solid thermal conductivity calculation model and gas–solid coupling model. Finally, we conclude and point out the aspects that still need further study of the existing thermal conductivity prediction models of micro/nano porous polymers. Such studies would provide a more accurate prediction of thermal conductivity and a broader application in energy conversion and storage systems.
2 Overview of the thermal conductivity predicted models
In order to analyze the heat transfer characteristic of the porous polymers, it ought to be noted that thermal convection plays a minor role in closed-pore materials with pore sizes below 4 mm in diameter [
22,
33] and in open-pore systems with pore sizes less than 2 mm [
25,
134]. Therefore, in microscale and nanoscale, heat transfer of porous polymers is mainly composed of two parts: thermal conduction and thermal radiation, while thermal convection at such a confined scale is negligible [
22,
25,
135]. That is, the effective total thermal conductivity of the micro/nano-porous polymers
is considered as composed of conductive and radiative components, as [
136-
139]:
where the conductive thermal conductivity , and the radiative thermal conductivity . To sum up, the overall calculation flow chart is shown in Fig. 3. Next, the prediction models of micro/nano-porous polymers would be discussed from thermal conduction and thermal radiation parts.
2.1 Conductive thermal conductivity models
Generally, while establishing the conductive thermal conductivity models, the gaseous and solid thermal conductivities would be calculated separately by selected predicted models. And then bring the results into the gas–solid coupling model to solve the overall conductive thermal conductivity [
14,
22,
25]. Next, these three kinds of predicted models would be introduced separately in the following sections.
2.1.1 Gaseous thermal conductivity models
According to the molecular kinetics theory [
140], the gaseous thermal conduction is carried out by the collision of gas molecules during thermal motion [
141]. The gas molecules move randomly, and the molecules at the high-temperature side have higher speed collide with molecules at the low-temperature side, causing the energy to be transferred to the lower speed molecules [
141,
142].
The gaseous thermal conductivity
inside a micro/nanocellular porous polymer was determined by taking the Knudsen effect [
122,
123,
143,
144] into account, which significantly lead to the gaseous thermal conductivity decrease [
25,
136]. This Knudsen effect implies that when cell size is comparable or smaller than the mean free path of the liquid or gas, the latter molecules collide more often with the molecules forming the surrounding solid part than among them, reducing the energy transfer through the gas molecules [
34,
143]. In a confined space, translation of gas molecules is governed by the Knudsen regime, in which the influence of cell size and mean free path on the efficiency of energy transfer is considered [
123,
136,
140],
where
is the Knudsen number [
145-
147], which refers to the ratio of the mean free path
of the gas to the cell size
d [
136].
β is the energy interaction between gaseous molecules and the solid surface [
25,
148], which is related to the thermal accommodation coefficient
[
149,
150] and heat capacity ratio
γ [
136,
151,
152]. The gaseous thermal conductivities are calculated for different polymer materials considering the Knudsen effect, compared with the experimental data in related reference, shown in Fig. 4 [
25,
122,
140,
153].
For the aerogel materials with large porosity and specific surface area, the gaseous thermal conductivity model proposed by Zeng
et al. [
154] also had a wide range of use in the prediction models of aerogel [
142,
155,
156], which is simplified as
where
p is pressure,
Ss is specific surface area,
T is temperature,
denotes porosity and
ρ represents the density of the porous material. Besides, considering the actual pore distribution in porous materials, Reichenauer
et al. [
153] proposed a calculation model for the gas phase thermal conductivity of porous materials based on double pore size distribution as well as the Gauss distribution. Meanwhile, Bi
et al. [
157] proposed a modified gaseous thermal conductivity model considering the randomness and inhomogeneity of pore distribution in aerogel materials.
What’s more, Fricke
et al. [
158] fitted the gaseous thermal conductivity of aerogel as a function of material density based on the experimental measurement results of aerogel thermal conductivity, which can be written as
It should be noted that the above equation is the contribution of the gaseous conductivity obtained directly from the experimental data to the overall thermal conductivity, not the gaseous thermal conductivity inside the pore. However, this empirical formula also needed to be refitted when the materials, densities, and boundary conditions changed.
2.1.2 Solid thermal conductivity models
Since the characteristic size of the porous polymer is at the microscale or even nanoscale, which is close to or smaller than the mean free path, leading to a noticeable size effect on the heat transfer progress of the solid framework [
140,
141]. As a result, the solid heat transfer of polymer materials is significantly reduced [
141].
In some research, the solid thermal conductivity
is theoretically extracted by considering the absence of gas in the porous polymers [
25,
136], which is directly equal to the thermal conductivity of the polymer skeleton.
Considering the influence of phonon scattering at the gas–solid interface as seen in Fig. 5, Wang
et al. [
122,
123] gave the thermal conductivity of the thin polymer skeleton based on the kinetic theory [
140,
159], as
where is the thermal conductivity of pure solid material, δ is the characteristic size of the polymer skeleton. is the phonon mean free path of the bulk polymer, which can be calculated based on the thermal conductivity of bulk polymer:
where
Cv is the specific heat per unit volume contributed by phonons, and
vg is the mean group velocity of phonons [
140,
159]. Besides, this kind of modified model, considering the size effect, also wildly used the thermal conductivity calculation of aerogel materials [
141,
142,
160]. Therefore, the thermal conductivity of the PS film based on Eq. (6) was calculated with a phonon mean-free path of 1.5 nm, which is in good agreement with those in the literature, as shown in Fig. 5.
Due to the size effect, the polymer film’s thermal conductivity reduces significantly with its thickness in the given range [
140,
159]. What’s more, Fricke
et al. [
157] fitted the solid thermal conductivity of aerogel as a function of material density based on the experimental measurement results of aerogel thermal conductivity, can be written as
It should be noted that this formula also needed to be refitted when the materials, densities, and boundary conditions changed.
2.1.3 Gas–solid coupling models
The existing models for the conductive thermal conductivity of porous polymers generally fall into four categories: (i) Empirical model; (ii) Equivalent circuit method; (iii) Equivalent fractal method; (iv) Numerical simulation method. In the following sections, we will introduce the simplified assumptions and application scope of these four methods.
Empirical models. Fricke
et al. [
157] fitted the thermal conductivity of aerogel as a function of material density based on the experimental measurement results. This method is used to adding the contribution values of thermal conductivity, which is simple in form and convenient to handle, as [
142,
158]
Alshrah
et al. [
161] identified the morphological features in an organic resorcinol-formaldehyde (RF) aerogel. They also correlated each mode of the heat transfer, assuming no coupling effect occurred within the solid–gas phase.
In subsequent studies, some scholars still use this method to calculate aerogel’s overall equivalent thermal conductivity, such as Reichenauer
et al. [
153], Lee
et al. [
162], Deng
et al. [
163], Spagnol
et al. [
164], and Swimm
et al. [
165]. In addition, combined with the effective medium method, this empirical method [
141,
166,
167] is often used to predict the equivalent thermal conductivity of polymer foam composites doped with carbon and other materials. However, this empirical prediction method usually has several empirical constants with no specific physical meaning in its equation, which is often changed with working conditions [
142]. Therefore, this kind of method needs empirical correlation comparing with experimental data each time, which is specifically dependent on the materials’ properties and cellular structure.
Equivalent circuit method. Based on Kirchhoff’s law, equivalent circuit method simplifies the complex structure into a regularized model and then connects the thermal resistance of each component in series and parallel. In the calculation of conductive thermal conductivity, this method is used as an analytical expression for the conductive thermal conductivity of the entire porous structure, determined by analyzing a network of thermal resistances for an ideal periodic unit consisting of cell struts and/or cell walls [
165-
169].
Wang
et al. [
122,
123] proposed the equivalent circuit conductive thermal conductivity models of both closed-cell and open-cell porous polymers, taking the phonon scattering effect into account, as shown in Fig. 6. This method has given mathematical expressions of the conductive thermal conductivity consist of the void fraction (
, the strut fraction (
and the gas-to-solid thermal conductivity ratio (
. The conductive thermal conductivities of cell walls
and struts
were determined by analyzing the thermal resistance circuit of heat conduction based on the parallel-series and series-parallel methods for cubic foams [
122,
123], as
where
and
are the conductive thermal conductivities of cell walls based on parallel–series and series–parallel methods, and
and
are the conductive thermal conductivities of struts based on parallel–series and series–parallel methods, respectively [
136].
Based on this, Buahom
et al. [
136] introduced a series of alternative mathematical expressions for the conductive thermal conductivity in terms of the void fraction and the gas-to-solid thermal conductivity, as summarized in Table 1. These expressions are derived from those existing expressions based on heat conduction through closed-cell and open-cell foams.
Zeng
et al. [
154] and Wei
et al. [
168] also used equivalent circuit method to calculate the overall conductive thermal conductivity of aerogel materials, while Lu
et al. [
155] used this method to calculate the conductive thermal conductivity of aerogel composite materials.
Equivalent fractal method. Fractal theory can describe the disordered and random generation process of numerous porous materials at the micro- and nano-scale. These porous materials are different from the Euclidean description of proportional scale [
170] but satisfy the fractal characteristics [
171] for their self-similarity property [
172,
173], which represents it can be simplified the whole process by repeating unit. Based on these, equivalent fractal method is a numerical solution [
174,
175] for the conductive thermal conductivity is based on fractal geometry theory [
170,
176], which presents a random structure and is independent of the scale.
The heat transfer properties of microscale and nanoscale porous materials were studied based on these dispersed fractal models, such as the Sierpinski gasket model [
164,
177,
178], Vicsek model [
179,
180], the Von Koch snowflakes model [
181,
182], and the fractal-intersecting sphere model [
141,
142]. To better represent the structure of the porous polymeric material, the fractal-intersecting sphere model was proposed by Xie
et al. [
141,
142] to calculate the conductive thermal conductivity of aerogel, as shown in Fig. 7.
Combined the equivalent circuit method and geometric of the fractal-intersecting sphere model, the conductive thermal conductivity of the fractal-intersecting model at the first stage
and second stage
respectively was calculated as follows [
141,
142]:
where
is the conductive thermal conductivity of the primary size, obtained from [
183]
and /d is the aspect ratio of the solid particle radius and the cubic skeleton structure.
Numerical simulation method. Sundarram
et al. [
184] used finite element analysis (FEA) [
184-
186] and molecular dynamics (MD) [
187-
190] to study the heat transfer properties of porous polymers (in particular, PMMA and PEI). Moreover, the thermal conductivity was predicted to reduce when the pore size decreased from 1 mm to 1 nm, mainly attributed to the phonon scattering effect in the solid polymer matrix. In these researches, a variety of pore configurations could represent the micro/nanoporous structure observed in polymer foams by FEA method. In contrast, MD method treats individual atoms in the model as point masses that interact via established interatomic potentials [
184,
191]. In addition, it has also been widely used to predict physical properties of bulk materials [
192-
196].
Han
et al. [
197] and Zeng
et al. [
198] used Lattice Boltzmann method (LBM) [
199-
203] to calculate the overall conductive thermal conductivity of aerogel materials. Meanwhile, the scale effect of structural parameters and the temperature dependence of the thermal conductivity of nanoparticles are also studied using this simulation method.
Although the numerical method is complicated to calculate, it can reflect the influence of various factors on thermal conductivity, and it can also consider the situation of complex structures. However, the calculations in the current literature are still concentrated on the situation of some simple structures [
141]. What’s more, it should be noted that not to repeat considering the gas–solid interface effect in choosing the solid thermal conductivity model and setting the boundary conditions in gas–solid coupling model.
2.2 Radiative thermal conductivity models
The radiant heat transfer inside the polymer material belongs to medium radiation because the polymer material is participatory media for radiant heat transfer. When incident radiant energy enters the material’s interior, the material will absorb, refract, and reflect the radiation energy. Namely, the material has attenuation effects such as absorption or scattering of radiation [
141]. Compared to the relatively mature development of the conductive thermal conductivity model, a model for the radiative thermal conductivity of porous polymers is still challenging. The radiative thermal conductivity becomes significant, especially in largely expanded porous dielectric materials having higher porosity [
12,
22,
204,
205] and a smaller pore size [
122,
206-
208]. Several existing approaches for modeling the radiative thermal conductivity are based on different assumptions on the micro- and nanoscales.
2.2.1 Radiative transfer equation approximation
Recently, an accurate model for the radiative thermal conductivity was introduced by solving the equivalent radiant thermal conductivity through the radiative transfer equation (RTE) for simple structures [
141,
209,
210].
Assuming that in a participating medium that emits, absorbs, and scattering, at position
s, the radiation energy transfer direction
Ω, according to the conservation of energy, the governing equation for the transfer of radiant energy in the medium can be derived-the radiation transfer equation, as [
211,
212]
where,
is the spectrum radiation intensity corresponding to space position
s and transmission direction
s,
and
are the spectrum absorption and scattering coefficients [
213,
214] in the medium, respectively, and the sum of the two is the spectrum extinction coefficient
,
is the scattering phase function [
215,
216]. The radiation transfer equation is integrodifferential and nonlinear in the three-dimensional radiant heat transfer progress for non-gray materials unless in some simple cases, it is difficult to obtain the exact solution of the equation [
14,
141]. Therefore, it is necessary to approximate the differential-integral radiation transfer equation with certain assumptions before solving it [
141].
Two-flux approximation method (also called Schuster-Schwarschild approximation) [
209,
217,
218] is obtained by simplifying the integral-differential radiation transfer equation into a set of linear differential equations, assuming the radiant intensity uniformly distributed in a specific angle range. However, this method is only accurate through optical thinness media, and this method is not suitable for radiation heat transfer problems with strong anisotropic scattering [
141].
For the radiation intensity can be divided into several discrete directions in the entire space, Chandrasekhar [
219] developed the two-flux approximation by solving the radiation intensity in these discrete directions, known as the discrete-ordinates method [
220]. Besides these methods, there are different approximate methods to solve the radiative transfer equation, such as the finite volume method [
221-
225], discrete transfer method [
226,
227], etc.
2.2.2 Ignore some physical processes
When the absorption or scattering of the material is unapparent, and the related term in radiation transfer equation (2.2.1) can be ignored [
141]. Likewise, the influence of the radiation term is well known in conventional and microporous materials, and this term is negligible for low density foams (<0.2 relative densities) [
34,
228]. However, the conventional models used to evaluate the radiation mechanism in porous polymers assumes that the wavelength of the infrared radiation is smaller than the pore size. It should note that this presumption is incorrect in nanoporous polymers [
34,
228].
2.2.3 Rosseland diffusion approximation
The diffusion approximation, known as the Rosseland diffusion approximation [
209,
230], is well established for optically thick media, which is widely used to calculate the radiative thermal conductivity of porous polymers. If the medium is optically thick, the medium has a strong attenuation effect on radiation, and the transmission distance of radiant energy is very short [
209]. Moreover, the radiation transfer equation can be simplified based on this assumption [
209].
The radiative thermal conductivity
at a mean temperature
is usually determined by the corresponding temperature-dependent Rosseland diffusion equation [
209,
230], as
where
is the Stefan–Boltzmann constant (5.670367 ×10
-8 W
m
K
and
is the Rosseland mean extinction coefficient of the foam, which indicates attenuation of the radiative energy.
is a harmonic mean of the spectral extinction coefficient
weighted by the temperature derivative of the blackbody spectral intensity for all wavelengths within Planck’s blackbody emission spectrum [
231,
232], as
where
is the spectral distribution of Planck’s law of blackbody emission. It indicates the fraction of the radiative energy at a wavelength
λ to the overall energy over the wavelength spectrum [
122,
140], as
where
is the blackbody spectral intensity at a wavelength
λ. Two auxiliary constants are defined as
and
.
c is the speed of light in vacuum (2.99792458 ×10
8 m
s
.
h is the Planck constant (6.62607015 ×10
-34 J
s) and
kB is the Boltzmann constant (1.380649 ×10
-23 J
K
[
209,
233].
As shown in Fig. 3,
can be measured by experiments [
232] or calculated by Mie theory [
234,
235]. Accordingly, the fraction of the radiative energy reflected, absorbed and transmitted at a wavelength
λ in the radiation spectrum can be defined as the spectral reflectivity (
, spectral absorptivity (
and spectral transmissivity (
of cell wall, respectively. Because
, the radiative thermal conductivity could be obtained by measuring or calculating the transmissivity
, as shown in Fig. 3. Therefore, Gong
et al. [
22,
25] used the transmissivity to derive the radiative thermal conductivity formula of the porous polymers. Furthermore, Wang
et al. [
122,
123] considered the interference of the micro/nano-porous polymers’ electromagnetic waves and modified the radiative thermal conductivity formula based on the transmissivity formula, as shown in Fig. 8. This method takes the size effect into account, leading to a more accurate thermal radiation prediction especially in the micro/nano-porous polymers. The detail of this microscale radiation method would be exhaustively introduced in the following section.
2.2.4 Microscale radiation method
Since Chen
et al. [
236,
237] found that the radiant energy of two very closely located objects is much larger than that predicted by the Planck Stefan–Boltzman law [
237,
238], the thermal radiation process needs to consider microscale effect [
239,
240]. Maxwell’s equations describe the propagation of electromagnetic waves and their interactions with matter [
160,
241]. Furthermore, the fluctuational electrodynamics combining the fluctuation-dissipation theorem with Maxwell’s equations [
140,
242,
243] fully describe thermal radiation’s emission, propagation, and absorption in both the near and far fields [
211,
244,
245].
Buahom
et al. [
136] determined the radiative thermal conductivity by analyzing the attenuation of radiative energy by absorption and scattering based on Mie’s theory, as well as interference of propagating waves and tunneling of evanescent waves in the radiative energy.
This method illustrates that, in the porous polymer, radiation in the form of electromagnetic waves propagates in space and interacts with gas and polymer molecules [
136]. Each photon carries a single discrete quantum of thermal radiation and affects the electrons in the material molecule to increase its energy level [
136]. Then some of these excited electrons release the electron’s excess energy and emit photons in the wavelength range according to Planck’s law of blackbody emission [
136,
209]. At the same time, multiple internal scattering, absorption, and re-emission occur in the porous polymer, resulting in tortuous transmission paths for radiant energy transfer, especially in such a complex structure, leading to a tortious transport path of radiative energy transfer [
22,
136,
209].
As shown in Eq. (22),
can be obtained by assuming randomly oriented concave-triangular struts and cell walls independently as [
136]
where
is the contribution of struts and
is the contribution of cell walls, the volume fraction of polymer located in struts, the strut fraction
, is used as the weighted factor in the calculation. The effective spectral extinction coefficient of cell wall
can be calculated as [
136]
where
is the number of cell walls per unit volume of foam,
,
is the spectral extinction cross-section of cell wall for the incident beam interacts with a wall tilted by the incident angle
. the spectral extinction efficiency factor, the ratio of the spectral extinction cross-section to the geometric cross-section of cell wall at normal incident
. Based on Fresnel formula,
and
in Eq. (25) can be written as [
246]
where the spectral reflection and transmission coefficients at the gas–solid (
,
and solid-gas (
,
interfaces are defined as [
246]
where and are the spectral complex refractive indices of polymer film and non-absorbing gas.
The spectral extinction coefficient of the strut can be expressed in the number density of struts and the geometric cross-section [
136]:
where
is the number of struts per unit volume of porous polymers, which is dependent on the porosity, the cell size, and the strut fraction.
is the normally projected geometric cross-section of a strut.
is the asymmetry factor of the strut.
and
are spectral scattering and extinction efficiency factors of the strut, separately split into the TM and TE incident modes [
136].
3 Summary and prospects
In this study, the existing thermal conductivity calculation models of micro/nano-porous polymers are introduced and reviewed. Firstly, it reviews the process of establishing the overall equivalent thermal conductivity calculation model for micro/nano-porous polymers. Then, the thermal conductivity prediction model is introduced separately based on thermal conduction and thermal radiation. In addition, the thermal conduction part is divided into the gaseous thermal conductivity calculation model, solid thermal conductivity calculation model and gas–solid coupling model. Finally, it is concluded that, compared with other porous materials, there are few studies on heat transfer of micro/nano-porous polymers, especially on the particular heat transfer mechanisms such as scale effects at the micro/nanoscale. In particular, the following aspects of porous polymers still need to be further studied:
(i) The study of characteristics and influences in micro/nano-porous polymers brought by microscale thermal radiation.
(ii) The research of the heat transfer characteristics of open/ closed-pores and particular morphologies at the nanoscale.
(iii) Heat transfer mechanism and impact factors of micro/nano-porous polymers. By analyzing the various influencing factors on heat transfer performance, the porous polymers can be further optimized for efficient thermal insulation materials, which are well used in energy conversion and storage systems.