CFD中文网

@李东岳 最近在将低版本OF的自定义边界条件移植到高版本OF的时候也发生了这个问题。串行求解没问题，并行求解报错。在不修改自定义边界条件源码的情况下，解决方案如下：1.在decomposePar前，将controlDict中引用的自定义库注释掉，例如

libs (
      "libOpenFOAM.so"
      //"libUserBC.so"
     );

2.执行decomposePar，3.取消注释，4.运行求解器，5.重新注释，6.执行reconstructPar。
以上，便可解决该报错。

如题，我在使用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是细导管的细节，可以看出网格表面呈现锯齿形状。

您好，我最近也遇到了相同的问题，请问您是怎么解决的呢？我已经在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;
}

@范准 已解决

感谢两位的讨论，请问平面的stl文件应该怎么生成呢？我只会生成3D的stl

楼主您好，请问您提到的“interFOAM一直有些问题，如交界面出现褶皱，交界面相对速度不正常的大”的信息是在哪里获取到的？我最近在自己的计算中也发现了这个问题，但是一直不知道如何解释以及如何避免。。。。谢谢~

我是这样做的：
在求解器中，添加这样的代码：

        IOdictionary meshProperties
        (
            IOobject
            (
                "meshProperties",
                runTime.constant(),
                mesh,
                IOobject::MUST_READ,
                IOobject::NO_WRITE
            )
        );
	scalar xnumber(int(readScalar(meshProperties.lookup("xnumber"))));
	scalar ynumber(int(readScalar(meshProperties.lookup("ynumber"))));

在算例文件夹的constant文件夹中，创建一个名为meshProperties的文件，内容为

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      meshProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

xnumber                       20;
ynumber                       10;

这样一来，参数xnumber和参数ynumber的值，就是通过算例文件夹里面的参数设定文件读取的了

@yfclark 刚刚开始苦逼的直博。。。。。。。

@东岳 举个例子，我希望把他从上面这样的网格，

变化成下面这样的网格

我试验过很多次，当实际底部边界的斜率比较大的时候，比如说0.2啊0.3啊这样的，画出来的网格边界就能够很好地保留实际边界的斜率，例如下图，

但是当我实际底部边界的斜率很小很小，例如此算例为0.002时，画出来就成了第一幅图那样了

@东岳 谢谢东岳前辈指导，我也发现随着网格的细化，这些台阶会越来越小，但是这也带来着计算量的迅速增加，怎样在保持原有网格数量的基础上，达到消除台阶的目的，这正是我目前死磕的地方。

是这样的，我画了一个长方体的背景网格，现在有一个斜率大概只有0.002的微微倾斜的边界面，是平面，我生成了STL格式的边界模型之后，用snappyHexMesh进行切割扣除，发现这个微微倾斜的斜面并不能够被很好地保留下这种几何特征，而是变成了下面这种，水平线-台阶-水平线-台阶-水平线，这样的形式，而不是一个连续的斜线。

不知道要怎么才能够完美保留微微倾斜变截面的几何特征，而避免这种台阶式的突变呢？谢谢各位前辈赐教~

@netbuger 哇！确实是这样的！我把这个选项取消掉就能够看得到了！谢谢前辈指点~~

@宝丁 我也是将边界设置成了不同的组！应该就是这个问题。多谢多谢~~~

不知道各位前辈有没有遇到过一种情况，就是在还没有进行计算的算例中（也就是只有0、constant、system文件夹。已经完成blockMesh和setFields等）打开ParaView查看初始条件，没办法加载出初始场的信息。
举个例子，在正常情况下，打开算例Apply的时候是这样的，我们可以在Cell Arrays中选择需要加载的场，

但是在某些算例中，Apply的时候是这样的

在Cell Arrays中并未显示任何可以加载的场。
不知道为什么会出现这种问题？

我目前找到的一个歪门邪道的解决办法是把0文件夹改名为1，然后再打开ParaVIew查看，查看完了把名字再改回0。

是的，现在我遇到的一个问题就是GCC的版本问题，还不知道怎么在单机状态下把GCC的版本更新了。。这个还需要我去探究

可以用一个名叫groovyBC的第三方插件，自由定义边界条件

@一二 哇，受教了受教了，谢谢前辈指点迷津~~~~~

在用fieldAverage工具的时候遇到了一个问题，比如说整个计算时间是0-10秒，假设这个流动5秒后才稳定，那么我只能对5-10秒内的流场进行时均，请问这应该怎么做呢？

我目前想到的非正式解决办法是，把整个计算分成两个算例，第一个算例正常计算0-5秒，然后第二个算例从5秒开始计算5-10秒并在计算过程中时均。不知道还有没有更巧妙些的办法呢？

希望各位前辈不吝赐教~

@东岳
哇，前辈好~~~