求解非线性反问题的大范围收敛梯度正则化算法

基于同伦映射的思想,改进了求解非线性反问题的梯度正则化算法。通过路径跟踪有效地拓宽了梯度正则化算法求解的收敛范围。对于正则化参数的修正,通过引入拟Sigmoid函数,提出了一种下降速率可调的连续化参数修正方法,在保证迭代稳定的条件下,得到较好的计算效率,同时保证该算法具有很好的抵抗观测噪声能力。实际算例表明,该方法收敛范围宽,计算效率高,在存在较强观测噪声的条件下也能得到很好的反演结果。     

求解非线性动力系统周期解大范围收敛方法

对于多自由度非线性动力系统,提出一种求解周期解的大范围收敛方法,这种算法对处理非线性动力系统有较强的功能。结合数值延拓算法,为求解具有系统参数的非线性动力系统在整个系统参数范围内的周期解提供了有效的方法。     

一种基于边界积分方程的隐式翼型反设计算法

本文由边界元方法出发,将适用于单连通空间Laplace问题的边界积分方程推广到带环量的多连通空间中,并对离散边界积分方程中的矩阵元积分式解析化,以避免在翼型尾缘处尖点附近直接利用数值积分计算矩阵元导致的数值振荡,对于以翼型表面压力分布为收敛目标的反设计问题,利用Newton- Raphson迭代求解满足该目标压力的非线...     

翼型跨音速流动计算和外形设计的变域变分有限元方法

在气体动力学问题研究中经常会碰到诸如激波、翼型设计等未知界面问题。未知界面的存在为该类问题的理论分析和数值求解带来了很大困难。刘高联针对未知界面问题发展了一种变域变分有限元方法,该方法将未知界面看作是一个变化区域的边界,采用变域变分将未知界面结合在变分泛函中,使其与求解流场的控制方程结合起来,从而将未知界面的求解和流场的求解完全耦合进行,因而是一种处理未知界面的独特工具,极适合于气动外形的设计求解。本文运用变域变分有限元方法对翼型跨音速流动正、反命题进行了数值研究。由于在跨音速翼型绕流中存在激波,所以为了得到压缩激波解,采用了“人工密度”办法。几个算例均得到了满意的计算结果和设计结果,证明了本文方法的有效性和优越性。     

平面跨音速叶栅正、反混合问题的一种高效率解法

本文在文献[1]的基础上提出了平面跨音速叶栅正、反混合问题的一种高效率解法,它保持了文献[1]中方法的优点,避免了打靶,从而达到了扩大适用范围,减少了计算时间。 文中导出了以Von Mises坐标为自变量的流线控制方程,简单而有效地处理了速度的双值问题,作为考核,计算了两个跨音速平面叶栅和一个拉伐尔喷管的流场,结果令人满意,本方法具有精度高、收敛快的特点,即使对于624所跨音涡轮动叶中部叶栅中比较复杂的流场,也只需迭代50次,费时不到1分钟(DPS-8机),即可获得满意的结果。     

基于气动特性翼型参数化方法的适用性研究

翼型的优化设计在飞机的研发生产过程中占据相当重要的地位,因为翼型的几何外形会直接影响到飞机的整体气动性能．合理选择恰当的翼型表达方式是翼型优化设计成功的关键之一,其中,参数化方法在近年来得到了广泛应用．针对翼型参数化方法,提出“几何鲁棒性”和“性能鲁棒性”作为评判方法优劣的重要参考标准．选取三种常用的翼型参数化方法--NURBS、CST和PARSEC方法,七种不同类型的翼型,采用四项指标作为系统的评判准则来研究各方法的拟合精度与多点扰动下的几何鲁棒性,并利用Foil Design软件计算气动性能的相关数据,研究各方法在多点扰动下的性能鲁棒性．经对比分析,这三种参数化方法对大部分翼型均具有良好的表达能力,且不同程度地体现了一定的几何鲁棒性与性能鲁棒性．     

