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

Problem with cyclicAMI and simpleFoam

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

Like Tree1Likes
  • 1 Post By taxalian

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 25, 2014, 09:29
Default Problem with cyclicAMI and simpleFoam
  #1
New Member
 
Matej Milavec
Join Date: Nov 2013
Posts: 3
Rep Power: 11
Bluewater is on a distinguished road
Dear Foamers,

i am learning how to calculate rotating machines in OpenFOAM. I have a problem with the cyclic boundary condition in the rotor domain. Attached a picture that shows the cross section at the midle of the rotor domain of pressure contour with very high or low (unrealistic) values close to the cycly BC.

Basically I have done the following steps.

1)The mesh I converted to OpenFoam using "Fluent3DMeshToFoam".

2)Then I defined for the cyclic boundary conditions "cyclicAMI" in the boundary file as below

...

p1-rotor
{
type cyclicAMI;
inGroups 1(cyclicAMI);
matchTolerance 0.0001;
transform rotational;
neighbourPatch p2_rotor;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
nFaces 5871;
startFace 6734366;
}
p2_rotor
{
type cyclicAMI;
inGroups 1(cyclicAMI);
matchTolerance 0.0001;
transform rotational;
neighbourPatch p1_rotor;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
nFaces 5871;
startFace 6751983;
}...

3) In the fvOptions the corresponding angular velocity on rotor domain.

MRF1
{
type MRFSource;
active true;
selectionMode cellZone;
cellZone rotor;

MRFSourceCoeffs
{
origin (0 0 0);
axis (0 0 1);
omega 77.5;
}
}

The residuals after 4000 iterations are for pressure just below 10-e3, velocities below 10-e4 and k and omega below 10-e5.

Does anyone knows why my strange results chose to the cyclicAMI boundary condition?

Thanks in advance,
Mat
Attached Images
File Type: jpg p_pMid_p01.jpg (15.6 KB, 126 views)

Last edited by Bluewater; September 30, 2014 at 02:06. Reason: I uploaded a .png file but something went wrong
Bluewater is offline   Reply With Quote

Old   September 26, 2014, 17:53
Default
  #2
Senior Member
 
Join Date: Nov 2010
Posts: 139
Rep Power: 14
taxalian is on a distinguished road
Send a message via Skype™ to taxalian
It seems that you are trying to use OF MRF module, as you are using the fvOptions dict. My suggestion would be to try setting your cyclicAMI patches as non-rotating ones. This might improve strange behaviour near the cyclicAMI boundaries.

Quote:
Originally Posted by Bluewater View Post
Dear Foamers,

i am learning how to calculate rotating machines in OpenFOAM. I have a problem with the cyclic boundary condition in the rotor domain. Attached a picture that shows the cross section at the midle of the rotor domain of pressure contour with very high or low (unrealistic) values close to the cycly BC.

Basically I have done the following steps.

1)The mesh I converted to OpenFoam using "Fluent3DMeshToFoam".

2)Then I defined for the cyclic boundary conditions "cyclicAMI" in the boundary file as below
...
p1-rotor
{
type cyclicAMI;
inGroups 1(cyclicAMI);
matchTolerance 0.0001;
transform rotational;
neighbourPatch p2_rotor;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
nFaces 5871;
startFace 6734366;
}
p2_rotor
{
type cyclicAMI;
inGroups 1(cyclicAMI);
matchTolerance 0.0001;
transform rotational;
neighbourPatch p1_rotor;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
nFaces 5871;
startFace 6751983;
}
...

3) In the fvOptions the corresponding angular velocity on rotor domain.

MRF1
{
type MRFSource;
active true;
selectionMode cellZone;
cellZone rotor;

MRFSourceCoeffs
{
origin (0 0 0);
axis (0 0 1);
omega 77.5;
}
}

The residuals after 4000 iterations are for pressure just below 10-e3, velocities below 10-e4 and k and omega below 10-e5.

Does anyone knows why my strange results chose to the cyclicAMI boundary condition?

Thanks in advance,
Mat
Bluewater likes this.
taxalian is offline   Reply With Quote

Old   September 29, 2014, 08:02
Default
  #3
