哈工大传输原理之相场模拟.pdfVIP

传输原理之相场法 6121900604 李晨星 本次作业意在了解相场法，并通过相场法完成晶粒长大的模拟。 模拟由本人与张九相共同完成。本人曾独自完成过问题分析，差分格式构建以及编程求解。 电脑故障又懒于重新编程，最终的程序由张九相编写，结果也来自张九相。张九相在编程过 程中也有一些独到的想法，这会在他的报告中得到体现。 一、对相场法的评价 本节涉及到本人对相场法的基本态度，因而放在此处进行说明。相场法通过能量的变分对微 观过程进行模拟，是一种很了不起的方法。但在相场法中有太多的唯象参数，而其本身又没 有很好的确定手段，可信性是有限的。此外，通过数值解法对方程进行求解，随着时间的推 移，解会逐渐偏离给定方程，也并不太理想。 相场法之所以会得到广泛的推崇，可能其一是在于其将物理问题数学化，其二是在于其得到 的结果与真实结果很“像”。不过在编程过程中觉得这只是一个简单的数学游戏，完全看不 到其将对生产生活带来的指导意义，这可能与本人的见识浅薄有很大的关系。随着日后的学 习，可能会对相场法有更进一步的了解，但在此持保留意见。 二、晶粒长大的相场法 晶粒长大的相场法由朗道给出，引入了序参量集来表示晶粒的取向。序参量集中的任一序参 量η应在-1 到1 之间。描述晶粒取向时，序参量集中应有一个序参量为1 或-1，而其他序参 量为0 。晶界处的序参量没有明显的特征。通过构造局域自由能密度函数以及进行变分，朗 道得到了如下方程 ∂η 3 2 2 = −L −aη + −aη + 2γη η −k∇ η , = 1,2, ⋯, p ∂t ≠j 其中L、a、γ 以及k 均为唯象参数。并未对方程的建立进行过详细的学习，因而对其真正想 表达的思想并不能完全琢磨透。本文所要做的只是对这一个方程进行编程求解。 三、唯象参数的确定 采用一种并不太严谨的方法对上述方程进行分析，希望得到一些对于唯象参数选择的启发。 这里将右端项拆开，分别考察其分别可能对动力学方程造成的影响。 3 ∂η ∂t = −L −aη + −aη 一项中，共有三个奇点，分别为0 、1 以及-1。其中0 是不稳定 的，而1 与-1 是稳定的，这说明其会使所有序参量稳定到1 或-1。 2 ∂η ∂t = −L 2γη η 一项并不太容易分析，不过可以轻易的看到其可能会使序参量达 ≠j 到晶粒的规定形式。一个序参量为1 或-1，而其他序参量为0 ，是一种稳定的解。 2 ∂η ∂t = −L −k∇ η 一项是一个扩散项，与传热传质的形式很相近，其保证了晶粒有长大 的趋势。这一项是极为关键的。 事实上，前两项描述的主要是晶粒自己的特性，而后一项描述的是晶粒间的交互作用。这样 的分析方法显然是不正确的，但起码会为唯象参数的调整提供一定的依据，不至于在唯象参 数调整过程中无从下手。 四、差分格式的构建 这里研究的是一个二维问题，其差分方法与一维问题有很大的不同。差分格式由自己构造， 构造方式如下。 2 2 ∂ η 1 ∂ η 1 = − 2 + = − 2 + ∂x 2 ℎ2 4 5 6 ∂y2 ℎ2 2 5 8 规定357 为方向n1，规定159 为方向n2，由方向导数的相关理论得到： 2 2 2 2 2 2 2 2