CFD中文网

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

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

大概就是，可以软件渲染（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上，
Screenshot from 2022-06-11 15-52-46.png
就是这个方法：
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 的

楼主，请问您解决这个问题了么，我最近做重叠网格入水也遇到这个问题，求助

@bestucan 真的非常感谢您的建议。对于二维的case尝试了一下，我设置一个不变形的界面，进入filter区域的粒子数量大约是等同于流通量。但我手上的计算资源和代码的优化方面并不乐观，耗时久。如果界面变成可变形的情况，计算时间应该会加倍，时间上等不起。

不好意思 数据在这里
0_1542843045909_48db69f6-3dde-40bc-9b8d-70c1b51caa28-image.png Dd [um] Number
0.00000220 17200.00000
0.00000439 71000.00000
0.00000659 97900.00000
0.00000878 84400.00000
0.00001100 53900.00000
0.00001320 35000.00000
0.00001540 23900.00000
0.00001760 15300.00000
0.00001980 12200.00000
0.00002200 9060.00000
0.00002410 6550.00000
0.00002630 5540.00000
0.00002850 4090.00000
0.00003070 3440.00000
0.00003290 2530.00000
0.00003510 2210.00000
0.00003730 1960.00000
0.00003950 1490.00000
0.00004170 1240.00000
0.00004390 938.00000
0.00004610 860.00000
0.00004830 743.00000
0.00005050 652.00000
0.00005270 550.00000
0.00005490 427.00000
0.00005710 339.00000
0.00005930 282.00000
0.00006150 268.00000
0.00006370 194.00000
0.00006590 208.00000
0.00006810 171.00000
0.00007030 131.00000
0.00007240 106.00000
0.00007460 125.00000
0.00007680 76.00000
0.00007900 66.10000
0.00008120 52.80000
0.00008340 63.80000
0.00008560 23.50000
0.00008780 42.40000
0.00009000 34.10000
0.00009220 33.50000
0.00009440 18.30000
0.00009660 21.60000
0.00009880 17.70000
0.00010100 14.00000
0.00010300 10.40000
0.00010500 10.30000
0.00010800 10.20000
0.00011000 6.71000
0.00011200 0.00000
0.00011400 3.29000
0.00011600 6.51000
0.00011900 0.00000
0.00012100 0.00000
0.00012300 0.00000
0.00012500 0.00000
0.00012700 0.00000
0.00013000 0.00000
0.00013200 0.00000
0.00013400 0.00000
0.00013600 2.98000
0.00013800 2.96000
0.00014100 0.00000
0.00014300 0.00000
0.00014500 0.00000
0.00014700 0.00000
0.00014900 0.00000
0.00015100 0.00000
0.00015400 0.00000
0.00015600 0.00000
0.00015800 0.00000
0.00016000 0.00000
0.00016200 0.00000
0.00016500 0.00000
0.00016700 0.00000
0.00016900 0.00000
0.00017100 0.00000
0.00017300 0.00000
0.00017600 0.00000
0.00017800 0.00000
0.00018000 0.00000
0.00018200 0.00000
0.00018400 0.00000
0.00018700 0.00000
0.00018900 0.00000
0.00019100 0.00000
0.00019300 0.00000
0.00019500 0.00000
0.00019800 0.00000
0.00020000 0.00000
0.00020200 0.00000
0.00020400 0.00000
0.00020600 0.00000
0.00020900 0.00000
0.00021100 0.00000
0.00021300 0.00000
0.00021500 0.00000
0.00021700 0.00000
0.00022000 0.00000
0.00022200 0.00000
0.00022400 0.00000
0.00022600 0.00000
0.00022800 0.00000
0.00023100 0.00000
0.00023300 0.00000
0.00023500 0.00000
0.00023700 0.00000
0.00023900 0.00000
0.00024100 0.00000
0.00024400 0.00000
0.00024600 0.00000
0.00024800 0.00000
0.00025000 0.00000
0.00025200 0.00000
0.00025500 0.00000
0.00025700 0.00000
0.00025900 0.00000
0.00026100 0.00000
0.00026300 0.00000
0.00026600 0.00000
0.00026800 0.00000
0.00027000 0.00000
0.00027200 0.00000
0.00027400 0.00000
0.00027700 0.00000
0.00027900 0.00000
0.00028100 0.00000
0.00028300 1.84000
0.00028500 0.00000
0.00028800 0.00000
0.00029000 0.00000
0.00029200 0.00000
0.00029400 0.00000
0.00029600 0.00000
0.00029900 0.00000
0.00030100 0.00000
0.00030300 0.00000
0.00030500 0.00000
0.00030700 0.00000
0.00031000 0.00000
0.00031200 0.00000
0.00031400 0.00000
0.00031600 0.00000
0.00031800 0.00000
0.00032100 0.00000
0.00032300 0.00000
0.00032500 0.00000
0.00032700 0.00000
0.00032900 0.00000
0.00033200 0.00000
0.00033400 0.00000
0.00033600 0.00000
0.00033800 0.00000
0.00034000 0.00000
0.00034200 0.00000

