# MRF solver

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

September 18, 2013, 14:04
MRF solver
#1
Member

Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 56
Rep Power: 5
Hello everyone，

I'm new to OpenFOAM. I have some experience with simpleFoam. Now I want to move to MRFSimpelFoam. What I'm simulating is an annular ring. In order to reduce the computational cost, I only simulate part of it. The domain is shown in the attachment. There are two walls, one is rotating and the other is stationary. When I select the rotating zone, I put all the cells in the rotaing zone. But the results are not correct. I know there must be something wrong, but I don't know where.

This is how I generate the mesh
convertToMeters 0.001;

vertices
(
(91.89382313 -9.658430006 0) //point 0
(103.7286337 -10.90231872 0) //point 1
(103.7286337 10.90231872 0) //point 2
(91.89382313 9.658430006 0) //point 3
(91.89382313 -9.658430006 0.127) //point 4
(103.7286337 -10.90231872 0.127) //point 5
(103.7286337 10.90231872 0.127) //point 6
(91.89382313 9.658430006 0.127) //point 7
);

blocks
(
hex (0 1 2 3 4 5 6 7) (20 27 10) simpleGrading (1 1 1)
);

edges
(
arc 1 2 (104.3 0 0)
arc 3 0 (92.4 0 0)
arc 5 6 (104.3 0 0.127)
arc 7 4 (92.4 0 0.127)
);

boundary
(
movingWall
{
type wall;
faces ((4 5 6 7));
}
fixedWall
{
type wall;
faces ((0 1 2 3));
}
inlet
{
type patch;
faces ((0 4 7 3));
}
outlet
{
type patch;
faces ((2 6 5 1));
}
side_half0
{
type cyclic;
neighbourPatch side_half1;
faces ((0 1 5 4));
}
side_half1
{
type cyclic;
neighbourPatch side_half0;
faces ((3 7 6 2));
}
);

Then I use topoSet to collect the cells. And it is my first time to use topoSet.
convertToMeters 0.001;

actions
(

// Get all cells
{
name rotor;
type cellSet;
action new;
source cylinderToCell;
sourceInfo
{
p1 (0 0 0); // start point on cylinder axis
p2 (0 0 0.127); // end point on cylinder axis
}
}
{
name rotor;
type cellZoneSet;
action new;
source setToCellZone;
sourceInfo
{
set rotor;
}
}

);
Since the simulation can be run, that means at least the mesh can provide enough information. And last is the U file.

dimensions [ 0 1 -1 0 0 0 0 ];

internalField uniform ( 0 0 0 );

boundaryField
{
inlet
{
}
outlet
{
}
movingWall
{
type fixedValue;
value uniform ( 0 0 0 );
}
fixedWall
{
type rotatingWallVelocity;
origin (0 0 0);
axis (0 0 1);
omega constant -104.72; //rad per sec (1000 rpm)
}
side_half0
{
type cyclic;
}
side_half1
{
type cyclic;
}
}

I know maybe this problem is simple and stupid. But I really need helps. Thank you so much.
Attached Images
 Capture1.JPG (21.9 KB, 29 views)

 September 19, 2013, 02:34 #2 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,195 Blog Entries: 1 Rep Power: 16 well 1- create a test case and upload it here 2- do you think, you need MRF method at all!, it seems only one patch is moving then you can assign it an appropriate BC, where is your MRFzone? __________________ Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog (http://openfoam.blogfa.com/) Training Course on OpenFOAM at (http://www.isme.ir/)

September 19, 2013, 08:24
#3
Member

Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 56
Rep Power: 5
Hi Nimasam，

Thank you for your reply. The test case is in the attachment. The reason why I have to use rotating reference frame is because later the geometry will become complex. There will be some grooves on the wall. As to MRF or SRF, later I will move to interFoam. And there is only MRFinterFoam, there is no SRFInterFoam. So I think I have to use MRF method.
Attached Files
 MRF_AA1.zip (12.0 KB, 2 views)

September 19, 2013, 08:29
#4
Member

Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 56
Rep Power: 5
The test case is in the attachment.

I first use blockMesh and then use topoSet to get the rotating zone.
Attached Files
 MRF_AA1.zip (12.0 KB, 11 views)

 Tags mrfsimplefoam, toposet

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Peter_600 OpenFOAM 4 August 2, 2014 09:52 Smaras FLUENT 2 February 19, 2013 07:58 sivakumar FLUENT 0 February 15, 2013 04:26 danvica OpenFOAM Running, Solving & CFD 16 December 22, 2012 03:09 bearcat CFX 6 April 28, 2008 14:08

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

 Contact Us - CFD Online - Top