Freeze Casting

Freeze-casting, also frequently referred to as ''ice-templating'', ''freeze casting'', or ''freeze alignment'', is a technique that exploits the highly anisotropic solidification behavior of a solvent (generally water) in a well-dispersed solution or slurry to controllably template directionally porous ceramics, polymers, metals and their hybrids. By subjecting an aqueous solution or slurry to a directional temperature gradient, ice crystals will nucleate on one side and grow along the temperature gradient. The ice crystals will redistribute the dissolved substance and the suspended particles as they grow within the solution or slurry, effectively templating the ingredients that are distributed in the solution or slurry.

Once solidification has ended, the frozen, templated composite is placed into a freeze-dryer to remove the ice crystals. The resulting green body contains anisotropic macropores in a replica of the sublimated ice crystals and structures from micropores to nacre-like packing between the ceramic or metal particles in the walls. The walls templated by the morphology of the ice crystals often show unilateral features. These together build a hierarchically structured cellular structure. This structure is often sintered for metals and ceramics, and crosslinked for polymers, to consolidate the particulate walls and provide strength to the porous material. The porosity left by the sublimation of solvent crystals is typically between 2–200 μm.

Overview

The first observation of cellular structures resulting from the freezing of water goes back over a century, but the first reported instance of freeze-casting, in the modern sense, was in 1954 when Maxwell et al. attempted to fabricate turbosupercharger blades out of refractory powders. They froze extremely thick slips of titanium carbide, producing near-net-shape castings that were easy to sinter and machine. The goal of this work, however, was to make dense ceramics. It was not until 2001, when Fukasawa et al. created directionally porous alumina castings, that the idea of using freeze-casting as a means of creating novel porous structures really took hold. Since that time, research has grown considerably with hundreds of papers coming out within the last decade.
Ice templating, freeze casting: Beyond materials processing

The principles of freeze casting are applicable to a broad range of combinations of particles and suspension media. Water is by far the most commonly used suspension media, and by freeze drying is readily conducive to the step of sublimation that is necessary for the success of freeze-casting processes. Due to the high level of control and broad range of possible porous microstructures that freeze-casting can produce, the technique has been adopted in disparate fields such as tissue engineering, tissue scaffolds, photonics, metal-matrix composites, dentistry,
, Freeze Casting of High Strength Composites for Dental Applications materials science,
, Dispersion, connectivity and tortuosity of hierarchical porosity composite SOFC cathodes prepared by freeze-casting
, Processing of Hierarchical and Anisotropic LSM-YSZ Ceramics
Lightweight and stiff cellular ceramic structures by ice templating and even food science.

There are three possible end results to uni-directionally freezing a suspension of particles. First, the ice-growth proceeds as a planar front, pushing particles in front like a bulldozer pushes a pile of rocks. This scenario usually occurs at very low solidification velocities (< 1 μm s⁻¹) or with extremely fine particles because they can move by Brownian motion away from the front. The resultant structure contains no macroporosity. If one were to increase the solidification speed, the size of the particles or solid loading moderately, the particles begin to interact in a meaningful way with the approaching ice front. The result is typically a lamellar or cellular templated structure whose exact morphology depends on the particular conditions of the system. It is this type of solidification that is targeted for porous materials made by freeze-casting. The third possibility for a freeze-cast structure occurs when particles are given insufficient time to segregate from the suspension, resulting in complete encapsulation of the particles within the ice front. This occurs when the freezing rates are rapid, particle size becomes sufficiently large, or when the solids loading is high enough to hinder particle motion.
To ensure templating, the particles must be ejected from the oncoming front. Energetically speaking, this will occur if there is an overall increase in free energy if the particle were to be engulfed ''(Δσ > 0)''.

$\Delta\sigma = \sigma_ - (\sigma_ + \sigma_)$

where ''Δσ'' is the change in free energy of the particle, σ_ps is the surface potential between the particle and interface, ''σ_pl'' is the potential between the particle and the liquid phase and ''σ_sl'' is the surface potential between the solid and liquid phases. This expression is valid at low solidification velocities, when the system is shifted only slightly from equilibrium. At high solidification velocities, kinetics must also be taken into consideration. There will be a liquid film between the front and particle to maintain constant transport of the molecules which are incorporated into the growing crystal. When the front velocity increases, this film thickness ''(d)'' will decrease due to increasing drag forces. A critical velocity ''(v_c)'' occurs when the film is no longer thick enough to supply the needed molecular supply. At this speed the particle will be engulfed. Most authors express v_c as a function of particle size where $v_c \propto \tfrac$. The transition from a porous R (lamellar) morphology to one where the majority of particles are entrapped occurs at ''v_c'', which is generally determined as:

