基于四边形积分子域的一致性高阶无网格法

近来提出的一致性高阶无网格法通过发展导数修正技术大幅度减少了所需积分点数目,并能够精确地通过分片试验,从而显著改善无网格计算的效率、精度和收敛性。然而,由于导数修正方程数目须与积分点数目相匹配,该方法仅限于使用三角形积分子域。在保留原有导数修正方程的基础上,提出了修正导数的共面条件,并据此建立补充方程,使得方程总数可匹配于所需的积分点数目,从而将一致性高阶无网格法方便地拓展到使用四边形积分子域。数值结果表明,本文方法精确地通过了分片试验,并展现出极好的计算精度、效率和收敛性。     

超声速气膜气动光学效应与Reynolds数相互关系实验研究

目前,随着相关项目研究的不断推进,如何在高Reynolds数下研究其对气动光学效应的影响成为重要命题.通过设计变Reynolds数气动光学效应实验平台,模拟的单位Reynolds数可以在7.2×106~2.2×108 m-1范围内变化.搭建的基于背景纹影（background oriented schlieren,BOS）的波前测试系统可以达到6 ns的时间分辨率.此系统测量的平凸透镜波前结果表明:实验测量结果与理论计算结果的误差在±4%以内.通过测量9种不同Reynolds数下的超声速气膜瞬态波前数据,分析结果表明:在高Reynolds数条件下,Reynolds数对于超声速气膜气动光学效应的影响比较明显,通过对实验数据进行函数拟合发现OPD_rms∝Re0.88,与推导结果OPD_rms∝Re0.9十分接近;利用小波分析方法研究了高Reynolds数条件下气动光学效应沿流向的分布特征,发现OPD_rms的低频部分（信号的主体）先降低后升高,但是高频部分的震荡幅度先升后降.分析认为OPD_rms的低频部分主要受到流场整体结构的影响,而高频部分更多地受到涡的空间分布影响.      

大型海上风力机尾迹区域风场分析

随着世界范围内海上风电场的不断兴起及海上风电场设计规模的日趋庞大,大型海上风力机尾迹区域风场特征的研究对于海上风电场的优化布局有着重要的指导意义。本文基于三维Navier-Stocks控制方程和适用于旋转流场分析的RNGk-ε湍流模型,采用滑移网格技术对美国可再生能源实验室(NREL)的5 MW海上风力机的性能及其尾迹区域的风场特征进行了较为系统的数值模拟。通过将不同风速下风力机输出功率的数值结果与NREL的设计参考数据进行对比,三维数值模型的有效性得到了很好地验证。此外,在此基础上进一步研究了大型海上风力机额定风速下及不同风轮转速下尾迹区域平均风速的分布特征,并得到了一系列具有参考价值的重要结论。     

适用于间断Galerkin方法的限制器研究

开发了一种适用于高精度间断Galerkin方法的斜率（多项式系数）限制器。与现有的斜率限制器不同,该限制器实施过程不考虑网格单元类型（三角形或四边形）,通过全微分方法构造新的多项式系数,因此,该限制器能够适用于各种类型网格——结构化网格、具有单一单元的非结构化网格和具有混合单元的非结构化网格。由于该限制器能够方便地应用于具有混合单元的非结构化网格,因此,本文使用的程序能够方便地求解具有复杂几何结构的流动问题。本文利用一些典型算例对其性能进行了验证,表明该限制器适用于不同类型的网格单元,能够在光滑解区保证高的精度,并能够在阊断区抑帛3非物理振荡。     

具有内部特征约束的四边形网格生成方法

提出一种具有内部特征约束的四边形网格自动生成方法,可以处理内部具有孔洞、约束线、约束点以及密度线、密度点等特征约束,满足数值分析中对网格生成的特殊要求。以区域分解法为基础,提出了确定最佳切割线的方法以及在切割线上生成过渡均匀网格节点的方法;将约束线、约束点以及密度线、密度点作为面积为零的孔洞,提出了对特征约束的处理方法;针对特征约束分隔区域的情况,提出了自动确定各子区域边界的方法。基于本文提出的方法已开发出具有特征约束的四边形网格自动生成程序,并成功应用于洪水分析系统和楼面设计分析系统。     

积分背景网格对自适应EFG法拓扑优化的影响研究

采用应力能量范数作为误差指标,探讨了EFG法中积分背景网格对计算精度的影响,得到了合理划分背景网格的建议;建立了以节点密度为设计变量、以最小化柔度为优化日标的拓扑优化模型。采用以节点密度值为加点判据的自适应规则加点方案,开展了连续体结构的拓扑优化研究,该加点方案能有效地减少设计变量的个数,探讨了背景网格对拓扑优化结果的影响。算例结果表明,采用合适的背景网格不仅能进一步减少设计变量的个数,而且能够改善拓扑优化结果的光滑性,使计算效率和精度得到提高。     

振动NACA0012翼型的移动网格数值模拟

利用流体力学方程的积分形式给出非结构移动网格上离散格式,利用自适应移动网格方法移动网格,进而得到网格速度.对振动Naca0012翼型问题,分三种类型确定网格速度,再结合Riemann问题的解法器构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高效、高分辨率的特点.     

使用AMR型背景网格的并行SPH方法

在分布式并行SPH方法中根据网格范围内粒子最大光滑长度调整网格尺度的AMR（Adaptive Mesh Refinement）背景网格,确保网格内粒子的相邻粒子位于同一个网格或者相邻的同尺寸网格范围内.与统一尺寸的背景网格相比,在光滑长度空间分布不均匀和随时间变化的情况下能够应用并减小搜索量,提高计算效率.数值结果表明,基于AMR型背景网格的并行SPH方法应用到粒子光滑长度变化很大的数值计算时能有效地降低搜索复杂度.     

液舱横荡非结构网格高阶格式数值模拟

提出一种液舱横荡数值模拟的方法,将气液两相交界面视为物理间断,通过高阶精度离散格式捕捉间断.根据NVD(Normalized Variable Diagram)实现非结构化网格上高精度离散格式,建立固定网格上自由表面运动模拟方法.在开发的非结构网格有限体积法求解器GTEA(General Transport Equation Analyzer)基础上,实现上述方法.首先对经典的溃坝过程进行模拟,并与文献结果对比验证方法和程序的可信度.对二维矩形液舱在不同激振频率时的横荡进行数值计算,并与实验以及商业软件CFX计算结果进行比较.结果表明方法和软件可以模拟自由面的翻卷、破碎运动现象,对距自由面较深点处流体载荷的计算结果与实验值符合较好,与商业软件CFX相比,在相同计算网格下,算法可以更好的计算次峰值,验证方法正确可行.     

振动桥梁CFD数值模拟的网格运动算法

针对非流线型桥梁断面模型,采用Delaunay背景网格方法、棱边弹簧及线性弹簧网格方法,研究非结构网格与结构网格条件下,计算域网格运动的特点及算法优劣.研究表明,在正常的桥梁断面振动振幅内,三种网格运动方法均能保持网格变形时不会失效,其中Delaunay背景网格方法耗时最少,线性弹簧网格方法获得网格质量最高;在大振幅状态下,结构化网格容易失效,线性弹簧网格方法可最大限度保证网格有效.