@yanoo 考虑一下看是不是编译器和依赖项的问题? OpenFOAM-8竟然也会出问题非常奇怪诶...

@chon01 哈哈，新年快乐
我不是前辈，今年刚博士一年级，大家互相交流学习！

求大佬修改一下，指出问题，假设每秒注射1ug的颗粒，粒径分布为RosinRammler，期望值为70e-6

injectionModels { model1 { type patchInjection; massTotal 0;//应该怎么设置呢，这里设置为总时间的注射量吗？ nParticle 10000; SOI 0; parcelBasisType fixed;//改成mass？ patch inlet; duration 100000; parcelsPerSecond 0; U0 (5 0 0); flowRateProfile constant 1; sizeDistribution { type RosinRammler;//假设 RosinRammlerDistribution { minValue 0.1e-6; maxValue 200e-6; d 70e-6; n 3;//这个含义是？ } } }

}

@Zephyr 我直接好家伙，楼上4年后回复你了

论文里面提到3D 玫瑰图是用来说明各向异性分布以及各向异性大小和主方向的。这些柱状图好像都是对称的，每个柱子的高度应该是这个方向的平均力的大小，颜色是根据数值在柱子不同高度处分别着色。（https://doi.org/10.1016/j.compgeo.2023.105692）

1-s2.0-S0266352X23004494-gr20_lrg (1).jpg

不仅如此，我还发现，tecplot有时候会乱码
左图tecplot右图paraview
ee7532199637e2fc3c875e7ed7bd5e8.png

你的网格是不是太大了导致内存爆了啊? 新开一个终端窗口运行

top

来监视你的硬件资源使用情况吧

@李东岳 谢谢李老师，联系卖的人了
@bestucan 老师，我写错了，是自动重启了，不是自动关机了。内存是256G，测试的网格是八十多万感觉不应该把内存爆了，因为测试96核和72核自动重启，但是对56核就没问题，单个CPU是64核的

@李东岳 在 关于openfoam与结构动力学的问题 中说：

楼上说的对。小变形比较简单，也可以自己写一个

好嘞，谢谢李老师！

建议运行

wmake libso $FOAM安装路径/src/finiteVolume

你这个文件夹里面就没有Make文件夹，正常来说就编译不了

@LeeSimon 在 请教一下，用interFoam计算一个有盖的渡槽（进出口均与明渠段相接），保证进口流量和水位一定的情况下，如何设置边界条件 中说：

请教一下解决了吗？

出口给个流速或者根据静压假定给压强

挺好，老铁

@dzw05 请问老师，您这边清楚interface并行边界如何处理吗？我想从openfoam中拿到完整矩阵数据，但是并行的话只能拿到各个线程对应的子矩阵，如何知道并行边界的关系将其组装成一个完整的矩阵

@李东岳 fluidfoam用于抽取数据进行后处理的，谢谢回复，刚刚问题解决了，将python文件放到上一级文件夹中运行即可

@benqing 你把红框那一串，通过下面类似的代码，起个名字

volScalarField("mingzi", ....);

然后在算例层面，通过输出cache的方法进行输出，参考 http://dyfluid.com/code.html

@ZY-GONG 好的，多谢老师，辛苦了

好的，感谢

刚才更新了一些在笔记上

代码里面为了减少代码复用会在基类中定义大量的变量。在衍生类中有一些用不到。每个算法都是一个衍生类。因此有一些算法会用不到一些变量。你自己编译状态方程，你就定义你自己需要的。

return neg(Re - 1000)*24*(1.0 + 0.15*pow(Re, 0.687)) + pos0(Re - 1000)*0.44*max(Re, residualRe_);

https://www.openfoam.com/documentation/guides/v2206/api/reactingEuler_2multiphaseSystem_2interfacialModels_2dragModels_2SchillerNaumann_2SchillerNaumann_8C_source.html

结合SchillerNaumann阻力模型的公式，你看一下，公式我笔记里面有