用间接边界元方法求解不可压缩孤立翼型绕流的一种库塔条件处理方法

本文对不可压二维翼型势流绕流的边界元法求解作了分析,在对基本方程与边界条件进行数值离散化时,将库塔条件代入基本方程。按本方法编制的计算机程序对若干算例进行了验算。结果表明,本文提出的方法是可靠的,该法计算简捷、方便,占用计算机内存少,具有实用意义。     

基于混沌多项式方法的层流超临界翼型稳健设计研究

开展层流超临界翼型稳健设计研究。针对传统的基于蒙特卡洛不确定分析方法的稳健设计系统计算量巨大、效率低下等问题,引入了混沌多项式方法,建立了新型稳健设计系统。采用RAE2822翼型稳健设计验证了设计系统的可靠性、高效性。结果表明,混沌多项式方法在保证设计精度的基础上设计效率提高了10倍左右。采用该系统进行了层流超临界翼型稳健设计研究,设计的结果表明:对于层流超临界翼型设计需要兼顾层流区长度、阻力发散特性;基于多目标稳健设计方法可以权衡设计点特性和阻力的波动特性;设计后翼型与原始翼型相比,设计点阻力减小了10cunts,阻力发散马赫数提高了0.01,在整个马赫数范围内维持了低的阻力特性。     

LES/RANS方法混合函数特性研究

将两方程k-ω SST湍流模型和Sagaut的混合尺度亚格子模型通过一个混合函数相结合, 构造出一种混合大涡/雷诺平均N-S方程模拟方法(hybird large eddy simulation/reynolds-averaged navier-stokes, Hybrid LES/RANS), 采用这种混合模拟方法结合5阶WENO格式对Ma=2.8平板湍流边界层进行了数值模拟, 并在计算区域上游入口处采用“回收/调节”方法生成湍流脉动边界条件, 通过考查RANS区域向LES区域的过渡参数及网格分辨率对这种混合模拟方法进行了评价. 计算结果表明: 该文采用的混合模拟方法可以捕捉到湍流边界层中的大尺度结构且入口边界层平均参数不会发生漂移, 混合函数应当将RANS区域和LES区域的过渡点设置在对数律层和尾迹律层的交界处, 而过渡应当迅速以获得正确的雷诺剪切应力分布, 在该文采用的模型及数值方法的条件下, 流向及展向的网格小至与Escudier混合长相当时, 能够获得可以接受的脉动速度的单点-二阶统计值.     

一种改进的位移反分析方法

为了发展可用于各种岩土工程及工程地质问题的一般的位移反分析系统方法, 本文在边界控制方法的基础上进行了改进。其目的是对于其各种复杂的系统问题确定最重要的力学参数。     

Navier-Stokes analysis of a circulation control airfoil

The two-dimensional, compressible, mass-averaged Navier-Stokes equations are used to investigate flows about a typical circulation control airfoil. The governing equations are solved using the implicit approximate-factorization algorithm of Beam-Warming with a modified algebraic eddy viscosity model. Results are compared with experimental data, and excellent agreement is obtained. The effects of different jet momentum coefficients and angles of attack on the flow are studied. The mechanism of genenating large lift by circulation control is discussed.     

Numerical investigation of transonic flow over deformable airfoil with plunging motion

In this article, the transonic inviscid flow over a deformable airfoil with plunging motion is studied numerically. A finite volume method based on the Roe scheme developed in a generalized coordinate is used along with an arbitrary Lagrangian-Eulerian method and a dynamic mesh algorithm to track the instantaneous position of the airfoil.The effects of different governing parameters such as the phase angle, the deformation amplitude, the initial angle of attack, the flapping frequency, and the Mach number on the unsteady flow field and aerodynamic coefficients are investigated in detail. The results show that maneuverability of the airfoil under various flow conditions is improved by the deformation. In addition, as the oscillation frequency of the airfoil increases, its aerodynamic performance is significantly improved.     

基于同伦技术的偶应力反问题求解

