Home | About | Journals | Submit | Contact Us | Français |

**|**Materials (Basel)**|**v.7(11); 2014 November**|**PMC5512629

Formats

Article sections

- Abstract
- 1. Introduction
- 2. Experimental
- 3. Morphological Characterization
- 4. Heat Conduction Modelling
- 5. Summary and Conclusions
- References

Authors

Related links

Materials (Basel). 2014 November; 7(11): 7173–7195.

Published online 2014 October 28. doi: 10.3390/ma7117173

PMCID: PMC5512629

Javier Narciso, External Editor

Received 2014 September 3; Revised 2014 October 13; Accepted 2014 October 15.

Copyright © 2014 by the authors;

licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/).

Reticulated porous ceramic (RPC) made of ceria are promising structures used in solar thermochemical redox cycles for splitting CO_{2} and H_{2}O. They feature dual-scale porosity with mm-size pores for effective radiative heat transfer during reduction and µm-size pores within its struts for enhanced kinetics during oxidation. In this work, the detailed 3D digital representation of the complex dual-scale RPC is obtained using synchrotron submicrometer tomography and X-ray microtomography. Total and open porosity, pore size distribution, mean pore diameter, and specific surface area are extracted from the computer tomography (CT) scans. The 3D digital geometry is then applied in direct pore level simulations (DPLS) of Fourier’s law within the solid and the fluid phases for the accurate determination of the effective thermal conductivity at each porosity scale and combined, and for fluid-to-solid thermal conductivity from 10^{−5} to 1. Results are compared to predictions by analytical models for structures with a wide range of porosities 0.09–0.9 in both the strut’s µm-scale and bulk’s mm-scale. The morphological properties and effective thermal conductivity determined in this work serve as an input to volume-averaged models for the design and optimization of solar chemical reactors.

Foam-type reticulated porous ceramics (RPC) structures are applied in a broad range of physical processes requiring enhanced heat and mass transfer [1,2]. Applications include microelectronics cooling [3], soil dynamics [4,5], catalytic reactors [6], radiant burners [7], tissue engineering [3,8] and volumetric heat exchangers for the conversion of concentrated solar energy [9,10,11]. Of special interest of the latter application is the solar-driven thermochemical redox cycle for splitting CO_{2} and H_{2}O [12,13,14,15,16], consisting of: (1) a high-temperature endothermic reduction, in which a metal oxide is thermally reduced and oxygen is evolved; and (2) a lower-temperature exothermic oxidation, in which the reduced oxide is re-oxidized with H_{2}O and CO_{2} to form H_{2} and CO (syngas), and further processed to liquid hydrocarbon fuels. Ceria-based oxides have emerged as highly attractive redox materials because of the rapid oxygen transport in the bulk [14,15,16,17,18]. Various porous structures made of ceria have been investigated for enhanced reaction rates [18,19,20], including structures with submicron-sized interconnected pores, but these are problematic to retain because of partial sintering at elevated temperatures [19]. Furthermore, their high optical thickness inhibits penetration of concentrated solar radiation, resulting in non-uniform heating and temperature distributions [14]. Most recently, Furler* et al.* [15] presented a unique and morphologically stable RPC structure featuring dual-scale porosity: mm-size pores with struts containing micron-size pores. The mm-size pores enable volumetric absorption of concentrated solar radiation and effective heat transfer during the reduction step, while the micron-size pores within the struts offer increased specific surface area leading to enhanced reaction kinetics during the oxidation step.

Optimization of solar reactors for thermochemical redox cycles requires computational models of heat transfer and fluid dynamics coupled to the reaction kinetics [14,21]. Since resolving the solar reactor at the pore scale would require tremendous computational demand, volume-averaging theory is often applied for solving the mass, energy, and momentum conservation equations using effective heat and mass transport properties [22,23,24,25]. These can be determined accurately by direct pore-level simulations (DPLS) using the detailed 3D digital geometry of the structure obtained by computer tomography (CT) [26,27,28]. For example, the Monte Carlo ray-tracing method has been applied at the pore level for solving the radiative heat transfer equations and determining the effective extinction coefficient and scattering phase function [29], and the finite volume (FV) technique has been applied at the pore level for solving the Navier-Stokes equations and determining the effective thermal conductivity, permeability, and heat transfer coefficient [30,31].

In this work we apply the tomography-based methodology to investigate RPC structures made of ceria with dual-scale porosity in the mm and µm scales. This structure is schematically depicted in Figure 1 [15]. The total and open porosity, pore size distribution, mean pore diameter and specific surface area are extracted from the CT-scans. The effective thermal conductivity is determined by DPLS for the RPC with non-porous struts and for the RPC with dual-scale porosity. We investigate the effect that the dual-scale porosity has on the morphological properties and on the conduction heat transfer across the RPC, and further compare the results to predictions by analytical models for structures with a wide range of porosities in both the strut’s µm-scale and bulk’s mm-scale.

The dual-scale RPC structure is manufactured using the Schwartzwalder foam replication method [32]. An organic foam template is coated with multiple slurry layers containing ceria particles and micron-sized carbon grains [15]. The carbon pore former content ranges from 10 to 50 vol%. After firing at high temperatures (>1800 K), the bearing carbon foam and grains are burned und the desired foam-type structure undergoes sintering.

