CFD中文网

icem网格划分导入fluent时检查失败，显示Error: interface zone 3 has two adjacent cell zones.
4f68eb65-132e-4314-b0f6-6964dc466673-image.png

@李东岳 明白了，谢谢李老师。排版排了好久，清楚一点方便讨论:papa:

https://www.cfd-online.com/Forums/openfoam-meshing/133850-snappyhexmesh-problems.html

teoria.png
如图所示，这是非常经典的Le Méhauté (1976)波浪理论适用图，但我有点迷糊，主要有两方面疑问：

图中红色部分看起来是Stokes三阶波/四阶波与五阶波的交接处，那么我应该选择用三阶波/四阶波还是五阶波呢？同样地，图中浅蓝色和五阶波也有交集，那么在既是浅蓝色又是五阶波的区域应该用cnoidal还是五阶波呢？ 我还注意到图中五阶波下面还有一个流函数理论，但我看olaFlow的用户手册，用流函数的话需要输入一大堆参数（如下图所示），这些参数应该怎么选取呢？本菜鸡在此请教一下各位:135:
2024-07-08 21-55-00 的屏幕截图.png
2024-07-08 22-01-59 的屏幕截图.png

李老师，我用的是可变时间步长，保证CFL数小于1，所以应该不是时间步长的问题，我注意到计算是计算到某一步之后rhoEqn的最小值突然变成0了，这可能是什么原因呢？

rhoEqn max/min : 1.62584 0

thermophysicalProperties中energy可以选择sensibleEnthalpy, sensibleInternalEnergy 和 absoluteEnthalpy。前两者与后者的区别，后者使用绝对能量表示包含热源的情况。而sensibleEnthalpy与sensibleInternalEnergy的区别，前者是利用焓值求解温度，后者是利用内能求解温度。热⼒学中分析闭式系统⽤内能，分析开式系统⽤焓值。那么在实际使用中，该如何选择这两种模式？例如高压离心风机。tutorials中的案例似乎找不出规律。
eb324645-1dd5-499a-bfc0-6d77667f7fa5-image.png

wall film 涉及到rebound、splash、adhesion等物理过程，必然是只对液体起作用。

1.我猜是symmetry这个边界条件，试一下改成wall，然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查，外部软件的网格节点编号和blockMesh的逻辑不一样。
2.可能是前后的empty边界问题，你可以在starccm里生成3D网格后，重新在openfoam里extrude来生成前后empty边界
类似这样
constructFrom patch;
sourceCase "../w3-d_hc1-3";
sourcePatches (symFront);

// If construct from patch: patch to use for back (can be same as sourcePatch)
exposedPatchName symBack;

// Flip surface normals before usage. Valid only for extrude from surface or
// patch.
flipNormals false;

//- Linear extrusion in point-normal direction
extrudeModel linearNormal;

nLayers 1;

expansionRatio 1.0;

linearNormalCoeffs
{
thickness 1;
}

// Do front and back need to be merged? Usually only makes sense for 360
// degree wedges.
mergeFaces false; //true;

// Merge small edges. Fraction of bounding box.
mergeTol 0;

@李东岳 感谢回复。我是想先写个简单的特定算例学习一下算法，格式边界条件啥的都直接写死，然后再慢慢做的通用一点，能搞出来肯定开源。

@李东岳 感谢李老师，我试试看能不能画出来

在构建类层次结构是，我们通常都是想要获得多态行为。我们需要对象一旦创建就能保持其类型和标识。但是在想复制一个对象时会发现C++中的对象复制语义与这一目标是冲突的，因为构造函数都是与类名相同的，也就是构造函数都是每个类特有的，复制构造函数不能设为虚拟的。
为了用基类指针（或引用）正确的将对象进行复制，更为通用的惯用法解决方案是，将基类的复制构造函数设为protected甚至是public，并改而依赖于一个虚拟的clone函数。例如particle类中

//- Construct as a copy particle(const particle& p); //- Construct a clone virtual autoPtr<particle> clone() const { return autoPtr<particle>::New(*this); } //Constructor Foam::particle::particle(const particle& p) : mesh_(p.mesh_), coordinates_(p.coordinates_), celli_(p.celli_), tetFacei_(p.tetFacei_), tetPti_(p.tetPti_), facei_(p.facei_), stepFraction_(p.stepFraction_), origProc_(p.origProc_), origId_(p.origId_) {} 在派生类KinematicParcel中 //- Construct as a copy KinematicParcel(const KinematicParcel& p); //- Construct and return a (basic particle) clone virtual autoPtr<particle> clone() const { return autoPtr<particle>(new KinematicParcel(*this)); }

通过虚拟函数clone就能正确的实现复制功能啦！

@King_RVM 你好，使用starCCM绘制网格，然后导入到openFoam里面跑算例，计算结果会出现质量不好的情况么？请教一下~

@史浩 他那个应该是个笔误。写反了。

我这周刚把原始的FCT算法植入进去了，OpenFOAM的算法跟原始的不太一样，主要在于

最大值和最小值的确定，这对应不同的方法，参考Zalesak1979 342页末

反扩散系数的确定（你贴的按个公式）

这也是他之所以叫做MULES的原因吧，不过思路都是一样的，Weller为了加快计算速度也优化的非常好。原始FCT已经完全可以对变量进行有界了，我对比原始的FCT和MULES差异非常小，原始FCT也不需要对$\lambda$进行迭代。我要用原始FCT验证几个算例分析一下。

在测试可压缩传输的时候，我发现可压缩FCT和不可压缩FCT的密度处理需要特殊处理。也就是说 http://www.cfd-china.com/topic/1029/mules算法中的const-rhotype-rho/3 帖子中的可压缩不能直接这么处理，要把rho和alpha分开。不能看做对rhoAlpha的传输。具体的我还在研究

下图FCT1 和FCT2对应原始FCT和MULES，差别微小。抽取的damBreak的alpha的最小值

捕获.JPG

居然还发公众号了

ffmpeg -i input.MP4 -c:v libx264 -preset slow -crf 40 -c:a copy output.mp4

是的，一般情况下你1）可以忽略correctPhiH，2）pcorr也可以不求解

需要写一个file文件，可以自定义颗粒物的状态。关于file文件怎么写，可以参照fluent帮助文件

大家好，请问paraview如何获得三维流线图？如下图这样的，用户指南上关于流线图的说明仅仅是寥寥几笔，请各位不吝赐教，谢谢！
feaadcd8-3566-4e3f-831c-5e9a11ca59e0-图片.png

今天收到一封邮件，邀请我审一本书。

Inked1_LI.jpg
这个是fluent自带的udf帮助文件里面的一个例子，也用到了NV_DS，可以对于比一下