一种基于MATLAB和paraview的fluent-DPM颗粒后处理方法技术

本发明专利技术公开了一种基于MATLAB和paraview的fluent

【技术实现步骤摘要】
一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法


[0001]本专利技术涉及计算流体力学领域仿真计算
，具体来说，涉及一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法。

技术介绍

[0002]随着流体力学理论的不断完善和高性能计算机技术的飞速发展，二者相互结合形成的计算流体动力学开始成为解决各种工程问题的重要手段。在实际的工程应用中，经常会遇到流体介质中含有固体颗粒的情况，商业软件Fluent成为解决这类含颗粒流动问题的主流方案之一。Fluent可以通过求解各种流体力学方程来对水利机械中的介质流动进行计算，将阀门、泵等内部的介质流动具体的呈现出来。在Fluent给出计算结果之后，还需要对结果进行处理，以得到水利机械内部介质的速度、压力、流线等，实现介质流动的可视化。这种可视化处理便被称为“后处理”。但是，Fluent自身对于计算结果的后处理并不是非常强大，经常出现不能满足用户需求的情况，特别是介质中含有颗粒时，Fluent自身并不能对计算的DPM颗粒进行深层次以及复杂的处理，就连一部分针对Fluent进行优化的后处理软件，对DPM颗粒的处理也不是非常友好。于是，Fluent计算出来的DPM颗粒的后处理成为仿真领域的难题之一。
[0003]ParaView是能够对二维和三维数据进行分析和可视化的非常强大的开源程序，其开源性和数据可视化为展示Fluent计算的DPM颗粒提供了可能性。但是ParaView并不能直接对FluentDPM颗粒进行处理，国内外文献资料也没有与之相关的案例。
[0004]针对相关技术中的问题，目前尚未提出有效的解决方案。

技术实现思路

