跨音速风扇全环叶片颤振特性的流固耦合分析

发展了求解叶片颤振问题的流固耦合计算方法和全环叶片振动的气动弹性模型,在每一时间步同步求解流体运动方程和叶片振动方程并交换边界信息;流体域求解了非定常雷诺平均N-S方程,得到每一步由于叶片变形而引起的流场变化;叶片变形则由积分叶片表面受到的气动力并求解结构动力学方程得到.颤振分析是在全环叶片模型上进行的,并解除了预先设定叶片间相位角的限制.此方法的显著特征是在一次气动弹性计算过程中,可同时分析叶片多个固有模态、多个节径下的气动弹性稳定性,大大提高了使用时域法进行叶片排气弹分析的计算效率.考察了NASA rotor 67风扇全环模型在堵塞点、最高效率点和近喘点3个气动工况下,节径变化对叶片气动弹性稳定性的影响,给出了不同模态下气弹最不稳定状态对应的叶片振动节径形式.结果表明,振动形式对于叶片气动弹性稳定性的影响很大.     

跨音速风扇全环叶片颤振特性的流固耦合分析

发展了求解叶片颤振问题的流固耦合计算方法和全环叶片振动的气动弹性模型,在每一时间步同步求解流体运动方程和叶片振动方程并交换边界信息;流体域求解了非定常雷诺平均N-S方程,得到每一步由于叶片变形而引起的流场变化;叶片变形则由积分叶片表面受到的气动力并求解结构动力学方程得到.颤振分析是在全环叶片模型上进行的,并解除了预先设定叶片间相位角的限制.此方法的显著特征是在一次气动弹性计算过程中,可同时分析叶片多个固有模态、多个节径下的气动弹性稳定性,大大提高了使用时域法进行叶片排气弹分析的计算效率.考察了NASA rotor 67风扇全环模型在堵塞点、最高效率点和近喘点3个气动工况下,节径变化对叶片气动弹性稳定性的影响,给出了不同模态下气弹最不稳定状态对应的叶片振动节径形式.结果表明,振动形式对于叶片气动弹性稳定性的影响很大.     

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

椭圆轨道编队构型的初始化控制研究

研究了椭圆参考轨道下的编队构型初始化问题,给出轨道面内绕飞构型的三种初始化策略.基于状态转移矩阵研究了构型初始化的两脉冲机动方法;基于高斯摄动方程,以降低燃料消耗为原则,研究了构型初始化的三迹向脉冲和四迹向脉冲方法,推导了速度增量的求解方程;对上述三种构型初始化策略进行了仿真分析,利用序列二次规划和遗传算法对计算结果进行优化.数值仿真表明,三种初始化策略均能完成构型初始化,三迹向脉冲初始化比两脉冲初始化节省约41%的速度增量,四迹向脉冲初始化需要的速度增量与三迹向脉冲初始化基本相等.      

降落伞开伞过程的多结点模型仿真

根据降落伞的结构和力学特征,在轴对称假设下创建了伞衣及回收物系统的多结点结构模型.通过考虑应力,重力和气动力的作用效果,建立了用于无气流攻角平面圆形降落伞充气模拟的多结点结构模型动力学方程组.对开伞过程中的流场变化引入准定常假设,利用simple算法数值模拟求解RNG(Renormalization Group)k-ε湍流模型下的雷诺平均N-S(Navier-Stokes)方程以获得选定时刻的伞衣表面压力分布.结合多结点模型动力学方程组的解算代码和计算流体力学程序,采用流固耦合的方法对选定的平面圆形降落伞模型的开伞过程进行了动态仿真,得到了开伞过程中降落伞外形和特性的变化.通过结果分析和比较,证明了多结点模型的可行性,发展出了一种用于降落伞流固耦合计算的新方法.      

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

非定常N－S方程的线化及其差分求解

