Finite Difference Heat Conduction Simulator

The core physics is governed by the 1D heat conduction equation (Fourier's Law):

Here, T(x,t) is temperature, t is time, x is position, and α is the thermal diffusivity (a material property). This equation says the rate of temperature change at a point is proportional to the spatial curvature (second derivative) of the temperature profile.

The Explicit Finite Difference Method discretizes this equation in space and time. The key update formula for each interior grid point i at time step n+1 is:

Where r = αΔt/Δx² is the stability parameter. The term \((T_{i+1}^n - 2T_i^n + T_{i-1}^n)\) is the finite difference approximation of the second spatial derivative ∂²T/∂x². For numerical stability, we must enforce r ≤ 0.5. This condition links the physical property (α), the numerical grid (Δx), and the time step (Δt).

Building Insulation Design: Engineers use this 1D analysis to model heat transfer through walls. By simulating different materials (changing α) and thicknesses (L), they can predict heat loss over time and optimize for energy efficiency, crucial for designing green buildings.

Electronic Cooling: Predicting temperature rise in a silicon chip or a heat sink fin is often modeled as 1D conduction. The explicit method allows for rapid simulation of transient "hot spots" after a processor starts a heavy computation, informing heatsink design.

Materials Processing: In quenching or annealing metals, controlling the cooling rate is vital to achieve desired material properties. Simulating the temperature history through a steel plate's thickness helps design the process to avoid cracks or residual stresses.

Geothermal Analysis: Modeling how temperature varies with depth in the Earth's crust over daily or seasonal cycles is a classic 1D heat conduction problem. This principle is used to design ground-source heat pump systems.

First, do not confuse "thermal diffusivity α" with "thermal conductivity k". While you directly manipulate thermal diffusivity α in the simulator, material data sheets in practice often list thermal conductivity k. Their relationship is $\alpha = k / (\rho c_p)$. For example, copper has high thermal conductivity (about 400 W/mK), so α is also large, meaning heat spreads quickly. Conversely, polystyrene foam has low thermal conductivity, so α is small, and heat tends to stay localized. When setting parameters, you must consider density ρ and specific heat $c_p$ together as a set.

Next, note that "grid size Δx and time step Δt cannot be chosen independently". The stability condition $r \le 0.5$ can be satisfied by either decreasing Δt or increasing Δx. For instance, if you refine Δx from 1mm to 0.5mm, the spatial resolution improves, but you must then reduce Δt to a quarter of its previous value; otherwise, r becomes too large. Consequently, the number of calculation steps needed to simulate the same real-time duration explodes, increasing computation time. In practice, you should always be mindful of the trade-off between "accuracy" and "computational cost".

Finally, be aware that the heat transfer coefficient h in the Robin boundary condition varies greatly depending on the situation. While you can change h to values like 10 or 100 in this tool, in reality, the order of magnitude differs significantly between natural convection (approx. 5–25 W/m²K in air), forced convection (approx. 25–250 W/m²K with a fan), or water cooling (approx. 500–10,000+ W/m²K). If you set this value too far from reality, your valuable simulation results will fail to reflect the real world. Therefore, it is essential to verify this carefully using literature or experimental data.