[Jason1114]

Hello All,

I hope you guys are doing great.

I'm a rookie on using OpenFOAM, I'm trying to simulate a propeller/fan by using MRF in simpleFOAM. Turbulence model is SST k-omega.

The case has two region, one is inner region(rotation zone), another is outer region. And the interface between two region's boundary are set as cyclicAMI.

The first problem is that when I start the calculation, there will be a maximum pressure value(about 10^7N) on the surface of the propeller, and my inflow is only 2m/s, rotation speed is 12 rps. I don't know if it is the wall function or the boundary condition is set wrong, so I attach the relevant files for the case.

The second problem is that when I set the rotation speed in MRFProperties as 0 rps, and there is only inflow in the flow field. I found that the fluid seems to deliberately bypass my rotation region, It's just like the fluid doesn't go into the rotation zone. It seems that there is a problem with the AMI setting.

Any help is appreciated

Thank you,
-Jason

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

7
(
    outlet
    {
        type            patch;
        nFaces          3152;
        startFace       10769020;
    }
    wall1
    {
        type            wall;
        nFaces          26746;
        startFace       10772172;
    }
    inlet
    {
        type            patch;
        nFaces          2978;
        startFace       10798918;
    }
    AMI_2
    {
        type            cyclicAMI;
        nFaces          94536;
        startFace       10801896;
		transform		noOrdering;
		neighbourPatch	AMI_1;
    }
    AMI_1
    {
        type            cyclicAMI;
        nFaces          137406;
        startFace       10896432;
		transform		noOrdering;
		neighbourPatch	AMI_2;
    }
    hub1
    {
        type            wall;
        nFaces          8680;
        startFace       11033838;
    }
    prop
    {
        type            wall;
        nFaces          68770;
        startFace       11042518;
    }
)

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

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{

    inlet
    {
        type            fixedValue;
        value           uniform (0.5 0 0);
    }
    
    outlet
    {
        type            zeroGradient;
    }


    wall1
    {
        type            noSlip;
    }
	hub1
    {
        type            noSlip;
    }
    prop
    {
        type            noSlip;
    }



    "(AMI_1|AMI_2)"
    {
        type            cyclicAMI;
    }
}

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

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{

    inlet
	{
        type            zeroGradient;
	}

    
    outlet
	{
        type            fixedValue;
		value			uniform 0;
	}

    "(wall1|prop|hub1)"
    {
        type            zeroGradient;
    }

    "(AMI_1|AMI_2)"
    {
		type			cyclicAMI;
	}
}

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

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

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

internalField   uniform 0.0459375;

boundaryField
{

    inlet
    {
        type            inletOutlet;
        inletValue      uniform 0.0459375;      
    }
    
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0.0459375;
    }

    wall1
    {
        type            inletOutlet;
        inletValue      uniform 0.0459375;
    }
	hub1
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
    prop
    {
        type            kqRWallFunction;
        value           uniform 0;
    }

    "(AMI_1|AMI_2)"
    {
        type            cyclicAMI;
    }
}


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

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 4593.75;

boundaryField
{
    inlet
    {
        type            inletOutlet;
        inletValue      uniform 4593.75; 
    }
    
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 4593.75;
    }

    wall1
    {
        type            inletOutlet;
        inletValue      uniform 4593.7;
    }


    prop
    {
        type            omegaWallFunction;
	      value		        $internalField;
    }
	
	hub1
    {
        type            omegaWallFunction;
	      value		        $internalField;
    }

    "(AMI_1|AMI_2)"
    {
        type            cyclicAMI;
    }

}


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

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

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

internalField   uniform 0;

boundaryField
{

    inlet
    {
        type            calculated;
        value           uniform 0;
    }
    
    outlet
    {
        type            calculated;
        value           uniform 0;
    }

    wall1
    {
        type            calculated;
        value           uniform 0;
    }

    prop
	{
	    // Mandatory entries (unmodifiable)
	    type            nutkRoughWallFunction;
	    Ks              uniform 0.0000384;
	    Cs              uniform 0.75;
	    value uniform 0;
	}
	
	hub1
	{
	    // Mandatory entries (unmodifiable)
	    type            nutkRoughWallFunction;
	    Ks              uniform 0.0000384;
	    Cs              uniform 0.75;
	    value uniform 0;
	}

    "(AMI_1|AMI_2)"
    {
        type            cyclicAMI;
    }
}


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

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

MRF1
{
    cellZone    rotor;
    active      yes;

    // Fixed patches (by default they 'move' with the MRF zone)
    

    origin    (0 0 0);
    axis      (-1 0 0);
    omega     0.0;
    nonRotatingPatches (AMI_1 AMI_2);
}

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

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         cellLimited Gauss linear 0.5;
    grad(U)			cellLimited Gauss linear 1;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss linearUpwindV grad(U);
    div(phi,k)      Gauss upwind;
    div(phi,omega)  Gauss upwind;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited 0.5;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.5;
}

wallDist
{
    method meshWave;
}


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

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0.05;
        smoother        GaussSeidel;
        nCellsInCoarsestLevel 20;
    }

    U
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-07;
        relTol          0.05;
    }

    k
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-07;
        relTol          0.1;
    }

    omega
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-07;
        relTol          0.1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 3;
    pRefCell                 0;
    pRefValue                0;
}

relaxationFactors
{
    fields
    {
        p               0.2;
    }
    equations
    {
        U               0.3;
        k               0.3;
        omega         	0.3;
    }
}


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