🔥 Semi-Infinite Solid Conduction

Calculate transient heat penetration, temperature profiles, and surface heat rates in semi-infinite media under constant temperature, heat flux, or convection boundary conditions.

⚡ Fortran 90 Engine Double Precision (IEEE 754) ✓ ISO / ASME Validated
Semi-Infinite Solid Conduction Heat Transfer
📊 Solver Telemetry ● ACTIVE
👁️ Views 155
⚡ Solves 123
💾 Downloads 481 📦 Fortran Code 3.4 KB
📅 Released Jun 2026
⏱️ Latency < 1 ms
⚡ TOOLS & REPORTS:
💾 Download Fortran 90

Configuration

Boundary Condition

Material Properties

Coordinates & Simulation Time

3600 s

Results & Visualization

Thermal Penetration Heat Field Map

T(x) x = 0 x_max ← Heat Conduction Path →

Transient Profile: T(x) at Multiple Times

Current Performance Summary

Penetration Depth
-
Surface Temp. T(0,t)
-
Surface Heat Flux
-
Temp. at depth T(x,t)
-
Total Energy Absorbed (Q)
-

Detailed text reports will appear here after clicking "Generate Engineering Report".

Methodology & Theory

Transient Thermal Penetration

A semi-infinite solid is an idealized body with a single plane surface, extending to infinity. This model represents real geometries (such as the earth's crust, thick concrete walls, or large metallic shafts) during the early stages of a thermal transient before the heat wave reaches the opposite boundaries.

The depth to which the temperature change propagates is defined as the **penetration depth**: $$\delta_p \approx 4\sqrt{\alpha t}$$ Beyond $\delta_p$, the temperature of the solid remains practically undisturbed at its initial value ($T_i$).

Academic Reference Formulas:

  • Case 1 — Constant $T_s$: $$T(x,t) = T_s + (T_i - T_s)\,\text{erf}\left(\frac{x}{2\sqrt{\alpha t}}\right)$$ $$Q(t) = 2\,k\,(T_s - T_i)\,\sqrt{\frac{t}{\pi \alpha}}$$
  • Case 2 — Constant $q''_s$: $$T(x,t) = T_i + \frac{2q''_s\sqrt{\alpha t/\pi}}{k}\exp\left(-\frac{x^2}{4\alpha t}\right) - \frac{q''_s x}{k}\text{erfc}\left(\frac{x}{2\sqrt{\alpha t}}\right)$$ $$Q(t) = q''_s \cdot t$$
  • Case 3 — Convection Boundary: $$T(x,t) = T_i + (T_\infty - T_i)\left[\text{erfc}\left(\frac{x}{2\sqrt{\alpha t}}\right) - \exp\left(\frac{hx}{k}+\frac{h^2\alpha t}{k^2}\right)\text{erfc}\left(\frac{x}{2\sqrt{\alpha t}}+\frac{h\sqrt{\alpha t}}{k}\right)\right]$$

Academic References:

  1. Incropera, F. P., DeWitt, D. P., Bergman, T. L., & Lavine, A. S. (2011). Fundamentals of Heat and Mass Transfer (7th Edition). John Wiley & Sons. Chapter 5.7 (The Semi-Infinite Solid).
  2. Carslaw, H. S., & Jaeger, J. C. (1959). Conduction of Heat in Solids (2nd Edition). Oxford University Press.
  3. Çengel, Y. A., & Ghajar, A. J. (2015). Heat and Mass Transfer: Fundamentals and Applications (5th Edition). McGraw-Hill. Chapter 4.3 (Transient Heat Conduction in Semi-Infinite Solids).

Worked Engineering Example

Problem Statement:
A thick concrete slab ($k = 1.4$ W/m·K, thermal diffusivity $\alpha = 7.0 \cdot 10^{-7}$ m²/s) is initially at a uniform temperature of $T_i = 20$°C. The surface is suddenly exposed to hot air at $T_\infty = 200$°C with a convection heat transfer coefficient of $h = 100$ W/m²·K. Calculate the temperature at a depth of $x = 50$ mm after $t = 1$ hour ($3600$ s) of exposure.

Step-by-step Solution:
1. Identify parameters and compute transient limits:
$$x = 0.05 \text{ m}, \quad t = 3600 \text{ s}, \quad k = 1.4 \text{ W/m·K}, \quad \alpha = 7.0 \cdot 10^{-7} \text{ m²/s}$$ $$\sqrt{\alpha t} = \sqrt{7.0 \cdot 10^{-7} \times 3600} = 0.0502 \text{ m}$$ $$\text{Dimensionless depth: } \eta = \frac{x}{2\sqrt{\alpha t}} = \frac{0.05}{2 \times 0.0502} = 0.498$$ $$\text{Convection parameter: } \frac{h\sqrt{\alpha t}}{k} = \frac{100 \times 0.0502}{1.4} = 3.586$$ 2. Calculate complementary error functions:
$$\text{erfc}(\eta) = \text{erfc}(0.498) \approx 0.482$$ 3. Evaluate Case 3 Convection boundary equation:
$$\frac{T(x,t) - T_i}{T_\infty - T_i} = \text{erfc}(\eta) - \exp\left(\frac{hx}{k} + \frac{h^2 \alpha t}{k^2}\right) \text{erfc}\left(\eta + \frac{h\sqrt{\alpha t}}{k}\right)$$ $$\frac{hx}{k} + \frac{h^2\alpha t}{k^2} = \frac{100 \times 0.05}{1.4} + (3.586)^2 = 3.571 + 12.86 = 16.43$$ $$\text{erfc}(0.498 + 3.586) = \text{erfc}(4.084) \approx 6.4 \cdot 10^{-9}$$ $$\exp(16.43) \times \text{erfc}(4.084) \approx 1.36 \cdot 10^7 \times 6.4 \cdot 10^{-9} = 0.087$$ $$\frac{T(x,t) - T_i}{T_\infty - T_i} \approx 0.482 - 0.087 = 0.395$$ 4. Find the local temperature $T(0.05\text{m}, 1\text{hr})$:
$$T(x,t) = T_i + 0.395 \times (T_\infty - T_i) = 20 + 0.395 \times (200 - 20) = 91.1 \text{ °C}$$
Final Result:
The concrete temperature at $50$ mm depth is 91.1 °C.

📘 Calculation Methodology: Semi-Infinite Solid Transient Heat Conduction

Mathematical Model & Theory

When a semi-infinite solid ($x \ge 0$) undergoes a sudden surface temperature change $T_s$, the temperature profile is governed by the Gaussian error function:

$$\frac{T(x,t) - T_s}{T_i - T_s} = \text{erf}\left(\frac{x}{2\sqrt{\alpha t}}\right), \quad q''_s(t) = \frac{k(T_s - T_i)}{\sqrt{\pi \alpha t}}$$
$$\text{Penetration Depth: } \delta_{th} \approx 3.64 \sqrt{\alpha t}$$

Assumptions

  • 1D transient conduction in semi-infinite domain.
  • Homogeneous material with constant diffusivity $lpha = k/( ho c_p)$.

Academic References

  1. Incropera, F. P. et al.: Fundamentals of Heat and Mass Transfer, Ch. 5.
  2. Holman, J. P.: Heat Transfer, McGraw-Hill.

Worked Engineering Example

Problem Statement:
A concrete wall ($\alpha = 7.0 \times 10^{-7}\text{ m}^2/\text{s}$) at $T_i = 20^\circ\text{C}$ is exposed to $T_s = 600^\circ\text{C}$. Calculate temperature at $x = 50\text{ mm}$ after $1\text{ hour}$ ($3600\text{ s}$).

Step-by-step Solution:
1. $2\sqrt{\alpha t} = 2\sqrt{7.0 \times 10^{-7} \times 3600} \approx 0.1004\text{ m}$.
2. $\eta = 0.05 / 0.1004 \approx 0.498 \implies \text{erf}(0.498) \approx 0.5205$.
3. $T = 600 + 0.5205 \times (20 - 600) = 600 - 301.9 = 298.1^\circ\text{C}$.
Final Result:
Temperature at 50 mm depth is $\mathbf{298.1^\circ\text{C}}$.