The Development of an Advanced Solid-Rocket Flow Simulation Program ROCFLO
作者:Alavilli P, Tafti D, Najjar F
简介
流体与固体有两方面的耦合:(1)燃烧导致的燃面推移(2)压力导致的装药大变形
The high gas pressures generated in the core fow load the propellant and the outer structure, thus deforming it. This deformation modifies the fow field which in turn alters the pressure field,and feeds back to the structural dynamics. The flow dynamics and the structural response are thus intimately coupled. Propellant combustion is also itself affected by the overlaid fluid dynamic field through the imposed pressure and temperature gradients. These aspects serve to illustrate the need for coupled multi-component simulation for the analysis and design of solid rocket motors.
使用动网格方法
where, U U U is the set of conserved variables. F F F and V V V are the inviscid and viscous flux vectors respectively. g g g is the grid speed vector and the term − g ⃗ U n ⃗ d A -\vec{g}U\vec{n}\mathrm{d}A −gUndA is the correction to the convective fluxes induced purely by mesh motion. S S S represents additional source terms that sometimes arise, for example in species equations with chemical reactions or in rotating systems, but is not considered here. Defining the modified convective flux vector F ⃗ = F ⃗ − g ⃗ U \vec{\mathcal{F}}=\vec{F}-\vec{g}U F=F−gU
网格运动可以分为三部分:药柱变形、燃面推移、网格自适应
有限体积法的半离散形式
使用带人工粘性的中心差分格式离散对流项,使用多阶段龙格库塔法离散时间1
非并行的计算流程如下:
并行的计算流程如下
针对航天飞机的固体火箭助推器开展仿真2,使用Gridgen生成网格
首先是稳态案例:网格不动,使用初始的几何模型开展仿真,得出结论
Note the huge difference caused to the chamber pressure by the increased burning surface provided by the star grain section of the motor. Accurate modeling of the effective burning surface is thus important. This portends the need for an accurate propellant surface regression model in order to achieve accurate long time burn simulations of solid rockcts.
然后是瞬态案例:只仿真了0.1s,出现了激波,这是由于假设初始时刻全部燃面都被点燃导致的。(an assumption of instantaneous
propellant ignition.),对于点火问题,见作者的其他论文3。
遗憾的是,这篇文章并没有给出复杂装药的动网格燃面推移(在作者其他文章中体现4),而且事实上并未结合装药的弹性变形,并未考虑侵蚀燃烧。
最后,给出了一些典型发动机5的计算结果。
A.Jameson, W. Schmidt,and E. Turkel. Numerical Solutions of the Euler Equations by Finite Volume Methods using Runge-Kutta Time Stepping Schemes. AIAA paper 81-1259, 1981. ↩︎
Thiokol Corporation. Design Data Book for Space Shuttle Reusable Solid Rocket Motor. TWR-16881 Revision A, Publi-cation no.930480, 1997. ↩︎
Alavilli,P.,Buckmaster,J.,Jackson,T.and Short,M. Ignition-Transient Model for Solid Propellant Rocket Motors. AIAA Paper.36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit,16-19 July 2000,Hunstville, AL. ↩︎
Alavilli,P.,Najjar, F.,Krier,H.,Tafti,D.and Fiedler, R. Dynamic Adaptive Simulation of Flows in Solid Propellant Rockets.AIAA Paper.36th AIAA/ASME/SAE/ASEE Joint Propulsion Conf.and Exhibit,16-19 July 2000, Hunstville,AL. ↩︎
Blomshield,F.S.,Crump, J.E.,Mathes,H.B.,Stalnaker, R. A.and Beckstead, M. W. Stability Testing of Full-Scale Tactical Motors. Journal of Propulsion 8 Power,13(3):349-355.1997. ↩︎