利用谱方法数值求解三维球坐标系下的Navier-Stokes方程

从三维球坐标系下的Navier-Stokes方程出发,采用谱配置方法建立了一个三维球坐标系下中层大气波动传播的全非线性动力学数值模式;构造了合理的纬向边界,并利用初始给定的重力波波包对该模式进行了数值检验．计算结果展示了重力波波包在大气中传播的三维图像,波包整体向上传播,波相关扰动速度随高度的增加而增长．     

软管锥套式空中加油系统建模与特性分析

针对自主空中加油研究中软管锥套运动模型过于简化,现有建模方法又存在成本高、计算量大、软管长度恒定等缺陷的问题,根据集中参数法原理,提出了一种长度可变的多级串联理想单摆系软管锥套运动模型.同时考虑软管收放、加油机牵连运动、重力、定常流、大气扰动、加油机尾流等内外部因素,推导了迭代形式的变长度软管锥套三维运动方程.由摆长约束导出了求解软管拉力的代数线性方程组,进而给出了模型稳定性证明和适用条件.通过数值仿真,测试了锥套阻力、平衡位置、软管收放等稳态特性,分析了加油机滚转运动、尾涡流场对软管锥套的动态影响以及软管甩鞭现象的产生机理.数值仿真结果验证了所建模型的有效性.      

三维球坐标系下波包非线性传播过程的一种数值求解方法

采用全隐欧拉（FICE）格式，将三维坐标系下的大气基本运动方程作为控制方程，建立起了三维球坐标系下可压大气中的全非线性数值模式；对纬向边界也进行了合理的处理，并对该模式进行了数值检验，计算结果显示出重力波波包在大气中非线性传播的三维图像：波包整体向上传播，波相关扰动速度随高度增加而增长。     

飞机连投货物实时仿真软件及应用

为了降低实时仿真对数字机性能的要求，提出了一个简化数学模型及程序结构的实时仿真软件，验证了可行性，所得结果与全量方程数字计算及实际一致，较好地解决了实时性和精度的要求，表明软件设计合理，性能良好，可用于飞机人机闭环连投飞行操稳特性和安全性的仿真研究和检验。     

某高空螺旋桨气动特性数值模拟与风洞试验

以某高空螺旋桨为研究对象,利用k-ε湍流模型,采用三维Navier-Stokes方程,并通过数值计算模拟该螺旋桨的复杂绕流,分析其气动特性.按照等雷诺数等前进比准则,进行了该螺旋桨地面缩比风洞试验.经过比较发现,计算结果与试验结果吻合良好.验证了数值模拟计算高空螺旋桨气动特性的正确性,为高空螺旋桨的设计和试验提供参考依据.     

人体温度场的边界元分析

针对人体温度场模型复杂,计算难度很大,需要简化计算,采用了边界元法.利用边界元法是以边界点上未知函数满足积分方程为基础进行近似计算的这一特点,对人体温度场模型进行了三维计算,并给出结果,与其它方法和实验比较,吻合较好.     

数字闭环光纤陀螺频率特性的计算和测试方法

建立数字闭环光纤陀螺FOG(fiber-optics gyroscopes)的数学模型,得到描写模型的差分方程,推导出与之对应的微分方程和传递函数.分析传递函数与FOG模型的误差,得出数字闭环FOG频率特性的计算方法.研究在普通速率转台测试数字闭环FOG频率特性的方法.并将测试数据与计算数据加以比较验证.计算和测试的结果证明了该计算的测试方法是准确和有效的.     

模拟三维重力波传播过程的并行数值模式

从三维大气运动的基本控制方程出发, 基于MPI消息传递接口和区域分解的思想, 建立了模拟三维可压缩大气中重力波传播过程的并行数值模式. 在对垂直方向进行区域分解的基础上, 针对跳点网格的特点, 以垂直速度主格点线为子区域的下边界, 状态变量主格点线为子区域的上边界进行区域分解. 利用MPI消息传递接口来传递计算各子区域中变量时所要用到的相邻子区域上相应变量, 从而顺利地完成对整个计算域的并行计算. 根据线性重力波理论, 通过模拟小振幅重力波的传播过程对所建模式进行了验证. 结果表明, 模式可以很好地模拟小振幅重力波在三维空间中的传播过程, 模拟的重力波振幅随着传播高度的增加以指数形式增长, 对能量传播路径以及有效扰动位能和扰动动能的模拟结果均与理论预测的结果吻合很好. 另外, 随着进程个数的增加, 完成相同的计算量所用的计算时间也显著减少. 这表明本文建立的并行数值模式不但能够很好的模拟重力波的传播过程, 而且能够有效地节约计算时间.      

