基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

一种适用于飞行器复杂外形的表面与空间网格生成方法

本文通过建立反映机身横截面外形特征的控制点概念，及在控制点间分配网格点数方式，给出了一种非常实用、交互式的表面网格生成方法，所生成的网格失真度小，物体保形好，生成适用于CFD计算需要的表面网络花费的时间少，适用性强。基于无限插值理论，引入本文提出的网格正交控制、物面法向量控制、加权平均光顺措施，给出了一种改进的空间网格生成方法，有效克服了传统无限插值网格生成方法在复杂外形网格生成方面的缺陷，有效改善了网格的法向疏密性、贴体性及周向网格的均匀性，对真实飞行器部件网格生成非常迅速。以某飞行器复杂机身表面与空间网格生成，及无人机翼身组合体空间网格生成为例，检验了本文网格生成方法的实用性和有效性。     

非定常预处理方法在倾转旋翼飞行器悬停状态气动干扰模拟中的应用

针对倾转旋翼飞行器的多种速度尺度并存的复杂流动问题，发展了能够统一求解同时包含不可压缩流动和可压缩流动的非定常预处理方法和一种考虑了流场物理信息传播方向的预处理远场边界条件。结合适用于倾转旋翼飞行器的复杂外形及大位移问题的非结构混合网格技术和重叠插值技术，建立了倾转旋翼飞行器的数值模拟方法。在此基础上，首先通过悬停算例验证了数值方法的有效性，确定了选择预处理截断参数的原则。然后探究了倾转旋翼飞行器直升机模式悬停状态的尾迹涡结构演化过程，研究了“喷泉效应”和“地面效应”对倾转旋翼气动性能的影响。计算结果表明，旋翼尾迹涡的耗散在不同方位角存在较大差异。“地面效应”抑制了旋翼尾迹的发展，提高了旋翼整个方位角区域的升力系数，而不仅限于机翼上方方位角区域。 “地面效应”作用强于“喷泉效应”作用，使旋翼相对于孤立旋翼拉力系数略微增大了5.26%，但同时整个倾转旋翼飞行器的总的向上的力只有孤立旋翼拉力的82.46%。结论表明预处理方法和预处理远场边界条件能够较好地模拟倾转旋翼飞行器悬停状态的干扰流场，兼顾了数值计算的耗散性和稳定性。气动干扰对倾转旋翼飞行器整体气动性能有重大影响。      

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

两种典型非结构网格变形方法特性对比研究

非定常流动仿真中常常遇到运动边界的情形,基于物理模型的线性弹簧法和基于数学方法的径向基函数(RBFs)插值法是实现计算网格随运动边界变化的两种主要的非结构网格变形方法。以NACA0012翼型的俯仰运动为例,对上述两种方法的变形特性进行较详细地定量分析研究,比较线性弹簧法和RBFs方法的CPU计算时间和最大变形能力,分析两种方法的网格质量随俯仰角度和运动步数的关系,并给出典型状态下网格单元质量云图和流场计算结果。结果表明:基于RBFs插值的网格变形方法计算效率较高,变形能力较强且能保持较高的网格质量,是一种高效的网格变形方法。     

基于线性弹性体的网格变形方法研究

网格变形方法在计算流体力学的各个方面都发挥着极其重要的作用,是各种应用问题解决方案的基本组成部分.研究和实现了基于线性弹性体的网格变形方法,对基于弹性体法的网格变形方法进行了综述,从弹性方程理论、高效求解方法等方面进行了介绍,给出了方法的最优参数选择,并从基本变形能力、多段翼型设计、非定常动态特性模拟等方面,展示基于线性弹性体网格变形方法的研究和应用进展.结果表明:实现的线弹性体网格变形方法特别善于处理几何大变形,尤其是网格扭转比较大的应用,并且变形后网格品质高、鲁棒性强.     

多面体网格在静气动弹性计算中的应用

介绍了多面体网格技术在静气动弹性问题中的应用.基于M6算例和径向基函数插值法,探究了多面体网格的CFD计算效率和变形能力,使用AGARD445.6机翼进行了静气动弹性数值计算.结果表明,多面体网格具有计算效率高、变形能力强、计算结果准确等优点,可适用于静气动弹性问题的数值计算.     

飞机着陆过程中非定常地面效应数值模拟方法

针对飞机着陆过程的大位移运动,提出适于多块结构动网格生成的改进谢别德-无限插值混合方法,改进谢别德方法准确计算各网格块的边界网格变形,无限插值方法快速计算内部网格点变形,结合局部重构保证动态网格的合理分布.混合方法有效解决了单一无限插值法存在的动网格交错问题.进而将二维非定常地面效应N-S方程计算方法发展至三维情形,应用于带增升装置的翼身组合体着陆过程中非定常地面效应计算分析.结果表明:模型的升力系数随离地高度的减小而减小,非定常和准定常计算方法的结果存在一定差异,而考虑了非定常效应的算法更为准确.     

一种交互式分区网格的生成方法和软件

介绍了一个具有交互性能的网格生成系统--BIGG系统。它是采用Windows环境下的C++语言与Fortran语言混合编程技术开发的,具有界面友好、显示灵活、交互性能强等特点。网格生成实例表明,采用BIGG软件可以缩短网格生成的周期,提高网格生成的质量。     

非结构化多重网格流场数值计算

本文对二维非结构网格无粘流场计算的多重网格方法作了探讨。研究了非结构化多重网格的生成方法,插值算子的构造以及流场计算．以绕流ＮＡＣＡ００１２翼型无粘流场为算例,研究了多重网格法在非结构网格上的实现方法及其加速效果表明本文的方法有显著的加速收敛效果。     

