Cellular Ceramics / p3
.1.pdf
3.1 Characterization of Structure and Morphology 245
tion of the equivalent pore diameter for samples 1, 3, 5, and 7; the mean pore diameter increases with increasing pore size of the original polyurethane foam.
Area number density |
|
|
5 |
|
7 |
|
|
|
|
|
|
|
|
||
1 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
200 |
|
400 |
600 |
800 |
1000 |
1200 |
1400 |
Pore diameter / µm
Figure 15 Area number density by differentiation of the cumulative curves as a function of the equivalent pore diameter for samples 1, 3, 5 and 7.
Assuming that the investigated samples are characterized by a random orientation of spherical pores with a relatively narrow cell size distribution, a simple stereological correction factor can be applied to determine the true three-dimensional mean pore diameter. It can be shown that the average sphere diameter Dsphere is larger than the average circular segment diameter Dcirc due to random truncation of the cells with respect to depth at the plane of the specimen surface. The relation can be expressed by the following equation [20, 22]:
Dcirc |
|
Dsphere ¼ 0:785 . |
(7) |
Figure 16 compares the mean pore diameters measured by the Visiocell method and the image analysis results. When a uniform distribution of spheres is assumed, the largest diameter in the pore size distribution corresponds to a section in the center of the sphere. Because of the limited number of measurements in this bin size, the d90 value of the cumulative curve is regarded as the largest diameter.
Applying the simple stereological correction factor from Eq. (7) results in a mean pore diameter that is in good agreement with the other methods. Differences between the d90 value and the corrected d50 value range between 1 and 7 %. This indicates that the basic assumptions for Eq. (7) are valid for these kinds of samples. For most samples, the difference between the d90 value from image analysis and the Visiocell method range between 3 and 5 % (except for samples 1 and 6, with differences of 9 and 17 %, respectively).
Alternatively, the two-dimensional slices from micro computer tomography (l-CT) can be used as input for the image analysis routine. As the difference in linear attenuation coefficient between air and ceramic material is fairly large, the image contrast is high. Figure 17a shows a slice of a l-CT scan of sample 8. The corresponding identification of the individual pores by the image analysis routine is presented in Fig. 17b.
246 |
Part 3 Structure |
|
|
|
|
|
|
|
|
|
µm |
1600 |
Visiocell |
|
|
|
|
|
|
|
1400 |
|
|
|
|
|
|||
|
Image analysis (d90) |
|
|
|
|||||
|
/ |
1200 |
Image analysis (d50/0.785) |
|
|
||||
|
diameter |
Image analysis (d50) |
|
|
|
||||
|
1000 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
800 |
|
|
|
|
|
|
|
|
Pore |
600 |
|
|
|
|
|
|
|
|
400 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
200 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
Sample number
Figure 16 Comparison of the equivalent pore diameter for the eight investigated samples, measured by the Visiocell method and image analysis: the d90 value, three-dimensional calculation according to Eq (7) (d50/0.785), and d50 value.
(a) |
(b) |
Figure 17 a) Slice by computerized X-ray microtomography of sample 8 and b) identification of the pore diameter in the imageanalysis routine.
Image analysis was performed on l-CT slices of samples 5, 6, 7, and 8. The measured data are the average of analyses on two randomly chosen l-CT slices. This corresponds to about 500 analyzed pores. Figure 18 compares the d90 and d50 values in the cumulative curves for analysis of images obtained by light microscopy and slices obtained by tomography. The l-CT-derived data are consistently lower for all investigated samples, both for the d90 and the d50 value. The differences in d90 value are about 4–8 % (17 % for sample 8). This might be due to difference in setting the threshold of the optical image and the thresholding in the microtomography procedure.
3)Shape factor:
The shape of the pores is of interest in many applications. Moreover, if stereological corrections are applied, information on the shape of the features is in most cases necessary. Depending on the expected form of the intersection (circular, elliptical, triangular, or rectangular), the aspect ratio and the degree of roundness must be
3.1 Characterization of Structure and Morphology 247
estimated. When performing image analysis on thin slices that are cut in x, y, and z directions, the differences in intersection forms provide information on the anisotropy and the general shape of the pore channel.
|
1600 |
Visiocell |
|
|
|
|
|
|
|
|
|
|
1400 |
IA on microscopy (d90) |
|
|
|
|
1200 |
IA on µ-CT slice (d90) |
|
|
|
/µm |
IA on microscopy (d50) |
|
|
||
1000 |
IA on µ-CT slice (d50) |
|
|
||
size |
800 |
|
|
|
|
|
|
|
|
|
|
Pore |
600 |
|
|
|
|
400 |
|
|
|
|
|
|
|
|
|
|
|
|
200 |
|
|
|
|
|
0 |
5 |
6 |
7 |
8 |
|
|
||||
Sample number
Figure 18 Comparison of the pore size determined by the Visiocell method, d90 values from image analysis on optical microscopy images and on microtomography (l-CT) slices and d50 values from image analysis on microscopy optical images and on l-CT slices.
As these samples can be described by a sphere-like pore model, a circularity factor of each pore in the sample can be calculated. The circularity factor Fcirc represents
the roundness of the region and can be defined as |
|
|
area |
. |
|
Fcirc ¼ 4p perimeter2 |
(8) |
|
For a perfect circle, the shape factor is 1. More elongated polygonal intersections lower the circularity factor. Figure 19 presents the shape factor for all pores in a micrograph of sample 1 (total number of pores 167). The mean circularity for this image is 0.74, with a standard deviation of 0.1.
500 |
|
|
|
1,0 |
|
400 |
|
|
|
0,8 |
factor Circularity |
|
|
|
|
||
|
|
|
|
0,6 |
|
300 |
|
|
|
0,4 |
|
200 |
|
|
|
||
|
|
|
0,2 |
||
diameterPoreµm/ |
|
|
|
|
|
|
|
|
|
0,0 |
|
0 |
40 |
80 |
120 |
160 |
|
Pore number
(a) |
(b) |
Figure 19 a) Identification of pores by image analysis on microscopic image of ceramic foam and b) ascending list of the equivalent diameter of each pore in the image of sample 1 and the corresponding circularity factor.
248 Part 3 Structure
3.1.2.2.4Capillary Flow Porometry
The structure of porous materials is often very complex because of the nature of the porosity (closed porosity, blind pores, and through-pores), the pore shape, the surface area, and the isotropy. For applications like filter media, the fluid or gas flow characteristics, the barrier properties, and the efficiency of the process are governed by the combination of all aspects of the pore structure.
Capillary flow porometry, also known as extrusion flow porometry, can be applied to determine some important characteristics of porous materials [38, 39]:
. |
Distribution of the constricted part of a pore channel. |
. |
Flow distribution curve (bubble point, mean diameter). |
. |
Liquid and gas permeability. |
. |
Surface area. |
Principle
The sample is placed in a wetting liquid to fill the pores of the sample. As the liquid/solid surface free energy cls is lower than the solid/gas surface free energy csg, pores are spontaneously filled by the liquid. When a nonreacting gas increasingly pressurises the sample, the liquid is removed from the pores and a gas flow through the pores results. To displace the liquid from the pores, the work done by the gas must be equal to the increase in surface free energy. This can be expressed by the following equation:
p dV ¼ ðcsg cls Þ dS |
(9) |
where p is differential pressure, dV the increase in volume of gas in the pore, and dS the increase in solid/gas surface area (or the corresponding decrease in solid/ liquid interfacial area).
The pore cross section can be quite complex. The diameter of a pore at any location along the pore path can be defined as the diameter of a cylindrical opening that has the same dS/dV ratio as the actual pore. This is illustrated in Fig. 20.
pore cross section
dS |
|
|
|
dV
dS |
4 |
|||
pore = |
|
|
cylindrical opening with diameter D = |
|
|
|
|||
dV |
|
D |
||
rectangular pore |
measured diameter |
12 x 6 |
8 |
Figure 20 Schematic representation of pore cross section and the definition of pore diameter in capillary flow analysis.
Taking into account the relation between the contact angle h and surface tension c, the relation between the pore diameter D and the differential pressure needed to displace the liquid from the pore can be expressed as:
p ¼ 4c |
cos h |
. |
(10) |
D |
3.1 Characterization of Structure and Morphology 249
For exact measurements, information on the contact angle is necessary. For low surface tension wetting fluids the contact angle approaches zero. Figure 21 gives a schematic representation of the principle of this method. When the pressure is gradually increased, the liquid is removed first from the pores with the largest diameter, and a gas flow through the pore results. Higher pressures are needed to empty the smaller pores.
The presence of pores is detected by measuring the flow rate at a given applied differential pressure with flow meters. The flow curve as a function of the differential pressure is measured for both the wet and the dry sample.
wetting liquid
sample flow
gas under pressure
Figure 21 Schematic view of the principle of capillary flow analysis.
Interpretation of measured characteristics
The measured differential pressures and flow rates are used to calculate a number of pore characteristics of the sample.
1) The measurement of the constricted pore diameter distribution is based on the complete removal of a wetting liquid from a pore with a certain diameter at a distinct differential pressure; it can be calculated from Eq. (10). Figure 22 shows the schematic view of a through-pore channel. To measure a gas flow through the pore, the differential pressure must be large enough to remove the liquid from the most constricted part of the pore. Until this pressure is reached, the gas displaces liquid in the pore up to the smallest diameter in the pore channel. Thus, the differential pressure that permits flow through a pore corresponds to the displacement of the liquid at the most constricted part of the pore. The flow meter detects the pore by sensing the increase in flow rate at a certain differential pressure. The calculated pore diameter is then the diameter of the through-pore at its most constricted part.
The bubble point, being the largest constricted through-pore diameter, corresponds to the lowest pressure at which gas flow through the wet sample is detected. The mean constricted through pore diameter is calculated by the intersection of the wet flow curve and the half-dry flow curve. This mean pore diameter corresponds to half of the flow through pores larger than the mean value. A flow distribution curve as a function of
the pore diameter can be calculated as well in terms of the function f:
!
f ¼ |
dfw |
100=dD |
(11) |
dfd |
250 Part 3 Structure
where fw and fd are the flow rates through wet and dry samples, respectively, at the same differential pressure. The area under the curve in a certain pore range is the percentage flow through that pore size range. A typical result of a capillary flow analysis of a ceramic foam is presented in Fig. 23.
Pore
Liquid
D2
Pore diameters |
|
|
measured by |
|
|
mercury intrusion |
|
|
porosimetry |
D1 |
Constricted pore |
|
||
|
|
|
Pore diameters |
|
diameter |
|
measured by |
|
measured by |
|
|
|
capillary flow |
|
liquid extrusion |
|
|
|
porometry |
|
porosimetry |
|
|
D2 |
|
|
|
|
|
|
Gas pressure |
|
Figure 22 Schematic view of a through-pore channel and the pore diameter measured by capillary flow analysis (in comparison with liquid extrusion and mercury intrusion porosimetry).
Flow rate / L min–1
80


