CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (http://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   Creating spherical domain for solving flow around a sphere (http://www.cfd-online.com/Forums/openfoam-meshing/104959-creating-spherical-domain-solving-flow-around-sphere.html)

jamessbriscoe July 19, 2012 11:36

Creating spherical domain for solving flow around a sphere
 
I am new to OpenFOAM and am using it to solve underwater implosion problems. I have been successful solving for flow around a spherical bubble of water in a large cubic domain using compressibleInterFoam. After a few oscillations instabilities start to creep in from the domain edge, I believe due to using a cubic domain outside of a spherical implosion.

I've been now trying to create a mesh of one octant of a sphere (with a hole in the center) using a hex block with two vertices collapsed onto other vertices. I am getting a fatal error saying "ill defined primitiveEntry".

Any help you all could give would be greatly appreciated! I haven't yet seen any posts about doing a spherical domain at all so maybe this is not the route to go.

James

Here is my blockMeshDict file:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices
(
    (0 10 0) //0
    (0 0.5 0) //1
    (0.5 0 0) //2
    (10 0 0) //3
    (0 0 10) //4
    (0 0 0.5) //5
);

blocks
(
    hex (0 1 2 3 4 5 5 4) (40 40 40) simpleGrading(1 1 1)
);

edges         
(
  arc 0 3 (7.071067 7.071067 0)
  arc 0 4 (0 7.071067 7.071067)
  arc 1 2 (0.353553 0.353553 0)
  arc 1 5 (0 0.353553 0.353553)
  arc 2 5 (0.353553 0 0.353553)
  arc 3 4 (7.071067 0 7.071067)
);

boundary
(
    sides
    {
        type symmetryPlane;
        faces
        (
            (0 1 5 4)
            (0 3 2 1)
        );
    }
    hole
    {
    type wall;
    faces
    (
        (1 2 5 5)
    );
    }
    out
    {
    type patch;
    faces
    (
        (0 4 4 3)
    );
    }
    frontAndBack
    {
    type empty;
    faces
    (
        (4 5 5 4)
    );
    }
);

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



All times are GMT -4. The time now is 22:53.