Strut samples are scanned using synchrotron submicrometer tomography with a voxel (3D pixel) size of *v*_{s} = 325 nm and a 0.832 × 0.832 × 0.702 mm^{3} field of view. The high-resolution CT is performed at the Swiss Light Source (SLS) of the Paul Scherrer Institute (PSI, Villigen, Switzerland) with the TOMCAT beamline for 40 keV photon energy, 400 µA beam current, a 100 µm thick aluminium filter, 40 µm thick copper filter, a 10 µm thick iron filter, 20× geometrical magnification, 1 s exposure time, 1001 projections. Figure 2 shows exemplary tomograms of strut samples manufactured with various concentrations of pore former ranging from 10 to 50 vol% and their corresponding 3D reconstructions of the pore space within isotropic strut regions. One tomogram contains 2560 × 2560 pixels. Numerous strut samples are scanned to verify reproducibility.

A ceria RPC sample with 10 pores per inch (ppi) is scanned by micrometer tomography with a voxel size of *v*_{s} = 35.7 µm and a 36.56 × 36.56 × 36.56 mm^{3} field of view. The low-resolution CT is performed with an unfiltered polychromatic X-ray beam at the Swiss Federal Laboratories for Materials Science and Technology (EMPA, Dübedorf, Switzerland) for 150 keV photon energy, 45 µA beam current, 6.272 s exposure time and 721 projections. Figure 3 shows a tomogram and its corresponding 3D digital reconstruction of the scanned RPC sample. One tomogram contains 1024 × 1024 pixels.

