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

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。也可以跟“像一篇文章一样长的”参数来定义自己的代码排版风格。

是不是模型的xyz搞翻了，可以先用paraview看看。

或者（16 2 1）试试

https://www.cfd-online.com/Forums/fluent/40645-non-newtonian-fluid-fluent.html

@深蓝 你的问题已经超出我的能力范围了
但是这里有个一个很完备的例子，应该对你有帮助。

一个简单而完整的UDS例子.zip

线性化，我能想到的方法就是泰勒级数展开了，不过要求值的变化范围不能太大，因为泰勒展开有个展开点，离展开点越近越准。这个过程应该是在软件外进行的，推出具体的计算方程再往UDF里写。

只能帮你到这了

@深蓝 这个网格数量……我猜应该是网格面积。这里刚好你的网格都一样大，不信你换个网格密度试试

这个材料密度，按说uds就是在场内运输个标量，跟其他的不太相关，但是一般应用中都会是某种物质，和其他物质有相互作用，然后就得有质量属性。这里材料密度能影响到uds的总量，应该是有个默认的1：1的关系，定义uds那里有个flux function，应该选的是mass flow rate。感兴趣可以找例子在udf写个flux function，兴许uds总量和材料密度的关系就没了。

这patch都解决问题了，怎么还倒腾udf

@深蓝 可以做这样一个测试，全部区域都给这个源，看浓度是不是线性升高的就知道了是哪种了。

如果只是放置个污染源看扩散，就在计算初始化的时候设定。这样应该比较简单。
大部分计算初始化都是allzone 设定为 0 。可以在patch里分区域设定计算初始状态。
面积乘以UDS值，就是初始污染物的量。想要单位1的污染物。就给单位面积上设定uds的初始值为1。

左上角的pipeline browser里的每个条目前面都有一个眼睛图标，挨个点击隐藏，看白箭头是哪个的

“Error: Invalid result name”

就是保存的太频繁了，电脑忙懵了。也可能是你把所有的核全用上了。按说应该留两个给电脑喘喘气。
或者你的c盘满了？

那个设置的值就是“值”，不是“通量”。就好比，设定一个一百度的铁球在那里，但是不是说把一百度铁球上的热量全都放到场内。

如果环境温度800度，这个100度的铁球还会吸收热量。吸收了应该升温，但是每次迭代都会又设定为100度。所以不停地吸收。

如果环境温度0度，这个100度的铁球持续释放热量。释放了应该降温，但是每次迭代都会又设定为100度，所以不停的释放。

是否吸收释放，以及吸收释放的速度取决于温度差。这里应该是浓度差。时间步越小，越能一直保持较高的浓度差，释放的就越多

就好比，给了个压力入口边界条件，却想得到固定的通量。

@conghongchuan 那多算一会儿呢。这个由出口向内的速度是压缩气泡时出现的。按说气泡大小稳定后，这个由出口向内的速度就该没了，否则，只进不出，不是炸了。然后，重力的作用就凸显了，大密度的液体会下沉，小密度的气泡上浮。

刚开始重力没什么作用，因为气液压差更大。

这样看起来，你的计算结果也是可能出现的，可能和论文中的计算条件并不一样。最后一张图，论文中的气泡体积明显大于模拟的气泡体积。可以多跑些时间，看看气泡整个变化过程。

@conghongchuan 那就剩重力加速度的方向设置可能有问题了
把g分别设置四个方向，看是不是每次气泡都逆着重力跑

还不行就很诡异了