火炸药学报    2018, Vol. 41 Issue (5): 506-511   DOI: 10.14077/j.issn.1007-7812.2018.05.015
0

引用本文  

乌岳, 李卓, 鲁荣. 固体推进剂花板浇注的数值模拟与实验研究[J]. 火炸药学报, 2018, 41(5): 506-511. DOI: 10.14077/j.issn.1007-7812.2018.05.015
WU Yue, LI Zhuo, LU Rong. Numerical Simulation and Experimental Study of Flower Plate Pouring System for Solid Propellant[J]. Chinese Journal of Explosives & Propellants, 2018, 41(5): 506-511. DOI: 10.14077/j.issn.1007-7812.2018.05.015

基金项目

中国航天科工集团第六研究院院管课题(GK201812);内蒙古工业大学科学研究项目(No.X201718)

作者简介

乌岳(1986-), 女, 博士研究生, 讲师, 从事固体推进剂研究。E-mail:569704028@qq.com

文章历史

收稿日期:2018-03-12
修回日期:2018-05-29
固体推进剂花板浇注的数值模拟与实验研究
乌岳1,3, 李卓1, 鲁荣2     
1. 内蒙古工业大学理学院, 内蒙古 呼和浩特 010051;
2. 中国航天科工集团 内蒙古航天红峡化工有限公司, 内蒙古 呼和浩特 010051;
3. 内蒙古工业大学化工学院, 内蒙古 呼和浩特 010051
摘要: 利用Herschel-Bulkley本构模型,通过CFD黏弹性流体软件Ansys-Polyflow对固体推进剂药浆浇注工艺过程进行数值模拟,研究了真空浇注系统不同孔径条件下药浆的流动速度、挤出胀大效应、黏度、剪切速率及流平特性,并通过实验对数值模拟结果进行了验证。结果表明,花板孔径越大,药浆流速越快,7 mm孔径条件下药浆流速最大可达0.10 m/s,而3 mm孔径流速最大值只有0.02 m/s;挤出胀大比随浇注孔长径比(L/D)的减小而增大,孔径为3、5、7 mm时的挤出胀大比在同一花板厚度下依次为1.05、1.31、1.44。HTPB四组元推进剂药浆为典型的屈服假塑性流体,当剪切应力低于某极限值时体系不流动,而高于该极限值时体系呈假塑性流体,黏度随剪切速率的增大而降低;当药浆流经花板和翼片时,流体剪切速率达到最大值,药浆黏度从最初的400 Pa·s降至300 Pa·s;药浆整体的流平性较好,未出现较大的孔洞和瑕疵;燃烧室芯孔右侧的微小空洞随着流动逐渐消失,直到燃烧室充满,云图显示孔洞的体积分数在0.5以下时,基本充满。
关键词: 固体推进剂药浆     本构模型     花板真空浇注     Ansys-Polyflow软件    
Numerical Simulation and Experimental Study of Flower Plate Pouring System for Solid Propellant
WU Yue1,3, LI Zhuo1, LU Rong2     
1. College of Science, Inner Mongolia University of Technology, Hohhot 010051, China;
2. Inner Mongolia Aerospace Hongxia Chemical Co., Ltd., China Aerospace Science & Industry Corporation, Hohhot 010051, China;
3. College of Chemical Engineering, Inner Mongolia University of Technology, Hohhot 010051, China
Abstract: The Herschel-Bulkley constitutive model was used to numerically simulate the process of solid propellant slurry pouring process through the CFD viscoelastic fluid software Ansys-Polyflow. The flow rate, extrusion swelling effect, viscosity, shear rate and leveling properties of slurry under the condition of different apertures in the vacuum pouring system were studied and the numerical simulation results were verified through experiments. Results show that:the bigger the aperture of the flower plate, the faster the flow rate of the slurry. Under the condition of 7 mm aperture, the maximum flow rate of the slurry can reach 0.10 m/s, while the maximum flow rate of 3 mm aperture is only 0.02 m/s. The extrusion swelling ratio increases with decreasing the pouring hole's length-diameter ratio (L/D). The extrusion swelling ratio at the aperture of 3, 5 and 7 mm under the same plate thickness is 1.05, 1.31 and 1.44, respectively. HTPB four-component propellant slurry is a typical pseudoplastic fluid of yield. When the shear stress is lower than a limit value, the system does not flow, while when the shear stress is higher than the limit value, the system is a pseudoplastic fluid, and the viscosity decreases with the increase of shear stress. When the slurry flows through the flower plate and the blade, the shear rate of fluid reaches the maximum, and the viscosity of the slurry drops from 400 Pa·s initially to 300 Pa·s. The overall fluidity of the slurry is good, without large holes and defects. The tiny cavity on the right side of core hole of the combustion chamber gradually disappears with the flow until the combustion chamber is full. When the volume fraction of the hole is below 0.5, the cloud map shows that the hole is basically full.
Keywords: solid propellant slurry     constitutive model     flower plate vacuum pouring     Ansys-Polyflow software    
引言

