相场模拟——合金，金属凝固模型，各向异性枝晶生长karma 合金凝固模型，选区激光熔融，激光增材制造，选择性激光熔融，SLM，定向凝固，熔铸 1matlab，实现合金各向异性枝晶生长，代码注释详细，包学包会，直接运行得到下面的演化过程，提供相关详细文献资料，包括控制方程，求解方法 2comsol偏微分方程，雪花生长模型，纯金属枝晶生长模型，详细资料

karmaAlloy：一个基于相场法的二元合金凝固过程模拟器

概述

相场模拟——合金，金属凝固模型，各向异性枝晶生长karma 合金凝固模型，选区激光熔融，激光增材制造，选择性激光熔融，SLM，定向凝固，熔铸 1matlab，实现合金各向异性枝晶生长，代码注释详细，包学包会，直接运行得到下面的演化过程，提供相关详细文献资料，包括控制方程，求解方法 2comsol偏微分方程，雪花生长模型，纯金属枝晶生长模型，详细资料

karmaAlloy 是一款面向科研与工程场景的二维二元合金凝固模拟器。它采用相场-溶质场耦合框架，在无需显式追踪固液界面的前提下，即可再现枝晶形貌演化、微观偏析及溶质分布等关键现象。整套代码以 MATLAB 语言编写，单文件即可运行，兼顾可读性与计算效率，适合桌面端快速实验，也方便向高性能平台移植。

核心能力

1. 各向异性枝晶生长
   通过引入四重对称的各向异性函数，自动诱发一次、二次乃至三次枝晶臂，真实再现“海藻状”形貌。
2. 溶质场-温度场耦合
   在相场方程中嵌入温度梯度与溶质再分配效应，可观察“温度梯度区熔”与“组分过冷”带来的形貌转变。
3. 守恒型溶质扩散
   采用有限体积格式离散守恒方程，天然保证全局溶质守恒，避免传统显式格式常见的质量漂移。
4. 可视化与数据落盘
   内置两套实时绘图窗口：相场界面与溶质浓度，可每 500 步自动刷新；同时提供帧缓存接口，一键输出 GIF 或位图序列，方便后期对比与汇报。
5. 零依赖开箱即用
   除 MATLAB 基础环境外，无需任何工具箱；默认 300×300 网格、8 万步迭代在主流笔记本上 10 分钟内即可完成。

技术路线

• 相场方程：采用 Karma-Rappel 型可逆-不可逆分解形式，界面厚度与界面能独立可控。
• 溶质场方程：基于保守 Allen-Cahn/Cahn-Hilliard 混合格式，界面处溶质分配系数通过局部相场值线性插值。
• 各向异性：在梯度方向上动态计算界面法向，实时更新各向异性模数与导数项，保证数值稳定。
• 时间推进：显式欧拉法，时间步长自适应上限由界面扩散率与网格尺寸共同约束。
• 空间离散：九点紧支模板近似 Laplacian，显著降低网格取向误差；溶质通量采用迎风-中心混合格式，兼顾稳定性与精度。
• 边界条件：四边零通量，通过镜像点隐式施加，无需额外矩阵扩展。

使用方式

1. 启动 MATLAB，将 karmaAlloy.m 置于工作目录。
2. 直接运行，默认参数即刻开始计算；若需调整抽拉速度、温度梯度或各向异性强度，只需修改文件顶部参数区。
3. 运行期间，两张实时图窗自动弹出，可中途暂停观察；计算结束后，命令行打印耗时统计。
4. 若需批量实验，可用 MATLAB 的 batch 或 parallel for 循环包裹主函数，配合参数数组实现高吞吐扫描。

性能与扩展

• 单核性能：在 3.2 GHz 的 i7 处理器上，300×300 网格 8 万步约 8–12 分钟，内存占用 < 500 MB。
• 向量化优化：核心循环已按 MATLAB 列优先原则展开，实测较原始四重 for 版本提速 4× 以上。
• GPU 移植：相场与溶质场更新均为逐点 stencil 操作，可直接改写为 arrayfun 或 CUDA kernel，预计再提速 10×。
• 三维扩展：将二维 stencil 升级为 19 点或 7 点模板，配合三维 FFT 求解温度场，即可形成 karmaAlloy-3D，代码改动量 < 15%。

典型输出解读

• 相场图：+1 代表固相，–1 代表液相，界面位于 0 附近。可直观读取一次枝晶间距、尖端速度及取向偏差。
• 浓度图：色彩映射溶质浓度，界面两侧浓度跃迁即为微观偏析；通过提取中心线数据，可进一步计算二次枝晶臂间距 λ2。
• 日志数据：命令行输出的“Compute Time”可用于衡量不同参数下的规模扩展性，为后续集群部署提供依据。

小结

karmaAlloy 以不到 400 行 MATLAB 代码，浓缩了相场法在合金凝固领域数十年的理论精华。它既是教学演示的理想范例，也是微观组织预测的快速原型。开发者可在其基础上继续嫁接流场、应力场或多组分热力学数据库，形成多物理场、多尺度材料计算流水线。