随笔分类 - 学习记录
用来记录一下研究生阶段的学习体会，和思考。主要mfix、fluent软件使用，为CFD理论等。
摘要：fortran对于浮点数异常默认是忽略的，如果要开启，有两种方法，一种可以通过调用c函数去开启这部分异常处理，并且需要将函数打包成库，链接到fortran中，另一个是通过编译选项设置。 第一种方法的代码： trapfpe.c /* * ref: * https://gcc.gnu.org/onlin
阅读全文
摘要：2022.5.31更新： 发现这是WSL 1特有的问题，遇到除0的情况，不会报错，只会直接终止程序，并且用文中提到的编译参数： gfortran -g main2.f90 -ffpe-trap=invalid,zero,overflow 也不会给出报错信息，导致无法排查错误。换成WSL 2，编译的时
阅读全文
摘要：paraview可以直接读入csv数据，并且显示为vtk的云图效果，和矢量图效果。 二维矢量图 导入csv数据（数据在文章最后有提供）： 成功导入： 关闭右侧数据，并选择Filters->Alphabetical->Table To Points，也可以直接Ctrl+空格调出搜索框，输入Table
阅读全文
摘要：之前在”DEM轨迹后处理“这篇文章中的第二种方法中介绍过一种方法，但是那种方法只适用于反应器内颗粒数量一定，没有新进入的颗粒的情况。后来在MFiX论坛询问了一下，解决了这个问题。具体方法见： https://mfix.netl.doe.gov/forum/t/how-to-plot-informat
阅读全文
摘要：1. 选中X和两个Y 2. 点击Double Y 3. 关掉gap to Symbol，否则Line+Symbol这种显示方式可能显示不出线条
阅读全文
摘要：基于MFiX-19.2.2 DEM并行程序中的颗粒循环 在DEM并行程序中，每个进程只循环该进程包含的颗粒，并且每个进程还有一层ghost cell，用来存放另一个进程发送过来的颗粒信息。 下面添加一些代码进行实验。在calc_force_dem.f这个文件里对颗粒做一次遍历，且输出颗粒的ID和位置
阅读全文
摘要：在计算颗粒碰撞的时候，需要进行neighbor颗粒的搜寻，只知道大概是基于网格与颗粒绑定的方式，但是具体的实现方式还是比较模糊。搜寻部分代码如下 (mfix-19.2.2)： 可以直接观察到的是，这里用到了两层do循环，外层循环是遍历所有颗粒的ID，LL，这个ID是每个颗粒的全局标记。内层循环用来遍
阅读全文
摘要：其实之前发现了这个问题，但是没有记录，过了好久又忘了。具体问题是，在做增量编译的之前，都会习惯性地删除多余文件再编译，随手就把*.mod和*.inc这类中间文件也删了，结果修改完代码执行 就会提示找不到usr模块中的变量，所以在增量编译的时候，不要删掉这些文件。
阅读全文
摘要：1. 导入数据 2. 选择一列，右键生成Frequent Count 3. 如果要显示相对频率，勾选Relative Frequency 4. 选择第一列和最后一列并生成柱状图 5. 双击生成的图形，调整一下spacing
阅读全文
摘要：更推荐的方法：https://www.cnblogs.com/Jay-CFD/p/12195294.html 2020.6.16更新 方法一：直接在paraview中显示 首先在输出颗粒信息的时候保存global ID： 然后在paraview中导入vtp数据（不要导入pvd），并使用Tempora
阅读全文
摘要：1. 下载Markdown Here源码包 网址：https://github.com/adam-p/markdown-here 2. 创建.xpi后缀文件 将文件夹 中的这几个文件放入同一个文件夹中，并压缩文件，压缩完毕后修改后缀为.xpi 导入.xpi文件到Zotero中 导入完成后重启即可。
阅读全文
摘要：首先分别做好点图和线图在两个graph中，然后选中其中一副图（点图或者线图），然后按下图选择： 如果有多个图，可以在弹出窗口中选择多个sheet。
阅读全文
摘要：基于mfix-19.1.2版本 方法一：直接输出差值网格固相速度 注：这种方式只适用于garg 2012颗粒差值格式在DEM中，默认是无法输出固相速度的网格值的： 但是通过搜索des文件夹下V_s关键字可以看到，实际上网格固相速度是会通过差值计算得到的： 并且在计算DEM颗粒脉动速度的地方也用到了网
阅读全文
摘要：更新： 之前提供的方法虽然没有什么问题，但是需要修改的代码太多，比较麻烦，后来又研究了一下，发现实际上数据自定义输出这一块是有现成接口的。 1）输出网格数据 参考https://www.cnblogs.com/Jay-CFD/p/10786299.html 这篇文章的方法二，利用 ReactionR
阅读全文
摘要：这里安装的时间节点为：2018.10.29，安装的是目前的最新版本CFDEM，支持到与OpenFoam-5.x的耦合。 1. 先安装openfoam:https://openfoam.org/download/source/不算很难，安装这个官网程序走就行，主要需要注意的是源码一定要放在~目录下编译
阅读全文
摘要：先转换为point data 提取空隙率在0.45-1.0之间的网格，为后面提取气泡内网格做准备 把free board部分去掉 然后积分 选择cell data后就得到气泡内所有网格的体积和，如果网格大小与extraction的阈值取得合适，这个体积可以大致代表气泡大小
阅读全文
摘要：具体步骤参考流沙的文章即可，如果python版本较高，可能有个地方需要小小的改动一下： tools文件夹下的compile_lib_edem_coupling.py文件中，导入模块有个地方需要修改 其他的，就按照流沙的教程走就行。
阅读全文
摘要：https://www.youtube.com/watch?v=DQYlnsAfs5s
阅读全文
摘要：在含有DEM反应的dat中，均相反应的block要在DEM反应之前，例如： 如果顺序反了，会报错找不到关键字("has no implicit type")
阅读全文
摘要：参考：https://www.youtube.com/watch?v=UjoSvWdxlTA
阅读全文
