做碲锌镉（CdZnTe，CZT）这类化合物半导体晶体生长模拟的人，早晚会面对同一个问题：凝固过程根本不是一个纯传热问题，也不只是一个流场问题，更不是单纯画一条固液界面移动线就能交代的几何变化——它是变形几何、固体传热和流场三个过程搅在一起的整体行为。市面上能摸到的晶体生长仿真方案很多，但真正肯用COMSOL Multiphysics把这三样东西同时放进一个模型里跑的人，都会经历一段"怎么搭都别扭、怎么调都发散"的时期。这篇文章不谈教科书推导，只讲我在CZT结晶过程模拟里实际做过的事情：模型怎么搭、三场耦合怎么设、网格怎么救，以及那些让COMSOL报错但帮助文档里根本不会写清楚的细节。给正在做晶体生长仿真、特别是准备用COMSOL动网格的朋友一份可以直接抄作业的参考。
1. 为什么碲锌镉结晶模拟必须同时盯着传热、流动和几何变形
1.1 碲锌镉晶体的实际应用与生长工艺背景
碲锌镉（CdZnTe，常简写CZT）是碲化镉基体中固溶一定比例锌形成的三元化合物半导体，典型组分为Cd₀.₉Zn₀.₁Te。它之所以在半导体材料里地位特殊，是因为平均原子序数高、密度接近6 g/cm³，对X射线和伽马射线有很强的阻止能力，是目前少有的能在室温下工作的核辐射探测器材料。同时它还是HgCdTe红外探测器的主流衬底，探测器级CZT晶体的市场价值一直很高，但生长难度也常年排在前列。
生长CZT晶体的主流工艺是垂直布里奇曼法（VBM）和移动加热器法（THM）。VBM的操作逻辑很直观：把多晶原料密封在高纯石英安瓿里，整只安瓿以1~3 mm/h的速度在具有特定温度分布的炉膛中移动，熔体从冷端开始定向凝固。这个速度慢到什么程度？一根150 mm长的晶锭往往要连续生长两三天。工艺窗口又窄，炉温梯度、安瓿移动速度、熔体对流强度只要有一点偏差，轻则界面形状不对，重则长出多晶或出现组分偏析。靠纯实验去试错，成本和周期都扛不住，所以数值模拟成了工艺前移的主要手段。
而这类模拟最麻烦之处在于：凝固界面本身在移动，移动的界面改变固体域和液体域的比例，进而改变热阻分布，热阻分布又反过来决定界面位置的推进速度。整个问题天然是强耦合的，想偷懒只算其中一个场，结果基本没有参考价值。
1.2 凝固过程中三个物理场的耦合逻辑
先把物理链条理清楚。热量从高温炉区经安瓿壁传入熔体，在固液界面处熔体凝固并释放潜热，这部分热量通过已凝固的晶体向冷端导出。与此同时，熔体内部因为轴向和径向都存在温度差，密度分布不均匀，浮力驱动出自然对流；如果熔体有自由表面，表面张力随温度变化还会叠加一层热毛细对流，也就是Marangoni对流。
流场一旦起来，熔体里的热量输运就不再只是导热，对流项会把温度场搅得面目全非。温度场一变，固液界面的位置和形状也跟着变。界面所在位置决定了高导热固体域和低导热液体域各占多大空间，这个几何比例又反过来改变整个计算域的传热路径。更麻烦的是，界面移动本身有自己的时间尺度，整个生长过程看似几天里缓慢推进，但实际上每一步都接近稳态、又永远不完全稳态，这种"准静态推进"对数值方法非常不友好。
这三个场的耦合强度并不对称。传热和界面几何之间是强耦合：界面位置由热平衡决定，热阻分布又由界面位置决定；流场与传热之间在格拉晓夫数Gr比较大时也极强；流场与几何之间的耦合相对弱一些，但在动网格框架下，熔体域收缩时流场边界形态改变，对涡结构和界面稳定性影响很直接。理解了这个主次关系，后面设计求解策略才有依据。
1.3 为什么选COMSOL Multiphysics而非自编程序
早期做晶体生长模拟的人大多自编有限体积或有限差分程序。一套二维轴对称代码要从N-S方程写起，再叠加热方程、固液界面追踪、潜热释放，没有大半年很难跑出像样的结果，而且每改一种工艺条件就要改一次边界条件，维护成本极高。COMSOL Multiphysics的优势在于把这些物理场封装成了标准接口：固体传热、层流、变形几何/移动网格都可以分别添加，再用多物理场节点把它们绑定在一起，初版模型一周内就能跑通。
对晶体生长模拟来说，COMSOL最有价值的其实是"变形几何"（Deformed Geometry）和"移动网格"（Moving Mesh）这一族工具。前者处理欧拉框架下计算域形状随时间变化的问题，允许材料流过变形网格；后者处理ALE框架下边界移动追踪问题，网格跟着材料边界走。CZT凝固过程里这两种思路都有用武之地，后面我会专门说怎么选。
当然COMSOL也不是没有脾气。全耦合求解器面对强非线性时经常不收敛，动网格畸变会毫无征兆地报错，而且很多报错信息根本不说人话。接下来这几章，重点就是怎么避开和解决这些问题。
2. 模型搭建：从材料参数到几何与边界条件
2.1 几何与坐标系选择
一上来就做全尺寸3D是新手最容易踩的坑。CZT安瓿直径几十毫米、高度两百毫米，界面附近又需要微米级分辨率才能分辨温度梯度和流动结构，3D模型随便一划就是几十万甚至上百万单元，再加上非稳态动网格，普通工作站根本转不动。我强烈建议从2D轴对称开始：垂直布里奇曼的安瓿绕轴旋转，几何和外场条件基本是轴对称的，用2D轴对称建模就能抓住主部物理，网格量省一个数量级。
几何上我习惯建三个域：熔体域（初始为液相）、晶体域（初始为固相）、石英安瓿壁。以直径D=80 mm、总长200 mm的安瓿为例，熔体初始高度取100 mm，晶体初始高度取100 mm。在COMSOL里画成矩形分区即可。注意2D轴对称模型的宽度是半径而不是直径，旋转对称轴放在左边。安瓿壁要单独画成一个域，后面给材料属性和外壁换热条件都方便。
这里有一个关键决定：初始的"熔体/晶体分界面"到底怎么处理。如果后续用等效热容法，两个域其实使用同一套温度相关材料属性，界面位置完全由温度场决定，不是几何上的真实边界。如果走显式移动网格追踪路线，则两个域在几何建模上就是不同的域，界面作为真实边界存在，固液属性分别指定。我的经验是：初版模型用等效热容法快速摸规律，确认流场和温度场都稳定后，再升到显式移动网格做精细化复核。
2.2 CdZnTe材料热物性参数的选取与处理
材料参数是CZT模拟最刁钻的一关。不同文献给出的热物性数据很分散，尤其是液态热导率、动力粘度和表面张力温度系数这三个，差别能达到两三倍。我通常先按主流文献值建一套基线参数，再对最不确定的参数做灵敏度扫描，绝不把一个孤立的数值当场真理。
| 参数 | 固相（晶体） | 液相（熔体） | 简要说明 |
|---|---|---|---|
| 密度（kg/m³） | 约5900 | 约5680 | 液相比固相约低4%，浮力项要配合参考密度使用 |
| 比热容（J/(kg·K)） | 约160 | 约200 | 文献分散，建议做灵敏度 |
| 热导率（W/(m·K)） | 约1.0~1.5 | 约1.0~2.0 | 液态有效热导率受对流影响，要区分表观值和本征值 |
| 动力粘度（Pa·s） | — | 约1.5×10⁻³~3×10⁻³ | 温度相关，常用指数律近似 |
| 热膨胀系数（1/K） | 约5×10⁻⁶ | 约2×10⁻⁴ | 体积膨胀系数，浮力项关键参数 |
| 液相线温度/熔点（℃） | 约1092 | 约1092 | Zn含量升高则液相线略升 |
| 熔化焓（kJ/kg） | 约200 | 约200 | 等效热容法必须用 |
| 表面张力温度系数（N/(m·K)） | — | 约-2×10⁻⁴ | 自由液面模型才需要 |
别小看固液密度那4%的差别。在布里奇曼生长中，熔体受热膨胀而晶体致密，密度差既是自然对流的驱动力，也是后期热应力的来源。做流动模拟时我在非等温流多物理场节点里勾选Boussinesq近似，参考密度取ρ₀=5680 kg/m³，浮力项写作ρ₀·β·(T-T₀)·g，β取2×10⁻⁴ K⁻¹。这个量级下，半径方向几十K的温差就足以驱动出显著的自然对流，流场不能忽视。
比热和热导率的温度依赖如果有实测数据，强烈建议做成插值函数而不是常数。CZT固态热导率随温度变化很明显，从高温降到室温能差好几倍。不过如果只关注凝固阶段的界面附近温度场，也就是1050~1150℃这个窄窗口，常数假设勉强能接受，前提是初版模型。
石英安瓿的参数也要认真给：热导率约2 W/(m·K)，比热约700 J/(kg·K)，密度约2200 kg/m³。很多人把注意力全放在CZT上，却忘了安瓿壁本身的热阻对界面温度场的影响。厚度几个毫米的石英壁，在热流路径上造成的影响经常比想象中大得多。
2.3 边界条件与初始条件设置
炉温分布是整个模型最重要的边界条件。实测的Bridgman炉轴向温度分布是一条典型S型曲线：高温区约1150℃，低温区约1000℃，中间绝热区温度梯度最大。在COMSOL里，我习惯用"定义→解析函数→分段"定义轴向温度曲线T_outer(z)，再把它作为安瓿外壁的温度边界条件。如果工艺是炉体移动而不是安瓿下降，这个温度曲线在空间上要随时间平移，写成T_outer(z - v_growth·t)即可。
安瓿外壁不建议一上来就用表面对表面辐射。辐射项是T⁴非线性，会直接把初次运行的稳态和瞬态求解器拖入泥潭。先用第一类给定温度或带对流系数的第三类边界把模型跑通，之后再开辐射做精细化验证，成功率会高很多。
熔体接触安瓿壁的内壁设置为无滑移边界。如果保留自由上表面，上表面设置为滑移边界，同时施加表面张力温度系数产生的Marangoni剪切应力。这个应力不能忽略，CZT熔体粘度虽然不算低，但自由液面附近的Marangoni流动对局部温度场影响很明显。
初始条件设置同样有讲究。我通常先单独求一次纯传热稳态，也就是只开导热不开流场和动网格，得到温度初值T_0；流场初值设0，压力初值设0。这比直接让一个非均匀温度场在耦合求解器里硬算要稳定得多。把整个模型初始温度都设为炉温曲线对应值，也比均匀设成1120℃再让边界慢慢拉平衡要靠谱。
2.4 潜热处理的注意点
潜热是能量方程里最凶的源项。最常见的做法是等效热容法，把相变潜热折算进一个等效比热容。核心公式就是：
cp_eq = cp + Lf * d(fL)/dT fL = 0.5 * (1 - tanh((T - Tm) / dT))其中fL是液相分数，Tm是熔点，dT是人为设定的相变过渡半宽度。用光滑的tanh函数代替阶跃，是为了避免两个网格节点之间发生温度跳变，否则求解器会把时间步压缩到无法接受。
dT的取值直接影响数值表现。ΔT取2~4 K是一个比较稳的起步区间。ΔT越小越接近真实相变面，但数值越硬；ΔT太大又会人为扩大两相区厚度，导致算出来的界面位置模糊。CZT这种需要精细获取界面形状的场景，我一般从ΔT=3 K开始，数值稳定后再逐步收紧。
等效热容法还有一个副作用：相变区间内等效比热会出现一个高尖峰，让能量方程变成强非线性，跟流场耦合时尤其敏感。所以正确的操作顺序是：先固定网格、固定流场，单独把温度场和相变收敛跑稳；再开流场；最后才引入动网格。分阶段逼近，比一次性全耦合省太多时间。
3. 三场耦合的实现细节：COMSOL接口设置与求解策略
3.1 物理场接口搭建
具体操作上，在"添加物理场"里选择以下接口：
- 传热→固体传热（Heat Transfer in Solids，ht），作用域选熔体域、晶体域和石英壁；
- 流体流动→单相流→层流（Laminar Flow，spf），作用域只选熔体域；
- 数学→变形网格→变形几何（Deformed Geometry，dg）或移动网格（Moving Mesh，ale），作用域选所有发生变形的域。
如果采用"等效热容+网格变形"的组合，物理本质上还是材料域在变形，固体传热和流场都定义在变形后的坐标里。这里有一个很容易漏掉的细节：用移动网格时，层流接口的物理模型记得勾选与网格速度相关的选项。COMSOL在移动网格下默认会用ALE形式的N-S方程，自动带出网格速度修正项，但如果你是从固定网格模型复制修改过来的，一定要去"物理模型"面板确认一遍，否则移动边界会在动量方程里产生虚假源项，结果怎么调都怪。
多物理场节点按这个思路搭：一是"非等温流动"（Nonisothermal Flow），把层流与流体传热耦合，它会在流场里加入浮力项、在能量方程里加入对流项；二是流体-固体边界传热，让熔体与晶体在界面上满足温度连续和热流连续。用等效热容法后，这个条件基本由材料属性隐式完成，但显式建一个节点方便后续查看界面热流。
材料节点里，固液两套属性分别指定到两个域。在等效热容法框架下，更常规的做法是只在"界面可能扫过的区域"启用含潜热尖峰的等效比热，晶体远端和熔体远端保持单一相，这样能避免不必要的非线性。
3.2 界面移动的驱动方式：Stefan条件与变形几何
这部分是核心中的核心。凝固界面移动的本构关系是Stefan条件，物理含义一句话：界面法向单位面积上释放的潜热，等于固相侧导走的热流与液相侧补进来的热流之差。公式写出来是：
ρ_s · Lf · v_n = k_s · ∂T_s/∂n - k_l · ∂T_l/∂n
界面法向速度v_n不是人为指定的，而是由界面两侧的温度梯度差决定。这是整个变形几何模型里最重要的"生长法律"。
等效热容法的好处就在于：不需要显式写出这个界面速度表达式，因为能量方程已经把潜热消耗吞进去了，温度场会自动让液相分数等值线移动。后处理时用fL=0.5等值线提取界面位置，再对时间求导，就能得到生长速度。这是集中参数式的界面追踪，实现难度低，适合先把整条物理链路打通。
如果必须显式追踪界面，则在移动网格接口里把界面边界上的法向速度写成上面那个公式的形式。这个速度不是常数，而是温度梯度的函数，需要在边界变量或边界偏微分方程里承载。COMSOL支持用边界偏微分方程来解决这类问题，但调试难度明显上了一个台阶，符号写反、法向方向取错都会直接发散成垃圾结果。所以我不建议一上来就硬刚显式Stefan条件，先拿等效热容法把物理趋势摸对，再用移动网格做复核更省钱。
这两种方法如果拿到的界面形状差在3%以内，说明模型不仅稳定而且正确。差很多的话，大概率不是数值方法的问题，而是材料参数或边界条件在两个模型里没有保持一致。
3.3 求解顺序与时间步控制
求解器配置是COMSOL动网格仿真里最大的坑之一。默认的"全耦合"求解器在三个强非线性物理场同时启动时，十次有九次不收敛。我的操作顺序是：
第一步，关闭流场和变形几何，单独解稳态传热，拿到初始温度场。 第二步，打开层流，固定网格，做瞬态非等温流，让流动达到准稳态。 第三步，打开变形几何或移动网格，保留上一步结果作为初值。 第四步，改用分离式（Segregated）求解器，把压力-速度放一组，温度场单独一组，网格位移单独一组。每个子问题都是相对成熟的小框架，整体反而比强耦合更快更稳。
时间步进用BDF（向后差分公式），同时设置最大时间步。凝固速率按1 μm/s估算，如果界面附近最小网格尺寸是0.1 mm，界面穿过一个单元需要100秒。时间步要是超过几十秒，界面状态可能一下跳过好几个单元，网格变形会突然剧烈起来。我习惯给BDF设"初始步长不超过0.1 s，最大步长不超过30~60 s"，让它在这个框架内做自适应。
容差建议收紧到1e-4，而不是默认的1e-3。默认容差对瞬态多物理场往往偏松，某些时刻已经偏离物理规律但求解器还在闷头跑。容差收紧后界面位置的锯齿明显减少，虽然会慢几倍，但翻车率低一个数量级。另外一定要设停止条件保护，比如界面移动超过设定值或出现非物理温度就提前终止，避免烧掉一整天算出一堆垃圾。
4. 网格畸变、数值振荡与收敛性踩坑记录
4.1 变形几何中网格畸变的典型表现与成因
动网格模型报错最吓人的一句话是"网格反移"或"倒置单元"，有时候直接报negative Jacobian。看到这些别急着改网格，先回去看几何为什么变成那样。CZT布里奇曼生长里网格畸变的主要导火索有三个：界面推进导致熔体域轴向变薄；温度梯度大、界面处网格被剧烈拉伸；流体连续性方程与动网格速度耦合产生压力尖峰。
界面向熔体方向推进几十毫米后，熔体域最初100 mm的高度最后可能只剩几个毫米，网格单元一边被压缩一边被拉长，长宽比飙升。COMSOL变形几何内置了Laplace和Winslow等网格平滑算法，但能力有限，网格相对厚度低于约3~5层时，界面附近的压力离散就会开始病态。
我的习惯是在整个模型里给变形几何设置一个"初始网格位移为零"的初始条件，然后在所有外边界用"滑移位移"而不是"固定位移"，让网格能沿安瓿壁滑移来容纳变形。很多人一出问题就把边界全设成固定，结果界面动一下，整片网格像拉一块硬板似的全被拽弯，倒置单元必然出现。
4.2 一次"不收敛"问题的完整定位过程
分享一次真实的排查过程，可以当作模板用。那轮参数扫描里，基线工况能跑通，但把冷端冷却速率从0.5 K/min改成1.5 K/min后，t≈3000 s时求解中止，报错信息是"最后步长小于最小步长，瞬态求解失败"。
我的第一步不是调网格，而是关掉流场，只保留传热+变形几何。结果这个组合顺利跑过了3000 s，说明问题不在动网格本身。第二步关掉变形几何，保留传热+流场，也能跑，说明问题不在流热耦合。问题就这样被锁定在"流场+动网格"的联合效应上。
第三步，打开流场和动网格，把时间步固定在一个极小值跑几步，看中间输出。果然，在熔体域厚度压缩到大约还剩8 mm时，界面附近压力出现高频振荡，紧接着最低层的网格单元开始反转。原因清楚了：熔体域变薄、界面又持续释放潜热，流场在越来越窄的通道里被"挤"，压力-速度耦合在动网格框架下没有足够阻尼，振荡超过阈值后网格直接塌掉。
解法分三路：第一，在熔体域顶部与气相接触的边界设置一个压力约束，给封闭空间一个泄压通道；第二，在界面边界上对网格位移施加额外的切向滑移自由度，让网格能顺着壁面滑而不被压扁；第三，检查层流接口的稳定化设置，COMSOL默认的流线扩散没有被自定义设置偷关掉就保留。改完后同一个扫描点顺利跑通，后面更极端的参数点只是变慢，没有再崩。
4.3 稳定收敛的几招储备
整理几条我在CZT模拟里沉淀下来的保命经验。
设置网格重剖分触发条件。COMSOL变形几何可以在网格质量低于阈值时自动重新划分网格，这个功能很诱人，但要慎用，任何重剖分都会带来历史场变量的插值误差，Crank-Nicolson类格式对插值尤其敏感。实在要用，就设一个"最大单元扭曲度"阈值，让重剖分尽量少触发。
给界面法向速度加限制。把界面推进速度限制在一个物理合理范围，比如0~5 μm/s，防止初期微小数值振荡被放大成不合理的界面位移。这个限制在COMSOL变量定义里用min和max函数就能实现，成本极低但效果明显。
用参数扫描接力逼近目标工况。把目标工艺参数拆成多个子工况，从简单到复杂逐个计算，每一步以上一步的解作为初值。这比直接怼一个极端工况稳定得多，本质上是用前一个解当了热启动。
定期检查能量守恒。每跑完一段，用COMSOL的报告功能统计"界面释放潜热总量"和"通过边界累计导热量"，两者偏差超过5%就优先怀疑潜热过渡区太宽或网格分辨率不足。能量账不平，后面所有结果都不能信。
5. 结果后处理与物理可信度检验
5.1 界面形态演化与生长速率判据
等效热容法的直接输出是温度场和液相分数场。界面位置可以通过派生值来定义：在每条轴向采样线上找fL=0.5的z坐标。把界面各径向位置的z坐标与中心位置的z坐标做差，就是界面凸度。经验值上，界面中心略微凸向熔体1~3 mm是比较理想的工艺状态，这个形状有利于把杂质推向晶体外周而不是滞留在中心，对探测器级CZT尤其重要。
生长速率后处理也很直接：记录中心点界面z坐标随时间的变化，取斜率就是局部生长速度。理想情况下，经过一段起步过渡后，生长速度应接近安瓿移动速度，比如1 mm/h。如果差太多，说明热平衡没有建立，或者冷却速率设置太快，界面在追赶温度场。这种"追赶"在实验中表现为界面过冷或局部回熔，模拟的价值就在这里——不用动炉子就能提前摸清工艺参数区间。
把界面推进速率、熔体最大流速和界面凸度三者放在同一张时间曲线上观察，几乎一眼就能看出某组参数是否"和谐"。如果流速随时间剧烈波动、界面推进速度跟着抖，那多半是对流进入了非稳态区间，需要进一步检查格拉晓夫数。
5.2 流场与温度场的交叉验证
解出流场后，先别急着截图，先做无量纲判断。格拉晓夫数Gr = g·β·ΔT·H³/ν²，以熔体高度H=50 mm、温差ΔT=30 K、体积膨胀系数2×10⁻⁴、运动粘度5.3×10⁻⁷ m²/s估算，Gr大概在10⁶量级，这已经是非常确定的自然对流，不是可以忽略的次级效应。如果你模拟出来的流速小到几乎看不见，第一时间怀疑Boussinesq项的符号或参考密度设置错了。
温度场和流场的交叉检查我从三个维度做。一是界面附近的等温线是否被对流拉变形，纯导热情况下等温线应该平滑地穿过界面，对流则会带来明显扭曲。二是对称轴附近是否出现规则环流，环流的方向应该与浮力方向一致：热区上升、冷区下降。三是通过界面的总热流是否等于"潜热消耗+两侧导热的收支差"，热流账不平基本就是潜热处理的问题。
我还喜欢看一个细节：熔体顶部或自由液面附近是否存在第二个反向涡。垂直布里奇曼里只出现一个主环流是常见情况，但如果工艺参数比较极端，可能出现上下层双涡结构，这种结构对溶质分布影响很大。做参数扫描时可以专门输出涡心位置随参数的变化，这比单看界面形貌更早暴露潜在不稳定。
5.3 参数扫描与工艺优化视角
等模型跑通，COMSOL的参数化扫描就变成了效率工具。用辅助扫描把炉温梯度、安瓿下降速度、熔体初始高度作为参数列表，一次性提交多个工况。每个工况都把上一个工况的解作为初值，相当于给每个极端参数点搭了一座桥，收敛率大幅提升。
扫描结果的判读要回到物理。提高炉温梯度会让界面更凸向熔体，但伴随而来的是更大的热应力；降低下降速度会让界面更接近稳态、流场更平静，但生长周期被拉长。用派生值表把每个参数点的界面凸度、最大流速、轴向温度梯度、界面过冷度整理成表，基本就能画出工艺窗口的轮廓。
有一点必须提醒：模拟得到的趋势比具体数值更可靠。比如"梯度增加10 K/cm，界面凸度增大约0.5 mm"这样的趋势，对工艺调整有直接指导意义；但"梯度30 K/cm时界面凸度是1.52 mm"这种绝对值只能当参考，因为材料参数和边界条件总有不精确的地方。我在对外输出报告时，永远把结果写成一个区间而不是一个点，这是做仿真的人最基本的专业素养。
6. 继续做下去的几个方向与我的实际体会
6.1 从二维到三维、从热流到场
2D轴对称模型能回答的问题已经很多，但它天然过滤掉了周向非对称扰动。真实安瓿在炉膛里不是完美对中，炉体也在缓慢旋转，这些因素会在晶体里留下螺旋状杂质条纹。想研究这些，就必须上3D模型并考虑安瓿旋转，无论是旋转坐标系还是滑移网格，复杂度都是2D的几倍，建议在2D物理彻底摸清后再碰。
更值得优先做的是把第三个场加进来——溶质场。CdZnTe熔体中的Zn和Cd在凝固时会发生再分配，造成组分偏析。用稀物质传递接口耦合浓度场，配合界面处的分凝系数，可以定量预测沿晶锭方向的电阻率均匀性。这个结果对探测器材料来说，比纯温度分布更接近用户真正关心的指标。
6.2 我的体会
回头做CZT结晶模拟最有价值的时刻，不是模型第一次跑通，而是拿模拟得到的温度曲线和实测热电偶数据对比、发现两者相差在2℃以内那次。从那以后我确定了一件事：材料参数校准和边界条件确认的优先级，高于任何高级物理算法。很多人一上来就追求相场、追求3D、追求辐射换热，结果连基线传热都没校出来，这是最划不来的。
我的最后一条建议是：先用最朴素的方法把能量守恒做对，再慢慢往模型里加细节。变形几何、固体传热、流场三者的耦合，本质上是把晶体生长最核心的物理搬进电脑里。基础骨架稳了，后面加什么细节都只是锦上添花。
