4. 圆柱坐标系下轴对称热传导建模
问题设定与几何简化
本例考察一个半径 R = 0.05 m、高度 H = 0.20 m 的实心圆柱体。圆柱外侧面维持恒温 T_w = 400 K，上下端面绝热，内部无内热源，材料导热系数 k = 45 W/(m·K)（典型碳钢量级）。由于几何、边界条件与材料均关于中心轴旋转对称，且端面绝热，温度只沿半径方向变化，即 T = T(r)，与 θ、z 无关。这正是采用二维轴对称接口的物理依据：用 (r, z) 平面上的一个矩形截面代表整个三维圆柱，计算量大幅下降而结果不失真。
轴对称成立的前提是：几何回转对称、材料各向同性、边界条件不随周向角变化。若圆柱侧面存在局部热源或周向不均匀冷却，则必须改用完整三维模型。
控制方程与解析解
在圆柱坐标系下，稳态、无内热源、仅径向变化的导热方程为：
d/dr ( r · dT/dr ) = 0积分两次得通解 T(r) = C₁·ln r + C₂。代入边界条件：
- 外壁
r = R：T = T_w = 400 K - 轴心
r = 0：由对称性dT/dr = 0，且温度有限
由于 ln r 在 r = 0 处发散，物理上要求 C₁ = 0，于是 T(r) = T_w 为常数。这说明：在无内热源、端面绝热、外壁恒温的实心圆柱中，稳态温度场是均匀的，径向不存在温度梯度。这个结论看似平凡，却是检验 COMSOL 模型正确性的关键基准——如果仿真出现明显径向温差，说明边界条件或对称性设置有误。
若圆柱为空心（内半径
r_i > 0），则 C₁ ≠ 0，温度呈对数分布 T(r) = T_w + (q·r_i/k)·ln(r/R)，径向温差显著。本例特意选用实心圆柱，正是为了获得一个"零梯度"的解析基准，便于快速判断模型是否搭建正确。COMSOL 建模关键设置
在 COMSOL 中选择 二维轴对称 空间维度，物理场选择 固体传热（Heat Transfer in Solids），研究类型为 稳态（Stationary）。几何只需绘制 r ∈ [0, 0.05]、z ∈ [0, 0.20] 的矩形。关键设置如下：
| 项目 | 设置值 | 说明 |
|---|---|---|
| 材料导热系数 k | 45 W/(m·K) | 各向同性，常数 |
| 外侧面边界（r = R） | 温度 T = 400 K | 恒温边界 |
| 上端面（z = H） | 热绝缘 | 默认边界，无需额外设置 |
| 下端面（z = 0） | 热绝缘 | 默认边界 |
| 对称轴（r = 0） | 轴对称 | 软件自动识别，无需手动添加 |
| 初始值 | T₀ = 400 K | 加速收敛 |
网格与求解
由于温度场为常数，网格密度对结果几乎无影响，但仍建议使用物理场控制网格（常规单元大小）以保证求解稳定。求解后，在结果中绘制 r 方向的温度分布曲线：沿 z = 0.10 m 处从 r = 0 到 r = 0.05 m 取一条截线，观察温度随半径的变化。
结果判读与验证
仿真完成后，沿径向截线提取温度数据。预期结果：
- 从
r = 0到r = 0.05 m，温度应稳定在400 K附近，波动幅度小于0.1 K（由数值离散误差引起）。 - 若出现明显径向温差（例如轴心温度低于外壁超过 1 K），应检查：外壁恒温边界是否施加在正确的边上；对称轴是否被误设为热绝缘以外的条件；材料导热系数是否被误设为随温度变化。
- 热流密度
q_r = -k·dT/dr在整个域内应接近零，可在结果中绘制r方向热流密度云图辅助确认。
本例的核心价值不在于获得复杂温度场，而在于建立"轴对称简化 → 解析基准 → 仿真验证"的闭环工作习惯。后续涉及空心圆柱、多层壁面或内热源的问题，都可以在此基础上扩展边界条件与材料属性，而几何与接口选择保持不变。