Novel Mixing Relations for Determining the Effective Thermal Conductivity of Open-Cell Foams.

This paper proposes a new approach to relate the effective thermal conductivity of open-cell solid foams to their porosity. It is based on a recently published approach for estimating the dielectric permittivity of isotropic porous media. A comprehensive assessment was performed comparing the proposed mixing relation with published experimental data for thermal conductivity and with numerical data from state-of-the-art relations. The mixing relation for the estimation of thermal conductivities based on dodecahedrons as building blocks shows good agreement with experimental data over a wide range of porosity.

1. Introduction

Effective permittivity, , and effective thermal conductivity, , of open-cell foams are important properties for the design and optimization of microwave-heated elements [1]. Determining both effective properties for foams with a reliable mixing relation based on foam porosity (ratio of the void volume to the total foam volume), bulk properties and skeletal morphology would greatly facilitate the design of systems using microwave-heated elements. Such “cross-property” relations are commonly used to interrelate changes in the effective values of various physical properties (e.g., thermal conductivity, elastic moduli, electrical conductivity and fluid permeability) [2,3,4,5] caused by microstructural features (e.g., pores and inclusions) of heterogeneous materials. Consequently, estimating effective properties turns out to be more complex than calculating a weighted volumetric mean. Predictive relations for or (or any cross-property) are based on models describing the microstructure of heterogeneous materials. Reliable predictions can only be obtained if the considered structural model resembles features of the foam microstructure. The geometrical representation of the foam morphology in the literature varies from the less accurate idealized assumptions (e.g., ordered, random and symmetrical distributions of solid and void phases) to more complex 2D structures (e.g., hexagonal honeycombs) and 3D unit cell morphology approaches (e.g., simple cubes, Weaire-Phelan unit cell and truncated tetrakaidecahedrons) [5,6,7,8,9,10,11,12].

Historically, the first relations formulated from idealized distributions of solid and void phases are based on the effective medium approximation (EMA) and belongs to the class of mean-field theories [5]. More recent relations are based on geometrical models considering details of the material structure, as well as structures modeled by using the finite element method (FEM).

In particular, for estimating the of open-cell foams, several relations can be found in the literature. These relations are mostly based on empirical data-fitting functions (e.g., Calmidi relation) [8], geometrical models (e.g., Bracconi, Dai and Yao relations) [9,10,11], probability distributions of parallel and series models (e.g., weighted arithmetic mean relation) [12,13,14], and those based on EMA assumptions, such as the Maxwell-type and self-consistent models (e.g., Differential Effective Medium relation) [12,15].

Recently, we proposed a numerical approach using FEM to derive predictive relations for the effective permittivity of open-cell foams based on two types of morphologies [16]. The first morphology corresponds to foams whose skeletons are based on Platonic solids (such as hexahedron, octahedron, icosahedron and dodecahedron) as building blocks, while the second morphology was reconstructed from micro-CT data of real open-cell ceramic foams. These relations based on Platonic solids and real open-cell foams are referred to as Platonic relation and OCF (open-cell foam) relation, respectively. Both relations agree well with the numerical data obtained from electromagnetic wave propagation calculations, as long as the foam behaves as an effective medium (effective medium approximation). In other words, on the macroscopic scale, foams behave similar to a homogenous medium. However, due to the lack of data, the relations are not yet validated at higher permittivity contrasts, where the EMA approach loses its applicability. Alternatively, the analogy of thermal and electrical networks can equally be used for the skeleton network of open-cell foams [9,17]. Here, extensive experimental and numerical data on the effective thermal conductivity of open-cell foams are available for validation. According to this analogy, it is reasonable to assume that any relation describing can be used to estimate , or vice versa as a cross-property relation [5,12,15,18]. Please note that only thermal conductivities are considered in this study and not convective heat-transfer coefficients or overall heat-transfer coefficients.

Following the thermal-electrical analogy, both Platonic and OCF relations can be used to estimate . The Platonic relation for estimating corresponds to the following: where is a correlation parameter that represents topological details of the skeleton morphology; is the porosity; and and are the thermal conductivities of the bulk materials, i.e., the solid skeleton and the medium that fills the voids of the skeleton, respectively. For complex-valued quantities, such as (), the calculation of also gives a complex-valued quantity () [16]. In contrast to , and, thus, also are real-valued quantities (). As a result, is calculated as follows: where and are expressed as follows: with = {,}. The respective coefficients, , are summarized in Table A1 of the Appendix A. Note that the Platonic geometry is no longer preserved if the size of the struts exceeds a certain limit and struts overlap each other, causing the closure of the open-cell faces. The porosity range in which foams exhibit the Platonic geometry is listed in Table A1. This corresponds to ideal Platonic skeletons and no to real foam skeletons where cell faces close at higher porosities (P ~ 0.5). The corresponding relation for estimating using the OCF relation is as follows:

