一种非结构网格上基于径向基函数重构的ENO格式

基于二维Euler方程,对非结构三角形网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个三角形单元构造插值径向基函数,而在计算交界面的流通量采用两点高斯积分公式以保证格式的整体精度。时间离散采用三阶TVD Runge—Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

任意多面体网格上的欧拉方程数值算法

发展并验证了一种新的支持多面体网格的欧拉方程离散算法,采用Fortran 95编写了支持任意网格拓扑及三维复杂外形的通用求解器.在空间离散上,基于径向基函数理论提出了一种新的梯度计算方法,并采用动能通量分裂格式来得到准确并且稳定的二阶精度重构.该方法不依赖于网格几何形状,因此对网格质量不敏感.由于在时间方向采用了点隐松...     

基于物理量梯度修正的RBF数据传递方法

径向基函数是复杂外形飞行器多场耦合计算中重要的数据传递方法。复杂飞行器物理量分布通常与局部流动特征密切相关且在不同区域变化剧烈,如何在具备各向同性特征的径向基函数中考虑实际物理量各向异性分布因素的影响成为制约其精度进一步提升的关键。针对径向基函数的问题,提出了一种基于局部物理量梯度对径向基三方向权重进行自适应修正的数据传递新方法。采用高超声速舵面、翼身组合体算例对该方法在单次数据传递中的可行性及效果进行了验证,同时基于该方法对压缩拐角外形开展了气动力/热/结构多场耦合数值模拟,并与风洞试验结果进行了对比分析,以验证该方法在实际耦合问题计算长时间高频次数据传递中的适用性及可靠性。研究结果表明该方法能够很好地提升如薄板样条插值（TPS）、多重二次函数插值（MQ）等全域基函数的单次插值效果,并且有效地改进了MQ方法对形态参数的鲁棒性,而对于梯度修正后的紧支基函数,其改进效果更佳,以较少的点达到了TPS、MQ等全域基函数的插值效果,实现了数据传递效率和精度的兼顾。采用该方法获得的耦合计算结果与试验对比一致,表明了该方法的可用性以及在复杂外形多场耦合问题中的良好工程应用前景。     

基于径向基函数的发动机瞬时推力计算方法

为了快速准确插值计算发动机的瞬时推力,提出了使用基于径向基函数(RBF)的推力计算方法。通过算例证明,采用径向基方法的插值结果精度高,算法成熟,使用计算机编程容易实现,可以快速求解发动机瞬时推力插值问题。该方法在飞机设计、特别是性能计算及参数优化过程中可以广泛使用。     

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

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

基于径向基函数的优化代理模型应用研究

针对已有径向基函数建模方法在小样本情况下存在的模型精度低,误差大的缺点,本文提出了基于交叉验证技术的径向基函数代理模型构造方法,即在计算径向距离时引入权值系数,通过对其优化以减小模型的预测误差平方和,从而提高模型的精度。采用两种测试函数对本文所提出的方法进行验证,结果证明了本方法相对于传统的径向基函数建模方法可以有效提高模型精度,减小模型误差。最后,将其应用于传感器飞机的总体优化设计之中,结果表明本文所提出的方法是可行的。     

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

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

非结构嵌套网格中的一种改进型径向基函数插值方法

非结构混合网格拓扑关系相对复杂,在嵌套边界区域进行流场信息插值传递容易有精度的损失。利用已知的解析函数作为插值样本,对线性插值(LINE)、距离倒数权重插值(WA)、逆二次径向基函数插值(IQ)、Wendland暞sC2径向基函数插值四种方法的精度进行对比,筛选出精度较高的C2径向基函数;提出一种改进的C2径向基函数插值方法,通过调整径向基函数的作用半径来控制插值矩阵条件数,进而消除使用径向基函数插值在嵌套边界区域出现的数值奇性问题;选用MD30P/30N 多段翼与三维AEDC典型外挂物分离模型对该插值方法进行验证。结果表明:改进后的径向基函数插值方法能够有效消除嵌套边界区域数值奇性,同时计算收敛速率更快,插值精度更高。     

两种高精度高效率静气弹数据传递函数

遵循静气弹问题中气动/结构网格之间的数据传递要求,引入了基于径向基函数的流固耦合界面数据传递方法,确定了紧支半径的选取原则,并编制相应的计算程序。以M6机翼为例,将不同径向基函数的插值精度及计算效率进行比较,结果表明体样条（Volume Spline）函数和C2函数具有良好的插值精度和计算效率。这两种流固耦合界面数据传递方法非常适用于工程应用。     

