@李东岳 感谢感谢，写得好！

@李东岳 谢谢东岳老师，我去研究下拉格朗日看下能不能满足我的要求

可以瞅瞅 blastFoam. 直接点爆炸:jingya:

@小白好海牛 https://zhuanlan.zhihu.com/p/315686722，量纲只要能计算本身并没有实际意义，你看看这个，量纲好像就是和Gammac（扩散系数）一起设置的，不是真正的浓度单位

打扰了……困迷糊了，结构错误，已解决:136:

https://www.cfd-china.com/topic/4790 看文章最后

@zhangxc0223 请问您使用的什么模型来做的粒子湍流扩散，边界条件是如何设置的？求教

单节点4核47秒，双节点每个2核72秒？

@vbcwl 老哥我这么理解对不对，就是速度设置成U+，也就是U/Ut，那这个按照自带的算例Ret=395的话，是不是就是U=Ubar/ut=0.1335/0.079呢

@尚善若水 感谢回复！无反射边界也在一些文献里看到过，想问不可压缩流体也可以使用无反射边界吗？

@李东岳 我在水槽底部5公分范围内加了颗粒，然后开始计算，水槽内是静水没有流速。开始颗粒正常沉降，最大相分数也随之增大，但很快就超过0.63了，之后颗粒的动量值还在一直增大，才算了0.2秒的过程就崩了:136: ，怎么都已经达到packing limit了颗粒动量和最大相分数还在持续增大，而且大得离谱，而且我流场是静态的啊，我不清楚是哪儿出的问题:143:
搜狗截图20220721121946.png

一开始的time step continuity errors就很大

稳态是可能得。瞬态要保证很小。

时间步长连续性误差会随步长增长，

这个看起来设置不太对。可以先用CPU跑没问题了之后，再跑GPU

@星星星星晴 org这面一直在整合。粒子那面自行will加入之后就变动很大。org那面湍流模型、热物理模型也是变动特别大。主要是架构方面的。

@yu_tian 这个是多计算域求解器里面的互相影响的两个边界条件

请问各位大佬，是否有phaseProperties设置教程呢？

of9中多相欧拉的算例，在constant/phaseProperties 中都要设置一个type；
已知的可选设置有:
basicMultiphaseSystem;
interfaceCompositionPhaseChangeMultiphaseSystem;
interfaceCompositionPhaseChangePopulationBalanceMultiphaseSystem;
populationBalanceMultiphaseSystem;
thermalPhaseChangeMultiphaseSystem;
thermalPhaseChangePopulationBalanceMultiphaseSystem.

目前比较关心的是什么时候使用interfaceCompositionPhaseChangeMultiphaseSystem，什么时候使用thermalPhaseChangeMultiphaseSystem呢？这两者具体有什么不同呢？

还有就是blending相里，算例有时会设置传热和曳力的离散/连续分数设置，有时又不会，这是与传热和曳力的模型有关吗
万分感谢:xiexie:

@xpqiu 好的，已经尝试成功，感谢老师指点啊！

79c9f160-dc4a-4ed1-8d71-3b9cf501958d-image.png
各位带佬上午好！:chouchou: :chouchou:
小弟最近在做一个二维圆柱绕流的算例，现在想绕着二维圆柱一圈提取数据，就是提取圆柱表面360°的数据，如上图的X轴所示，纵坐标提取U和P啥的。有带佬知道这该咋样实现嘛？:chitang:

请教，目前计算的是一个三维域，视线方向指的是Z方向。OpenFOAM如何实现，对某一标量沿着Z方向的积分计算？

如果颗粒刚度比较大的话，可以用重叠网格（overset mesh），也就是背景网格+颗粒周边网格，两套网格重合的地方会进行插值处理。这样可以回避网格重画的难点，不过这东西算二维还行，三维想算快算好比较折磨。

各位大佬大家好，我最近准备在rhoSimpleFoam的基础上加入自己的算法，算法中比较关键的点是有另一组速度、压力和温度场（暂且叫它们Uf，pf，Tf）。

之前在simpleFoam中顺利地植入了，但是在可压缩这里我发现：thermo.rho()这些函数在调用时，读取的物理量是求解器原来的U、p，如果我再另外定义Uf，pf的话，好像不能直接使用thermo.rho()函数进行另一组速度压力场的密度修正？

附上一段src/thermophysicalModels/basic/fluidThermo/fluidThermo.C中的代码，麻烦大佬们看下是不是这样的

Foam::fluidThermo::implementation::implementation ( const fvMesh& mesh, const word& phaseName ) : p_(lookupOrConstruct(mesh, "p")) {} Foam::fluidThermo::implementation::implementation ( const fvMesh& mesh, const dictionary& dict, const word& phaseName ) : p_(lookupOrConstruct(mesh, "p")) {}