RBFs-MSA Hybrid Method for Mesh Deformation

Simulating unsteady flow phenomena involving moving boundaries is a challenging task,one key requirement of which is a reliable and fast algorithm to deform the computational mesh.Radial basis functions(RBFs) interpolation is a very simple and robust method to deform the mesh.However,the number of operations and the requirement of memory storage will be increased rapidly as the number of grid nodes increases,which limits the application of RBFs to three-dimensional(3D) moving mesh.Moving submesh approach(MSA) is an efficient method,but its robustness depends on the method used to deform the background mesh.A hybrid method which combines the benefits of MSA and RBFs interpolation,which is called RBFs-MSA,has been presented.This hybrid method is proved to be robust and efficient via several numerical examples.From the aspect of the quality of deforming meshes,this hybrid method is comparable with the RBFs interpolation;from the aspect of computing efficiency,one test case shows that RBFs-MSA is about two orders of magnitude faster than RBFs interpolation.For these benefits of RBFs-MSA,the new method is suitable for unsteady flow simulation which refers to boundaries movement.     

基于动态嵌套网格的飞行器外挂物投放的数值模拟

运用Euler方程和六自由度动力学方程数值模拟了外挂物投放的气动特性与运行姿态和轨迹。为了模拟外挂物从飞行器分离后的绕流特性和相对运动，采用了动态嵌套网格技术，并发展了一种洞边界建立和外边点的插值单元搜寻蛇行法，提高了计算效率。通过算例，证明了本文方法的有效性和合理性。     

Numerical simulation of aerodynamic interaction for a tilt rotor aircraft in helicopter mode

A rotor CFD solver is developed for simulating the aerodynamic interaction phenomenon among rotor, wing and fuselage of a tilt rotor aircraft in its helicopter mode. The unsteady Navier–Stokes equations are discretized in inertial frame and embedded grid system is adopted for describing the relative motion among blades and nacelle/wing/fuselage. A combination of multi-layer embedded grid and ‘‘extended hole fringe" technique is complemented in original grid system to tackle grid assembly difficulties arising from the narrow space among different aerodynamic components, and to improve the interpolation precision by decreasing the cell volume discrepancy among different grid blocks. An overall donor cell searching and automatic hole cutting technique is used for grid assembly, and the solution processes are speeded up by introduction of Open MP parallel method. Based on this solver, flow fields and aerodynamics of a tilt rotor aircraft in hover are simulated with several rotor collective angles, and the corresponding states of an isolated rotor and rotor/wing/fuselage model are also computed to obtain reference solution.Aerodynamic interference influences among the rotor and wing/fuselage/nacelle are analyzed,and some meaningful conclusions are drawn.     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

飞翼气动优化中参数化和网格变形技术

几何参数化和网格变形是飞行器气动外形数值优化迭代过程中的两个关键技术。基于非均匀有理B样条(NURBS)的自由型面变形(NFFD)技术对几何表示形式具有普适性,距离权函数(DWF)网格变形技术具有计算快速和网格拓扑无关性,两者广泛应用于曲面优化。基于NFFD技术改进了反求参数的Newton迭代算法,并通过提高物面附近网格刚度改进的距离权函数(IDWF)技术使其适用于更大程度的网格变形。还提出了改进后的参数化和网格变形两种技术并行计算的具体实现算法。结合离散伴随方法,使用参数化和网格变形技术,实现了由NACA0012初始翼型到飞翼标准翼型EH1590的反设计;针对某飞翼标模完成了单点全机升阻比优化,升阻比提高约18%。数值结果表明,建立的NFFD和IDWF动网格技术可满足飞翼气动外形优化参数化和快速网格变形的需求。     

应用Delaunay图映射与FFD技术的层流翼型气动优化设计

采用非均匀有理B样条(NURBS)基函数属性建立了任意空间的自由式变形(FFD)翼型参数化方法,进一步结合基于Delaunay图映射技术建立了结构对接网格变形模式,通过粒子群优化(PSO)算法进行参数化方法、网格变形模式以及计算流体力学(CFD)数值模拟技术之间的整合,研究、构建了气动优化设计系统,并对某型层流理念设计的高空长航时(HALE)飞机基本翼型进行气动优化设计。气动特性目标函数评估方法中,边界层转捩数值模拟技术采用γ-Re_θt转捩模型耦合剪切应力输运(SST)模式湍流模型。优化设计后翼型气动特性表明:采用相关技术建立的层流翼型气动优化设计系统对于层流理念设计的HALE飞机翼型的设计具备较高的优化效率。     

一种用于机动目标跟踪的新自适应卡尔曼滤波算法

在“当前”模型的概念下,从工程实现的背景出发,提出了一种用于机动目标跟踪的新自适应卡尔曼滤波算法。基本思想是通过对加速度项引入加权因子来进一步突出“当前”信息的作用,为“当前”模型提供更加准确的“当前”信息。蒙特卡罗模拟结果表明,算法不仅克服了“当前”模型自适应卡尔曼滤波算法的缺陷,而且使跟踪性能得到进一步的提高。     

高效气动优化设计方法

基于CFD方法开展气动外形优化设计通常计算量较大,采用离散共轭方法计算目标函数梯度,建立了高效的飞机气动优化设计系统。采用NURBS方法对翼剖面几何外形的扰动量进行参数化,避免了对原始外形的拟合,建立了基于NURBS方法的机翼剖面参数化方法。在复杂外形的优化设计中,动网格方法是关键技术之一,采用无限插值方法生成变形后的网格,并提出采用无限插值方法处理部件之间相贯线发生变化的情况。最后开展了翼身组合体+吊舱+挂架等复杂外形的优化设计,成功地减小了阻力,证明了方法的有效性。