PISO算法中的对流项fvm::div(phi,U)是如何迭代求解的呢？

@李东岳 李老师，对于周期边界下的压力场，piso循环中有pEqn.setReference(pRefCell,pRefValue);用来确定压力场。那在计算周期边界下的被动标量时，应该用什么来对温度场进行修正呢？
:135:

ff12dfb1-d78b-4f7d-812d-db36e6749698-图片.png
请问各位老师，根据参差图 --Initial residual , Final residual
这个能判定收敛吗？
初始残差没到设定的值，且计算到后边开始波动，红色方框所示；
Final residual虽然可以到设定值，但是在一部分计算中波动很大，红色圆圈所示。
请问是什么原因呢。
期待大佬回答解惑。

@五好青年 可以，就是easycae的网盘，你注册个用户应该就有，

@李东岳 哈哈哈，物理隔绝了:mihu:

图一梯度算子的分量$\frac{\p }{\p x}$，不是某个物理量的梯度分量

$\frac{\p }{\p x}$不是一个函数么，你要写一个函数？

@李东岳 收到，谢谢李老师。

液滴尺寸比网格大，蒸发源项(组分、能量)加到网格里面会不会出问题，有没有老师知道OpenFoam这方面是怎么处理这个问题的

@李东岳 谢谢李老师。

网格流动性是什么..

在WallSpringSliderDashpot.C中，有这样几句：

scalar etaN = alpha_*sqrt(p.mass()*kN)*pow025(normalOverlapMag); vector fN_PW = rHat_PW *(kN*pow(normalOverlapMag, b_) - etaN*(U_PW & rHat_PW)); ... scalar etaT = etaN; ... (if条件：颗粒不滑移的条件下) fT_PW = - kT*tangentialOverlap_PW - etaT*USlip_PW;

这里的颗粒碰撞中的切向阻尼系数，etaN，等于法向值，然后分别应用到求碰撞力的切向和法向分量，fN, fT，中。是不是意味着，我在颗粒配置文件中通过alpha设置的回弹系数值e, 就是同时设定了切向和法向的回弹系数，eT, eN, 都等于e？

请问这个理解对不对？

最近在做一个污染物浓度扩散的算例，需要选定一条线并将其设置为污染物浓度扩散的线源，但在setFields中好像没有选定线的设定，所以哪位大佬要是知道请告知，万分感谢。

你说的把Q的插值变成比值，应该是Ω方法吧？？？
Liutex-Ω的话好像比较麻烦

@西湖冷月 你具体编程过吗？

@yfleung @荷兰后腰阿贾克斯 @波流力
大佬朋友，请问你们处理400万以上网格，速度很慢，无法出结果的问题，解决了吗？
感谢

@队长别开枪 明白了，谢谢老师！

@oitocfd 好的，我尝试tecplot打开了其他格式是binary的算例，发现可以打开，那么就说明是之前的算例的问题，我再修改试试。谢谢您。

@星星星星晴 谢谢，我再看看

@西极若木
这个算例中在Stokes2里的angle是不是可以设置波浪来的角度？也就是模拟斜浪

@tens 谢谢您的热心指点，虽然还有点晕，我后面仔细体会体会:quwan: ，再次感谢:xinxin:

@sungda 这是二维的自适应网格吗？是怎么实现的