![](https://img-blog.csdnimg.cn/20190927151132530.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
后处理
后处理
beidou111
研究图形学流体仿真。个人学习笔记,仅作记录,侵权必删。曾经是搞CFD流化床生物质的。现北航计算机博士生在读。联系请发邮件。私信不看。邮箱:li_cl@foxmail.com
展开
-
ffmpeg图片序列转视频或者GIF
图片序列转换为视频ffmpeg -f image2 -i ~/Videos/anime/anime.%04d.png movie.mp4视频转换为GIFffmpeg -i test.mp4 test.gif图片序列转换为视频ffmpeg -f image2 -r 50 -i ~/Videos/anime/anime.%04d.png movie.gif-r 调整帧率-i表示输入文件-y表示覆盖原有文件...原创 2021-05-18 19:07:20 · 4651 阅读 · 0 评论 -
图片序列转换为AVI视频或GIF的方法(MATLAB)
转换为AVI的代码(MATLAB)framesPath = 'E:\study\biomass\MFiX_Cases\DEMPIC\Analysis postprocess\PRJ DEMPIC_0619_HTandRe\DEMPIC_DEM_0619_SBR1\post\anime\';%图像序列所在路径,同时要保证图像大小相同videoName = 'E:\study\biomass\MFiX_Cases\DEMPIC\Analysis postprocess\PRJ DEMPIC_0619_HTa原创 2021-03-27 17:07:14 · 613 阅读 · 0 评论 -
画出某一高度截面上的颗粒速度分布线图
目标:画出某一高度截面上的颗粒速度分布线图如下图时间:2021-3-11原始算例:MFiX官方算例2D DEM方法1:借鉴Yin博客,先把颗粒速度插值到网格再画图借鉴Yin 博客上将颗粒速度插值到网格上的方法。利用reactionRates这个数组。插值以后再取时均,最后再截面。https://www.cnblogs.com/Jay-CFD/p/10786299.html结果:这种方法是可行的,但是最新版本由于更改了后处理源码会导致错误,需要退回到1931版本如图这种方法需要更改源代码原创 2021-03-12 10:39:27 · 1490 阅读 · 1 评论 -
【paraview】采用tracer和macro宏来记录paraview的操作过程用于重复性工作
录制宏的意义在于重复性工作,相当于记录了手动操作的步骤。注意要点录制tracer的paraview版本不能高于 读取的paraview版本,也就是向下兼容但是不能向上兼容(兼容问题可能出现但不一定出现,尽量采用同一个版本)。录制的脚本语言是python录制的操作可以选默认的第二项any *modified* properties,也可以选第三项,only-user*modified*properties。但是推荐选第三项,因为记录的操作更少,更不容易出问题。基本步骤1. 读取源数据2. 点原创 2021-01-15 21:12:17 · 1511 阅读 · 4 评论 -
用于批量处理气体产量的bash脚本
csvName='GAS_YIELD_AVG_ALL.csv'prjDir="/home/ubuntu1/disk_DATA/Works/MFiX_cases/save_prj_1121/DEMPIC/PRJ DEMPIC_rerun_90s"cd 'newCoeff_save_failed/rerun1028_DEM_newCoeff(30sFailed)'cp '/home/ubuntu1/intcsv.py' .python intcsv.py $csvNamecd "$prjDir"原创 2021-01-29 19:24:00 · 88 阅读 · 0 评论 -
【python后处理】用python对excel数据进行积分,用以得到四种气体产量分数
无异常版本需要预处理,先手动删掉中间出现的非数字行(前三行表头不算)这里说的“异常”是中间的非数字行例如重启了算例的时候会出现下面类似行,需要手动删除import openpyxl# 创建一个工作簿wb = openpyxl.Workbook()# 创建一个test_case的sheet表单# wb.create_sheet('test_case')# 保存为一个xlsx格式的文件# wb.save('cases.xlsx')# 读取excel中的数据# 第一步:打开工作簿wb原创 2021-01-14 23:00:43 · 1133 阅读 · 1 评论 -
【Paraview】决定Paraview后处理速度的因素是CPU
以前一直以为paraview后处理大量数据的时候限制速度的主要因素是硬盘的传输速度。今天做了个测试,发现居然可能是CPU决定的。此测试只是初步的,因为我自己没有其他测试主机了,所以只能测试一下不同硬盘。结果如下测试对象是二维流化床的算例,1800个网格全选后plot over time比较令人惊诧的是,机械硬盘和固态硬盘几乎无差距(甚至有时固态更慢)而CPU占用率一直保持在30%多。Paraview这个软件并行做的非常好,每个线程的负载几乎一样,都是30%多。不亏Paraview(para应该是pa原创 2020-12-18 11:01:00 · 1013 阅读 · 0 评论 -
【ParaView】使用Paraview 的comparitive视图同时显示不同时间步的结果
来源paraview 5.7 user guide chapter 11问题描述一种比较常见的应用场景是:瞬态的流场,要动态地显示其过程,除了每个时间点截屏一次,还可以用本方法快速完成。操作流程以官方算例can.ex2为例从examples文件夹即可加载这个算例就是挤压一块铁皮1. 打开compartive view出现如下视图2. 打开comparative view inspector3 设置参数comparative view inspector中Layout是选原创 2020-11-26 16:28:26 · 1610 阅读 · 1 评论 -
【Paraview】【转】Paraview显示颗粒轨迹
https://www.cnblogs.com/Jay-CFD/p/11066776.html转载 2020-05-01 17:06:53 · 3499 阅读 · 1 评论 -
gnuplot画图初步探究和命令讲解(处理变密度 的reactive chemistry 得到的算例结果)
原理:命令含义plot ‘POST_Xs_D.dat’根据数据文件绘图注释u 1:2using 的简写,表示使用第一列和第二列的数据w lpwith linespoints的简写,表示使用点线图set title ‘mass frac of 4 species’设置图片的标题,如果在同一个图画出多条曲线,要用逗号分割,除了plot其他的都要重写一遍换行 表示下行实际上和上一行是同一行t ‘XD DEM’title 的简写,这里的title其原创 2020-11-05 17:11:32 · 570 阅读 · 0 评论 -
【paraview】vtkwrite函数:利用matlab写入vtk流场输入
参考https://www.cnblogs.com/dalanjing/p/10696804.html目的:利用vtkwrite函数写出符合vtk格式的文件(写出文件保存为*.vtk并可以被paraview读入)MATLAB函数vtkwrite使用方法这里解释一下原博主给出的一个例子前6个参数是必须要给的分别是写出的文件名,数据类型,x坐标,y坐标,z坐标注意这里的X Y Z必须都是三维的数据。即X为354115的张量,Y也是354115的张量,Z也是354115的张量。点的顺序是先X方向上排原创 2020-11-02 09:50:55 · 3918 阅读 · 1 评论 -
paraview显示自己写的流场(入门)
在YouTube上看到Sharcnet HPC这个人的介绍paraview写入流场格式的视频。很简单,基本上不用什么额外的格式。在此尝试了一下。操作1 自己准备一个流场其实就一列,就是数组a,用c++ 写出到二进制文件aaa.raw#include <math.h>#include <stdlib.h>#include <stdio.h>int main(){ const int n=20; const int dim=2*n+1;原创 2020-10-31 21:21:23 · 1137 阅读 · 0 评论 -
MFiX中的时间对应的模块、变量、子程序和输出方法wall_time()和time
代码讲解模块名:time_cpu源文件:time_cpu_mod.f这里面定义了三个函数:wall_time获取当前现实时间 主要调用的是SYSTEM_CLOCKget_cpu_time获取当前CPU时间, 主要调用的是cpu_timeremaining_wall_time_estimate计算距离算完大约还剩余多少时间get_cpu_time先看第二个函数get_cpu_time仅仅就用了下fortran内置函数cpu_time这个函数是在MFiX-20版本新加的其实在19原创 2020-10-29 13:29:12 · 761 阅读 · 0 评论 -
MFiX中后处理输出雷诺数Re的方法
操作方法废话少说,直接上操作方法从\legacy_tests\dem-tests\evaporation\usr_rates_des.f中复制代码如下!----------------------------------------------------------------------!! Function: calc_NRe(M) !!原创 2020-10-25 10:52:28 · 700 阅读 · 0 评论 -
MFiX利用ReactionRates数组把颗粒变量存储到网格上
思路受到这篇Jay的文章的启发,修改了一下https://www.cnblogs.com/Jay-CFD/p/10786299.html借鉴原文的方法二但是这个方法修改了vtk_out.f这个文件,由于19版本和20版本的MFiX在这个文件上改动很大,所以在20上无法运行成功。于是做了如下修改:我发现原博主更改vtk_out.f这个文件主要其实就是用来call des_granular_temperature的,所以何不在别的地方call des_granular_temperature呢?由原创 2020-09-19 09:08:56 · 958 阅读 · 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 · 3001 阅读 · 2 评论 -
【paraview】paraview做时均图并去掉前20s(temporal statistics过滤器)【2020-12-2更新】
Filter名: temperalStatistics原创 2020-09-04 21:59:22 · 1366 阅读 · 0 评论 -
paraview输出高清晰度带速度箭头颗粒分布图
打开.pvd文件导入颗粒(这个文件本身只包含时间信息,用来和.vtp文件关联的)用速度着色(着色选jet,表示形式选point gausssian),并添加时间标注(control+空格搜索 annote time filter ), 如图调整下时间标注的字体等等为颗粒加下速度箭头使用glyph稍微调整下参数,如图输出file->save screenshot随便给个名字,一般学术出版,格式选tiff或者jpeg像素给大一些看看结果...原创 2020-08-24 22:03:05 · 6523 阅读 · 3 评论 -
【MFiX后处理】MFiX中使用des_usr_var输出颗粒的传热量(2020-11-27更新变更输出变量名方法)
参考之前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自带的ed原创 2020-08-21 00:54:16 · 1037 阅读 · 0 评论 -
【paraview进行MFiX后处理】如何计算某一气体组分的质量流量/气体产量
MFiX版本:20.1.2paraview版本:5.8目的:在一个反应容器中,发生反应并产生H2,我们需要算出H2的质量流量,并且测定出气体产量。当前MFiX只能输出H2气体质量分数,但是没有办法输出总质量或者质量流量,所以问题归结于求质量流量,并且进行时间积分。求总产量,相当于求对某一切面流过的质量流量积分。导入流场。做切面。点Y normal可以直接把切面与Y轴垂直如果没显示切面,勾选show plane在 origin中输入切面所在位置对切面积分在filter &g原创 2020-05-21 17:49:09 · 2637 阅读 · 0 评论 -
【Barracuda】后处理使用plot manager画线图的方法总结(包括残差,探针,通过某个面的组分质量流量等)(持续更新)
使用plot manager 可以画出二维线图包括但不限于各种残差(残差图其实就藏在这里),某一点的压力、温度值随时间的变化(相当于fluent或mfix里的probe探针), 还有气体通过某一平面(Flux Plane)的各组分质量分数、质量流量等。今天先开个头、残差图点击plot manager点add plot然后填好各种标题(subtitle非必须)name是该图在barracuda里面的一个名称,类似给BC起名字,比如我们给residualtitle就是显示在图里的 我们也给re原创 2020-05-11 00:08:44 · 1706 阅读 · 3 评论 -
【Barracuda】Barracuda用GMV后处理的使用方法总结(持续更新2020-5-22)
读取文件或者改变所显示的时间read GMV file -> same simulation选择文件。因为我是每隔0.1s保存一次的,所以100就代表第10s的流场读取流场display->cells勾选Faces, Shaded和Refine点击Color by->node fiield下的 new field-> 选择你要展示的流场文件(如H2O....原创 2020-05-08 16:41:31 · 3476 阅读 · 6 评论