固体推进剂药浆兼具黏性体和弹性体的双重特性,是典型的黏弹性非牛顿流体。固体推进剂低黏度药浆的浇注通常采用真空花板浇注工艺[1],推进剂药浆流动的驱动力为料斗和浇注缸之间的压差以及料斗内推进剂药浆自重[2]。其本构方程和流动边界条件的复杂性,很难采用理论解析方法分析其流动[3]。目前随着计算机辅助技术的迅速发展,数值模拟已成为分析黏弹流体流动的主要手段。

姜爱民等[4]对固体推进剂药浆在混合釜内的流动过程进行了数值模拟,结合网格叠加技术和网格加密技术,对混合过程中药浆流动进行了有限元模拟,详细研究了捏合机混合过程中混合釜内剪切应力、各点的速度、混合指数、黏度的分布和极值关系;Connelly C K等[5]利用粒子示踪法,模拟了行星式混合内玉米糖浆的三维混合流动;Shimada T等[6]利用Polyflow软件将混有铅示踪颗粒的推进剂药浆浇注到一个透明型腔中,并在周围架设了两个方向的X射线透析仪器,通过计算得到药浆内颗粒的三维迹线,从而模拟得出比较直观的药浆内部的流场。

本研究利用黏弹性流体软件Polyflow提供的广义牛顿流体的本构模型,对复合推进剂药浆真空花板浇注流动的全过程进行数值模拟,针对真空花板浇注系统不同孔径的药浆流动速度、挤出胀大效应、药浆黏度和剪切速率及流平性做了一系列研究,同时通过实验对模拟结果进行了验证。

1 数值模拟

计算过程采用标准的CFD模拟分析流程:模型导入-划分网格-设置模型及参数-设置边界条件-利用CFD-POST对结果进行后处理。

1.1 推进剂浇注流动过程的基本假设

由于Polyflow软件本身的局限性,本研究依据实际情况作了以下假设:

(1) 流动为低雷诺数的不可压缩流体的层流流动[7],且流动过程不考虑推进剂药浆与外界环境的热交换(实际的浇注工房有循环水作保温);

(2) 推进剂药浆与壁面之间作无滑移边界条件设置;

(3) 由于Polyflow软件本身的设定,体积分数云图中若体积分数在0.5以上,则认为充满。

1.2 计算模型参数及工作原理

真空花板浇注工艺装置如图 1所示。

图 1 真空花板浇注工艺装置图 Figure 1 Process equipment of vacuum flower plate pouring

工作原理:将混合好的可浇注药浆定量放入料斗中保温(温度50℃),开启胶管阀,药浆按指定下料速度通过胶管阀进入真空度符合要求的燃烧室,将已除气的药浆通过真空压差和自重流入发动机燃烧室内壁和芯模之间的空隙内[8]

计算模型在图 1的基础上作了必要的简化设置,如图 2所示。

图 2 计算模型 Figure 2 Calculation model

选取各个装置的参数如下:料斗内径90 mm,花板直径90 mm,孔简化为36个,孔的直径依次为3、5、7 mm(实际中同一型号的花板孔径一致,为分析不同孔径的流速及挤出胀大效应,假设同一花板有不同孔径的分布),浇注缸内径200 mm,燃烧室内径100 mm,燃烧室采用后翼、环型槽结构[9],12个翼片。药柱芯孔最大处直径约40 mm,前机口直径约15 mm,后机口直径40 mm,燃烧室总长200 mm,直径100 mm。