For the analysis of the strut structures, isotropic regions within submicrometer tomograms are cropped with the size of 501 × 501 × 501 voxels (0.163 × 0.163 × 0.163 mm^{3}). For the analysis of the RPC structures, micrometer tomograms are cropped to 500 × 500 × 500 voxels (17.85 × 17.85 × 17.85 mm^{3}). The cropped 8 bit tomograms obtained from measurements are pre-processed with a 3D Gaussian blurring filter to remove unwanted image noise derived from the photon sensor. Histograms computed from 3D tomogram stacks show bimodal character representing two grey scale pixel classes, with the threshold found by Otsu’s method of intra-class variance minimization [33,34]. Finally, based on the threshold, each pixel is assigned to be either void or solid. Porosity, ε = *V*_{f}/*V*, is defined as the ratio between the void space volume and the total cube volume. It is calculated by counting void and solid voxels of the 3D stack. The representative elementary volume (REV) defines the minimum volume containing a porous zone for which the continuum assumption is valid. It is determined from incrementally growing cubic subvolumes until their calculated porosities convergences within a certain band, ±γ. The conditions for the minimum edge length of the REV are [27]:

(1)

where *V _{L*}* is the sample subvolume and

Dual-scale porosity, ε_{dual}, is calculated from strut-scale porosity, ε_{strut}, and RPC-scale porosity, ε_{RPC}, as:

ε_{dual} = ε_{RPC} + (1 − ε_{RPC}) ε_{strut}

(2)

ε_{strut} is linearly fitted to the pore former concentration, ϕ, as ε_{strut} = 0.008707 ϕ. RPC structures with varying strut thicknesses are generated by altering the original segmented tomography scans through a dilation process with 3D spherical elements of a certain diameter, *d*, with . Table 1 lists the dilation radius, digital porosity, and mean pore diameter of the original RPC reconstruction and of the digitally altered RPC for 3 increasing strut thicknesses, and the corresponding digital section cut and 3D rendering. As expected, porosity and mean pore diameter of the RPC decrease with increasing strut dilation because the thicker struts consume void space.

Dilation radius, digital porosity, and mean pore diameter of the original RPC reconstruction and of the digitally altered RPC for increasing strut thickness, and the corresponding digital section cut and 3D rendering.

Of special interest is the connectivity of the µm-sized pores within the struts. The pore connectivity scales directly with the specific surface area reachable by reacting gases, and thus scales with the fuel production rates [15]. Open porosity, ε_{open}, is defined as the pore space accessible from one of the 6 cube sample surfaces. An iterative routine starts searching from one cube side for connected neighbour void voxels in order to find all pores connected to this side. This reconstruction is performed to detect closed pores within the structure that account for the porosity but are not exposed to the gaseous reactants. Table 2 shows exemplary results of total and open porosity for strut samples with various pore former concentrations.

Strut and open porosity, mean pore diameter, and the corresponding 3D rendering for strut samples manufactured with various pore former concentrations.

Figure 4 shows all results of the porosities collected from selected 3D tomography reconstructions. ε_{open} is presented with error bars because 6 evaluations are obtained per sample (one from each side). For ϕ ≤ 20 vol%, there is no pore connectivity observed as seen graphically in Table 2. For ϕ ≥ 30 vol%, the majority of the pores are connected and the pore network passes through the entire cube sample. For 50 vol%, practically every pore is connected to the pore network since the open porosity and total porosity are almost the same ε_{open} ≈ ε_{strut}. For high RPC porosity, e.g., ε_{RPC} = 0.825, the dual-scale porosity changes by less than 0.1 for ϕ = 50 vol% because only 17.5% of the volume is filled with µm-sized pores within the thin struts (see Figure 4). For lower RPC porosities, e.g., ε_{RPC} = 0.459, there is a stronger trend for increasing ϕ since there is more solid to be filled with µm-sized pores. The porosities of the model fit for ϕ = 0 vol% represent ε_{RPC} as listed in Table 1.

Total and open porosity as a function of the pore former concentration used to manufacture the µm-sized strut pores, total dual-scale porosity of RPC (original and digitally altered strut thickness).

To investigate the transition from closed pores to interconnected pores in a systematic way, spheres with a small Gaussian size distribution of 2–3 voxels are randomly placed within the 3D volume. Pores of such artificially generated structures get connected between ε_{strut} = 0.2 and 0.3. Such connectivity behavior is consistent with empirical correlations of the effective gas diffusivity within porous carbon [35]. Figure 4 also shows the comparison between the numerically determined ε_{open} and the experimentally measured values by mercury intrusion porosimetry (MIP) [36,37], as reported by Furler* et al.* [15]. The agreement is reasonable well. For ϕ = 30 vol%, the experimentally measured ε_{open} is lower than the numerically determined one due to poor pore connectivity across coating layers, as seen in Figure 5. This anisotropic region results from the 2-step coating applied during the fabrication process, and it is not considered in the determination of morphological properties and effective thermal conductivity within isotropic regions.

Pore size distribution is determined by applying a morphology-altering algorithm to the digitally segmented 3D structure consisting of an inversion of the solid and void space, followed by erosion and dilation with successively increasing spherical elements of diameter *d* with . In a last step, the algorithm inverts back the solid and void space. The cumulative pore size distribution 1 − *F*(*d*) is defined as the ratio of the opening-closing porosity, ε_{oc}(*d*), and the original porosity [27]:

(3)

The pore size distribution is then calculated as: *f*(*d*) = d*F*(*d*)/d*d*. Figure 6 shows the cumulative pore size distribution (left *y*-axis) and pore size distribution (right *y*-axis) as a function of *d*.·*f*(*d*) agrees qualitatively well with the values obtained by mercury intrusion porosimetry measurements. The mean pore diameter is then defined in such a way to split the area under the pore size distribution curve into two equal areas [27].

(4)

Cumulative pore size distribution (**left**
*y*-axis) and pore size distribution (**right**
*y*-axis) obtained by morphology operations with spherical structuring elements of diameter *d*.

The mean pore diameter for porous strut samples with various concentrations of pore former in the range 10–50 vol% is listed in Table 2. As expected, *d*_{mean} remains independent of ϕ at around 10 µm because the same pore forming material is used to manufacture all samples.

The specific surface area (SSA) is determined in three different ways: (1) using statistical two-point correlation function computed on the 3D segmented structures with an in-house Fortran code; (2) resampling of the phase interface area with a surface mesh-based algorithm using the open source software ImageJ (version 1.47v, Java 1.6.0_20 (64-bit)) [38] extended with the free BoneJ plugin (version 1.3.11) [39]; and (3) using the actual phase interface area of tetrahedral 3D meshes generated with an in-house Fortran code [40]. The in-house mesh generator covers the void and solid domain with tetrahedral elements for unstructured body-fitted grids and subsequently refines the elements at the phase boundary. These 3D meshes are later used to solve the steady-state energy conservation equation to determine the effective thermal conductivity. The two point correlation *s*_{2}(*r*) is a statistical function that indicates the probability of two arbitrary points Ψ(*r*) and Ψ(**r** + *r***ŝ**) separated by the distance *r* to be in the void phase [41]:

(5)

where Ω is the solid angle and *V* the cube volume. Porosity and specific surface area are then calculated using the following expressions [27,30]:

(6)

(7)

(8)

A code loops once through the entire structure (*x*-, *y*-, *z*-direction), counts each void voxel, and computes *s*_{2}(*r* = 0) in a digitally exact manner. Additionally, the code counts for each void voxel the number of direct neighbour void voxels in all 6 directions (-*x*, *x*, -*y*, *y*, -*z*, *z*) and computes 6·*s*_{2}(*r* = 1). The volumetric specific surface area is then calculated as:

(9)

The specific surface area is presented in two different units: per total volume (fluid+solid phases), *A*_{0} [m^{2} · m^{−3}], and per ceria mass, *ssa* [m^{2} · g^{−1}]. The density of ceria is ρ_{CeO2} = 7.22 g/cm^{3} [21,42]. Dual-scale *A*_{0} is calculated by multiplying *A*_{0} determined for the struts by the solid volume fraction of the RPC: *A*_{0,dual} = (1 − ε_{RPC}) · *A*_{0,strut}. *A*_{0} of RPC with non-porous struts is converted to *ssa* by: *ssa*_{RPC} = *A*_{0,RPC}/[(1 − ε_{RPC}) · ρ_{CeO2}]. Table 3 lists *ssa* of RPC with non-porous struts obtained from the original tomography scans and digitally dilated struts.

Because *ssa* is defined per unit mass, the strut *ssa* is necessarily equal to the *ssa* of the entire structure, including mm and µm sized pores. Thus, *A*_{0,strut} is used to calculate *ssa* of porous struts and dual-scale porous structures:

(10)

*A*_{0} and *ssa* are plotted as a function of ϕ in Figure 7a,b, respectively. Of special interest is the open *ssa*, which is directly related to the surface area reachable by the reacting gases for conversion of CO_{2} and H_{2}O to CO and H_{2}. The open *ssa* is calculated for the struts without closed pores. Open *A*_{0} and *ssa* are presented with error bars (standard deviation) for the data evaluated. Surface areas calculated from two-point correlation and phase interface area of the 3D meshes (black and white symbols) lie within close proximity to one another, whereas those calculated using ImageJ (grey symbols) are higher. This is because the two-point correlation is based on a statistical model leading to smoothing effects and the 3D meshes actually contain a smoothing algorithm for the phase interface, whereas the resampled phase interface of ImageJ incorporates fine mesh surface irregularities. For ϕ ≤ 20 vol%, the open *ssa* is smaller than the total *ssa* because the majority of the pores are not connected. For ϕ = 30 vol% most pores are connected and the open *ssa* approaches the total *ssa*. For ϕ = 50 vol%, open and total *ssa* are nearly identical because at this point there are practically no closed pores (see Table 2). Trend lines are plotted defining the mean value of all different calculation methods for total (solid) and open (dashed) *A*_{0} and *ssa*.

The governing steady-state heat conduction equations within the solid phase and the stagnant fluid phase are given by:

Solid phase: (*k*_{s}*T*_{s}) = 0

(11)

Fluid phase: (*k*_{f}*T*_{f}) = 0

(12)

where *k*_{s} and *k*_{f} are the thermal conductivity of the solid and fluid, respectively. The cubic domain is schematically shown in Figure 8. The boundary conditions are given in Equations (13)–(17). An inlet and outlet temperature is set (*T*_{hot} > *T*_{cold}) to provide a steady heat flux through the two phases with length *L*. Lateral walls of the sample cube are adiabatic. Local thermal equilibrium is assumed at the phase interface. Heat flux across the interface is driven by the temperature gradient and the thermal conductivities in each phase at the interface.

Schematic of the steady state heat conduction simulation setup with hot inlet temperature, *T*_{hot}, at *z* = 0 and cold outlet temperature, *T*_{cold}, at *z* = *L*.

Inlet temperature: *T*_{s}(*z* = 0) = *T*_{f}(*z* = 0) = *T*_{hot}

(13)

Outlet temperature: *T*_{s}(*z* = *L*) = *T*_{f}(*z* = *L*) = *T*_{cold}

(14)

Adiabatic lateral walls: **n** · = 0

(15)

Local thermal equilibrium at phase interface: *T*_{s} = *T*_{f}

(16)

Heat flux across phase interface: **n**·*k*_{s}*T*_{s} = **n**·*k*_{f}*T*_{f}

(17)

The governing volume-averaged steady-state equation for effective heat conduction within the isotropic porous structure reduces to one equation [22,23]:

(*k*_{eff}*T*) = 0

(18)

The effective thermal conductivity is calculated using the 1D Fourier’s law and the heat flux determined by DPLS:

(19)

where *k*_{eff} is the effective thermal conductivity of the cubic porous structure, the effective heat flux at the inlet or outlet, and *A*_{flux} = *L*^{2} is the inlet or outlet area constraint with *T*_{hot} or *T*_{cold}, respectively. The methodology for determination of *k*_{eff} for dual-scale porous structures is schematically shown in Figure 9. In a first step, *k*_{eff} of the strut with µm-sized pores (*k*_{eff,strut}) is determined according Equation (19). In a second step, this *k*_{eff,strut} serves as an input for the solid domain of a further simulation performed with the mm-sized pores of the RPC.

Methodology for the determination of the effective thermal conductivity of the RPC with dual-scale porosity.

Numerical DPLS are performed for RPC with non-porous/porous struts with/without digital strut dilation for different fluid-solid thermal conductivity ratios ranging from 10^{−5} up to 1. The cases covered include 4 RPC with non-porous struts (original scan and digitally dilated struts with 2, 5, and 10 voxels) and 16 RPC with porous struts (*i.e.*, 4 RPC, each with 4 different strut porosities ϕ* =* 10, 20, 30, 50 vol%). Simulations are performed using a commercial computational fluid dynamics (CFD) software (ANSYS^{®} Academic Research, release 14.0). Initially, grid resolution study is performed, indicating convergence for structures containing element sizes between 0.57 µm at fluid-solid interface to 2.28 µm within bulk for the µm-size pores within the struts, and between 62.0 µm at fluid-solid interface to 247.9 µm within bulk for the mm-size pores of the RPC. Typical number of elements is 20 million, yielding an error of less than 1% compared to the finest mesh tested of 40 million. The correctness of the DPLS was verified by solving simple geometrical cases with exact analytical solutions, while its accuracy was fine-tuned by grid refinement.

Figure 10 shows the ratio of the effective thermal conductivity to the solid thermal conductivity *vs.* the ratio of the fluid-to-solid thermal conductivity for a single porous strut (ϕ = 50 vol%), a RPC with non-porous struts, and a RPC with porous struts (ϕ = 50 vol% and 0 mm strut dilation). The analytical curves for serial and parallel slabs are indicating the maximum and minimum possible heat flux [43,44]. These exemplary simulation results correspond to ε_{strut} = 0.410 and ε_{RPC} = 0.825, leading to ε_{dual} = 0.897. *k*_{eff} decreases with increasing porosity and decreasing *k*_{f}/*k*_{s}. For *k*_{f}/*k*_{s} < 10^{−3}, *k*_{eff} does not significantly change anymore, indicating heat conduction dominated by the solid domain (e.g., for vacuum applications). In that range, the ratio of *k*_{eff} for RPC with non-porous struts to *k*_{eff} for RPC with porous struts is 2.4. As expected, this ratio approaches 1 for increasing *k*_{f}/*k*_{s} as the thermal conductivities of the fluid and solid phases approach each other. The serial and parallel heat conduction mode of the lumped fluid and solid material bracket the minimum and maximum possible heat flux (also called Wiener lower and upper bound), respectively [43,44]. As expected, the simulation results (black symbols) are between the minimum (white symbols with dashed line) and maximum (white symbols with solid line) possible heat flux for each porosity.

Ratio of the effective thermal conductivity to the solid thermal conductivity* vs.* ratio of the fluid-to-solid thermal conductivity for a single porous strut (ϕ = 50 vol%), a RPC with non-porous struts, and a RPC with porous struts (ϕ = **...**

Table 4 lists various analytical models for *k*_{eff} [31,45]. For simplicity, analytical equations are given in terms of η = *k*_{f}/*k*_{s}, and ς_{eff} = *k*_{eff}/*k*_{s}. Several models allow fitting with geometrical shaping parameters. They were least-squares fitted to three different sets of simulation data: (1) *k*_{eff,strut} for a single porous strut (black squares shown in Figure 10); (2) *k*_{eff,RPC} for a RPC with non-porous struts (black circles shown in Figure 10); and (3) *k*_{eff,dual} for a RPC with porous struts (black triangles shown in Figure 10). To identify the model which agrees best with all data sets, an overall least-squares approximation, *k*_{eff,all}, was fitted for all simulation data.

The root-mean-square error (RMS) is defined to compare *k*_{eff} calculated by the analytical models with that determined by our simulation:

(20)

where *n* is the number of data points per data set over the entire range of *k*_{eff}*/k _{s}* and

