CFD Online Discussion Forums - fvscheme setting in Actuator Line simulation

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.

Code:

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!