$v_c = \frac \left (\frac \right)^z$

where ''a₀'' is the average intermolecular distance of the molecule that is freezing within the liquid, ''d'' is the overall thickness of the liquid film, ''η'' is the solution viscosity, ''R'' is the particle radius and ''z'' is an exponent that can vary from 1 to 5. As expected, ''v_c'' decreases as particle radius ''R'' goes up.

Waschkies et al. studied the structure of dilute to concentrated freeze-casts from low (< 1 μm s⁻¹) to extremely high (> 700 μm s⁻¹) solidification velocities. From this study, they were able to generate morphological maps for freeze-cast structures made under various conditions. Maps such as these are excellent for showing general trends, but they are quite specific to the materials system from which they were derived. For most applications where freeze-casts will be used after freezing, binders are needed to supply strength in the green state. The addition of binder can significantly alter the chemistry within the frozen environment, depressing the freezing point and hampering particle motion leading to particle entrapment at speeds far below the predicted ''v_c''. Assuming, however, that we are operating at speeds below v_c and above those which produce a planar front, we will achieve some cellular structure with both ice-crystals and walls composed of packed ceramic particles. The morphology of this structure is tied to some variables, but the most influential is the temperature gradient as a function of time and distance along the freezing direction.

Freeze-cast structures have at least three apparent morphological regions.
Morphological instability in freezing colloidal suspensions At the side where freezing initiates is a nearly isotropic region with no visible macropores dubbed the Initial Zone (IZ). Directly after the IZ is the Transition Zone (TZ), where macropores begin to form and align with one another. The pores in this region may appear randomly oriented. The third zone is called the Steady-State Zone (SSZ), macropores in this region are aligned with one another and grow in a regular fashion. Within the SSZ, the structure is defined by a value λ that is the average thickness of a ceramic wall and its adjacent macropore.

Initial zone: nucleation and growth mechanisms

Although the ability of ice to reject suspended particles in the growth process has long been known, the mechanism remains the subject of some discussion. It was believed initially that during the moments immediately following the ice nucleus, nucleation of the ice crystals, particles are rejected from the growing planar ice front, leading to the formation of a constitutionally super-cooled zone directly ahead of the growing ice. This unstable region eventually results in perturbations, breaking the planar front into a columnar ice front, a phenomenon better known as a Mullins-Serkerka instability. After the breakdown, the ice crystals grow along the temperature gradient, pushing ceramic particles from the liquid phase aside so that they accumulate between the growing ice crystals. However, recent in-situ X-ray radiography of directionally frozen alumina suspensions reveal a different mechanism.

Transition zone: a changing microstructure

As solidification slows and growth kinetics become rate-limiting, the ice crystals begin to exclude the particles, redistributing them within the suspension. A competitive growth process develops between two crystal populations, those with their basal planes aligned with the thermal gradient (z-crystals) and those that are randomly oriented (r-crystals) giving rise to the start of the TZ.

There are colonies of similarly aligned ice crystals growing throughout the suspension. There are fine lamellae of aligned z-crystals growing with their basal planes aligned with the thermal gradient. The r-crystals appear in this cross-section as platelets but in actuality, they are most similar to columnar dendritic crystals cut along a bias. Within the transition zone, the r-crystals either stop growing or turn into z-crystals that eventually become the predominant orientation, and lead to steady-state growth. There are some reasons why this occurs. For one, during freezing, the growing crystals tend to align with the temperature gradient, as this is the lowest energy configuration and thermodynamically preferential. Aligned growth, however, can mean two different things. Assuming the temperature gradient is vertical, the growing crystal will either be parallel (z-crystal) or perpendicular (r-crystal) to this gradient. A crystal that lays horizontally can still grow in line with the temperature gradient, but it will mean growing on its face rather than its edge. Since the thermal conductivity of ice is so small (1.6 - 2.4 W mK⁻¹) compared with most every other ceramic (ex. Al₂O₃= 40 W mK⁻¹), the growing ice will have a significant insulative effect on the localized thermal conditions within the slurry. This can be illustrated using simple resistor elements.

