修改湍流模型，定义表达式的问题

[King_RVM]

@Cr 是的，我定义I1 I2的时候调用了转置函数，但是我现在发现实际上我不需要，于是删去了转置计算。关于量纲的问题，你说的很对，虽然sij omegeij、标量I1 I2和张量V1 V2 V3没有量纲，但是定义过程中S 与 Omega有量纲，为了减少量纲的问题，我定义如下

sij 与 omegaij的无量纲化是通过除以omega_实现的。现在的报错是与定义的张量I 有关，在sij中，场离散的项symm(tgradU())不能直接与单位对角矩阵I相加减，该怎样把I转换成一样的类型呢？我想通过volTensorField定义I，但是报错。我定义单位对角阵是因为定义sij中1/3*partU/partx需要乘以狄拉克函数，实际计算中就是单位对角阵。
目前报错如下:

[King_RVM]

@Cr 我发现问题来自于单位矩阵的相乘，当我用更简洁的形式去表达sij来避免单位对角矩阵I的表达，sij没有报错，但是V3关于I 的部分报错了。其实我想添加的表达式是

其中sij为除以omega_来做到无量纲化，而图中omegeij为
除以omega_来做到无因次化
我现在修改后的定义为

V3中的I不可避免，因为需要将那一项从标量转为张量，现在的报错如下

[Cr]

@King_RVM 这个错误还是类型不统一，sij是tmp<volTensorField>类型， I1是一个volTensorField::Internal或DimensionedField<tensor, volMesh>类型，OpenFOAM没有定义两种类型的减法，我建议你将所有的vol##Type##Field::Internal都换成vol##Type##Field试试。

[Cr]

@King_RVM 这个编译错误是OpenFOAM本身源码上的bug，需要对DimensionField的源代码进行改动才可以使编译通过。OpenFOAM-v2006 bug Issues #1868

[King_RVM]

@Cr 我尝试过改成一样的类型，所有都改成vol##Type##Field，但是这样sij会报错

在这里 dev(symm(tgradU()()))是张量除以一个标量还是张量，还是报错。
dev(symm(tgradU()()))的灵感来自于李老师的东岳流体网站http://www.dyfluid.com/docs/tensor.html

[King_RVM]

@Cr 非常感谢这部分资料的提供，我没想过源代码会有问题

[Cr]

@King_RVM dev(symm(tgradU()()))改为 dev(symm(tgradU()))，同时将sij的类型声明为tmp<volSymmTensorField>。
原因：tgradU()为tmp<T>的operator()重载返回T类型的引用，在这里为GeometricField类型，GeometricField的operator()的重载返回GeometricField::Internal类型，即dimensionedField<Type, GeoMesh>，所以类型不匹配，场运算报错。

[King_RVM]

@Cr 非常非常感谢你的分析和回复，我的编译通过了，编译通过之后运行算例有新的问题，原因是tmp这个内存管理小工具用得不好，我最新的通过编译的部分修改如下，供大家参考

[Chen_hao]

你好，请问你修改湍流模型的雷诺应力项（跟你一样修改boussnesiq 假设）是在哪修改的？是在相应的komegaSST.C中，还是我看他们说的linearViscousStress.C？还有一直有个疑惑，相应的动量方程中这个雷诺应力项用修改吗？ 不好意思打扰你了，希望能不吝赐教。

[李东岳]

@chen_hao

1. linearViscousStress.C

2. UEqn不用动