读MFiX源代码
读MFiX源代码
beidou111
研究图形学流体仿真。个人学习笔记,仅作记录,侵权必删。曾经是搞CFD流化床生物质的。现北航计算机博士生在读。联系请发邮件。私信不看。邮箱:li_cl@foxmail.com
展开
-
【读MFiX源代码】0 读代码准备工作:如何利用gdb调试mfix程序?并且用Visual studio code将debug可视化(使用vscode配置linux下Fortran编译调试环境)
文章目录前提准备方法1. 编译你的case2. 把你编译的程序装入gdb3. gdb中运行和查看源代码常见错误总结进阶:用VS code 可视化gdb1. 配置launch.json2 配置task.json前提准备安装了gdb的linux/windows没有的话ubantu下apt-get 下载一个,一般都自带window要安装mingw64(注意不要弄错了64位和32位不能通用)本...原创 2020-04-21 01:09:37 · 1274 阅读 · 0 评论 -
MFiX-19.3.1的调用关系图(用doxygen和grapviz生成的)
对model文件夹生成的效果不是很好,凑合看吧MFiX-19.3.1版本使用方法解压后打开html文件夹下面的index.html下载链接:https://pan.baidu.com/s/15YP9QLBJ-FDtUjM4AlGXeQ提取码:mdxu复制这段内容后打开百度网盘手机App,操作更方便哦例子效果一览(mfix.f中的调用关系)...原创 2020-11-04 17:19:57 · 320 阅读 · 0 评论 -
【读MFiX源代码】MFiX中四种传热方式全面详解(对流、导热、辐射、反应热)并且输出以供后处理(2020-12-15更新)
文章目录1 目标2 代码结构2.1 des_time_step外层总控制2.2 颗粒导热2.2.1 颗粒-颗粒导热:calc_force_dem.f和des_thermo_cond_mod.f2.2.2 颗粒-墙壁导热:calc_dem_thermo_with_wall_stl.f和des_thermo_cond_mod.f2.3 除导热外其他三种传热的控制:calc_thermo_des2.3.1 CONV_GS_DES1 对流传热量2.3.2 DES_RADIATION计算辐射量2.3.3 RXNS_G原创 2020-09-13 22:51:19 · 3286 阅读 · 2 评论 -
MFiX存储ReactionRates的模块
正文位置:/model/rxns_mod.f第五行! reaction rates DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: ReactionRates在usr_rates.f和usr_rates_des.f中导入rxns模块即可直接使用ReactionRates数组USE rxns后记这个模块rxns还有许多其他有用的变量。该模块基本作用就是变量定义。比如把组分定义成一个结构体,把反应定义成结构体。以后再挖掘这里列原创 2020-08-19 23:33:41 · 398 阅读 · 0 评论 -
MFiX-PIC中计算颗粒固相应力的相关代码
源文件名:calc_ps_pic.f90源文件位置\model\des\pic\calc_ps_pic.f9081-90行为核心,如下 ! Particle stress :: Snider (Eq 33) ! tau = (PsubS*ThetasubP^Beta)/max(ThetasubCP-ThetasubP,small#*Theatsubf) DO IJK = IJKSTART3, IJKEND3 IF(VALID_PS_AT(IJK)) THEN Ps(原创 2020-08-19 14:15:25 · 500 阅读 · 1 评论 -
【MFiX源代码】sutherland公式计算粘度
MFiX采用sutherland公式计算粘度只能适用于空气,因为源代码里只给了空气源代码位置\mfix-19.3.1\model\calc_mu_g.f其中F2O3为2/3最主要的语句是136行MU_G(IJK) = to_SI*1.7D-4 * & (T_G(IJK)/273.0D0)**1.5D0 * (383.D0/(T_G(IJK)+110.D0))动力粘度以273K 时 1.7E-4为基础,依照该网格温度变化μg(IJK)=1.7×10−4×(Tg273)1.5×(3原创 2020-06-12 17:22:14 · 5556 阅读 · 1 评论 -
【MFiX源代码】MFiX中的颗粒-流体换热关系式为Ranz-Marshall关系式
目前为止应该只有这一个关系式源代码所在位置:calc_gamma_des.f中第111行尽管源代码前面给出了 select case, 貌似是要给出多种换热关系式,然而目前代码里植入了这一种参考J. Musser博士毕业论文的P72 (pdf实际页面是第94页)...原创 2020-05-18 16:49:47 · 1006 阅读 · 0 评论 -
【读MFiX源代码】4 run_dem(des_time_march)计算颗粒相(主要步骤DEM_TIME_STEP)(源文件在des_time_march.f)
位置被谁调用(入口)run_mfix(dt_loop)(mfix.f)mfix.f文件中第146行,注意观察输出输出仅为这三行run_dem进入run_dem(仍然在mfix.f)可以发现主要步骤有三步:des_time_initdes_time_stepdes_time_end运行完des_time_init会输出第一行运行完des_time_step会输出第二行...原创 2020-05-08 17:22:44 · 1206 阅读 · 3 评论 -
MFiX中的颗粒辐射传热是如何计算的?
一直以为流化床中颗粒辐射传热会极其复杂,谁能想到MFiX中采用的是一个极其简陋的辐射传热模型。主要内容简化假设所有颗粒都是漫灰体某个颗粒的辐射传热简化为颗粒对无限大空间的辐射该无限大空间(环境)的温度为该颗粒周围所有颗粒温度的算术平均值计算公式其中TenvT_{env}Tenv 是环境温度上标或者下标iii代表该颗粒Np,ΩN_{p,\Omega}Np,Ω...原创 2020-04-30 19:31:49 · 798 阅读 · 1 评论 -
【读MFiX源代码】1程序入口mfix.f文件
2020/4/11 第一天本文以19.3.1版本为例主要的源代码都在mfix的源代码文件夹(一般在anaconda安装位置的envs文件夹下能找到)下的model文件夹我的mfix安装在D:\App\Anaconda3\envs程序的入口是mfix.f找到PROGRAM MFIX这就是程序的入口END PROGRAM就是程序的终点在这个源代码里,定义了多个subroutine,有的...原创 2020-04-11 14:47:30 · 949 阅读 · 0 评论 -
【读MFiX源代码】1.5 GET_DAT读取mfix.dat并且检查BC IC 反应配平 热物性(源文件在main.f)
所在位置调用位置mfix.f第127行子程序所在位置/model/main.f第469行-557行所在module: module mainSUBROUTINE GET_DATA(MFIX_DAT)作用! Read input data, check data, do computations for IC and BC locations! and flows, and s...原创 2020-04-23 15:20:57 · 652 阅读 · 0 评论 -
【读MFiX源代码】2 读time march loop(dt_loop)(仍在mfix.f)
目录上次总结Time march looprun_fluidrun_dem上次总结mfix.f是程序入口mfix.f的入口在PROGRAM MFIX 而program mfix的主要作用是调用了run_mfixrun_mfix的主要步骤CALL GET_DATA(MFIX_DAT) 读取deck文件,包括设定几何,BC,IC位置CALL INITIALIZE(MFIX_DAT) ...原创 2020-04-13 20:57:55 · 553 阅读 · 0 评论 -
【读MFiX源代码】3 run_fluid计算流体相 (主要步骤DO_ITERATION) (源文件iterate.f)
MODULE ITERATE文件层次ITERATE_INIT(118行-200行)DO_ITERATION合理的创建标题,有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能,丰富你的文章UML 图表FLowchart...原创 2020-04-24 22:29:44 · 1152 阅读 · 1 评论 -
【读MFiX源代码】番外:calc_thermo_des.f
目录概述主体内容(36-66行)概述该文件定义了一个模块MODULE CALC_THERMO_DES_MOD这个模块的内部只有一个子程序SUBROUTINE CALC_THERMO_DESMODULE CALC_THERMO_DES_MODSUBROUTINE CALC_THERMO_DES先看下注释! Module name: CALC_THERMO_DES ...原创 2020-04-20 17:17:51 · 394 阅读 · 0 评论