[gxy200992243]

Hi, everyone!

I am now trying to use pisofoam combined with les to simulate wind turbine wake with the help of actuator line method. I have modified the pisoFoam solver. The modified solver works and gives a result, but this result is quite different from what I did with Fluent, in which the actuator line is achieved by UDF. The comparison of the results from openfoam and fluent is attached. The x direction is streamwise direction.

I think the difference may be brought by the scheme I pick in Openfoam.

gradSchemes
{
    default         cellMDLimited leastSquares 1;
    grad(p)         cellMDLimited leastSquares 1;
    grad(U)         cellMDLimited leastSquares 1;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss filteredLinear2 1 0;
    div(phi,k)      Gauss filteredLinear2 1 0;
    div(phi,B)      Gauss filteredLinear2 1 0;
    div(phi,nuTilda) Gauss filteredLinear2 1 0;
    div(B)          Gauss filteredLinear2 1 0;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

Can anyone help me with the fvscheme setting?

Billions of thanks!