mfix中输出DEM颗粒的固相速度到网格
基于mfix-19.1.2版本
方法一:直接输出差值网格固相速度
注:这种方式只适用于garg 2012
颗粒差值格式
在DEM中,默认是无法输出固相速度的网格值的:
但是通过搜索des文件夹下V_s
关键字可以看到,实际上网格固相速度是会通过差值计算得到的:
并且在计算DEM颗粒脉动速度的地方也用到了网格固相速度:
说明网格固相速度只是默认没有写出,因此只需要人为写出一下即可。
首先找到网格输出的代码位置:
由于MMAX
代表TFM的固相数,因此在循环范围部分需要修改,默认只写出TFM的网格固相速度。修改如下:
最后在dat文件里手动添加一下VTK_U_S(VTK, PHASE)
这个关键字(且确保颗粒的差值格式为garg 2012
)例如:
结果如下:
方法二:保存网格内颗粒的平均速度到网格
网格颗粒速度还可以仿照DES_GRANULAR_TEMPERATURE
函数内计算网格拟颗粒温度的方式,对网格内颗粒速度进行统计平均。
按照J. Musser的建议,可以把需要后处理的变量存到ReactionRates
这个数组中,虽然这个数组本来是用来后处理反应速率的,不过由于可以直接在gui里设置该数组的个数以及输出频率,而且该变量仅仅用来保存输出,计算中未被使用,因此可以灵活使用,来保存其他需要保存的变量。
例如这里我要保存网格内颗粒Y方向的平均速度V_s
,可以先设置一个ReactionRates
数组,以备后面使用:
在des_granular_temperature
中USE这个数组:
定义一个局部变量用于后面保存网格内颗粒速度的和:
将颗粒速度累加到对应网格内的M
相内:
除以网格内M
相颗粒个数,得到网格内颗粒的平均速度,并保存到ReactionRates(ijk,1)
中。
注意拟颗粒温度由于其定义跟维数有关,因此要多除以一个维数,但是速度不需要:
添加代码如下:
现在,只要调用这个函数,就会往ReactionRates
中写入固相速度,因此只要在输出的部分,每次输出之前调用一次DES_GRANULAR_TEMPERATURE
就可以了:
重新编译计算可以看到输出Y方向速度:
mfix中输出DEM颗粒的固相速度到网格的更多相关文章
- shell单引号中输出参数值
因为在shell的单引号中,所有的特殊字符和变量都会变成文本,那么如果需要在字符串中输出变量值怎么办呢? 这里记录以下两种方法: 使用双引号 shell> X='parameter' shell ...
- python3中输出不换行
python2中输出默认是换行的,为了抑制换行,是这么做的: print x, 到了python3中,print变成一个函数,这种语法便行不通了.用2to3工具转换了下,变成这样了: print(x, ...
- 在 ASP.NET MVC Web 应用程序中输出 RSS Feeds
RSS全称Really Simple Syndication.一些更新频率较高的网站可以通过RSS让订阅者快速获取更新信息.RSS文档需遵守XML规范的,其中必需包含标题.链接.描述信息,还可以包含发 ...
- 在JSP页面中输出JSON格式数据
JSON-taglib是一套使在JSP页面中输出JSON格式数据的标签库. JSON-taglib主页: http://json-taglib.sourceforge.net/index.html J ...
- MVC中如何在controller的action中输出JS到页面上
MVC中如何在controller的action中输出JS到页面上 可以通过Http上下文对象(httpContext)就可以了,在Action中的HttpContext就是这个Action所指向的页 ...
- makefile中使用echo向文件中输出版本号和编译时间
@echo "#define BUILD_TIME" `date +"%F_%H:%M:%S"` > buildTime_svnVer.h @echo & ...
- 在jsp中用一数组存储了数据库表中某一字段的值,然后在页面中输出其中的值。
List<String> list = new ArrayList<String>(); String sql = "select userName from us ...
- 在Spring Boot中输出REST资源
前面我们我们已经看了Spring Boot中的很多知识点了,也见识到Spring Boot带给我们的各种便利了,今天我们来看看针对在Spring Boot中输出REST资源这一需求,Spring Bo ...
- 在Gazebo中使用DEM構建起伏地形環境
所需資料下載地址: 1. https://bitbucket.org/osrf/gazebo_tutorials/raw/default/dem/files/ 数字高程模型(致謝谷歌翻譯)概述数字高程 ...
随机推荐
- [SoapUI]怎样获取上一个Test Step的名字
def currentStepInd = context.currentStepIndex def previousStep = testRunner.testCase.getTestStepAt(c ...
- 何时使用[self release]
这样的语句 [self release]; 乍看上去让人很困惑. 从release方法本身的作用上来说,就是给self的引用技术减一,就像release对其他对象所做的一样. 一般来说,唯一用到, ...
- Java之RandomAccessFile小结
今天跟大家分享一下javase中的关于I/O的操作: 有时我们需要在文件的末尾追加一些内容,在这时用RandomAccessFile就很好. 这个类有两个构造方法: RandomAccessFile( ...
- [转]快速入门系列--WebAPI--01基础
本文转自:http://www.cnblogs.com/wanliwang01/p/aspnet_webapi_base01.html ASP.NET MVC和WebAPI已经是.NET Web部分的 ...
- leetcode-8-String to Integer (atoi) (已总结)
8. String to Integer (atoi) Total Accepted: 93917 Total Submissions: 699588 Difficulty: Easy Impleme ...
- 转发一篇能看懂的关于ingress的说明
最近发现好多人问 Ingress,同时一直也没去用 Nginx 的 Ingress,索性鼓捣了一把,发现跟原来确实有了点变化,在这里写篇文章记录一下 https://mritd.me/2017/03/ ...
- Zero Clipboard js+swf实现的复制功能使用方法
开发中经常会用到复制的功能,在 IE 下实现比较简单.但要想做到跨浏览器比较困难了.本文将介绍一个跨浏览器的库类 Zero Clipboard .它利用 Flash 进行复制,所以只要浏览器装有 Fl ...
- django drf 改变retrive的pk查询字段
lookup_filed可以改变retrive查询时默认以pk查询的逻辑 from django.shortcuts import render from rest_framework import ...
- JS 面向对象详解
面向对象详解1 OO1.html <!DOCTYPE html> <html> <head> <meta charset="utf-8" ...
- TortoiseGit revert failed - unable to revert local changes
TortoiseGit revert failed - unable to revert local changes 下载完按照提示一步一步安装就可以了 成功之后就是windows.1版本,如图