低重环境航天器贮箱内三维液体晃动数值模拟.pdfVIP

全国第十一届空间及运动体控制技术学术年会论文 低重环境航天器贮箱内 三维液体晃动数值模拟 周宏李俊峰王天舒 清华大学工程力学系，北京100084 摘要： 低重环境下液体的晃动具有与常重环境不同的动力学特性．该文采用ALE 描述对煎重环境下航天嚣贮箱内带自由液面边界的粘性不可压液体运动进行了有限 垂坌堑．选用分步法求解Navier—s矗五i。茅石，仅需用同阶插值函数离散控制方程即 可满足数值计算的精度要求；引入最小二乘数值逼近方法拟合自由液面平均曲率，给 出了舍表面张力的动力学边界条件的数值计算形式．模拟得到贮箱壁面位置随时间变 化的历程和自由液面晃动的空问区域，获得了低重环境下的液体晃动不同于常重晃动 下的主要特征． 关键词：航天器，ALE方法，表面张力，低重力 在航天器姿态控制精度要求更加严格，充液量增大的工程背景下，研究低重环境 (10w—gravity environment)的充液贮箱内液体的大幅晃动，由于更接近航天器在外太空的工 作环境，而显现出了更加重要工程应用价值。采用数值方法研究低重环境的液体晃动，能弥补 使用实验研究实验机会少，外载环境实现困难，实验初始条件不易保证等条件的限制…，同时 也能避免理论研究中建立数理模型多重假定条件的限制，以及研究范围受限等影响。 针对液体大幅晃动中液体所占区域会发生大范围的运动和变形，采用任意拉格朗日欧拉描 述方式(ArbitraryLagrange—EulerMethods，简称ALE)对液体的运动进行描述，能很好的 克服单独使用以物质域变量为描述参数的Lagrange描述格式，在处理流场区域大范围变形的网 格畸变导致计算不能进行的情况。同时也避免了在以空间域变量为描述参数的Euler描述格式 有强烈的非线性对流项给求解方程带来的困难，同时又兼有这两种运动描述方式的优点”””。 行了离散，在液体不可压的假设下，采用了分步法处理有限元离散方程“’，使得可以将液体速 度变量和液体压力变量采用同阶插值方案离散。 另外，在低重环境下，表面张力效应对于液体晃动响应的影响将远大于环境重力对液体运 动的影响。研究低重环境下的晃动问题就必须处理含有表面张力作用效应的动力学边界条件， 这是数值模拟低重环境下液体晃动的一个技术难点。本文提出了用最4、乘数值逼近方法拟合 自由液面平均曲率的方法，将表面张力效应引入自由液面边界条件的方法模拟低重环境的液体 晃动。 1 ALE描述下的Navier-Stokes基本变量方程 在ALE描述下任意物理量f的隧体导数可以写为： 堂dt一翌I+q丢 (f-1】2，3) (1) c，为』￡F参考坐标下的对流速度，是物质速度叶和网格速度wJ之差： 32 全盈第十一届空问及运动体控制技术学术年会论文 Cj=叶一叶 (2) 不可压粘性流体的Navier-Stokes基本变量方程和连续方程，在ALE描述格式下可以写为： 』+Cj等一≯1·州％，M。一厶 ㈤ IHu=0， 其中虬为流体质点的运动速度，Pi为流场压力，f，J为质点的自由度数，P，v分别是流体的密 度和运动粘性系数，，是流体所受的外力。 现假设自由面的曲面方程为F(x，f)=0，在ALE格式描述下，其随体导数为 u：l，2，3)u21，2，川 性’(4) 譬：iOF+Cj_OF：o i2百i。o 为保证自由液面边界上的网格点始终在自由液面上．(2)式中的当地导数项应该为零，所以在 自由面上网格点的速度变量应满足： (Uj--M)豢_0’㈣，2，3)．