我来分享一哈我的配置方法。总之就是把用户自定义的和官方的分开。省的乱了。至于命名，如果修改版本特别多，官方名加日期命名比较好。代码最好用版本管理器，还能加加些笔记标注每个版本做了什么修改，是从哪个版本来的。

alias ofd='source $HOME/OpenFOAM/OpenFOAM-dev/etc/bashrc;
export FOAM_USER_SRC=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/src;
export FOAM_USER_APP=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications;
export FOAM_USER_SOLVERS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/solvers;
export FOAM_USER_UTILITIES=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/utilities;
export FOAM_USER_TUTORIALS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/run/tutorials'

alias of22x='source $HOME/OpenFOAM/OpenFOAM-2.2.x/etc/bashrc;
export FOAM_USER_SRC=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/src;
export FOAM_USER_APP=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications;
export FOAM_USER_SOLVERS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/solvers;
export FOAM_USER_UTILITIES=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/utilities;
export FOAM_USER_TUTORIALS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/run/tutorials'

alias of30x='source $HOME/OpenFOAM/OpenFOAM-3.0.x/etc/bashrc;
export FOAM_USER_SRC=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/src;
export FOAM_USER_APP=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications;
export FOAM_USER_SOLVERS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/solvers;
export FOAM_USER_UTILITIES=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/utilities;
export FOAM_USER_TUTORIALS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/run/tutorials'

alias of4x='source $HOME/OpenFOAM/OpenFOAM-4.x/etc/bashrc;
export FOAM_USER_SRC=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/src;
export FOAM_USER_APP=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications;
export FOAM_USER_SOLVERS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/solvers;
export FOAM_USER_UTILITIES=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/utilities;
export FOAM_USER_TUTORIALS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/run/tutorials'

alias of5x='source $HOME/OpenFOAM/OpenFOAM-5.x/etc/bashrc;
export FOAM_USER_SRC=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/src;
export FOAM_USER_APP=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications;
export FOAM_USER_SOLVERS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/solvers;
export FOAM_USER_UTILITIES=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/utilities;
export FOAM_USER_TUTORIALS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/run/tutorials'

alias of6='source $HOME/OpenFOAM/OpenFOAM-6/etc/bashrc;
export FOAM_USER_SRC=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/src;
export FOAM_USER_APP=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications;
export FOAM_USER_SOLVERS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/solvers;
export FOAM_USER_UTILITIES=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/utilities;
export FOAM_USER_TUTORIALS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/run/tutorials'

alias of7='source $HOME/OpenFOAM/OpenFOAM-7/etc/bashrc;
export FOAM_USER_SRC=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/src;
export FOAM_USER_APP=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications;
export FOAM_USER_SOLVERS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/solvers;
export FOAM_USER_UTILITIES=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/applications/utilities;
export FOAM_USER_TUTORIALS=$HOME/$WM_PROJECT/$USER-$WM_PROJECT_VERSION/run/tutorials'

# default version
of6

#对应OpenFOAM自带的src lib app sol util tut.
alias msrc='cd $FOAM_USER_SRC'
alias mlib='cd $FOAM_USER_LIBBIN'
alias mapp='cd $FOAM_USER_APP'
alias msol='cd $FOAM_USER_SOLVERS'
alias mutil='cd $FOAM_USER_UTILITIES'
alias mtut='cd $FOAM_USER_TUTORIALS'

那个地址是运行程序临时分配的地址，也不是完整的物理地址，应该是某个段地址的偏移地址。如果要看，得用gdb调试，设置断点再看。
除非是相当底层的调试，一般用不到这种调试方法。

这是个ddtFoam()函数引起的浮点溢出错误，就排查排查程序参数有没有不合理的。

因为你肯定没回头去把手册或者视频看一遍

从哪算取决于starttime的设定和有没有那个时间节点的文件夹。

• 看是不是只能从0.8算，那就mapfields过来后，从0看看能不能算；算了之后0.8有没有更新。

• 至于看0用没用上，把0里的边界条件瞎胡改看看抱不报错不就知道了