When ice crystals are aligned with their basal planes parallel to the temperature gradient (z-crystals), they can be represented as two resistors in parallel. The thermal resistance of the ceramic is significantly smaller than that of the ice however, so the apparent resistance can be expressed as the lower R_ceramic. If the ice crystals are aligned perpendicular to the temperature gradient (r-crystals), they can be approximated as two resistor elements in series. For this case, the R_ice is limiting and will dictate the localized thermal conditions. The lower thermal resistance for the z-crystal case leads to lower temperatures and greater heat flux at the growing crystals tips, driving further growth in this direction while, at the same time, the large R_ice value hinders the growth of the r-crystals. Each ice crystal growing within the slurry will be some combination of these two scenarios. Thermodynamics dictate that all crystals will tend to align with the preferential temperature gradient causing r-crystals to eventually give way to z-crystals, which can be seen from the following radiographs taken within the TZ.

When z-crystals become the only significant crystal orientation present, the ice-front grows in a steady-state manner except there are no significant changes to the system conditions. It was observed in 2012 that, in the initial moments of freezing, there are dendritic r-crystals that grow 5 - 15 times faster than the solidifying front. These shoot up into the suspension ahead of the main ice front and partially melt back. These crystals stop growing at the point where the TZ will eventually fully transition to the SSZ. Researchers determined that this particular point marks the position where the suspension is in an equilibrium state (i.e. freezing temperature and suspension temperature are equal). We can say then that the size of the initial and transition zones are controlled by the extent of supercooling beyond the already low freezing temperature. If the freeze-casting setup is controlled so that nucleation is favored at only small supercooling, then the TZ will give way to the SSZ sooner.

Steady-state growth zone

The structure in this final region contains long, aligned lamellae that alternate between ice crystals and ceramic walls. The faster a sample is frozen, the finer its solvent crystals (and its eventual macroporosity) will be. Within the SSZ, the normal speeds which are usable for colloidal templating are 10 – 100 mm s⁻¹ leading to solvent crystals typically between 2 mm and 200 mm. Subsequent sublimation of the ice within the SSZ yields a green ceramic preform with porosity in a nearly exact replica of these ice crystals. The microstructure of a freeze-cast within the SSZ is defined by its wavelength ''(λ)'' which is the average thickness of a single ceramic wall plus its adjacent macropore. Several publications have reported the effects of solidification kinetics on the microstructures of freeze-cast materials. It has been shown that ''λ'' follows an empirical power-law relationship with solidification velocity ''(υ)'' (Eq. 2.14):

$\lambda = A\nu^$

Both ''A'' and ''υ'' are used as fitting parameters as currently there is no way of calculating them from first principles, although it is generally believed that ''A'' is related to slurry parameters like viscosity and solid loading while ''n'' is influenced by particle characteristics.

Controlling the porous structure

There are two general categories of tools for architecture a freeze-cast:
# Chemistry of the System - freezing medium and chosen particulate material(s), any additional binders, dispersants or additives.
# Operational Conditions - temperature profile, atmosphere, mold material, freezing surface, etc.

Initially, the materials system is chosen based on what sort of final structure is needed. This review has focused on water as the vehicle for freezing, but there are some other solvents that may be used. Notably, camphene, which is an organic solvent that is waxy at room temperature. Freezing of this solution produces highly branched dendritic crystals. Once the materials system is settled on however, the majority of microstructural control comes from external operational conditions such as mold material and temperature gradient.

Controlling pore size

The microstructural wavelength (average pore + wall thickness) can be described as a function of the solidification velocity v (λ= Av⁻ⁿ) where ''A'' is dependent on solids loading. There are two ways then that the pore size can be controlled. The first is to change the solidification speed that then alters the microstructural wavelength, or the solids loading can be changed. In doing so, the ratio of pore size to wall size is changed.
It is often more prudent to alter the solidification velocity seeing as a minimum solid loading is usually desired. Since microstructural size ''(λ)'' is inversely related to the velocity of the freezing front, faster speeds lead to finer structures, while slower speeds produce a coarse microstructure. Controlling the solidification velocity is, therefore, crucial to being able to control the microstructure.

Controlling pore shape

