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

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

人工释能法求解线性代数方程组

本文建立求解线性代数方程组的人工释能法迭代格式,并给出该迭代格式收敛的充分必要条件和最佳释能时机,同时给出实现人工释能法的数值计算方法及其数值稳定性条件。     

Automatic differentiation based discrete adjoint method for aerodynamic design optimization on unstructured meshes

A systematic methodology for formulating, implementing, solving and verifying discrete adjoint of the compressible Reynolds-averaged Navier-Stokes (RANS) equations for aerodynamic design optimization on unstructured meshes is proposed. First, a general adjoint formulation is con-structed for the entire optimization problem, including parameterization, mesh deformation, flow solution and computation of the objective function, which is followed by detailed formulations of matrix-vector products arising in the adjoint model. According to this formulation, procedural components of implementing the required matrix-vector products are generated by means of auto-matic differentiation (AD) in a structured and modular manner. Furthermore, a duality-preserving iterative algorithm is employed to solve flow adjoint equations arising in the adjoint model, ensur-ing identical convergence rates for the tangent and the adjoint models. A three-step strategy is adopted to verify the adjoint computation. The proposed method has several remarkable features:the use of AD techniques avoids tedious and error-prone manual derivation and programming;duality is strictly preserved so that consistent and highly accurate discrete sensitivities can be obtained; and comparable efficiency to hand-coded implementation can be achieved. Upon the cur-rent discrete adjoint method, a gradient-based optimization framework has been developed and applied to a drag reduction problem.     

A coupling VWM/CFD/CSD method for rotor airload prediction

A coupling fluid-structure method with a combination of viscous wake model(VWM),computational fluid dynamics(CFD) and comprehensive structural dynamics(CSD) modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method' is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.     

适于直升机配平计算的混合遗传算法

结合传统的数值计算如拟牛顿法的局部快速收敛特性,以及遗传算法的全局收敛与群体搜索的优点,发展出一种新的混合遗传算法,使之具备全局快速收敛特性。通过两个算例及某型直升机在前进比为0.2状态下的配平,验证了该算法的可行性及高效性,从而为直升机的配平计算提供了一种新的成功的算法,同时也在传统数值分析方法和现代智能进化算法之间的有效结合做了一次比较成功的尝试。     

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

非线性边界条件的Burgers－KdV方程的静态解

对几类非线性边值条件，讨论了Ｂｕｒｇｅｒｓ－ＫｄＶ方程的边值问题的静态解的存在性和解的个数．借助不动点定理，证明了该边值问题至少有一个，至多有有限个静态解．     

The use of Gaussian equations of motions of a satellite for local gravity anomaly recovery

The orbital elements of a low Earth orbiting satellite and their velocities can be used for local determination of gravity anomaly. The important issue is to find direct relations among the anomalies and these parameters. Here, a primary theoretical study is presented for this purpose. The Gaussian equations of motion of a satellite are used to develop integral formulas for recovering the gravity anomalies. The behaviour of kernels of the integrals are investigated for a two-month simulated orbit similar to that of the Gravity field and steady-state ocean circulation explorer (GOCE) mission over Fennoscandia. Numerical investigations show that the integral formulas have neither isotropic nor well-behaved kernels. In such a case, gravity anomaly recovery is not successful due to large spatial truncation error of the integral formulas. Reformulation of the problem by combining the orbital elements and their velocities leads to an integral with a well-behaved kernel which is suitable for our purpose. Also based on these combinations some general relations among the orbital elements and their velocities are obtained which can be used for validation of orbital parameters and their velocities.     

Friction moment analysis of space gyroscope bearing with ribbon cage under ultra-low oscillatory motion

This paper presents the model of calculating the total friction moment of space gyroscope ball bearings which usually work under ultra-low oscillatory motion and are very sensitive to the friction moment. The aim is to know the proportion of the friction moment caused by each frictional source in the bearing＇s total friction moment, which is helpful to optimize the bearing design to deduce the friction moment. In the model, the cage dynamic equations considering six degree-of-freedom and the balls dynamic equations considering two degree-of-freedom were solved.The good trends with different loads between the measured friction moments and computational results prove that the model under constant rate was validated. The computational results show that when the speed was set at 5 r/min, the bearing＇s maximum total friction moment when oscillation occurred was obviously larger than that occurred at a constant rate. At the onset of each oscillatory motion, the proportion of the friction moment caused by cage in the bearing＇s total friction moment was very high, and it increased with the increasing speed. The analyses of different cage thicknesses and different clearances between cage pocket and ball show that smaller thickness and clearance were preferred.     

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。