能否通过直接设置扩散率为0来实现关闭扩散项?
理论上可以。但是要注意对流格式。容易不稳定。
射流雾化这个用VOF算,网格量不小。不太好debug。你这个问题,只能看出来是发散了。我看你最好先用不可压缩整收敛,然后再上弱可压缩。
你首页那个不是真正的报错。那个是intelmpi的log。你可以贴一下openfoam的log。不过类似这种问题,对于一个没有接手你们代码的任何人来说,都不容易debug,
一般情况下不会有差异。但是你这个不知道是不是就是本地就报错了。那肯定放云上也报错。
如果本地不报错,云上报错。那估计是并行的事情。比如本地3核不报错,超算300核报错,这也会有可能。
如果本地300核不报错,云上300核报错。这种很少见。但是总体,应该跟超算那面没什么关系。主要是openfoam
老铁你这啥都算啊,前几天还记得你弄IBM。为啥不用不可压缩呢。变量越多越容易出问题。这个错误就是发散了。
才看到这个回复。现在A100/V100太贵了,之前我问忘了是A100还是V100来着,要10万一个。好像国内有个A800稍微便宜一些。GPU计算这个如果硬件能够解决的话。目前来看GPU都是第三方写的外挂。还需要各大vendor做适配。也需要一段功夫。
不过你这整的挺厉害,都跑起来了,我对第三方写的包总是不敢尝试,经常出错。目前上亿网格我偏向petsc
我就是boxTurb添加的,然后mapField做的非均匀网格。你这个问题,看起来不像是初始场的问题,看起来像是压力梯度没设置好
setFields这种情况出现,一般就是那个box没写好,如果写好的话,0文件夹里面的U会出现内部场,而不是uniform,这个非常好debug
东岳老师今年是线下吗,北京还是杭州?上次在杭州上课受益匪浅,希望今年还是在杭州,直接从学校去听课,再跟着老师学习一下
下回大概率在北京。2020年之后,我在北京就没搞了。下回我会回国专门干这个事,同时把我博士国外导师叫过来,给大家线下做开场。然后我还会找几个助理,我在上面写代码,助理在下面处理各种编译不通的问题。
很激动。好久没线下讲课了!!
这个收敛特别慢,最后形成了正激波,choke了。我再减小出口压力的时候,也可以算出来,但是出口压力有一层异常值。这个我下周可能完事后再给你debug。现在没时间了。
“现在需要慢慢向你的设置靠拢”是啥意思呢
在一个正确的算例基础上,逐渐趋向你的设置
我直接跑了个你那个算例,确实这个样子。逐渐改了一下,正常了一点,现在需要慢慢向你得设置靠拢。就是一个算例算的时间太长了。下面这个图算一个要1000多秒。等我搞完了发你。
最近按你算例的codeStream代码加入了初始脉动速度
这个codeStream写的扰动,可以保证$\nabla\cdot\bfU=0$么
@MingfengWang 在 snappyHexMesh并行后导致对topoSet建立的faceZone流量监测不准! 中说:
就会出现对faceZone流量监测不准的情况
差多少。1%?5%?还是20 30%?
监测内部面的流量
怎么监控的
这个扩展磁盘不太好弄,经常出问题,因此我都默认搞成了100G。我也没有什么好办法可以简单的扩展硬盘。
如果开不开机,你可以再下一个虚拟机,用新虚拟机,挂上开不开机的虚拟机硬盘,单做一个硬盘用。
@Rachel0096 isoFoam那面使用的isoAdvector算法对alpha求解的。interFoam默认的MULES
@疏影横斜水清浅 矩方法在壁面上可以类似拉格朗日那种处理碰撞类似的速度。如果简单考虑颗粒流的话,一些小众的矩方法可以实现。你们要用矩方法算的话,这个坑挺大的,适合科研,不适合工程。
@疏影横斜水清浅 你说的确实有道理。
目前我能想到的只是在输运方程上植入重力导致的沉降。但是更细的,斜面分速度,还是颗粒速度,这个应该会影响结果,感觉拉格朗日可以算的很好。欧拉方法可能会确实细节
Airflow and Particle Transport in the Human Respiratory System
你可以看一下这个文章。可以处理沉积率。欧拉方法的。