@dzw05
k：

dimensions      [0 2 -2 0 0 0 0];
internalField   uniform 0.081;
boundaryField
{ 
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
    }
    frontAndBack
    {
        type            symmetry;
     }
    inter2
    {
        type            cyclicAMI;
     }
    inter1
    {
        type            cyclicAMI;
     }

    blades
    {
        type            kqRWallFunction;
        value           $internalField;
    }
}

nut

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

internalField   uniform 0;

boundaryField
{
        
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type             calculated;
        value           uniform 0;

    }
    frontAndBack
    {
        type            symmetry;
     }
    inter2
    {
        type            cyclicAMI;
     }
    inter1
    {
        type            cyclicAMI;
     }
     blades
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
}

omega

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

internalField   uniform 0.095;

boundaryField
{
    

    outlet
    {
        type            zeroGradient;
    }
    blades
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    inlet
    {
        type            zeroGradient;
    }
        frontAndBack
    {
        type            symmetry;
     }
    inter2
    {
        type            cyclicAMI;
     }
    inter1
    {
        type            cyclicAMI;
     }
}

p

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            fixedValue;
        value           uniform 0;
   }
        frontAndBack
    {
        type            symmetry;
     }
    inter2
    {
        type            cyclicAMI;
     }
    inter1
    {
        type            cyclicAMI;
     }


    blades
    {
        type            zeroGradient;

    }
}

u

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

internalField   uniform (14 0 0);

boundaryField
{
    
    inlet
    {
        type            fixedValue;

        value      uniform (14 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }
        frontAndBack
    {
        type            symmetry;
     }
    inter2
    {
        type            cyclicAMI;
     }
    inter1
    {
        type            cyclicAMI;
     }


     blades
    {
        type            rotatingWallVelocity;

       origin           (0 0 0);

       axis             (0 0 1);

       omega           0;
    }
}

boundary

inlet
    {
        type            patch;
        nFaces          400;
        startFace       6377766;
    }
    outlet
    {
        type            patch;
        nFaces          400;
        startFace       6378166;
    }
    frontAndBack
    {
        type            symmetry;
        inGroups        1(symmetry);
        nFaces          4800;
        startFace       6378566;
    }
    blades
    {
        type            wall;
        inGroups        1(wall);
        nFaces          80334;
        startFace       6383366;
    }
    inter1
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          13994;
        startFace       6463700;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  inter2;
    }
    inter2
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          11572;
        startFace       6477694;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  inter1;
    }

fvSchemes

ddtSchemes
{
    default Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)         Gauss linear;
}

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

laplacianSchemes
{
    default         Gauss linear limited corrected 0.5;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}

fvSolution

solvers
{
    "pcorr.*"
    {
        solver           GAMG;
        tolerance        0.02;
        relTol           0;
        smoother         GaussSeidel;
    }

    p
    {
        $pcorr;
        tolerance        1e-7;
        relTol           0.01;
    }

    pFinal
    {
        $p;
        tolerance        1e-7;
        relTol           0;
    }

    "(U|k|omega)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-06;
        relTol          0.1;
    }

    "(U|k|omega)Final"
    {
        $U;
        tolerance       1e-06;
        relTol          0;
    }

    cellDisplacement
    {
        solver          GAMG;
        tolerance       1e-5;
        relTol          0;
        smoother        GaussSeidel;
    }
}

PIMPLE
{
    correctPhi          yes;
    nOuterCorrectors    2;
    nCorrectors         2;
    nNonOrthogonalCorrectors 5;
    pRefCell  0;
    pRefPoint 0;
    pRefValue 0;
}

relaxationFactors
{
    fields
    {
        p               0.2;
    }
    equations
    {
        "(U|k|omega)"   0.7;
        "(U|k|omega)Final" 0.7;
    }
}

cache
{
    grad(U);
}

CFD中文网

lv1995