ection的方式进行自定义颗粒的初始位置这个命令十分方便在CFDEM中也有类似的命令不过CFDEM中的命令更加强大我们不仅可以定义颗粒的初始位置而且还可以定义颗粒的初始速度,唯一的缺点就是CFDEM学起来资料比较少。

。

。

anyway这里我们介绍如何在OpenFOAM中实现这个操作。

其核心步骤是在cosntant文件夹下创建kinematicCloudPositions文件其大致内容如下

于是自定义粒子的位置的思路就非常简单我们利用MATLAB创建上面的文件内容即可。

MATLAB编程

MATLAB程序初始化OpenFOAM颗粒位置颗粒位置为均匀分布

clc;clear;format

生成初始颗粒位置每个点处速度的标准形式为(0.00.00.0)

position

linspace(Lx/(2*NumXofP),Lx-Lx/(2*NumXofP),NumXofP);

linspace(Ly/(2*NumYofP),Ly-Ly/(2*NumYofP),NumYofP);

[X,Y]

reshape(X,NumXofP*NumYofP*NumZofP,1);postion_y

reshape(Y,NumXofP*NumYofP*NumZofP,1);postion_z

zeros(NumXofP*NumYofP*NumZofP,1);

elseif

linspace(Lz/(2*NumZofP),Lz-Lz/(2*NumZofP),NumZofP);

X3d

[Z3d;z(i).*ones(NumXofP,NumYofP)];endpostion_x

reshape(X3d,NumXofP*NumYofP*NumZofP,1);postion_y

reshape(Y3d,NumXofP*NumYofP*NumZofP,1);postion_z

reshape(Z3d,NumXofP*NumYofP*NumZofP,1);disp(done)

end%%

zeros(NumXofP*NumYofP*NumZofP,1);

position_xyz

1:NumXofP*NumYofP*NumZofPposition_xyz(i)

strcat((,num2str(postion_x(i)),

end

writematrix(position_xyz,position_file);二维粒子初始位置如下(采用scatter函数)

我们把上面的坐标粘贴到kinematicCloudPositions文件里即可。