能量粒子进入磁层的数值模拟研究

来自行星际的能量粒子进入磁层的过程一直是磁层研究的热点,也是空间天气研究关注的焦点之一,而以往的工作大都集中在讨论不同行星际条件下粒子进入磁层的统计规律.本文在T89磁层模式的基础上,建立了太阳能量粒子进入磁层运动的程序,计算了几种不同能量的质子从不同方位角入射到磁层的运动轨迹.模拟结果表明,能量粒子沿着磁力线方向才可能穿越磁层到达地球表面,越偏离这个方向,则越早被反弹.能量越高的粒子进入到地球磁层上空的角度范围越大,但仍然只有沿着磁力线入射的粒子才能到达地球表面.这些结果与理论预言是一致的.      

耀斑激波传播的数值研究

采用一维磁流体力学模型和激波装配法,分析耀斑激波在行星际空间的传播特性和激波下游区的波动结构,并就激波装配法和激波捕捉法的模拟结果的精度和可靠性进行比较.     

带肋变截面回转通道内流动与换热的数值模拟

开发了三维流动换热的通用计算程序,数值研究了带肋变截面回转通道内流动与换热的特性.湍流模型采用低雷诺数k-ε模型.通道肋间距为25mm,肋高分别为1mm,1.5mm,2mm,冷气进口雷诺数Re分别为7500,12500,18500,25000.计算结果表明:①通道的平均努赛尔数均随进口雷诺数的增大而增大;②对于Re=7500和12500,肋高越高,换热越强;对于Re=18500和25000,肋高为1.5mm的通道换热最强;③局部雷诺数的不同和离心力的影响导致通道内各区域的局部换热随肋高的变化趋势并不一致;在进口段,肋高越高,换热越强;在出口段,当Re=7500和12500时,肋高越高,换热越强,而当Re=18500和25000时,存在最佳肋高1.5mm.     

板料成形数值模拟接触搜索模型的研究

基于参数曲面描述法和有限元网格描述法,提出模具几何形状模型的曲面单元描述法,采用了逐级更新的接触搜索算法,推导出求交迭代初值优化选择方案及满足无穿透原理的接触判断准则,并建立了板料成形数值模拟接触搜索模型,此模型不仅改善了数值模拟的收敛性,而且提高了计算效率和稳定性.     

重力波波包相互作用时空演变的数值研究

本文设计了一种研究二维空间中重力波波包共振相互作用演变的数值模式.利用此模式获得二个大振幅重力波波包通过碰撞而发生的完整的相互作用过程.数值结果表明,能量上行重力波可以通过波-波相互作用激发能量下行重力波;在相互作用的时空演变过程中,3个波包的能量收支情况和相互作用的强度具有局地性;由于考虑了空间传播效应,相互作用过程不再具有周期性,而是出现一个具有不同物理含义的相互作用特征时间,这个新的相互作用特征时间和相互作用最终的发展程度不仅由相互作用系数和初始波包振幅决定,还受到波包空间尺度和波包群速度的相对大小控制.     

太阳风--磁层相互作用的磁流体力学数值模拟研究

磁层位于地球空间的最外层, 太阳风与磁层的相互作用是空间天气变化因果链中承上启下的关键环节, 是揭示地球空间天气基本规律的关键科学课题. 地球空间由于时变、多成分、多自由度的关联相互作用,使得传统的理论分析变得非常困难. 数值模拟作为近几十年发展起来的一个新的研究手段,对地球空间的理论和应用研究产生了深刻的影响. 国际上磁层的全球MHD数值模拟工作开始于20世纪70年代末, 最初的研究局限于二维空间. 由于磁层内在的三维特性, 20世纪80年代三维MHD数值模拟工作兴起. 本文简要说明了三维全球磁层MHD (磁流体力学)研究的特点及现状, 给出了三维全球磁层模型的基本框架, 综述了行星际激波与磁层相互作用、大尺度电流体系、重联电压和越极电位、磁层顶K-H不稳定性等方面的太阳风--磁层相互作用的MHD数值模拟的研究进展.     

发动机内剩余推进剂排放对卫星表面污染的研究