Additives can prove highly useful and versatile in changing the morphology of pores. These work by affecting the growth kinetics and microstructure of the ice in addition to the topology of the ice-water interface. Some additives work by altering the phase diagram of the solvent. For example, water and NaCl have a eutectic phase diagram. When NaCl is added into a freeze-casting suspension, the solid ice phase and liquid regions are separated by a zone where both solids and liquids can coexist. This briny region is removed during sublimation, but its existence has a strong effect on the microstructure of the porous ceramic. Other additives work by either altering the interfacial surface energies between the solid/liquid and particle/liquid, changing the viscosity of the suspension, or the degree of undercooling in the system. Studies have been done with glycerol, sucrose, ethanol, acetic acid and more.

Static vs. dynamic freezing profiles

If a freeze casting setup with a constant temperature on either side of the freezing system is used, (static freeze-casting) the front solidification velocity in the SSZ will decrease over time due to the increasing thermal buffer caused by the growing ice front. When this occurs, more time is given for the anisotropic ice crystals to grow perpendicularly to the freezing direction (c-axis) resulting in a structure with ice lamellae that increase in thickness along the length of the sample.

To ensure highly anisotropic, yet predictable solidification behavior within the SSZ, dynamic freezing patterns are preferred. Using dynamic freezing, the velocity of the solidification front, and, therefore, the ice crystal size, can be controlled with a changing temperature gradient. The increasing thermal gradient counters the effect of the growing thermal buffer imposed by the growing ice front. It has been shown that a linearly decreasing temperature on one side of a freeze-cast will result in near-constant solidification velocity, yielding ice crystals with an almost constant thickness along the SSZ of an entire sample. However, as pointed out by Waschkies et al. even with constant solidification velocity, the thickness of the ice crystals does increase slightly over the course of freezing. In contrast to that, Flauder et al. demonstrated that an exponential change of the temperature at the cooling plate leads to a constant ice crystal thickness within the complete SSZ, which was attributed to a measurably constant ice-front velocity in a distinct study. This approach enables a prediction of the ice-front velocity from the thermal parameters of the suspension. Consequently, if the exact relationship between the pore diameter and ice-front velocity is known, an exact control over the pore diameter can be achieved.

Anisotropy of the interface kinetics

Even if the temperature gradient within the slurry is perfectly vertical, it is common to see tilting or curvature of the lamellae as they grow through the suspension. To explain this, it is possible to define two distinct growth directions for each ice crystal. There is the direction determined by the temperature gradient, and the one defined by the preferred growth direction crystallographically speaking. These angles are often at odds with one another, and their balance will describe the tilt of the crystal.

The non-overlapping growth directions also help to explain why dendritic textures are often seen in freeze-casts. This texturing is usually found only on the side of each lamella; the direction of the imposed temperature gradient. The ceramic structure left behind shows the negative image of these dendrites. In 2013, Deville et al. made the observation that the periodicity of these dendrites (tip-to-tip distance) actually seems to be related to the primary crystal thickness.

Particle packing effects

Up until now, the focus has been mostly on the structure of the ice itself; the particles are almost an afterthought to the templating process but in fact, the particles can and do play a significant role during freeze-casting. It turns out that particle arrangement also changes as a function of the freezing conditions. For example, researchers have shown that freezing velocity has a marked effect on wall roughness. Faster freezing rates produce rougher walls since particles are given insufficient time to rearrange. This could be of use when developing permeable gas transfer membranes where tortuosity and roughness could impede gas flow. It also turns out that z- and r-crystals do not interact with ceramic particles in the same way. The z-crystals pack particles in the x-y plane while r-crystals pack particles primarily in the z-direction. R-crystals actually pack particles more efficiently than z-crystals and because of this, the area fraction of the particle-rich phase (1 - area fraction of ice crystals) changes as the crystal population shifts from a mixture of z- and r-crystals to only z-crystals. Starting from where ice crystals first begin to exclude particles, marking the beginning of the transition zone, we have a majority of r-crystals and a high value for the particle-rich phase fraction. We can assume that because the solidification speed is still rapid that the particles will not be packed efficiently. As the solidification rate slows down, however, the area fraction of the particle-rich phase drops indicating an increase in packing efficiency. At the same time, the competitive growth process is taking place, replacing r-crystals with z-crystals. At a certain point nearing the end of the transition zone, the particle-rich phase fraction rises sharply since z-crystals are less efficient at packing particles than r-crystals. The apex of this curve marks the point where only z-crystals are present (SSZ). During steady-state growth, after the maximum particle-rich phase fraction is reached, the efficiency of packing increases as steady-state is achieved. In 2011, researchers at Yale University set out to probe the actual spatial packing of particles within the walls. Using small-angle X-ray scattering (SAXS) they characterized the particle size, shape and interparticle spacing of nominally 32 nm silica suspensions that had been freeze-cast at different speeds. Computer simulations indicated that for this system, the particles within the walls should not be touching but rather separated from one another by thin films of ice. Testing, however, revealed that the particles were, in fact, touching and more than that, they attained a packed morphology that cannot be explained by typical equilibrium densification processes.
Particle-scale structure in frozen colloidal suspensions from small angle x-ray scattering