Root-mean-square (RMS) error of analytical models compared to three simulation data sets and to all simulation data.

High and low resolution computer tomographic scans were performed on complex reticulated porous ceramics (RPC) structures to capture the 3D digital representations of their dual-scale porosity in the mm and µm range. The CT scans were processed with a Gaussian blurring filter for a clustering-based image thresholding of the void and solid phases using Otsu’s method. The struts containing µm-size pores were digitally dilated with spherical structuring elements generating structures with different thickness and porosity. The morphological properties analyzed include porosity, pore size distribution, specific surface area, and pore connectivity within representative sample volumes of the isotropic strut regions and of the RPC. The total strut porosity was linearly dependent on the concentration of pore forming agent, and no pore connectivity was observed for concentration less than 20 vol%, consistent with mercury intrusion porosimetry measurements. A well-connected pore network results in high specific surface area and penetration of reactant gas for high fuel production. The effective thermal conductivities of a single porous strut, a RPC with non-porous struts, and a RPC with porous struts (dual scale) were determined by direct pore level simulations of the heat conduction equation with a CFD code. Values were compared to predictions by analytical models over a wide range of porosities. Models without shaping parameters were generally inaccurate (overall RMS > 10%). Miller’s model with two shaping parameters predicted *k*_{eff} with RMS error below 2.1% and the modified three-resistor model with three empirical fitting parameters predicted *k*_{eff} with a RMS error below 1.5%. These analytical correlations are applicable to RPC with porosities in both the strut’s µm-scale and bulk’s mm-scale ranging from 0.09 to 0.9.