近壁区理论耗散率模型的改进

提出了在边界层的近壁区,采用共振三波的理论模型描述湍流边界层相干结构,由此解决了对称三维波理论模型不能处理三阶相关量的问题,根据新的理论模型对ε方程逐项进行计算和分析.与原理论模型计算结果相比,本文理论结果与直接数值模拟(DNS)符合更好.     

基于相干结构的大涡模拟壁面模型

首先对湍流边界层近壁区相干结构群进行研究,根据共振三波理论使用沿流向周期性发展的多个三维波构造其初始结构,通过直接数值模拟得到该相干结构群的演化特性.在此基础上,提出了基于近壁区相干结构流动参数的新型大涡模拟壁面模型.确定了此壁面模型与主流大涡模拟衔接的原理方案.计算结果表明:该壁面模型在减少计算量的同时,能够为主流区大涡模拟提供重要的壁面大尺度结构信息,使计算取得准确的结果,实现了壁面模型的基本功能.      

非等间距显式算法在重力波非线性传播数值模拟中的初步应用

采用3阶TVD型显式时间积分方法,建立了2维可压缩大气中重力波非线性传播的非等间距显式数值模式.对小振幅重力波传播过程的模拟结果表明,该模式能够很好地再现小振幅重力波的传播过程,并能保持能量守恒关系和各扰动量之间的相位关系,与线性重力波理论预测的结果吻合很好.对有限振幅重力波的模拟结果表明,与格点数相同的等间距网格中的模拟结果相比,采用非等间距网格能够以更高的分辨率模拟重力波的不稳定传播直到破碎的整个过程;在非等间距网格中的模拟结果与在加密一倍的等间距网格中的模拟结果相同,说明了用非等间距网格模拟重力波的饱和与破碎过程是正确的,并且采用非等间距网格能够大大减少计算量.      

Autonomous navigation for lunar final approach based on gravity gradient measurements

Lunar final approach navigation is critical for pin-point lunar landing in future missions. This study investigates the use of lunar gravity gradient measurements for autonomous navigation of a lunar probe during the final approach phase. As the spacecraft approaches the Moon, the strength of gravity gradient signals improves. A spaceborne gravity gradiometer can precisely measure local gravity gradients, and the latest lunar gravity model GL1500E is used to provide reference values. The employed truncation degree and order of the gravity model are increased stepwise considering the decreasing altitude of the spacecraft in order to reach a compromise between computational costs and model accuracy. An iterative Kalman filter is developed for coupled orbit and attitude estimation using gravity gradient measurements and attitude quaternions obtained from star sensors. A simulated spacecraft with a gradiometer noise level of 0.01 E is considered. Simulation results show that the spacecraft’s position converges rapidly and achieves an accuracy of less than 100 m at the last epoch.     

一种空间相变换热器热设计与仿真分析及其改进

基于相变计算方法即焓法处理相变材料凝固/融化模型,对一个以水为主动冷却介质,内填充石蜡类相变蓄热材料的板式相变换热器的换热进行数值模拟.得到了不同重力条件下冷却面的温度分布,相变材料在融化过程中的动态温度场分布、相变界面分布、融化时间等结果,验证了该相变换热器的可行性.对比该相变换热器在重力与微重力不同条件下的性能差异,利用添加强化传热肋片与泡沫复合相变材料方法,提高了微重力条件下该类相变换热器的效率,可为空间相变蓄热装置的设计及实验研究提供重要参考.      

不同重力条件下单气泡池沸腾现象的数值研究

对不同重力条件下常压饱和水中单气泡池沸腾现象的气泡生长过程及传热特性进行了数值模拟. 采用简化的润滑流模型计算生长气泡底部微液膜的贡献, 而其他宏观区域的气液两相介质则用连续界面模型统一处理. 气液界面形状和加热面上接触线的运动分别采用Level Set方法和固定的表观接触角来近似刻画. 计算结果表明, 气泡生长过程中, 当量直径近似与生长时间的1/3～1/2次方成正比, 重力对相关趋势的影响不大, 但强烈影响着气泡脱落直径和生长时间, 其中脱落直径反比于重力的1/3次方, 生长时间反比于重力的4/5次方. 在固定的核化点数密度条件下, 加热面平均热流密度近似与壁面过热度的3/2次方成正比, 该趋势并不随重力的减弱而改变.      

Evaluation of using two-phase frictional pressure drop correlations for normal gravity to microgravity and reduced gravity