1.3 网格划分

本研究中计算域形状不规则,为了提高花板网格质量将整个模型分割为38个体,分别对其进行结构化网格划分。总网格数为442 798个。网格划分情况如图 3所示。

图 3 网格划分 Figure 3 Mesh division
1.4 本构方程

研究表明[9-10],固含量较高的推进剂药浆往往表现出非线性宾汉流体特征:当剪切应力低于屈服值时体系不流动,而高于屈服值时,剪切应力与剪切速率呈非线性幂律关系,且黏度随剪切速率增大而降低,即屈服假塑性流体[11],亦称为Herschel-Bulkley塑性流体。

本研究针对广义牛顿流体等温流动模型Herschel-Bulkley进行计算,形式如下[11]

$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\eta = \frac{{{\tau _0}}}{{\dot \gamma }} + K{\left( {\frac{{\dot \gamma }}{{{{\dot \gamma }_0}}}} \right)^{n - 1}}, \dot \gamma >{{\dot \gamma }_0}\\ \eta = {\tau _0}\frac{{\left( {2 - \dot \gamma /{{\dot \gamma }_0}} \right)}}{{{{\dot \gamma }_0}}} + K\left[ {\left( {2 - n} \right) + \left( {n - 1} \right)\frac{{\dot \gamma }}{{{{\dot \gamma }_0}}}} \right], \dot \gamma <{{\dot \gamma }_0} \end{array} $

式中:τ0为屈服应力;${\dot \gamma _0}$为临界剪切速率;K为流体黏度系数;n为非牛顿指数。当n=1时,流体表现为线性牛顿流动,称为宾汉塑性流体;n<1时,流体表现为非线性剪切稀变流动,称为屈服假塑性流体;n>1时为非线性剪切增稠流动,称为屈服胀流体,如图 4所示。

图 4 广义宾汉流体的3类塑性流动 Figure 4 Three types of plastic flow of generalized Bingham fluids

文献[11]表明,该固体质量分数为88%的四组元HTPB推进剂为典型的屈服假塑性流体,当施加的应力超过屈服应力后,体系发生流动,此时流动符合n<1的幂律定律。

将文献[11]中$\eta - \dot \gamma $曲线的点利用Polymat拟合及外推后,得到拟合曲线。选定非牛顿指数n和屈服应力分别为0.64和47.1 Pa[11],拟合出黏稠系数和临界剪切速率分别为3 788和0.005 s-1

1.5 边界条件

本研究在数值计算中涉及的边界条件有:流动入口边界条件、流动出口边界条件和壁面边界条件。

工程文件任务一的流动入口采用作用法向应力fn和切向应力fs, 由于入口边界大气压强为0.1 MPa,即fn=1.0×105fs=0;流动出口边界采用作用法向速率vn和作用切向速率vs,且均为零;其余为壁面边界条件;由于选用有限元瞬态体积分数VOF(volume of fluid)模型,系统自动生成任务二Fluid Fraction Transport task,该任务边界条件入口采用Inflow边界条件,其余为壁面设置。

2 数值模拟结果与分析 2.1 花板孔径对药浆流动速度的影响

当黏弹性流体流经不同直径的花板孔时,其流动速率会有很大不同。圆截面不同孔径速度云图如图 5所示。

图 5 圆截面不同孔径速度云图 Figure 5 The cloud map with different orifice speeds in the circle section

图 5可以看出,孔径越大,流体的流速越大,且最大值出现在圆孔的中心轴线上:7 mm孔径条件下,流速最大值约0.10 m/s;5 mm孔径流速的最大值约0.04 m/s;3 mm孔径流速最大值约为0.02 m/s,这是由于孔径越大,流动阻力越小,质量流量越大导致流速越大,其次考虑到3种不同孔径,压强的相互竞争作用也会导致出现速度差异。

一般认为,与圆孔管壁处接触的流体层可视为静止的,因此圆孔管壁处流体流速接近于0[12]