The morphological properties and effective thermal conductivity determined in this work serve as an input to volume-averaged models for the design and optimization of solar chemical reactors.

We gratefully acknowledge the financial support by the Swiss Competence Center Energy & Mobility, the Helmholtz–Gemeinschaft Deutscher Forschungszentren (Virtuelles Institut SolarSyngas), and the European Research Council under the European Union’s ERC Advanced Grant (SUNFUELS—No. 320541).

Nomenclature

*A*_{0,(.)}- Volumetric specific surface area (m
^{−1}) *A*_{flux}- Cross sectional inlet/outlet area of cubic sample for heat flux (m
^{2}) *ssa*_{(.)}- Physical specific surface area (Index: strut, RPC, dual) (m
^{2}g^{−1}) *d*- Diameter of spherical structuring elements (m)
*d*_{mean}- Mean pore diameter (m)
*f*(*d*)- Pore size distribution (–)
*F*(*d*)- Cumulative pore size distribution (–)
*k*_{s}- Solid thermal conductivity (W m
^{−1}K^{−1}) *k*_{f}- Fluid thermal conductivity (W m
^{−1}K^{−1}) *k*_{eff,(.)}- Effective thermal conductivity of porous structure (W m
^{−1}K^{−1}) *l*_{REV}- Cube edge length of representative elementary volume (m)
*L*- Cube edge length (m)
*n*- Number of simulation data points (–)
- Heat flux through sample (W m
^{−2 }) *s*_{2}(*r*)- Two point correlation function (–)
*T*_{cold}- Cold side temperature (K)
*T*_{f}- Fluid temperature (K)
*T*_{hot}- Hot side temperature (K)
*T*_{s}- Solid temperature (K)
*v*_{s}- Voxel size (m)
*V*- Total sample cube volume (m
^{3}) *V*_{f}- Void volume (m
^{3}) - ε
_{(.)} - Porosity (–)
- γ
- Error band of porosity (–)
- ς
_{eff} - Ratio of effective to solid thermal conductivity (–)
- η
- Ratio of fluid to solid thermal conductivity (–)
- ϕ
- Pore former concentration (vol%)

