[求助]FLUENT模拟爆震燃烧温度和压力过大的问题

luoj303

我这段时间用FLUENT模拟了一下爆震燃烧的问题：
非稳态的，但是计算时温度一直就是上限值：5000K，而压力随着时间的推移从30atm（高能区初始压力）一直增大到48.4（压力上限值），整个计算过程燃烧的激波从左传到右。
我的迭代步长为1e-5步数为150，对solution中的参数进行过调整，但是都得不到稳定的结果，压力要不是一直增加要不是一直减小。
为什么阿？要怎么解决阿？

eric

时间步长太大，要降到1e-6以下，有时候需达到1e-8

luoj303

eric兄，昨天按照你的方法将时间步减为1e-6了，但是温度和压力还是一直以上限值进行推进的，搞不懂为何会一直是上限值，我的初始高能区给的条件是：30atm，3000K.
将时间步长减为1e-8后，由于机子原因，有点慢，至今没有结果，但是算了100个时间步左右，高能区的温度和压力一直还是上限，不知道为何阿？？
上限，为何？

eric

你把高能区的反应物初始浓度改低，初始条件增加部分产物

luoj303

[这个贴子最后由luoj303在 2006/08/10 08:48am 第 2 次编辑]

Eric兄，我回去试验了一下，但是还是不行，我写了一下我的整个计算过程，你帮我
看看到底是什么地方的问题啊。cas文件我也传上来了。
目的：爆震波传播模拟
模型：分离式求解器，非稳态，组分模型（燃烧模型是EDC模型）,辐射模型曾经加入过
燃料: CH4-AIR
初始条件：CH4:0.3,O2:0.6,压力为1atm,温度300K
高能区：CH4:0.3,O2:0.6,压力30atm，温度3000K
        CH4:0.2,O2:0.4,CO2:0.1,H2O:0.1，压力温度一样
时间步:1e-5,1e-6,1e-7
结果：产生一道激波从左向右传，四种情况：1、压力逐渐升高到48.7atm（限制值为50），温度一直在
      最高5000K；2、压力逐渐减小到4，5个压力，最后传到管口或者在中途整个流场混乱了，激波没
      有，温度从3000K左右一直减小到几百K;3、根本产生不了激波，压力和温度计算中会出现最低值，流场中出现0atm和0K；；4、从开始压力和温度一直就是最高
。看过组分的变化，的确有CO2和H2O产生，CH4和O2的量减少。
-------------------传说中的分割线-----------------------------
温度初始图

附：曾经有一次，竟然超过上限值，我并未调过。达到60多

howyie

用的是一阶精度还是二阶精度啊？燃烧非稳态计算精度最少要在二阶以上的，一阶有时候就是捕捉不到激波