场计算类型的问题

[fangyuanaza]

请教各位老师，我的问题是：读取了一个对称张量，用于计算的时候说计算类型不一致：

以上表达式中的aij_是已经计算出来了的一个无量纲的对称张量，放在0文件夹中，我读取了这个量，也能在屏幕上打印出来，湍流模型编译也能通过，但是调用求解器运算模拟会报错：

为了解决这个问题，采取的措施是：
1、强制转换变量，改成symm(aij_), 还是报错
2、猜测是不是因为这个变量只是对称张量，而其它量是volField的量，也就是说这个量没有带上网格信息？我尝试用volSymmTensorField强制转换也不对，可能也不能这样转换，那么请教各位大神，应该怎么办呢？

[xpqiu]

可能因为初始 mut_ 是0，所以上面公式中除以 mut_ 这里就出现 floating point exception了。

[fangyuanaza]

@xpqiu 谢谢回复我~ 针对这个建议，我给了个很小的量加在mut 上：

dimensionedScalar small_val
    (
        "small_val",
        dimensionSet(0, 2, -1, 0, 0, 0 ,0),
        1e-25
    );

用mut_+small_val 还是会报同样的错误。而且其实我的初始场mut不是0，是DNS算的结果，我是放到RANS计算中，想得到新的湍流频率，因为这个两种模拟的湍流频率不是一个数量级的

[xpqiu]

@fangyuanaza

从上面报错截图看，floating point exception 是因为除法导致的，可能是其他位置有除法运算中出现除以0的情况。

[fangyuanaza]

@xpqiu 我还以为这里是在说类型不一致的问题，我是一行行排查错误的，把其他的改动都注释掉了，现在排查定位到了这一行。

[fangyuanaza]

@xpqiu
问题就出现在这里，我把这个注释掉就能用

[fangyuanaza]

@xpqiu 其他行代码也调试完了，就剩这个地方了，是我读取信息有问题么？

    aij_
    (
	 IOobject
        (
	   "aij",
	    runTime_.timeName(),
	    mesh_,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE		
	),
        mesh_		
    ),

我是这样读取的，我用Info是可以打印出来aij各个分量的值的

[fangyuanaza]

@xpqiu 您是对的，是我犯低级错误了，我定义的小量small_val在这个要修改的式子后面，放到前面就可以了。谢谢

[xpqiu]

@fangyuanaza