高中学习和大学考试还是那种“我一恼”“我一生气”就做到了的模式。往后就是“临渊羡鱼不如退而结网”。重点是“退”，死磕的话。。。。不是有个成语叫“磕（客）死他乡”

• 17页2.1.4.5：Max Propagation长度为0.5 Maximum Step Length of 0.5

• 22页倒数第二段: Attribute Mode : Attribute Type。这个英文原版也错了。

• 比如22页倒数第四段开始。这段最后一句，说取消勾选。实际5.4中movingwall、fixedwalls、frontandback默认取消勾选。

• 22页倒数第二段。“X Axis Data我们选择Use Data Array以及arc_length，这样x轴代表了离底距离”。
  实际操作。“X Axis Data”现在叫“X Axis Parameters”(22页倒数第一段的"Line Series"也变了了"Series Parameters")。没有了“Use Data Array”选项（默认了），多了个“X Array Name”选项。也是默认了arc_length。

• （这个可忽略）然后22页倒数第一段。这里说“可以看见默认设置就已经有了速度的模以及三个分量”。如果按教程，操作到这里已经执行了postprocess命令。无法判断看到的U_x U_y是默认的还是命令生成的。（实际，默认或者命令生成的都是一个效果。命令并不改变效果。只是在每个时间节点文件夹多了三个文件）当然上面提到postprocess -func "components(U)"时说了“例如”，没说让操作。

• 22页倒数第一段，最后一句，把英文原文的“双击鼠标”翻译丢了。"The user can edit this most easily by a double click of the mouse over that selection"

• 23页第二段，改图中字体、颜色、名称不在图形左上方。在pipeline Browser 下方

• 23页第三段，“图2.11我们使用了下述选项：采用标准格式，指定了坐标轴范围”。原文是“Notation”选“Standard”，并选中“Specify Axis Range”。是软件名词，翻译成中文了。然后。英文版的也不对。没有Standard了。而且四个轴可以分别设置了。mixed、scientific、fixed。最后一句：enable line series，没有这个选项了。只要在line style里选none就可以不显示线。然后在marker style里选circle或者 square能实现图2.12效果。设置的时候还要把要调整的三条线都选中。

• 还有一个细节。英文原版在代码部分用的是“等宽字体”。一般代码都用这种字体，上下对的齐，变量名太长不至于看起来很稠。正常用来阅读的英文不等宽，i最窄ZW这些占的宽，符合书写和阅读习惯。

• 61页：3.2.2.1中“1.”“2.”中有三处出现“InInclude”。应该该是"lnInclude"（lninclude\LNINCLUDE）英文版也错了。

• 61页：3.2.2.1中"4."中"$WM_ARCH"前多了个"-"。英文原版那有个"-"是因为在行尾。

• 61页：倒数第二行，“EXE”和“INC”中间有个下划线丢失了。

• 62页：第二个"3."，“EXE”和“INC”中间有个下划线丢失了。

• 66页：倒数第三行，“FOAM”和“RUN”中间有个下划线丢失了。

• 66页：倒数第二行，“FOAM”和“SOLVERS”中间有个下划线丢失了。

• 67页：“来进行编译。编译过程会产生下列信息”下面的信息好像还是openfoam2.3的。英文版的手册这里已经变了。再下面的再次编译操作的信息也和英文版不一样。我操做的和英文版的也不一样。现在重复编译，不返回任何信息。

• 92页：量纲表第5个量是“物质的量”，单位“mol”。“摩尔质量”是另外一个单位kg/mol，表示1mol物质的质量。

• 120页：第二行“需要兴义网格单元形状”->“需要定义网格单元形状”

• 150页：倒数第二行，“首先需要创建一个文件系统”，英文版是创建个算例“create a new OpenFOAM case”

106页：量纲表第4个量，名称应该是“温度”，USCS单位应该是"。R"(Rankine是读音的一部分，不是单位的一部分)；
106页：量纲表第5个量，SI单位应该是“mol”，USCS单位也是"mol"

• epsilon是湍流耗散，而湍流耗散是 粘性产生的摩擦力 引起的 粘性耗散。一滩流体，需要考虑粘性耗散对流体的影响；一坨流体，粘性耗散大多数情况下可以忽略。而湍流充分发展就是个不可忽略的情况。

• 粘性耗散在大多数情况下，是固体以粘性为桥梁向流体施加影响（流体和别人摩擦）的副作用：动能转为热能。湍流耗散是自施加影响（流体自己摩擦）的副作用，是个内耗的过程（不团结一致向前冲，而把动能浪费在内耗上【注1】）。

• 这个内耗过程越严重，湍流耗散就越严重。那 内耗的量 取决于“摩擦力”（和粘性有关），“摩擦距离”（和湍流涡的大小和速度有关（决定涡微元的线速度）），“涡变形”（决定摩擦发生的频次）

• 那越粘的东西，涡微元转的越猛，涡变形越大，epsilon越大，我猜是这样的

【注1】自然界不是有个原理，耗能最小。反映在生物上，粘菌吃食物自规划最有效路线。反映在物理上，光的折射路线即是光经过折射面到达目的地的最快路线。那么流体转捩有点不符合这个原理，还是流动的目的不是向前冲，所以通过湍流达到目的。湍流好像是个熵增的过程，也是自发的。（一瓶子不满，半瓶子咣当中）

这更像是用代码格式化软件统一生成的。高级编辑器自带的代码格式化一般只改行前缩进，代码格式化软件是全部重新排版。

一般大型代码库都有非常细致的格式要求，代码格式化软件可以非常细致的定义每一个格式化要求：每行最多字符；算符前后加不加空格、换不换行，前换 后换还是都换……

然后用被处理过的代码非常满足每一个细致要求，所以某些地方就看着有一点点别扭了。 indent是GNU出的格式化软件。里面有一些常用的风格模板：GNU、K&R、linux。也可以跟“像一篇文章一样长的”参数来定义自己的代码排版风格。

@zhe DPMFoam.C 这个的108行在算粒子，下面就开始用算过的信息更新alphac了。

108行的定义在这里 KinematicClud.C 673行。这里用到的solve函数就在这个文件的91行。

然后这solve里一大堆追着函数名找，大概就是particle distribution function的过程了。我猜的

@zhe 我说错了，本来想说alpha1+alpha2不等于1。就是括号里的对立面

（如果alpha1 + alpha2 = 1，那么rho = alpha1*rho1 + alpha2*rho2，未找到alpha1和alpha2的关系。）

这个确实不能说 alpha1 和 alpha2 是未归一化的。
但是这个不影响后面的思路，虽然思路很模糊。
其实绝对正确的平均密度公式是：(alpha1 * rho1 + alpha2 * rho2) / (alpha1 + alpha2)。无论alpha1、alpha2有没有归一化都成立。

还有，刚才又回看了一下。你的第一个帖子说

密度rho應該已經經過alpha1*rho1 + alpha2*rho2得到的整體流體密度

如果这个表达式表达整体流体的密度 ，那是不是就意味着aplha1+alpha2=1?
如果 alpha1+alpha2\=1，那么剩下的体积分数代表的密度不是没算进来。但是代码中并未有限制alpha1和alpha2关系的地方。alphac倒是限制为1-xxx。

不考虑求解器。多相流的平均密度应该是连续相和离散相的平均alphac * rhoc + alphad * rhod，其中alphac + alphad = 1。

MPPIC方法在OpenFOAM里计算好像用的是particle parcel？所有都没有见 alphad。

我说的太外行就不用管了，我没研究过，就是插一句嘴。
附上一个某个帖子的总结：

忘了哪个帖子了，是这个论坛的。

引用自己的两段微信聊天记录：

结构网格和非结构网格的区别在于是否可以根据网格的索引编号确定相邻网格。 结构网格的索引号可以判断网格的相邻网格和面归属，非结构网格不可以，要另存储。是不是完全垂直只会影响算面通量要不要乘以面向量

最后一句换为是不是完全垂直只会影响算面通量要不要乘以带夹角的面向量。

