dev(grad(U)) 与 dev2(grad(U)) 的区别

OpenFOAM
1

在 foam-extend-3.1 中，interPhaseChangeFoam 的 U 方程里面

    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
      - fvm::laplacian(muEff, U)
      - (fvc::grad(U) & fvc::grad(muEff))
    //- fvc::div(muEff*(fvc::interpolate(dev2(fvc::grad(U))) & mesh.Sf()))
    );

而在有的求解器中，比如，MRFInterFoam 的 U 方程

    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      - fvm::laplacian(muEff, U)
      - (fvc::grad(U) & fvc::grad(muEff))
    //- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
    );

虽然这两行都被隐去，但是这两行代表什么意思呢？

在 cavitatingFoam 中，

    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(phi, U)
      - fvm::laplacian(muEff, U)
    //- (fvc::grad(U) & fvc::grad(muf))
      - fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
    );

这一行又启用了。

有什么区别或用意吗？