@浪迹天大 老师您好，按照这个方法，我需要更改controlDict中的maxDeltaT。但是会发现文件更新不起作用。

/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application interFlowVb1; startFrom latestTime; startTime 0; stopAt endTime; endTime 0.3; deltaT 0.00005; writeControl adjustableRunTime; writeInterval 0.05; purgeWrite 0; writeFormat binary; writePrecision 8; writeCompression off; timeFormat general; timePrecision 8; runTimeModifiable yes; adjustTimeStep yes; maxCo 0.5; maxAlphaCo 0.5; maxCapillaryNum 10; // stability criteria is 1 maxDeltaT 2e-5;//controlDict.a; // maxDeltaT 1e-5;//controlDict.b; // maxDeltaT 5e-6;//controlDict.c; functions { fileUpdate1 { type timeActivatedFileUpdate; libs ("libutilityFunctionObjects.so"); writeControl timeStep; writeInterval 1; fileToUpdate "$FOAM_CASE/system/controlDict"; timeVsFile ( (-1 "$FOAM_CASE/system/controlDict.a") (1e-4 "$FOAM_CASE/system/controlDict.b") (2e-4 "$FOAM_CASE/system/controlDict.c") ); } }

三个文件只有maxDeltaT有差别，请老师捞捞@李东岳 @浪迹天大

运行后结果提示：

timeActivatedFileUpdate: copying file "~/orifice3Da/system/controlDict.c" to: "~/orifice3Da/system/controlDict" Courant Number mean: 2.0171174e-08 max: 0.3369988 Interface Courant Number mean: 2.978204e-10 max: 0.22052861 Capillary Number: 2.4526254 deltaT = 2e-05 Time = 0.00022

文件有更新，并且maxDeltaT也更新为了1e-5; 但是它没有发挥作用，可以看到文件更新后的deltaT 依然是 2e-05，如果是我自己运行过程中手动改小maxDeltaT，它会马上从deltaT = 2e-05，便为deltaT = 1e-05。

@JC 在 在interfoam两相流中，如何计算相对渗透率 中说：

@JC 在 在interfoam两相流中，如何计算相对渗透率 中说：

有一个标准里面写有计算方法，记录每个时刻流出的两相体积，以及压力梯度，应该就可以算出来了。对了还需要这个多孔介质的初始渗透率

《SY/T 5345-1999》一个很老的标准，可以参考下

谢谢，后面我也是直接用出口流量来算每一相渗透率再计算相对渗透率，部分文献里也存在根据饱和度计算相对渗透率的模型，不过不好确定误差。

@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 好的，多谢老师，辛苦了

好的，感谢

刚才更新了一些在笔记上