CFD中文网

bestucan

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

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'

bestucan

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

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

bestucan

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

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

看是不是只能从0.8算，那就mapfields过来后，从0看看能不能算；算了之后0.8有没有更新。
至于看0用没用上，把0里的边界条件瞎胡改看看抱不报错不就知道了

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

bestucan

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”

bestucan

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

bestucan

epsilon是湍流耗散，而湍流耗散是粘性产生的摩擦力引起的粘性耗散。一滩流体，需要考虑粘性耗散对流体的影响；一坨流体，粘性耗散大多数情况下可以忽略。而湍流充分发展就是个不可忽略的情况。
粘性耗散在大多数情况下，是固体以粘性为桥梁向流体施加影响（流体和别人摩擦）的副作用：动能转为热能。湍流耗散是自施加影响（流体自己摩擦）的副作用，是个内耗的过程（不团结一致向前冲，而把动能浪费在内耗上【注1】）。
这个内耗过程越严重，湍流耗散就越严重。那内耗的量取决于“摩擦力”（和粘性有关），“摩擦距离”（和湍流涡的大小和速度有关（决定涡微元的线速度）），“涡变形”（决定摩擦发生的频次）
那越粘的东西，涡微元转的越猛，涡变形越大，epsilon越大，我猜是这样的

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

bestucan

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

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

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

bestucan

一个是锥子喷射，一个是锥子孔喷射。

可能是：

锥子喷射能保证喷着时是锥形分布。

锥子孔空射，是保证出口形状是锥子孔，至于能不能喷射成锥形，得看喷射压力和喷射流量了。

大口径小流量，是会单边喷射的，就像洗澡的花洒没水了。。。。

bestucan

应该和左下角的 tolerance 相关。
这tolerance都那么小了，还能12点覆盖

估计是把一个几何体修复成一个点了，所以崩了

bestucan

@dxl 不用放时刻文件夹。文件名中的001 002 003就被当作顺序了

bestucan

@dxl

这个tableToPointsFilter可以区分时间保存吗

可以呀。每个时间步存一个txt，只要txt的命名有规则，比如

abc001.txt
abc002.txt
abc003.txt

paraview就能批量读进来。

如果不是这样规范命名，桌和代表时间的数字前面没补零，长度不一样，可以用rename之类的命令批量修改。

另外，你这时在linux里直接读写win的文件啊？我没怎么用过wsl，我记得很早的时候。在win上用wsl互相交换文件。要在linux里挂载win的盘符，从linux的根目录出发去找文件，比如：
/mnt/F

因为两个系统的“文件系统”不一样，一个是ext4，一个是NTFS。权限、命名规则不一样（win不接受一些特殊字符、不区分大小写）

用挂载的方式去读写文件，会有个不同文件系统转换的过程。

直接读写也能操作，但容易出问题。

bestucan

@Kirknwpu 怎么增大的计算域？整体放大还是拉长？

按说出现回流，加长流域大部分情况下能解决。出口那里本就该有个涡，被出口边界一分为二。加长，就可以躲开那个本就应该存在的涡。

要么是那里本不该有涡，但是因为精度或者还没收敛，不稳定。

试试先用压力出口稳一稳，阻碍一下回流。算到大概稳定，再切换到原本的出口边界条件。

bestucan

@dxl
开一个终端，用这个命令监视内存

watch -n 1 free -h

看看是不是“内存渐增，崩的时候满了”。用system monitor 看也行。

保存成视频应该有两种方法，一种是直接出视频，一种是连续图片。按说出连续图片的应该不怎么占内存。

不知道怎么存的坐标。按说应该是文本文件。这样的话就可以读取存储坐标的文本文件，然后用table To Points filter 转成坐标系的点。点的形状尺寸也能调。

bestucan

报错就是那个文件不存在

two_compiler_particle.0393.nut

paraview可以直接读openfoam的数据呀。

也可以试试更新以下paraview试试

bestucan

初始化的时候有几个选项，有混合、标准；相对值、绝对值。

可以试试查查不同

bestucan

https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1compressible_1_1humidityTemperatureCoupledMixedFvPatchScalarField.html

https://www.cfd-online.com/Forums/openfoam-solving/147129-condensation-using-openfoam.html

https://www.researchgate.net/publication/279643607_Modeling_of_Direct_Contact_Condensation_With_OpenFOAM

https://jhmtr.semnan.ac.ir/article_467.html

https://github.com/NimaSam/phaseChangeHeatFoam

bestucan

仿真结束的时候室温还是20度么？

有没有考虑温度越高，加热过的空气密度越低，对流会强烈。如果房间被加热，还会向外辐射。

实际中，房间向外辐射多少且不说。钢筋混凝土也会吸热。
比如说，背阴的房子“阴凉”，就是房体、地基（山里的阴凉就是山体很难暖热）温度低，进去很难把房间暖热。

仿真中的模型应该是理想绝热的。墙体、地基有没有设定材料，算比热容什么的。地基的吸热应该比墙体吸的更多。墙体还有个厚度，地面是无限厚的。应该设定个固定温度之类的。

可以检查下散热片附近的流场有没有加强，墙体、地面的温度升高多少。

bestucan

@bestucan

bestucan 发布的最佳帖子

bestucan 发布的最新帖子