圆周石墨密封流体温度场MATLAB分析（微分求积法求解二维能量方程）
问题描述
我们需要使用微分求积法(DQM)分析圆周石墨密封装置中流体温度场的二维分布，求解能量控制方程。
数学模型
二维稳态能量方程（不考虑内热源）：
λ(∂²T/∂r² + (1/r)∂T/∂r + ∂²T/∂z²) = ρc_p(v_r ∂T/∂r + v_z ∂T/∂z)
其中：
- T(r,z)为温度场
- λ为导热系数
- ρ为密度
- c_p为比热容
- v_r, v_z为流速分量
MATLAB实现代码
clear; clc; close all;
