灵敏度问题

SuperDJ

问题如图，当压力从21.2bar开始时，21.3bar是OK的，但是我把值压力范围改为21.3bar开始时，21.3bar却是error，这种情况持续到22.3bar才结束，有时候还会出现在中间某个值，或者是模块的值和灵敏度分析值相同时，出现同一个值一个对一个错的情况，不清楚这是怎么回事，麻烦各位指点一下！

                               
登录/注册后可看大图

该贴已经同步到 SuperDJ的微博

adog8207

造成不收敛的原因是：在20~22bar时，计算的时候出现了冷凝。
如图1所示下半部分，在20bar的时候，C3H8的露点为57.0123℃，此时气相的密度d=46.7083kg/m3，如要达到收敛目标d=50kg/m3，则需要降温，降温则会冷凝，必然不能收敛。


如图2所示，在22bar的时候，C3H8的露点为61.5941℃，此时气相的密度d=52.5322kg/m3，如要达到收敛目标d=50kg/m3，则需要升温，从此图的上班部分可以看到，升温至66。7306℃时，密度d=50.0000kg/m3。


建议将收敛精度降低，比如0.01，干嘛要整到0.00001呢。或是换一种物性方法。

jingwan1999

是不是高压状态热力学方程有点问题呀？坐等高人解答

SuperDJ

jingwan1999 发表于 2014-4-7 15:50
是不是高压状态热力学方程有点问题呀？坐等高人解答

试换过其他物性，也存在以上问题，搞的我都无法理解！

SuperDJ

adog8207 发表于 2014-4-28 23:06
造成不收敛的原因是：在20~22bar时，计算的时候出现了冷凝。
如图1所示下半部分，在20bar的时候，C3H8的露 ...

感谢你的帮忙，我对密度进行灵敏度进行分析，得出的结论也是在21.2bar压力下达不到要求的密度，但是我不明白的是为什么会出现一时可以收敛，一时又收敛不了这种情况。

adog8207

SuperDJ 发表于 2014-4-29 22:05
感谢你的帮忙，我对密度进行灵敏度进行分析，得出的结论也是在21.2bar压力下达不到要求的密度，但是我不 ...

能否收敛跟初值有很大的关系，尤其是要求高精度的时候。