大气边界层浮力项问题

[Vortex]

大气边界层风力机的控制方程为上图（比较前沿的文章基本都采用上式）。我觉得OpenFOAM中的boussinesq浮力项与文章的并不一致，因此我并不敢直接用，有知道其中原理的朋友可以解释一下吗？假如要自己写一个浮力项，<potential temperature>这个带了“<>”的势温按照文献中的意思为"horizontal average"，在求解器应如何实现？是同一水平面的势温取平均还是邻近网格取平均？