二维非结构网格Euler方程的LU-SGS算法

研究了将二维非结构三角形网格重新排序后，用LU-SGS隐式算法计算Euler方程的方法。分析了网格不平衡情况下该方法的可行性和应用效果。计算了不同翼型在不同流动条件下的二维Euler方程，并与四步Runge-Kutta显式方法进行了比较，结果验证了此方法具有良好的计算和收敛效果，以及在粘性计算方面的潜力。运用此方法对某预研型号的头部外形进行了设计和选择。该方法克服了以往隐式方法大量耗费内存的弱点，达到了计算耗时短和占用内存少的统一。     

迟滞非线性二元机翼颤振特性分析

 采用多项式迟滞非线性模型建立二元机翼气动弹性运动方程，并用数值积分法进行求解。通过系统响应振幅随来流速度变化的分叉图和频谱分析发现，俯仰方向由于含有非线性因素，振动中的高阶分量随速度提高不断增加，并引起高次分叉。重点研究“机翼/空气”质量比以及“沉浮/俯仰”两个自由度的自然振动频率比对非线性颤振速度边界的影响，并提出可以通过提高自然振动频率比来减小迟滞非线性因素的不利影响。     

时域有限体积法的格式稳定性分析

以模型方程为例,采用MUSCL插值的S-W通量分裂格式,分别结合二步、四步龙格一库塔时间推进格式建立数值求解方法.用Von Neumann傅里叶级数分析方法,通过对不同格式的放大因子及色散因子的计算,分析了时域有限体积法(FVTD)格式的稳定性.     

圆柱形霍尔推力器内等离子体数值模拟

利用一维磁流体动力学模型,对圆柱形霍尔推力器的放电等离子体进行了数值模拟。考虑了等离子体的电离、中和、碰撞、玻姆扩散及阳极鞘层的影响,由龙格-库塔方法得到离子速度、离子数密度、电子的温度等分布,其与实验结果有很好的一致性。经分析可知,此分布与电磁场的分布、粒子碰撞及电子阻抗等因素有关,并分析了电子温度分布、电子速度与电磁场的关系。结果表明,离子数密度沿通道方向增加,但在出口附近略有下降;而中性粒子数密度逐渐降低;离子速度在出口达到最大值,电子速度在下游有较大的梯度分布。     

点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

Simulation study of motion of charged particles trapped in Earth’s magnetosphere

This article aims to understand the motion of the charged particles trapped in the Earth’s inner magnetosphere. The emphasis is on identifying the numerical scheme, which is appropriate to characterize the trajectories of the charged particles of different energies that enter the Earth’s magnetosphere and get trap along the magnetic field lines. These particles perform three different periodic motions, namely: gyration, bounce, and azimuthal drift. However, often, the gyration of the particle is ignored, and only the guiding center of the particle is traced to reduce the computational time. It is because the simulation of all three motions (gyro, bounce, and drift) together needed a robust numerical scheme, which has less numerical dissipation. We have developed a three-dimensional test particle simulation model in which the relativistic equation of motion is solved numerically using the fourth and sixth-order Runge-Kutta methods. The stability of the simulation model is verified by checking the conservation of total kinetic energy and adiabatic invariants linked with each type of motion. We found that the sixth-order Runge-Kutta method is suitable to trace the complete trajectories of both proton and electron of a wide energy range, 5 keV to 250 MeV for L = 2 – 6. We have estimated the bounce and drift periods from the simulations, and they are found to be in good agreement with the theory. The study implies that a simulation model with sixth-order Runge-Kutta method can be applied to the time-vary, non-analytical form of magnetic configuration in future studies to understand the dynamics of charged particles trapped in Earth’s magnetosphere.     

卫星推进系统发动机启动过程数值仿真

以一维管道瞬变流理论和特征线方法为基础,建立了仿真系统部件流体动力学模型。计算方法采用带内插的特征线方法和四阶龙格库塔法。对燃烧室模型,考虑燃烧时滞的影响,计算了卫星推进系统发动机在启动过程中各参数的变化。仿真结果与卫星推进系统热试车试验结果基本吻合。结果表明计算模型较好地描述了发动机启动过程中的水击和管流振荡现象,采用的仿真方法符合精度要求,在工程上可为卫星推进系统的设计与试验提供指导。      

AUSM+方法在黏性内流模拟中的应用

简要阐述了三维雷诺平均Navier-Stokes(RANS)方程求解程序中使用的对流迎风分裂格式改进形式(AUSM+)以及Spalart-Allmaras(S-A)一方程湍流模型,然后对Sajben收缩扩散通道内跨声速流动和高速透平静叶排内部流动进行了数值模拟,并将数值结果与实验结果进行了比较,结果表明自主开发的三维R...     

轻负荷磁头的动态特性

本文采用有限差分法求解非定常修正Ｒｅｙｎｏｌｄｓ方程，采用四阶Ｒｕｎｇｅ－Ｋｕｔｔａ法求发磁头滑块的动力方程组，分别计算出磁头滑块受到冲击或飞越各种障碍物时的动态响应，较全面地反映了磁头滑块的动态时域特性。