在微幅简谐振荡的假设下，对计算平面内非定常Ｎ－Ｓ方程进行了线比，并在ＭａｃＣｏｒｍａｒｋ于１９８０年提出的半隐格式的基础上，发展了一种适用于求解计算平面内线化Ｎ－Ｓ方程的差分格式．用该差分格式求解计算平面内的线化Ｎ－Ｓ方程时，可以避免求解五对角方程；且时间步长可以取得较大，计算效率高．     

非定常N—S方程的线化及其差分求解

在微幅简谐振荡的假设下，对计算平面内非定常Ｎ－Ｓ方程进行了线化，并在ＭａｃＣｏｒｍａｒｋ于１９８０年提出的半隐格式的基础上，发展了一种适用于求解计算平面内线化Ｎ－Ｓ方程的差分格式，用该差分格式求解计算平面内的线化Ｎ－Ｓ方程时，可以避免求解五对角方程，且时间步长可以取得较大，计算效率高。     

高可靠性航空电子设备热分析中的有限体积法

阐述了热对系统可靠性的影响,对比了数值法求解温度场分布的优缺点.针对航空电子设备热分析中同时要求满足高的分析精度和快的计算速度的问题,提出用有限体积法求解温度场分布的计算方法.根据典型航空电子设备的结构特点,建立了温度场的数学模型,给出了边界条件.在此基础上,用有限体积法对待求解方程离散化,给出了离散化过程中网格划分及边界的处理方法.结合实例进行了温度场求解,计算结果与热测量结果进行了对比,得到了较高的分析精度,验证了方法的有效性.最后对计算误差进行了分析,并提出了进一步热设计的合理化建议.      

矩形齿磁流体密封装置的磁场计算

 考虑了不同磁介质磁化性能的非线性,并采用介质边界上磁场传播关系对介质边界漏磁的精确描述,利用不等间距网格有限差分数值方法,对矢量磁位描述的轴对称场偏微分方程进行求解,计算出磁流体密封装置中具有非线性、稳态特征的强磁场分布,得出了不同齿数密封装置的磁场分布,克服了磁路法计算中对磁介质磁化性能和漏磁理想化估计带来的误差,计算结果准确、合理、可靠,为密封装置的结构设计提供了设计依据.     

磁场重联的二维粒子模拟研究

利用二维全粒子模拟方法研究了无碰撞等离子体中的磁场重联过程，得到了不同区域的离子和电子速度分布．计算结果表明，电子和离子在扩散区中的不同动力学特性产生的Hall电流使磁场的y分量By呈现四极形分布．离子和电子的速度分布偏离了初态时的Maxwell分布，呈现非局域的多重分布．同时由于磁场重联而产生的电场使电子在X点附近得到加速和加热，因而在电子的能谱分布中形成-高能尾。     

行星际结构与垂直无碰撞激波的相互作用

应用一维混合模拟方法数值研究了两种行星际结构──反向磁场和高密度等离子团与垂直无碰撞激波的相互作用．结果表明，随着激波上游区磁场的反向，下游区磁场将逐渐改变符号，且等离子体密度和速度分别呈现较强的湍动．激波上游和下游的物理量依然满足Rankine-Hugonoit关系．当高密度的等离子体团通过垂直无碰撞激波时，部分质子被激波反射，部分质子被加速并进入下游区域．由于质子速度分布为非Maxwell分布，在激波下游也激发出较强的湍动．     

磁层顶通量传输事件的经验重构

在地球磁层顶附近观测到的通量传输事件（Flux Transfer Event,FTE）一般被认为是瞬态局域磁重联的产物,是太阳风质量、动量和能量进入地球内磁层的重要通道.重构FTE的磁场结构可促进对其形成、演化过程及其与周围等离子体环境相互作用的理解.Grad-Shafranov重构法和磁通量绳拟合法等传统磁场重构方法适用于满足特定物理条件的磁场结构.基于平面线性插值原理,设计了一种不受具体物理条件限定的二维FTE磁场结构重构法.模型测试以及对THEMIS和Cluster卫星簇分别观测到的两个FTE的实际应用表明,在合适的多卫星位形条件下,该方法能快速有效重构出FTE的磁场空间分布,有助于推测FTE的磁场线位形,理解卫星测量数据的时间变化,以及分析等离子体物理量相对于FTE的磁场空间分布特征.      

