太阳表面水平运动驱动的磁圈与开放磁场重联的模拟研究

太阳大气的诸多观测事件(如耀斑、喷流等)均被归因于磁重联产生的能量转换. 近年来, 关于太阳风起源, 有研究提出了磁重联使闭合磁圈开放为太阳风供应物质的新模式. 在该模式中, 闭合磁圈被光球超米粒组织对流携带, 向超米粒边界运动, 与位于边界的开放磁场相碰撞进而发生磁重联. 该模式中磁重联的驱动及其效应是本文的研究目标. 磁流体力学(MHD)数值模拟是研究太阳大气磁重联物理过程的重要途径. 本文建立了一个二维MHD数值模型, 结合太阳大气温度和密度的分层分布, 在超米粒组织尺度上模拟了水平流动驱动的闭合磁圈与开放磁场的重联过程. 通过对模拟结果的定量分析, 认为磁重联确实能够将闭合磁圈的物质释放, 进而供应给新的开放磁结构并产生向上流动. 该结果为进一步模拟研究太阳风初始外流奠定了基础.      

2.5维自适应磁场重联MHD模式

磁雷诺数(R_m)是影响磁场重联的重要因素. 真实的物理环境中R_m往往很高, 例如, 在行星际空间和太阳日冕中R_m通常大于104量级. 高R_m条件下的磁重联表现出很多异常特性, 然而高R_m条件下的磁场重联数值模拟需要很高的时空分辨率, 否则很难分辨出重联过程中形成的薄电流片. 本文基于自适应软件包PARAMESH将并行自适应网格技术引入磁场重联数值模拟, 建立了一个2.5维自适应磁场重联MHD模式, 研究高磁雷诺数条件下重联的动态演化过程, 进而将不同磁雷诺数的参数进行对比研究. 结果表明, 该模式可以自动捕捉到磁场重联产生的奇性电流片, 高磁雷诺数条件下产生的慢激波结构可提供一种快速磁能释放机制.      

基于数值分析的空速-动压信息处理方法

从数值计算方法的角度对大气数据校验系统设计中空速、指示空速与动压、静压转换关系进行研究,提出了一种新的转换方法。在PC计算机和ARM嵌入式计算机S3C2410上进行了模拟仿真,给出了计算程序和分析结果,与传统方法相比,具有转换速度快、程序简捷、精度高的特点,该方法用软件容易实现,可以很方便地应用到其它传感器的研究与设计中。     

空间平台发射拦截器动力学与最优脉冲拦截

研究了空间平台发射拦截器的两体耦合动力学,以及拦截器拦截目标星的最优脉冲控制问题。平台首先与目标星形成绕飞关系,保持其发射筒轴线始终瞄准目标星。接到发射指令后,拦截器从发射筒中射出,本文采用拉格朗日第二类方程建立了发射过程平台-拦截器两体动力学模型。因为两体耦合影响,平台姿态偏转,拦截器出筒时的速度已经不能瞄准目标星。通过小型火箭发动机给其施加速度脉冲,使其进入拦截轨道,保证拦截的同时,将脉冲速度最小化以节省燃料,本文将其归结为一个非线性规划问题,采用三级优化的策略来求解。在拦截飞行时间相较于平台绕飞目标星的周期是小量的条件下,可以视绕飞平均角速度为小参数,采用正则摄动方法求出非线性规划的一阶近似解,然后以此为迭代初值,寻找最优真解。最后进行了数值仿真验证。     

移15°双十二相隐极同步整流发电机数学模型的研究

分别在abc坐标系和dq0坐标系下建立了移15°双十二相隐极同步整流发电机的基本数学方程,推导了电机电压、磁链、输出功率、电磁转矩和转子运动方程。在MATLAB/Simulink中建立了移15°双十二相隐极同步整流发电机的仿真模型。对两路输出进行了不平衡负载仿真和试验。仿真与试验结果相比基本一致,验证了移15°双十二相隐极同步整流发电机数学模型的正确性。这为双十二相同步整流发电机运行性能和励磁控制系统的研究奠定了基础。     

基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法,建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题,旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法,通过添加冗余约束,避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块,对基于运动嵌套网格 的空间流场进行求解,湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块,基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解,验证了网格变形模块和CFD模块的有效性,然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上,计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷,并与飞行测试数据进行了对比。计算结果表明,文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析,提高了旋翼气弹载荷的预测精度。     

燃烧室污染物生成大涡模拟模型的研究

采用亚网格动力学模型和联合概率密度函数(PDF)亚网格模型对污染物NO的生成进行大涡模拟.计算中采用的数学模型有:k方程亚网格尺度模型来估算亚网格湍流粘性;亚网格EBU燃烧模型估算化学反应速率;在交错网格体系下, 采用SIMPLE算法和混合差分格式来求解差分方程.数值分析两种不同亚网格污染物生成模型和两种进口气流温度对NO排放的影响, 计算结果与实验数据比较表明:两种亚网格污染物生成模型都能较好的预估燃烧室污染物NO的生成, 但亚网格动力学模型的计算工作量要比联合概率密度函数亚网格模型少得多, 更适合工程应用.      

新型航空润滑油油膜拖动力计算研究

在自制的试验装置上测量了新型航空润滑油油膜拖动力;给出了高剪切率弹流润滑下的Ery-ing特征应力、极限剪切应力、剪切弹性模量的新型航空油流变参数拟合计算式;确定了该油本构方程;建立了该油油膜拖动力计算式, 并进行了试验验证.试验结果表明:该油在试验条件下表现为粘弹性;拖动力的理论计算值与试验测量值相当吻合.给出的拖动力计算方法可用于指导润滑分析领域的工程计算.      

一类高阶有理差分动力系统

讨论了高阶有理差分方程xn 1=(αxn-k 1)/(-1 xnxn-1…xn-k 1),for n=0,1,…,其中α<0,k∈{1,2,…},x-k 1,…,x0都是实数,并且使得x-k 1…x0≠1。给出了这类方程解的明确表达式,并对其解的稳定性和收敛性也做了进一步讨论。     

内埋式弹舱舱门气动载荷计算分析研究

采用中心有限体积法进行空间离散和四步Runge-Kutta显式推进方法进行时间推进数值求解了非定常欧拉方程;采用Delaunay方法生成内埋式弹舱舱盖打开时的非结构网格.数值模拟了内埋式轰炸机弹舱单开舱门在分别打开8°和45°时的复杂流场.通过对上述复杂流场的分析,发现作用于舱门的气动载荷具有使其趋向于关闭的特性,且有随打开角度的增大而减小的趋势.