The calculation of two-phase frictional pressure drop (TPFPD) is required by two-phase systems operating under microgravity and reduced gravity. There are a large number of correlations for the TPFPD in tubes under normal gravity. However, it is hard to find out a TPFPD correlation obtained from microgravity and/or reduced gravity conditions, and thus people have to use TPFPD correlations for normal gravity to calculate TPFPD under microgravity and reduced gravity. It is necessary to evaluate the feasibility of such practice. This paper offers a comprehensive review of the TPFPD correlations for normal gravity and an up-to-data survey of the TPFPD experimental study under microgravity and reduced gravity. There are 23 TPFPD correlations for normal gravity reviewed and 135 experimental data under microgravity obtained from the literature. These experimental data are used to evaluate the reviewed TPFPD correlations. It is found that the smallest mean absolute relative deviation (MARD) of the correlations is greater than 34%. Using TPFPD correlations for normal gravity to reduced gravity and microgravity may be acceptable for the first approximation, but correlations intended for microgravity and reduced gravity are needed and more experiments are desired to obtain more data with high accuracy.     

模拟重力波非线性传播的并行数值模式

基于MPI消息传递接u和区域分解的思想,对跳点网格上的串行数值模式进行并行化处理,建立了模拟重力波非线性传播的并行数值模式.针对跳点网格的特点,详细论述了跳点网格系统下区域分解和各子区域间的数据通信.对小振幅重力波传播过程的模拟结果表明,并行数值模式可以很好地模拟小振幅重力波的传播过程,也能保持重力波传播过程中的能量守恒关系.并行数值模式的并行效率可以达到0.65,在理想情况下可以达到最大值1.0.作为与串行数值模式的比较,采用不同网格分辨率以保证计算时间相同,分别用串行和并行数值模式模拟有限振幅蕈力波的非线性传播过程,结果表明,通过引入更多的进程参与计算,并行数值模式可以有效地提高模式的分辨率,模拟出重力波在发生翻转之后,而存演化成湍流之前的这段时间内有Kelvin-Helmholtz billows出现,但是,在相同的计算时间内,串行模式的分辨率较低,不能刻画这些精细的现象.     

Directional solidification of complex-forming eutectic melt of the PbCl2-AgCl dielectric system under conditions of zero, normal and increased gravity

First results in the area of studying the influence of gravity field on the degree of ordering of directionally solified eutectict are presented. As the model complex-forming melts, PbCl₂ and its eutectic with AgCl was chosen. Directional solidification of the mentioned eutectic was realised under conditions of zero, normal and increased gravity. It was shown that increasing level of gravity results in a marked improvement of the degree of spacial ordering of the lamelar structure.     

微重力下胶原蛋白纤维化及羟基磷灰石结晶的初步研究

在长期空间飞行过程中, 骨质丢失是一个严重问题. 羟基磷灰石(HAP)晶体是骨骼的主要成分, 骨骼中的胶原蛋白纤维在HAP生长结晶过程中起到关键作用. 研究了胶原蛋白纤维化过程在模拟微重力和常重力条件下的变化, 对以胶原 蛋白纤维作为模板生长出的HAP晶体形貌进行了观察. 结果表明, 不同浓度胶原蛋白溶液中形成的胶原蛋白纤维, 其内部孔隙数量和尺寸在模拟微重力条件下要明显大于常重力条件下, 胶原蛋白纤维内部孔隙的分布也不同于常重力条 件下的结果. 以模拟微重力条件下形成的胶原蛋白纤维为模板生长出的HAP 晶体主要为立方体状, 而以常重力条件下形成的胶原蛋白纤维为模板生长出的 HAP晶体形貌主要为板状. 该结果有助于未来进一步阐明空间骨质丢失的机理.      

Possible actions of gravity on the cellular machinery.

Since the first flight of the ESA Biorack on the German Spacelab Mission D1 in 1985 evidence has been obtained that biological cells and small unicellular organisms function differently under conditions of microgravity. However, there is still lack of scientific proof that these effects are caused by a direct influence on the cells in the weightlessness condition. The question how normal gravity may play a role in cellular activity is being addressed and the results show that gravity may provide important signals during certain state transitions in the cell. These would be gravity-sensitive windows in the biological process. Also, by amplification mechanisms inside the cell, the cell may assume a state that is typical for normal gravity conditions and would change in microgravity. Experimental tools are discussed that would provide the conditions to obtain evidence for direct action of gravity and for the possible existence of gravity-sensitive windows.