臨界断熱厚
Theory and Physics
Concept of Critical Insulation Thickness
Professor, are critical insulation thickness and critical insulation radius the same thing?
They are essentially the same phenomenon. The critical insulation thickness $t_{cr}$ is the thickness at which the heat dissipation becomes maximum when insulation is wrapped around the outer surface of a cylinder, expressed as $t_{cr} = r_{cr} - r_i = k/h - r_i$.
Flat plates don't have a critical insulation thickness, right?
Correct. For a flat plate, adding more insulation does not change the outer surface area, so the total thermal resistance increases monotonically. This phenomenon occurs only for cylinders and spheres where the area changes.
Total Thermal Resistance of a Cylinder
The total thermal resistance as a function of insulation thickness $t$ is
Setting $dR/dt = 0$ yields $r_i + t = k/h$. Thus, $t_{cr} = k/h - r_i$. If $r_i > k/h$, then $t_{cr} < 0$, meaning the critical insulation thickness issue does not arise.
Case of a Sphere
In spherical coordinates, the total thermal resistance is
The critical radius is $r_{cr} = 2k/h$, which is twice that of a cylinder.
Is the critical radius larger for a sphere because the way the area increases is different?
Yes. The surface area of a sphere is $4\pi r^2$, proportional to the square of the radius, while the lateral area of a cylinder $2\pi r L$ is proportional to the radius to the first power. The area increase effect is greater for a sphere, hence the larger critical radius.
Design Considerations
| Shape | $r_{cr}$ | Practical Implication |
|---|---|---|
| Flat Plate | None | Adding insulation is always effective |
| Cylinder | $k/h$ | Caution for small-diameter pipes/wires |
| Sphere | $2k/h$ | Caution for insulation design of spherical vessels |
Adding Insulation Can Be Counterproductive
For flat plate insulation, increasing thickness always reduces heat loss. However, for insulation wrapped around a cylindrical pipe, the increased outer surface area enhances convection, creating a paradox where heat loss actually increases until the "critical radius" rcrit = λ/h is exceeded. This phenomenon was first mathematically demonstrated by Lord Rayleigh in an 1896 paper discussing heat radiation from cylinders.
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 heats up slowly and cools down slowly, while an aluminum pot heats up quickly and cools down 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 experience 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】A metal spoon placed in a hot pot gets hot all the way to the handle—metals have high thermal conductivity $k$, so heat transfers quickly from the high-temperature side to the low-temperature side. A wooden spoon doesn't get hot because its $k$ is small. Insulation materials (e.g., glass wool) have extremely small $k$, making heat transfer difficult even with a temperature gradient. This term mathematically expresses the natural tendency of "heat flowing 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 under a fan is because the wind (fluid flow) carries away warm air near the body surface and supplies fresh, cool 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 in a PC's CPU cooler also uses forced convection for heat dissipation. Convection is an order of magnitude more efficient heat transport method 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 through microwave absorption inside the food (volumetric heating). The heater wire in an electric blanket warms up through Joule heating ($Q = I^2 R / V$). Heat generation during charging/discharging of lithium-ion batteries and friction heat from brake pads are also considered as heat sources in analysis. Unlike boundary conditions where heat is supplied from the "surface" externally, 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, etc.)
- 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: Phase change (melting/solidification) requires consideration of latent heat. Extreme temperature gradients necessitate thermal-stress coupling
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Temperature $T$ | K (Kelvin) or Celsius | Be careful not to confuse absolute temperature and Celsius. 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
Choosing a Calculation Method
What is an appropriate method to calculate the critical insulation thickness?
Choose based on the complexity of the problem.
| Condition | Manual Calculation | Numerical Analysis |
|---|---|---|
| Simple cylinder, constant $k$, constant $h$ | Sufficient | Unnecessary |
| Temperature-dependent $k(T)$ | Approximately possible | Recommended |
| Includes radiation | Approximation via linearization | Recommended |
| Non-circular cross-section | Not possible | Required |
| Multi-layer insulation | Possible but cumbersome | Recommended |
Exact Calculation of Heat Dissipation
Heat dissipation per unit length as a function of insulation thickness:
Creating a table of this in Excel makes it easy to identify the critical point. You can find the maximum value numerically without taking the derivative.
So Excel is sufficient.
For basic cases, Excel's Goal Seek or Solver is sufficient. However, if temperature dependence is included, using Python's SciPy.optimize or FEM is more reliable.
Verification with FEM
Here is the verification procedure in Ansys Mechanical.
1. Create a 2D axisymmetric model (PLANE55, KEYOPT(3)=1)
2. Model an annular cross-section with inner radius $r_i$
3. Vary the outer radius as a parameter using an APDL DO loop
4. Obtain the Total Heat Flow for each case
5. Compare with the theoretical value $r_{cr} = k/h$
What about meshing for thin insulation layers?
At least 3 elements in the radial direction. For an insulation thickness of 1mm, an element size of about 0.3mm is sufficient. Aim for an element aspect ratio of 5 or less.
Economic Calculation of Optimal Insulation Thickness
The optimal insulation thickness is an economic optimization problem determined by "incremental insulation cost = reduction in energy cost due to reduced heat loss." The piping insulation design standard (JSME S010) established by the Japan Society of Mechanical Engineers in the 1980s specifies a procedure for calculating thickness based on Life Cycle Cost (LCC) minimization, and it is still used as a standard in chemical plant design.
Linear Elements vs. Quadratic Elements
In heat conduction analysis, linear elements often provide sufficient accuracy. For areas 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 Problem
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 small relative to 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 strong nonlinearities like radiation.
Steady-State Analysis Determination
Convergence is determined 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
The explicit method is like "a weather forecast predicting the next step using only current information"—fast to compute but unstable with large time steps (misses storms). The implicit method is like "a prediction that also considers 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 Scenarios in Practice
Please tell me about cases where critical insulation thickness actually becomes a problem.
Let me list some typical cases.
Electric Wire Sheathing Design
For an AWG24 copper wire (outer diameter 0.56mm) with PVC sheathing ($k = 0.16$ W/(m K)), with natural convection $h = 10$ W/(m$^2$ K), $r_{cr} = 16$ mm. Heat dissipation improves as the sheathing is thickened until the outer diameter reaches 32mm.
Getting cooler by thickening the sheathing is counterintuitive.
The allowable current calculation for power cables (IEC 60287) considers the heat dissipation improvement effect of the sheathing. However, practical sheathing thickness is determined by mechanical protection and insulation requirements.
Insulation for Small-Diameter Refrigerant Pipes
For a refrigerant pipe (outer diameter 6.35mm) with rubber-based insulation ($k = 0.04$), $r_{cr} = 4$ mm. This is close to the pipe radius of 3.2mm, so thin insulation has little effect.
| Insulation Thickness [mm] | Heat Dissipation Ratio | Surface Temperature |
|---|---|---|
| 0 | 1.00 | Low |
| 1 | 1.01 | Slightly lower |
| 5 | 0.92 | Begins to rise |
| 15 | 0.70 | Sufficient insulation effect |
So you need at least 5mm to see an effect.
In practice, insulation thickness of 10mm or more is common, so it's not a problem. However, it's important not to make the simplistic judgment that "thin insulation is sufficient."
Economic Insulation Thickness
In practice, the optimal thickness is determined by comparing the investment cost of insulation with energy savings payback. JIS A 9501 "Standard for Thermal Insulation Work" specifies a method for calculating the economic thickness.
So it's a different optimization problem from critical insulation thickness.
Critical insulation thickness is a physical limit, while economic insulation thickness is a cost optimization. Design decisions should be made understanding both.
Insulation Standards for Process Plants
For steam pipes common in petrochemical plants (150°C, 100 mm diameter), glass wool insulation is typically designed with a thickness of 50–75 mm. JIS A 9504 (2018 edition) specifies minimum thicknesses for different operating temperatures. Falling below this can lead to high surface temperatures, violating the "surfaces that may be touched should be below 60°C" clause of the Industrial Safety and Health Act.
Analogy for Analysis Flow
Think of the thermal analysis flow as "designing a bath reheating system." Decide the bathtub shape (analysis target), set the initial water temperature (initial conditions) and outside air temperature (boundary conditions), and adjust the reheater output (heat source). Predicting "will it become lukewarm after 2 hours?" through calculation—this is the essence of transient thermal analysis.
Common Pitfalls for Beginners
"Can I ignore radiation?" — Usually OK around room temperature. But for high temperatures...
Related Topics
なった
詳しく
報告