应用真空物理的相关理论,分析了在高空超高真空环境中,发动机二次启动前推进剂排放及其在空中迁移的物理过程,给出模型假设和模拟方程,并针对此真空环境及物理过程应用合适的计算方法对方程进行求解.最后计算了推进剂的排放对卫星造成的污染,给出计算域内浓度变化的非定常过程及卫星表面不同时刻的最大污染量,为推进剂排放方案的实施与否提供理论依据.     

1998年5月2日日冕亮度观测图的数值研究

采用数值手段模拟了1998年5月2日日冕亮度观测图.计算模式改进为球坐标系下特殊的二维理想磁流体(MHD)模型,即把(γ,φ)坐标建立在SOHO观测日冕亮度的子午面上,消除了子午面极区的几何奇异性.根据SOHO日冕观测布置磁极子得到初始磁场位形,内边界条件采用自洽的投影特征线边界条件,计算迭代出稳态的多冕流磁场结构,得到了与观测基本一致的亮度图.计算结果表明在太阳表面附近磁场位形对太阳风等离子参数分布起控制作用.     

基于一维模型的高频电波加热电离层F层数值研究

基于电离层一维仿真加热模型,详细介绍了模型中电子的动量方程、连续性方程、能量方程和各类参量表达式,利用对角矩阵追赶法数值求解电离层F层加热过程,分析了不同时次电子数密度和电子温度的变化,讨论了不同频率和不同功率电波加热的情形.结果表明：当高频电波加热高电离层时,电子温度迅速增加,并很快趋近于饱和状态;电子密度的变化较为迟缓,但在加热过程中其变化幅度却越来越大;电子密度变化量在沿磁场方向上形成空洞和上下稠团两峰一谷构型;频率越大、功率越高的电波加热时,电子密度的变化也越大,但存在一适值频率的电波对电子温度的影响最小.     

磁尾等离子片中偶极化锋面的数值模拟研究

利用守恒型TVD格式对8波模型磁流体方程组进行数值模拟, 对磁尾中偶极化锋面的物理和演化特性进行研究. 构建了由BBF类型通量管机制产生的偶极化锋面数值模拟模型, 该模型由磁尾平衡模型、亚暴增长相模型和亚暴触发及BBF形成模型三部分组成. 数值模拟结果很好地再现了磁尾中BBF类型通量管机制产生的偶极化锋面特性. 伴随着高速 流的出现, 磁场B_z分量呈非对称双极变化结构, 即锋面前减小为负值, 在锋面上急剧增大. 当B_z增大到极大值后回落并趋于稳定. 随着偶极化锋面伴随地向高速流向地球运动, 偶极化锋面上B_z的变化越来越小.      

对旋轴流风机非定常流场数值模拟

利用全隐格式的非稳态压力修正求解方法,分析和研究了对旋叶轮级间流场非定常流动情况,上游尾迹与下游叶片位势作用对流场的影响,并与动态实验结果进行了有益的对比.结果表明:非定常流场的计算更能反映叶轮机实际运行性能,如果能很好地组织对旋叶轮的非定常流动,还可以进一步提高其压比和效率.     

Euler方程数值模拟绕悬停旋翼桨叶的流动

本研究应用格心格式有限体积龙格－库塔时间推进法求解三维欧拉方程,模拟绕旋翼桨叶的流动,民流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角,选择Ｏ－Ｈ型式生成绕桨叶贴体的弦向和展向网格,通过算为进一步寻求有效的Ｅｕｌｅｒ或ＮＳ方程数值解积累了经验。     

对旋轴流风机的稳态流场数值模拟

在任意曲线转动坐标系及重整化群(RNG)方法改进的k-ε两方程湍流模型基础上,提出了"通用进出口边界条件"的给法,并利用改进的强隐式(MSIP)方法进行计算迭代,提高了对多排叶轮机数值模拟的稳定性和准确性.通过对某对旋轴流风机定常全三维流场的计算,并与设计参数和实验性能参数进行对比分析后,验证了在风扇/压气机的设计中,不宜在叶尖区分配最大负荷,这种设计对对旋轴流压缩系统尤其不利.     

定常压力梯度边界层相干结构的直接数值模拟

根据流动稳定性理论,将不稳定波的一个周期作为相干结构的初值,采用直接数值模拟方法对有压力梯度湍流边界层中相干结构的演化进行了研究,得出其各种特性的变化与实验观测到的结果一致.     

