做wedge网格做了一个下午

[李东岳]

openfoam wedge网格真是弱爆了

convertToMeters 1;

r  0.025;

Xpoint  #calc "$r*cos(1.5*3.1415/180)";
yPos  #calc "$r*sin(1.5*3.1415/180)";
yNeg  #calc "-$r*sin(1.5*3.1415/180)";

vertices
(
    (0   0    -0.02)//0
    (0   0    0)//1
    (0   0    0.7)//2
    ($Xpoint $yNeg  -0.02)//3
    ($Xpoint $yNeg  0)//4
    ($Xpoint $yNeg  0.7)//5
	(#calc "0.05*cos(1.5*3.1415/180)" #calc "-0.05*sin(1.5*3.1415/180)" 0)//6
	(#calc "0.05*cos(1.5*3.1415/180)" #calc "-0.05*sin(1.5*3.1415/180)" 0.7)//7

    ($Xpoint $yPos  -0.02)//8
    ($Xpoint $yPos  0)//9
    ($Xpoint $yPos  0.7)//10
	(#calc "0.05*cos(1.5*3.1415/180)" #calc "0.05*sin(1.5*3.1415/180)" 0)//11
	(#calc "0.05*cos(1.5*3.1415/180)" #calc "0.05*sin(1.5*3.1415/180)" 0.7)//12
);

blocks
(
    hex (0 1 4 3 0 1 9 8) (10 10 1) simpleGrading (1 1 1)
    hex (1 2 5 4 1 2 10 9) (10 10 1) simpleGrading (1 1 1)
    hex (4 5 7 6 9 10 12 11) (10 10 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    inlet
    {
        type patch;
        faces
        (
            (0 8 3 0)
        );
    }
	inletF
    {
        type patch;
        faces
        (
            (4 9 11 6)
        );
    }
    outlet
    {
        type patch;
        faces
        (
            (2 10 5 2)
			(10 5 7 12)
        );
    }
    walls
    {
        type wall;
        faces
        (
            (8 3 4 9)
            (11 6 7 12)
        );
    }
    front
    {
        type wedge;
        faces
        (
			(0 3 4 1)
            (1 4 5 2)
			(4 6 7 5)
        );
    }
    back
    {
        type wedge;
        faces
        (
			(0 8 9 1)
            (1 9 10 2)
			(9 11 12 10)
        );
    }
    axis
    {
        type empty;
        faces
        (
            (0 1 1 0)
			(1 2 2 1)
        );
    }
);

mergePatchPairs
(
);

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

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

[李东岳]

这个是我画了一下午的wedge网格，先给所需，在这个基础上改吧........