引入Bregman距离构造同伦函数,建立偶应力反问题的一种求解模式,可以对偶应力理论相关参数进行识别.利用有限元技术,建立了偶应力正/反问题数值求解模型,该模型不仅考虑了非均质的影响,而且也便于反问题的敏度分析.文中给出了相关的数值算例,并对信息误差和不同的函数形式计算结果做了初步探讨,得到满意的结果.数值算例表明该模...     

Unsteady transonic aerodynamic loadings on the airfoil caused by heaving,pitching oscillations and control surface

An implicit upwind finite volume solver for the Euler equations using the improved flux-splitting method is established and used to calculate the transonic flow past the airfoils with heaving, pitching oscillations and the control surface. Results are given for the NACA64A-10 airfoil which is in harmonic heaving and pitching oscillation and with the control surface in the transonic flow field. Some computational results are compared with the experiment data and the good agreements are shown in the paper.     

基于混合网格的低声爆反设计方法研究

低声爆设计方法已成为新一代军民用超声速飞机研制过程中必须解决的关键难题之一。针对传统SGD低声爆外形反设计方法无法对声爆近场非线性效应进行描述和分析的缺点,提出了利用CFD方法求解得到的声爆近场压力分布代替F函数进行低声爆反设计的方法。声爆近场预测采用点-点对接的结构/非结构混合网格,充分利用非结构网格对复杂外形适应性强和结构化网格计算效率高的优点。结果分析表明,基于改进后的低声爆反设计方法得到的方案在声爆超压以及感觉噪声级等方面都比基于原始SGD方法得到的方案有较大改善。     

A new Navier–Stokes inverse method based on mass‐averaged tangential velocity for blade design

This paper presents a novel viscous inverse method for blade design. In this inverse design method the mass‐averaged tangential velocity and the blade thickness are prescribed, and the corresponding blade profile is sought. The blade profile is then computed iteratively using the discrepancies between the prescribed mass‐averaged tangential velocity distribution and its calculated distribution on an initial blade. The re‐design of an axial rotor blade, starting from an initial arbitrary profile in subsonic flow regimes, demonstrates the merits and robustness of this approach. Copyright © 2008 John Wiley & Sons, Ltd.     

Robust design of natural laminar flow supercritical airfoil by multi-objective evolution method

A transonic, high Reynolds number natural laminar flow airfoil is designed and studied. The γ-θ transition model is combined with the shear stress transport(SST)k-w turbulence model to predict the transition region for a laminar-turbulent boundary layer. The non-uniform free-form deformation(NFFD) method based on the non-uniform rational B-spline(NURBS) basis function is introduced to the airfoil parameterization.The non-dominated sorting genetic algorithm-II(NSGA-II) is used as the search algorithm, and the surrogate model based on the Kriging models is introduced to improve the efficiency of the optimization system. The optimization system is set up based on the above technologies, and the robust design about the uncertainty of the Mach number is carried out for NASA0412 airfoil. The optimized airfoil is analyzed and compared with the original airfoil. The results show that natural laminar flow can be achieved on a supercritical airfoil to improve the aerodynamic characteristic of airfoils.     

Nonlinear effects caused by pulsed periodic supply of energy in the vicinity of a symmetric airfoil in a transonic flow

Changes in the structure of a transonic flow around a symmetric airfoil and a decrease in the wave drag of the latter, depending on the energy-supply period and on localization and shape of the energy-supply zone, are considered by means of the numerical solution of two-dimensional unsteady equations of gas dynamics. Energy addition to the gas ahead of the closing shock wave in an immediate vicinity of the contour in zones extended along the contour is found to significantly reduce the wave drag of the airfoil. The nature of this decrease in drag is clarified. The existence of a limiting frequency of energy supply is found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 64–71, May–June, 2006.     

A remark on the projection-3 method

We show that the continuous (in time) form of the projection-3 scheme proposed in Reference 2 is not a proper approximation of the unsteady Navier-Stokes equations. Hence, the projection-3 scheme and its variants are not appropriate for the numerical computation of the Navier-Stokes equations.