CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [blockMesh] Zero thickness plate simulation (https://www.cfd-online.com/Forums/openfoam-meshing/240603-zero-thickness-plate-simulation.html)

adrian.delser January 14, 2022 13:00

Zero thickness plate simulation
 
Hi everyone, I'm trying to do a thin plate simulation, and having looked around the internet I have found references to createBaffles and topoSet in order to achieve this. None of those posts actually explain how to set up and run the simulation nor provide examples I could test (I am still quite new to openfoam). I am using a modified blockmesh file (based off the movingCone tutorial) that fills in the gap left by the movingCone with a new block. The thin plate should go where the slanted wall of the moving cone was before. I am therefore looking for a simulation of a slanted flat plate moving from rest.



It is possible that blockMesh is not the best utility for this and that snappyHexMesh would be more adapted, although I have not used it before. I also intend to move this plate vertically as well as rotate it around a certain point on its length during the simulation, although these are to be worked out later. This is in effect a zero thickness airfoil simulation.


The blockMeshDict file is as follows (the dynamicMeshDict is so far unchanged from the movingCone tutorial):


Code:

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

convertToMeters 0.001;

vertices
(
    (-7.5 0 0)
    (-7 0 0)
    (-3.5 0 0)
    (0 0 0)
    (-7.5 0.75 -0.032745707)
    (-7 0.75 -0.032745707)
    (-3.5 2 -0.087321886)
    (0 2 -0.087321886)
    (-7.5 2.5 -0.10915236)
    (-7 2.5 -0.10915236)
    (-3.5 2.5 -0.10915236)
    (0 2.5 -0.10915236)
    (-7.5 0.75 0.032745707)
    (-7 0.75 0.032745707)
    (-3.5 2 0.087321886)
    (0 2 0.087321886)
    (-7.5 2.5 0.10915236)
    (-7 2.5 0.10915236)
    (-3.5 2.5 0.10915236)
    (0 2.5 0.10915236)
);

blocks
(
    hex (0 1 5 4 0 1 13 12) (15 15 1) simpleGrading (1 1 1)
    hex (2 3 7 6 2 3 15 14) (20 15 1) simpleGrading (2 1 1)
    hex (4 5 9 8 12 13 17 16) (15 15 1) simpleGrading (1 1 1)
    hex (5 6 10 9 13 14 18 17) (50 15 1) simpleGrading (1 1 1)
    hex (6 7 11 10 14 15 19 18) (20 15 1) simpleGrading (2 1 1)
    /**/
    hex (1 2 6 5 1 2 14 13) (50 15 1) simpleGrading (1 1 1)
   
);

edges
(
);

boundary
(
    movingWall
    {
        type wall;
        faces
        (
            /*the line below should become the flat plate ideally*/
            /*(5 6 14 13)*/
        );
    }

    farFieldMoving
    {
        type patch;
        faces
        (
            (9 17 18 10)
        );
    }
    fixedWall
    {
        type wall;
        faces
        (
            (3 7 15 3)
            (7 11 19 15)
        );
    }
    axis
    {
        type empty;
        faces
        (
            (0 1 1 0)
            (2 3 3 2)
            /**/
            (1 2 2 1)
           
        );
    }
    left
    {
        type patch;
        faces
        (
            (0 0 12 4)
            (4 12 16 8)
        );
    }
    farField
    {
        type patch;
        faces
        (
            (8 16 17 9)
            (10 18 19 11)
        );
    }
    back
    {
        type wedge;
        faces
        (
            (0 4 5 1)
            (2 6 7 3)
            (4 8 9 5)
            (5 9 10 6)
            (6 10 11 7)
            /**/
            (1 5 6 2)
           
        );
    }
    front
    {
        type wedge;
        faces
        (
            (0 1 13 12)
            (2 3 15 14)
            (12 13 17 16)
            (13 14 18 17)
            (14 15 19 18)
            /**/
           
            (1 2 14 13)
           
        );
    }
);

mergePatchPairs
(
);

Thanks!


All times are GMT -4. The time now is 03:18.