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