OpenFOAM手动网格分块

OpenFOAM手动网格分块

1. 通过OpenFOAM的格式自编程来对各个cell命名，确定他们所属的processor

/*---------------------------------------------------------------------------*\

  =========                 |

  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox

   \\    /   O peration     |

    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation

     \\/     M anipulation  |

-------------------------------------------------------------------------------

License

    This file is part of OpenFOAM.

 

    OpenFOAM is free software: you can redistribute it and/or modify it

    under the terms of the GNU General Public License as published by

    the Free Software Foundation, either version 3 of the License, or

    (at your option) any later version.

 

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT

    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License

    for more details.

 

    You should have received a copy of the GNU General Public License

    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.

 

Application

    scalar2Label

    codeBy CloudBird, HIT,cloudbird7@foxmail.com

 

\*---------------------------------------------------------------------------*/

 

#include "fvCFD.H"

 

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

 

 

int main(int argc, char *argv[])

{

    #include "setRootCase.H"

    #include "createTime.H"

    #include "createMesh.H"

 

    const volVectorField& center = mesh.C();

 

    labelIOList finalDecomp

    (

        IOobject

        (

            "manualDecomposition",

            mesh.facesInstance(),

            mesh,

            IOobject::NO_READ,

            IOobject::AUTO_WRITE,

            false

        ),

        labelList (center.size())

    );

 

    Info<< "Making manualDecomposition file\n" << endl;

    scalar x_=0,y_=0;

label num0 = 0, num1 = 0, num2 = 0, num3 = 0;

//这里是将网格分为四块，分块的依据是下边if的判断语句，也就是先挖一个小方框出来，然后对剩下的部分分成三块，并输出分块之后每一块的cell的数量。

forAll(center, cellID)

    {

x_ = center[cellID].x();

y_ = center[cellID].y();

if (x_ > -0.75 && x_ < 0.75 &&y_>-0.75&&y_ < 0.75)

{

finalDecomp[cellID] = 0;

num0++;

}

else

{

if (x_ < 0)

{

finalDecomp[cellID] = 1;

num1++;

}

else if (x_ < 3)

{

finalDecomp[cellID] = 2;

num2++;

}

 

else

{

finalDecomp[cellID] = 3;

num3++;

}

}

 

    }

Info << "num0 = " << num0 << nl << "num1 = " << num1 << nl << "num2 = " << num2 << nl << "num3 = " << num3 << nl;

    finalDecomp.write();

 

    Info<< "Done\n" << endl;

 

    return 0;

}

 

 

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

这样会在constant文件夹下边生成一个名为manualDecomposition的文件，这个文件本质上是一个标量文件，用于标记各个cell所属的processor。

2. 在system文件夹下的decomposeParDict中修改分割方式

numberOfSubdomains 4;

method           manual;

manualCoeffs

{
    dataFile        "manualDecomposition";
}

其中dataFile的名字要改成我们生成的文件的名字，也就是maualDecomposition，还需要注意的一点是numberOfSubdomains的数量要和你在maulDecomposition中定义的数量一致。

3.之后直接decomposePar就是按照手动分割网格的方式对你定义的网格进行分割。
processor0

processor1

processor2

processor3

文件的编译按照一般OpenFOAM文件编译的方式进行，参考http://dyfluid.com/docs/book/_book/

参考文章：

https://blog.csdn.net/CloudBird07/article/details/105349031

https://www.dazhuanlan.com/2019/08/16/5d55fc7975c16/