Morphological instabilities

In an ideal world, the spatial concentration of particles within the SSZ would remain constant throughout solidification. As it happens, though, the concentration of particles does change during compression, and this process is highly sensitive to solidification speed. At low freezing rates, Brownian motion takes place, allowing particles to move easily away from the solid-liquid interface and maintain a homogeneous suspension. In this situation, the suspension is always warmer than the solidified portion. At fast solidification speeds, approaching VC, the concentration, and concentration gradient at the solid-liquid interface increases because particles cannot redistribute soon enough. When it has built up enough, the freezing point of the suspension is below the temperature gradient in the solution and morphological instabilities can occur. For situations where the particle concentration bleeds into the diffusion layer, both the actual and freezing temperature dip below the equilibrium freezing temperature creating an unstable system. Often, these situations lead to the formation of what are known as ice lenses.

These morphological instabilities can trap particles, preventing full redistribution and resulting in inhomogeneous distribution of solids along the freezing direction as well as discontinuities in the ceramic walls, creating voids larger than intrinsic pores within the walls of the porous ceramic.

Mechanical properties

Most research into the mechanical properties of freeze casted structures focus on the compressive strength of the material and its yielding behavior at increasing stresses. According to Ashby, the mechanical properties of a freeze-casted, open pore structure can be approximately modeled with an anisotropic, cellular solid. These include naturally occurring materials such as cork and wood that have properties that have anisotropic structures, and thus mechanical properties that are directionally dependent. Donius et al. have investigated the anisotropic nature of freeze-casted aerogels, comparing their mechanical strength to isotropically freeze casted aerogels. They found that the Young's modulus of the anisotropic structure was significantly higher than that of the isotropic aerogels, particularly when tested parallel to the freezing direction. The Young's modulus is several orders of magnitude higher in the parallel direction as compared to the direction perpendicular to freezing, demonstrating the anisotropic mechanical properties.

The mechanical behavior of the freeze casted structure can be classified into distinct regions. At low strains, the lamallae follow a linear elastic behavior. Here, the lamellae bend under a compressive stress, and thus deflect. According to Ashby, this deflection can be calculated from single beam theory, in which each of the cellular sections are idealized to be cubic shaped where each of the cell walls are assumed to be beam-like members with a square base. Based on this idealization, the amount of bending $\delta$ in the cell walls under a compressive force $F$ is given by $\delta = \frac$ where $l$ is the length of each cell, $I$ is the second moment of area, $E_s$ is the Young's modulus of the cell wall material and $C_1$ is a geometry dependent constant. Furthermore, we find that the Young's modulus of the entire structure $E$ is proportional to the square of the relative density: $\frac = C_2 (\frac)^2$. This shows that the density of the material is an important factor when designing structures that can withstand loads, and that the Young's modulus of the structure is heavily determined by the porosity of the structure. Past the linear region, the lamellae start to buckle elastically and deform non-linearly. In a stress-strain curve, this is shown as a flat plateau. The critical load at which buckling begins is given by: $F_ = \frac$ where $n^2$ is a constant dependent on the boundary constraints of the structure. This is one of the main failure mechanisms for freeze casted materials. From this, the maximum compressive stress that an anisotropic porous solid can maintain is given by $\sigma^ = \sigma_\frac$ where $\sigma_$ is the fracture stress for the bulk material. These models demonstrate that the bulk material selection can drastically impact the mechanical response of freeze casted structures under stress. Other microstructural features such as the lamellar thickness, pore morphology and degree of macroporosity can also heavily influence the compressive strength and Young's modulus of these highly anisotropic structures.