是的，fluent用的也是非结构网格。如果用icem画过图，一般在导出到fluent时都会转换成非结构网格才能导出。因为非结构避不开，结构的可以转换成非结构的，所以基本都是基于非结构网格的。结构网格比较配有限差分，但是现在有限元更常见。约翰。安德森的《计算流体力学基础及其应用》主要讲结构网格配有限差分，要建贴体坐标系，要坐标系变形，主要例子是算飞机翼型，应该说应用范围没有有限元广泛。openfoam用blockmesh生成的网格，或者其他格式转过来能用的网格，在constant/polymesh里，就有定义面归属（owner）和相邻网格（neighbour）的文件，所以是非结构的。总之，fluent和openfoam都是基于非结构的，来了结构网格，先转成非结构才能算。

可能有理解错的，但应该能解除部分疑惑。

说两个猜想。
如果是alphc(fake)*rho = alpha1*rho1 + alpha2*rho2。
那得是alpha1、alpha2都是未归一化的值。
（如果alpha1 + alpha2 = 1，那么rho = alpha1*rho1 + alpha2*rho2，未找到alpha1和alpha2的关系。）
这种情况下rho = (alpha1*rho1 + alpha2*rho2)/alphac(fake)即密度=质量/体积

实际上代码里的rho(code) = alpha1*rho1 + alpha2*rho2，即上面的质量。
即rho(code)是体积为alphac(fake)的混合物的质量
即rho(code)是单位alphac(fake)的质量

下面开始猜：
alphac(code)是归一化过的，在代码里有显示，是1-xxx。
那么猜测：与alpha1和alpha2类似的数值alphac(real)隐含在alphac(code)中。
即alphac(code)=alphac(real)/alphac(fake)

那么rho(code) * alphac(code)代入以上定义就等于：(alpha1*rho1 + alpha2*rho2)/alphac(fake)*alphac(real)
即rho*alphac(real)
(rho(code) * alphac(code) -> rho*alphac(real))

上面是猜rho可能的定义。下面是猜有的地方乘alphac，有的地方不乘。

上面说了alphac是归一化过的，则rho(code)是单位体积内的，这种单位内的量容易搭积分的便车补上计算。对流项没有考虑alphac，因为对流计算过程涉及体积(速度即单位时间内的体积)。

至少把你的现有进展说一下，进行到哪一步了。particle distribution function 是代码里的还是文献里的，在代码哪里。你看了什么文献、哪个代码却没有找到。好歹让想帮你的人把你踩过的坑给避开了。

用 sudo fdisk -l 看看 swap 分区多大。如果很小或者没有，可能是因为内存。如果很大可以排除内存问题。swap是虚拟内存，物理内存不够用会把硬盘当内存用，只会很慢不会崩。

也可以在编译过程中开着 system monitor 看着内存和 swap 使用状况。

当初装linux的时候分区选择了lvm方式，还能调整swap分区大小。如果不是lvm而是直接分区。只能重装系统修改分区大小。

然后你的源代码放在OpenFOAM1，1??? 不知道会不会影响。

估计你的源代码是用apt下载的，或者是哪位大佬的祖传代码，所以默认安装位置改成了/opt。

各种教程的编译安装位置基本都在$HOME\OpenFOAM，不会有权限问题。
调整 OpenFOAM 的安装位置还是挺麻烦的。

@RebelYoung 这种错误可能是数据本身有问题，不一定是paraview的问题。

如果想确定，可以去paraview的官网下载最新的paraview。如果连它都崩，那就是数据的问题。

https://cfd-china.com/topic/3453/运行pisofoam的时候-卡住不动了

里面有个可以尝试的方法。

在求解器的各处随机插入
cout << "run in line xx"
重新编译这个求解器，看看具体在哪一步之后就不输出了。

根据描述，断点续算可以，但是直接突破不过去，而且只是不输出而不是不计算。应当是某个变量积累超过输出语句的判断范围。所以还在算但不输出。所以上面提到的具体哪一步就是修改了这个变量超出了范围。

更精准的定位错误就要看输出残差的函数怎么写的了。

OF 的 debug 模式编译有更强大和方便的调试功能。