MATLAB程序初始化OpenFOAM颗粒位置

原创

已于 2023-08-10 19:02:44 修改 · 1k 阅读

3 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #开发语言

于 2023-08-10 19:01:28 首次发布

本文介绍了如何在OpenFOAM的颗粒两相流模拟中，使用MATLAB编程自定义颗粒的初始位置和速度，涉及创建kinematicCloudPositions文件并提供示例代码。

该文章已生成可运行项目，

问题引入

在OpenFOAM的颗粒两相流求解器中，我们可以采用manualInjection的方式进行自定义颗粒的初始位置，这个命令十分方便，在CFDEM中也有类似的命令，不过CFDEM中的命令更加强大，我们不仅可以定义颗粒的初始位置，而且还可以定义颗粒的初始速度,唯一的缺点就是CFDEM学起来资料比较少。。。
anyway，这里我们介绍如何在OpenFOAM中实现这个操作。
其核心步骤是在cosntant文件夹下创建kinematicCloudPositions文件，其大致内容如下：
在这里插入图片描述
于是自定义粒子的位置的思路就非常简单，我们利用MATLAB创建上面的文件内容即可。

MATLAB编程

这里不对程序做解说，直接放在这里。

%% write by rqli 2023/8/10
%% MATLAB程序初始化OpenFOAM颗粒位置（颗粒位置为均匀分布）
clc;clear;format long

%% 各个方向颗粒的数目
NumXofP = 5;NumYofP = 5;NumZofP = 1;
%% 计算域参数
Lx = 0.278;Ly = 0.278;%这里是二维计算域，因此只有高度和宽度
Lz = 0.278;
%% 生成初始颗粒位置，每个点处速度的标准形式为(0.0，0.0，0.0)
position = zeros(NumXofP,NumYofP);
x = linspace(Lx/(2*NumXofP),Lx-Lx/(2*NumXofP),