blockMesh楔形体网格划分

Joker

QQ截图20240505123429.jpg
用blockMesh划分这种楔形体网格的时候，最左侧那条轴线上的边界该如何定义呢？像这种三棱柱类型的体又该如何定义呢？

李东岳

正儿八经的wedge左边是一条线，定义成axis。但interFoam那面有个算例叫climbing rob，左边不是一条线，你看一下他定义成了什么

Joker

@李东岳他定义成了一个面，然后用了一个extrudeMeshDict文件对网格进行转化，我正在尝试，谢谢老师

李东岳

我记得好像也是axis。你看看。手头没开虚拟机。

Joker

@李东岳 ```
/--------------------------------- C++ -----------------------------------\

\ / F ield	OpenFOAM: The Open Source CFD Toolbox
\ / O peration	Website: https://openfoam.org
\ / A nd	Version: 10
\/ M anipulation

*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class dictionary;
object extrudeProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

constructFrom patch;
sourceCase "$FOAM_CASE";

sourcePatches (front);
exposedPatchName back;

extrudeModel wedge;

sectorCoeffs
{
axisPt (0 0 0);
axis (0 -1 0);
angle 5;
}

flipNormals false;
mergeFaces false;

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

老师，这是我的extrudeMeshDict文件，我用这个文件进行拉伸，它提示我的wedge不是一个平面，为什么会这样呢？

/--------------------------------- C++ -----------------------------------\

\ / F ield	OpenFOAM: The Open Source CFD Toolbox
\ / O peration	Website: https://openfoam.org
\ / A nd	Version: 10
\/ M anipulation

*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.001;

vertices
(
(0 0 0)//0
(7 0 0)//1
(10 0 0)//2
(1010 0 0)//3
(2010 0 0)//4

(0 5 0)//5
(5 5 0)//6
(10 5 0)//7
(1010 5 0)//8
(2010 5 0)//9

(0 40 0)//10
(7 40 0)//11
(10 40 0)//12
(1010 40 0)//13
(2010 40 0)//14

(0 1640 0)//15
(7 1640 0)//16
(10 1640 0)//17


(0 0 10)//0
(7 0 10)//1
(10 0 10)//2
(1010 0 10)//3
(2010 0 10)//4

(0 5 10)//5
(7 5 10)//6
(10 5 10)//7
(1010 5 10)//8
(2010 5 10)//9

(0 40 10)//10
(7 40 10)//11
(10 40 10)//12
(1010 40 10)//13
(2010 40 10)//14

(0 1640 10)//15
(7 1640 10)//16
(10 1640 10)//17

);

blocks
(
hex (0 1 6 5 18 19 24 23) (20 50 1) simpleGrading (1 ((1 1 50)(0 0 0)) 1)//1

hex (1 2 7 6 19 20 25 24) (30 50 1) simpleGrading (((1 1 0.05)(0 0 0)) ((1 1 50)(0 0 0)) 1)//2

hex (2 3 8 7 20 21 26 25) (300 50 1) simpleGrading (((1 1 1500)(0 0 0)) ((1 1 50)(0 0 0)) 1)//3

hex (3 4 9 8 21 22 27 26) (40 50 1) simpleGrading (1 ((1 1 50)(0 0 0)) 1)//4



hex (5 6 11 10 23 24 29 28) (20 100 1) simpleGrading (1 1 1)//5

hex (6 7 12 11 24 25 30 29) (30 100 1) simpleGrading (((1 1 0.05)(0 0 0)) 1 1)//6

hex (7 8 13 12 25 26 31 30) (300 100 1) simpleGrading (((1 1 1500)(0 0 0)) 1 1)//7

hex (8 9 14 13 26 27 32 31) (40 100 1) simpleGrading (1 1 1)//8



hex (10 11 16 15 28 29 34 33) (20 400 1) simpleGrading (1 ((1 1 50)(0 0 0)) 1)//9

hex (11 12 17 16 29 30 35 34) (30 400 1) simpleGrading (((1 1 0.05)(0 0 0)) ((1 1 50)(0 0 0)) 1)//10

);

boundary
(
inlet
{
type patch;
faces
(
(15 33 34 16)
(16 34 35 17)
);
}

outlet
{
    type patch;
    faces
    (
        (4 9 27 22)
        (9 14 32 27)
    );
}

upper
{
    type symmetry;
    faces
    (
        (12 30 31 13)
        (13 31 32 14)
    );
}

axis
{
    type symmetry;
    faces
    (
        (0 18 23 5)
        (5 23 28 10)
        (10 28 33 15)
    
    );
}

bottom
{
    type wall;
    faces
    (
        (0 1 19 18)
        (1 2 20 19)
        (2 3 21 20)
        (3 4 22 21)

    );
}

wall
{
    type wall;
    faces
    (
        (12 17 35 30)

    );
}

back
{
    type symmetry;
    faces
    (
        (0 5 6 1)
        (1 6 7 2)
        (2 7 8 3)
        (3 8 9 4)
        
        (5 10 11 6)
        (6 11 12 7)
        (7 12 13 8)
        (8 13 14 9)
        
        (10 15 16 11)
        (11 16 17 12)
        
        
        
    );
}

front
{
    type symmetry;
    faces
    (
        
        
        
        (18 19 24 23)
        (19 20 25 24)
        (20 21 26 25)
        (21 22 27 26)
        
        (23 24 29 28)
        (24 25 30 29)
        (25 26 31 30)
        (26 27 32 31)
        
        (28 29 34 33)
        (29 30 35 34)
    );
}

);

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

这是我的blockMesh文件，根据你之前的提示，我看那个算例里面是先定义了一个矩形出来，然后进行拉伸的，想仿照人家的来做

CFD中文网

blockMesh楔形体网格划分