CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

interDymFoam skew AMI surfaces

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 19, 2017, 11:26
Default interDymFoam skew AMI surfaces
  #1
New Member
 
Arne
Join Date: Oct 2013
Posts: 16
Rep Power: 11
Arne87 is on a distinguished road
Dear foam experts,

I am trying to run a similar case as the mixerVesselAMI2D tutorial. Differences are that I want to make it a 3D case by adding more cells in the z dimension and the front and back surfaces have to be physical walls (no slip). Furthermore, I want to use pointwise to generate the grid as this is our main meshing software. If I create the mesh using blockMesh, increase the number of cells in the third dimension and change the empty surfaces to walls everything is fine. The case basically runs in the same way as in the tutorial.

However, if I change the geometry slightly by importing the pointwise grid I run into problems. The geometry is visible in the pictures. It is basically a circle with a spared out square in its middle. This square is supposed to rotate. To make this possible I created (as in the tutorial) two AMI patches which distribute the domain in a rotating part and a static part. The center square rotates as intended and on the inside of the domain the rotation is performed as planned. In the close up you can see that in this layer the cells are actually sliding along the AMI surface. However, on the outside of the domain in the front and the back surfaces problems arise. In this layer the cells are not sliding on the AMI, hence solving the rotation by deforming the cells. This leads to a crash of the simulation at some point.

In my opinion my mesh is set up in a similar way as the tutorial case and I don't understand why I get a different behavior. If somebody can help with this issue I would be very glad. I attached all the files which I considered relevant. However, if you need any further information please let me know.

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

8
(
    AMI1
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          1280;
        startFace       141528;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  AMI2;
    }
    AMI2
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          1280;
        startFace       142808;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  AMI1;
    }
    back1
    {
        type            wall;
        inGroups        1(wall);
        nFaces          2540;
        startFace       144088;
    }
    back2
    {
        type            wall;
        inGroups        1(wall);
        nFaces          2432;
        startFace       146628;
    }
    front1
    {
        type            wall;
        inGroups        1(wall);
        nFaces          2540;
        startFace       149060;
    }
    front2
    {
        type            wall;
        inGroups        1(wall);
        nFaces          2432;
        startFace       151600;
    }
    rotor
    {
        type            wall;
        inGroups        1(wall);
        nFaces          1280;
        startFace       154032;
    }
    stator
    {
        type            wall;
        inGroups        1(wall);
        nFaces          1480;
        startFace       155312;
    }
)

// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      topoSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
actions
(
    {
        name AMI;
        type faceSet;
        action new;
        source patchToFace;
        sourceInfo
        {
            name "AMI.*";
        }
    }

    {
        name c1;
        type cellSet;
        action new;
        source cylinderToCell;
        sourceInfo
        {
            p1 (0 0 -100);
            p2 (0 0 100);
        radius 0.999999;
        }
    }

    {
        name rotor;
        type cellZoneSet;
        action new;
        source setToCellZone;
        sourceInfo
        {
            set  c1;
        }
    }


);

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

dynamicFvMesh   solidBodyMotionFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solidBodyMotionFvMeshCoeffs
{
    cellZone        rotor;

    solidBodyMotionFunction  rotatingMotion;
    rotatingMotionCoeffs
    {
        origin        (0 0 0);
        axis          (0 0 1);
        omega         6.2832; // rad/s
    }
}


// ************************************************************************* //
I am looking forward to you answers,

Best

Arne
Attached Images
File Type: png geometry.png (58.1 KB, 13 views)
File Type: png innerdoaminclose.png (8.5 KB, 12 views)
File Type: png innerdomain.png (80.1 KB, 14 views)
File Type: png outsidedeformed.png (69.6 KB, 14 views)
Arne87 is offline   Reply With Quote

Reply

Tags
cyclicami, interdymfoam, mixervessel2d, pointwise

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
AMI speed performance danny123 OpenFOAM 21 October 24, 2020 05:13
Compressor Simulation using rhoPimpleDyMFoam Jetfire OpenFOAM Running, Solving & CFD 107 December 9, 2014 14:38
AMI interDyMFoam for mixer nu problem danny123 OpenFOAM Programming & Development 8 September 6, 2013 03:34
mixerVesselAMI2D's mass is not balancing sharonyue OpenFOAM Running, Solving & CFD 6 June 10, 2013 10:34
Problem with inconsistent patches (AMI, interDyMFoam) jrrygg OpenFOAM Running, Solving & CFD 3 February 1, 2013 04:23


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