机翼-机身-短舱-挂架外形气动优化设计方法

采用基于径向基函数的无限插值方法进行了复杂外形的动网格生成,并针对其存在的可能导致壁面附近出现负体积与网格质量下降的问题,提出对不与物面直接相连的边的位移采用径向基函数插值,对与物面相连的边的位移进行线性插值获得位移量的方法解决该问题.采用无限插值方法建立了一种简单有效的物面相贯线处理方法.利用离散共轭方法计算目标函数梯度,对DLR-F6 机翼-机身-短舱-挂架外形的机翼和短舱进行了几何外形参数化与气动外形优化设计.结果表明:考虑短舱安装方式的优化设计较不考虑短舱安装方式的优化设计可降低大约0.0001的阻力系数.通过全机优化设计,将全机阻力系数降低了0.00153.      

Nonlinear deforming of laminated composite shells of revolution under finite deflections and normal’s rotation angles

The discretization of the boundary value problem for laminated composite shells is based on the finite difference approach using the regular mesh with the constant grid step and the difference operators of the second order of accuracy. The dynamic relaxation method is proposed for the solution of the nonlinear problem. The evolutionary equations of the dynamic relaxation are constructed, and the optimum parameters of the converging linear iterative process are estimated.     

基于四面体单元的三维裂纹自动化建模

三维裂纹前缘布置结构化网格是一项耗时且困难的工作。以二次四面体单元虚拟裂纹闭合法（VC‐CM）为基础,提出基于几何特征的三维裂纹建模方法,阐述应用四面体单元在ABAQUS上进行三维裂纹建模的实现细节;应用Python编程语言结合NumPy和SciPy开源科学计算库开发参数化三维裂纹分析程序包,实现三维裂纹问题分析的全流...     

Solid rocket motor propellant grain burnback simulation based on fast minimum distance function calculation and improved marching tetrahedron method

To efficiently compute arbitrary propellant grain evolution of the burning surface with uniform and non-uniform burning rate for solid rocket motor, a unified framework of burning surface regression simulation has been developed based on minimum distance function. In order to speed up the computation of the mini-mum distance between grid nodes of grain and the triangular mesh of burning surface, a fast distance querying method based on the equal size cube voxel structure was employed. An improved marching tetrahedron method based on piecewise linear approximation was carried out on second-order tetrahedral elements, achieved high-efficiency and adequate accuracy of burning surface extraction simultaneously. The cases of star grain, finocyl grain, and non-uniform tube grain were studied to verify the proposed method. The observed result indicates that the grain burnback computation method could realize the accurate simulation on unstructured tetrahedral mesh with a desirable performance on computational time.     

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

三维非结构网格的生成及优化

叙述了一种新的生成三维非结构网格的Delaunay方法,该方法能够高效地生成优质四面体网格。由于采用了局部剖分的方法来保证边界的完整性,因而不需要添加新的节点到原始表面三角形中。该方法在四面体外接球的球心处插入新的节点,在插入全部内点后,采用删除边、边/面交换、网格光顺等方法对生成的网格进行了优化,基本能删除质量差的四面体。     

基于高阶谐波平衡的跨声速颤振高效预测方法

跨声速流场激波及其诱导的附面层分离等非线性因素导致跨声速颤振边界很难被准确预测,尤其是目前工程常用的偶极子格网法,在跨声速时该方法的预测精度大幅下降。在雷诺平均Navier-Stokes方程流场求解器的框架内,利用结构模态建立广义结构运动方程,利用径向基函数建立模态振型的插值方法,结合径向基函数和无限插值两种网格变形方法的优点实现高效高鲁棒性网格变形方法,从而实现颤振时间推进分析流程,利用国际颤振标模AGARD445.6机翼验证程序在跨声速颤振边界预测中的可靠性。然而,时域方法在气动/结构反复迭代,需要耗费大量的计算资源和时间。为了提高颤振预测效率,基于高阶谐波平衡（HOHB）方法快速获得广义力影响系数矩阵,利用该矩阵建立频域模态位移和气动力之间的关系,实现高效颤振频域分析方法。通过二维翼型和三维机翼算例进行验证,结果表明：在不对计算精度产生明显影响的前提下,HOHB方法能够提高颤振预测效率约6倍。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

Delaunay graph-based moving mesh method with damping functions

The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.