Posts made by bestucan

@李东岳 最新版的 paraview 已经支持 “不重组数据” 就 ”可以查看了“

估计李老师看成“最新的已经不支持了“

@hy1112006 看不懂

哦，对了，最新版的 paraview 已经支持不重组数据就可以查看了。可以到官网下载

加载好OpenFOAM环境后，运行下面两句，

sol
find . -name 'createFields.H' |xargs grep 'NO_READ' .

然后，总有一个能抄对...

另外

which mysolverFoam | xargs ls -la 

核验时间戳是不是刚编译的

信息量过于少没法判断，看我签名里的链接

先用串行，不开多线程试试。如果不行，可能是导出的格式有问题。如果行，那就是并行换镜没配置好

这个属于后处理了吧。可以在paraview里画，plot data over time

报错可能没复制完，可以在你运行的命令后面加个导出log的尾巴

reconstructPar |& tee reconst.log

然后把这个log上传上来

不阔气不阔气，我也不是很懂，我猜的

@forest-wang 我们这讨论有点对不上呀
我说可能是网格导致的扩散形成了假粘度。
你把问题又更详细的描述了一遍。

你可以测试一下是否是“网格-扩散-粘度”造成的，或者明确不是这个造成的。

气泡上浮速度明显变慢了

这和液体中粘度越高，气泡越慢是通的。

我做的其他两种形状的气泡都没有出现这种问题

就是第一个帖子中的另外两种？那M1和M2也有一段不重合呀。兴许M3和他们不一样的原因找到了，那段不重合的原因也找到了。

三个可以大概有这样的结论：网格越密，雷诺数越有可能小。

现在得不到明确的结论是三个之间缺少个”闭环“。

如果M3里的默认网格密度是M1，自适应网格密度是M2。然后三个实验两两比较会有很多结论。然后增减网格密度验证结论。

现在M3里多出了一种网格密度（自适应的那个），比较难比了

雷诺数的意义是惯性力和黏性力的比值。

考虑可能网格的划分影响了黏度的计算。不是有个数值扩散么，数值扩散相当于人工粘度了，数值扩散好像是和网格划分相关的

@vbcwl 是的。那么问题就变成“of6可能因为openmpi无法并行了” 和无法切换版本无关了。

按说按照官方指导应该用的是源里的openmpi，不是不自己编译的。

如果确实用了自己编译的openmpi，那么也可以多版本openmpi共存。这是我的配置文件，可以参考修改成你的

#--------------------------------------- openfoam
alias openmpi_1.10.2='export PATH="$HOME/.local/share/OpenFOAM/ThirdParty-4.x/openmpi_1.10.2/bin:$PATH" \
export LD_LIBRARY_PATH="$HOME/.local/share/OpenFOAM/ThirdParty-4.x/openmpi_1.10.2/lib:$LD_LIBRARY_PATH"'

alias openmpi_1.6.5='export PATH="$HOME/.local/share/OpenFOAM/ThirdParty-2.3.x/openmpi_1.6.5/bin:$PATH" \
export LD_LIBRARY_PATH="$HOME/.local/share/OpenFOAM/ThirdParty-2.3.x/openmpi_1.6.5/lib:$LD_LIBRARY_PATH"'

alias of23x='source $HOME/.local/share/OpenFOAM/OpenFOAM-2.3.x/etc/bashrc; openmpi_1.6.5; export WM_CC='gcc-4.6' export WM_CXX='g++-4.6''
alias of4x='source $HOME/.local/share/OpenFOAM/OpenFOAM-4.x/etc/bashrc; openmpi_1.10.2;'

openmpi的多版本加载得在编译OpenFOAM之前。

openmpi编译安装的步骤

@vbcwl 我说的是就在你的终端输入这四个字母：f o a m
然后回车键。

你就会跳转到你用的版本的OpenFOAM的文件夹了，你就知道你的版本环境了。

至于求解器运行不起来，再找原因，但先把“是不是因为版本环境没切换导致的”这个原因排除掉

@vbcwl 那怎么定义没效果呢？

切换完环境后运行
foam
就知道有没有效果了

@vbcwl 按你的路径写，我们安装的位置不一样。
写了alias之后就不用些source ...开头的了

alias就是给source ...起个别名。运行of6就相当于运行source ...

如果你要使用6版本的 OpenFOAM，就在终端运行命令：
of6

换版本就of6h者of9

@strive_dut

通过将我的求解器求出的压力和温度传递给函数FXc(float temperature, pressure)，然后在函数中可以得到具体的Xc或者Yc值

是的。

因为求解过程中每个网格单元的温度和压力都是在不断变化的，那我的拟合函数应该是一个曲面的数据吧？

是的，就是你第一个帖子里的那两个曲面（其实就是用函数表达那两个曲面图）

往.bashrc里添加类似于这些的东西

