各位老师,早上好!小白最近在做openfoam燃烧时发现得到的模拟结果过小,比如我在化学反应机理中添加激发态自由基的反应,模拟得到的激发态自由基的结果量级在10^(-9),这在实验中结果是可信的,因为这种自由基不稳定,很快会跃迁到基态,但是从数值角度怎么解释呢?也有一些基态自由基比如CH,它的量级在10^(-5)。在openfoam中很多方程,它的残差到10^(-6),我们就认为它收敛了。10^(-9)量级的值我们怎么才能够保证其数值模拟结果的可信性呢,而不是计算得到的误差?是对一些不大可能发生的反应乘上了一个很小的系数吗?还是经过什么处理?文献中也只是给出了一个结果,但并没有对结果的可信性作出解释,恳请各位老师解惑。
风云5091
帖子
-
openfoam中收敛值过小是否可信,怎么从数值角度解释? -
reactingFoam模拟甲烷空气扩散燃烧,结果不收敛@BlancheDH 好嘞,我去试试,谢谢大佬
-
reactingFoam模拟甲烷空气扩散燃烧,结果不收敛各位大佬们,早上好!小白最近在对一篇关于甲烷空气层流扩散燃烧的文献进行复现,文献是采用fluent进行模拟的,我用的是reactingFoam求解器进行模拟。二维燃烧器结构如下:
文章中说先不开启化学反应,只通甲烷空气,形成一个冷态流场再进行点火燃烧,但并没有说在哪点火的,这是我自己设置的一个点火区域:
模拟中成功实现了点火燃烧,反应是发生了的,但结果总是不收敛,发现是压力方程的起始残差过大(在0.4-0.7之间),也增加了网格数量但是并没有什么效果,想问问各位大佬们这是什么原因导致的?
压力方程残差:
离散方程求解:
-
icoFoam源代码解读各位老师,大佬晚上好!openfoam小白最近在阅读icoFoam源代码时遇到了一点问题,根据李老师对icoFoam求解器的理论推导,它的压力泊松方程如下:
在源代码中首先定义了phiHbyA,应该是HbyA的面通量:
但是,它的压力泊松方程里面的用的是phiHbyA而不是HbyA,这里我不是很能理解:
另外按理说piso算法每算一次U、P,按理来说应该要更新HbyA一次:
但是源代码里面更新的好像是phiHbyA(速度变了,且是个通量)
在我的理解里,这段代码每次进行piso算法时,HbyA似乎没变,毕竟它是从之前的速度方程那里得到的?(对UEqn.A这个函数不太了解)
但是piso算完以后的U用的又是和理论推导的结果一致,用的是HbyA,而不是phiHbyA:
恳请各位老师,大佬解惑,谢谢各位! -
reactingFoam甲烷空气预混燃烧后处理@wangfei9088 好滴,感谢大佬
-
reactingFoam甲烷空气预混燃烧后处理各位老师,大佬们晚上好!openfoam小白最近做了一个关于甲烷空气预混燃烧的数值模拟,想请教一下应该怎么计算层流火焰的燃烧速度以及火焰厚度,另外reactingFoam似乎是一个不可压缩求解器,因此结果文件中并没有密度文件,也想问一下怎样将得到的物质的质量分数分布转换成摩尔浓度?是有相关的指令吗,还是得自己设置一个函数,以下是我速度的结果图:
红圈那块感觉不太合理?
速度边界条件如下:
-
chemkin中的therm.dat文件理解还有一个问题就是我在chemkin反应机理中添加了这个反应:
对于hv这个光子应该在哪里定义呀?我在chem.inp以及therm.date文件中没有定义光子,也不清楚如何定义,但是采用chemkinToFoam时显示了如下错误:
hv改成大写以后就成了这样:
恳请各位老师,大佬解惑,感谢各位了! -
chemkin中的therm.dat文件理解@han1112 好滴,了解了,谢谢大佬,另外激发态的定压比容和熵与激发态相同我是看了一篇文献的假设,然后基于这篇文献做的
-
chemkin中的therm.dat文件理解还有就是对于某种激发态分子,比如CO2*,它和CO2除了焓不同,其他都设置成一样就行了吗?因为之前看到离子和基态分子之间要加上E表明是离子,所以不太确定要不要加上什么符号表示它是激发态,比如下面这个:
找chemkin机理文件也没找到相关的例子,希望各位大佬解惑,谢谢 -
chemkin中的therm.dat文件理解@han1112 原来是这样,了解了,感谢大佬
-
chemkin中的therm.dat文件理解各位老师,大佬们晚上好!openfoam小白最近在阅读chemkin机理中的therm.dat这个文件出来点问题就是下面红圈这一块我不是很明白:
按照chemkinToFoam之后的结果文件来看这一块描述的似乎是分子的摩尔质量以及组成元素,请问各位大佬们它这个到底描述的是什么?在哪本资料里面对这块有所描述呢,谢谢各位! -
reactingFoam求解器修改化学反应机理@LiuWenchao 好嘞,我去看看,谢谢大佬!
-
reactingFoam求解器修改化学反应机理各位老师,大佬们好!openfoam萌新想请问一下在reactingFoam求解器中添加新的反应机理,比如CO+O+M=CO2*+M这个反应是直接在reactingsGRI这个文件里面添加就行了吗?我看了一下reactingsGRI这个文件里面的化学反应可以分为以下这三种反应类型:
我这个反应应该属于上述哪一种呀?还有就是下面的两种反应类型里面分别添加了这两行代码,也不是很清楚它是什么意思
希望各位老师、大佬们解惑!