[Sunxing]

Hi foamers,

I'm simulating a wind farm case. However, the velocity result shows that the wake effect is too strong. In this simulation, I used the simpleFoam solver and the standard k-epsilon model. The U and Fvschemes files are displayed following. Can anyone give me some hints please~

#include        "include/initialConditions"
#include "include/ABLConditions"
dimensions      [0 1 -1 0 0 0 0];
internalField   uniform $flowVelocity;
boundaryField
{
    south
    {
        type            zeroGradient;
    }
    north
    {
        type            atmBoundaryLayerInletVelocity;
        Uref            $Uref;
        Href            $Href;
        n               $windDirection;
        z               $zDirection;
        z0              $z0;
        value           $internalField;
        zGround         $zGround;
     kappa       $kappa;
    }

    terrain
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
     top
    {
        type           slip;
    }
     east
    {
        type            slip;
    }
 
  west
    {
        type            slip;
    }
}

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
    default            steadyState;
}
gradSchemes
{
    default              Gauss linear;
    grad(p)              Gauss linear;
    grad(U)              Gauss linear;
}
divSchemes
{
    default             none;
    div(phi,U)           Gauss upwind;
    div((nuEff*dev(T(grad(U)))))     Gauss linear;
    div(phi,epsilon)     Gauss upwind;
    div(phi,k)           Gauss upwind;
}
laplacianSchemes
{
    default              Gauss linear limited 0.333;
}
interpolationSchemes
{
    default             linear;
}
snGradSchemes
{
    default             limited 0.333;
}
fluxRequired
{
    default             no;
    p;
}