2.2 花板孔径挤出胀大效应分析

挤出胀大效应(BARUS效应)是非牛顿流体的典型流变特征[13]图 6为流体出口挤出截面,其大小用挤出胀大比(B=Dmax/D)表示。

图 6 圆孔挤出截面 Figure 6 Export extrusion cross section

图 6可以看出,由于花板厚度为5 mm,所以3个不同孔径(3、5、7 mm)的长径比依次为5:3、5:5、5:7。通过Matlab软件编写程序拟合圆孔挤出截面直径,小圆(3 mm)、中圆(5 mm)、大圆(7 mm)挤出截面直径拟合结果分别为3.137 1、6.554 6和10.08 mm。即挤出胀大比依次为1.05、1.31和1.44。

由此可见,花板孔径的长径比强烈影响挤出胀大比。挤出胀大比随孔径长径比减小而增大,当孔径长度一致时,孔径越大挤出胀大比越大。

2.3 黏度和剪切速率云图

HTPB丁羟四组元推进剂为典型的非线性屈服假塑性流体,其主要特征是在某一剪切速率范围内,黏度随剪切速率的增大而降低,即剪切稀变流体,从黏度和剪切速率云图中可以得到证实。

图 7为某一YZ截面的黏度云图和剪切速率云图,图 8为某一XY圆截面翼片位置药浆的黏度云图和剪切速率云图。所采用药浆初始黏度为400 Pa·s。

图 7 某一YZ截面药浆的黏度云图和剪切速率云图 Figure 7 The viscosity and shear rate clouds for a YZ section slurry
图 8 某一XY圆截面翼片位置药浆的黏度云图、剪切速率云图和圆截面不同孔径速度云图 Figure 8 The viscosity and shear rate cloud maps of the fin position for a certain XY circular section

图 7图 8可以看出,药浆以较高的剪切速率流经花板和燃烧室翼片时,黏度降至最小值300 Pa·s,当离开喷孔和翼片时,所受的剪切速率迅速衰减。产生以上结果的原因为:

(1) 流体流经花板不同孔径时,剪切应力和剪切速率大幅度提升,这是因为药浆从大面积料斗流道进入小截面圆孔时,流动状况发生巨大变化,入口区药浆会受到拉伸作用,出现明显的流线收敛现象,进而导致药浆不同层流之间的摩擦生热加剧,使得药浆连续相分子间热运动动能增大,从而提高流体的流动性,导致黏度下降[14]

(2) 流体流经燃烧室翼片时,固体推进剂连续相的黏合剂在翼片流动剪切作用下其高分子链段发生取向,这些高分子链时而相互缠绕,时而解开,从而分子之间的平衡力受到破坏,在高剪切速率下这种破坏难以在短时间内恢复,进而破坏了药浆内部黏滞结构,从而导致层流之间的运动黏滞度下降,黏度变小[15]

此外,从计算模型的角度分析,所选用的本构方程为H-B方程,当流体超过临界剪切速率时,流动过程遵从假塑性流体幂律变化规律。除花板和翼片的其他位置的剪切应变速率很低,进一步证实药浆经真空花板浇注系统时,燃烧室中的流动是在低剪切应力和低剪切应变速率下的流动。

2.4 药浆流平性分析

将药浆离开花板喷口至充满燃烧室空间并淌平的过程称为药浆流平过程。不同时刻药浆的体积分数云图如图 9所示。

图 9 不同时刻药浆的体积分数云图 Figure 9 Volume fraction cloud chart of the slurry at different time

图 9可知,药浆在燃烧室均匀缓慢地从中间向四围展开,在真空压差和自重的综合作用力下逐渐充满燃烧室空间,这一过程与实际流动规律极为相似。药浆在4 s左右通过花板流入燃烧室,整个过程流平性良好,直到t=46 s时,燃烧室芯孔右侧出现一微小空洞,但其随着流动逐渐消失,直到燃烧室充满,孔洞的体积分数在0.5以下,可以忽略不计,基本充满。

本算例所采用药浆的屈服值为47.1 Pa,从体积分数云图可以看出该药浆的流平性很好,表明推进剂在流平过程中屈服值越小流平性越好。

