DNS of turbulence in a channel flow
 

Now I am using channelOodes (turn off the turbulence model) to run a DNS for a channelflow problem. However, there is slightly difference (5%) compared with other DNS result. The difference might caused by the linear discretization scheme.

Now I am considering using higher order, like fourth order scheme. But how to set the numerical schemes? like interpolationSchemes,gradSchemes,divSchemes, laplacianSchemes, timeScheme, etc.

I just look into the tutorials, boxTurb16. Only the divSchemes is set to cubic(fourth order), Is there any reason for the setting?

Thank you