- 2pc
- 2-point correlation
- mesh
- 3D-mesh generated from digitally segmented structures
- ImageJ
- Calculated using open source software ImageJ
- strut
- Morphological property of porous strut
- RPC
- Morphological property of RPC with non-porous struts
- dual
- Morphological property of RPC with porous struts (dual-scale)
- open
- Open morphological property excluding closed pores
- oc
- Morphological property after applying opening-closing algorithm

Conflicts of Interest

The authors declare no conflict of interest.

1. Bear J., Buchlin J.M. Modelling and Applications of Transport Phenomena in Porous Media. Kluwer Academic Publishers; Dordrecht, The Netherlands: 1991. p. 380.

2. Ho C.K., Webb S.W. Gas Transport in Porous Media. Springer; Dordrecht, The Netherlands: 2006. p. 44.

3. Vadász P.T. Emerging Topics in Heat and Mass Transfer in Porous Media: From Bioengineering and Microelectronics to Nanotechnology. Springer; Dordrecht, The Netherlands: 2008. p. 328.

4. Verruijt A. An Introduction to Soil Dynamics. Springer; Dordrecht, The Netherlands: 2010. p. 433.

5. Bear J., Verruijt A. Modeling Groundwater Flow and Pollution: With Computer Programs for Sample Cases. D. Reidel Pub. Co.; Dordrecht, The Netherlands: 1987. p. 414.

6. Dhamrat R.S., Ellzey J.L. Numerical and experimental study of the conversion of methane to hydrogen in a porous media reactor. Combust. Flame. 2006;144:698–709. doi: 10.1016/j.combustflame.2005.08.038. [Cross Ref]

7. Pantangi V.K., Mishra S.C., Muthukumar P., Reddy R. Studies on porous radiant burners for lpg (liquefied petroleum gas) cooking applications. Energy. 2011;36:6074–6080. doi: 10.1016/j.energy.2011.08.008. [Cross Ref]

8. Zermatten E., Vetsch J.R., Ruffoni D., Hofmann S., Müller R., Steinfeld A. Micro-computed tomography based computational fluid dynamics for the determination of shear stresses in scaffolds within a perfusion bioreactor. Ann. Biomed. Eng. 2014;42:1085–1094. doi: 10.1007/s10439-014-0981-0. [PubMed] [Cross Ref]

9. Fend T., Hoffschmidt B., Pitz-Paal R., Reutter O., Rietbrock P. Porous materials as open volumetric solar receivers: Experimental determination of thermophysical and heat transfer properties. Energy. 2004;29:823–833. doi: 10.1016/S0360-5442(03)00188-9. [Cross Ref]

10. Hischier I., Leumann P., Steinfeld A. Experimental and numerical analyses of a pressurized air receiver for solar-driven gas turbines. J. Sol. Energy Eng. 2012;134:1–8.

11. Hischier I., Poživil P., Steinfeld A. A modular ceramic cavity-receiver for high-temperature high-concentration solar applications. J. Sol. Energy Eng. 2011;134 doi: 10.1115/1.4005107. [Cross Ref]

12. Romero M., Steinfeld A. Concentrating solar thermal power and thermochemical fuels. Energy Environ. Sci. 2012;5:9234–9245. doi: 10.1039/c2ee21275g. [Cross Ref]

13. Miller J.E., McDaniel A.H., Allendorf M.D. Considerations in the design of materials for solar-driven fuel production using metal-oxide thermochemical cycles. Adv. Energy Mater. 2014;4 doi: 10.1002/aenm.201300469. [Cross Ref]

14. Furler P., Scheffe J., Gorbar M., Moes L., Vogt U., Steinfeld A. Solar thermochemical CO_{2} splitting utilizing a reticulated porous ceria redox system. Energy Fuels. 2012;26:7051–7059.

15. Furler P., Scheffe J., Marxer D., Gorbar M., Bonk A., Vogt U., Steinfeld A. Thermochemical CO_{2} splitting via redox cycling of ceria reticulated foam structures with dual-scale porosities. PCCP. 2014;16:10503–10511. doi: 10.1039/c4cp01172d. [PubMed] [Cross Ref]

16. Furler P., Scheffe J.R., Steinfeld A. Syngas production by simultaneous splitting of H_{2}O and CO_{2} via ceria redox reactions in a high-temperature solar reactor. Energy Environ. Sci. 2012;5:6098–6103. doi: 10.1039/c1ee02620h. [Cross Ref]

17. Ackermann S., Scheffe J.R., Steinfeld A. Diffusion of oxygen in ceria at elevated temperatures and its application to H_{2}O/CO_{2} splitting thermochemical redox cycles. J. Phys. Chem. C. 2014;118:5216–5225. doi: 10.1021/jp500755t. [Cross Ref]