3 实验验证

本实验研究对象为某固体质量分数88%的端羟基聚丁二烯(HPTB)四组元推进剂,配方(质量分数)为:黏合剂,10.5%;高氯酸铵,60%;含能炸药,10%;金属粉末,18%;其他功能助剂,1.5%。密度为1.78×103 kg/m3

药浆在混合出料后约20 min内通过落球黏度计进行黏度测定。推进剂药浆起始黏度值较低,为400 Pa·s,模拟中初始黏度的取值与其对应。1 h后药浆黏度升至759.12 Pa·s,之后每1 h测定一次,药浆黏度增长平稳,说明药浆整体性能良好。

3.1 模拟实验器浇注前准备

为尽可能真实地还原真空花板浇注实验,按照发动机燃烧室结构尺寸作如图 10的模拟实验器。内径100 mm,高200 mm,体积为0.001 57 m3。与模拟选用的燃烧室尺寸一致。

图 10 模拟实验器 Figure 10 Simulation experimental device

该模拟实验器采用6 mm厚度普通结构钢板焊接而成。装配前,首先用清理剂清理内表面,晾干后在内表面涂敷防黏剂,之后将试验器各个零件按照先底板后侧板的顺序进行装配,零部件之间加装密封垫。

3.2 模拟实验器浇注

浇注工装装配从下而上,依次为模拟实验器、套筒、花板、胶管阀和料斗等。花板采用Φ3 mm×36 mm规格,进行一次除气。

在浇注过程中,需有效控制胶管阀门开度,保证药浆下料速度,通过观察事先安装的线性速度传感器,显示流体经过3 mm孔径时,速度为0.015 8 m/s,与模拟结果误差在允许范围内。记录四组元配方保温温度为50℃,记录浇注时U型压力机的读数为101 000 Pa,浇注缸内的余压为4 000 Pa,最终浇注历程总时长为2 min。

实际浇注历程是模拟浇注时长的1.66倍,分析原因认为:(1)药浆非常黏稠,在流动过程中会在料斗和套筒等形成挂壁现象,挂壁越久越厚,而且在实际浇注过程中药浆会阻塞花板孔,延迟了总的浇注时间;(2)由于实验过程和数值模拟二者所选用花板规格的差异性,导致总浇注时间的差异。

综上,误差在一定的范围内,可以认定该模拟方法一定程度上正确反映了药浆的流动规律。

3.3 药柱质量检测

浇注工序结束后,将模拟实验器进行固化,并立即进行保温。检查其外表面无残药,并按照装配过程的逆顺序进行脱模。将药柱分解为Φ20 mm×20 mm的长条药柱,如图 11所示。分解过程中,未发现药柱存在较大孔洞等质量问题。

图 11 固体药柱的截面 Figure 11 Cross section of solid grain
4 结论

(1) 花板孔径越大,药浆流速越快。挤出胀大比随孔径长径比(L/D)减小而增大,当孔径长度一致时,孔径越大,挤出胀大比越大。

(2) 推进剂药浆为典型的屈服假塑性流体,一定剪切速率范围内,黏度随剪切速率的增大而降低,且当流经花板和翼片时,流体剪切速率最大,黏度最低。药浆黏度和剪切速率在流动过程中发生变化。

(3) 目前固体推进剂药浆的本构模型无法精确拟合,未来需补充实验数据,将本构模型精确化,如引入弹性项,建立线性黏弹性模型或微分形式和积分形式的非线性黏弹性模型,从而更好地描述药浆的流动状态。