压气机间隙流与处理机匣作用的三维数值分析

利用Numeca CFD对某一压气机静子叶栅的间隙流动进行流场计算,并将其与具有圆弧斜槽处理结构的间隙及流场计算结果进行对比、分析.详细揭示了叶栅顶部间隙区及处理槽内的流动特征.结果表明,通过采用机匣处理,改变了压气机气流流路,形成叶尖漏流的通道,减少漏流下洗对叶片通道造成的阻塞.叶尖附近主流在叶盆尾缘气流高压作用下进入斜槽,而后气流在叶背前缘以高速由斜槽射入主流,该高速射流有效地扫除叶尖易失速的附面层,从而延迟气流分离,扩大压气机的失速裕度并减少二次流损失.这种籍助于动量交换而形成高速射流对主流的作用可能是机匣开槽结构改善失速裕度的主要原因.      

水洞收缩段流场的数值模拟及优选

介绍了一种水洞非对称收缩段的数值模拟方法和程序,计算采用有限体积预估修正格式求解非定常欧拉方程.结合具体水洞设计对几种典型收缩曲线的三维流场进行了计算和分析,给出了沿轴向的压力分布,出口截面速度分布及分布的均匀性,选出最佳曲线,为水洞的设计提供依据.结果表明:本文的计算方法和程序对于水洞或风洞的三元收缩段,是有价值的分析和辅助设计手段.     

一种板式推进剂管理装置（PMD）性能的数值仿真

板式推进剂管理装置（PMD）是板式表面张力贮箱最重要的部件之一．为分析板式PMD管理推进剂能力,以某板式贮箱为研究对象,采用三维非定常流体（VOF）体积函数的两相流模型,对不同填充比时板式PMD的推进剂管理性能进行数值仿真,得到微重力下推进剂在PMD上及贮箱内部的分布情况,仿真结果与国外空间搭载试验结果比较吻合,数值仿真结果为板式贮箱的性能提供验证依据     

板料冲压成形回弹的数值模拟

回弹是板料成形过程中不可避免的现象.回弹的数值模拟作为解决回弹问题的有效方法,越来越受到工业界和学术界的重视.完整的回弹数值模拟分为2步,即冲头加载成形过程模拟和冲头卸载回弹过程模拟.本文比较了静力隐式和动力显式2种算法在计算回弹问题时的优劣,采用静力隐式算法建立了计算冲头加载成形过程模拟和冲头卸载回弹过程模拟的有限元模型,并利用自行开发的软件给出了计算实例.      

微重力板式贮箱推进剂重定位的数值仿真

板式表面张力贮箱内推进剂重定位对确定推进剂分布情况、研究晃动影响、提高控制精度等具有重要意义．为研究板式贮箱内推进剂重定位的规律,对微重力下板式贮箱内液体重定位问题进行数值仿真．计算时使用三维非定常两相流动流体（VOF）模型,对某一板式贮箱寿命末期在不同微重力加速度情况下各种重定位过程进行数值仿真,得到各种工况下重定位的全过程,以及定位后推进剂的分布情况．数值仿真结果为板式贮箱的设计提供有利依据．     

阻尼网特性数值模拟

依靠经验公式和工程估算等传统方法,无法对阻尼网性能进行精确评估.为确定阻尼网压力损失系数和降湍性能,采用计算流体力学方法,结合适当的边界条件,对阻尼网进行了模拟.使用数值模拟能够得到不同开孔率阻尼网在不同雷诺数、不同入射方向、不同目数的损失系数,与试验结果更为接近,在30°~45°大角度入射时得到的损失系数更为精确;在入射气流与阻尼网平面呈一定夹角时,阻尼网后的气流压力和速度呈现脉动趋势,传播距离大约为100d;在雷诺数小于40时流动保持层流状态,扰动传播距离为50d,此时降湍效果最好,随雷诺数增加,扰动传播距离增加至400d;在流动未失稳时,开孔率越低,降湍效果越明显,开孔率低于0.5时流动容易失稳;开孔率保持不变,随目数增加阻尼网损失系数增加明显,降湍能力提升.因此可以根据数值模拟结果选择阻尼网的最优参数.