关于几个运算符运算后的量纲疑问

[Chen_hao]

定义如下：

     
    tensor I(1,0,0,0,1,0,0,0,1); 
    volSymmTensorField Sij(twoSymm(fvc::grad(this->U_)));
    volTensorField Wij(2.0*skew(fvc::grad(this->U_)));

运算如下：

     Sij&Sij  - (1.0/3.0)*magSqr(Wij)*I 

Sij应该是[0 0 -1 0 0 0 0]，为啥 &内积 后还是[0 0 -1 0 0 0 0]？
而 magSqr(Wij) 模的平方后乘以I为啥 [0 0 -2 0 0 0 0] ？

[Chen_hao]

@chen_hao 在 关于几个运算符运算后的量纲疑问 中说：

定义如下：

     
    tensor I(1,0,0,0,1,0,0,0,1); 
    volSymmTensorField Sij(twoSymm(fvc::grad(this->U_)));
    volTensorField Wij(2.0*skew(fvc::grad(this->U_)));

运算如下：

     Sij&Sij  - (1.0/3.0)*magSqr(Wij)*I 

Sij应该是[0 0 -1 0 0 0 0]，为啥 &内积 后还是[0 0 -1 0 0 0 0]？
而 magSqr(Wij) 模的平方后乘以I为啥 [0 0 -2 0 0 0 0] ？

显示错误如下，我通过排查也确实是这个地方，Sij & Sij 量纲显示为[0 0 -1 0 0 0 0 ]，解释不通 难道要强行在乘一个量纲数？

--> FOAM FATAL ERROR: 
LHS and RHS of - have different dimensions
     dimensions : [0 0 -1 0 0 0 0] - [0 0 -2 0 0 0 0]

[Chen_hao]

@chen_hao
解决了，改成了无量纲，就避免出问题：

 volSymmTensorField Sij(twoSymm(tgradU())/omega_);
 volTensorField Wij(2.0*skew(tgradU())/omega_);