@King_RVM 我解决了我的问题，是乱用了tmp这个内存管理工具，删掉就编译通过了

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

@悬浮的猪 我想请问你解决了你的问题么？我是修改了komegaSST湍流模型，编译后运行自己的算例，结果报错，跟你的错误很类似

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

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

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

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

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

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

@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需要乘以狄拉克函数，实际计算中就是单位对角阵。
目前报错如下:

@Cr 谢谢你的回复！我把跟sij的量都修改成了：：Internal的形式，这几个量编译倒是没有报错，但是出现了新的错误，提示是openfoam里面源代码的错误，如下

这是怎么回事呢？DimensionedField是定义单位的，我定义的标量矢量I1 I2 V1 V2 V3都是无量纲的量，最后的量纲是乘以密度与湍动能决定的

@bestucan 谢谢跟我讨论！divU的定义是标量

S跟I 都是张量，我觉得应该可以相加减啊