Table A1

Parameters used to calculate and .

Valid Porosity Range	Variable	a_1	a_2	a_3	a_4	a_5	a_6
	Hexahedron
0.058 to 1	gm′	0.5229	−0.9951	2.4460	−4.7673	3.9566	−1.1630
	g0′	0.2775	4.7271	−14.1449	15.0204	−6.6742	0.7939
	Octahedron1
0.415 to 1	gm′	0.4008	0.6699	−3.2970	3.9169	−2.2435	0.5528
	g0′	2.1613	−7.0162	17.1386	−28.2326	23.5875	−7.6392
	Octahedron2
0.093 to 1	gm′	0.6546	−2.1904	6.0160	−9.7020	7.2485	−2.0266
	g0′	0.3698	7.6377	−26.3063	34.1470	−20.6464	4.7987
	Dodecahedron
0.078 to 1	gm′	0.4617	0.3546	−3.9082	7.6496	−7.0236	2.4661
	g0′	1.4799	−10.1720	50.1347	−109.2526	103.7717	−35.9641
	Icosahedron
0.267 to 1	gm′	0.2213	4.1649	−17.8917	29.0420	−21.8497	6.3132
	g0′	2.4320	−21.2239	84.7916	−152.1911	124.4845	−38.2943

The novelty of this work is the use mixing relations for predictions developed by using FEM electromagnetic wave propagation calculations to estimate the effective permittivity of foams. Therefore, Equations (1) and (4) represent novel mixing relations for estimating .

In the following assessment, predictions for the of open-cell foams using Platonic and OCF relations, as well as predictions from selected relations from the literature relations are analyzed and compared with experimental and numerical data.

2. Materials and Methods

The state-of-the-art relations [8,9,10,11,12,13,14,15] considered in this study are summarized in Table 1. They were selected because of their excellent prediction capability [8,9,10,11,12,13,14,15,19]. The considered skeleton materials of the open-cell foams and involved filling medium (as well as their bulk thermal conductivities, and ), and the open-cell foams data from which was obtained are summarized in Table 2 and Table 3, respectively [8,11,13,15,19,20,21,22,23,24,25,26,27,28].

Table 1

Mixing relations applied for estimating the effective thermal conductivity of open-cell foams.

Relation	Expression	Remarks
Bracconi [9]	k_eff=k_s(23(1−P)^2+13(1−P))	This relation ignores k_f because it was derived from a correlation for the tortuosity in periodic ordered structures, which is defined by the skeleton structure only.
Weighted arithmetic mean (WAM) [12]	k_eff=Ψ_arithm(Pk_f+(1−P)k_s)+1−Ψ_arithm(Pk_f+1−Pk_s)	Weighted arithmetic mean of the Wiener bounds using the following: Ψ_arithm=0.35 as proposed by Bhattacharya et al. [14] Ψ_arithm=0.49 as proposed by Dietrich et al. [13]
Calmidi [8]	k_eff=k_fP+k_s0.181(1−P)^0.763
Dai [10]	k_eff=22(R_A+R_B+R_C+R_D)	R_A=4d(2e^2+πd(1−e))(k_s−k_f)+4k_f
		R_B=(e−2d)e^2(k_s−k_f)+2k_f
		R_C=2(2−2e)πd^22(k_s−2k_f)+2k_f
		R_D=2ee^2(k_s−k_f)+4k_f
		d=(8(1−P)−3e^32)π(3−e(32+1))
		e=0.198
Differential Effective Medium (DEM) [15]	k_f−k_effk_f−k_s(k_sk_eff)^1/3=1−P	Known also as the Bruggeman relation, non-symmetric.
Yao [11]	k_eff=1(λ/k_E+(1−2λ)/k_F+λ/k_G)
		k_E=26πλ(3−4λ)1+b^2b^2(k_s−k_f)+k_f
		k_F=26πλ^21+b^2b^2(k_s2−k_f)+k_f
		k_G=26πλ^21+b^2b^2(k_s−k_f)+k_fλ is calculated (implicit method) from P=1−22πλ^2(3−5λ)1+b^2b^2,where a is a geometrical parameter (b=2.01 as recommended by Yao et al. [11]).

