吸气式高超声速飞行器耦合运动数值模拟

为了研究吸气式高超声速飞行器在俯仰/滚转两自由度耦合运动下的动稳定性问题,基于气动/运动耦合数值模拟方法并结合理论分析,针对一种类似SR-72构型的吸气式高超声速飞行器开展了进气道通流状态下强迫俯仰/自由滚转耦合运动数值模拟研究。结果表明:强迫俯仰/自由滚转耦合运动下,滚转通道的动力学方程可以简化为有阻尼的Mathieu方程形式,并且可以求得稳定性边界。理论分析表明:滚转通道的动稳定性与俯仰振荡角频率相关,在飞行器滚转振荡固有角频率附近存在2个临界角频率,当俯仰振荡角频率位于2个临界角频率之间时,滚转通道是动不稳定的。在俯仰振荡振幅较小时,数值模拟结果与理论符合较好,但实际的临界角频率与理论分析求解的值有一定偏差;数值模拟结果表明随着俯仰振荡振幅增大,导致滚转发散的角频率范围变得更宽,且向更高频率偏移。      

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

关于拟线性椭圆型方程在广义Sobolev空间中解的正则性

本文讨论了下面方程Ｄｉｒｉｃｈｌｅｔ问题在广义Ｓｏｂｏｌｅｖ空间中解中的正则性－ｄ／ｄｘ１ａ１（ｘ，ｕ，Ｄｕ）＋ａ（ｘ，ｕ，Ｄｕ）＝０，ｘ∈Ω其中Ω∈Ｒ是有界区域，证明了上述问题在Ｗ（Ω）和Ｗ（Ω）在存在有界广义解。     

一类非线性椭圆型方程非平凡解和分歧点的存在性

本文给出了所谓集中到列紧原理的另一形式,并利用它讨论了一类非线性椭圆型方程广义银和分歧点的顾在性。     

旋转对管内流体流动的影响

本文描述了旋转管流的基本运动方程。在笛卡尔坐标系和柱坐标系中,针对两种基本的几何流动,推导出旋转情况下管内的流动方程,定义了影响旋转管流的流动和换热的两个无量纲准则数J和Ra,分析了二次流产生的原因,从而说明旋转对流动的影响以及旋转管流的运动特征。     

预移相型恒温热线(膜)流速计的动态响应方程

分析传统恒温热线热膜流速计(简称 HWFA)动态方程后,提出了一种新型HWFA 线路模型——主电桥预移相模型,推导出真实条件下小扰动动态响应方程,并按此原理设计制造了崭新的高性能 IFV—900型 HWFA 样机。该机全部革除了传统HWFA 的三个调节参量,大大地简化了调节步骤;不自激振荡的范围明显增宽,极大地减少了使用上的麻烦;频带宽度比当今世界流行的 HWFA 增宽了2～5倍;加上一系列智能化功能后使自动化程度达到了新的高度。这一新型 HWFA 的问世将可能导至 HWFA技术的重大改进和新的变革。首先阐明主电桥预移相模型的提出和真实条件下小扰动动态响应方程的推导。随后将用若干姐妹篇阐明动态方程频率响应调节元物理值的解析解,频率特性实验分析手段的电扰动试验响应情况的理论分析等内容。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。