Novel freeze-casting techniques

Freeze-casting can be applied to produce aligned porous structure from diverse building blocks including ceramics, polymers, biomacromolecules, graphene and carbon nanotubes. As long as there are particles that may be rejected by a progressing freezing front, a templated structure is possible. By controlling cooling gradients and the distribution of particles during freeze casting, using various physical means, the orientation of lamellae in obtained freeze cast structures can be controlled to provide improved performance in diverse applied materials. Munch et al. showed that it is possible to control the long-range arrangement and orientation of crystals normal to the growth direction by templating the nucleation surface. This technique works by providing lower energy nucleation sites to control the initial crystal growth and arrangement. The orientation of ice crystals can also be affected by applying electromagnetic fields as was demonstrated in 2010 by Tang et al. in 2012 by Porter et al., and in 2021 by Yin et al. Using specialized setups, researchers have been able to create radially aligned freeze-casts tailored for biomedical applications and filtration or gas separation applications. Inspired by nature, scientists have also been able to use coordinating chemicals and cryopreserved to create remarkably distinctive microstructural architectures.

Freeze cast materials

Particles that are assembled into aligned porous materials in freeze casting processes are often referred to as building blocks. As freeze casting has become a widespread technique the range of materials used has expanded. In recent years, graphene and carbon nanotubes have been used to fabricate controlled porous structures using freeze casting methods, with materials often exhibiting outstanding properties. Unlike aerogel materials produced without ice-templating, freeze cast structures of carbon nanomaterials have the advantage of possessing aligned pores, allowing, for example unparalleled combinations of low density and high conductivity.

Applications of freeze cast materials

Freeze casting is unique in its ability to produce aligned pore structures. Such structures are often found in nature, and consequently freeze casting has emerged as a valuable tool to fabricate biomimetic structures. The transport of fluids through aligned pores has led to the use of freeze casting as a method towards biomedical applications including bone scaffold materials. The alignment of pores in freeze cast structures also imparts extraordinarily high thermal resistance in the direction perpendicular to the aligned pores. The freeze casting o
aligned porous fibres
by spinning processes presents a promising method towards the fabrication of high performance insulating clothing articles. In addition, materials with aligned pores produced via freeze casting from sintered nickel powder have gained significant attention in phase-change systems, such as loop heat pipes (LHPs), due to their excellent thermal properties. In these systems, wicks play a critical role in maintaining liquid-vapor equilibrium, enabling efficient circulation of the working fluid. Traditional wicks are often manufactured separately and integrated later, creating an interface that limits liquid transfer efficiency. To address this limitation, a porous wick with a gradient structure was developed in a single operation using freeze casting. This innovative approach eliminates interfacial resistance, ensuring seamless liquid transport while maintaining the high thermal conductivity and efficient capillary action required for optimal LHP performance.

Another emerging and promising application of freeze casting is the production of porous foams for green hydrogen generation through advanced thermochemical processes like Chemical Loop Combustion (CLC) and the Steam Iron Process (SIP). These processes leverage the unique properties of porous metal structures, such as optimized reaction kinetics, enhanced thermal efficiency, and sustainability. In Chemical Loop Combustion (CLC), foams made from materials like iron oxides act as oxygen carriers, enabling fuel combustion without direct air contact, separating CO₂ for capture while producing high-purity hydrogen. Similarly, in Steam Iron Process (SIP), dendritic pore structures ensure efficient water vapor distribution and maximize hydrogen yield. The precise control over porosity and thermal properties afforded by freeze casting, along with the use of eco-friendly solvents like camphene, positions these foams as a vital innovation for scalable and sustainable hydrogen production, contributing to the fight against climate change.

See also

*Freeze gelation

Further reading

*
*J. Laurie, ''Freeze Casting: a Modified Sol-Gel Process'', University of Bath, UK, Ph.D. Thesis, 1995
*M. Statham, ''Economic Manufacture of Freeze-Cast Ceramic Substrate Shapes for the Spray-Forming Process'', Univ. Bath, UK, Ph.D. Thesis, 1998
*S. Deville, "Freezing Colloids: Observations, Principles, Control, and Use." Springer, 2017
*

External links

*A website with large dataset, allowing creation of graph

References

{{Reflist, colwidth=30em

Casting (manufacturing)
Ceramic engineering
Colloids
Water ice