Transient Conduction: Semi-Infinite Solid Solutions

Understanding how heat penetrates thick materials during brief heating or cooling events is essential for designing processes from metal quenching to geothermal systems. The semi-infinite solid model provides a powerful analytical shortcut for these scenarios, allowing engineers to predict temperature changes without solving complex boundary value problems. Mastering this concept enables you to analyze transient thermal responses in everything from concrete slabs to the Earth's crust.

What is a Semi-Infinite Solid?

In heat transfer, a semi-infinite solid is an idealization where one boundary of a body is defined, but the other extends to infinity. This approximation is valid when a thermal disturbance at the surface has not had enough time to penetrate to the far boundary. Physically, this means the material is so thick that for the duration of the event, the interior remains at its initial temperature. The key criterion is that the thermal penetration depth is much smaller than the actual thickness of the object. You encounter this condition in short-duration processes, such as the initial seconds of quenching a thick steel block or the diurnal heating of a deep concrete wall.

The model relies on several simplifying assumptions: the material is homogeneous and isotropic, its properties are constant, and the initial temperature $T_i$ is uniform. Most importantly, the surface experiences a sudden change to a new temperature $T_s$, which is then maintained. This step change is a common boundary condition, though other surface fluxes can be analyzed. By adopting this idealized geometry, the complex partial differential equation of transient conduction reduces to a form solvable with a similarity variable, combining depth and time into a single parameter.

The Complementary Error Function Solution

The governing equation for one-dimensional, transient conduction without heat generation is the heat diffusion equation: $$\frac{\partial T}{\partial t}=\alpha \frac{{\partial }^{2}T}{\partial x^2}.$$ Here, $T$ is temperature, $t$ is time, $x$ is the distance from the surface, and $\alpha$ is the thermal diffusivity, a material property that measures how quickly heat diffuses through a substance. For a semi-infinite solid with a constant surface temperature change, the solution is expressed using the complementary error function, denoted as erfc.

The temperature at any depth $x$ and time $t$ is given by: $$T(x,t)=T_i+(T_s-T_i)\text{erfc}\left(\frac{x}{2\sqrt{\alpha t}}\right).$$ The term $\frac{x}{2\sqrt{\alpha t}}$ is the dimensionless similarity variable, often represented by $\eta$. The complementary error function erfc($\eta$) is defined as 1 - erf($\eta$), where erf is the standard error function. This function decreases monotonically from 1 at $\eta =0$ (the surface) to 0 as $\eta$ approaches infinity (deep within the solid). This mathematical form directly shows that temperature change propagates into the material in a predictable, smooth manner dependent on both position and time.

Interpreting Temperature Profiles and Key Parameters

The solution reveals that temperature distribution is not linear but follows the erfc curve, which has a steep gradient near the surface that gradually flattens with depth. The argument of the erfc function, $\eta =x/(2\sqrt{\alpha t})$, is crucial. It tells you that depth $x$ and time $t$ are not independent; doubling the depth has the same effect on temperature as quadrupling the time, due to the square root dependence. This is why heat penetration slows down considerably as time progresses.

Thermal diffusivity $\alpha$ is the engine of this process. Materials with high $\alpha$, like metals, allow heat to penetrate quickly, leading to a shallower temperature gradient for a given time. Materials with low $\alpha$, such as plastics or soil, cause heat to linger near the surface, creating a steeper gradient. You can use the solution to answer practical questions: For example, to find the time it takes for a specific temperature to reach a certain depth, you would set $T(x,t)$ to the desired value and solve for $t$, which often requires consulting erfc tables or using computational tools.

A related useful concept is the penetration depth, often defined as the depth where the temperature change is 1% of the surface change. This can be approximated as $x\approx 4\sqrt{\alpha t}$. This rule of thumb helps you quickly estimate how far a thermal wave has traveled, confirming whether the semi-infinite assumption holds for your specific material thickness and process duration.

Engineering Applications of the Model

The semi-infinite solution is not just theoretical; it directly applies to numerous engineering fields where thick bodies experience short-term thermal events.

