C++实现一维热传导方程的有限差分数值解法

1. 项目概述：一维热方程的有限差分数值求解

热传导方程（Heat Equation）是描述热量在介质中传递过程的经典偏微分方程，在工程热物理、材料科学、地质勘探等领域有广泛应用。这个项目展示如何用C++实现一维热传导方程的有限差分数值解法，适用于随时间变化的热传导问题模拟。

有限差分法（Finite Difference Method）的核心思想是用离散的网格点代替连续空间，用差分近似代替微分运算。对于一维热方程 ∂u/∂t = α·∂²u/∂x²，我们将分别用前向差分处理时间导数，用中心差分处理空间二阶导数，构建显式迭代格式。这种方法计算效率高，代码实现直观，特别适合作为数值计算的入门案例。

注意：显式格式虽然简单，但存在稳定性限制。时间步长Δt必须满足CFL条件：α·Δt/Δx² ≤ 0.5，否则计算结果会发散。

2. 数学建模与算法设计

2.1 热方程离散化处理

考虑一维热传导方程：

code复制∂u/∂t = α·∂²u/∂x²,  0 ≤ x ≤ L, t ≥ 0

其中α是热扩散系数，u(x,t)表示温度分布。

采用均匀网格离散：

• 空间步长Δx = L/N，网格点xi = i·Δx, i=0,1,...,N
• 时间步长Δt，时间层tn = n·Δt, n=0,1,2,...

用差分近似微分：

• 时间导数：前向差分 ∂u/∂t ≈ (u_i^{n+1} - u_i^n)/Δt
• 空间二阶导数：中心差分 ∂²u/∂x² ≈ (u_{i+1}^n - 2u_i^n + u_{i-1}^n)/Δx²

得到显式迭代公式：

code复制u_i^{n+1} = u_i^n + r·(u_{i+1}^n - 2u_i^n + u_{i-1}^n)

其中r = α·Δt/Δx² 是稳定性参数。

2.2 边界条件处理

常见边界条件类型：

1. 狄利克雷边界（固定温度）：

   cpp复制u[0] = left_temp;  // 左端固定温度
   u[N] = right_temp; // 右端固定温度
   

2. 诺伊曼边界（绝热/热流）：

   cpp复制// 绝热边界示例（热流为0）
   u[0] = u[1]; 
   u[N] = u[N-1];