在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 这是二维的自适应网格吗？是怎么实现的

单核计算：
ff8a928a-d108-454d-8548-a662114c6a34-image.png
并行计算：
2674cfa8-2f86-410a-9ebc-144e98f89800-image.png

@mania 你用RAS计算得到U、R、L和points，直接放到主程序计算就可以。of2006有自带的算例

@tens 好的谢谢:146: 我试一试

@jqbu 基于OpenFOAM的甲烷湍流燃烧LES数值研究_龙文增;不知道这篇文章能不能帮助到你

最近几天调程序老是对不上，今天才发现只有标准的结构性网格才是这样，如果是非结构网格，或者SHM生成的网格在过渡层处这个数量不是固定的
还是要细心点:136:
d0883aaa-8a7f-4939-9a9f-cf30552dfc9a-image.png

额 好像发现个问题，质量守恒方程应该是这个
7973d26e-efef-4af5-9e55-686a447c4996-image.png

其中𝑆(𝜌)与 𝑆(𝛾)形式一样，只不过自变量由相体积分数变成了密度，不知道这样对不对？

感觉有点弄乱了，求大佬指点，不胜感激:135:

@小考拉 自问自答来解决这个问题。

首先定义一个全为0的流场命名为c_flow

然后利用下面的命令

c_flow.boundaryField()[patchID] = c_flow_origin.boundaryField()[patchID]

最后输出c_flow流场即可:chitang:

@鲸落 请问你解决了吗~动态读取cellZone~

众所周知，当我们在进行并行计算的时候，需要把网格进行剖分，OF中可以使用第三方库进行区域分解，可是我们知道，OF中的区域分解是不带重叠层的（或者说没有halo单元），也就是，同一个网格仅可能属于一个处理器。
那么问题来了，我现在想实现带有halo层的重叠型的区域分解，我该怎么实现呢？或者说，有没有那位大佬做过类似的工作或者相应的资料吗？跪谢了！

@hongjiewang 在 温度边界条件 中说：

1632279639(1).png
气液界面是不断变化的，请问红色方框内的边界条件如何添加~

可以识别出气液界面，但是想知道如何在同一个边界里使用第二类和第三类边界条件