• Ground Temperature Response: In geothermal engineering and building science, the temperature variation of the Earth's surface due to daily or seasonal cycles is modeled using this approach. The soil acts as a semi-infinite medium for diurnal cycles because the thermal wave from the sun only penetrates a meter or two. The model helps design ground-source heat pumps and assess thermal storage in the shallow subsurface.
• Quenching of Thick Parts: In metallurgy, rapidly cooling a hot metal part by immersion in a fluid (quenching) is critical for achieving desired material properties. During the initial seconds, the core of a thick forging remains hot while the surface cools abruptly. Using the semi-infinite model, you can estimate the temperature gradient and cooling rate at different depths, which informs predictions about residual stress and hardening.
• Short-Duration Heating of Walls and Slabs: In fire safety engineering, assessing the temperature rise in a concrete wall during a brief fire exposure often relies on this model. For short times, the heat from the fire hasn't reached the interior side of the wall, validating the semi-infinite assumption. This analysis is key for determining fire resistance ratings and structural integrity under thermal load.

Consider a step-by-step example: A thick steel slab ($\alpha =1.2\times 10^{-5}{\text{m}}^2/\text{s}$) initially at $300^{\circ }\text{C}$ has its surface suddenly cooled to $50^{\circ }\text{C}$. To find the temperature 2 cm deep after 30 seconds:

1. Calculate the similarity variable: $\eta =\frac{0.02\text{m}}{2\sqrt{(1.2\times 10^{-5})\times 30}}\approx 0.527$.
2. Look up erfc(0.527) from a table or calculator: erfc(0.527) $\approx$ 0.416.
3. Apply the solution: $T=300+(50-300)\times 0.416=300-250\times 0.416\approx 196^{\circ }\text{C}$.

This quick calculation shows the interior is still relatively hot, confirming the short-time behavior.

Common Pitfalls

1. Applying the Model When Boundaries Are Not Semi-Infinite: The most frequent error is using this solution for times long enough that the thermal disturbance reaches the far boundary of the actual object. For a slab of finite thickness $L$, the model is only valid for times approximately $t<L^2/(16\alpha )$. Beyond this, you must use solutions for finite geometries. Always check the penetration depth estimate against the physical thickness.

• Correction: Before applying the solution, verify that $4\sqrt{\alpha t}\ll L$, where $L$ is the material thickness. If not, switch to a transient solution for a plane wall.

1. Confusing erf and erfc: The complementary error function (erfc) is used for the constant surface temperature problem, while the error function (erf) appears in solutions for other boundary conditions, like a constant surface heat flux. Using the wrong function will invert your temperature profile.

• Correction: Remember the mnemonic: For a changed surface temperature, you use the complementary error function (erfc). The solution form $T=T_i+(T_s-T_i)\text{erfc}(\eta )$ is standard for this specific case.

1. Ignoring the Assumption of Constant Properties: The solution assumes thermal diffusivity $\alpha$ is constant. In reality, for large temperature ranges, properties like conductivity and heat capacity change, making the problem nonlinear.

• Correction: For processes involving extreme temperatures (e.g., quenching from very high heat), use property values averaged over the expected temperature range or employ numerical methods that account for property variations.

1. Misinterpreting the Surface Condition: The classic solution requires a sudden, step-change in surface temperature. In practice, the surface temperature might change gradually, or the boundary condition might involve convection (e.g., quenching in a fluid bath).

• Correction: For a convective boundary condition, you can still use a semi-infinite model but with a slightly more complex solution involving erfc and an additional parameter, the Biot number. Ensure you select the solution that matches your actual surface thermal interaction.

Summary

• The semi-infinite solid model is a vital tool for analyzing transient heat conduction in thick bodies where the thermal pulse has not yet reached the interior boundary, valid for short-time or large-thickness scenarios.
• The temperature distribution is governed by the complementary error function solution $T(x,t)=T_i+(T_s-T_i)\text{erfc}(x/(2\sqrt{\alpha t}))$, linking depth and time through the similarity variable and material thermal diffusivity $\alpha$.
• Key applications include predicting ground temperature response to environmental cycles, analyzing quenching gradients in metallurgy, and evaluating short-duration heating in structures for fire safety.
• Always verify that the thermal penetration depth $4\sqrt{\alpha t}$ is much less than the actual object thickness to ensure the semi-infinite assumption holds.
• Avoid common errors by correctly identifying boundary conditions, distinguishing between erf and erfc, and using constant material properties only when justified by the temperature range.