参考之前Yin的那一篇博客
https://www.cnblogs.com/Jay-CFD/p/10621274.html
输出一下颗粒传热量
使用的算例是fluid_bed_dem_2d
- 打开des scalar tracking
前面都和Yin那篇一样，打开enable user scalar tracking
输入2
因为第一个用来测试Y加速度了
第二个输出传热量
然后在output vtk中勾选好
-
更改源代码
源代码为des_thermo_newvalues.f
在mfix-GUI自带的editor里搜索
然后copy到自己的文件夹中
添加如下代码
先在导入模块处导入des_usr_var变量然后在Q_source(: ) 也就是颗粒传热量清零之前复制到des_usr_var中
这里可以看出来，这个变量是每个时间步都清零的，所以必须在这里输出，要不然只能得到零。
-
编译运行
结果如下，可以看到传热过程了