对于不同的初始场internalField，使用Simplefoam计算收敛结果不同

[Gallon]

各位老师们好，研究ABL的问题使用SimpleFoam，监测下游距离入口处1/3L处的速度和湍动能。边界条件boundaryField相同，为什么对0文件夹下的internalField进行修改，会对最后收敛的计算结果产生影响。一个internalField设置为0，另外一个internalField设置为入口边界条件的剖面数值

包括后续使用冠层模型中同一组源项去影响结果，稳态计算的收敛结果也是不同的，这是什么原因造成的呢？

[李东岳]

估计还没收敛吧

植入树木冠层么？

[Gallon]

@李东岳 李老师好，是植入从树木冠层改进的一个源项ECD模型，都是从动量方程右端添加源项，但是我用foamMonitor看残差判断已经收敛了吧。不知道为什么对于同一组源项或者不添加源项，对于不同的internalField算出来的结果不一样。我用的singleGraph后处理监测的线

[李东岳]

不知道为什么对于同一组源项或者不添加源项，对于不同的internalField算出来的结果不一样

感觉是你源项加的有点不太好

如果你不加源项，应该没问题。你可以测试下。

[Gallon]

@李东岳
又做了一次测试，测试的内容是采用Openfoam原生的simplefoam求解器（未作任何修改，也没添加源项），唯一的改变在于0文件夹下的U和k以及omega的internalField。如图目标位置处收敛后的监测同一位置处的湍动能还是有一些差距。

还有就是之前添加源项的问题是我添加一组源项可能在一个初始场能够实现很好的自保持效果，但是换了个初始场internalField之后就不能够实现自保持了，怀疑是计算对于源项太敏感了？

[李东岳]

又做了一次测试，测试的内容是采用Openfoam原生的simplefoam求解器（未作任何修改，也没添加源项），唯一的改变在于0文件夹下的U和k以及omega的internalField。如图目标位置处收敛后的监测同一位置处的湍动能还是有一些差距。

如果网格少的话，把这2个算例发上来我看看

[Gallon]

@李东岳 testSimplefoam.zip 谢谢老师，我这边上传了，分别是internalField为边界条件和为0的两个0文件，然后我用的是singleGraph3监测x=2m位置处的速度和湍动能信息。网格数量不大。