18. Chueh W.C., Falter C., Abbott M., Scipio D., Furler P., Haile S.M., Steinfeld A. High-flux solar-driven thermochemical dissociation of CO_{2} and H_{2}O using nonstoichiometric ceria. Science. 2010;330:1797–1801. doi: 10.1126/science.1197834. [PubMed] [Cross Ref]

19. Venstrom L.J., Petkovich N., Rudisill S., Stein A., Davidson J.H. The effects of morphology on the oxidation of ceria by water and carbon dioxide. J. Sol. Energy Eng. 2012;134 doi: 10.1115/1.4005119. [Cross Ref]

20. Gibbons W.T., Venstrom L.J., de Smith R.M., Davidson J.H., Jackson G.S. Ceria-based electrospun fibers for renewable fuel production via two-step thermal redox cycles for carbon dioxide splitting. Phys. Chem. Chem. Phys. 2014;16:14271–14280. doi: 10.1039/c4cp01974a. [PubMed] [Cross Ref]

21. Keene D.J., Davidson J.H., Lipiński W. A model of transient heat and mass transfer in a heterogeneous medium of ceria undergoing nonstoichiometric reduction. J. Heat Transfer. 2013;135 doi: 10.1115/1.4023494. [Cross Ref]

22. Kaviany M. Principles of Heat Transfer in Porous Media. 2nd ed. Springer-Verlag; New York, NY, USA: 1995. p. 708.

23. Whitaker S. The Method of Volume Averaging. Kluwer Academic; Dordrecht, The Netherlands: 1999. p. 219. (Series: Theory and Applications of Transport in Porous Media)).

24. Bodla K.K., Weibel J.A., Garimella S.V. Advances in fluid and thermal transport property analysis and design of sintered porous wick microstructures. J. Heat Transf. 2013;135 doi: 10.1115/1.4023569. [Cross Ref]

25. Krishnan S., Murthy J.Y., Garimella S.V. Direct simulation of transport in open-cell metal foam. J. Heat Transf. 2006;128:793–799. doi: 10.1115/1.2227038. [Cross Ref]

26. Petrasch J., Meier F., Friess H., Steinfeld A. Tomography based determination of permeability, dupuit-forchheimer coefficient, and interfacial heat transfer coefficient in reticulate porous ceramics. 2008;29:315–326.

27. Petrasch J., Wyss P., Stämpfli R., Steinfeld A. Tomography-based multiscale analyses of the 3D geometrical morphology of reticulated porous ceramics. J. Am. Ceram. Soc. 2008;91:2659–2665. doi: 10.1111/j.1551-2916.2008.02308.x. [Cross Ref]

28. Bodla K.K., Murthy J.Y., Garimella S.V. Microtomography-based simulation of transport through open-cell metal foams. Numer. Heat Transf. A Appl. 2010;58:527–544. doi: 10.1080/10407782.2010.511987. [Cross Ref]

29. Petrasch J., Wyss P., Steinfeld A. Tomography-based monte carlo determination of radiative properties of reticulate porous ceramics. J. Quant. Spectrosc. Radiat. Transf. 2007;105:180–197. doi: 10.1016/j.jqsrt.2006.11.002. [Cross Ref]

30. Haussener S., Coray P., Lipinski W., Wyss P., Steinfeld A. Tomography-based heat and mass transfer characterization of reticulate porous ceramics for high-temperature processing. J. Heat Transf. 2010;132 doi: 10.1115/1.4000226. [Cross Ref]

31. Petrasch J., Schrader B., Wyss P., Steinfeld A. Tomography-based determination of the effective thermal conductivity of fluid-saturated reticulate porous ceramics. ASME J. Heat Transf. 2008;130 doi: 10.1115/1.2804932. [Cross Ref]

32. Karl S., Somers A.V. Method of Making Porous Ceramic Articles. 3090094 A. U.S. Patent. 1963 May 21;

33. Otsu N. A threshold selection method from gray-level histograms. IEEE Trans. Syst. Man Cybern. 1979;9:62–66. doi: 10.1109/TSMC.1979.4310076. [Cross Ref]

34. Sezgin M., Sankur B.L. Survey over image thresholding techniques and quantitative performance evaluation. ELECTIM. 2004;13:146–168.

35. Zamel N., Li X., Shen J. Correlation for the effective gas diffusion coefficient in carbon paper diffusion media. Energy Fuels. 2009;23:6070–6078. doi: 10.1021/ef900653x. [Cross Ref]

36. Giesche H. Mercury porosimetry: A general (practical) overview. Part. Part. Syst. Charact. 2006;23:9–19. doi: 10.1002/ppsc.200601009. [Cross Ref]

37. Webb P.A. An Introduction to the Physical Characterization of Materials by Mercury Intrusion Porosimetry with Emphasis on Reduction and Presentation of Experimental Data. Micromeritics Instrument Corp.; Norcross, GA, USA: 2001.

38. Abràmoff M.D., Magalhães P.J., Ram S.J. Image processing with imageJ. Biophoton. Int. 2004;11:36–43.

39. Doube M., Kłosowski M.M., Arganda-Carreras I., Cordelières F.P., Dougherty R.P., Jackson J.S., Schmid B., Hutchinson J.R., Shefelbine S.J. Bonej: Free and extensible bone image analysis in imageJ. Bone. 2010;47:1076–1079. doi: 10.1016/j.bone.2010.08.023. [PMC free article] [PubMed] [Cross Ref]