参考文献
[1]
Waldy M, Gabl R, Seibl J, et al. Alternative methods for the implementation of trash rack losses in the 3D-numerical calculation with FLOW-3D[J]. Österreichische Wasser-und Abfallwirtschaft, 2015, 67(1): 64–69.
[2]
汪志清. 定量插管浇注工艺装置[J]. 推进技术, 1984, 5(3): 11–18.
WANG Zhi-qing. Quantitative intubation pouring process device[J]. Journal of Propulsion Technology, 1984, 5(3): 11–18.
[3]
Gao D, Morley N B, Dhir V. Numerical simulation of wavy falling film flow using VOF method[J]. Journal of Computational Physics, 2003, 192(2): 624–642. DOI:10.1016/j.jcp.2003.07.013
[4]
姜爱民, 于胜春, 邹常青. 立式捏合机混合釜内推进剂药浆混合过程数值仿真[J]. 海军航空工程学院学报, 2011, 26(4): 381–386.
JIANG Ai-min, YU Sheng-chun, ZOU Chang-qing. Numercial simulation of the mixing of composite solid propellant slurry in a kneading mixer[J]. Journal of Naval Air Engineering Institute, 2011, 26(4): 381–386. DOI:10.3969/j.issn.1673-1522.2011.04.006
[5]
Connelly Kay Robin, Jordan Valenti-James. Mixing analysis of a Newtonian fluid in a 3D planetary pin mixer[J]. Chemical Engineering Research and Design, 2008, 86: 1434–1440. DOI:10.1016/j.cherd.2008.08.023
[6]
Shimada T, Habu H, Seike Y, et al. X-ray Visualization measurement of slurry flow in solid propellant casting[J]. Flow Measurement & Instrumentation, 2007, 18(5): 235–240.
[7]
王小英, 尹欣梅, 吴倩, 等. 温度对RDX/PET/NPBA推进剂药浆流变特性的影响[J]. 火炸药学报, 2014, 37(3): 52–55.
WANG Xiao-ying, YIN Xin-mei, WU Qian, et al. Effect of temperature on rheological properties of RDX/PET/NPBA slurry[J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2014, 37(3): 52–55. DOI:10.3969/j.issn.1007-7812.2014.03.012
[8]
邓凯, 阳建红, 陈飞, 等. HTPB复合固体推进剂本构方程[J]. 宇航学报, 2010, 31(7): 1815–1818.
DENG Kai, YANG Jian-hong, CHEN Fei, et al. HTPB composite solid propellant constitutive equation[J]. Journal of Aeronautics, 2010, 31(7): 1815–1818. DOI:10.3873/j.issn.1000-1328.2010.07.018
[9]
Phan-Thien N. A nonlinear network viscoelastic model[J]. Journal of Rheology, 1987, 22: 259–283.
[10]
Phan-Thien N, Tanner R I. A new constitutive equation derived from network theory[J]. Journal of Fluid Mechanics, 1997, 2(4): 353–363.
[11]
庞爱民. 固体火箭推进剂理论与工程[M]. 北京: 中国宇航出版社, 2014.
[12]
苏昌银, 姚谦, 史旭辉, 等. 固体火箭发动机捆绑加压成型装药工艺研究[J]. 固体火箭技术, 2006, 29(6): 432–434.
SU Chang-yin, YAO Qian, SHI Xu-hui, et al. Solid rocket engine bundle research on pressure forming charge process[J]. Journal of Solid Rocket Technology, 2006, 29(6): 432–434. DOI:10.3969/j.issn.1006-2793.2006.06.010
[13]
肖立群, 樊学忠, 王晗, 等. 含铝固体推进剂燃烧过程中铝粉团聚现象研究进展[J]. 火炸药学报, 2018, 41(1): 7–15.
XIAO Li-qun, FAN Xue-zhong, WANG Han, et al. Research progress on the agglomeration phenomenon of aluminum powder in the combustion of aluminized solid propellants[J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2018, 41(1): 7–15.
[14]
李勇.黏弹性聚合物熔体流动和挤出胀大的数值模拟研究[D].上海: 华东理工大学, 2012.
LI Yong. Numercial simulations of flow and extrudate swell for viscoelastic polymer melt[D].Shanghai: East China University of Science and Technology, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10251-1012516422.htm
[15]
韩先洪.成型填充过程中非等温非牛顿黏弹性流动数值模拟[D].大连: 大连理工大学, 2007.
HAN Xian-hong. Nonisothermal non-Newtonian viscoelastic flow simulation of filling process[D].Dalian: Dalian University of Technology, 2007. http://cdmd.cnki.com.cn/Article/CDMD-10141-2007207212.htm