1 问题描述

本案例演示挂载从机翼上投放后的运动行为以及流场分布。

2 Fluent设置

以3D、Double Precision模式启动Fluent
读取网格文件delta.msh.gz

2.1 General设置

选择使用Density-Based
激活Gravity，设置重力加速度为Z方向9.807 m/s2

注：这里先计算稳态，目的是确保计算稳定性。

2.2 Models设置

选择使用Inviscid无粘模型

2.3 Materials设置

指定空气的Density为ideal-gas

2.4 边界条件设置

指定farfield的Mach Number为1.2
指定流动方向为X轴-1

指定温度为216.65 k

指定出口边界outflow回流温度为216.65 k

2.5 操作条件设置

指定参考压力为20646 Pa

注：20646 Pa指的是11600 m高度处的压力。

2.6 初始化计算

初始化计算

2.7 迭代计算

指定迭代次数1000
激活选项Solution Steering，指定Flow Type为supersonic，如下图所示

可以如下图所示修改求解控制参数

计算求解。

2.8 稳态计算结果

马赫数分布

2.9 编译并加载UDF

编译UDF并加载

案例UDF比较简单，主要定义了挂载的质量及外力信息。

#include "udf.h"
 DEFINE_SDOF_PROPERTIES(delta_missile, prop, dt, time, dtime)
 {
    prop[SDOF_MASS] = 907.185;
    prop[SDOF_IXX] = 27.116;
    prop[SDOF_IYY] = 488.094;
    prop[SDOF_IZZ] = 488.094;
 
    /* add injector forces, moments */
    {
      register real dfront = fabs (DT_CG (dt)[2] - (0.179832*DT_THETA (dt)[1]));
      register real dback = fabs (DT_CG (dt)[2] + (0.329184*DT_THETA (dt)[1]));
 
      if (dfront <= 0.100584)
        {
          prop[SDOF_LOAD_F_Z] = 10676.0;
          prop[SDOF_LOAD_M_Y] = -1920.0;
        }
 
      if (dback <= 0.100584)
        {
          prop[SDOF_LOAD_F_Z] += 42703.0;
          prop[SDOF_LOAD_M_Y] += 14057.0;
        }
    }
    printf ("ndelta_missile: updated 6DOF properties");
 }