[jamessbriscoe]

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:

/*--------------------------------*- 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)
    );
    }
);

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