我很好奇是否有人有关于如何将 Crank-Nicolson（具有近似分解）应用于

$$\nabla \cdot (\nu (\nabla \mathbf{u} + \left(\nabla \mathbf{u}\right)^{T}))$$

我不确定的是如何处理混合导数。目前我正在使用 Crank-Nicolson

$$\nu_{Laminar}\nabla^{2} \mathbf{u}$$

这是直截了当的，因为我可以逐个维度地解决它。