口语有一定的随意性，这种随意性会和书面语有细微的容易忽视的区别。比如。脖子涡，其实“子”是常见后缀，也可以叫脖涡，但是无论叫什么，这脖子涡是手脖子还是脚脖子？万一以后有“手脖子涡”、“脚脖子涡”，那“脖子涡”将何以自处！
“颈项”准，但是“项”多余。其实叫“颈涡”就行，但是不方便口语交流。加上“项”让人更容易联想到“颈”，而不是“净、“睛”之类的。
参考诗分段：
首联，颌联，颈联，尾联
首涡，颌涡，颈涡，尾涡
如果是一堆小涡聚集可以这样叫
如果是单个涡拉长，可以叫
涡首，涡颌，涡颈，涡尾。
看看，那些口口声声“涡头”“涡尾”的人多么的粗俗！应该叫“涡首”！
别急，等你当导师的以后，你也可以对你的学生说：
https://www.openfoam.com/documentation/user-guide/standard-boundaryconditions.php
这是所有的边界条件。fluent的这个uds，可以把它类比成某种物质，比如空气，这样设固定值就是在设浓度，设流量就是在设速度。挺方便。OpenFOAM就要自己按需配制了，比如固定值的速度就是固定流量的空气，但是固定值的浓度就是固定流量，固定梯度的浓度是固定流量。
做一个点线图，横轴是孔隙率，纵轴是最大液相分数，看看到哪里突变了。二分法找这个突变点，应该很快。
出图的时候用同一标尺（最大值最小值相等），这样不同图之间也有可比性。看看那个刚好导致突变点的孔隙率的case，液相分数是怎么消失的，是跑出边界了还是在流场内被计算过程吃了
@starDust 已经成功了？在tecplot里插值就可以显示了？
按说tecplot和paraview应该都可以，快捷的可视化数据。其他的方法就麻烦了，用各种画图工具比如gnuplot去画
应该是可以不一样的。甚至就应该不一样，否则很容易混乱和覆盖。
这里改了，用到这个库的地方也得改。
这里不应当是个错误或者警告。可能有些后续操作没跟上所以出错或警告。但这个“不同命名”的操作方向是对的。
网格变形成这样，肯定算不对啊。
所以找问题得找：为什么网格变形成这样。
如何不让网格变形成这样。
至于发散、库朗数之类的，得在网格正常的基础上讨论。
看起来是圆柱的移动逼近了极限。所以要么限制圆柱的运动（离深蓝色和浅蓝色交界处原点），要么扩大运动的范围（就是可以变形的深蓝色区域）。
可以找已有的代码对比一下，除非特别新的东西，常用的东西一般都有模板，还有各种代码的模板。
http://www.pudn.com/Download/item/id/3745439.html
https://github.com/Singh-T/TDMA-2D
这些八成是对的，把他们的代码一点一点的替换成你的代码，看看替换到哪就不对了
The Lagrangian scale-similarity model:
https://www.cfd-online.com/Forums/openfoam-solving/90351-dynlagrangian-sgs-model.html
The Lagrangian scale-dependent model:
https://core.ac.uk/download/pdf/30828294.pdf
对着引用文献应该能找到很多
颗粒属于多相流了吧。组分一般是同相的，都是气体，或者都是液体，然后比如空气，百分之七十氮气，百分之二十多氧气。他们相溶在一起，无法分开。但颗粒和周围并没有相溶。即使是同相的，油水混合物，颗粒内的物质和颗粒外完全不同。组分好像没有这个特征。
即使颗粒特别小，计算中用相分数表示，感觉和组分还是差别挺大的
可以找一下openfoam的旋转机械案例。还有船舶的模拟，螺旋桨转着推着船跑的那种，用的重叠网格，overset
还有这个帖子
https://www.cfd-china.com/topic/3921/openfoam中被动变形下的主动控制技术问题-柔性翼形俯仰/6
这是在计算时（fluent或者openfoam）设置的，边界条件。画网格的时候形状对了就行了。
前后上下左右都是周期性边界条件的也有啊，这个当然可以。
应该是这样的。
can@x260-TC ~/.local/share/OpenFOAM/OpenFOAM-8 %echo $WM_OPTIONS
linux64GccDPInt32Opt
按说没加载openfoam的环境，应该调用不了求解器。可是加载了环境就应该有这个环境变量。
这重的也太猛了……重了五个，如果说解决问题是开锁，这个锁眼都快堵死了
这是长期乱套导致的，越往后越纠缠不清。这问题太细致，一个操作都能写出来几个分叉，操作多了就说不全了。
好习惯建议：
https://cfd-china.com/topic/3190/分享一个在任意文件夹中编译库或求解器的办法
这个里面我也跟了个帖子，是写在.bashrc里的，现在又改良了，就是下面这个代码，单独成一个文件，在openfoam里的etc/bashrc里加载。可以把代码和算例分开放，备份方便（被屏蔽的前三行）。还增加了openfoam里没有的快捷命令，可以把用户的东西（代码、二进制、算例）与自带的完全分开。
#\$export FOAM_INST_DIR=\$HOME/.local/share/\$WM_PROJECT
#export WM_PROJECT_USER_DIR=\$HOME/Documents/openfoam/\$USER-\$WM_PROJECT_VERSION
#export WM_PROJECT_USER_CODE_DIR=\$HOME/Documents/code/C/openfoam/\$USER-\$WM_PROJECT_VERSION
#-------------------------------------------------------------------------------
export FOAM_USER_SRC=\$WM_PROJECT_USER_CODE_DIR/src;
#-------------------------------------------------------------------------------
export FOAM_USER_APP=\$WM_PROJECT_USER_CODE_DIR/applications;
export FOAM_USER_SOLVERS=\$WM_PROJECT_USER_CODE_DIR/applications/solvers;
export FOAM_USER_UTILITIES=\$WM_PROJECT_USER_CODE_DIR/applications/utilities;
#-------------------------------------------------------------------------------
export FOAM_USER_TUTORIALS=\$WM_PROJECT_USER_DIR/run/tutorials
alias msrc='cd \$FOAM_USER_SRC'
#-------------------------------------------------------------------------------
alias mapp='cd \$FOAM_USER_APP'
alias msol='cd \$FOAM_USER_SOLVERS'
alias mutil='cd \$FOAM_USER_UTILITIES'
#-------------------------------------------------------------------------------
alias mtut='cd \$FOAM_USER_TUTORIALS'
alias mappb='cd \$FOAM_USER_APPBIN'
alias mlibb='cd \$FOAM_USER_LIBBIN'
加载上面的代码的位置放在下面这个地方，我是连安装位置一起修改了，所以看起来比较乱（只看第三行，前两行是另一个修改，可以定位）。上面的代码我放在my_enviroment文件夹里，每个版本都有各自的，大体相同，不同处可以分开修改
-export WM_PROJECT_USER_DIR=\$HOME/\$WM_PROJECT/\$USER-\$WM_PROJECT_VERSION
+export WM_PROJECT_USER_DIR=\$HOME/.local/share/\$WM_PROJECT/\$USER-\$WM_PROJECT_VERSION
+source \$FOAM_INST_DIR/my_environment/\$USER-$WM_PROJECT_VERSION
openfoam自带的快捷命令在etc/config.sh/alisas，自定义的环境变量在etc/config.sh/setting。但是并不是每个环境变量都有带user版本的，快捷命令更没有。
最好是不要有重名的库，不重名，其实放不放user都行，但是分开放更好。无论是备份自己的工作还是不弄乱原有的代码。
还有一个好处就是，解决这种重复。即使库名字不冲突，但库内容冲突。比如不同的库中重复定义了同一个类。这时候把编译时要链接的库地址删了带冲突库的路径，这样就不冲突，甚至会缺。缺的单独拉到user里就行。
这，还带顶帖的…… 你一顶，大家一看有回复了，好，下一题。
contour 的范围是正负七百多。可是2里的 contour levels 才正负15？
@hongjiewang 我也编译过 debug 模式的 of ，就是比正常慢点，没出现过这种错误。这种错误看起来很初级不像是源代码带的。我搜我的代码的这个文件也没这句。
单这个错误，是的，加个分号就行。就是一个语句忘了结尾。
@hongjiewang 第二个错误信息我看错了，确实不是都和thermo有关。
可是你的第二个错误你说 g 出的错，
Info<< "Need to store rho for ddt(rhoCp, U)\n" << endl;
这一句后面只跟一个 “1”，在代码里不是这里么
Info<< "1" << endl;
#include "createMRF.H"
#include "createFvOptions.H"Info<< "1" << endl;
#include "readGravitationalAcceleration.H"
createMRF.H 或者 creat FvOptions.H
输出第二个“1”才到g
250的网格都错，和内存无关了。这种面向对象的程序一般写的时候都有脚手架，就像官网的代码有很多框图，但又比那些框图要详细。这留个窟窿那钻个眼都有用，可是只看代码非常费功夫才知道有什么用。不知道有什么用就很容易错。互相之间的依赖关系，先后顺序。那些不报错的是没语法错误了，不一定没语义错误。所以一般都在原有的类似的求解器上增删改；改变结构造个新的，还是挺麻烦的。of的路慢慢熬吧。
如果还不行，用gdb调试，有gdb的of版，gdbof
这两篇高引和一本书应该挺足的了
这个 ALE 用谷歌搜还是很多资料的，如果不会翻墙可以用 bing 的国际版搜 Arbitrary Lagrangian and Eulerian
openfoam的过程，可以看东岳流体网站上的求解器解析。
曾经，我排过一个很狗血的错，一台电脑的某个软件也是这样捉摸不定的出错，不知道在哪里，也没有一点点防备。找了半天发现是磁盘满了。还剩一丢丢空间，所以一开软件，这个空间什么时候被临时文件占满了，就崩了。错误特征完全无法捉摸。
这个，可以du -h看看各挂载点的空间。然后，开两个窗口。一个开top，监视内存；一个开求解器，看内存的实时占用，不过有捕捉不到的风险。可以把网格尽量缩小再试试。这些变量占的空间基本都依托于网格。
还有一个点，这两个错误都跟 thermo 有关。thermo好像也和湍流有关
https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1ThermalDiffusivity.html
按这个图，肯定是难求的放内循环，难求的求不出来循环自己就行了。
如果好求的放内循环，外面难求的每次求不成功，再来一圈还得把好求的带上再求一次。
https://www.openfoam.com/documentation/guides/latest/doc/index.html
https://cpp.openfoam.org/v8/
在这搜，每个版本的都不一定一样，先找对版本，然后搜。其实在自己电脑上也能用grep搜，但是这个看着跳转更方便。
里面重名很多，按求解器或者后处理工具可以辨别你用的哪个。
入口处的压力小，因为出口在“泄压”。极限条件，出口处堵死，入口处压力就会增的特别大。
所以当设定为固定压力时，会大大增加流速。
总之就是阀门对流体的阻滞作用弱，所以导致压力上不去。
而你的模型里，阻滞流体的部分主要依靠狭窄的通道。（极限条件，那狭窄的通道特别宽，以至于出入口根本就是联通的，等压的）
而狭窄的通道可以阻滞流体，直接分析，这流体得足够粘，通道的壁面得足够粗糙。阻滞的作用才大，才能在入口处把压力憋上去。可能模拟相对于实验：壁面太滑、粘度太低。也可能是网格在那狭窄的地方太糙，以至于没有反映出阻滞效果。
边界设置好气温不就行了。传热就三种途径，辐射、对流、传导，空气的传导热作用比较弱。鸭绒就是通过绒毛阻止空气流动形成空气墙保暖的。边界处的空气除非进入流场参与热交换，也不会有辐射（假设太阳在头顶……），传导又弱。
没发现解释就自己解释啊，解释清楚了往论文里写
@Jacob-Z 搞不明白，或者说不能分开搞明白，这是个知识网，交错分布，不能一小块一小块的搞明白，否则老得回头看回头搜，时间成本更大。
OpenFOAM 并行和通常意义上的并行有点区别。一般并行的部分都是写在代码里的。OpenFOAM 的并行部分不在代码里，严格说，都不是并行的。并行的任务分配是通过 decomposePar 完成。并行计算是通过 runParallel 完成的。所谓并行就是一下开几个求解器进程，同时算不同的区域，几个求解器的计算区域的交接处信息交换也是 runParallel 完成的。
这样求解器中的变量就不能“特殊”，要“通识”。不能直接“全场”信息，要用“我计算的所有网格”信息。
你的第一个帖子说要用每个时间步计算的到的数据，这个数据的调用直接写进求解器。那么在各个区域上计算的进程有的快，有的慢。全场的信息（或者说某个时间步计算的到的结果）对于单个进算进程来说是未知的。对于单个进程：我只计算左上角的区域，然后下一步计算需要引用全区域变量。
runParallel 应该是没有这样的调度的。进程有快慢，所以偶尔会互相等待。等待也是等待交界面的信息。而这个全场信息是等不来的，得合并各个进程的计算结果才有。
所以有些求解器它就是不能并行，施主又何必强求
https://www.dataphysics-instruments.com/us/knowledge/understanding-interfaces/surface-tension/
这里有一句
Thus, for liquids the surface tension and the surface energy are equal
至于动能，k就是湍动能啊，如果不只是要湍动能，那k怎么来的肯定要涉及动能啊。
wiki是个大类，维基百科是一种实现形式。也有轻型的，比如 MediaWiki。要说记录知识，wiki就是为此而生的
一般 wiki 都是网页形式的、服务器模式，因为一般用wiki这么规范的都是很多人一起编辑的，给很多人看的；所以跟服务器网页的模式很搭。个人笔记用wiki有点重。
wsl是win模拟linux的api，就像wine一样。wsl2就是个虚拟机，不过经过微软优化，没其他虚拟机那么嗝应人。
如果是升级失败，系统里应该有还原点，还原一下更新前的系统状态就好，不会擦除其他个人数据。
putty连ssh，filezilla传文件，xming接图形界面，都是开源的。
Xming 也能接本地的图形界面，比如wsl的。
Xming在软件下载网站softonic下载。
Xming安装用默认的。然后运行,它会显示即将接收到的图形信号的标识。标识就是display number 的编号，默认的是0。
Xming之所以能够接收到DISPLAY的信号 , 得让wsl知道往这里发送信号。让wsl往这里发送信号的操作就是在wsl系统中的 ~/.bashrc 文件中添加下面一句。
export DISPLAY=:0
然后保存 .bashrc，重启wsl系统。这时，当wsl系统运行图形界面程序时
(如xclock)，Xming 会接收到信号。进而显示。
@merlinizlf 当然是用画网格的软件检查，检查完不是有个质量系数，应该也有个柱状图，网格质量的区间分布，应该能设定显示某个质量区间的网格，比如icem，点质量最低的柱子，就能显示了。
按说网格质量质量应该靠预期控制的，比如在某些结构就用某类型网格，否则就会变差，比如圆弧上的O形网格，分叉上的Y形网格。如果用了还质量差，那就只能加密了。
加密会提高质量，所以，让网格变稀疏些，会暴露可能的低质量网格区域
@zhe 原来你是疑惑怎么从网格里提取 parcel 啊。
https://www.openfoam.com/documentation/guides/latest/api/PairCollision_8C_source.html
你放的代码是这个89行吧，84行创建了cellOccupancy，
那个创建容积率的函数在多个模型下有同名的，我猜这里这个是
https://www.openfoam.com/documentation/guides/latest/api/KinematicCloudI_8H_source.html
这个的371行，
这个函数就是判断是否存在，否则创建。就是这个函数buildCellOccupancy()，在
https://www.openfoam.com/documentation/guides/latest/api/KinematicCloud_8C_source.html
这里140行
你给的代码里 pA_ptr, pB_ptr 都是通过 cellOccupancy 的返回值得到 parcel 的指针，就是 那个371行，如果存在就返回
return *cellOccupancyPtr_;
而这个 cellOccupancyPtr_ 也有多个定义，在这里应该是这个
https://www.openfoam.com/documentation/guides/latest/api/KinematicCloud_8H_source.html
173行。
173行的注释大概解释了cell 和parcel怎么联系起来的。
好了
模板继承用的太多，我不太熟悉这部分，同名函数定义太多，剩下的我也没招了
"时间相关系数的积分"，这个，在Durbin那本书里应该是 integral correlation time-scale，而不是 Integral of time correlation coefficient 。
书中这句：
The integral correlation time-scale is regarded as that on which the turbulence equilibrates and is of order k/ε.
再去 bing 或者 google 上搜 correlation time scale，有很多讨论。
刚添加点，就发现又回复了……我已经被榨干了……我不是搞这个的，
请洒潘江，各倾陆海云尔
然而目前学者提出的各种积分尺度计算方法，个人觉得基本不是遵循原始定义计算，这样产生的问题是不同方法计算的结果区别很大，更像是一个“算不准”的参数。用一个“算不准”的参数去分析问题，总觉得这是有点奇怪的。
这个不叫算不准，叫计算标准不统一。就像同样是表征温度，有华氏，有摄氏，有开氏。温度的原始定义是分子内能，这也不能直接测。但也不妨碍使用各种度分析问题。
一直没明白空间相关系数的积分怎么能代表涡的平均尺度大小？
不是“积分”代表了“尺度”，是“积分”反映了“幕后主使”的状况，“幕后主使”决定了“尺度”。就像看面相，不是某些特征代表他是什么样的人，是他的为人性格决定了他的面部特征和他是什么样的人。
当然这个"幕后主使"是什么，就靠大家挖掘了……
“时间相关系数的积分”这个统计量，在数学书中未看到过，好像基本只在湍流的书有提及。
某某的积分在数学书中有定义吧，比如“未知数x的积分”。现在只是把“未知数x”换成了“时间相关系数”，他们都是“数”，如果还是单独定义"时间相关系数的积分"，那“空间想关系数的积分”、“愉悦感相关系数的积分”、……能说三天三夜不带重样的。
1 2 有很多我不懂的名词，插不上嘴
这里不是看不起中文教材，主要是，不是特别特别特别经典的著作且发展并不迅速；同一个英文词汇很难有统一的中文翻译。就很难在pdf里搜，用眼找太慢了还容易错过。上面的图片是 Pope的 turbulent flows 的，还有一本 Durbin 的statistical theory and modeling for turbulent flow 有很多相关内容。
分两部分显示。用 paraview 里的 find data 把大于Ystoic的值的网格筛选出来，再把小于 Ystoic 的网格筛选出来。然后只对大于 Ystoci 的那部分网格做 calculation 处理。然后同时显示两部分网格。
这段一看 A B 就是网格。但是把 A B 理解成 parcel 相关也不算没有依据，parcel 好像就是从网格上划分出来的。比如，进口面有10个网格，那么每个网格喷出来1个 parcel。一次喷出来10个parcel。
如果把 parcel 理解成单位时间内单位空间上的粒子物理量。那理解成网格还是parcel的区别不就是欧拉和拉格朗日的区别。
我猜的