【MFix】用particle_input.dat自定义指定初始颗粒位置的要点

最新推荐文章于 2021-05-20 22:06:03 发布
最新推荐文章于 2021-05-20 22:06:03 发布
阅读量571 收藏 3
点赞数
分类专栏： # MFiX
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43940314/article/details/109790321
版权

文章目录

旧版particle_input.dat（MFiX-20.3以前，新版本也兼容）

particle_input.dat中的信息

根据官方指南描述
在这里插入图片描述从左到右每列的含义是（对于2D)
x坐标
y坐标
半径（注意是半径，GUI里给的是直径）
密度
速度u分量
速度v分量

单位都是SI

要点

要点：

1 particle_input.dat要和GUI solids里指定的密度和直径相对应

MFiX就是靠这个分辨颗粒的相的。

例如我给了两个颗粒：一个是生物质，一个是沙子。生物质密度562，沙子密度2640

在这里插入图片描述

在particle_input.dat文件中本身是没法区分哪个是生物质，哪个是沙子的。

但是根据GUI中给定的密度和直径，MFiX可以在particle_input.dat找到与其相对应的值

例如
在这里插入图片描述这样就可以分辨哪个是沙子，哪个是生物质了。

注意，如果没有完全对应密度直径，运行时会报错。

报错如下
在这里插入图片描述

PS:我想大概这也是为什么GUI里会在上面显示Diamter 和 Density两项。

2 要在solid中指定颗粒数目，这个数目可以小于particle_input中的数目。

GUI中指定颗粒数目
在这里插入图片描述或mfix.dat中关键字指定颗粒数目
在这里插入图片描述

3 要关闭自动颗粒生成

4 要在IC中指定组分分数

尤其注意，你的颗粒在哪个IC区域，它就会根据哪个IC区域的组分分数来指定。

5 注意颗粒位置 不能让在particle_input中指定的颗粒重叠

否则会运行时报错

报错如下
在这里插入图片描述

待续。。。。

新版particle_input.dat（MFiX-20.3以后）

