@sxz0823 楼上大佬做湍流干了1篇JCP 1篇POF 相当流弊

你加的那些私有成员，声明了么

@李东岳 我在KinematicParcel.H文件，以及这个文件的前面，我都对我新加变量声明了，不然不声明这些东西源代码都无法进行编译，会有很多error。（目前只不过编译时候有一些warning，连着源代码以前自带的一些东西都在warning）。
目前这个问题是我开始编译求解器时候，求解器识别源代码时候出现的情况。

我再看看我代码声明问题吧，我也不知道和这个后没有关系

另：之前我代码是这样

template<class ParcelType> inline Foam::KinematicParcel<ParcelType>::KinematicParcel ( const polyMesh& owner, const vector& position, const label cellI, const label tetFaceI, const label tetPtI ) : ParcelType(owner, position, cellI, tetFaceI, tetPtI), active_(true), typeId_(-1), nParticle_(0), d_(0.0), dTarget_(0.0), U_(vector::zero), rho_(0.0), age_(0.0), tTurb_(0.0), dRt_(0.0), //- myadd Ro_(0,0), //- myadd pgo_(0.0), //- myadd Pc_(0.0)， //- myadd UTurb_(vector::zero), rhoc_(0.0), Uc_(vector::zero), muc_(0.0) {}

他求解器会报错

/home/zly/OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/intermediate/lnInclude/KinematicParcelI.H:105:12: error: expression list treated as compound expression in mem-initializer [-fpermissive] muc_(0.0)

看起来只报错最后一个地方，我也不知道为什么。这些代码的声明和修改，一部分是我看网上的一篇资料修改，另一部分是我参考着原始代码格式，进行改写。

@李东岳 李老师您好，之前确实没考虑到压力波。不过我这里的泄压过程是缓慢过程，是通过密闭空间周围的多孔介质渗流漏出去的，给的流速很低，大概在10e-4m/s量级，也会有这种问题么？然后考虑到空间内的热对流，我可能更倾向于使用buoyantPimpleFoam求解器？

这个数值单位是m（米）吗，如果是，这个距离远大于建立的模型的尺寸，那么这两个关键词的意义是什么。

是米。如果过大，就没有意义了。

@tidedrinker 在 服务器计算速率问题咨询。 中说：

其他没看出来，内存是不是插满比较好？
另外，有些时候，同时算两个算例平均时间会比仅算一个快

也在琢磨扩充内存的问题，感谢。

老师您好，我想请问一下关于字典中关键词的含义
innerDistance 100;
outerDistance 101;

再找到的一个介绍中这样解释：innerDistance 内的任何内容直接将网格节点作为刚体移动。在innerDistance 和outerDistance 之间，网格节点变形。在outerDistance 之外，不会发生变形。outerDisatance永远大于innerDistance。

想请问一下老师，这个数值单位是m（米）吗，如果是，这个距离远大于建立的模型的尺寸，我在很多算例都看见是这么设置的，那么这两个关键词的意义是什么。

之前snappyHexMesh画的3D网格进行重叠网格计算没有问题。后来为了节省计算成本用的2D网格，用的snappyHexMesh画的网格，然后使用extrudeMesh将3D网格变成2D网格后就出现问题了

李老师您好，这是两个文件里的，第一个是定义C4H6质量分数，第二个是定义边界速度

@江山不如凉竹 最近也看到一篇类似的文献，感觉思路挺好的，请问这个问题解决了吗？很好奇如何实现。

截图_20240705140532.png

@bestucan 在 求助：使用interDyMFoam计算船舶浮态 中说：

是不是水面起浪了？看上飘和下沉会不会持续，还是回归

太感谢了，确实是流速导致的起浪了，我把流速调为0以后就不会出现这种问题了。

@coolhhh 对的，如果湍流入口可能旋转比较方便，对于目前这个问题，乘以角度就可以了。谢谢大佬:146:

想请教一下有没有人用过openfoam中旋流入口边界条件swirlInletVelocity，我用该边界条件计算了一个冷态流场如图所示为轴向速度的云图。
1.png
我查看了源代码中关于该边界条件的使用方法如下
2.png
只需要指定轴向速度 径向速度以及切向速度，我想知道这个边界条件是如何计算旋流数的，因为我貌似在源代码中没有查看到如何通过旋流数去确定这三个方向的速度值。不知道有没有人知道openfoam中是怎么可以通过旋流数去计算得到这个三个速度的量？

@huangyuhui723 我尝试了可以通过U.boundaryField().oldTime()进行调用，可供参考。

我建议你不要上来先植入这个方程。而是做一个直径与时间的分段函数，现把这个路走通。再植入你的方程。

那可能就是影响太小了，主动把y方向的速度增大10倍看一下

@李东岳 有的老师 dynamicmeshdict里面我漏了一行设置 现在解决了

在使用dsmcFoam+求解器时，定义了一个旋转边界条件，来流气体分子均加上了一个绕z轴旋转的角速度 omega,当计算稳定之后，使用tecplot在入口处画Streamtraces，结果如图所示，貌似加上的角速度与预期不一致，有同行遇见过类似的问题吗？

725a22b4-d894-400e-9a70-cf5d266706aa-image.png

请问您的这个问题解决了，我在改进SST的时候一直提示我重复定义函数

最新进展：我将0文件中的pd边界条件进行了修改，可以实现中断后继续运算，但原问题依旧不知该如何解决。

代码如下：

/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | foam-extend: Open Source CFD | | \\ / O peration | Version: 4.0 | | \\ / A nd | Web: http://www.foam-extend.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object pd; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { out { //type zeroGradient; type fixedValue; value $internalField; } "(sym1|sym2)" { type symmetryPlane; } in { //type fixedFluxPressure; //value $interField; type zeroGradient; } cylinder { //type fixedFluxPressure; //value $interField; type zeroGradient; } "(front|back)" { type empty; } } // ************************************************************************* //

@chapaofan 好的，谢谢