Sintering, densification and microstructure control
**Guiding question:** When does heating remove porosity, and when does it merely change grain size?
§3.1 Why this unit matters to ceramic engineering
Guiding question: When does heating remove porosity, and when does it merely change grain size?
Two kiln programs produce the same apparent mean relative density but different grain size tails and strength scatter. This is a question about pore elimination and grain growth and the distinction between densification versus grain coarsening. A defensible answer identifies feedstock composition, chemistry, processing path and microstructure before comparing measured properties. It also names the potentially controlling failure mode: a recipe that optimizes density, dielectric response or simple bend strength can still worsen thermal cycling, porosity topology, lifetime or chemical compatibility. The technical sections below connect simplified mathematics to the kinds of measurements that would justify a practical claim.
The shrinking compact evolves through interconnected and then isolated pores while grains grow and the driving force changes. Surface diffusion can create visible necks without moving particle centers closer.
Separate shrinkage, mass loss, pore topology and grain-size distributions. The temperature at the sample core may differ from a furnace set point, especially during rapid cycles.
§3.2 Measurement and conservation checkpoints
Write down the starting composition, defect-sensitive atmosphere and powder or glass history; define whether reported fractions are by area, mass, moles or volume; and record the process history that could alter furnace atmosphere, time and geometry. For any local observation, describe why its sampling plan represents the part. Treat trapped gases and nonuniform heating as an explicit rival explanation and plan dilatometry plus volumetric imaging as an independent check.
§3.3 Early-stage necking and transport mechanisms
Scientific foundation
During initial sintering, atoms redistribute to form necks between particles. Surface diffusion can enlarge a neck without reducing pore volume, while grain-boundary and lattice diffusion paths may transport matter so that particle centers approach and densification occurs. At a fixed temperature, more time is not synonymous with constant densification rate; geometry and the dominant mechanism evolve. Idealized two-sphere formulas apply only over their derived small-neck regime.
Governing relationship or decision rule: In a restricted two-sphere viscous-flow example, x/r proportional to [gamma t/(eta r)]^(1/2), with model-specific constants. Read every term in the equation as part of a declared specimen condition. In early-stage necking and transport mechanisms, a valid calculation requires the state of the ceramic, unit system and physical or geometrical meaning of each quantity. A nominally identical oxide composition can conceal different dopant distributions, porosity, grain-boundary glass and processing damage. Equal bulk density does not guarantee equivalent fracture, transport or thermal stability. Interpret measurements only after identifying which feature of the state they can support.
Worked ceramic investigation
Use the early-stage viscous-flow teaching expression for neck size, with fictional surface energy, viscosity and radius, to compare two times. State explicitly that it represents a different physical mechanism from diffusion-driven densification in a crystalline alumina compact. Sketch where mass comes from and where it is deposited for a surface-diffusion mechanism versus a densifying mechanism.
A companion calculator called Illustrative viscous neck ratio supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Isothermal hold duration = 214.2857 s, the declared formula returns 0.17928429 neck/radius; at Isothermal hold duration = 428.5714 s, it returns 0.25354628 neck/radius. On this interval the modeled value increases, and the two-point slope is 0.00034655593 neck/radius per s. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
The browser expression is Math.sqrt(3t/(21e10*1e-6)). Model validity is limited to: Two equal spheres r=1 µm, viscosity 1e10 Pa·s, surface energy 1 J/m²; early-stage toy viscous-flow only, not crystalline solid-state sintering. Reproduce the first result with a calculator using the expression and units stated here. Before accepting the number, check the units analytically and identify what is being held constant. If a result looks realistic, that coincidence does not establish physical validation. Try an extreme permitted slider input and ask whether the computed change is physically interpretable or an artifact of the simplifying assumptions.
Evidence, uncertainty, and a competing mechanism
Real particle contacts vary in curvature and coordination number. The illustrative viscosity may be inappropriate for a non-glassy specimen, and extrapolating the formula beyond small necks is physically invalid.
Exercises and self-assessment
Exercise 3.1A — Physical mechanism. How can a micrograph show larger particle necks without any measurable overall shrinkage? State one plausible explanation, one competing explanation and the material-state observation that separates them. Relate the claim to the governing relationship written above rather than merely repeating definitions.
Exercise 3.1B — Numerical cross-check. Recalculate 0.17928429 neck/radius from the companion formula at 214.2857 s and compare with 0.25354628 neck/radius at 428.5714 s. State the sign of the sensitivity and explain whether the changed input would reasonably be controlled independently in a real ceramic specimen.
Exercise 3.1C — Test design. Propose a preparation or test matrix that addresses trapped gases and nonuniform heating. Give the specimen orientation, spatial sampling scheme, temperature-time record and principal instrument calibration. Explain which observation would falsify your first hypothesis and name a condition under which your model should not be used.
§3.4 Densification and isotropic shrinkage
Scientific foundation
Relative density is bulk density divided by the theoretical density for the relevant phase mixture. If specimen mass is conserved and shrinkage is isotropic, linear dimensions scale with the cube root of the initial-to-final volume ratio. In a porous ceramic, these assumptions can fail due to chemical loss, closed pores, anisotropic forming texture, constrained substrates or density variation along a part. Reporting only percentage linear shrinkage hides directional distortion and possible weight loss.
Governing relationship or decision rule: With constant mass and isotropic dimensions, L_f/L_0=(rho_0/rho_f)^(1/3), where rho values are absolute or consistently normalized densities. Read every term in the equation as part of a declared specimen condition. In densification and isotropic shrinkage, a valid calculation requires the state of the ceramic, unit system and physical or geometrical meaning of each quantity. A nominally identical oxide composition can conceal different dopant distributions, porosity, grain-boundary glass and processing damage. Equal bulk density does not guarantee equivalent fracture, transport or thermal stability. Interpret measurements only after identifying which feature of the state they can support.
Worked ceramic investigation
Starting from a hypothetical green relative density of 0.55, calculate isotropic linear shrinkage at final relative densities 0.90 and 0.98 and compare with a length-gauge measurement plan. Distinguish a simple conservation prediction from a fitted dilatometer trace. Include a mass-loss check to test whether binder residue or volatile reactants invalidate the ideal model.
A companion calculator called Ideal isotropic linear shrinkage supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Final relative density = 0.8185714 fraction, the declared formula returns 12.413864 % length decrease; at Final relative density = 0.8528571 fraction, it returns 13.603636 % length decrease. On this interval the modeled value increases, and the two-point slope is 34.701692 % length decrease per fraction. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
The browser expression is 100*(1-Math.pow(0.55/rhof,1/3)). Model validity is limited to: Start relative density=0.55, conserved mass, unchanged reference solid density, spatially uniform isotropic shrinkage. Reproduce the first result with a calculator using the expression and units stated here. Before accepting the number, check the units analytically and identify what is being held constant. If a result looks realistic, that coincidence does not establish physical validation. Try an extreme permitted slider input and ask whether the computed change is physically interpretable or an artifact of the simplifying assumptions.
Evidence, uncertainty, and a competing mechanism
The density ratio must be based on consistent skeleton chemistry before and after firing; an actual phase reaction can change theoretical density. Porosity may not remain spatially uniform.
Exercises and self-assessment
Exercise 3.2A — Physical mechanism. Why can a sintering shrinkage curve look smooth while internal closed pores grow? State one plausible explanation, one competing explanation and the material-state observation that separates them. Relate the claim to the governing relationship written above rather than merely repeating definitions.
Exercise 3.2B — Numerical cross-check. Recalculate 12.413864 % length decrease from the companion formula at 0.8185714 fraction and compare with 13.603636 % length decrease at 0.8528571 fraction. State the sign of the sensitivity and explain whether the changed input would reasonably be controlled independently in a real ceramic specimen.
Exercise 3.2C — Test design. Propose a preparation or test matrix that addresses trapped gases and nonuniform heating. Give the specimen orientation, spatial sampling scheme, temperature-time record and principal instrument calibration. Explain which observation would falsify your first hypothesis and name a condition under which your model should not be used.
§3.5 Grain growth and pore pinning
Scientific foundation
Grain boundaries lower total interfacial energy when their network coarsens, but moving boundaries interact with pores and second-phase particles. Pore drag may delay boundary motion while isolated closed pores can become trapped inside fast-growing grains. Excessive grain growth sometimes reduces mechanical strength through flaw enlargement or anisotropic expansion. Grain growth kinetics can be fitted with D^n-D0^n=Kt, but the exponent is not a universal constant across temperatures and mechanisms.
Governing relationship or decision rule: A common empirical growth form is D^n-D_0^n=Kt; n and K must be validated for that processing regime. Read every term in the equation as part of a declared specimen condition. In grain growth and pore pinning, a valid calculation requires the state of the ceramic, unit system and physical or geometrical meaning of each quantity. A nominally identical oxide composition can conceal different dopant distributions, porosity, grain-boundary glass and processing damage. Equal bulk density does not guarantee equivalent fracture, transport or thermal stability. Interpret measurements only after identifying which feature of the state they can support.
Worked ceramic investigation
Compare two hypothetical firing schedules achieving nearly identical final bulk density yet differing in measured grain-size distribution. Use microscopy to count both interior and boundary pores, and then select a mechanical test sensitive to the dangerous tail of the flaw population. Ask whether the schedule with the larger average grain is necessarily inferior for electrical transport, where grain-boundary resistance may matter.
A companion calculator called Illustrative viscous neck ratio supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Isothermal hold duration = 642.8571 s, the declared formula returns 0.3105295 neck/radius; at Isothermal hold duration = 857.1429 s, it returns 0.35856858 neck/radius. On this interval the modeled value increases, and the two-point slope is 0.00022418238 neck/radius per s. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
Evidence, uncertainty, and a competing mechanism
Measurements from one etched section cannot determine three-dimensional pore topology. Abnormal grain growth may be controlled by rare seed grains and transient liquid films not captured by a mean-grain model.
Exercises and self-assessment
Exercise 3.3A — Physical mechanism. Under what conditions can additional soaking reduce strength after porosity seems already minimal? State one plausible explanation, one competing explanation and the material-state observation that separates them. Relate the claim to the governing relationship written above rather than merely repeating definitions.
Exercise 3.3B — Numerical cross-check. Recalculate 0.3105295 neck/radius from the companion formula at 642.8571 s and compare with 0.35856858 neck/radius at 857.1429 s. State the sign of the sensitivity and explain whether the changed input would reasonably be controlled independently in a real ceramic specimen.
Exercise 3.3C — Test design. Propose a preparation or test matrix that addresses trapped gases and nonuniform heating. Give the specimen orientation, spatial sampling scheme, temperature-time record and principal instrument calibration. Explain which observation would falsify your first hypothesis and name a condition under which your model should not be used.
§3.6 Liquid-phase and reactive sintering
Scientific foundation
A transient liquid may speed material redistribution through dissolution and reprecipitation, particle rearrangement or viscous flow. Wettability and liquid volume determine whether liquid bridges particles or creates continuous weak intergranular films. Reactive sintering simultaneously changes chemistry, phase fractions and volume, so a dimensional shrinkage measurement is not evidence of pore elimination by itself. Reaction-induced expansion can partially compensate densification shrinkage.
Governing relationship or decision rule: Liquid-phase sintering requires liquid formation, wetting and a transport pathway; conservation of each element constrains the final phase assemblage. Read every term in the equation as part of a declared specimen condition. In liquid-phase and reactive sintering, a valid calculation requires the state of the ceramic, unit system and physical or geometrical meaning of each quantity. A nominally identical oxide composition can conceal different dopant distributions, porosity, grain-boundary glass and processing damage. Equal bulk density does not guarantee equivalent fracture, transport or thermal stability. Interpret measurements only after identifying which feature of the state they can support.
Worked ceramic investigation
For a hypothetical silicon-nitride system with a small oxide sintering aid, propose a heating cycle that maps wetting and residual intergranular composition. Compare fired density, high-temperature creep and oxidation instead of selecting the maximum as-fired bend strength alone. Use a mass and phase ledger to track whether a volatile loss changes reference density.
A companion calculator called Ideal isotropic linear shrinkage supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Final relative density = 0.8871429 fraction, the declared formula returns 14.731287 % length decrease; at Final relative density = 0.9214286 fraction, it returns 15.802279 % length decrease. On this interval the modeled value increases, and the two-point slope is 31.237267 % length decrease per fraction. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
Evidence, uncertainty, and a competing mechanism
Real sintering-aid phase chemistry is multicomponent. Changing sintering additive concentration may also change grain aspect ratio and crystallization; the contribution of each must be separated experimentally.
Exercises and self-assessment
Exercise 3.4A — Physical mechanism. Why might a ceramic densify faster with a liquid phase but become less reliable at elevated temperature? State one plausible explanation, one competing explanation and the material-state observation that separates them. Relate the claim to the governing relationship written above rather than merely repeating definitions.
Exercise 3.4B — Numerical cross-check. Recalculate 14.731287 % length decrease from the companion formula at 0.8871429 fraction and compare with 15.802279 % length decrease at 0.9214286 fraction. State the sign of the sensitivity and explain whether the changed input would reasonably be controlled independently in a real ceramic specimen.
Exercise 3.4C — Test design. Propose a preparation or test matrix that addresses trapped gases and nonuniform heating. Give the specimen orientation, spatial sampling scheme, temperature-time record and principal instrument calibration. Explain which observation would falsify your first hypothesis and name a condition under which your model should not be used.
§3.7 Pressure-assisted sintering and defect inspection
Scientific foundation
Hot pressing, spark-plasma-assisted sintering and hot isostatic pressing use stress and thermal history to change diffusion paths, contact forces and pore closure. Pressure can reduce open porosity yet fail to heal flaws isolated by chemistry, contamination or trapped gas. Temperature measurement in fast electric-current processing is especially sensitive to tooling geometry and sensor location; labels such as spark plasma do not establish a unique microscopic transport mechanism.
Governing relationship or decision rule: The applied mechanical work and heat flow affect densification kinetics, but a universal pressure-to-density equation is not justified. Read every term in the equation as part of a declared specimen condition. In pressure-assisted sintering and defect inspection, a valid calculation requires the state of the ceramic, unit system and physical or geometrical meaning of each quantity. A nominally identical oxide composition can conceal different dopant distributions, porosity, grain-boundary glass and processing damage. Equal bulk density does not guarantee equivalent fracture, transport or thermal stability. Interpret measurements only after identifying which feature of the state they can support.
Worked ceramic investigation
Define a matched comparison of pressureless and assisted routes using the same starting powder and final geometry. Record temperature gradients and tool-contact contamination. Inspect density spatially with Archimedes measurements plus an imaging technique capable of detecting large defects. If one process is faster, evaluate the total energy and yield before calling it inherently greener.
A companion calculator called Illustrative viscous neck ratio supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Isothermal hold duration = 1071.429 s, the declared formula returns 0.40089186 neck/radius; at Isothermal hold duration = 1285.714 s, it returns 0.43915503 neck/radius. On this interval the modeled value increases, and the two-point slope is 0.00017856146 neck/radius per s. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
Evidence, uncertainty, and a competing mechanism
The measured relative density can be near unity while rare large flaws still dominate flexural strength. Applied pressure may introduce texture that matters in anisotropic ceramics.
Exercises and self-assessment
Exercise 3.5A — Physical mechanism. What evidence would be needed to claim that applied pressure closed the critical cracks rather than only smaller pores? State one plausible explanation, one competing explanation and the material-state observation that separates them. Relate the claim to the governing relationship written above rather than merely repeating definitions.
Exercise 3.5B — Numerical cross-check. Recalculate 0.40089186 neck/radius from the companion formula at 1071.429 s and compare with 0.43915503 neck/radius at 1285.714 s. State the sign of the sensitivity and explain whether the changed input would reasonably be controlled independently in a real ceramic specimen.
Exercise 3.5C — Test design. Propose a preparation or test matrix that addresses trapped gases and nonuniform heating. Give the specimen orientation, spatial sampling scheme, temperature-time record and principal instrument calibration. Explain which observation would falsify your first hypothesis and name a condition under which your model should not be used.
§3.8 Integrated investigation: Choose a reliable sintering schedule
The evidence problem
Two kiln programs produce the same apparent mean relative density but different grain size tails and strength scatter. The problem cannot be solved by selecting the equation that gives the most convenient number. First define which property must be achieved, the physical constraint that could prevent it, and the microstructural condition of the specimen being evaluated. List known facts separately from inferred mechanisms, uncertain constants and quantities not yet measured. The central decision is not simply which powder process or ceramic composition seems attractive: it is what observations would allow a researcher or an engineer to prefer one explanation over competing causes.
Integrated analysis and model workflow
Contrast pore closure, abnormal grain growth, trapped gas and near-surface versus interior temperatures. Build a minimal processing–structure–property map linking the material history with observed state and relevant response. Carry the distinction between pore elimination and grain growth and densification versus grain coarsening through the map. For phase equilibria, show composition conservation and do not infer processing kinetics from a phase diagram. For a densification relation, separate pore disappearance from neck growth and grain coarsening. For electrical function, distinguish ion migration, electronic leakage, electrode polarization and nonlinear ferroelectric switching. Use the actual specimen time–temperature history instead of only a furnace setting. For damage, compare local flaw driving force with the appropriate fracture resistance.
Validation plan, report, and counterfactual
Collect dilatometry, 3D defect imaging and strength distributions from held-out production geometries. Include an alternative interpretation capable of explaining the same first observation and choose a different instrument or held-out processing condition that separates it. State which measurement would cause you to reject the favored hypothesis. Retain raw counts, image scales, processing logs and sample identifiers. If comparing ceramic families, write down the constraints that matter for weight, stiffness, temperature, oxidation and availability rather than reporting only one maximum strength metric.
Known limits: Final relative density alone does not rank reliability. A reader should be able to reconstruct the entire inference from original specimen provenance and measured inputs. The synthetic datasets intentionally contain no physical replicate scatter or hidden model mismatch; those omissions would make a real-world confidence interval far wider than a formula plot suggests.
Five cumulative exam questions
- Draw a cause-and-effect diagram linking composition, process history, pore elimination and grain growth and the decision-controlling property; mark one arrow that cannot yet be established from the observations.
- Write a conservation or constitutive relation for densification versus grain coarsening with all units and physical validity assumptions; interpret one limiting case.
- Propose two mutually exclusive experimental signatures for the favored and competing mechanisms and name the required independent instruments.
- Show how trapped gases and nonuniform heating could reverse a seemingly beneficial result, and design the minimum additional test needed to detect that reversal.
- Write a one-paragraph evidence statement specifying what has been demonstrated, what is merely predicted, and which industrial decision is unjustified until qualification.
Step-by-Step Statistical Simulation and Monte Carlo Solutions
Three original questions connect calculations, definitions and interpretation. Open each solution to follow the reasoning.
Use the early-stage viscous-flow teaching expression for neck size, with fictional surface energy, viscosity and radius, to compare two times. State explicitly that it represents a different physical mechanism from diffusion-driven densification in a crystalline alumina compact. Sketch where mass comes from and where it is deposited for a surface-diffusion mechanism versus a densifying mechanism.
A companion calculator called Illustrative viscous neck ratio supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Isothermal hold duration = 214.2857 s, the declared formula returns 0.17928429 neck/radius; at Isothermal hold duration = 428.5714 s, it returns 0.25354628 neck/radius. On this interval the modeled value increases, and the two-point slope is 0.00034655593 neck/radius per s. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
The browser expression is Math.sqrt(3t/(21e10*1e-6)). Model validity is limited to: Two equal spheres r=1 µm, viscosity 1e10 Pa·s, surface energy 1 J/m²; early-stage toy viscous-flow only, not crystalline solid-state sintering. Reproduce the first result with a calculator using the expression and units stated here. Before accepting the number, check the units analytically and identify what is being held constant. If a result looks realistic, that coincidence does not establish physical validation. Try an extreme permitted slider input and ask whether the computed change is physically interpretable or an artifact of the simplifying assumptions.
Starting from a hypothetical green relative density of 0.55, calculate isotropic linear shrinkage at final relative densities 0.90 and 0.98 and compare with a length-gauge measurement plan. Distinguish a simple conservation prediction from a fitted dilatometer trace. Include a mass-loss check to test whether binder residue or volatile reactants invalidate the ideal model.
A companion calculator called Ideal isotropic linear shrinkage supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Final relative density = 0.8185714 fraction, the declared formula returns 12.413864 % length decrease; at Final relative density = 0.8528571 fraction, it returns 13.603636 % length decrease. On this interval the modeled value increases, and the two-point slope is 34.701692 % length decrease per fraction. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
The browser expression is 100*(1-Math.pow(0.55/rhof,1/3)). Model validity is limited to: Start relative density=0.55, conserved mass, unchanged reference solid density, spatially uniform isotropic shrinkage. Reproduce the first result with a calculator using the expression and units stated here. Before accepting the number, check the units analytically and identify what is being held constant. If a result looks realistic, that coincidence does not establish physical validation. Try an extreme permitted slider input and ask whether the computed change is physically interpretable or an artifact of the simplifying assumptions.
Compare two hypothetical firing schedules achieving nearly identical final bulk density yet differing in measured grain-size distribution. Use microscopy to count both interior and boundary pores, and then select a mechanical test sensitive to the dangerous tail of the flaw population. Ask whether the schedule with the larger average grain is necessarily inferior for electrical transport, where grain-boundary resistance may matter.
A companion calculator called Illustrative viscous neck ratio supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Isothermal hold duration = 642.8571 s, the declared formula returns 0.3105295 neck/radius; at Isothermal hold duration = 857.1429 s, it returns 0.35856858 neck/radius. On this interval the modeled value increases, and the two-point slope is 0.00022418238 neck/radius per s. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
For a hypothetical silicon-nitride system with a small oxide sintering aid, propose a heating cycle that maps wetting and residual intergranular composition. Compare fired density, high-temperature creep and oxidation instead of selecting the maximum as-fired bend strength alone. Use a mass and phase ledger to track whether a volatile loss changes reference density.
A companion calculator called Ideal isotropic linear shrinkage supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Final relative density = 0.8871429 fraction, the declared formula returns 14.731287 % length decrease; at Final relative density = 0.9214286 fraction, it returns 15.802279 % length decrease. On this interval the modeled value increases, and the two-point slope is 31.237267 % length decrease per fraction. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.
Define a matched comparison of pressureless and assisted routes using the same starting powder and final geometry. Record temperature gradients and tool-contact contamination. Inspect density spatially with Archimedes measurements plus an imaging technique capable of detecting large defects. If one process is faster, evaluate the total energy and yield before calling it inherently greener.
A companion calculator called Illustrative viscous neck ratio supplies a separate, deliberately limited numerical experiment. Its controls isolate one effect so a reader can check dimensions and sensitivity instead of guessing the direction of change. At Isothermal hold duration = 1071.429 s, the declared formula returns 0.40089186 neck/radius; at Isothermal hold duration = 1285.714 s, it returns 0.43915503 neck/radius. On this interval the modeled value increases, and the two-point slope is 0.00017856146 neck/radius per s. These are worked calculations, not fabricated experimental observations or certified properties of any real ceramic formulation.