臨界断熱半径
Theory and Physics
What is Critical Insulation Radius?
Professor, I thought that adding insulation always reduces heat loss, but is it true that it can sometimes increase it?
It's true. When insulation is wrapped around the outside of a cylinder or sphere, the thermal resistance due to conduction increases, but the external surface area also increases. Since convective thermal resistance is inversely proportional to area, up to a certain radius, the decrease in convective resistance outweighs the increase in conductive resistance, and the overall heat loss actually increases.
That's counterintuitive.
The radius at which heat loss is maximized is called the critical insulation radius $r_{cr}$. For a cylinder, the total thermal resistance is
Differentiating this with respect to $r$ and setting it to zero gives
For a sphere, it becomes $r_{cr} = 2k/h$.
So the smaller $k$ is and the larger $h$ is, the smaller the critical radius becomes.
Correct. For natural convection ($h \approx 5$ W/(m$^2$ K)) and epoxy insulation ($k \approx 0.2$ W/(m K)), $r_{cr} = 0.04$ m = 40 mm. If the original pipe outer diameter is less than 40mm, there exists a region where adding insulation will have the opposite effect.
Physical Interpretation
It's easier to understand by breaking down the total thermal resistance.
| Radius $r$ | Conduction Resistance | Convection Resistance | Total Resistance | Heat Loss |
|---|---|---|---|---|
| $r < r_{cr}$ | Increase (small) | Decrease (large) | Decrease | Increase |
| $r = r_{cr}$ | — | — | Minimum | Maximum |
| $r > r_{cr}$ | Increase (large) | Decrease (small) | Increase | Decrease |
So it's a battle between whether the decrease in convection resistance outweighs the increase in conduction resistance.
Exactly. In practice, this becomes an issue with wire insulation and pipe insulation. During design, you must confirm that $r_i > r_{cr}$ before determining the insulation thickness.
Intuitive Meaning of rcrit = λ/h
The critical insulation radius for cylindrical insulation, rcrit = λins/ho, is the point where the rate of increase in the insulation's conductive thermal resistance balances the rate of increase in the external surface's convective heat transfer. For typical outdoor conditions (ho = 10 W/m²·K) and glass wool (λ=0.04 W/m·K), rcrit = 4 mm. The outer radius of most industrial pipes far exceeds this, so it's not a practical issue, but it becomes important in the design of wire insulation and medical tubing.
Physical Meaning of Each Term
- Heat Storage Term $\rho c_p \partial T/\partial t$: Rate of thermal energy storage per unit volume. 【Everyday Example】 An iron frying pan is slow to heat up and cool down, while an aluminum pot heats and cools quickly—this is due to the difference in the product of density $\rho$ and specific heat $c_p$ (Heat Capacity). Objects with large heat capacity have slower temperature changes. Water has a very high specific heat (4,186 J/(kg·K)), which is why temperatures near the ocean are more stable than inland. In transient analysis, this term determines the rate of temperature change over time.
- Heat Conduction Term $\nabla \cdot (k \nabla T)$: Heat conduction based on Fourier's law. Heat flux proportional to the temperature gradient. 【Everyday Example】 Putting a metal spoon in a hot pot makes the handle hot—metal has a high thermal conductivity $k$, so heat transfers quickly from the hot side to the cold side. A wooden spoon doesn't get hot because its $k$ is small. Insulation materials (like glass wool) have extremely small $k$, making heat transfer difficult even with a temperature gradient. This term mathematically expresses the natural tendency of heat to flow where there is a temperature difference.
- Convection Term $\rho c_p \mathbf{u} \cdot \nabla T$: Heat transport accompanying fluid motion. 【Everyday Example】 Feeling cool in front of a fan is because the wind (fluid flow) carries away the warm air near your skin and supplies fresh, cold air—this is forced convection. The ceiling area of a room becoming warm with heating is due to natural convection where heated air rises due to buoyancy. The fan on a PC's CPU cooler also dissipates heat via forced convection. Convection is an order of magnitude more efficient at heat transport than conduction.
- Heat Source Term $Q$: Internal heat generation (Joule heat, chemical reaction heat, radiation absorption, etc.). Unit: W/m³. 【Everyday Example】 A microwave oven heats food via microwave absorption inside the food (volumetric heating). The heater wire in an electric blanket warms up via Joule heating ($Q = I^2 R / V$). Heat generation during lithium-ion battery charging/discharging and friction heat from brake pads are also considered as heat sources in analysis. Unlike boundary conditions that supply heat to the "surface" from the outside, the heat source term represents energy generation "inside" the material.
Assumptions and Applicability Limits
- Fourier's Law: Linear relationship where heat flux is proportional to temperature gradient (non-Fourier heat conduction is needed for extremely low temperatures or ultra-short pulse heating)
- Isotropic Thermal Conductivity: Thermal conductivity is independent of direction (anisotropy must be considered for composite materials or single crystals)
- Temperature-Independent Material Properties (Linear Analysis): Assumption that material properties do not depend on temperature (temperature dependence is needed for large temperature differences)
- Treatment of Thermal Radiation: Surface-to-surface radiation uses the view factor method; for participating media, the DO method or P1 approximation is applied
- Non-Applicable Cases: For phase change (melting/solidification), latent heat must be considered. For extreme temperature gradients, thermal-stress coupling is essential
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Temperature $T$ | K (Kelvin) or Celsius | Be careful not to confuse absolute and Celsius temperatures. Always use absolute temperature for radiation calculations. |
| Thermal Conductivity $k$ | W/(m·K) | Steel: ~50, Aluminum: ~237, Air: ~0.026 |
| Heat Transfer Coefficient $h$ | W/(m²·K) | Natural Convection: 5–25, Forced Convection: 25–250, Boiling: 2,500–25,000 |
| Specific Heat $c_p$ | J/(kg·K) | Distinguish between specific heat at constant pressure and constant volume (important for gases) |
| Heat Flux $q$ | W/m² | Neumann condition as a boundary condition |
Numerical Methods and Implementation
Analytical Solution
The critical insulation radius problem has a neat analytical solution, doesn't it?
Yes. An exact solution is obtained for steady-state heat conduction in 1D cylindrical coordinates. The temperature distribution is
The heat loss $q$ is
The critical radius is derived by treating $r_o$ as a variable, differentiating $q$ with respect to $r_o$, and setting it to zero.
Verification via Numerical Analysis
Does FEM give the same result?
Of course it does. In fact, comparison with the theoretical solution is used as a benchmark for solver verification. In Ansys Mechanical, you can create a cylindrical mesh with SOLID70, apply a temperature constraint on the inner surface, and set a convection condition on the outer surface.
| Parameter | Value |
|---|---|
| Inner Radius $r_i$ | 5 mm |
| Insulation $k$ | 0.2 W/(m K) |
| External $h$ | 10 W/(m$^2$ K) |
| Theoretical $r_{cr}$ | 20 mm |
Performing a parametric analysis by varying insulation thickness from 5mm to 50mm confirms that heat loss is maximized at $r_o = 20$ mm. The error compared to the theoretical value is less than 0.1%.
It's a perfect problem for verification.
You can automate all cases with a DO loop in an APDL macro. Similarly, in Abaqus, you can create a study with $r_o$ as a parameter using Python scripting.
When Including Temperature Dependence
When the insulation's $k$ is temperature-dependent, the question becomes at which temperature to evaluate $k$ in $r_{cr} = k(T)/h$. A self-consistent solution must be found using Newton-Raphson iteration.
In reality, insulation $k$ often increases with temperature, right?
Glass wool becomes about twice its room temperature value at 200°C. For high-temperature pipe insulation design, ignoring temperature dependence can lead to underdesign.
Visual Confirmation by Plotting Heat Loss Curves
To truly grasp the existence of the critical radius, plot heat loss Q on the vertical axis against insulation outer radius on the horizontal axis. The critical point can be found by differentiating Q(r)=2πLΔT/[ln(r/ri)/λ + 1/(h·r)] and setting dQ/dr=0. The heat transfer textbook by Incropera & DeWitt from the 1970s (now in its 7th edition) popularized this graph as a standard example problem in thermal engineering education worldwide.
Linear Elements vs. Quadratic Elements
In heat conduction analysis, linear elements often provide sufficient accuracy. For regions with steep temperature gradients (e.g., thermal shock), quadratic elements are recommended.
Heat Flux Evaluation
Calculated from the temperature gradient within an element. Smoothing may be required, similar to nodal stresses.
Convection-Diffusion Problems
When the Peclet number is high (convection-dominated), upwinding stabilization (e.g., SUPG) is needed. Not required for pure heat conduction problems.
Time Step for Transient Analysis
Set a time step sufficiently smaller than the characteristic thermal diffusion time $\tau = L^2 / \alpha$ ($\alpha$: Thermal Diffusivity). Automatic time step control is effective for rapid temperature changes.
Nonlinear Convergence
Nonlinearity due to temperature-dependent material properties is often mild, and Picard iteration (direct substitution method) is often sufficient. Newton's method is recommended for the strong nonlinearity of radiation.
Steady-State Analysis Convergence Criterion
Convergence is judged when the temperature change at all nodes falls below a threshold (e.g., $|\Delta T| / T_{max} < 10^{-5}$).
Analogy for Explicit and Implicit Methods
Explicit method is like "predicting the next step using only current information, like a weather forecast"—calculation is fast but unstable with large time steps (misses storms). Implicit method is like "prediction considering future states"—stable even with large time steps but requires solving equations at each step. For problems without rapid temperature changes, using the implicit method with larger time steps is more efficient.
Practical Guide
Application to Design
How is the critical insulation radius used in practical work?
There are two main scenarios.
1. Insulation Design: When determining insulation thickness for pipes or ducts, confirm that $r_i > r_{cr}$ before optimizing thickness.
2. Heat Dissipation Design: When determining coating thickness for electrical wires, if $r_i < r_{cr}$, then increasing coating thickness improves heat dissipation.
Using it inversely for heat dissipation design is interesting.
Wire insulation is a classic example. For an AWG24 copper wire (outer diameter 0.56mm) with PVC coating ($k = 0.16$ W/(m K)), with natural convection $h = 10$ W/(m$^2$ K), $r_{cr} = 16$ mm. Heat dissipation improves as the coating is thickened until the coated outer diameter reaches 32mm.
Properties of Typical Insulation Materials
| Insulation Material | $k$ [W/(m K)] | Temperature Range | $r_{cr}$ ($h$=10) |
|---|---|---|---|
| Glass Wool | 0.04 | Up to 450°C | 4 mm |
| Rock Wool | 0.04 | Up to 700°C | 4 mm |
| Polyurethane Foam | 0.02 | Up to 100°C | 2 mm |
| Silica Aerogel | 0.015 | Up to 650°C | 1.5 mm |
| Ceramic Fiber | 0.08 | Up to 1200°C | 8 mm |
Higher-performance insulation materials have smaller $r_{cr}$, so they can be used safely even on thin pipes.
Aerogel has $r_{cr} = 1.5$ mm, so for pipes with an outer diameter above 3mm, it's hardly ever a problem. It's expensive but has a track record in aerospace and LNG fields.
Result Verification
Key points for verifying critical insulation radius analysis results are as follows.
- Comparison with Theoretical Value: Does it match $r_{cr} = k/h$ (cylinder) or $2k/h$ (sphere)?
- Heat Loss Curve: Does it peak at $r = r_{cr}$ and change monotonically before and after?
- Energy Balance: Does the heat input at the inner surface match the heat loss at the outer surface?
I see. It's important to check both the peak position and the overall trend.
Related Topics
なった
詳しく
報告