板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析.pdf

中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.4 8 5-9 2C h i n e s e S p a c e S c i e n c e a n d T e c h n o l o g yI S S N1 0 0 0-7 5 8 X C N1 1-1 8 5 9/Vh t t p:z g k j.c a s t.c nD O I:1 0.1 6 7 0 8/j.c n k i.1 0 0 0-7 5 8 X.2 0 2 3.0 0 5 6板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析戴炜1,张钟元1,李光昱1,韩伟2,*1.南昌航空大学 飞行器工程学院,南昌3 3 0 0 6 32.军事科学院 国防科技创新研究院,北京1 0 0 0 7 1摘 要:为了验证2.6 7L卫星贮箱性能,通过建立贮箱仿真模型,分别对贮箱5 0%与6 0%填充率下推进剂的重定位过程进行仿真验证。为了提高仿真效率,通过引入G C I指数这一概念,对比2 9 5万、6 2 5万、7 6 7万、9 5 5万网格数量的仿真结果,对仿真模型的网格收敛性进行分析。通过两个工况仿真结果的对比验证,认为当网格数量为6 2 5万时,数值模拟的误差能够收敛,具有较高的仿真效率。基于该仿真模型,对比分析了贮箱在受到轴向、周向扰动时与不受扰动时推进剂重定位过程的质心坐标变化,发现在扰动持续1 s情况下,推进剂质心变化小于11 0-3m,因此认为该板式贮箱具有较好的抗扰动能力。关键词:卫星贮箱;仿真模型;C F D;重定位过程;G C I指数中图分类号:V 4 1 1.8 文献标识码:A收稿日期:2 0 2 2-0 4-1 1;修回日期:2 0 2 2-0 6-2 9;录用日期:2 0 2 2-0 9-0 2;网络出版时间:2 0 2 2-0 9-2 3 1 0:4 9基金项目:国家自然科学基金(5 2 0 6 5 0 4 4)*通信作者.E-m a i l:h a n w e i 1 9 8 4h o t m a i l.c o m引用格式:戴炜,张钟元,李光昱,等.板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析J.中国空间科学技术,2 0 2 3,4 3(4):8 5-9 2.D A I W,Z HAN G Z Y,L I G Y,e t a l.G r i d c o n v e r g e n c e a n a l y s i s f o r n u m e r i c a l s i m u l a t i o n o f p r o p e l l a n t r e p o s i t i o n p r o c e s s i n v a n e-t y p e s u r f a c e t e n s i o n t a n kJ.C h i n e s e S p a c e S c i e n c e a n d T e c h n o l o g y,2 0 2 3,4 3(4):8 5-9 2(i n C h i n e s e).G r i d c o n v e r g e n c e a n a l y s i s f o r n u m e r i c a l s i m u l a t i o n o f p r o p e l l a n t r e p o s i t i o n p r o c e s s i n v a n e-t y p e s u r f a c e t e n s i o n t a n kD A I W e i1,Z H A N G Z h o n g y u a n1,L I G u a n g y u1,H A N W e i2,*1.S c h o o l o f A i r c r a f t E n g i n e e r i n g,N a n c h a n g H a n g k o n g U n i v e r s i t y,N a n c h a n g 3 3 0 0 6 3,C h i n a2.N a t i o n a l I n n o v a t i o n I n s t i t u t e o f D e f e n s e T e c h n o l o g y,C h i n e s e A c a d e m y o f M i l i t a r y S c i e n c e,B e i j i n g 1 0 0 0 7 1,C h i n aA b s t r a c t:I n o r d e r t o v e r i f y t h e p e r f o r m a n c e o f 2.6 7L s a t e l l i t e t a n k,a t a n k s i m u l a t i o n m o d e l w a s e s t a b l i s h e d.T h e r e p o s i t i o n p r o c e s s e s o f p r o p e l l a n t a t 5 0%a n d 6 0%f i l l i n g r a t e s w e r e s i m u l a t e d.F o r i m p r o v i n g t h e s i m u l a t i o n e f f i c i e n c y,t h e c o n c e p t o f G C I i n d e x w a s i n t r o d u c e d.T h e s i m u l a t i o n r e s u l t s o f 2.9 5m i l l i o n,6.2 5m i l l i o n,7.6 7m i l l i o n a n d 9.5 5m i l l i o n g r i d s w e r e c o m p a r e d t o a n a l y z e t h e g r i d c o n v e r g e n c e o f t h e s i m u l a t i o n m o d e l.T h r o u g h t h e c o m p a r i s o n a n d v e r i f i c a t i o n o f t h e s i m u l a t i o n r e s u l t s,i t w a s c o n s i d e r e d t h a t t h e e r r o r o f n u m e r i c a l s i m u l a t i o n c o u l d c o n v e r g e a n d h a d h i g h s i m u l a t i o n e f f i c i e n c y w h e n t h e n u m b e r o f g r i d s w a s 6.2 5m i l l i o n.B a s e d o n t h e r e s u l t s,t h e c a p a b i l i t y o f t h e t a n k p r o p e l l a n t m a n a g e m e n t d e v i c e t o s u p p r e s s l i q u i d s l o s h i n g w a s a n a l y z e d.A n d i t i s s h o w n t h a t t h e 2.6 7L s a t e l l i t e t a n k h a s g o o d p e r f o r m a n c e i n a n t i-d i s t u r b a n c e c a p a b i l i t y.K e y w o r d s:s a t e l l i t e t a n k;s i m u l a t i o n m o d e l;C F D;r e p o s i t i o n i n g p r o c e s s;G C I i n d e x8 6 中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.41 引言微重力环境下,气液高度混合,需要卫星贮箱内的推进剂管理装置对推进剂进行有效管理,实现气液分离,为推进系统提供不夹气的推进剂。板式表面张力贮箱以其质量轻、可靠性高的优点被广泛应用于多型号卫星1-3。在板式表面张力贮箱的设计过程中,需要对其性能进行验证。目前落塔试验有效时间短4、微重力飞机试验成本过高4-6,都难以对板式表面张力贮箱性能进行大量、长周期的验证。随着计算流体力学(C F D)和计算机技术的迅猛发展,通过数值模拟技术对卫星贮箱性能进行仿真成为热点7,8,极大提高了板式表面张力贮箱的设计和优化效率。根据计算流体力学理论可知,在合理的参数设置下,数值仿真的误差随着网格数量的增加而逐渐减小直至收敛9。长期以来,工程人员常通过提高网格密度来提高仿真结果的精确度。然而,随着网格数量的增加,所消耗的计算成本以及计算时间也增加,因此需要对网格收敛性即网格无关性进行验证。在满足精度要求的情况下,较少地消耗计算资源。针对网格收敛性验证问题,早在1 9 9 4年,R o a c h e1 0-1 1基于R i c h a r d s o n外推法,提出了网格收敛指数(G r i d C o n v e r g e n c e I n d e x,G C I)用以评估网格模型。郑秋亚1 2等通过将网格逐步细分或粗化,采用G C I指数对网格无关性进行评估,选择一个理想的网格,做到既经济又能达到一定的精度要求。刘厚林1 3等通过对比得出G C I指数随着网格的加密会不同程度的减小,且由于网格越密离散结果越接近精确解。安恩科1 4等通过对比 实 验 亦 得 出 了 同 样 的 结 论。C h u n g1 5等通过G C I指数分析液氢贮箱模型的收敛情况,获得了效率较高的仿真模型。以上研究均表明通过G C I指数能够较好地验证C F D数值仿真中网格的收敛性。目前部分仿真研究,对仿真模型并未进行网格收敛性分析,仿真效率较低。本文以板式表面张力贮箱的应用为背景,首先建立了2.6 7L板式表面张力贮箱的几何模型,并对该模型离散化,生成了不同数量网格;然后对贮箱在5 0%、6 0%填充率时的重定位过程进行了数值仿真。基于仿真结果,通过计算G C I指数,对不同网格数量的仿真模型进行网格无关性分析选出兼顾计算速度与精度的仿真模型;最后基于该模型计算不同方向扰动下贮箱质心的变化,以检验其抗外部扰动的能力。2 板式表面张力贮箱仿真模型2.1 物理模型贮箱模型为球冠加圆柱段箱体(如图1所示),球面直径为1 2 0mm,圆柱段长度1 6 0mm,总高度为2 8 0mm,容积为2.6 7L。箱体内部安装有6片内导流板和6个外导流板,使得该贮箱具有较强的推进剂蓄留及导引能力1 6。内外导流板成3 0 夹角均匀分布在中央支撑柱周围。贮箱底部有2 4片叶片组成的蓄液器,用于卫星机动过 程 中 蓄 留 推 进 剂,为 推 进 系 统 提 供 推进剂。仿真中贮箱内选择的工质为液态水,贮箱顶部留有气垫,选择工质为空气,其密度为1.2 2 5 k g/m3,粘度为1.7 8 9 41 0-5k g/ms。液态水的相关参数如表1。图1 2.6 7L板式贮箱模型F i g.1 2.6 7L v a n e-t y p e t a n k m o d e l表1 工质的相关参数T a b l e 1 R e l a t e d p a r a m e t e r s o f w o r k i n g m e d i u mP a r a m e t e r sV a l u e sW o r k i n g m e d i u mW a t e r(l i q u i d)D e n s i t y/(k gm-3)9 9 8.2V i s c o s t y/(k gm-1s)0.0 0 1 0 0 3S u r f a c e t e n s i o n c o e f f i c i e n t/(Nm-1)0.0 6 6 7C o n t a c t a n g l e/()5戴炜,等:板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析8 7 2.2 网格划分根据贮箱物理模型,建立仿真分析模型。采用非结构四面体网格划分方式(如图2所示)对贮箱内部流域进行离散化。图2 2.6 7L板式贮箱网格模型F i g.2 G r i d m o d e l o f 2.6 7L p l a t e t a n k通过调整网格密度及边界层的数量(如图3所示),得到网格密度不同的四个模型。如表2所示图3 4种不同密度的模型在X Y平面(z=0)的体网格截图F i g.3 V o l u m e m e s h s c r e e n s h o t s o f f o u r m o d e l s w i t h d i f f e r e n t d e n s i t i e s i n X Y p l a n e(z=0)表2 4组网格的基本参数T a b l e 2 B a s i c p a r a m e t e r s o f 4 g r o u p s o f g r i d sN u n m b e r o f g r i d sN u m b e r o f n o d e s2 9 5 0 0 0 06 0 6 2 5 96 2 5 0 0 0 01 1 2 0 0 4 87 6 7 0 0 0 01 3 6 4 1 8 89 5 5 0 0 0 01 6 6 9 8 5 32.3 仿真参数设置为模拟流体在贮箱内的流动情况,设置为非定常流动,不考虑重力的影响。各个壁面均定义为无滑移边界。选用VO F模型,设定初始状态下液体保持在贮箱底部,气体保持在贮箱顶部。将空气作为基本项,第二相为液相。湍流模型为L a m i n a r模型。不考虑液体流动过程中的能量传递和热交换。设置压力-速度耦合方程求解算法为P I S O方 式,梯 度 插 值 方 案 采 用G r e e n-G a u s e C e l l B a s e d法,压 力 插 值 算 法 选 择B o d y F o r c e W e i g h t e d,动量方程选用QU I C K方法。本文分别按照推进剂5 0%和6 0%的填充比,设定各相组分的初始条件,建立气液交界面并对其运动过程进行监测。监控并记录气液交界面最高点的位置,以推进剂爬升的最高点为仿真结果分析仿真模型的收敛情况。3 网格收敛指数假设数值仿真过程中存在仿真收敛解f与真实解fe x a c t。则f可表示为:f=fe x a c t+g1h+g2h2+g3h3+式中:h为网格间距,gi为与步长无关的常数。令g1=0,得到:f=fe x a c t+(h22f1-h21f2)(h22-h21)+o(h2)令网格间距比r=h2/h1,hi为不同网格数模型对应的网格间距。又因在同一组对比网格中,模型未发生改变,模型总体积不变,故:r=h2h1=3N2N18 8 中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.4fe x a c tf1+(f1-f2)/(r2-1)=f2-f1f1E=1-rp式中:Ni为节点数;p为收敛率;为收敛误差;fi为不同网格模型得到的收敛解。基于以上计算,R o a c h e引入了一个安全系数Fs,通常Fs被设置为1.2 53之间1 7。定义网格收敛指数为:G C I=Fs|E|由于本次仿真网格数量采用非定常数量增加,则其收敛率p可由以下方程解出:2 3rp2 3-1=rp1 21 2rp1 2-1(1)式中:1 2、2 3为选定考察项目的网格收敛误差;r为网格增长率,r的下标为3套不同的网格,编号分别为1、2、3;r1 2、r2 3为对应网格数量增长率。则网格收敛指数(G C I)定义为:G C I=Fsrp-14 讨论与分析4.1 仿真结果不同于地面环境,卫星运行在微重力环境下,表面张力占主导地位。贮箱内的推进剂在表面张力的作用下,将沿导流板及贮箱壁面形成的内角流动。当贮箱推进剂达到最小势能状态时,气液交界面按一定的曲率半径附着在导流板板壁周围。此时贮箱内液体呈连续分布,气体被包裹在贮箱顶部。当贮箱受到扰动时,例如卫星调姿时,贮箱气液界面被破坏,当扰动停止时,推进剂在表面张力的作用下进行重定位,恢复到稳定的气液界面。为了保证卫星两次调姿时推进系统正常运行,需要获得贮箱推进剂在扰动停止后的重定位时间。图4的数据来源于第二套模型(网格数为6 2 5万的模型),显示了贮箱6 0%填充率时重定位过程的仿真情况。图4(a)为仿真初始状态下气液交界面的形状,当重力消失,推进剂在表面张力的作用下沿罐壁和导流板爬升,气液交界面由平面逐渐收缩成弯液面。图4(d)为重定位完成时气液交界面的形状。图4 6 0%填充率时液面爬升趋势图F i g.4 L i q u i d c l i m b i n g t r e n d a t 6 0%f i l l i n g r a t e4.2 不同网格数量结果汇总4组不同网格的计算结果,绘制推进剂液面最高点的轨迹图。如图5所示,流动趋势保持一致,气液分界面的最高点相差较小。相较于网格数为2 9 5万的仿真模型,其余3组模型的重定位时间和液面高度曲线与文献1 6 中的理论计算结 果 基 本 吻 合 且 重 定 位 时 间 相 同,均 为0.7 5 s左右。为保证仿真过程中流场内的计算精度大致相同,要求网格的疏密程度相合理。对比图2网格剖面,网格数为2 9 5万的模型在贮箱边缘附近网格单元质量较差,且各导流板间的流域内网格疏密程度和排列方式均不相同。通过对数据后处理分析,该模型难以对接近壁面处推进剂的流动情况进行较高精度的仿真。图5 6 0%填充率下4种网格模型液面爬升过程对比F i g.5 C o m p a r i s o n o f l i q u i d c l i m b i n g p r o c e s s o f f o u r g r i d m o d e l s u n d e r 6 0%f i l l i n g r a t e通过观察图4中推进剂重定位过程中的气液交界面知,贮箱导流板具有良好的导流能力,戴炜,等:板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析8 9 且推进剂的最高点在内导流板的根部附近。为评估仿真模型对计算结果的影响,对仿真结果进行网格收敛性的分析。在评估过程中以气液分界面最高点的值作为对比量,为消除仿真计算过程中舍入误差的影响,每组工况均重复计算3次,分析数据并对差别过大数据予以剔除并求取计算平均值。汇总相关计算数据。如表3所示:表3 不同填充率下仿真结果T a b l e 3 S i m u l a t i o n r e s u l t s u n d e r d i f f e r e n t f i l l i n g r a t e sN o.T o t a l g r i d s T o t a l n o d e sL1/m L2/m12 9 5 0 0 0 06 0 6 2 5 90.0 7 8 3 0 0.0 7 0 9 526 2 5 0 0 0 01 1 2 0 0 4 80.0 7 8 7 4 0.0 7 3 9 937 6 7 0 0 0 01 3 6 4 1 8 80.0 7 9 0 8 0.0 7 5 0 449 5 5 0 0 0 01 6 6 9 8 5 30.0 7 9 9 4 0.0 7 3 1 1 其中L1为6 0%填充率液面最高点,L2为5 0%填充率液面最高点。4.3 网格无关性对比分析为检验贮箱重定位模型的网格收敛性,基于表3得到的仿真结果进行计算。由网格收敛指数定义知,收敛率p与网格模型有关。因为在建模时采用非结构化网格的网格划分方式,故网格间距的比值r为非定常量。因此在按式(1)计算收敛率p时,无法统一计算收敛率。所以将4组模型将数据按编号分为(1,2,3)、(2,3,4)、(1,2,4)和(1,3,4)共计4组对比组。R o a c h e 1 6指出,当对比组内存在3项数值时,此时安全系数Fs=1.2 5,带入公式(2)分别计算每组网格的网格收敛指数G C If i n e 和G C Ic o a r s e。计算结果如表4所示。表4 6 0%填充率G C I指数计算结果T a b l e 4 C a l c u l a t i o n r e s u l t s o f G C I i n d e x o f 6 0%f i l l i n g r a t eG r o u p1 22 3r1 2r2 3pG C If i n eG C Ic o a r s e1,2,30.0 0 0 4 40.0 0 0 3 41.2 2 71.0 6 87.0 4 90.0 1 70.0 7 22,3,40.0 0 0 3 40.0 0 0 8 61.0 6 81.0 6 91 3.7 80.0 2 9 80.0 7 1 31,3,40.0 0 0 7 80.0 0 0 8 61.3 1 01.0 6 91 0.6 90.0 0 5 80.1 0 31,2,40.0 0 0 4 40.0 0 1 21.2 2 71.1 4 28.9 7 40.0 1 0 40.0 6 5 4 通过对比四组对比组的数据,当加注量在6 0%时四组网格的网格收敛指数(G C I)均在1%以内,符合小于3%的要求1 4。且随着网格数量的增加,网格收敛指数呈下降的趋势,且网格收敛指数趋于0。A l i等1 8在R o a c h e的基础上,定义了收敛比R的概念:R=f2-f1f3-f2 若0R1,则网格单调收敛,R1则网格发散。经计算,当充液量为6 0%时第1对比组的网格收敛比R1,其余3个对比组的网格收敛比均在区间(0,1)内。符合收敛条件。结合G C I指数计算结果。判定由网格数量分别为6 2 5万、7 6 7万和9 5 5万的三组网格组成的网格组收敛。为验证上述结论,对5 0%填充率下,不同网格数量的仿真结果进行分析。如表5所示。在5 0%加注量时,网格数为2 9 5万的模型的计算结果相较剩余3组模型有较大偏差,使得其在网格收敛指数计算结果偏离正常值,甚至对比组(1,2,4)无正数解。且当充液量为5 0%时第1、3、4对比组的网格收敛比R1,仅有第2对比组的网格收敛比在区间(0,1)内。符合收敛条件。综合以上计算,当充液量为5 0%时,仅有第二对比组的网格符合收敛条件。从而得到网格数量分别为6 2 5万、7 6 7万和9 5 5万的三组网格组成的网格组收敛。表5 5 0%填充率G C I指数计算结果T a b l e 5 C a l c u l a t i o n r e s u l t s o f G C I i n d e x o f 5 0%f i l l i n g r a t eG r o u p1 22 3r1 2r2 3pG C If i n eG C Ic o a r s e1,2,30.0 0 3 0 40.0 0 1 0 51.2 2 71.0 6 80.5 3 23.3 0 53.6 8 52,3,40.0 0 1 0 50.0 0 1 9 31.0 6 81.0 6 98.9 6 71.6 3 32.9 51,3,40.0 0 4 0 90.0 0 1 9 31.3 1 01.0 6 94.1 2 80.7 6 12.4 9 61,2,40.0 0 3 0 40.0 0 0 8 81.2 2 71.1 4 2-4.6 69 0 中国空间科学技术A u g.2 5 2 0 2 3 V o l.4 3 N o.4 通过比较4组不同数量网格的仿真模型,结果表明,网格数量分别为6 2 5万、7 6 7万和9 5 5万的仿真模型仿真收敛,能够满足目标精度。其中数量为6 2 5万时,不同工况的仿真结果已具有较高精度。经对比,网格数为6 2 5万、7 6 7万和9 5 5万的仿真模型计算耗时分别为8小时、1 0小时和1 3小时。综合考虑网格数为6 2 5万的模型计算时间较少,具有较高性价比。4.4 贮箱受扰动情况下的重定位过程卫星在轨运行时,根据其任务不同,需要进行在轨机动及轨道保持。为验证贮箱在轨运行时 的 抗 扰 动 能 力,针 对 卫 星 脉 冲 的 典 型 工况1 9,2 0,分别对贮箱5 0%填充率时受到周向和轴向扰动的质心晃动情况进行仿真分析。由于贮箱为轴对称结构,故仅模拟沿-X方向及-Z方向受到外部扰动的情况。设置扰动时长为1 s,扰动大小为1.41 0-2m/s2。分析并记录扰动过程整体质心的变化,用以对比不同扰动对重定位过程的影响,以检验贮箱抗扰动的能力。图6显示了贮箱质心坐标在不同扰动下的变化情况。当受到-X方向扰动时,X轴质心偏移量随时间逐渐加大,最大值为71 0-4m,此时Y轴质心偏移量经历小幅度波动后,逐渐恢复基本与无扰动时相一致。对比Z方向的质心坐标,无论是受到-X向还是-Z向扰动,其坐标 曲 线 基 本 重 合,最 大 偏 移 量 仅 为51 0-5m。如表6所示。由图可知,贮箱周向扰动对贮箱质心的影响较轴向扰动更大,但总体变化小于1 0-3m,因此单方向的扰动对其他方向的质心坐标影响亦十分有限。通过与无扰动时对比,在施加外部扰动下,质心各坐标轴的变化较小,小于1 0-3m。因而认为该贮箱有较好的抗扰动能力。该结果符合设计时同时采用内外导流板的初衷。表6 各坐标轴重心最大偏移量T a b l e 6 M a x i m u m c e n t r o i d d e v i a t i o n o f e a c h a x i sD i r e c t i o n-X s l o s h i n g-Z s l o s h i n gX a x i s/m71 0-421 0-5Y a x i s/m2.51 0-521 0-5Z a x i s/m51 0-54.51 0-5图6 质心坐标在不同扰动下的对比图F i g.6 C o m p a r i s o n o f c e n t r o i d i n a x i s u n d e r d i f f e r e n t d i s t u r b a n c e5 结论论文以板式表面张力贮箱的典型工况仿真为背景,建立了2.6 7L板式表面张力贮箱的几何模型以及不同网格数量的仿真模型,通过计算不同贮箱模型在填充率为5 0%和6 0%时重定位仿真结果的G C I指数,进行网格收敛性分析。结果表明:在2.6 7L贮箱推进剂重定位仿真过程中,对于网格数为6 2 5万、7 6 7万和9 5 5万的戴炜,等:板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析9 1 仿真模型对比组,其数值模拟的误差逐渐收敛,数值模拟的计算值结果不会因网格数量不足而失真。综合考虑各方面的因素,选择网格数为6 2 5万的网格具有较高计算效率,用于板式贮箱后续工况的仿真计算。基于该仿真模型,分别对贮箱受到轴向、周向1.41 0-2 m/s2加速度扰动以及不受扰动情况下的质心变化进行分析,质心变化小于11 0-3m,认为该贮箱具有较好的抗扰动能力。通过监测外加扰动下贮箱质心的变化,得到在加速度为1.41 0-2 m/s2的外加扰动下,贮箱内的液面管理装置能有效的抑制液面晃动。参考文献(R e f e r e n c e s)1 T AM W,B A L L I N G E R I,J A E K L E D.S u r f a c e t e n s i o n PMD t a n k f o r o n o r b i t f l u i d t r a n s f e rC 4 4 t h A I AA/A S ME/S A E/A S E E J o i n t P r o p u l s i o n C o n f e r e n c e&E x h i b i t.H a r t f o r d,C o n n e c t i c u t:A I AA,2 0 0 8:5 1 0 5.2 N E T T E R G,R E NN E R U,G E R S TMANN J,e t a l.D e s i g n o f a n e w r e f i l l a b l e r e s e r v o i r f o r a s t a n d a r d s u r f a c e t e n s i o n t a n kC 3 6 t h A I AA/A S ME/S A E/A S E E J o i n t P r o p u l s i o n C o n f e r e n c e a n d E x h i b i t.H u n t s v i l l e,A l a b a m a:A I AA,2 0 0 0:3 4 4 6.3 T AM W,HE R S H M,B A L L I N G E R I.H y b r i d p r o p e l l a n t t a n k s f o r s p a c e c r a f t a n d l a u n c h v e h i c l e sC 3 9 t h A I AA/A S ME/S A E/A S E E J o i n t P r o p u l s i o n C o n f e r e n c e a n d E x h i b i t.H u n t s v i l l e,A l a b a m a:A I AA,2 0 0 3:4 6 0 7.4 魏延明,潘海林.全管理圆柱形表面张力贮箱的微重力实验验证 静平衡 与 重 定 位J.控 制 工 程,1 9 9 7,(5):1 4-1 9.5 苗楠,刘战合,王晓璐,等.飞机抛物线飞行试验中的液体行为仿真J.飞行力学,2 0 2 0,3 8(2):8 1-8 7.M I AO N,L I U Z H,WANG X L,e t a l.S i m u l a t i o n o f l i q u i d b e h a v i o r i n a i r c r a f t p a r a b o l i c f l i g h t t e s t J.F l u i d M a c h i n e r y,2 0 2 0,3 8(2):8 1-8 7(i n C h i n e s e).6 张晨辉,芮伟,段俐,等.基于板式表面张力贮箱的落塔实验C第九届全国实验流体力学学术会议.杭州,2 0 1 3.Z HAN G C H,RU I W,D UAN L.e t a l.T h e t o w e r d r o p e x p e r i m e n t s b a s e d o n t h e v a n e t y p e s u r f a c e t e n s i o n t a n kC T h e 9 t h N a t i o n a l C o n f e r e n c e o n E x p e r i m e n t a l F l u i d M e c h a n i c s.H a n g z h o u,2 0 1 3(i n C h i n e s e).7 胡齐,李永,梁军强,等.微重力下板式贮箱内液体晃动性能研究J.航天控制,2 0 1 7,3 5(5):1 9-2 4.HU Q,L I Y,L I ANG J Q,e t a l.L i q u i d s l o s h i n g p e r f o r m a n c e i n v a n e t y p e t a n k u n d e r m i c r o g r a v i t yJ.A e r o s p a c e C o n t r o l,2 0 1 7,3 5(5):1 9-2 4(i n C h i n e s e).8 胡齐,李永,潘海林,等.微重力环境下大叶片板式贮箱内流体行为的数值仿真与试验验证J.空间控制技术与应用,2 0 1 3,3 9(2):5 8-6 2.HU Q,L I Y,P AN H L,e t a l.N u m e r i c a l s i m u l a t i o n a n d e x p e r i m e n t v e r i f i c a t i o n o f f l u i d b e h a v i o r i n t h e v a n e t y p e t a n k w i t h b i g v a n e s i n m i c r o g r a v i t y e n v i r o n m e n t J.A e r o s p a c e C o n t r o l a n d A p p l i c a t i o n,2 0 1 3,3 9(2):5 8-6 2(i n C h i n e s e).9 张德良.计算流体力学教程M.北京:高等教育出版社,2 0 1 0:9 0-9 2.Z HANG D L.A c o u r s e i n c o m p u t a t i o n a l f l u i d d y n a m i c s M.B e i j i n g:H i g h e r E d u c a t i o n P r e s s,2 0 1 0:9 0-9 2(i n C h i n e s e).1 0 R OA C HE P J.P e r s p e c t i v e:a m e t h o d f o r u n i f o r m r e p o r t i n g o f g r i d r e f i n e m e n t s t u d i e sJ.J o u r n a l o f F l u i d s E n g i n e e r i n g,1 9 9 4,1 1 6(3):4 0 5-4 1 3.1 1 R OA C HE P J.Q u a n t i f i c a t i o n o f u n c e r t a i n t y i n c o m p u t a t i o n a l f l u i d d y n a m i c s J.A n n u a l R e v i e w o f F l u i d M e c h a n i c s,1 9 9 7,2 9(1):1 2 3-1 6 0.1 2 郑秋亚,王宝圆.网格密度对流场解的影响J.航空计算技术,2 0 0 0,(4):1 3-1 6,2 2.Z HE N G Q Y,WAN G B Y.T h e i n f l u e n c e o f g r i d d e n s i t y t o f l o w f i e l d s o l u a t i o n J.A e r o n a u t i c a l C o m p u t i n g T e c h n i q u e,2 0 0 0,(4):1 3-1 6,2 2(i n C h i n e s e).1 3 刘厚林,刘明明,白羽,等.基于G C I的离心泵网格的收敛性J.江 苏 大 学 学 报(自 然 科 学