dry flow curve


wet flow curve
60 


half-dry flow curve
40
20
0 |
|
|
|
|
0,00 |
0,05 |
0,10 |
0,15 |
0,20 |
Pressure / bar
Figure 23 Capillary flow analysis of sample 1, showing the dry flow curve, the wet flow curve, and the half-dry flow curve.
2) The gas flow curve through the dry sample permits the calculation of the gas permeability using Darcy’s law [40]:
F ¼ k |
A |
|
ðpi p0 Þ |
(12) |
|
l |
|
l |
|||
|
|
|
|
|
|
where F is the volume flow rate at the average pressure per unit time, k the permeability, l the viscosity of the fluid, A the cross-sectional area of the porous sample, l the thickness of the porous sample, pi the inlet pressure, and po the outlet pressure.
Analogously, liquid permeability for a variety of liquids and solutions can be measured by using the corresponding wet flow curve.
3.1 Characterization of Structure and Morphology 251
3) The flow rate data enable calculation of the surface area by using the Kozeny– Carman relation [41, 42]. Accurate measurements are possible when the specific surface area is less than 10 m2 g–1 and the distribution of pore diameters is relatively narrow. The agreement with gas adsorption measurements is reasonably good if the sample is free of blind pores.
Specifications
In this study, a PMI Capillary flow porometer, type CFP-1200-A was used (Porous Materials Inc., USA). The specifications according to the manufacturer are summarized in Tab. 3. The measurable pore size stretches over a broad range from 13 nm to 500 mm. The low pressure needed to remove the liquid enables the characterization of mechanically weaker structures (e.g., polymer foams) and avoids distortion of the pore structures.
Table 3 Specification of the capillary flow analysis equipment.
Property |
Specification |
|
|
Pore size range |
13–500 mm |
Sample diameter |
1–6 cm |
Pressure range |
0–35 bar |
Pressurizing gas |
dry, compressed air, nonflammable, noncorrosive |
Pressure transducer range |
0–35 bar |
Resolution |
1 in 20 000 |
Accuracy |
0.15 % of reading |
Mass flow transducer range |
10 ml min–1 to 500 l min–1 |
Results
Figure 24a presents pore diameter and filter flow (ratio of the wet and dry flow curves) as a function of the applied pressure for sample 4. From these raw data, the bubble point can easily be identified as the pore diameter that corresponds to the
µm/diameterPore |
|
Bubble point |
|
|
|
100 |
distributionsizePore |
|
|
|
|
|
|
|
|
300 |
|
|
|
|
|
%/flowFilter |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pore diameter |
80 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Filter flow |
|
|
|
|
|
|
|
|
|
|
|
|
200 |
|
|
|
|
60 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
|
|
|
|
100 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
0,00 |
0,05 |
0,10 |
0,15 |
0,20 |
|
|
0 |
50 |
100 |
150 |
200 |
250 |
300 |
350 |
Differential Pressure / bar |
Average diameter / µm |
(a) |
(b) |
Figure 24 a) Relation between pore diameter and filter flow as a function of differential pressure and b) size distribution of the most restricted part in the pore channels for sample 4.
252 Part 3 Structure
lowest pressure needed to initiate flow through the wet sample. For this sample, the bubble point is 338 mm, and 90 % of the filter flow occurs for pores larger than 49 mm. A pore size distribution is obtained by the differentiation of the filter flow as a function of the pore diameter (Fig. 24b).
Figure 25 presents an overview of the filter flow curves for the eight samples. Although the pore diameters of the different samples vary over a relatively broad range (from 200 to 1700 mm), the pore sizes as measured by capillary flow analysis are distributed over a much smaller range (ca. 20 to 700 mm). This can be explained by the presence of semiclosed windows and other restricted pore channels that cannot be detected by image analysis or other techniques based on microscopy.
Filter flow / %
100 |
|
|
|
|
|
|
|
80 |
|
|
|
|
|
|
|
60 |
5 |
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
|
20 |
1 |
2 |
|
|
|
|
|
|
3 |
4 |
7 |
8 |
|
|
|
|
|
|
|
||||
0 |
|
|
|
|
|
|
|
0 |
150 |
300 |
|
450 |
600 |
750 |
|
Diameter / µm
Figure 25 Filter flow curves measured by capillary flow analysis for the eight samples.
Comparing the information of the capillary flow analysis with the pore diameter determined by the Visiocell measurements clearly shows the difference between the two techniques regarding the measured property (Fig. 26). As capillary flow analysis only measures the most restricted part of the pore channel, mean pore flow diameter and the bubble point are lower than the average diameter of the pore size determined by the Visiocell technique.
Although a general increase in the restricted pore diameter is seen with increasing cell size of the polyurethane foam, this trend is not followed to the same degree. When coating with more viscous suspensions, the formation of semiclosed cell windows cannot be avoided, as exemplified by the lowering of the restricted pore channel diameter for samples 6 and 8 compared to samples 5 and 7, respectively.
The mean flow pore diameter and the bubble point are critical parameters of porous materials that are used in applications such as filtration or porous bone replacements. For the latter, bone ingrowth requires a minimum pore channel diameter that enables migrating biological material like cells or veins to penetrate the porous implant. The combination of the mean pore diameter and the most restricted part of the pore channel play a critical role in the ability for complete ingrowth.
3.1 Characterization of Structure and Morphology 253
Diameter / µm
1500 |
|
|
|
|
|
|
|
|
Visiocell |
|
|
|
|
|
|
1200 |
Bubblepoint |
|
|
|
|
|
|
|
Mean flow pore |
|
|
|
|
||
900 |
|
|
|
|
|
|
|
600 |
|
|
|
|
|
|
|
300 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
Sample number
Figure 26 Comparison of the pore diameter by the Visiocell method, bubble point, and mean flow pore diameter measured by capillary flow analysis.
3.1.2.2.5Micro Computer Tomography
Knowledge of the structure and morphology of highly porous ceramic materials is essential for various applications. Structural properties such as fractional density, which corresponds to the overall porosity, the mean size and size distribution of the foam cells and the surrounding struts, the shape of the cells, and the isotropy of the individual cells influence many properties. For example, gas permeability is increased by large windows interconnecting cells and decreased by small cell diameters.
Several approaches have been used to evaluate the structure, such as magnetic resonance imaging and optical analysis. Tomography facilitates three-dimensional imaging of complex structures. Optical tomography was used to evaluate transparent aqueous foams: light was transmitted through a sample and detected by a camera with a plane detector [43, 44]. The drawback of light optical analysis is that it cannot be used for nondestructive evaluation of the three-dimensional structure of nontransparent ceramics. However, X rays transmit through most ceramic materials, and thus X-ray tomography can be used to characterize porous ceramic materials [45].
Principle
X-ray micro tomography is based on the evaluation of stacked two-dimensional layers of the sample that are scanned sequentially [46] (Fig. 27a). An X-ray source, usually a microfocus X-ray tube, emits a fan-shaped beam of intensity I0 that passes through the specimen. The intensity I is detected by a line array of length x. The relative intensity I/I0 of the transmitted beam depends on the amount and nature of the material being analyzed. Between two recordings the sample is slightly rotated. After a full rotation up to several hundred recordings have been taken. The intensity data can be displayed as a sinugram (Fig. 27b).
254 Part 3 |
Structure |
|
I(L, ) |
|
Measured Data: |
|
|
Sinugram (single layer) |
|
L |
0˚ |
|
|
|
|
y |
|
|
x |
|
Usually up |
x-ray source |
to 500 |
|
layers |
180˚ |
|
|
(a) |
(b) |
Figure 27 a) Measurement schematic of X-ray tomography and b) raw intensity data displayed as sinugram.
The intensities I and I0 with and without material and can be described by a line integral
R |
lðx; yÞ dx dy! |
|
I ¼ I0 exp L |
(13) |
where l is the material-dependent X-ray absorption coefficient and x is the position on the line L between the X-ray source and the detector; l(x,y) is called the object function. During the measurement the sample rotates around the center axis with an angle h
s |
|
cosh |
sinh |
|
x |
|
(14) |
t |
¼ sinh |
cosh |
y |
||||
|
|
|
|
|
|
|
|
where s and t are the coordinates of the rotating system, and x,y the coordinates of the fixed system. The X-ray beam is parallel to the s-axis in the fixed system. If t0 is the distance between the beam and the s-axis, the points that are on the beam can be described according to
x sinh þ y cosh ¼ t0 . |
(15) |
Equation (13) can be transformed using the Dirac d(x), which has the value of infinity for x = 0, zero elsewhere, and a total integral of one. Equation (15) can be described as a line integral after normalization by I0 and inverting:
þ ¥ þ ¥ |
I0 |
|
l ðh; t0 Þ ¼ Ð¥ Ð¥ |
lðx; yÞ dðx sinh þ y cosh t0 Þ dx dy ¼ ln Iðh; t0 Þ. |
(16) |
This projection of l(x,y) on the line integrals l(h,t0) is equivalent to the Radon transformation. A projection Ph is defined as the line integrals (–¥ < t < +¥) for a