Table 2

Thermal conductivities of filling media and skeleton bulk materials (* average value [23]).

Skeleton Material	k_s/Wm^−1K^−1	Filling Medium	k_f/Wm^−1K^−1
Aluminum [8]	218	Air [11]	0.0265
Alumina [13]	25.9	Paraffin [11]	0.305
AlSi7 [15]	167	Vacuum * [23]	0.003
Cupper [11]	401	Water [11]	0.613
FeCr-alloy [28]	16
Mullite [13]	4.4
Nickel [25]	91.4
OBSiC [13]	8.1
Stainless steel (SS) [19]	15
Zirconia [28]	2.5
Polyurethane (PU) [19]	0.2

Table 3

List of references with published values for different skeleton and filling media combinations and corresponding thermal conductivity contrast values.

Skeleton-Fluid	k_s/k_f	Skeleton-Fluid	k_s/k_f
Al-Air [8,19,20,21,22]	8.2×10^3	Nickel-Air [19]	3.4×10^3
Al-Water [8,19,22]	355.6	Nickel-Water [19]	149.1
Al-Paraffin [19]	714.7	Nickel-Paraffin [19,27]	299.7
Al-Vacuum [23,24]	7.3×10^4	Polyurethane-Air [19]	7.5
Alumina-Air [13]	977.4	Polyurethane-Water [19]	0.3
AlSi7-Air [15]	6.3×10^3	Polyurethane-Paraffin [19]	0.7
Cu-Air [11,19]	1.5×10^4	OBSiC-Air [13]	305.7
Cu-Paraffin [15,22,25]	1.3×10^3	Stainless steel-Air [19]	566.0
Cu-Water [11,19]	654.2	Stainless steel-Water [19]	24.5
FeCr-alloy-Air [26]	603.8	Stainless steel-Paraffin [19]	49.2
Mullite-Air [13,26]	166.0	Zirconia-Air [26]	94.3

For easier data comparison, the of open-cell foams is normalized as follows:

This way, all values are scaled between 0 to 1, allowing us to describe as a function of , which depends on the foam microstructure and thermal conductivity contrast. The difference between estimated from relations and those from experimental or simulated data is quantified by the root-mean-square error (RMSE).

3. Results and Discussion

Figure 1 shows the normalized data (symbols) found in the literature for open-cell foams against estimations (lines) using the literature mixing relations (upper figures), the Platonic relations (lower figures using Equation (1) based on different Platonic solids) and the OCF relation (lower-right subfigure using Equation (4)). As increases, decreases for both experimental data and estimates from relations (as illustrated in the right figure of the upper row for the prediction of the Dai relation [10]). In addition, Table 4 summarizes the RMSE of each relation at three different porosity ranges ([0.5, 1.0], [0.85, 1.0] and [0.9, 1.0]) to characterize the deviations from the experimental data.

Figure 1

Normalized experimental values (represented by symbols corresponding to foams listed in Table 2) compared with those estimated from relations (represented by lines corresponding to relations listed in Table 1 and Equations (1) and (4)). The embedded subplots provide an enlarged view for the porosity ranging from 0.85 to 1.0.

Table 4

Comparison of the RMSE for of the considered relations (* weighted arithmetic mean approach with corresponding arithmetic coefficient proposed by Bhattacharya or Dietrich; see remarks in Table 1).

Relation	RMSE
	P≥0.50	P≥0.85	P≥0.90
Bracconi	1.55	0.78	0.68
Calmidi	7.80	0.98	0.81
Dai	84.05	0.80	0.73
DEM	2.10	1.77	1.72
WAM/Bhattacharya *	5.38	0.68	0.58
WAM/Dietrich *	3.58	2.32	2.13
Yao	3.02	0.62	0.53
Hexahedron	1.85	1.29	1.14
Octahedron1	1.89	1.34	1.18
Octahedron2	1.90	1.32	1.18
Dodecahedron	1.65	0.87	0.73
Icosahedron	1.67	1.07	0.92
OCF	2.30	1.55	1.52