New Member
 
Matej Milavec
Join Date: Nov 2013
Posts: 3
Rep Power: 11
Bluewater is on a distinguished road
Hi Taxalian,
in fact the cyclic boundary conditions have to rotate with the domain defined as rotational domain. The model consist of one blade (see picture posted in my first post) that spans 72 degrees (so the real whole model has 5 blades). As you suggested I have defined in "fvOptions" file at "noRotatingPatches()" field the two periodic patches (they are defined in "boundary" file as "patches" but with a "createPatchDict" I tranform them to "cyclicAMI", correct me if this is a wrong procedure, because I do not know how to address directly "cyclicAMI" in "noRotatingPatches()").

Code:
 
MRF1
{
    type            MRFSource;
    active          true;
    selectionMode   cellZone;
    cellZone        fluid_rotor;

    MRFSourceCoeffs
    {
        nonRotatingPatches (p1-fluid_rotor p2-fluid_rotor);
        origin      (0 0 0);
        axis        (0 0 1);
        omega       77.5;
    }
}
The "createPatchDict" looks:

Code:
patches
(
    {
        name            p1_rotor;
        patchInfo
        {
            type            cyclicAMI;
            matchTolerance  0.0001;
            neighbourPatch  p2_rotor;
            transform       rotational;
            rotationAxis    (0 0 1);
            rotationCentre  (0 0 0);
        }
        constructFrom patches;
        patches (p1-fluid_rotor);
    }

    {
        name            p2_rotor;
        patchInfo
        {
            type            cyclicAMI;
            matchTolerance  0.0001;
            neighbourPatch  p1_rotor;
            transform       rotational;
            rotationAxis    (0 0 1);
            rotationCentre  (0 0 0);
        }
        constructFrom patches;
        patches (p2-fluid_rotor);
    }
);
I do not understand why for example under 0/* I have to address the cyclicAMI with the names that I use in "constructFrom patches;" of "createPatchDicht", since from the moment I have redefined them as "cyclicAMI" it has already a new name (i. e. in my case p1_rotor and p2_rotor).

If I look at the results now as expected the two cyclicAMI do not rotate, they behave as walls (velocity is 0.) I would expect an easy way to get the cyclicAMI rotate in the proper way), there are any idea? I am doing something wrong?

Quote:
Originally Posted by taxalian View Post
It seems that you are trying to use OF MRF module, as you are using the fvOptions dict. My suggestion would be to try setting your cyclicAMI patches as non-rotating ones. This might improve strange behaviour near the cyclicAMI boundaries.

Last edited by Bluewater; September 29, 2014 at 11:46. Reason: to clarify better content
Bluewater is offline   Reply With Quote

Old   October 1, 2014, 04:54
Default
  #4
New Member
 
Matej Milavec
Join Date: Nov 2013
Posts: 3
Rep Power: 11
Bluewater is on a distinguished road
after some trials I do not see how to get that the periodic bondary conditions (defined as "cyclicAMI") in the rotational domain do what is expected, look the two pictures attached, where it is clear that they rotate in the contrary direction. This clarifies also the strange pressure contours closed to the periodic b.c. I attached in my first post. Does anyone in this forum work with this kind of problems, or "cyclicAMI" inside a "cellZone" that rotates and defined as such in "fvOptions" with "MRFSource"?
Attached Images
File Type: jpg p_pMid_veloVector01.jpg (39.2 KB, 85 views)
File Type: jpg p_pMid_veloVector02.jpg (73.3 KB, 58 views)
Bluewater is offline   Reply With Quote

Reply

Tags
cyclicami, frozen rotor, simplefoam

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
cyclicAMI + simpleFoam crashes in parallel jmf OpenFOAM Running, Solving & CFD 6 June 28, 2017 16:30
Law of the wall: k-omega vs. Spalart-Allmaras (simpleFoam) MrFourier OpenFOAM Running, Solving & CFD 17 June 17, 2014 05:08
CyclicAMI with translation and pressure jump across inlet and outlet santoshgoku OpenFOAM Running, Solving & CFD 5 January 15, 2014 12:37


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