alias of6='source $HOME/.local/share/OpenFOAM/OpenFOAM-6/etc/bashrc;'
alias of7='source $HOME/.local/share/OpenFOAM/OpenFOAM-7/etc/bashrc;'
alias of8='source $HOME/.local/share/OpenFOAM/OpenFOAM-8/etc/bashrc;'
alias of9='source $HOME/.local/share/OpenFOAM/OpenFOAM-9/etc/bashrc; source $HOME/.local/share/OpenFOAM/blastfoam/etc/bashrc'

然后用of6、of7之类的命令加载想要的版本的环境

@strive_dut 不是识别呀，就是，你给出的流程图中，读取完$X_c$之后就开始$Z_c$与$X_c$做比较了。哪里用到$X_c$，就把$X_c$替换成其拟合出的函数。

if (Zc >= Xc) {
...
} else {
...
}

换成

if (Zc >= FXc(temperature, pressure)) {
...
} else {
...
}

温度和压力要取该网格的？还是该区域的？看具体情况。

至于这个FXc()函数怎么定义

float FXc(float temperature, pressure) {
result = ...// matlab拟合出来的函数，如何用温度和压力计算质量分数

return result;
}

@sea 没有解决

但是可以保存当前设置然后读取：
file -> save state

开paraview之后进行的操作都可以保存住

@strive_dut 关于那两张图，如果没有原始数据，只有数据点。那得选差值方法，可以看看这个 。

这样插值了之后再读入OpenFOAM（按我第一回帖说的方法），就是查表法获得数值了。如果是在OpenFOAM里插值，没做过，估计是C++手动实现那些插值算法。比较麻烦。

那两张图可能就是用实验测得的数据点拟合出来的（因为比较光滑）。

也可以用现有的数据点，把那两张图拟合成函数表达式。可以看看这个。

这样把求解器里的 $X_c$出现的位置换成拟合出的函数就行了： FX_c(temperature, pressure)

@strive_dut 你问的第一个问题是怎么把外部数据导入OpenFOAM参与计算：

利用 C++程序语言将正庚烷-氮气混合物气液平衡数据库转化为 OpenFoam 可以调用和读取的脉谱数据

我说的可以解决这个问题。

你的第二个问题是求解过程中如何使用脉谱数据：

每执行一个时间步，我都需要调用这个脉谱数据来对比计算域的每个单元的数据，然后对于需要部分单元的数据采用两相参数进行重新求解

这个是关于如何 选用求解器 或者 修改求解器 的问题了。选用求解器得看需求。

另外，我又看了看你的图，大概想象到，你想要在计算过程中随时根据当前点的温度和压力得到那个质量分数？

如果是这样，可以把那个面拟合成个函数写到求解器里需要用的位置。
或者把散点数据输入OpenFOAM，使用时再插值，得到任意一点的值。

不知道哪个更快

@臭臭爱趴趴 如果操作比较复杂，可以搜 cfd post 如何进行相关操作。
cfd post 也没有，就导出成ASCII，然后用 matlab

@strive_DUT 随便找个算例，计算一两个时间不。时间节点文件夹里有温度（标量）、速度（矢量）的文件，看看那些文件的格式。然后无论用什么语言，把数据按那个格式输出

view -> light inspector

然后在右边的面板把 light Kit 的对勾取消

@李东岳 摩尔定律之东岳流体分律

你看你图中圈着的上一个Phase-sum volume fractionmin 是正的，圈着的地方变成负的了，当然max那上下都负的。max先数值溢出，min跟着溢出。

说明这是计算本身的不稳定。

速度快了，能够不发散，是把不稳定多出来的那部分赶紧吹出去了；慢了，吹不出去就原地爆炸了。

我看了格式，里面有不少迎风格式，好像迎风格式不能考虑回流？
这种管子冲快了，刚出狭管上下两个角落应该会有回流涡吧？
可能因为这个

我对大型计算缺乏想象力
要是往气象、核爆那边想，好像占多少内存都理直气壮

要是普通的，那可能泄漏，或者后“处理方法”本身比较吃内存

可以试试扩大swap分区。
如果装系统的时候没有用LVM，后面想改分区大小就比较麻烦了。

可以把 swap 当成虚拟内存，可以用

free -h

查看

@李东岳 佛系计算，不给服务器压力，不push

这个问题看起来像机器翻译的。。。
请问你是中国大陆的北方网友么

流体比较细？我理解为颗粒流中颗粒比较小？

网格比流体颗粒大，每个网格存在变量：相分数（包含流体颗粒的量）
网格比颗粒流体小，被流体占据的网格堆一起，他们的边界是相界面。

这是多相流的两种方案

工作中常做好备份工作，否则很容易：添加新功能成没成功且不说，旧功能无法实现了。

如果有备份拉出来，用 diff 对比不同就很容易排查了

是否可以通过一个时间步后粒子的坐标相对于界面的位置来判断粒子是否穿过界面

