一、一维热传导方程
MATLAB实现（显式法）：
% 参数设置
L = 1.0; % 杆长 (m)
Nx = 50; % 空间节点数
dx = L/(Nx-1); % 空间步长 (m)
alpha = 0.01; % 热扩散率 (m²/s)
T_left = 100; % 左端温度 (℃)
T_right = 0; % 右端温度 (℃)
dt = 0.001; % 时间步长 (s)
Nt = 1000; % 时间步数% 初始化温度分布
T = zeros(Nx, Nt);
T(:,1) = 20; % 初始温度 (℃)
T(1,:) = T_left; % 左边界条件
T(end,:) = T_right;% 右边界条件% 显式差分迭代
r = alpha*dt/dx^2;
for n = 1:Nt-1for i = 2:Nx-1T(i,n+1) = T(i,n) + r*(T(i+1,n) - 2*T(i,n) + T(i-1,n));end
end% 可视化
x = linspace(0, L, Nx);
surf(linspace(0, dt*Nt, Nt), x, T');
xlabel('时间 (s)'); ylabel('位置 (m)'); zlabel('温度 (℃)');
title('一维热传导显式法模拟');
关键点：
-
显式法稳定性条件：r≤0.5
-
可视化使用
surf函数展示温度场随时间演化。
二、二维热传导方程
MATLAB实现（交替方向隐式法，ADI）：
% 参数设置
Lx = 1.0; Ly = 1.0; % 板长 (m)
Nx = 50; Ny = 50; % 网格数
dx = Lx/(Nx-1); dy = Ly/(Ny-1);
alpha = 0.01; % 热扩散率
h = 10; % 对流系数
T_inf = 20; % 环境温度% 初始化温度场
T = 20*ones(Nx,Ny);% 边界条件（对流边界）
for n = 1:1000T(1,:) = T_inf + (T(2,:) - T_inf)*exp(-h*dx/k);T(end,:) = T_inf + (T(end-1,:) - T_inf)*exp(-h*dx/k);T(:,1) = T_inf + (T(:,2) - T_inf)*exp(-h*dy/k);T(:,end) = T_inf + (T(:,end-1) - T_inf)*exp(-h*dy/k);
end% ADI迭代
r_x = alpha*dt/(2*dx^2); r_y = alpha*dt/(2*dy^2);
for iter = 1:1000% x方向隐式for j = 2:Ny-1A = diag(1 + 2*r_x*ones(Nx-2,1)) + diag(-r_x*ones(Nx-3,1), 1) + diag(-r_x*ones(Nx-3,1), -1);b = T(2:end-1,j) + r_x*(T(3:end,j) - 2*T(2:end-1,j) + T(1:end-2,j));T(2:end-1,j) = A\b;end% y方向隐式for i = 2:Nx-1A = diag(1 + 2*r_y*ones(Ny-2,1)) + diag(-r_y*ones(Ny-3,1), 1) + diag(-r_y*ones(Ny-3,1), -1);b = T(i,2:end-1) + r_y*(T(i,3:end) - 2*T(i,2:end-1) + T(i,1:end-2));T(i,2:end-1) = A\b;end
end% 可视化
[X,Y] = meshgrid(linspace(0,Lx,Nx), linspace(0,Ly,Ny));
contourf(X,Y,T',50); colorbar; title('二维热传导ADI法模拟');
关键点：
-
ADI法结合显式与隐式格式，无条件稳定
-
对流边界采用指数格式处理。
三、三维热传导方程
MATLAB实现（有限元法，PDE工具箱）：
% 使用PDE工具箱求解三维热传导
model = createpde('thermal','transient');
geometryFromEdges(model,@(p) createCubeGeometry(1)); % 创建1×1×1立方体
generateMesh(model,'Hmax',0.1);% 材料属性
thermalProperties(model,'ThermalConductivity',1, 'Density',1, 'SpecificHeat',1);% 边界条件
thermalBC(model,'Face',1:6, 'Temperature', 20); % 全固定温度% 初始条件
thermalIC(model, 0);% 时间求解
tlist = linspace(0, 1, 100);
result = solvepde(model, tlist);% 可视化
pdeplot3D(model, 'ColorMapData', result.Temperature(:,end));
title('三维热传导PDE工具箱模拟');
关键点：
-
PDE工具箱支持复杂几何和自动网格划分
-
适用于工程复杂结构（如电子元件散热）。
参考代码 一维，二维，三维热传导方程的数值求解 www.youwenfan.com/contentcnq/63240.html
四、通用优化策略
-
稳定性控制：
-
显式法需满足 r≤0.5（一维）、r≤0.25（二维）
-
隐式法（如ADI）无条件稳定，适合长时模拟。
-
-
网格自适应：
% 在PDE工具箱中设置自适应网格 generateMesh(model, 'GeometricOrder', 'linear', 'Hmax', 0.05); -
并行计算：
% 使用parfor加速三维迭代 parfor iter = 1:1000% 并行计算各节点温度 end
五、应用案例对比
| 场景 | 维度 | 方法 | MATLAB工具 | 结果精度 |
|---|---|---|---|---|
| 金属杆瞬态加热 | 一维 | 显式/隐式法 | 自定义代码 | 高 |
| 电子芯片散热 | 二维 | ADI法 | 自定义代码 | 中高 |
| 地热储层温度场 | 三维 | PDE工具箱 | PDE工具箱 | 高 |
