DEM反应不收敛问题
之前算DEM反应一直不收敛,后来把计算规模减少到两个固定颗粒,也就是仿照reactive-chemistry和evaperation这两个算例。
目前reactive-chemistry这个算例,把各种输运方程都打开以后,DES_ONEWAY_COUPLED设置为False,就会出现不收敛的情况,而evaperation换成同样的网格和几何,就没有这个情况,于是索性在evaperation上改。
经过尝试发现,反应焓值应该是导致容易发散的一个原因。我设置的反应如下:
如果把异相反应的反应焓设为0,则反应可以顺利进行,不会导致发散(虽然这样并不合理)。
接下来还需要继续寻找解决办法,感觉应该是反应速率编写的不够合理。
下面列表对结果进行记录:
修改内容 | 结果记录 | 备注 |
DH = 0.0 fracDH(1) = 1.0 |
不发散 | |
DH = -52834.0 ! (cal/moles-reacted) |
不发散 |
paraview内检测不到气相反应速率, 但是能正常产生CO2 (注意ReactionRates(IJK, :)是从1到NRR按顺序储存的,DES_RATES(:)和RATES(:) 各自分别是从1开始,所以这里要注意数组的对应关系,比如这里的赋值分别应该是 ReactionRates(IJK, 1) = RATES(1)和ReactionRates(IJK, 2) = DES_RATES(1) |
DH = -52834.0 ! (cal/moles-reacted) fracDH(1) = 0.8 ! HoR assigned to solid phase 1 |
发散 |
最新文章
- 【转】利用反射快速给Model实体赋值
- Eclipse中一个Maven工程的目录结构
- Android4.4访问外部存储
- MongoDB学习笔记~ObjectId主键的设计
- 【C++】 struct结构自动对齐的问题
- 我的iOS之路2
- [AX 2012] Woker user request
- Heap(堆)和stack(栈)有的区别是什么。
- 慕课网-安卓工程师初养成-1-3 使用记事本编写Java程序
- E - 最短的名字
- 【Linux学习笔记】用nc实现两台主机间的文件传输(不需要输密码)
- EasyUI篇のico
- Extjs中Chart利用series的tips属性设置鼠标划过时显示数据
- PHP5的类与对象
- Redis的基本知识
- python timeit模块简单用法
- Java中随机数生成的问题
- JMeter-java.lang.OutOfMemoryError: PermGen space错误
- Replication主要配置项
- hadoop搭建伪分布式集群(centos7+hadoop-3.1.0/2.7.7)