李东岳 发布的帖子

能否通过直接设置扩散率为0来实现关闭扩散项？

理论上可以。但是要注意对流格式。容易不稳定。

射流雾化这个用VOF算，网格量不小。不太好debug。你这个问题，只能看出来是发散了。我看你最好先用不可压缩整收敛，然后再上弱可压缩。

你首页那个不是真正的报错。那个是intelmpi的log。你可以贴一下openfoam的log。不过类似这种问题，对于一个没有接手你们代码的任何人来说，都不容易debug，

一般情况下不会有差异。但是你这个不知道是不是就是本地就报错了。那肯定放云上也报错。

如果本地不报错，云上报错。那估计是并行的事情。比如本地3核不报错，超算300核报错，这也会有可能。

如果本地300核不报错，云上300核报错。这种很少见。但是总体，应该跟超算那面没什么关系。主要是openfoam

很奇怪的问题。
如果你可以提交一下算例，我可以给你看一下（网格数量少的话）

你可以在本地尝试，本地跑通了，再弄到集群上测试。

老铁你这啥都算啊，前几天还记得你弄IBM。为啥不用不可压缩呢。变量越多越容易出问题。这个错误就是发散了。

才看到这个回复。现在A100/V100太贵了，之前我问忘了是A100还是V100来着，要10万一个。好像国内有个A800稍微便宜一些。GPU计算这个如果硬件能够解决的话。目前来看GPU都是第三方写的外挂。还需要各大vendor做适配。也需要一段功夫。

不过你这整的挺厉害，都跑起来了，我对第三方写的包总是不敢尝试，经常出错。目前上亿网格我偏向petsc

这个是rhoCentralFoam算的，我加密了一下网格，看起来可以，分别是速度跟压力

我就是boxTurb添加的，然后mapField做的非均匀网格。你这个问题，看起来不像是初始场的问题，看起来像是压力梯度没设置好

这玩意不太好整...也不能给出啥建设性意见

setFields这种情况出现，一般就是那个box没写好，如果写好的话，0文件夹里面的U会出现内部场，而不是uniform，这个非常好debug

@upc_ngh 还来啊，那多谢支持！北京见！

东岳老师今年是线下吗，北京还是杭州？上次在杭州上课受益匪浅，希望今年还是在杭州，直接从学校去听课，再跟着老师学习一下

下回大概率在北京。2020年之后，我在北京就没搞了。下回我会回国专门干这个事，同时把我博士国外导师叫过来，给大家线下做开场。然后我还会找几个助理，我在上面写代码，助理在下面处理各种编译不通的问题。

Make some noise!

很激动。好久没线下讲课了！！

这是我用rhoPimpleFoam跑的，跟你跑的fluent的差不多

CDV2.tar.xz

用0文件分网格出来，和0.0005分网格出来，确实后者会少一个“point”文件

你把0.0005改成0试一下

在楼上

laval.tar.xz

这个收敛特别慢，最后形成了正激波，choke了。我再减小出口压力的时候，也可以算出来，但是出口压力有一层异常值。这个我下周可能完事后再给你debug。现在没时间了。

“现在需要慢慢向你的设置靠拢”是啥意思呢

在一个正确的算例基础上，逐渐趋向你的设置

codedfixedvalue怎么写的

decomposePar的log怎么提示？

我直接跑了个你那个算例，确实这个样子。逐渐改了一下，正常了一点，现在需要慢慢向你得设置靠拢。就是一个算例算的时间太长了。下面这个图算一个要1000多秒。等我搞完了发你。

可压缩那面就是比不可压缩容易发散。你把算例发上来我给你调一下。顺便当做oks课程算例。

看起来猜测边界条件搭配设置的问题

U的outletVelocity人为指定的么

@杨英狄 这个跟interFoam非常像，就没写，可以等我课程之后写一下

可能粒径算的不对，导致流场没对

flowRate采用哪种方式监控的？

@SHUKK 这个里面有算例，你直接下载，allrun就可以

https://www.cfd-china.com/topic/6231/中性大气环境湍流动能的自保持-附有算例下载/24

我那个算例好像没有添加symmetry类型。这个是你自己加的么。

最近按你算例的codeStream代码加入了初始脉动速度

这个codeStream写的扰动，可以保证$\nabla\cdot\bfU=0$么

颗粒是质点。如果你要看球那种，表面温度高，里面温度低那种的，应该工作量不小。

@aiweimo 也行，有链接不

@MingfengWang 在 snappyHexMesh并行后导致对topoSet建立的faceZone流量监测不准！ 中说：

就会出现对faceZone流量监测不准的情况

差多少。1%？5%？还是20 30%？

监测内部面的流量

怎么监控的

这个扩展磁盘不太好弄，经常出问题，因此我都默认搞成了100G。我也没有什么好办法可以简单的扩展硬盘。

如果开不开机，你可以再下一个虚拟机，用新虚拟机，挂上开不开机的虚拟机硬盘，单做一个硬盘用。

是yf么？不应该啊。你现在应该openfoam玩的很溜啊。

还是说你同学借了个号发帖 :-)

可以看一下compressibleInterIsoFoam

• 你把tgradU.clear();去掉，
• 然后使用的时候用tgradU()()

前面的tgradU怎么用的，把涉及到tgradU的代码贴一下

tgradU是不是tmp，前面怎么用的

你最好debug到具体某一行（比如这一行去掉就没问题，这一行放进去就有问题）

@Rachel0096 isoFoam那面使用的isoAdvector算法对alpha求解的。interFoam默认的MULES

这个问题是某个tmp类型，已经销毁掉了，但是又被继续调用了。tmp类型使用的使用可以后面加()试一下

@aiweimo 不知道后来这个文章怎么样了？

这个网格按道理来说 好像是负体积啊

只要不到大气以上那种太大的尺度，多大都行。

@疏影横斜水清浅 计算域多大都行。为什么会问这个，CFD这面计算域大，网格大了就行了啊？

@疏影横斜水清浅 矩方法在壁面上可以类似拉格朗日那种处理碰撞类似的速度。如果简单考虑颗粒流的话，一些小众的矩方法可以实现。你们要用矩方法算的话，这个坑挺大的，适合科研，不适合工程。

@疏影横斜水清浅 你说的确实有道理。

目前我能想到的只是在输运方程上植入重力导致的沉降。但是更细的，斜面分速度，还是颗粒速度，这个应该会影响结果，感觉拉格朗日可以算的很好。欧拉方法可能会确实细节

Airflow and Particle Transport in the Human Respiratory System

你可以看一下这个文章。可以处理沉积率。欧拉方法的。

我还不知道VMware vSphere 我这面服务器跑的都是普通的vmware 我让学生研究下