CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[Other] Creating spherical domain for solving flow around a sphere

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 19, 2012, 12:36
Default Creating spherical domain for solving flow around a sphere
  #1
New Member
 
James Briscoe
Join Date: Nov 2010
Posts: 3
Rep Power: 15
jamessbriscoe is on a distinguished road
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)
    );
    }
);

// ************************************************************************* //
jamessbriscoe is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] How to correctly mesh a spherical domain with an .stl objet in it ? TerenceA OpenFOAM Meshing & Mesh Conversion 3 September 14, 2017 10:00
Periodic Pressure drop cfd_begin CFX 10 May 25, 2017 08:09
Domain Imbalance HMR CFX 5 October 10, 2016 06:57
Pressure distribution on a wall darazsbence CFX 17 October 6, 2015 11:38
OpenFOAM : flow around a sphere at low Reynold number maxou1993 OpenFOAM Running, Solving & CFD 1 July 7, 2015 06:13


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