40. Friess H., Haussener S., Steinfeld A., Petrasch J. Tetrahedral mesh generation based on space indicator functions. Int. J. Numer. Methods Eng. 2013;93:1040–1056. doi: 10.1002/nme.4419. [Cross Ref]

41. Berryman J.G., Blair S.C. Use of digital image analysis to estimate fluid permeability of porous materials: Application of two-point correlation functions. J. Appl. Phys. 1986;60:1930–1938. doi: 10.1063/1.337245. [Cross Ref]

42. Mogensen M., Sammes N.M., Tompsett G.A. Physical, chemical and electrochemical properties of pure and doped ceria. Solid State Ion. 2000;129:63–94. doi: 10.1016/S0167-2738(99)00318-5. [Cross Ref]

43. Bruggeman D.A.G. Berechnung verschiedener physikalischer konstanten von heterogenen substanzen. I. Dielektrizitätskonstanten und leitfähigkeiten der mischkörper aus isotropen substanzen. Ann. Phys. 1935;416:636–664. doi: 10.1002/andp.19354160705. (In German) [Cross Ref]

44. Deissler R.G., Boegli J.S. An investigation of effective thermal conductivities of powders in various gases. Trans. ASME. 1958;80:1417–1425.

45. Wang M., Pan N. Predictions of effective physical properties of complex multiphase materials. Mater. Sci. Eng. R Rep. 2008;63:1–30. doi: 10.1016/j.mser.2008.07.001. [Cross Ref]

46. DeVera A.L., Strieder W. Upper and lower bounds on the thermal conductivity of a random, two-phase material. J. Phys. Chem. 1977;81:1783–1790. doi: 10.1021/j100533a014. [Cross Ref]

47. Hashin Z., Shtrikman S. A variational approach to the theory of the effective magnetic permeability of multiphase materials. J. Appl. Phys. 1962;33:3125–3131. doi: 10.1063/1.1728579. [Cross Ref]

48. Woodside W., Messmer J.H. Thermal conductivity of porous media. II. Consolidated rocks. J. Appl. Phys. 1961;32:1699–1706. doi: 10.1063/1.1728420. [Cross Ref]

49. Russell H.W. Principles of heat flow in porous insulators. J. Am. Ceram. Soc. 1935;18:1–5. doi: 10.1111/j.1151-2916.1935.tb19340.x. [Cross Ref]

50. Loeb A.L. Thermal conductivity: VIII, a theory of thermal conductivity of porous materials. J. Am. Ceram. Soc. 1954;37:96–99. doi: 10.1111/j.1551-2916.1954.tb20107.x. [Cross Ref]

51. Maxwell J.C. A Treatise on Electricity and Magnetism. 3rd ed. Dover Publications; New York, NY, USA: 1954.

52. Choy T.C. Effective Medium Theory: Principles and Applications. Oxford University Press; Oxford, UK: 1999. p. 182.

53. Sattler K.D. Handbook of Nanophysics. Nanoparticles and Quantum Dots. Taylor & Francis; Boca Raton, FL, USA: 2011.

54. Kamiuto K. Cellular and Porous Material Thermal Properties Simulation Prediction. Wiley-VCH; Weinheim, Germany: 2008. Modeling of composite heat transfer in open-cellular porous materials at high temperatures; pp. 165–198.

55. Coquard R., Loretz M., Baillis D. Conductive heat transfer in metallic/ceramic open-cell foams. Adv. Eng. Mater. 2008;10:323–337. doi: 10.1002/adem.200700331. [Cross Ref]

56. Bhattacharya A., Calmidi V.V., Mahajan R.L. An analytical-experimental study for the determination of the effective thermal conductivity of high porosity fibrous foams. ASME Appl. Mech. Divis. Pupl. AMD. 1999;233:13–20.

57. Boomsma K., Poulikakos D. On the effective thermal conductivity of a three-dimensionally structured fluid-saturated metal foam. Int. J. Heat Mass Transf. 2001;44:827–836. doi: 10.1016/S0017-9310(00)00123-X. [Cross Ref]

58. Hamilton R.L., Crosser O.K. Thermal conductivity of heterogeneous two-component systems. Ind. Eng. Chem. Fundam. 1962;1:187–191. doi: 10.1021/i160003a005. [Cross Ref]

59. Miller M.N. Bounds for effective electrical, thermal, and magnetic properties of heterogeneous materials. J. Math. Phys. 1969;10:1988–2004. doi: 10.1063/1.1664794. [Cross Ref]

60. Calmidi V.V., Mahajan R.L. The effective thermal conductivity of high porosity fibrous metal foams. J. Heat Transf. 1999;121:466–471. doi: 10.1115/1.2826001. [Cross Ref]

61. Matsushita M., Monde M., Mitsutake Y. Predictive calculation of the effective thermal conductivity in a metal hydride packed bed. Int. J. Hydrog. Energy. 2014;39:9718–9725. doi: 10.1016/j.ijhydene.2014.04.072. [Cross Ref]

62. Dul’Nev G.N., Zarichnyak Y.P. A study of the generalized conductivity coefficients in heterogeneous systems. Heat Transf. Sov. Res. 1970;2:89–107.

Articles from Materials are provided here courtesy of **Multidisciplinary Digital Publishing Institute (MDPI)**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |