请教关于进口条件的设置的问题

CFD中文网

质量严重不守恒，为何不把网格做到10万以下，轻量级跑跑试试？

小葱小虫

@cfd-china 好的，我这么设置周期边界没有问题吧，两个面的距离是22mm

| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

8
(
    hub
    {
        type wall;
        nFaces          30748;
        startFace       5380555;
    }
    inlet
    {
        type patch;
        nFaces          5546;
        startFace       5411303;
    }
    outlet
    {
        type patch;
        nFaces          4543;
        startFace       5416849;
    }
    per1
    {
          type            cyclicAMI;
        nFaces          15340;
        startFace       5421392;
        neighbourPatch   per2;
        matchTolerance    0.01;
	transform         translational;
	separationVector (0 0.022 0);
    }
    per2
    {
         type            cyclicAMI;
       nFaces          15340;
       startFace       5436732;
       neighbourPatch   per1;
       matchTolerance    0.01;
       transform         translational;
       separationVector (0 -0.022 0);
    }
    ps
    {
        type            wall;
        nFaces          10266;
        startFace       5452072;
    }
    ss
    {
        type            wall;
        nFaces          11151;
        startFace       5462338;
    }
    sym
    {
        type symmetryPlane;
        nFaces          30748;
        startFace       5473489;
    }
)

// ************************************************************************* //

CFD中文网

为何不用cyclic？

per1
    {
          type            cyclic;
        nFaces          15340;
        startFace       5421392;
        neighbourPatch   per2;
    }
    per2
    {
         type            cyclic;
       nFaces          15340;
       startFace       5436732;
       neighbourPatch   per1;
    }

小葱小虫

@cfd-china 我感觉这个我理解的可能有问题，是应该给平移距离的大小（0.022m），还是（0,1,0）这么给。当时从这里看见的

separationVector (0 -0.022 0);

小葱小虫

@cfd-china 使用cyclic会有这个问题，因为两侧网格并不完全匹配

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-119cac7e8750
Exec   : simpleFoam
Date   : Mar 25 2017
Time   : 05:33:37
Host   : "ubuntu"
PID    : 48984
Case   : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/cynasa
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0



--> FOAM FATAL ERROR: 
face 47 area does not match neighbour by 10.4554% -- possible face ordering problem.
patch:per1 my area:1.63789e-10 neighbour area:1.47515e-10 matching tolerance:0.1
Mesh face:5421439 fc:(0.0411928 0.00446869 5e-07)
Neighbour fc:(0.040347 0.0264296 5e-07)
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with cyclic debug flag set for more information.

    From function cyclicPolyPatch::calcTransforms()
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 221.

FOAM exiting

小葱小虫

@cfd-china 我的算例可以算了，但是不知道为什么，壁面的分离算的特别小

小葱小虫

@小葱小虫我只是换了一套网格，其他没动

小葱小虫

@cfd-china 老师您好，由于我的算例进口马赫数为0.7，所以要考虑可压缩性，但是算着算着就发散了，用的rhosimplefoam，请问，能否帮我调试一下，由于我的上面是对称面，下面是壁面，所以网格10万很难满足壁面要求，能否帮忙调试，谢谢

Izumi

@小葱小虫你好！我在计算压缩机时，也出现了time step continuity errors过大导致发散的情况，参考苏军伟的博客，我调整了非正交修正次数，目前没发散，得到了结果。你可以参考一下。
主要就是修改fvsolution文件。

SIMPLE
{
    //nNonOrthogonalCorrectors 0;
    nNonOrthogonalCorrectors 1;
    residualControl
    {
        p               1e-2;
        U               1e-3;
        "(k|epsilon)"   1e-3;
    }
}

小葱小虫

@Izumi 是增加nNonOrthogonalCorrector这个数吗，谢谢我试试

Haining LUO

@小葱小虫
 @Izumi 应该说的就是nNonOrthogonalCorrectors

另外你两个网格之间差别是什么，我也是用pointwise画网格，出现aspect ratio那么大可能是手误也有可能是中间有breakpoint，你用examine aspect ratio看看那些畸形的网格，用connector->distribution->delete breakpoint去掉。因为如果一个方向的网格意外的小，虽然你的求解器是simpleFoam没有CFL数，但也对稳定性有影响吧？

宝丁

我也跟你遇到了相同的问题，不知道是不是计算设置的问题，总压入口是不是还有别的需要注意的地方

宝丁

你好，你现在解决这个问题了吗？我现在也遇到了这个问题，不过是因为pressureDirectedInletVelocity 这个进口条件的问题。
我也在做类似的研究，可以交流一下吗

CFD中文网

请教关于进口条件的设置的问题