[mpl]

Hello!

I'm experimenting with the MRF capability in the simpleFoam solver (turbulence turned off). I'm trying to compare rotatingWallVelocity vs. fvOptions to establish the rotating walls. As a validation and trying to understand how to set up the problem, I'm just using a simple rotating cylinder and I'm trying to see the pressure rise at the outer boundary (like a centrifuge).

My problem is that this doesn't seem to converge. It should be a simple problem that readily converges.

I've attached my blockMeshDict, 0/p and 0/U.

Any suggestions would be appreciated!

Best regards,
Mark


0/p:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}

dimensions [ 0 2 -2 0 0 0 0 ];

internalField uniform 101.325;

boundaryField
{
ID
{
type fixedValue;
value uniform 101.325;
}
OD
{
type zeroGradient;
}
TB
{
type empty;
}
}

0/U:
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}






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

internalField uniform ( 0 0 0 );

boundaryField
{
ID
{
type rotatingWallVelocity;
origin (0 0 0);
axis (0 1 0);
omega 837;
// type fixedValue;
// value uniform (0 0 0);
}
OD
{
type rotatingWallVelocity;
origin ( 0.0 0.0 0.0 );
axis ( 0.0 1.0 0.0 );
omega 837;
// type fixedValue;
// value uniform (0 0 0);
// type surfaceNormalFixedValue;
// refValue uniform 0;
}
TB
{
type empty;
}
}


blockMeshDict:
/*--------------------------------*- 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 0.001;

vertices
(
(25 0 0)
(100 0 0)
(0 0 100)
(0 0 25)
(-25 0 0)
(-100 0 0)
(0 0 -25)
(0 0 -100)
(25 5 0)
(100 5 0)
(0 5 100)
(0 5 25)
(-25 5 0)
(-100 5 0)
(0 5 -25)
(0 5 -100)
);

blocks
(
hex (0 1 9 8 3 2 10 11) rotor (20 1 30) simpleGrading (1 1 1)
hex (3 2 10 11 4 5 13 12) rotor (20 1 30) simpleGrading (1 1 1)
hex (4 5 13 12 6 7 15 14) rotor (20 1 30 ) simpleGrading (1 1 1)
hex (6 7 15 14 0 1 9 8) rotor (20 1 30) simpleGrading (1 1 1)
);

edges
(
arc 1 2 (70.7106781202421 0 70.7106781202421)
arc 2 5 (-70.7106781202421 0 70.7106781202421)
arc 5 7 (-70.7106781202421 0 -70.7106781202421)
arc 7 1 (70.7106781202421 0 -70.7106781202421)
arc 3 0 (17.6776695300605 0 17.6776695300605)
arc 3 4 (-17.6776695300605 0 17.6776695300605)
arc 4 6 (-17.6776695300605 0 -17.6776695300605)
arc 6 0 (17.6776695300605 0 -17.6776695300605)
arc 9 10 (70.7106781202421 5 70.7106781202421)
arc 10 13 (-70.7106781202421 5 70.7106781202421)
arc 13 15 (-70.7106781202421 5 -70.7106781202421)
arc 15 9 (70.7106781202421 5 -70.7106781202421)
arc 8 11 (17.6776695300605 5 17.6776695300605)
arc 11 12 (-17.6776695300605 5 17.6776695300605)
arc 12 14 (-17.6776695300605 5 -17.6776695300605)
arc 14 8 (17.6776695300605 5 -17.6776695300605)
);

boundary
(
TB
{
type empty;
faces
(
(11 10 9 8)
(13 10 11 12)
(13 12 14 15)
(15 14 8 9)
(0 1 2 3)
(3 2 5 4)
(4 5 7 6)
(6 7 1 0)
);
}
OD
{
type patch;
faces
(
(1 9 10 2)
(13 5 2 10)
(15 7 5 13)
(7 15 9 1)
);
}
ID
{
type patch;
faces
(
(8 0 3 11)
(11 3 4 12)
(12 4 6 14)
(14 6 0 8)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //

[crixman]

Hi Mark, did you make any progress?
I am using simpleFoam with fvOptions for MRF (OF 2.2.2.), making the geometry with Salome and meshing with cfMesh - I find it much easier to setup!