beidou111
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
particle_filter_demo.rar_DEMO_particle filter_particle_filter_de
09-24
particle_filter_demo
MFIX使用指南
11-12
CFD模拟仿真的开源软件MFIX使用方法官方英文指南，主要应用于DEM模型模拟
【MFiX】使颗粒静止（常常用于测试）
beidou111的博客
11-13 598
说明 在测试过程中为了防止颗粒飞出，有时需要保持颗粒静止 或者做固定床的模拟的时候。 操作 参考reactive-chemistry算例 添加如下文件即可 usr1_des.f !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv! ! ! ! Module name: URS1_
mfix输出自定义数据
weixin_30449453的博客
04-04 388
更新： 之前提供的方法虽然没有什么问题，但是需要修改的代码太多，比较麻烦，后来又研究了一下，发现实际上数据自定义输出这一块是有现成接口的。 1）输出网格数据 参考https://www.cnblogs.com/Jay-CFD/p/10786299.html这篇文章的方法二，利用ReactionRates这个变量来储存和输出。 2）输出颗粒数据 利用软件提供的user scalar ...
MFiX中后处理输出雷诺数Re的方法
beidou111的博客
10-25 716
操作方法 废话少说，直接上操作方法 从\legacy_tests\dem-tests\evaporation\usr_rates_des.f中复制代码如下 !----------------------------------------------------------------------! ! Function: calc_NRe(M) ! !
画出某一高度截面上的颗粒速度分布线图
beidou111的博客
03-12 1551
目标：画出某一高度截面上的颗粒速度分布线图 如下图 时间：2021-3-11 原始算例：MFiX官方算例2D DEM 方法1：借鉴Yin博客，先把颗粒速度插值到网格再画图 借鉴Yin 博客上将颗粒速度插值到网格上的方法。利用reactionRates这个数组。插值以后再取时均，最后再截面。 https://www.cnblogs.com/Jay-CFD/p/10786299.html 结果：这种方法是可行的，但是最新版本由于更改了后处理源码会导致错误，需要退回到1931版本 如图 这种方法需要更改源代码
Particle_For_UnlineOneDiv.zip_particle filter_粒子滤波
07-15
粒子滤波（Particle Filter）是一种基于蒙特卡洛模拟的非线性、非高斯状态估计方法，广泛应用于导航、视觉跟踪、信号处理等领域。在标题提到的"Particle_For_UnlineOneDiv.zip"压缩包中，包含了一个名为"Particle_...
particle_builder.zip
09-20
《cocos粒子特效：particle_builder深度解析》 在游戏开发领域，粒子特效是不可或缺的一部分，它们为游戏场景增添了视觉震撼力和生动性。Cocos2dx作为一款强大的跨平台游戏开发框架，内置了丰富的粒子系统，让...
Particle_GS.zip_高斯粒子滤波
07-14
高斯粒子滤波（Gaussian Mixture Particle Filter, GM-PF）是一种在滤波理论中的高级算法，主要用于解决非线性、非高斯状态估计问题。在这个“Particle_GS.zip”压缩包中，包含了一个名为“Particle_GS.m”的Matlab...
JohannesPfeifer-Particle_Filtering.zip
最新发布
08-18
粒子滤波（Particle Filtering）是一种在非线性、非高斯状态估计问题中广泛使用的算法，由Gordon、Salmond和Smith于1993年提出，它属于贝叶斯滤波框架下的一种重要方法。Johannes Pfeifer的资料集"JohannesPfeifer-...
【MFiX】MFiX中变量总结 【持续更新 2020-12-21】
beidou111的博客
05-22 1465
名称 变量名 模块 气体密度 RO_g(IJK) physprop 和 fldvar 想到就更新
学生成绩管理系统纯C（仿木鱼不会水）
he___H的博客
06-06 1599
在这里插入代码片 ```#include <stdio.h> #include <string.h> #include <stdlib.h> #include<conio.h > #define MAX 3 struct score{ char id[10];//学号 char name[10]; int chinese; int math; int english; int sum; int rank;//排
ccs读取dat文件c语言程序,详解CCS中的.dat文件
weixin_42298032的博客
05-20 2970
CCS支持的.dat文件的格式为：文件头为定数　　 数据格式　　起始地址　　页类型　　数据块大小1651 1 80000000 0 10固定标识 数据格式 基地址页类型 长度数据格式：1-十六进制 2-十进制 3-十进制长整型 4-十进制浮点型...
【转】MFiX输出自定义颗粒或网格信息
beidou111的博客
08-19 654
转自：https://www.cnblogs.com/Jay-CFD/p/10621274.html JayYin 本文仅作自己方便查看用，建议查看原文 输出颗粒信息 总结一下 在solid选项卡 DEM页面下打开user scalar tracking并输入变量个数 在output中勾选DES user scalar 1 修改源代码（先把待改源码拷贝到自己的文件夹）。如输出颗粒的Y方向加速度，需要在cfnewvalues.ffnewvalues.f中添加 DES_USR_VAR(1,:) = FC(:
纠正mfix19.1 user guide里的一个错误
beidou111的博客
09-03 342
安装mfix19.1的时候死活安装不上conda环境 最后发现是这个命令错误了： conda env create -n mfix-19.1 mfix/mfix-win64 如果输入这个命令，会出现下面的报错 (base) C:\Users\Administrator>conda env create -n mfix-19.1 mfix/mfix-win64 SpecNotFound: m...
【MFiX后处理】MFiX中使用des_usr_var输出颗粒的传热量（2020-11-27更新变更输出变量名方法）
beidou111的博客
08-21 1095
参考之前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
Tensorflow学习笔记（8）——input_data.py解析
热门推荐
lhanchao的博客
05-03 3万+
这里学习一下前面用到的读取mnist数据库文件的代码。其实并没有用到Tensorlfow的东西，但是读取数据库文件是使用Tensorflow编程实现功能的基础，因此归到Tensorflow的学习笔记中。 这里需要注意的主要有以下几点：
MFiX利用ReactionRates数组把颗粒变量存储到网格上
beidou111的博客
09-19 997
思路 受到这篇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呢？ 由
loongson2f 64位内核编译
涓涓溪流的专栏
07-13 1472
在文件/etc/apt/sources.list 添加源 http://dev.lemote.com/debian lenny main contrib non­free http://dev.lemote.com/debian­loongson loongson main 安装编译工具: apt­get install gcc   下载内核源码 git clone git://de
function [GlobalBest, particle] = initialize(problems, params) nVar = problems.nVar; VarMin = problems.VarMin; VarMax = problems.VarMax; nPop = params.nPop; VarSize = [1 nVar]; empty_particle.Position = []; empty_particle.Velocity = []; empty_particle.Cost = []; empty_particle.Best.Position = []; empty_particle.Best.Cost = []; particle = repmat(empty_particle, nPop, 1); GlobalBest.Cost = inf; for i=1:nPop particle(i).Position = unifrnd(VarMin, VarMax, VarSize); particle(i).Velocity = zeros(VarSize); particle(i).Cost = CostFunction(particle(i).Position); particle(i).Best.Position = particle(i).Position; particle(i).Best.Cost = particle(i).Cost; if particle(i).Best.Cost < GlobalBest.Cost GlobalBest = particle(i).Best; end end end帮我解释一下这段程序
05-11
这段程序是用来初始化粒子群优化算法的粒子的位置、速度、成本等信息的函数。下面对每一行代码进行解释： 1. `function [GlobalBest, particle] = initialize(problems, params)`：定义了一个名为 `initialize` 的函数，输入参数为 `problems` 和 `params`，输出参数为 `GlobalBest` 和 `particle`。 2. `nVar = problems.nVar;`：获取问题的变量个数。 3. `VarMin = problems.VarMin;`：获取问题中每个变量的最小值。 4. `VarMax = problems.VarMax;`：获取问题中每个变量的最大值。 5. `nPop = params.nPop;`：获取粒子群的个体数量。 6. `VarSize = [1 nVar];`：定义一个大小为 1 行 nVar 列的矩阵 VarSize，用来表示每个粒子的位置和速度的大小。 7. `empty_particle.Position = [];`：定义一个空粒子，包含位置、速度和成本。 8. `empty_particle.Velocity = [];` 9. `empty_particle.Cost = [];` 10. `empty_particle.Best.Position = [];`：定义空粒子的历史最优位置和最优成本。 11. `empty_particle.Best.Cost = [];` 12. `particle = repmat(empty_particle, nPop, 1);`：使用 repmat 函数将空粒子复制 nPop 次，得到一个大小为 nPop 行 1 列的粒子矩阵。 13. `GlobalBest.Cost = inf;`：初始化全局最优成本为无穷大。 14. `for i=1:nPop`：对于每个粒子进行循环。 15. `particle(i).Position = unifrnd(VarMin, VarMax, VarSize);`：将第 i 个粒子的位置初始化为在 VarMin 和 VarMax 之间均匀分布的随机值。 16. `particle(i).Velocity = zeros(VarSize);`：将第 i 个粒子的速度初始化为大小为 VarSize 的零向量。 17. `particle(i).Cost = CostFunction(particle(i).Position);`：计算第 i 个粒子的成本，即将其位置传入成本函数中计算得到。 18. `particle(i).Best.Position = particle(i).Position;`：将第 i 个粒子的历史最优位置初始化为其当前位置。 19. `particle(i).Best.Cost = particle(i).Cost;`：将第 i 个粒子的历史最优成本初始化为其当前成本。 20. `if particle(i).Best.Cost < GlobalBest.Cost`：如果第 i 个粒子的历史最优成本小于全局最优成本。 21. `GlobalBest = particle(i).Best;`：将全局最优位置更新为第 i 个粒子的历史最优位置。 22. `end`：结束 for 循环。 23. `end`：结束函数。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值