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

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

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

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

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

非结构弹性体动网格方法应用策略及改进

针对3种代表性的经典弹性体动网格方法，从网格最大变形能力、质量评价指标和空间分布趋势3个方面综合分析了各方法的网格变形特性，重点总结了各方法的适用范围.结合工程实际，研究了特定问题的弹性体动网格应用策略.最后，整合各经典弹性体动网格方法的优势提出了改进方法，并通过NACA 0006翼型的超声速俯仰振荡加以验证.研究结果表明:基于面积的弹性体法适用于加密型网格，而基于距离的弹性体法适用于均匀型网格，该应用策略可为工程问题提供指导.改进方法效果明显，在大变形时相对基于距离的弹性体法的各项网格质量增幅均在10%以上，相对基于面积的弹性体法的最小网格质量增幅几乎达到4倍，并且在数值模拟中具有较高的精度，可显著增强弹性体法的鲁棒性.      

基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

基于非结构运动对接网格的旋翼前飞流场数值模拟

发展了一种基于线性弹簧理论的非结构运动对接网格方法,通过求解非定常Euler方程模拟了旋翼的前飞流场。利用运动区网格的转动模拟旋翼的旋转运动,并运用线性弹簧动态网格技术处理旋翼的挥舞和变距运动,方程求解中采用一种精度较高的插值和重构方法。算例表明,本文的方法能够有效地模拟旋翼的前飞流场。     

非结构动网格在三维可动边界问题中的应用研究

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法并通过耦合求解基于(Arb itrary Lagrangian-Eu lerian ALE)描述的Eu ler方程,模拟了做谐和振动的刚性无后掠机翼及柔性M6机翼的跨音速绕流,计算结果与实验结果吻合良好。     

非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

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

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

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.     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

大幅运动复杂构形扑翼动态网格生成的一种新方法

 基于Delaunay图映射的动态网格生成方法无需迭代计算,效率高,稳定性好。但对大幅运动复杂构形的动边界问题,背景图容易交叉,重新生成背景图和重新定位网格点信息不仅费时而且会导致网格质量的严重下降。提出一种双重Delaunay图映射的动网格生成方法,通过在初始背景图中添加辅助点,生成虚实两种背景图和虚实两种映射关系。分别根据虚映射关系和虚背景图、实映射关系和实背景图,移动辅助点和网格点。几个复杂构形的扑翼算例表明,双重图映射方法多付出极少的内存代价即可避免背景图交叉引起的问题,提高了动网格生成的效率和质量,增强了处理大变形复杂动边界问题的强健性。     

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

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

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

一种基于径向基函数和峰值选择法的高效网格变形技术

基于径向基函数的网格变形方法因其具有诸多优点,而被广泛应用于气动外形优化设计等领域。对于大规模网格或复杂构型,该方法所需计算量是难以承受的。为了提高网格变形效率,可以通过减少建立插值模型所需支撑点数目来实现。为此,提出一种高效的选点算法--峰值选择法。该算法在选点过程的每个迭代步中对边界节点处的误差进行分析,从物面节点中选取多个峰值点来更新支撑点集,减少迭代步数,提高选点效率。在该算法的基础上,实现了网格的高效变形。三段翼型的网格变形算例证明：该方法可以在保证网格质量的同时实现复杂网格的变形。以DLR-F6复杂模型（约1 000万网格）的刚性运动和弹性大变形为算例对该方法的变形效率和变形后网格质量做了进一步评估：当相对误差设置为5.0×10^-7时,在保证变形后网格质量的前提下,该方法变形效率最快比传统贪婪算法提高了13倍,其中在选点效率方面最快提高了31倍。