数字旗手

写在前面——如何阅读OpenPhase源码 OpenPhase的源码真的是挺难读的，它是一整套完全的很成熟的求解相场模型的框架，其模块化完成度很高(尤其设计的数据结构面面俱到，号称能处理所有类型的数据)，所以用它来解方程时写的代码很少，但带来一个问题是你必须先读懂它：高度集成的代码让初学者“一进去就出不来”(从而迷失在各个类中)或“压根进不去”(一行代码就完成一个大功能，根本不知道怎么入手)。 我的建议是： 1. 善用Doxygen OpenPhase支持Doxygen，其注释也较全面，使用Doxygen生成类之间的关系图，方便理清思路。 附录有Doxygen生成的dot图

持续更新中~~~ 条件断点 为断点设置条件表达式，只有条件满足时才激活断点 如： break 670 if sum == 3 转自这里 出现”value optimized out”错误 gdb调试程序的时候打印变量值会出现value optimized out情况,可以在gcc编译的时候加上-O0参数项,意思是不进行编译优化,调试的时候就会顺畅了,运行流程不会跳来跳去的,发布项目的时候记得不要在使用-O0参数项,gcc默认编译或加上-O2优化编译会提高程序运行速度. 转自这里。

Update 2016-11-16更新： 在Ubuntu10.10上安装gcc4.8.0时，出现错误： 1 ‘CHAR_BIT’ was not declared in this scope 解决方法： 1 unset LIBRARY_PATH CPATH C_INCLUDE_PATH PKG_CONFIG_PATH CPLUS_INCLUDE_PATH INCLUDE 如果在/etc/profile 或者 ~/.bashrc中设置了ccache的话,也暂时去掉ccache的设置. 参考自这里。 另外： 单纯修改gcc的指向后，会造成其他编译器的不匹配，比如gcc虽然指

写在前面：本文转自原AFEPack论坛afepack.org。 AFEPack教学系列转载自北京大学数学学院李若教授之个人博客。 如果需要转载本教学系列，请保持文章完整性，并注明转载出处。 李老师的话 学校的服务器老是有问题，希望这个地方能够稳定一些。我把AFEPack的使用写一些系列的小例子，从最简单的开始，希望能够为我这个软件包的使用者提供一些方便。 我在大年三十的晚上更新了AFEPack 的主页上的一点内容，现在主要是照抄下来，给想用AFEPack的人一些方向性的指导。在过去的几年中间，我多少花了一些力气来推广这个包，希望能够为计算数学的发展做一点事情。由于一些令人失望的经验，我对

模型概况 该模型是Appolaire等人于2008年在MSE上发表。 B. Appolaire, H. Combeau, G. Lesoult. Modeling of equiaxed growth in multicomponent alloys accounting for convection and for the globular/dendritic morphological transition, Mater. Sci. Eng. A 487 (2008) 33-45. 该模型是一个半解析半数值的模型，一方面对于枝晶的细微结构如一次枝晶尖端半径、二次枝晶间距及直径，对流效应

安装AFEPack软件包，需要提前安好的软件有gcc、g++、doxygen、automake(1.11版本)、dx、emacs、vim。 准备： 我会有一个用户名，比如为qixinbo。在home/qixinbo下建立files，include，lib三个目录。其中include目录是用于存放程序编译时调用的.h头文件，lib目录是用于存放动态链接库.so文件 然后将 boost_1_37_0.tar.bz2, deal.nodoc-6.3.1.tar.gz, AFEPack-snapshot.tar.gz, easymesh.c.gz 四个文件放在files里。 安装boost： 安

纯铁在不同条件下的平衡态计算 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 r t steel1 //read TDB steel的缩写，读取steel1.TDB文件 fe //选择Fe元素 l data //list data的缩写，列出所有的热力学数据 set c t=1000 p=1e5 n=1 //set condition的缩写，设置温度K，压力Pa和摩尔数mol c e //calculate equilibrium的缩写，计算当前条件下的系统平衡态 l,,,, //list的缩写

本文是对徐祖耀老师所著的《相变原理》(科学出版社2000年出版)第一章的读书笔记。 相变的定义 在均匀单相中，或在几个混合相中，出现了不同成分或不同结构(包括原子、离子或电子位置及位向的改变)、不同组织形态或不同性质的相，就称为相变。 相变何时发生 当母相失稳而新相具有较高的稳定性时，就会发生相变。 母相失稳可以需要跨过一定的能垒(需要一定的形核驱动力，即先形核)，也可以连续失稳(没有能垒的连续型相变，不需要进行形核过程，如Spinodal分解)。 在一些固态相变中，除温度因素外，压强和应变也引起母相的失稳。 固态相变热循环时往往发生较大的热滞，即冷却时的相变温度不等于再次加热时的相变温度

传统Ginzburg-Landau自由能泛函形式 其中$\gamma$是传统明锐界面模型中的表面张力，$\epsilon$是界面宽度。 多相体系的混合能形式 其中$(\frac{\psi-1}{2})^2$这一项是为了保证两个不同相(标记为$\phi=1,\psi=-1$和$\phi=-1,\psi=-1$)之间的相互作用不直接影响第三相(标记为$\psi=1$)。 多相运动体系的总能量 通过在系统中加入流体方程，得到整个流体动力学体系的总能量，其是动能和混合能的加权之和： 这里$\lambda$表示两种能量之间的竞争。 应力张量 通过虚功原理求得Ginzburg-Landau能中的应

VTK文件格式概况 VTK，全称为Visualization Toolkit，即“可视化工具箱”，制订了一个统一的文件输入输出格式，这样就可在不同的软件间方便通信。 VTK文件格式包括五个基本部分： (1)第一部分是文件版本说明： # vtk DataFile Version 3.0 (2)第二部分是文件头，是一个由”\n”结尾的字符串，最大为256个字符，如： PhaseField (3)第三部分说明文件的格式，是ASCII或BINARY，两者必填其一。 (4)第四部分是数据集的结构。 该部分起始是关键词DATASET加上数据集的类型，数据集的类型，即几何/拓扑的类型，包括()：STRUCT