Figure 1 (upper row) illustrates that the literature relations estimate better for , and this is expected, because they were preferably developed for foams with higher porosity. The best predictions are obtained by the relation from Bracconi (see Table 4). However, it should be mentioned that the Bracconi relation accuracy decreases as increases, since the relation does not include a term.

The relations inspired by the Platonic solids (bottom row of Figure 1) provide excellent predictions with only minor deviations from the experimental data (see Table 4). The results suggest that the dodecahedral structure (followed by the icosahedral) best mimics the skeleton of real foams. It is important to highlight that other researchers [10,22] have proposed the Kelvin tetrakaidecahedron (known to reproduce packings with low surface area) as the best geometrical element to mimic foams. However, the modeling of this polyhedron has been performed only in two dimensions, using ligaments corresponding to struts and thus lacks important features, such as the effect of the geometrical shape of joints. Accordingly, a better description of tetrakaidecahedron-inspired foams would be required for a fair comparison with the dodecahedral and the icosahedral relations from this study.

The predictions from the OCF relation (see lower right graph of Figure 1) are similar to those of the DEM relation, which agrees with our previous study on the estimation of [16] but deviates—to a certain extent—from the experimental data of .

Figure 2 shows the estimated from the relation of Bracconi [9] and from dodecahedron and OCF relations in comparison with those reported from numerical simulations [19].

Figure 2

estimated from Bracconi, dodecahedron and OCF relations compared with those from numerical simulations of Al, Cu and Ni foams (left) and SS and PU foams (right). The errors bars correspond to the standard deviation as reported by August et al. [19]. Lines from Bracconi’s relation are overlapping, as it does not consider the filling medium.

Figure 2 reveals that, in contrast to the experimental values, the best predictions for the simulated data are obtained from the OCF relation (RMSEOCF = 0.41 < RMSEBracconi = 1.27 < RMSEDodecahedron = 1.77). This is consistent with previous electromagnetic wave propagation calculations used for computing [16] and with calculations via the diffuse interface representation of the phase-field model used for computing [19]. The structural models for computing were reconstructed from tomographic scans of open-cell foams, while the structural models for computing correspond to synthetic foam structures using the algorithm proposed by August et al. [29]. Although the open-cell structural models used for the simulations are different, both are well represented by the OCF relation. Thus, it can be concluded that the simulation models are significantly different from real foams for the following reasons:

For the simulations performed to calculate , the skeleton morphology was reconstructed from µCT-scans of samples with porosities of 0.90 ± 0.01. At this porosity, both experiments and simulations are well estimated by the OCF relation (see Figure 1 and Figure 2). However, 3D erosion and dilation filters [16] were applied for generating models of different porosity, which—depending on the mesh resolution—may produce significant differences compared with the microstructure of real foams.

The synthetic foam structures generated by August et al. [19] are not morphologically identical despite having the same porosity. Numerically computed and measured experimental values of (as reported by August et al. [19]) only agree if the standard deviation is considered. This indicates that only a few of the synthetic structures are morphologically consistent with real foams.

4. Conclusions

Recently, relations have been derived to estimate the effective permittivity of open-cell foams based on two approaches: (1) using Platonic solids as building blocks of foam skeletons and (2) using the morphology of foam samples extracted from tomographic scans. Based on the thermal-electrical analogy, these relations can be used to estimate cross-properties, such as the effective thermal conductivity. In this work, an assessment of the predictions of the effective thermal conductivity of open-cell foams from these new relations has been presented. The relations have been compared with experimental and numerical data from the literature as well as with predictions from available mixing relations. It has been shown that foam’s thermal conductivity can be well estimated from the Platonic relation by using dodecahedrons (which describes foams based on dodecahedral building blocks).

The following recommendations can be derived from this work to properly select the most suitable relation:

For foams with porosities ranging from 0.5 to 1.0 and low bulk thermal conductivities of the filling medium, the Bracconi relation [9] is recommended.

Novel mixing relations are recommended for porosities ranging from 0.5 to 0.9, except for the Platonic relation based on dodecahedrons, which can be applied for porosities from 0.5 to 1.0. In addition, the Platonic relation based on dodecahedrons is recommended over the Bracconi relation for thermal conductivities of the filling medium higher than  .

The relations of Yao [11] and the weighted arithmetic mean approach with arithmetic coefficient proposed by Bhattacharya [14] are recommended if the porosity is higher than 0.85.

Finally, the simulated numerical data were well met by using the OCF relation. However, a significant difference was identified between the numerically predicted values and those from experiments.