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

Mixer,ICEM and AMI

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

Like Tree1Likes
  • 1 Post By sharonyue

Reply
 
LinkBack Thread Tools Display Modes
Old   March 4, 2013, 04:27
Default Mixer,ICEM and AMI
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Hi foamers;

I am using pimpleDyMFoam,the mesh was generated by ICEM.I find my result is a little wired.Anyone can help me?I dont know where is going wrong.

first I copy the fluent.msh into /incompressible/pimpleDyMFoam/mixerVessel2DAMI



1,Stationary Mesh.

The image is attached.

2,Rotate mesh.

The image is attached.

3,Merge mesh

I use ICEM to merge mesh,the image is attached.
AMI1 and AMI2 is actually the same.but there are two boundaries.

4,fluentMeshToFoam -writeZones fluent.msh

After this I get the polyMesh.and there are two parts in cellzones dictionary.


5,Change boundary dictionary.

Code:
5
(
    AMI1
    {
        type            cyclicAMI;
        nFaces          116;
        startFace       6944;
neighbourPatch  AMI2;
        transform       noOrdering;
    }
    IMPELLER
    {
        type            wall;
        nFaces          116;
        startFace       7060;
    }
    AMI2
    {
        type            cyclicAMI;
        nFaces          156;
        startFace       7176;
neighbourPatch  AMI1;
        transform       noOrdering;
    }
    WALLS
    {
        type            wall;
        nFaces          156;
        startFace       7332;
    }
    frontAndBackPlanes
    {
        type            empty;
        nFaces          7216;
        startFace       7488;
    }
)
6,dynamicMeshDict
Code:
dynamicFvMesh   solidBodyMotionFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solidBodyMotionFvMeshCoeffs
{
    cellZone        ROTATE;

    solidBodyMotionFunction  rotatingMotion;
    rotatingMotionCoeffs
    {
        CofG        (0 0 0);
        radialVelocity (0 0 5000); // deg/s
    }
}
7,p
Code:
internalField   uniform 0;

boundaryField
{
    AMI1
    {
        type		cyclicAMI;
	value		$internalField;
    }

    AMI2
    {
       type		cyclicAMI;
	value		$internalField;
    }

    IMPELLER
    {
        type            zeroGradient;
    }

    WALLS
    {
       type            zeroGradient;
    }

   
}
8, U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{


    AMI1
    {
        type		cyclicAMI;
	value		$internalField;
    }

    AMI2
    {
       type		cyclicAMI;
	value		$internalField;
    }

    IMPELLER
    {
        type            movingWallVelocity; // fixedValue;
        value           uniform (0 0 0);
    }

    WALLS
    {
        type            movingWallVelocity; // fixedValue;
        value           uniform (0 0 0);
    }
}
9,pimpleDyMFoam


the result is here.Looks like the velocity in the stationry mesh and rotate mesh is totally different, in the rotate mesh the velocity is big. in the stationary is small.

Looks like the velocity can not penetrate the AMI boundary.

Anyone know how to handle this?
Attached Images
File Type: jpg Stationary Mesh.jpg (54.0 KB, 40 views)
File Type: jpg rotate mesh.jpg (93.5 KB, 36 views)
File Type: jpg MERGE mesh.jpg (82.6 KB, 44 views)
File Type: jpg u.jpg (37.7 KB, 48 views)
kiddmax likes this.
sharonyue is offline   Reply With Quote

Old   March 4, 2013, 22:48
Default
  #2
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
but in interDyMFoam, it looks fine. check it out here:ICEM and AMI unfinished tutorial
sharonyue is offline   Reply With Quote

Old   July 24, 2013, 10:06
Default
  #3
Senior Member
 
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 112
Rep Power: 8
be_inspired is on a distinguished road
Hi,

Did you solve the issue with pimpleDyMFoam?

I am facing a similar problem. I have merged the two zones using Fluent. The AMI interfaces were set as "wall" inside Fluent and then, after conversion to OF with fluent3DMeshtoFoam, I change the type to cyclicAMI.

I have run simpleFoam previously to have a field initialized with no problems in the AMI patches. Then I changed to pimpleDyMFoam and it fails. The max Courant goes to 5.86138e+07. If I start from time 0, the solver runs but after some pimple iterations, the number of iterations required to converge Ux goes up to 1001 and it crashes.

Comparing my case with the propeller case(tutorial) I do not know if it is needed to do all those steps if previously I have already generated the AMI patchs with ICEM and I have already two regions.

Thank you very much
Attached Images
File Type: jpg picture_AMI.JPG (82.7 KB, 17 views)
Attached Files
File Type: txt log.pimple.txt (5.3 KB, 1 views)
be_inspired is offline   Reply With Quote

Old   July 24, 2013, 10:18
Default
  #4
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by be_inspired View Post
Hi,

Did you solve the issue with pimpleDyMFoam?

I am facing a similar problem. I have merged the two zones using Fluent. The AMI interfaces were set as "wall" inside Fluent and then, after conversion to OF with fluent3DMeshtoFoam, I change the type to cyclicAMI.

I have run simpleFoam previously to have a field initialized with no problems in the AMI patches. Then I changed to pimpleDyMFoam and it fails. The max Courant goes to 5.86138e+07. If I start from time 0, the solver runs but after some pimple iterations, the number of iterations required to converge Ux goes up to 1001 and it crashes.

Comparing my case with the propeller case(tutorial) I do not know if it is needed to do all those steps if previously I have already generated the AMI patchs with ICEM and I have already two regions.

Thank you very much
Hi,

If it works well in simpleFoam and crashes with pimpleDyMFoam, it's probably incurred by the AMI mesh, in your pimpleDyFoam, run this solver:

moveDynamicMesh

and see if the mesh is rotating.

Regards.
sharonyue is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem using AMI vinz OpenFOAM Running, Solving & CFD 255 June 12, 2015 10:18
Trouble with simple AMI case Horus OpenFOAM 0 February 26, 2013 10:31
Problem with inconsistent patches (AMI, interDyMFoam) jrrygg OpenFOAM Running, Solving & CFD 3 February 1, 2013 04:23
problems running AMI MUZ OpenFOAM 6 November 20, 2012 07:18
Problem when case with AMI A.A. OpenFOAM Meshing & Mesh Conversion 1 October 4, 2012 11:54


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