板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析.pdf
《板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析.pdf》由会员分享,可在线阅读,更多相关《板式表面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析.pdf(8页珍藏版)》请在咨信网上搜索。
1、中国空间科学技术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.南昌航空大学 飞行器工程
2、学院,南昌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万时,数值模拟的误差能够收敛,具有较高的仿真效率。基于该仿真模型,对比分析了贮箱在受到轴向、周向扰动时与不受扰动时推进剂重定位过程的质心坐标变化,发现在扰动
3、持续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引用格式:戴炜,张钟元,李光昱,等.板式表
4、面张力贮箱推进剂重定位过程数值仿真的网格收敛性分析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.
5、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
6、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
7、 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
8、 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
9、 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
10、 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
11、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
12、 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 引言微重力环境下,气液高度混合,需要卫星贮箱内的推进剂管理装置对推进剂进行有效管理,实现气液分离,为推进系统提供不夹气的推进剂。板式表面张力贮箱以
13、其质量轻、可靠性高的优点被广泛应用于多型号卫星1-3。在板式表面张力贮箱的设计过程中,需要对其性能进行验证。目前落塔试验有效时间短4、微重力飞机试验成本过高4-6,都难以对板式表面张力贮箱性能进行大量、长周期的验证。随着计算流体力学(C F D)和计算机技术的迅猛发展,通过数值模拟技术对卫星贮箱性能进行仿真成为热点7,8,极大提高了板式表面张力贮箱的设计和优化效率。根据计算流体力学理论可知,在合理的参数设置下,数值仿真的误差随着网格数量的增加而逐渐减小直至收敛9。长期以来,工程人员常通过提高网格密度来提高仿真结果的精确度。然而,随着网格数量的增加,所消耗的计算成本以及计算时间也增加,因此需要对
14、网格收敛性即网格无关性进行验证。在满足精度要求的情况下,较少地消耗计算资源。针对网格收敛性验证问题,早在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等通过对
15、比 实 验 亦 得 出 了 同 样 的 结 论。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指数,对不同网格数量的仿真模型进行网格无关性分析选出兼顾计算速度与精度的仿真模型;最后
16、基于该模型计算不同方向扰动下贮箱质心的变化,以检验其抗外部扰动的能力。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
17、 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
18、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通过调整网格
19、密度及边界层的数量(如图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
20、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模型。不考虑液体流动过程中的能量传递和热交换。设置压力-速
21、度耦合方程求解算法为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+g
22、1h+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
23、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 仿真结果不同于地面环境,卫星运行在微重力环境下,表面张力占主导地位。贮箱内的推进剂在表面张力的作用下,将沿导流板及贮箱壁面形成的
24、内角流动。当贮箱推进剂达到最小势能状态时,气液交界面按一定的曲率半径附着在导流板板壁周围。此时贮箱内液体呈连续分布,气体被包裹在贮箱顶部。当贮箱受到扰动时,例如卫星调姿时,贮箱气液界面被破坏,当扰动停止时,推进剂在表面张力的作用下进行重定位,恢复到稳定的气液界面。为了保证卫星两次调姿时推进系统正常运行,需要获得贮箱推进剂在扰动停止后的重定位时间。图4的数据来源于第二套模型(网格数为6 2 5万的模型),显示了贮箱6 0%填充率时重定位过程的仿真情况。图4(a)为仿真初始状态下气液交界面的形状,当重力消失,推进剂在表面张力的作用下沿罐壁和导流板爬升,气液交界面由平面逐渐收缩成弯液面。图4(d)为
25、重定位完成时气液交界面的形状。图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左右。为保证仿真过程中流场内的计算精度大致相同,要求网格的疏密程度
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 板式 表面张力 推进 定位 过程 数值 仿真 网格 收敛性 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。