https://www.baidu.com/s?wd=判断点在线的某侧

给粒子在界面某侧声明个变量和一个flag，每次修改变量的时候给flag扶植1，然后每次计算完统计flag为1的粒子数量，同时给flag赋值0.

然后看每个时间步统计的粒子数量差就知道通量了

搜 openQBMM

@李东岳 我就记得我之前搞过，后来坑太深了，就没搞了。

这些涉及到图形处理的各个层，又是平台不同，又是商业公司夹私货，那些个组件的名字本来就乱七八糟，缩写规则还不一样。想短时间了解就更乱。

大概就是，可以软件渲染（OSMesa），可以硬件渲染（EGL）。
其实都是基于Open GL，就这三个英语词空格隔开一起搜都一大堆讲更多术语关系的文章。

paraview 默认用软件渲染，想要硬件渲染把 EGL 编译选项打开就行了，这就有个忘了打开报错的：
https://discourse.paraview.org/t/build-paraview-for-nvidia-gpu/8309
但这个是单GPU的，不知道多GPU的怎么搞，

我的GPU是当显卡用的，一开paraview就调用GPU，当然是为“图形界面渲染"的,不是为“paraview渲染”的，这俩还都是opengl下面的，那些名词太乱了。

我试着编了一下 pvserver能多线程跑在GPU上，

就是这个方法：
https://github.com/Kitware/ParaView/blob/master/Documentation/dev/build.md

当然编译之前已经安装好nvidia自带的驱动了。

nvidia 早些年自己出过 paraview 的插件：
https://www.nvidia.com/en-us/data-center/index-paraview-plugin/
不过是针对 paraview 5 的

@李东岳 我电脑上有个GPU，有空我试试

@李东岳 下载的paraview里有个 bin ，bin里有 pvserver

mpirun -np 4 pvserver

然后会显示pvserver的ip（或者localhost）和端口（-p 1000 可以修改端口）。这个开pvserver的终端不能关。

pvserver --help能看到更多选项。

然后开 paraview 窗口，file->connect->add server。把上面的ip和端口填上。

pvserver有挺多模式，什么数据服务、渲染服务，共同服务，加参数还能指定GPU什么的，但是接口老变，得现用现搜了。

要是编译的paraview，得在编译前有个修改编译选项，没搞过。

paraview 可以并行，可以用GPU，但是需要些麻烦的操作（开并行服务然后链接或者修改paraview编译选项打开并行）。不这样，靠单核上限有些贵且上限低。

启动的时候在启动命令后面加上后缀试试

-driver x11

不对，你这是 wsl 啊，那不知道了，wsl 的图形界面比较麻烦，wsl2 相对于1改动还挺大，不清楚怎么搞出来图形界面。

可以试试

$ sudo apt install xclock
$ xclock

如果这个都打不开，那是系统的图形界面没配置好

按说apt安装的肯定没问题，你是不是照着编译paraview的教程添加了什么 export QT_SELECT=5 之类的东西？

实在不行
https://www.paraview.org/download/
这里下载，解压，进入 bin 目录

./paraview

paraview 的操作可以写 python 脚本。直接写比较麻烦，可以宏录制之后复制粘贴改参数

tools -> start trace

不知道

如果不是私有代码的话，可以在官网的手册里搜，右上角可以搜：

https://www.openfoam.com/documentation/guides/v2112/api/classFoam_1_1phasePair.html#ad38615aa705d601149b6ac5d346ae4f2

@zy-gong

看多相流分类：

https://cfd-china.com/topic/1152/多相流与矩方法的数学模型

找函数出现的位置（终端路径切换到 OpenFOAM 文件夹）：

grep -r 'dispersed()' .

@vbcwl 压力驱动一般都有功率吧？比如 P=FV，
改了g变成恒力驱动了，这种驱动效果，是改整体的势能。

就像，自由落体不能理解为压力驱动吧，我猜是这样的

@李东岳 国内计算机系的用 linux 的都没有多少。相关软件生产商因为生态和“在固定规则下博弈”的能力不足，也天然喜 win 厌 linux，搞些改不掉的默认，卸不掉的组件。都穷怕了。

好在现在国家采购在硬推这些。

mac 至少在普及，mac 上的生态对 linux 也有助力（appimage 格式的软件），至少打破了电脑一定是 win 的固有印象。
虽然很多人买了 mac，在 mac 上装 win。。。

把 return 挨个屏蔽，看看屏蔽了哪个报错。。。
就知道调用哪个了

我想知道它是什么意思呢？

好像是用 u 和 dim 初始化了一个新的 fvVectorMatrix 类型的变量

这个截图太错位了，可以上传 log看看

icoFoam > run.log

很多问题就是：
试试就知道了
代码不被卡脖子，相对于芯片的根本性优势就是：试错成本低

@ljxneu军祥 是的