---
title: "读 Karma 2001 的合金相场论文，并复现其中的枝晶图"
description: "这篇四页短文提出了反截留流。本章带你读：结论在哪里，怎样检验收敛，以及我们复现了哪些图中的数值。"
canonical_url: "https://sandboxphysics.com/zh/stories/crystal-defects-microstructure/reading-karma-2001/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/crystal-defects-microstructure/reading-karma-2001/"
last_modified: "2026-10-30"
---

[晶体缺陷与显微组织](https://sandboxphysics.com/zh/stories/crystal-defects-microstructure/)

L4 · 第 7 章

# 读 Karma 2001 的合金相场论文，并复现其中的枝晶图

Sandbox Physics · 2026-10-30

Karma 在 2001 年的这篇短文中提出：加一股溶质流，稀合金相场模型就能用足以算出枝晶的宽界面，同时界面处仍满足局部平衡。证据是两张图里的收敛检验。本章带你一边开着本实验，一边读这篇论文：每个结论在哪里，这两张图能说明什么、不能说明什么，复现时哪些要自己选，以及我们的运行与哪些数值吻合。

## 开始前需要什么

第 6 章：尖锐界面问题、宽界面造成的溶质截留、反截留流和参数关系。能读懂画出的曲线和简短的渐近论证。

- 在论文中找到模型、参数选择、收敛检验和溶质检验，并说明每张图支持什么。
- 在较粗的比值下复现图 1 和图 2 中的平台尖端速度和溶质剖面，并说明这个比较检验了什么、没检验什么。

## 一段话说清结论

早先的合金相场模型，为了避免虚假的截留、表面扩散和界面拉伸，界面不能比毛细长度厚多少，因此模拟慢速生长代价太高。Karma 加入反截留流，并选好插值函数，让这三个效应在薄界面极限下全部消失，模型于是回到带局部平衡的经典自由边界问题。支持这个结论的有两样：一段薄界面分析的梗概，以及二维枝晶模拟：界面宽度改变时，模拟的尖端速度和溶质剖面保持不变。

来源：A. Karma

## 先看哪里

短文前半部分，式 (1)—(12)，讨论液相线与固相线平行的理想合金，说明这股流怎样消除化学势跳变，并给出常数 $a_1$ 和 $a_2$。稀合金从式 (13) 开始：式 (13)—(15) 是要逼近的尖锐界面问题，式 (16)—(18) 是模型，式 (20) 是插值函数，式 (22) 是四重各向异性。

结尾附近有一段，列出了复现所需的全部数字：$\Delta x=0.4$、$\Delta t=0.008$、$W=\tau=1$、$\epsilon_4=0.02$、$k=0.15$、$\Omega=0.55$，以及半径 $22\,d_0$ 的晶核。图 1 画出两种模型在 $d_0/W=0.277$ 和 $0.544$ 下尖端速度 $Vd_0/D$ 随 $tD/d_0^2$ 的变化。图 2 画出沿枝晶轴留在固相中的溶质。薄界面推导本身，文中说将另文发表；Echebarria 等人在 2004 年发表了它。

来源：A. Karma · B. Echebarria

## 论文怎样检验收敛

枝晶没有精确解可以对照。Karma 改为让每个模型在两种界面宽度下各跑一次：如果模型已经到达薄界面极限，两条曲线就必须重合。Karma 的新模型（下文称反截留模型）通过了检验，两个比值下的速度和溶质剖面都重合。标准模型没有通过，溶质剖面最明显，即使在较细的比值下也远未收敛。

另一个独立的检验是自洽性。用运行中测得的尖端半径 $\rho$，局部平衡可以预测尖端后方固相应含多少溶质。Karma 报告，反截留模型与预测相差在百分之几以内。再看代价的论证意味着什么：$d_0/W$ 减半，一次运行便宜约 $32$ 倍，所以收敛检验同时也在检验网格能粗到什么程度。

$$
V\,\frac{d_0}{D}\Big|_{d_0/W=0.277}\;\overset{?}{=}\;V\,\frac{d_0}{D}\Big|_{d_0/W=0.544},\qquad \frac{c_s(x)}{c_l^0}\;\overset{?}{=}\;k\Bigl[1-(1-k)\frac{d_0}{\rho}\Bigr]
$$

图 1 和图 2 背后的两项检验：两种界面宽度之间相互吻合，以及固相浓度与按测得尖端半径算出的局部平衡值吻合。

$$
\text{cost}\propto\Bigl(\frac{d_0}{W}\Bigr)^{5}:\qquad \frac{\text{cost}(0.544)}{\text{cost}(0.277)}\approx2^5=32
$$

Karma 对一次运行的代价如何随毛细长度与界面宽度之比增长的估计。

来源：A. Karma

## 哪些要自己选

短文没有给出方盒大小、边界条件、晶核位置，也没说尖端速度和半径怎样测量。我们的运行只算晶体的四分之一：晶核放在角上，边界取镜像，方盒边长 $128\,W$，在较粗的比值下相当于 $462\,d_0$。尖端取在轴线上 $\phi$ 过零处；半径由轴线附近四行上零等值线拟合的抛物线给出；速度在 $5000$ 到 $6000\,d_0^2/D$ 之间取平均。

图 1 和图 2 只以曲线给出数值，所以下面每个参照数都是我们目测读出的，速度约精确到 $0.001$，浓度约精确到 $0.005$。溶质平台在曲线平坦的 $150$ 到 $400\,d_0$ 之间读取。

来源：A. Karma

## 完整例题：在较粗的比值下复现图 1 和图 2

题目。在 Karma 的网格上，以 $d_0/W=0.277$ 运行两种模型，把平台速度和固相浓度与他的图比较。

反截留模型。尖端速度为 $0.0180$，图 1 约为 $0.018$。固相平台为 $0.1447\,c_l^0$，图 2 约为 $0.145$。尖端半径 $21.3\,d_0$，局部平衡据此预测 $0.144$，运行的自洽程度在 $0.5\%$ 以内。

标准模型。尖端以 $0.0244$ 前进，图上约为 $0.025$。固相含 $0.2247\,c_l^0$，图上约为 $0.225$；峰值 $0.3185$ 出现在 $50\,d_0$，图上约 $0.315$，也在 $50\,d_0$ 附近。尖端半径为 $16.9\,d_0$，局部平衡预测 $0.142$，所以这里的固相多含了 $58\%$。

较细的比值，较粗的网格。$d_0/W=0.544$ 我们只在 $\Delta x=0.8$ 的网格上运行，网格间距是 Karma 的两倍。标准模型的固相平台为 $0.181$，图上约为 $0.183$；峰值 $0.262$ 在 $30\,d_0$，图上约 $0.265$，在 $27\,d_0$ 附近。反截留模型的速度为 $0.0166$，比图上低 $8\%$；在这个网格上速度没有收敛。

结论。在 Karma 的网格上、$d_0/W=0.277$ 时，两种模型都在读数误差内复现了图中的值；反截留模型与局部平衡自洽，标准模型则在固相中留下过多溶质。

来源：A. Karma

## 这次复现没有说明什么

这些都是确定性的单次运行，没有种子间的离散可报；基准文件中的离散只是平均时间窗内的起伏。参照值是目测读数，不是论文给出的数字。较细的比值没有在 Karma 的网格上重跑，他的收敛检验本身因此没有复现。复现的是较粗的那一半，外加较细那一半在粗网格上的一次考察。

为了能交互运行，本实验在 $\Delta x=0.8$ 的网格上跑较粗的比值。那里尖端速度为 $0.0161$，比 Karma 网格上慢 $11\%$，而固相仍含 $0.1446\,c_l^0$。溶质分配经得起粗网格，尖端选择经不起。这里没有与任何实验比较，也没有重新推导薄界面渐近分析，后者 Echebarria 等人给出了完整推导。运行是二维的、等温的，不含噪声。

来源：A. Karma · B. Echebarria

## 到实验里试一试

### [把本实验和图 1 对照](https://sandboxphysics.com/zh/labs/alloy-dendrite-antitrapping/)

1. 选“四臂枝晶”，让它长。先按图 1，再按本实验较粗网格的路线基准，预测运行后期尖端速度的读数。
2. 把界面模型切换为“标准模型，无反截留流”。预测尖端速度和溶质读数相对 Karma 的标准模型曲线会往哪里走。
3. 解释为什么收敛的模型在两种界面宽度下必须给出相同的曲线，以及为什么本实验的粗网格能让尖端变慢，却几乎不改变溶质读数。

## 练习与解答

1. 在以 $W$ 和 $\tau$ 为单位、网格间距和时间步长固定的条件下，解释 Karma 关于运行代价按 $(d_0/W)^5$ 增长的估计。

   展开解答

   方盒以 $d_0$ 为单位固定，所以边长方向的格点数按 $d_0/W$ 增长，二维格点数按 $(d_0/W)^2$ 增长。以 $\tau$ 为单位的时间单位 $d_0^2/D$ 等于 $(d_0/W)^2/(a_2a_1W/d_0)$，按 $(d_0/W)^3$ 增长，步数也一样。两者相乘，按 $(d_0/W)^5$ 增长。

2. 对同样的方盒和时间，本实验的网格 $\Delta x=0.8$、$\Delta t=0.032$ 比 Karma 的便宜多少？精度上付出了什么？

   展开解答

   格点少四倍，步数也少四倍，大约便宜 $16$ 倍。尖端速度降低 $11\%$，为 $0.0161$；固相浓度仍为 $0.145$。

3. 两条曲线在两种界面宽度下重合。这证明了模型的什么，又没证明什么？

   展开解答

   它说明结果已经不依赖界面宽度，这正是薄界面极限所承诺的。它不说明结果与实验相符，也不说明结果没有与宽度无关的误差，比如以毛细长度计的网格间距、方盒大小，或者各向异性的离散方式。用独立预测做的局部平衡检验因此才重要。

[上一章：合金相场模型中的溶质截留与反截留流](https://sandboxphysics.com/zh/stories/crystal-defects-microstructure/thin-interface-alloy-solidification-and-antitrapping/)

## 来源与复算范围

1. [A. Karma, Phase-field formulation for quantitative modeling of alloy solidification, arXiv:cond-mat/0103289v1 (2001)](https://arxiv.org/abs/cond-mat/0103289v1)

   第 1 版，全文。式 (16)—(18) 与 (22)：含反截留流的等温稀合金模型。式 (19)—(20)：平衡剖面，以及消除表面扩散的 $q(\phi)$ 取法。参数段落：$\Delta x=0.4$、$\Delta t=0.008$、$k=0.15$、$\Omega=0.55$、$\epsilon_4=0.02$、$a_1=0.8839$、$a_2=0.6267$，晶核半径 $22\,d_0$。图 1：尖端速度随时间的变化。图 2：沿枝晶轴固相中的溶质。

2. [B. Echebarria, R. Folch, A. Karma and M. Plapp, Quantitative phase-field model of alloy solidification, arXiv:cond-mat/0404164v2 (2004)](https://arxiv.org/abs/cond-mat/0404164v2)

   第 2 版。第 I B 节，式 (2)：界面厚度取介观尺度时，溶质截留、表面扩散和界面拉伸都与 $W$ 成正比地放大。第 I D 节：把反截留流作为非变分的修正。第 I E 节，式 (7)—(13)：完整模型，式 (10) 的反截留流，系数 $a=1/(2\sqrt2)$，以及 $d_0=a_1W/\lambda$ 和式 (13) 的 $\beta$。第 II A 节，式 (17)—(21)：等温生长的尖锐界面问题，包括 Stefan 条件和 Gibbs–Thomson 关系。第 IV 节给出薄界面分析；本课程只略读，没有重新推导。