🌐 Anisotropic Thermal Conduction

Analyze directional conductivity tensor rotation for composites, PCBs, and wood. Calculates effective thermal conductivity and directional heat fluxes.

⚡ Fortran 90 Engine Double Precision (IEEE 754) ✓ ISO / ASME Validated
Anisotropic Thermal Conduction Heat Transfer
📊 Solver Telemetry ● ACTIVE
👁️ Views 106
⚡ Solves 85
💾 Downloads 415 📦 Fortran Code 3.4 KB
📅 Released Jun 2026
⏱️ Latency < 1 ms
⚡ TOOLS & REPORTS:
💾 Download Fortran 90

Parameters Input

Directional Tensor Rotation: $$k_{eff}(\theta) = k_x\cos^2\theta + k_y\sin^2\theta$$ $$q_x = k_{eff}\frac{T_{hot} - T_{cold}}{L_x}$$ $$Q_x = q_x A$$

Results & Analysis

Results and visualizations will be displayed here upon completion of the computation.

Calculation Methodology

Mathematical Model & Theory

Fourier's Law for heat conduction in an isotropic solid links heat flux directly to temperature gradient. However, in anisotropic materials (like composites, wood, or graphitic crystals), conductivity is a tensor quantity. The directional heat flux in coordinate systems rotated relative to principal material coordinates in 2D is resolved as:

$$k_{eff}(\theta) = k_x \cos^2(\theta) + k_y \sin^2(\theta)$$ $$k_{xy}' = (k_x - k_y)\sin(\theta)\cos(\theta)$$

Variable Definitions & Units:

  • $k_x, k_y, k_z$: Thermal conductivity along principal axes [W/m·K]
  • $\theta$: Fiber or layer rotation angle [rad or deg]
  • $k_{eff}$: Rotated effective thermal conductivity [W/m·K]
  • $q_x$: Thermal heat flux along x-axis [W/m²]
  • $Q_x$: Heat transfer rate [W]

Academic References:

  1. Incropera, F. P., & DeWitt, D. P. (2011). Fundamentals of Heat and Mass Transfer. Wiley.
  2. Özisik, M. N. (1993). Heat Conduction. Wiley.
  3. Carslaw, H. S., & Jaeger, J. C. (1959). Conduction of Heat in Solids. Oxford University Press.

Worked Engineering Example

Problem Statement:
A Carbon Fiber Reinforced Polymer (CFRP) sheet ($k_x = 7.0$ W/m·K along fibers, $k_y = 0.8$ W/m·K perpendicular to fibers) is 3 mm thick ($L_x = 0.003$ m) and has an area of 0.001 m². The hot face is at 120°C and the cold face is at 25°C. If the fibers are rotated by 30° relative to the thermal system, calculate the effective conductivity and heat rate.

Step-by-step Solution:
1. Calculate the effective rotated conductivity $k_{eff}$:
$$k_{eff} = 7.0\cos^2(30^\circ) + 0.8\sin^2(30^\circ)$$ $$k_{eff} = 7.0(0.75) + 0.8(0.25) = 5.25 + 0.20 = 5.45 \text{ W/m·K}$$ 2. Calculate directional thermal heat flux $q_x$:
$$q_x = k_{eff} \frac{T_{hot} - T_{cold}}{L_x} = 5.45 \frac{120 - 25}{0.003} = 172,583.3 \text{ W/m}^2$$ 3. Calculate heat transfer rate $Q_x$:
$$Q_x = q_x A = 172,583.3 \times 0.001 = 172.58 \text{ W}$$
Final Result:
Effective rotated conductivity is 5.45 W/m·K and total heat rate is 172.6 W.