@liuxiaodong les计算圆管管道流动的，尝试了很多都最后变成层流了，这是为什么呢

谢谢版主。

从代码来看，e表示回弹系数，mu表示一种损耗。e越大，回弹越完全，mu越大，速度被消耗的越大。这些值影响你颗粒撞击壁面的回弹行为。
你可以手动计算一下，看看U有什么变化

case PatchInteractionModel<CloudType>::itRebound: { keepParticle = true; active = true; vector nw; vector Up; this->owner().patchData(p, pp, nw, Up); // Calculate motion relative to patch velocity U -= Up; scalar Un = U & nw; vector Ut = U - Un*nw; if (Un > 0) { U -= (1.0 + e_)*Un*nw; } U -= mu_*Ut; // Return velocity to global space U += Up; break; }

fluent自定义过饱和函数时应该用那个压力数据来计算？我的静压只有不到100pa，太小了

这是用离散法求解PBE的增长过程。如果没有增长过程，如果不是离散法，就不需要这么做。a是波速，Minmod是limiter。如果Minmod目前你还不熟悉，我大体知道你的CFD基础在哪，你需要补CFD基本理论，再来看PBE:138:

@东岳 总压减小是因为速度变小，还是静压变小，或二者都有？无粘不可压的时候就是动能与压力势能的转化，是守恒的。

请问在OpenFOAM中，如何定义一个轴对称模型来仿真Taylor Couette flow？如何设定wedge边界？谢谢大神！

@leo_chou 大批量的工业算例（本身可能不复杂，但就是工况多），LBM是不是也需要大规模网格和较长时间？工业应用CFD对时间的要求基本上是小时级别的。

另外好奇问下，复杂几何下的LBM如何，需要湍流的LBM呢？我觉得大部分案例，复杂几何+湍流模型是标配，没有湍流的话可能有多相或者化学反应啥的，这些LBM是不是会更强一些？

图片
我最近在使用reactingParcelFilmFoam计算水滴与结构表面撞击并形成水膜的过程中，我想往模型中添加一个水滴的收集系数，这个系数需要知道水滴在撞击到结构时的位置坐标和初始的位置坐标，请问有人知道该怎么确定这两个坐标吗？

怎么能够对计算完的case设定单元区域，取得这个区域内自己想要的值

我在用interFoam时，笛卡尔网格可以得到合理的计算结果。但换成非结构网格（prism之类）后，计算虽没有发散，但界面有问题，请问有什么需要改进的地方吗？@东岳 @tidedrinker

sed -e 's|^mirrorlist=|#mirrorlist=|g' \ -e 's|^#baseurl=http://mirror.centos.org|baseurl=http://mirrors.tuna.tsinghua.edu.cn|g' \ -i.bak \ /etc/yum.repos.d/CentOS-*.repo yum clean all && yum makecache

@东岳 i其实是第一个网格，n是最后一个网格
也就是
先求每一个网格中的alpha * x 再对整个计算域求和
然后求和所有计算域中alpha的值
两者做比值

缩小模型只能保证一个无量纲数相同，一般的流动考虑雷诺数就行，缩放会带来一些非真实的物理效应，但是计算量不够就不要讲究了