[0005]针对相关技术中的问题，本专利技术提出一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法，以克服现有相关技术所存在的上述技术问题。
[0006]为此，本专利技术采用的具体技术方案如下：
[0007]一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法，包括以下步骤；
[0008]S1：通过Fluent软件对DPM颗粒进行计算；
[0009]S2：在Fluent软件中将DPM颗粒的基本数据(颗粒坐标、颗粒速度、颗粒直径、颗粒ID)导出后缀名为“.dpmrpt”格式的文件；
[0010]S3：在MATLAB中，对存储DPM颗粒数据的“.dpmrpt”文件进行处理，并将结果写入“.vtk”文件；
[0011]S4：在ParaView中打开“.vtk”文件，对其进行后处理。
[0012]作为优选的，所述步骤S3中；
[0013]S31：利用MATLAB创建“.m”文件，将颗粒坐标、颗粒速度、颗粒直径、颗粒ID储存到矩阵中；
[0014]S32：利用MATLAB重新定义储存着颗粒坐标、颗粒速度、颗粒直径、颗粒ID信息的矩
阵；
[0015]S33：利用MATLAB中的循环结构对矩阵进行重新排列；
[0016]S34：利用MATLAB将必要信息和标记信息写入到“.vtk”文件；
[0017]S35：利用MATLAB创建“.vtk”文件，将重新排列的矩阵写入“.vtk”文件。
[0018]作为优选的，所述步骤S31的关键在于：
[0019]通过fid＝fopen('filename'，'permission')命令打开“.dpmrpt”文件，并通过lines＝textscan(fid，'％f％f％f％f％f％f％f％f％f％[^\n]'，'Headerlines'，17)命令读取“.dpmrpt”文件中的颗粒基本数据(颗粒坐标、颗粒速度、颗粒直径、颗粒ID)，将其转化为“M行
×
1列”的矩阵数据。
[0020]作为优选的，所述步骤S32的关键在于：
[0021]利用MATLAB重新定义储存着颗粒基本数据(颗粒坐标、颗粒速度、颗粒直径、颗粒ID)信息的矩阵，重新定义之后的矩阵仅仅作为处理过程的中间量，而不作为写入“.vtk”文件的最终结果，并在重新定义的同时将颗粒直径转化为颗粒半径。
[0022]作为优选的，所述步骤S33的关键在于：
[0023]通过for循环将储存颗粒基本数据信息的矩阵重新进行排列，将“M行
×
1列”的矩阵数据转化为“N行
×
9列”的矩阵；
[0024]作为优选的，所述步骤S34的关键在于：
[0025]利用MATLAB将能够标记颗粒基本数据(颗粒坐标、颗粒速度、颗粒直径、颗粒ID)的标记信息，通过fprintf(fileID，'format'，'vtkData')命令写入到“.vtk”文件，以在ParaView后处理步骤中能更方便快捷的识别、标记和区分各种不同的颗粒信息，同时也将构成“.vtk”文件的必需信息写入“.vtk”文件中。
[0026]作为优选的，所述步骤S35的关键在于：
[0027]通过fprintf(fileID，'format'，'vtkData')命令将转化后的矩阵写入到“.vtk”文件中，并且放置于所述步骤S34中的必要信息和标记信息之后，保证颗粒标记信息和与之相关矩阵能够相互对应。
[0028]作为优选的，所述步骤S3中使用的MATLAB程序代码为(双百分号％％后面的内容为对代码的解读和注释)：
[0029][0030][0031][0032][0033][0034]本专利技术的有益效果为：在MATLAB中使用前文所述的程序，将Fluent计算出来的DPM
颗粒文件转化成为ParaView可识别的文件，使得Paraview能够对Fluent计算出来的DPM颗粒进行后处理。与现有的FluentDPM颗粒后处理方法相比较，本专利提供的程序和方法能够将商业软件Fluent计算出来的DPM颗粒放置在开源环境中，使得DPM颗粒的后处理、可视化等更加方便快捷，在原有的基础上获得更加多变的数据以及进行更加复杂的操作。
附图说明
[0035]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0036]图1是根据本专利技术实施例的一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法的流程示意图；
[0037]图2是根据本专利技术实施例的操作示意图；
[0038]图3是根据本专利技术实施例的提及的颗粒信息图；
[0039]图4是根据本专利技术实施例的通过ParaView对DPM颗粒进行部分后处理的示意图。
具体实施方式
[0040]为进一步说明本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法，其特征在于，包括以下步骤；S1：通过Fluent软件对DPM颗粒进行计算；S2：在Fluent软件中将DPM颗粒的基本数据(颗粒坐标、颗粒速度、颗粒直径、颗粒ID)导出后缀名为“.dpmrpt”格式的文件；S3：在MATLAB中，对存储DPM颗粒数据的“.dpmrpt”文件进行处理，并将结果写入“.vtk”文件；S4：在ParaView中打开“.vtk”文件，对其进行后处理。2.根据权利要求1所述的一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法，其特征在于，所述步骤S3中；S31：利用MATLAB创建“.m”文件，将颗粒坐标、颗粒速度、颗粒直径、颗粒ID储存到矩阵中；S32：利用MATLAB重新定义储存着颗粒坐标、颗粒速度、颗粒直径、颗粒ID信息的矩阵；S33：利用MATLAB中的循环结构对矩阵进行重新排列；S34：利用MATLAB将必要信息和标记信息写入到“.vtk”文件；S35：利用MATLAB创建“.vtk”文件，将重新排列的矩阵写入“.vtk”文件。3.根据权利要求2所述的一种基于MATLAB和paraview的fluent
‑
DPM颗粒后处理方法，其特征在于，所述步骤S31的关键在于：通过fid＝fopen('filename'，'permission')命令打开“.dpmrpt”文件，并通过lines＝textscan(fid，'％f％f％f％f％f％f％f％f％f％[^\n]'，'Headerlines'，17)命令读取“.dpmrpt”文件中的颗粒基本数据(颗粒坐标、颗粒速度、颗粒直径、颗粒ID)，将其转化为“M行
×
1列”的矩阵数据。4.根据权利要求2所...

【专利技术属性】
技术研发人员：林哲，孙希望，朱祖超，郑旭，陈德胜，
申请(专利权)人：浙江理工大学，
类型：发明
国别省市：