@coolhhh 厉害，测试代码copy了。不过，这里直接操作就又可以了:wocao:

OpenFOAM-v2206/tutorials/incompressible/pimpleFoam/LES/planeChannel/setups.orig/common/system/controlDict，有一个计算Cf场的函数，可以参考下：

Cf { type coded; libs (utilityFunctionObjects); name Cf; timeStart $/timeStart; writeControl writeTime; codeExecute #{ auto* CfPtr = mesh().getObjectPtr<volScalarField>("Cf"); if (!CfPtr) { Info<< "Create skin-friction coefficient field" << nl; CfPtr = new volScalarField ( IOobject ( "Cf", mesh().time().timeName(), mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh(), dimless ); regIOobject::store(CfPtr); } auto& Cf = *CfPtr; Info<< "Computing skin-friction coefficient field\n" << endl; const auto& tau = mesh().lookupObject<volVectorField>("wallShearStress"); const dimensionedScalar Ubulk(dimVelocity, 17.55); Cf = mag(tau.component(0))/(0.5*sqr(Ubulk)); #}; }

compressibleInterFoam其实自带了一个，就在tutorials里面，好像交bubbleCollapse

@Jasper-0 您好，请问这个问题您现在解决了吗？

@chenboyao 可以试试把fvsolution中pimple的动量预测关了

边界条件应该没啥问题，就是接触角边界条件，我用的时候可能会搞反0°和180°，你自己最好确定一下。
这种网格，大网格和小网格之间的变化是2倍的变化，有些地方生成的可能也不是很好。反正你现在是二维，要不就全加密到最小的尺寸？

我在CFD-Online中找到了解决办法，Particle tracking and AMI patches error for parallel processing，供大家参考

感谢，这个回答救命了

使用的是1912版本

各位老师、专家，请问我在计算物体移动过程中想要了解他局部位置的受力情况该怎么办？能从parafoam中提取吗

如题，我在使用snappyHexMesh根据stl文件对连接两个区域的细导管进行网格划分时，在非加密区贴合良好，而在加密区出现了网格与表面不贴合，即锯齿的情况，请各位不吝赐教，谢谢大家。
我查看了论坛以往的网格贴合问题的帖子，基本都是全局不贴合，而非局部不贴合，与我的情况不太一样。我的面贴合参数snapControls如下，

snapControls { nSmoothPatch 5; nSmoothInternal 5; tolerance 3; nSolveIter 500; nRelaxIter 5; nFeatureSnapIter 10; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap false; }

已经尝试过增加面贴合迭代次数和面平滑迭代次数等，几乎没有效果，是按照https://zhuanlan.zhihu.com/p/295681514的指导进行的调整。整体算例文件在附件中给出，./Allrun一键运行即可。
snapTest.zip
图1是光滑的原始stl文件，图2是使用snappyHexMesh生成的贴体网格，图3是细导管的细节，可以看出网格表面呈现锯齿形状。
stl.png
whole.png
local.png

我是在多孔介质里面模拟气液两相流，为什么我的速度流线是这样的呢，蓝色气泡位置的流速很大
image.png

最近看《Holzmann T. Mathematics, numerics, derivations and OpenFOAM》，里边也提到大部分情况忽略k，对结果影响不大：
05a3730d-81ea-47f7-93d1-d3e7f8814605-image.png

@疏影横斜水清浅 你好，可以咨询一下这个的相关问题吗

另外，CFDEM如何实现壁面变形呢，像下面这样逐渐撑开流动通道
image.png

@Tens 请问下 如何在算例中调用 src\thermophysicalModels\thermophysicalProperties\liquidProperties的物性

http://dyfluid.com/datadrivenles.html 李老师方便提供完整编译的文件么？最近在编译数据驱动RANS模型，但是编译的时候，总是报错，libtorch与openfoam的typename冲突，不知道咋解决。

就是速度入口，柱子面积和入口面积的比×表观速度

您好，我最近也遇到了相同的问题，请问您是怎么解决的呢？我已经在snappyHexMeshDict中指定了面贴合参数

snapControls { nSmoothPatch 5; nSmoothInternal 5; tolerance 2; nSolveIter 300; nRelaxIter 5; nFeatureSnapIter 10; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap false; }

@DY大世界 您好，我最近也遇到了相同的问题，请问您是怎么解决的呢？我已经在snappyHexMeshDict中指定了面贴合参数

snapControls { nSmoothPatch 5; nSmoothInternal 5; tolerance 2; nSolveIter 300; nRelaxIter 5; nFeatureSnapIter 10; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap false; }