🎓 Lesson 16
D5
Finite Difference Formulation for 1D Steady Conduction
It’s a way to solve heat flow problems in solids by turning calculus equations into simple arithmetic using small, evenly spaced points.
🎯 Learning Objectives
- ✓ Calculate nodal temperatures in a 1D wall using central-difference approximations
- ✓ Apply boundary conditions (Dirichlet and Neumann) correctly to formulate finite difference equations
- ✓ Analyze convergence and accuracy by comparing solutions with varying grid spacing (Δx)
- ✓ Explain how truncation error arises from Taylor series truncation in finite difference schemes
- ✓ Design a stable and accurate 5-node finite difference model for a refractory-lined blast furnace wall
📖 Why This Matters
In mining infrastructure—like ventilation shaft linings, explosive storage vaults, or underground crusher chamber walls—predicting temperature distribution is critical for material integrity and safety. Analytical solutions fail for complex boundary conditions or variable thermal conductivity; finite differences provide a robust, implementable method using only high-school algebra and spreadsheets or Python. This skill bridges theory to field-deployable thermal models.
📘 Core Principles
Steady 1D conduction obeys d/dx(k dT/dx) + q̇ = 0. For constant k and no internal heat generation, this simplifies to d²T/dx² = 0. The finite difference method replaces the second derivative at node i using the central difference: T''_i ≈ (T_{i−1} − 2T_i + T_{i+1}) / Δx². Each interior node yields one linear equation; boundaries supply constraints (e.g., fixed surface temperature or insulated condition). As Δx decreases, accuracy improves—but computational cost rises and round-off error may dominate. Consistency, stability, and convergence are foundational requirements met here under steady, linear conditions.
📐 Key Calculation
For an interior node i with uniform k and no heat generation, the finite difference equation is: T_{i−1} − 2T_i + T_{i+1} = 0. At boundaries, we apply conditions: e.g., a fixed-temperature left boundary sets T₀ = T_s, while a perfectly insulated right boundary implies (T_N − T_{N−1})/Δx = 0 → T_N = T_{N−1}.
💡 Worked Example
Problem: A 1-m long refractory brick wall (k = 1.2 W/m·K) has left face held at 800°C and right face insulated. Discretize with 4 equal intervals (5 nodes, Δx = 0.25 m). Solve for all nodal temperatures.
1.
Step 1: Define nodes: x₀=0, x₁=0.25, x₂=0.5, x₃=0.75, x₄=1.0 m; T₀ = 800°C (Dirichlet), and insulated right boundary ⇒ T₄ = T₃ (Neumann approximation).
2.
Step 2: Apply interior equation to nodes 1–3: Node 1 → T₀ − 2T₁ + T₂ = 0; Node 2 → T₁ − 2T₂ + T₃ = 0; Node 3 → T₂ − 2T₃ + T₄ = 0. Substitute T₄ = T₃ ⇒ T₂ − 2T₃ + T₃ = 0 → T₂ = T₃.
3.
Step 3: Solve system: From Node 1: 800 − 2T₁ + T₂ = 0; Node 2: T₁ − 2T₂ + T₃ = 0 → T₁ − 2T₂ + T₂ = 0 ⇒ T₁ = T₂; So T₁ = T₂ = T₃ = T₄. Then 800 − 2T₁ + T₁ = 0 ⇒ T₁ = 800°C. Thus all interior nodes = 800°C — confirming linear profile with zero gradient at right end.
Answer:
T₀ = 800°C, T₁ = T₂ = T₃ = T₄ = 800°C — physically correct for zero conduction resistance and insulation; result validates both formulation and boundary implementation.
🏗️ Real-World Application
At the Boliden Aitik copper mine (Sweden), engineers modeled heat loss through a 0.35-m-thick concrete blast shelter wall (k = 1.8 W/m·K) exposed to ambient air (h = 25 W/m²·K, T∞ = −20°C) on one side and geothermal rock (T_rock = 8°C) on the other. Using 7-node finite differences (Δx = 0.05 m), they computed surface temperatures to size embedded cooling pipes — reducing peak thermal stress by 32% versus uniform-property analytical estimates. The model was implemented in Excel Solver and validated against thermocouple data from pilot instrumentation.
📋 Case Connection
📋 Thermal Management System for EV Traction Inverter
Peak junction temps >175°C causing derating and reliability concerns
📋 Thermal Design of Satellite Payload Radiator for Lunar Orbit Mission
Extreme radiative environment: solar flux up to 1360 W/m², albedo up to 0.12, IR emission from hot regolith (~390 K)