CFD数学模型的线性化方法及其应用

计算流体力学(CFD)方法不仅仅起到数值模拟的作用,它本身是一个复杂的非线性系统。在流动稳定性分析、气动弹性分析、优化设计以及流动控制等领域,从系统的角度出发,对CFD数学模型线性化后,可以对模型的系统矩阵进行定量分析获得更多的系统特性。但是CFD数学模型往往非常复杂且阶数很高,因此其线性化系统矩阵的获得比较困难。鉴于此,采用人工编程和自动微分相结合,构造有限体积法并行CFD模型的线性化系统矩阵。其中自动微分只被用来得到每个界面通量的局部雅可比矩阵,而采用人工编程方法来实现并行环境下的稀疏雅可比矩阵的组装。线性化系统的并行求解采用了块雅可比预处理的广义最小残量法,每个并行进程内部则采用零填充不完全LU分解预处理。为了验证这种线性化方法,上述方法被用于:① NACA 0012翼型的非定常绕流线性系统构造与求解;② NACA 0012翼型稳态流动的伴随方程构造与求解;③ AGARD wing 445.6机翼颤振问题降阶建模。上述三个算例的结果与CFD模拟的吻合一致。     

Nonlinear electron-acoustic rogue waves in electron-beam plasma system with non-thermal hot electrons

The properties of nonlinear electron-acoustic rogue waves have been investigated in an unmagnetized collisionless four-component plasma system consisting of a cold electron fluid, non-thermal hot electrons obeying a non-thermal distribution, an electron beam and stationary ions. It is found that the basic set of fluid equations is reduced to a nonlinear Schrodinger equation. The dependence of rogue wave profiles on the electron beam and energetic population parameter are discussed. The results of the present investigation may be applicable in auroral zone plasma.     

基于多项式拟合SDRE的三维导引律设计

应用多项式拟合的SDRE方法结合改进的极坐标系(MPC)设计了三维次优导引律。介绍了 SDRE方法与多项式拟合的SDRE方法,后者是前者十分优越的逼近;推导了MPC下的弹目相对 运动方程,将球坐标下的六个状态方程减少到了三个并且满足多项式拟合SDRE方法的应用前 提;在此基础上,推导出了三维拟合SDRE导引律(nSDRE)。仿真显示,nSDRE是一种有效的导 引律,较广义理想比例导引律(GIPN)具有更好的导引品质,特别在目标机动时,nSDRE能更 好地应对目标机动引起的视线转率发散而导致脱靶的问题。
     

低温传感器测量数据基于MATLAB的 处理与分析

文章介绍基于MATLAB的低温温度传感器测量数据的处理及分析方法,通过线性最小二乘法和无约束非线性优化法中的拟合模型工具对数据进行拟合,其中包括多项式、双曲线、幂函数、Steinhart-Hart经验公式等拟合模型工具,以低温温度传感器电阻与温度数据的曲线拟合为例,开展了各种模型的曲线拟合比较分析。最后通过剩余标准差和相关系数两项指标评价不同模型的拟合效果。     

Hardware-in-the-loop simulations of GPS-based navigation and control for satellite formation flying

A relative navigation and formation control algorithm for satellite formation flying was developed, and a hardware-in-the-loop (HIL) simulation testbed was established and configured to evaluate this algorithm. The algorithm presented is a relative navigation estimation algorithm using double-difference carrier-phase and single-difference code measurements based on the extended Kalman filter (EKF). In addition, a state-dependent Riccati equation (SDRE) technique is utilized as a nonlinear controller for the formation control problem. The state-dependent coefficient (SDC) form is formulated to include nonlinearities in the relative dynamics. To evaluate the relative navigation and control algorithms developed, a closed-loop HIL testbed is configured. To demonstrate the performance of the testbed, a test formation flying scenario comprising formation acquisition and keeping in a low earth orbit (LEO) has been established. The relative navigation results from the closed-loop simulations show that a 3D RMS of 0.07 m can be achieved for position accuracy. The targeted leader–follower formation flying in the along-track separation of 100 m was maintained with a mean position error of approximately 0.2 m and a standard deviation of 0.9 m. The simulation results show that the HIL testbed is capable of successful demonstration of the GPS-based satellite autonomous formation flying mission.     

HTPB复合固体推进剂本构方程

本文将Schapery非线性粘弹本构模型和Perzyna非线性粘塑本构模型组合建     

基于Khokhlov-Zabolotskaya-Kuznetsov方程的声爆频域预测法

建立了一种快速预测声爆传播特性的频域方法,基于传统Khokhlov-Zabolotskaya-Kuznetsov(KZK)方程,描述声爆沿激波波阵面法线方向的传播.为了验证模型正确性,以NASA TD N-161中试件C为对象,首先基于不同网格,利用计算流体动力学(CFD)方法得到超声速流场;将CFD结果作傅里叶变换后代入模型方程,快速求解声爆传播特性.预测结果与NASA实验结果符合很好,研究表明:预测方法能够捕捉声爆的非线性传播,声衍射项使得声场在远场趋于轴对称分布,远距离传播后声能量集中于低频分量.      

泵间管气液两相流压力波传播速度数值研究

采用ANSYS CFX热相变模型数值模拟了泵间管气氧和液氧的混合冷凝过程.并应用整体平均两流体模型建立了气液两相流压力波的色散方程,系统研究了泵间管两相压力波的传播速度特性.计算表明:泵间管气液两相波速随空泡份额的增加先减小而后在弹状流区域略微增大,且两相波速随扰动角频率增加而增大.增加液氧流量或降低进口气氧温度,气液的冷凝效果增强,沿管流的空泡份额减小,气液两相波速增大.      

基于KZK方程的声爆频域预测法

建立了一种快速预测声爆传播特性的频域方法,基于传统Khohklov-Zabolotskaya-Kuznetsov (KZK)方程,描述声爆沿激波波阵面法线方向的传播。为了验证模型正确性,以NASA TD N-161中试件C为对象,首先基于不同网格,利用计算流体力学(CFD)方法得到超音速流场;将CFD结果作傅立叶变换后代入模型方程,快速求解声爆传播特性。预测结果与NASA实验结果符合很好,研究表明：预测方法能够捕捉声爆的非线性传播,声衍射项使得声场在远场趋于轴对称分布,远距离传播后声能量集中于低频分量。