有没有大佬知道怎么求颗粒的剪切应力？类似流体相的UPrime2Mean？

cpu型号：单路9654，主板超微
内存：16g*12 ddr5
系统：Ubuntu20.04
OpenFOAM版本：wsl-openfoam-9

#cores Wall time (s):
32 165.68
16 194.41
8 191.77
4 289.18

有一个比较基础的问题想要请教：仿照其他函数格式构造下列函数，目的是调用Tv温度场，必须加入const Thermo& thermo才能实现thermo.Tv() 调用，但不了解const Thermo& thermo是什么变量？Tvenewdefinition该如何调用？
scalar Tvenewdefinition
(
const Thermo& thermo
) const
{
return thermo.Tv();
}

@李东岳 李老师感谢您的及时回复
在OpenFOAM中液体属性的adiabaticPerfectFluid模型，使用的是Tait状态方程这个模型
bdfe6710-0a9a-4fa8-a960-c1d2b345e7e1-image.png
其中B的值通过实验测得为3e8，在这种情况下如果绝对负压不超过B的值，就不会出现负密度，这样应该也不会影响到计算吧？我不知道这样理解对不对

也有相关领域的液相绝对负压的文献：CAUPIN, F. & HERBERT, E. 2006 Cavitation in water: a review. C. R. Phys. 7 (9–10), 1000–1017.

好的老师，我再调试调试:144:

压力梯度项前面是否包含体积分数，应该跟采用的模型是A模型还是B模型有关

multiphaserulerfoam是存在界面压缩项。这个东西同时存在drag，同时存在界面压缩。

Wardle那个最开始就是跟Weller开发的multiphaseeulerfoam

@capillaryFix 下午好，铁子。您那边有新的解决方案和想法吗？

当前化学反应速率均只包含一个控制温度，想更改成与两个温度相关，通过插值读表获取反应速率。编写新函数inline scalar HoffertLienReactionRate::interpolate(double T) 实现上述功能，代码可以成功编译，但输出Tv检查时发现其值一直为0，这是不是表示未能成功调用？具体代码如下。向各位大佬求助！！！HoffertLienReactionRate.H HoffertLienReactionRateI.H

@dzw05 因为limitSum仅仅对anti diffusion flux进行限制，在限制完后再和低阶格式组装成总通量:134:

@Imech 在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：

@学流体的小明 老哥，这个问题你解决了吗，我最近算完美滑移边界发现算的结果与参考对不上

最近还在搞，你是什么问题？算例是什么？

@Zhy2022 您好，frictionalStressModel和granularPressureModel里面的模型我都试了一下,并且有的跟fluent选的一样，但还是这个问题。我能否把算例分享给您，你帮我看一下，我的模型网格是20倍粒径，单核跑起来很快，不知道怎么分享给您合适？还有一个问题就是，我如果不把物性参数改成water，用官方算例air，结果就是正常的，能达到紧密堆积0.6.不确定是不是我物性参数改的有问题？faac38a7-acef-406c-8ae4-47b7d5bbfd64-image.png

@队长别开枪 感谢大佬回复，这就去学习一下！

compressibleInterFoam自带的tutorial里面就有

为什么在controlDict中把writeInterval调大，计算就发散？？？为什么会有这么诡异的bug
writeInterval不就是控制写入频率的么，跟计算有关系吗？
有没有大佬遇见过
7f83b6e1550cdbc87dc895138ddbea5.png 1748766340208.png

@队长别开枪 好的，感谢！

解决问题了，主要是openFoam的版本问题，切换到低版本（2206以下），就可以实现边界位移了

image.png

@李东岳 好的李老师，搞定了