修改下回复格式:
- 满足第2条,大于50币;
- 需要Turbulence modelling for CFD;
- 921703893@qq.com;
修改下回复格式:
@李东岳 李老师,满足大于50,想要Turbulence modelling for CFD这本书,邮箱921703893@qq.com
@King_RVM 我解决了我的问题,是乱用了tmp这个内存管理工具,删掉就编译通过了
@Cr 非常非常感谢你的分析和回复,我的编译通过了,编译通过之后运行算例有新的问题,原因是tmp这个内存管理小工具用得不好,我最新的通过编译的部分修改如下,供大家参考
@悬浮的猪 我想请问你解决了你的问题么?我是修改了komegaSST湍流模型,编译后运行自己的算例,结果报错,跟你的错误很类似
@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需要乘以狄拉克函数,实际计算中就是单位对角阵。
目前报错如下: