2020年那道“带相变材料的低温防护服御寒仿真模拟”A题从数学建模到程序求解我把整个过程完整跑了一遍。说实话这道题表面上是个传热问题真正做起来才发现它同时在考材料科学直觉、数值稳定性判断和工程简化能力。如果你正在备战数学建模竞赛尤其是对偏物理/工科的A题没底这篇复盘应该能给你一条可以抄作业的完整链路从控制方程怎么列、相变潜热怎么处理到有限差分程序怎么写、参数怎么标定再到结果怎么分析才能拿分。1. 为什么是相变材料题目背后的工程痛点1.1 从普通防护服到“智能控温”需求2020年的特殊背景相信大家还有印象医护人员长时间穿着防护服工作闷热和不透气问题非常突出。但低温场景下问题刚好反过来——冷链运输、北方冬季户外检疫、寒区作业等场合普通防护服保温能力有限穿厚了影响操作灵活性穿薄了又扛不住低温。传统做法是堆砌隔热棉层但隔热层只能“减缓”热量散失温度还是会持续下降时间一长人就扛不住了。这里就引出了相变材料Phase Change MaterialPCM的核心价值。PCM在一定温度下会从固态熔化成液态熔化过程会吸收大量潜热而温度几乎保持不变。把它缝进防护服里就相当于内置了一个“恒温缓冲池”当外界低温试图把人体热量带走时PCM先熔化吸热让服装内表面温度在相当长一段时间内维持在相变点附近。这个逻辑用大白话讲就是普通保温杯只能减慢水变凉的速度而PCM相当于在杯壁里藏了一块冰水温降到冰点时会被冰的熔化“撑住”很久。竞赛题把这样一个真实工程需求抽象成数学问题本质上是让你回答三个层次的问题第一PCM层到底能维持多久的防护时间第二不同环境温度和材料参数下温度场怎么演化第三怎么优化PCM的厚度、熔点等参数来满足设计指标。1.2 三层结构模型的物理图景低温防护服在实际产品里可能有四五层但建模阶段没必要全堆上去。我当时的做法是简化为三层结构最外层是防水耐磨面料中间是PCM功能层最内层是贴近人体的舒适/隔热层然后皮肤表面单独处理。这样既保留了对热行为影响最大的结构特征又避免因为层数过多导致参数冗余、数值计算量失控。每一层的热物性差异非常大外层导热系数不算低但很薄热阻贡献很小PCM层导热系数低、比热容大而且存在相变潜热内层通常是比较蓬松的纤维材料密度小、导热系数低类似羽绒服的保暖原理。这三层串联起来相当于三个热阻和两个热容PCM还有个潜热热容构成的动态热系统。把题目翻译成工程语言之后你就会发现它不是纯粹让你做热力学计算而是考察“如何把多层结构、相变过程、人体热调节、环境对流换热”这些物理因素统一写进一组方程组里。这也是建模题最关键的环节——物理图像不清晰后面程序写得再漂亮都是空中楼阁。1.3 竞赛题的正确打开方式先定指标再建模我见过很多队伍拿到题就开始查文献、套公式这是最要命的。建模的第一步应该是定义“输出指标”这道题用什么量来评价防护服的性能通常可选的指标包括皮肤表面温度降到某个阈值比如10°C或15°C的时长、PCM完全熔化所需的时间、最大热流密度、温度均匀性等。我当时定的核心防护指标是“PCM层完全熔化前皮肤表面温度不低于18°C的持续时间”。为什么用18°C而不是更低因为皮肤长时间接触低于这个温度的表面会明显感觉冰冷甚至产生冻伤风险。这个指标同时串联了材料参数潜热、熔点、厚度和环境条件环境温度、对流强度后续灵敏度分析和优化设计都围绕它展开。指标先行还有一个好处它帮你排除了大量不必要的物理细节。比如你不需要精确模拟人体血管舒张收缩只需要给皮肤表面一个合理的边界条件不需要考虑相变材料在熔化后的自然对流因为服装夹层内空间封闭对流强度远弱于导热。建模不是越复杂越好而是把所有因素按“对目标指标的影响程度”排队只保留影响大的。2. 传热模型的推导从傅里叶定律到“带内热源的三层微分方程”2.1 每一层的控制方程到底怎么写核心工具就是傅里叶导热定律加能量守恒即固体/静止流体中的导热微分方程[ \rho c_p \frac{\partial T}{\partial t} \frac{\partial}{\partial x}\left(k \frac{\partial T}{\partial x}\right) \dot{q} ]其中 (\rho) 是密度(c_p) 是比热容(k) 是导热系数(\dot{q}) 是内热源项。这里我把问题简化为一维模型方向为从皮肤向外表面因为防护服的典型厚度只有一到两厘米而横向尺寸是米量级边缘效应可以忽略。人体侧的内热源来自代谢产热。更精确的做法是使用生物传热领域的Pennes方程[ \rho_t c_t \frac{\partial T}{\partial t} k_t \frac{\partial^2 T}{\partial x^2} \rho_b c_b \omega_b (T_{art} - T) \dot{q}_{met} ]这里引入了血液灌注项(\rho_b c_b \omega_b (T_{art} - T))表示动脉血以一定灌注率向组织传热。但在防护服问题里人体组织我们未必需要建出来通常只需要在皮肤表面给定边界条件即可。这个取舍我在后面会细说。对于PCM层问题就复杂了在相变过程中材料会吸收潜热等效比热容不再是常数而是会在相变温度附近出现一个很大的“尖峰”。如果你直接用固定比热容去计算温度会在相变温度附近产生不合理的平台或突变数值上很难收敛。2.2 相变潜热怎么“骗”进方程焓法与显热容法的抉择处理相变潜热有两大流派焓法和显热容法又叫等效比热容法。焓法的思路是用焓 (H) 代替温度 (T) 作为基本未知量把潜热自然地包含在焓的变化里。控制方程变成[ \frac{\partial H}{\partial t} \frac{\partial}{\partial x}\left(k \frac{\partial T}{\partial x}\right) ]然后再通过温度和焓的转换关系包括显热区间和相变区间把温度恢复出来。焓法的最大优点是守恒性强不容易出现潜热“丢失”的问题缺点是实现时要处理分段函数。显热容法则是把相变潜热折算成相变温度区间内的一个附加比热容[ c_{eff} c_s \frac{L}{T_l - T_s} \quad (T_s \le T \le T_l) ]其中 (L) 是潜热(T_s) 和 (T_l) 分别是固相线和液相线温度。这样方程形式上和其他层一致但注意这里有个隐含假设相变发生在一个温度区间而不是某个单一温度点上。纯物质的相变温度是一个点实际材料比如石蜡混合物通常有一个几摄氏度的熔化区间所以这个方法在工程上完全可行。我最终选择的是焓法原因有两个第一显热容法在相变区间附近温度和焓关系曲线比较陡非线性迭代容易振荡而焓法更稳定第二焓法能很方便地追踪“液相分数”——我可以随时知道PCM层已经熔化百分之多少这对结果分析非常有价值。2.3 边界条件与初始条件最容易丢分的地方很多队伍把微分方程列对了结果画出来的曲线一看就有问题十有八九是边界条件写错了。低温防护服问题有三个重要的边界/初始条件外表面边界条件最外层直接接触低温空气存在对流换热写成 [ -k_{out} \left. \frac{\partial T}{\partial x} \right|{xL} h (T{surf} - T_{\infty}) ] 这里 (h) 是对流换热系数(T_{\infty}) 是环境温度。皮肤表面边界条件皮肤侧可以给第一类边界条件固定温度如 (T_{skin}33°C)也可以给第二类/第三类边界给定热流或与血液温度耦合。我选择第二类给一个稳定的生理热流密度理由是这样能反映人体实际产热能力而不是假设皮肤温度恒定——现实中皮肤温度会随着热量流失而下降。初始条件假设穿着者从室温环境25°C左右进入低温环境初始时刻整个防护服温度均匀。这里有个细节值得提醒如果你直接让初始时刻服装温度为25°C但环境温度是-20°C那么从 (t0^) 开始外表面就承受巨大的温差驱动力第一秒的热流可能高达几千瓦每平方米这在物理上是不符合实际的。真实的场景是人体先穿着防护服在室温里待了一段时间服装内部已经接近体温水平进入低温环境是“突然”发生的所以初始条件应该是让内部温度基本均匀且接近体温。可以做一个短暂的“预热平衡”模拟来生成更合理的初始温度分布而不是简单赋一个常数。2.4 为什么选择一维模型精度与效率的权衡有人可能会问为什么不直接建三维模型用现成的有限元软件如COMSOL或ANSYS跑仿真竞赛场景首先要考虑的是可控性和可解释性。你在论文里写“我用COMSOL做了仿真”评委看不到代码细节而且商业软件的黑箱很难展示你真正理解了物理过程。手写一维有限差分程序每层网格几十个节点时间步长一秒以内几分钟就能跑完一轮参数扫描效率远比三维模型高。更重要的是一维模型已经足以回答这道题的核心问题——沿厚度方向的温度分布和时间演化。你想知道的“PCM层撑多久”本身就是一维问题。三维模型引入的额外自由度比如边缘效应、衣服褶皱、空气间隙不均匀对数值结果的影响远小于参数不确定性带来的误差。建模原则永远是“够用就好”在竞赛时间有限的情况下更是如此。3. 数值求解我为什么选有限差分而不是软件仿真3.1 空间离散与非均匀网格的取法确定了控制方程之后数值离散是下一个关键决策。我用的是经典有限差分法空间上采用非均匀网格。为什么不用均匀网格因为三层材料厚度差异很大外层可能只有0.5mmPCM层可能达到5-8mm内层2-4mm。如果整体用均匀网格为了分辨薄层至少要几千个节点大部分算力浪费在不需要精细分辨的地方。我的做法是每一层单独划分网格层内网格均匀层与层之间通过界面通量连续性衔接。比如外层用8个节点PCM层用30个节点内层用10个节点总共不到50个节点精度已经非常可观。网格独立性的验证我放在后面说但先剧透一个经验节点数翻倍目标指标防护时间变化控制在1%以内就可以认为网格收敛了。很多人在这上面堆节点数完全没必要。3.2 显式时间推进的稳定性约束与步长选择显式格式前向欧拉实现简单但有稳定性条件。一维导热问题显式格式的稳定性要求是[ \Delta t \le \frac{\Delta x^2}{2\alpha} ]其中 (\alpha \frac{k}{\rho c_p}) 是热扩散系数。细网格和低密度的内层材料会强烈限制步长。比如内层材料导热系数0.05、密度80、比热1500算出来 (\alpha \approx 4.2 \times 10^{-7}) (\text{m}^2/\text{s})如果网格尺寸 (10^{-4}) m允许的最大步长约为0.012秒。这意味着模拟两小时需要60万步虽然每步计算量不大但参数扫描时会比较耗时。这里有两种应对思路一是缩小内层网格尺寸把它并入PCM层用较大网格处理二是直接用隐式格式。隐式格式无条件稳定步长可以取到0.5秒甚至更大但每步需要解三对角方程组。我最终选择了隐式格式用Python的scipy.linalg.solve_banded解三对角矩阵每个时间步的计算开销比显式大但总步数从60万降到1万左右整体提速十分可观。而且隐式格式对非线性源项相变区间的焓-温度关系的迭代收敛也更宽容。3.3 焓法迭代的细节实现焓法的核心是焓场与温度场的相互转换。每一步的流程是用当前焓场 (H^n) 和上一轮温度场 (T^n)计算新的焓场 (H^{n1})通过离散导热方程。根据 (H^{n1}) 反算温度 (T^{n1})。如果节点处于相变区间需要由液相分数插值得到温度。如果节点处于相变区间且温度变化剧烈需要迭代1-2次保证收敛。焓与温度的关系可以写成[ H(T) \begin{cases} \rho c_s (T - T_{ref}) T T_s \ \rho c_s (T_s - T_{ref}) \rho L \cdot \frac{T - T_s}{T_l - T_s} T_s \le T \le T_l \ \rho c_s (T_s - T_{ref}) \rho L \rho c_l (T - T_l) T T_l \end{cases} ]实际操作中我把相变区间的焓-温度关系拆成了两步先算液相分数 (f_l \frac{H - H_s}{H_l - H_s})限制在 ([0,1])再用 (T (1-f_l) T_s f_l T_l) 恢复温度。这样数值实现非常干净而且在输出液相分数用于结果分析时也是现成的。3.4 程序模块结构核心代码骨架我使用Python完成的整个求解框架大致分为四个模块参数定义、网格生成、时间推进、后处理。核心时间推进部分的骨架如下import numpy as np from scipy.linalg import solve_banded # 参数定义 # L1, L2, L3: 各层厚度 # k1, k2, k3: 导热系数 # rho2, L_latent, Ts, Tl: PCM密度、潜热、固液相线 # h_conv, T_inf: 对流系数、环境温度 # 网格生成每层单独划分界面处共用节点 N1, N2, N3 8, 30, 10 x np.concatenate([ np.linspace(0, L1, N1, endpointFalse), np.linspace(L1, L1L2, N2, endpointFalse), np.linspace(L1L2, L1L2L3, N31) ]) # 注意端点去重形成按厚度方向排列的节点数组 # 构建三对角系数矩阵 A隐式格式 # A 的对角元素取决于各节点的 k, rho, cp_eff或 dh/dT # 时间推进 for n in range(int(t_total / dt)): # 1. 组装系数矩阵和右端项含边界条件 # 2. 解三对角线性方程组得到 H^{n1} H_new solve_banded((1, 1), ab, rhs) # 3. 由 H_new 反算 T_new 和液相分数 f_l T_new enthalpy_to_temperature(H_new) f_l liquid_fraction(H_new) # 4. 更新边界记录皮肤表面温度、PCM液相分数等 T T_new这段代码虽然简化了细节但整体结构就是竞赛时的实际程序没有用任何现成有限元包。这么做的好处是每一步你都知道计算机在算什么调试和论文解释都很有底气。程序里最需要小心的两个地方界面处的热流连续性处理以及焓反算温度时的边界分段判断。4. 材料物性标定与关键参数所有数字都要有出处4.1 PCM选型石蜡和其他候选材料当时查阅文献和材料手册可供选择的相变材料主要集中在三类石蜡类正十八烷、正二十烷、无机水合盐如六水氯化钙、脂肪酸类如癸酸、月桂酸。它们的物性差异很大材料类型熔点范围潜热 (kJ/kg)导热系数 (W/mK)优点缺点石蜡正十八烷27-29°C220-2450.15-0.25潜热大、化学稳定、无过冷导热低、可燃六水氯化钙29-30°C170-1900.5-1.0导热好、价格低过冷明显、易腐蚀癸酸30-32°C150-1800.15-0.25可再生、无毒价格较高、有气味从防护服场景出发最合适的是相变点在28-32°C的石蜡类材料因为这个温度刚好比人体皮肤表面温度33°C左右低几度又远高于环境低温PCM熔化时能“吸走”足够多的热量而不至于让皮肤感觉过冷。潜热取200 kJ/kg是一个合理标杆值实际竞赛中你也可以对多个候选材料做对比分析这本身就是加分项。4.2 对流换热系数的取值逻辑外表面和环境之间的对流换热系数 (h) 不是随便拍脑袋定的。无风环境下自然对流的典型值在5-10 W/(m²K)如果有风强制对流可以到20-50 W/(m²K)甚至更高。我从传热学的经验关联式出发按空气掠过大平板、风速1-3 m/s的场景估算取了一个基准值 (h15) W/(m²K)然后把它作为参数扫描变量之一范围取5-30。这里有一个很多人忽略的点对流换热系数对防护时间的影响是高度非线性的。h从5增加到10防护时间可能缩短40%但从20增加到30缩短幅度就没那么大了。原因很简单当h增大到一定程度外表面温度已经非常接近环境温度整体传热热阻由服装本身主导环境侧的对流热阻占比很小。这个“热阻串联”的思路在结果分析时非常有用它解释了为什么在强风环境下单纯加厚PCM层对提升防护时间的边际收益在递减。4.3 人体产热模型与皮肤边界人体代谢产热在不同活动强度下差别很大静坐时约60-80 W轻度活动100-150 W中重度活动可达300 W以上。如果按体表面积1.8 m²折算单位面积产热大约在40-180 W/m²。但要注意这些热量并不是全部通过皮肤表面传导出去的相当一部分通过呼吸蒸发、出汗蒸发散失。因此在防护服建模中有效通过服装传导的热流密度通常在20-50 W/m²。我是这样设置皮肤边界条件的不直接给定皮肤温度恒定而是在皮肤表面给定一个恒定的有效热流密度 (q_{skin} 30) W/m²作为基准工况反映人体静坐状态下的净散热需求。这样皮肤温度会自然波动如果服装保温好皮肤温度会稳定在舒适区间如果环境太冷皮肤温度会逐渐下降直到逼近危险阈值。这种处理方式比“皮肤恒温33°C”更符合真实生理而且不会出现“不管外部多冷皮肤永远33°C”这种明显违背常识的结果。材料参数方面我还把PCM和纤维材料的导热系数随温度变化考虑为分段常数固态和液态分别取值并在PCM层增加了一个很小的等效辐射导热修正。事实上对于夹层结构如果你只考虑导热熔融液态PCM内部的自然对流和辐射换热会低估热传递。工程上常用“等效导热系数”方法在液相区把 (k) 放大1.5-3倍。我做了一组敏感性测试发现这个修正对防护时间的影响大约在5-10%不算主导因素但引入它可以体现你对物理过程的理解深度。5. 仿真结果怎么分析温度场、液相分数与防护时间5.1 典型温度时间曲线的三个特征阶段程序跑通之后第一件事不是画一堆漂亮的云图而是盯着“皮肤表面温度-时间”曲线和“PCM层平均温度-时间”曲线看。这两条曲线能直观告诉你整个热过程经历了几个阶段第一阶段快速降温期。防护服刚进入低温环境时服装整体温度梯度迅速建立皮肤表面温度以一个较快的速率下降。这个过程通常持续几分钟取决于服装初始温度和环境温差。内层材料越薄、导热越快这个阶段降温斜率越陡。第二阶段相变平台期。当PCM层的温度降到熔点附近材料开始熔化大量潜热被吸收相当于在热量传递路径上“耗掉”了一大部分能量。这个阶段皮肤表面温度下降速率明显放缓曲线出现一个明显的“膝弯”之后进入近似线性缓慢下降的阶段。平台期的长度直接取决于PCM的潜热和厚度以及通过服装的总热流。第三阶段失效期。PCM完全熔化后潜热缓冲耗尽材料变成单一的液态导热体温度再次加速下降。如果在熔化完全之前皮肤温度还没有跌破阈值那么从这个时刻起防护服的性能就迅速恶化。这也说明一个设计原则PCM层的设计目标必须是“在预定防护时长内不要完全熔化”否则后段会非常被动。这三个阶段的划分对论文写作特别重要它们直接对应了PCM防护服“延迟降温”的物理机制而且你可以很方便地画出液相分数随时间的曲线来佐证。5.2 液相分数推进速度揭示的设计瓶颈另一个关键输出是PCM层液相分数的空间分布和时间演化。我发现一个非常有意思的现象相变并不是从PCM层的某一侧均匀推进的。由于内外两侧温差不同PCM层靠近外层的那部分先熔化然后相变界面逐步向内层推进。当PCM层较厚时内层甚至还能保持固态很久。这意味着“PCM总熔化量”不能简单地用“平均温度达到熔点”来估算。这个现象带来的设计启示是如果你只关心防护时间那么不一定需要让整个PCM层都熔化。设计时可以定义一个“有效PCM利用率”在防护时间结束时实际熔化的PCM体积占总PCM体积的比例。我做的仿真显示PCM层厚5mm时两小时内利用率大约只有40-60%而厚度10mm时利用率进一步下降。换句话说盲目增加PCM厚度前50%的厚度在提供显热隔热后面50%则只是“备用缓冲”并没有被充分利用。这个结论对后续优化非常关键与其无限加厚PCM层不如优化内层隔热层的热阻让PCM在更长时间内保持在相变点附近。这也是为什么后来我做优化时把内层厚度和PCM厚度作为双变量联合优化而不是单调地增加PCM。5.3 防护时间的定义方式不能只看皮肤温度防护时间的定义有几种口径选择不同结论可能完全相反。我比赛时用的是“皮肤表面温度不超过18°C的持续时间”。但后来复盘发现更稳健的做法是同时给出两个指标一是“皮肤温度低于18°C的时间点”二是“PCM层液相分数达到90%的时间点”。两者交集才是真正的“安全窗口”。为什么这么说因为只看皮肤温度可能会被初始温度场误导。如果初始温度设得比较高比如30°C皮肤需要一段时间才能降到18°C这段时间即使PCM已经熔化了一半皮肤表面温度可能还很高。只看温度会高估防护时间看液相分数能揭示PCM寿命的真实剩余。两个指标一起看才能对防护服状态有一个客观的评估。我最终的评判标准定为防护时间 皮肤温度降到18°C对应的时刻且要求此时PCM液相分数不超过80%。这样既保证了热舒适性又留出了安全裕度也更符合工程设计的理念。5.4 灵敏度分析与参数扫描结论我做参数扫描时的变量包括环境温度、对流换热系数、PCM厚度、PCM潜热、内层厚度每组变量扫描区间设为基准值的±50%。结论可以用一句话概括对防护时间影响最大的三个因素依次是PCM潜热、PCM厚度、环境温度对流换热系数的影响排第四内层厚度的影响最小。这个排序不是拍脑袋出来的而是通过归一化灵敏度系数算出来的[ S \frac{\Delta t_{protection} / \bar{t}_{protection}}{\Delta p / \bar{p}} ]其中 (p) 是某个参数。我发现PCM潜热增加20%防护时间大约提升25-30%PCM厚度增加20%防护时间提升约10-15%而内层厚度增加20%防护时间只提升3-5%。这说明在基准工况下内层隔热已经接近“瓶颈”——它的热阻虽大但受到本身厚度的限制可调空间有限。真正的设计杠杆在PCM层。这个灵敏度分析结论值得写进论文的“优化设计”部分因为它不仅给出了最优参数方向还告诉了读者“为什么是这个方向”。评委看到这种级别的分析比看到你堆了十个三维温度分布云图要有用得多。6. 复盘验证方法、常见坑与延伸思考6.1 如何用极限情形验证程序程序写完之后必须做验证否则直接交付的结果很容易被评委挑刺。我的验证策略有三个层次第一层退化验证。把PCM的潜热设成0相变温度范围缩到无限小模型应当退化为普通三层导热问题。这时我可以手动计算稳态温度分布和程序结果比较。如果差异超过1%说明相变处理模块有问题。第二层能量守恒检查。计算整个模拟时间内进入系统的总能量和离开系统的总能量之差应该等于系统内焓的总变化。我用这个指标追踪每一个时间步的累计误差结果表明隐式格式的能量守恒误差在0.1%以内非常理想。第三层网格独立性验证。连续加密网格观察防护时间的变化。当网格尺寸从2mm缩到0.5mm时防护时间变化仍在1%以内说明结果已经收敛。我把这三层验证结果全部写进了论文附录这既是学术严谨性的体现也在答辩时省了很多口舌。6.2 竞赛现场常见的四个翻车细节第一初始条件处理不当。直接给整件衣服设一个统一温度导致第一秒的热流爆表。解决方法是先跑一个持续的常温平衡阶段再切换边界条件模拟进入低温环境。第二相变区间太窄导致迭代振荡。如果熔点在28-30°C之间而 (T_s) 和 (T_l) 只差0.1°C焓-温度关系曲线几乎垂直任何数值求解器都会颤抖。解决办法是适当放宽相变区间到2-4°C或者把PCM建模为混合物多组分依次熔化。第三外表面边界条件用了固定温度。如果你设置外表面温度恒等于环境温度相当于假设衣服外表面没有对流热阻这会导致热量散失被严重高估防护时间可能被低估20-30%。正确做法是对流边界条件。第四时间步长选择拍脑袋。隐式格式虽然无条件稳定但步长过大时时间离散误差会拖垮精度。我的建议是先做一组时间步长收敛测试找到“步长再缩小结果变化不超过0.5%”的临界值用这个值跑正式模拟。6.3 从竞赛题到实际产品的延伸思考坦白讲竞赛模拟和真实产品之间还有一条不小的鸿沟。真实防护服中PCM往往被封装修复成微胶囊掺杂在纤维内部而不是一整块材料。这种情况下PCM的分布是离散的微胶囊壳层的导热热阻也会影响吸放热速率这些在竞赛模型里都没有考虑。不过从另一个角度看竞赛模拟恰恰给了你一个做“系统级设计”的机会。你可以不考虑微胶囊怎么造但需要知道PCM潜热多大、放多少克、放在哪一层最优。这个从材料微观设计到系统热管理的视角正是实际工程中热设计工程师每天都在做的事情。如果你对这方面的兴趣被勾起来了可以找个时间用同样的框架把微胶囊PCM的等效热物性也用公式估算出来那是另一道很有意思的题目。最后分享一个个人经验拿到题不要急着写代码先花半天时间把“模型里到底包含什么过程、忽略什么过程”想清楚。竞赛A题和B题最大的区别在于A题一旦物理图像错了后面程序写得再花哨也拿不到分。把傅里叶定律、对流换热、相变潜热这三块搞透再加上有限差分的稳定性概念这道题你就已经掌握了对付它的全部武器。至于代码本身反而是整个过程中最不稀缺的部分。