首页/文章/ 详情

Fluent VOF水下固体火箭发射仿真(一)

2天前浏览32

本案例利用VOF模型对水下固体火箭(10m水深)发动机点火初期的流场进行了仿真。该案例所用模型为假设模型,且缺少相关燃气参数,仅作计算设置参考。通过此案例后续跨可以对不同水深、不同模型的工况展开类似仿真计算。

1 假设说明

为简化计算流程做出如下假设:

(1)对于水下高温燃气,忽略发动机推进剂的烧蚀产物,将其视为具有温、压的理想气体,考虑粘连特性及可压缩性;

(2)对于环境水介质,由于其密度相对于气体十分稳定,将其考虑为有粘不可压缩流体;

(3)对于气液边界,不考虑两相间的掺混与相变,但考虑两相间的热量传递;

(4)对于浮力的影响,由于发动机喷管出口竖直向上在静水中工作,同时两相密度差大,不忽略浮力

2 workbench 设置

本案例以二维轴对称展开计算,因计算模型简单,计算网格采用结构化网格进行划分。需要选择几何+icem+fluent三个软件,相关的workbench设置如下图:

3 SCDM 设置

3.1 导入几何

计算域包括发动机内流域(喷管)和外界计算域两个部分。模型中Laval喷管的喉部直径为20mm,出口直径为38mm。计算域大小为1250 mm × 350 mm。相关几何模型如下图。

4 ICEM 设置

ICEM为结构化网格划分软件,由于模型简单,不再单独介绍网格划分流程,相关的网格划分结果如下图所示,因为模型具有对称性,数值仿真模型采用二维轴对称模型,网格划分仅需考虑轴对称中一侧的网格即可。

5 FLUENT 设置

5.1 General设置与网格导入

由于采用二维轴对称模型,此处需要勾选轴对称选项。  

5.2 表达式设置

本案例需要对喷管入口的相关参数进行设置,主要进行压强和温度的设置,根据大量文献和实验的结果,温度和压强在达到稳定前可以拟合为随时间线性变化的方程,此处可以使用UDF和表达式两种方式进行设置,此处仅考虑了压强和温度为达到稳定前的计算。所以喷管入口的压强和温度的表达式设置比较简单。相关设置如下图所示。  

5.3 边界条件设置

首先将下图所示的对称轴设置为对称轴边界,其余边界设置如下图所示。

需要注意入口边界设置,需要调用表达式的相关参数。

出口边界为压力出口,为水深环境压力,本案例采用水深10m的环境,相关环境压力为0.2MPa

5.4 材料设置

由于是对水下发射问题展开仿真,因此需要采用水和燃气两种材料,此处用空气(理想气体)直接代替燃气,具体的燃气参数应该根据燃烧室相应的物性参数估算获得。具体的操作如下图。

5.5 模型设置

选择理想气体后,能量方程默认打开,此处需要进一步打开VOF模型,因为不考虑相变等其他问题,此处操作简单,不进行Lee模型设置。VOF模型相关设置如下图所示。


5.6 初始化设置

VOF初始化设置一般采用标准初始化,需要对不同的区域展开初始化。分配不同的材料。

此处对发动机内流域采用3.2MPa进行初始化,外部水域采用环境压力进行初始化。

相关初始化设置如下图。

5.7 计算设置

此处计算时间步长设置为1e-610000步。

此处建议选取更小的时间步长进行相关计算。

6 后处理结果

6.1 后处理云图结果

受计算时间与计算资源影响,此处未进行完整计算,仅计算了00000s

除此之外,还可以展开对发动机推力、尾流压力云图、燃汽泡演化等后处理结果进行分析,本文不再进一步阐述。

来源:CFD仿真库
FluentWorkbench燃烧UDF材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-03-21
最近编辑:2天前
CFD仿真库
硕士 分享知识,互帮互助,共同进步。
获赞 1粉丝 2文章 55课程 0
点赞
收藏
作者推荐

Fluent 重叠网格+UDF 扑翼机气动仿真(一)

本案例利用Fluent重叠网格与UDF,对扑翼机的气动特性展开仿真。该案例所用模型为假设模型,仅作计算设置参考。通过此案例后续可以对进一步添加udf代码与更换模型,实现更为复杂的扑翼机运动,对其展开气动仿真计算。1UDF说明在本研究中采用重叠网格模型对扑翼机扑翼运动进行模拟。本案例选择DEFINE_CG_MOTION进行定义,omega[0]代表z轴旋转方向,本案例设计翼型上下摆动18°,相关的UDF代码如下:#include"udf.h"#include"mem.h"#include"dynamesh_tools.h"DEFINE_CG_MOTION(pyj,dt,vel,omega,time,dtime){NV_S(vel,=,0.0);NV_S(omega,=,0.0);omega[0]=0.314*cos(2*3.14*time);}2workbench设置本案例需要设置如下三个模块的计算,其中包括背景网格区域、前景网格区域与fluent计算三个部分,具体设置如下图:3SCDM设置3.1导入几何整体几何结构如下图:扑翼机翼型采用NACA0012,具体的几何结构如下图,x轴正向为压力出口,负轴位速度入口,扑翼机表面为壁面,其余面位对称面。重叠网格区域为内部圆柱区域。扑翼机运动翼型命名为naca,静止区域命名为bird。4FluentMeshing设置4.1网格设置采用SCDM进行网格划分,背景网格与前景网格皆采用六面体网格划分,并划分相对应的边界层网格。背景网格如下图所示:前景网格如下图所示:5FLUENT设置5.1General设置与网格导入首先导入背景网格,其次通过下图所示的方法将前景导入。由于本文涉及到NACA0012翼型运动,因此需要探讨瞬态计算结果,此处的设置比较简单,勾选为瞬态计算。5.2边界条件设置由于本案例采用的计算模型为半模,对称面命名为axi用于区分。首先进行视图设置:其余各边界条件的设置如下:速度入口设置为5m/s,具体设置如下图:5.3动网格设置首先将本案例采用的UDF导入,具体的导入方式如下图。其次对运动区域进行设置,即整个运动区域到网格都要添加udf,具体设置如下:5.4初始化设置相关初始化设置如下图。5.5计算设置此处进行的计算设置如下:来源:CFD仿真库

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