火星磁场产生的原因及其分布

构建了一个可以得到火星赤道面上磁场分布的模型. 模型根据卫星观测数据, 提出了火星电离层、磁层顶和磁尾电流片上都各自通有电流的假设. 由电流的连续性条件可知, 这三种背景条件下的电流之间满足一定关系, 即火星磁层顶上的总电流是电离层上的总电流与磁尾电流片上的总电流之和. 这些电流产生的磁场与太阳风磁场共同构成了火星赤道面上磁场分布. 通过计算发现, 采用这种磁场模型得到的结果与目前卫星所观测的结果以及与采用其他方法得到的结果符合得较好.      

基于连续性条件的体PIV后处理技术

对粒子图像测速技术(PIV,Particle Image Velocimetry)体空间内三维速度场的后处理技术进行研究,引入不可压缩连续性方程的约束对实验测得流场进行修正,使得修正后流场严格满足差分形式的连续性方程.通过流场修正的后处理技术实现拉格朗日乘数法求解全流场修正速度场二范数的极小值.理论推导严格证明了修正后的速度场在数学上更接近满足连续性方程的真实流场.对解析流场引入高斯噪音来进行人工模拟的测试,其结果进一步验证了修正后流场能更趋近于真实流场,对流场中出现的PIV常见的高斯分布噪音有一定的消弱作用.该技术被成功应用于实验数据的后处理,结果表明其能改善测量速度场的平滑度.     

Velocity correlation functions of charged particles derived from the Fokker–Planck equation

An important ingredient in theories for diffusion of charged particles across a mean magnetic field are velocity correlation functions along and across that field. In the current article we present an analytical study of these functions by investigating the two-dimensional Fokker–Planck equation. We show that for an isotropic pitch-angle Fokker–Planck coefficient, the parallel velocity correlation function is an exponential function in agreement with the standard model used previously. For other forms of the pitch-angle diffusion coefficient, however, we find non-exponential forms. Also a new, velocity correlation function based, approach for deriving the so-called Earl-relation is presented. This new derivation is more systematic and simpler than previous derivations. We also discuss higher-order velocity correlations and the applicability of the quasi-normal hypothesis in particle diffusion theory. Furthermore, we compute velocity correlation functions across the mean field and develop an alternative theory for perpendicular diffusion.     

流动对太阳黑子精细结构的影响

考虑磁扩散率为温度的函数以及绝热条件,线性扰动磁流力学方程组的色散方程是频率的三次方程。由于沿磁场流动的速度u₀小于声速u_s,即ε=u₀2/u_s2是一个小量,在零级近似下解色散方程,得到了磁扩散不稳定性。本文进一步讨论在一级近似下小流动对磁扩散不稳定性的影响。结果表明,小流动不破坏磁扩散不稳定性的存在,并且,它对波数空间中的不稳定区也影响不大。      

基于经验分布的区间数据分析方法

现有区间数据分析的方法通常假设数据在某一区间上服从均匀分布,这在实际数据分析中通常是不成立的.针对此问题,在原始数据来源于连续分布的简单假设下,利用经过分布函数变换后的随机变量服从(0,1)上的均匀分布,分别采用经验分布函数和核估计对原始数据的分布函数进行估计.基于此设计变换,对变换后的数据进行均匀分布的假设检验,通过检验后进行后续的区间数据分析,使得均匀分布的假定得以成立,保证了统计理论上的严谨性.数据模拟结果表明,将经验分布函数变换后的数据作为研究对象,进行区间数据分析,所得到的统计建模结果更加合理且具有较强的解释力.