毛细泵环路系统运行特性的数值研究

给出了描述毛细泵环路(CPL)系统运行特性的非稳态模型,并运用该模型模拟了系统的启动过程,考察了储液罐中温度波动以及蒸发器上热负荷的突然改变对系统稳定性的影响;从数值模拟的角度更好地预测和分析CPL系统运行的过程和特点,该模型的数值结果与实验现象和实验数据相比有较好的一致性。通过该模型还可以进一步揭示CPL系统运行过程中启动困难、储液罐温度不能波动过大以及热流密度突然改变引起系统失效等内在原因,为设计出性能优良的CPL系统提供了理论依据。     

毛细泵回路蒸发器的数值研究

文章针对毛细泵回路(CPL)蒸发器内的多孔芯结构, 建立了描述换热和流动的非稳态数值模型。该模型采用焓模型来处理多孔介质内工质的相变问题, 通过分区计算既避免了采用耗时复杂的自适应网格技术又能方便地跟踪汽液界面, 为处理多孔介质中的毛细力带来方便; 在处理流动问题时, 通过对源项的处理考虑了毛细抽力的作用, 体现了毛细抽力是CPL的动力源, 从给定的算例来分析, 该数值模型较好地反映了毛细芯的工作特性。     

耀斑激波传播的数值研究

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

限制器对高超声速体表面热流数值模拟的影响

采用Roe的FDS（Flux Difference Splitting）差分格式结合min mod,van Leer和Osher-C限制器,对高超声速体外部流场进行数值模拟.研究了限制器对高超声速体气动力/热数值模拟的区别,对比分析了3种限制器对三维模型表面热流数值模拟结果的影响.研究表明,限制器对数值模拟精度有重要影响,且对气动热数值模拟的影响大于对气动力的影响;对简单或较复杂的气动外形,采用Roe的FDS格式结合min mod限制器获得的表面热特性比它结合van Leer和Osher-C限制器时要好,可为热防护设计提供可靠结果.      

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

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

上游双扰动对圆柱气动特性影响的数值模拟

通过用商用计算流体力学软件FLUENT 6.0中的SIMPLE方法求解二维非定常层流模型,对雷诺数Re=200时上游双扰动体对下游圆柱气动特性影响进行了数值模拟研究.计算结果表明当双扰动体相隔较近(H/d=1.1,1.33)时,其后形成单一的Karman涡街,它们对下游圆柱的干扰作用类似于单扰动体时的情况,但是流动结构发生变化的临界间距Lcri较之有效宽度相同的单扰动体的要大,因此可以推断出在该种布置下双扰动体和圆柱绕流会出现2种流动模式.在空穴流动模式下面,下游圆柱的阻力和升力脉动最多减小94%和79%,系统阻力(包括圆柱和扰动体)也能减小33%.而当双扰动体尾流为双Karman涡街时(双扰动体相隔较远H/d=3.33,5.0),圆柱迎风面部分区域仍受到来流的冲击作用,气动力减小有限.     

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

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

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

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

背景压强对电推进羽流场影响的数值模拟

真空舱内背景压强是电推力器地面试验过程中影响工作性能评估和羽流场参数诊断的重要参数。针对LIPS-200型离子推力器羽流场参数的数值仿真中采用的背景压强建立方法进行了仿真分析。仿真中采用混合粒子网格(PIC)方法和直接模拟蒙特卡罗(DSMC)方法处理羽流场中等离子体运动和粒子间碰撞,分别采用虚拟粒子和计算粒子建立压强的方式,对电推进羽流场进行了数值模拟,并与绝对真空环境进行对比分析。结果表明:背景压强的存在导致中性粒子和电荷交换离子数密度较绝对真空环境高1个量级以上。虚拟粒子可大幅提高计算效率,获得的流场中电荷交换离子分布与计算粒子结果相近,但中性粒子分布相差较大,虚拟粒子无法表征壁面及真空泵的影响。      

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

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

用全球原始方程半谱模式研究QBO对行星波传播的影响

本文建立了一个全球原始方程半谱模式,模拟了赤道上空风场准两年振荡(QBO)及其相应的副热带急流大小对冬季半球行星波向上传播及平流层突然增温的影响。结果表明,波数1的行星波在QBO东风相比西风相更易向上传播,平流层增温更快更强。波数2则相反。QBO对低纬对流层里的行星波上传的影响限制在低纬低平流层,对中高纬影响不大。     

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

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

微吹减阻技术影响因素的数值模拟

以NASA格伦中心的PN3和PN23型微孔壁板为原形建立"1排微孔"、"4排微孔"、"8排微孔"和"16排微孔"4种计算模型,针对不同的几何、物理参数进行了微吹技术MBT(Micro-Blowing Technique)降低平板摩擦阻力的数值模拟参数研究.结果显示:微量吹气使主流边界层近壁面流动减速,从而改变了壁面局部摩擦力的分布;不同几何、物理参数对MBT技术减阻能力的影响满足一定的规律.参数研究的结果可为该技术的减阻应用提供一定的指导.     

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

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

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

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

微重力量值对气液两相流相分布及液相湍流统计量影响的数值研究

气液两相流在空间领域具有广阔的应用前景, 深入理解微重力量值(微重力大小)对相分布和液相湍流的影响十分必要. 采用欧拉elax-elax拉格朗日双向耦合模型深入研究了不同微重力量值对相分布和液相湍流的影响. 液相速度场通过直接数值模拟求解, 气泡的运动轨迹由牛顿运动方程跟踪. 研究表明, 气泡分布和液相湍流与微重力量值均具有直接联系. 在低微重力量值下, 气泡近似均匀分布在槽道内, 且对液相湍流统计量几乎没有影响; 然而当微重力量值较高时, 大量气泡聚集在壁面附近, 液相湍流由于气泡的注入受到极大调制.     

非均质壁面对液滴俘获能力的数值模拟研究

为研究不同非均质壁面对液滴的俘获能力,采用界面追踪法（FTM）结合广义滑移边界建立了接触角模型,对液滴在非均匀润湿的非均质壁面上的运动过程进行了数值模拟研究。液滴在倾斜壁面上受到重力作用由均匀润湿部分下滑至非均匀润湿部分,通过改变Bo数、Oh数、非均匀润湿程度研究了液滴在非均匀润湿区域的运动规律。研究表明:Bo数越大,液滴运动受壁面阻力影响越小,液滴下滑的速度越快,液滴越难以被俘获;Oh数越大,液滴运动受壁面阻力影响越小,液滴越难以被俘获;非均匀润湿程度越大,非均质壁面对液滴的阻力越大,液滴越易被俘获。      

环境压力对胶体推力器喷雾过程的影响

对应用于卫星微推进胶体推力器的喷雾过程进行模拟研究,采用拉格朗日粒子跟踪法计算带电液滴的运动轨迹,获得抽取极板前后液滴速度和密度的空间分布．分析推力器工作环境压力对喷雾扩散角、液滴轴向速度和推力器性能的影响．模拟结果表明：随着推力器工作环境压力的降低,液滴轴向速度增大,喷雾扩散角变化较小,推力器性能提高．当工作环境压力低于一定值时,其对喷雾过程的影响可以忽略,而推力器性能基本保持不变．     

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

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

应用机－电模拟理论研究测力传感器的动态特性

在某些特定情况下，某些电器系统或电路网络的运动微分方程在数学上与力传感器的运动微分方程是相似的，因而说测力传感器的机械模型与电器元件的网络模型也是相似的。从传递函数的角度讨论了这种相似性并指出了它在应变式测力传感器动态特性研究中的应用。