数字旗手

电气化、自动化、数字化、智能化、智慧化

0%

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

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

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虽然指
Read more »

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

模型概况 该模型是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. 该模型是一个半解析半数值的模型,一方面对于枝晶的细微结构如一次枝晶尖端半径、二次枝晶间距及直径,对流效应
Read more »

安装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: 安
Read more »

纯铁在不同条件下的平衡态计算 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的缩写
Read more »

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

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

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