基于格子Boltzmann方法的可压缩翼型绕流模拟

采用能够恢复可压缩Navier-Stokes方程的耦合双分布(DDF)格子Boltzmann方法(LBM)模拟了贴体网格下的NACA0012翼型绕流.首先在Ma＝0.5,α＝0.0,Re＝5 000状态下通过模拟低雷诺数流动检验了该方法;然后分别在Ma＝0.5,0.85和1.2,α＝-0.05,1.0和0.0状态下进行了模拟,并在Ma＝0.85时通过加密网格得到了更好的结果.为了避免转捩问题采用了无粘边界条件进行模拟,但得到的压力系数与实验结果吻合良好,证明了该方法的可靠性并显示了其模拟可压缩的潜力.该方法在翼型绕流上的应用为湍流流动的数值模拟提供了基础.     

改进的CFD/CSD耦合系统设计方法及其计算

论述了界面上的数据交换方法和时间推进格式在流体动力学（CFD）和结构动力学（CSD）耦合求解非线性气动弹性问题中的关键作用;改进了CFD／CSD耦合系统．包括基于边界元方法设计了CFD和CSD耦合界面的数据转换方法,该方法可在同一映射矩阵处理结构响应和非定常气动载荷转换,保证了耦合边界上的能量守恒;改进了一种松耦合方法流程,在保持模块化基础上提高了计算精度和效率。最后将本文方法应用于二维位移外插算例和A—GARD445．6机翼的动响应分析中,结果表明本文方法具有较高精度。     

航天器推进剂晃动与控制系统之间耦合效应的混合型估算方法(英文)

基于等效力学模型理论与计算流体力学技术,提出了一种用于估算航天器中推进剂晃动的力学影响的方法.以计算流体力学技术为依托,建立了晃动的等效力学模型并将其整合为晃动模块,以便在航天器导航控制回路的仿真分析中应用.该模块体现了航天器运动特征参数与推进剂晃动动力响应之间的映射关系,使其在用于分析晃动与控制系统之间的耦合效应时,比直接应用计算流体力学分析软件更为方便.该方法的正确性与优势通过检验算例得到了验证.通过在实例中应用,证明了该方法可方便地用于确定晃动的影响和检验数种防晃设计的有效性.     

一种模拟不可压流的节点间断有限元格子波尔兹曼方法

基于节点间断有限元方法,本文提出了一种求解格子波尔兹曼方程(Lattice Boltzmann equation,LBE)的新方法,即节点间断有限元-格子波尔兹曼方法(Nodal discontinuous Galerkin-Lattice Boltzmann method,NDGLBM)。在该方法中,LBE的碰撞过程和迁移过程被拆分成了两步:碰撞过程用LBM多松弛时间(Multiple relaxation time,MRT)模型进行求解,迁移过程则写成对流方程并采用节点间断有限元方法进行求解。其中,空间离散采用了非结构网格,时间离散采用了四阶、五步龙格-库塔格式。为了验证NDG-LBM的可行性,文中模拟了顶盖驱动方腔流、静止单个圆柱绕流、旋转-静止双圆柱绕流以及高雷诺数NACA0012翼型绕流。计算所得的数值结果与其他文献中的结果吻合度很好。     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

自适应计算网格的生成和优化方法

本文提出一种简单通用的自适应计算网格生成方法。配合文[1]的网格优化方法,能在生成自适应网格的同时,保证网格具有较好的光滑性和正交性。本文还阐述了单元加权函数W的取值和边界点的确定对自适应网格性质的影响。作为应用实例,本文将自适应网格用于二维无粘跨音速叶栅流场计算。结果表明,采用自适应网格能明显改善流场的计算精度。     

格子Boltzman方法的微间隙旋转流动研究

为了全面了解动压轴承的内部流动特性,使用了格子玻尔兹曼方法对内壁高速旋转的圆环间隙内的二维流动进行研究,在偏心率0.2～0.9,间隙比0.005～0.01,转速3×104～1.8×105 r/min的范围内讨论了流动和压力周向分布.通过研究发现,内壁的转速与压力极值的大小几乎成线性关系,而大偏心率和小间隙比则会带来压力...     

复杂固体边界流场的三维 PIV 测试技术:任意三维边界识别算法

针对复杂固体边界三维流场的 PIV 测试应用,以及流固耦合实验研究中流场和固体结构特征的瞬态同步测试需求,发展了一种基于双相机布置形式的任意三维边界识别算法以精确获取三维表面几何信息;并以基于MLOS-SMART 三维粒子场重构的 Tomo-PIV 算法计算三维速度矢量场,可同步获取三维表面结构运动/变形信息和三维瞬态速度场。这一边界识别算法基于 SURF(加速稳健特征)模式识别算法进行三维曲面重构,可以确定流场中三维物体结构的边界特征。论文采用双相机布置方式获取了三种不同曲率的圆柱曲面图像,验证了所发展的三维边界识别算法的准确性。最后以圆柱绕流 Tomo-PIV 数字合成粒子图像序列为验证对象,采用所发展的边界识别算法和 Tomo-PIV 算法分别高质量地计算出圆柱曲面信息和三维速度场。     

圆转方非轴对称喷管无粘流动的多层网格法求解

本文利用多层网格法求解三维无粘亚音和跨音速圆转方非轴对称喷管的流动,基本格式是把Ni[1]的二维有限体积积分格式推广到三维流场。边界条件采用文[2]中的“预测一修正”处理方法。为了进一步提高格式稳定性,在计算流动变量的一阶时间变化量时,采用单元体积平均法而不是文[1]中的算术平均法。引入当地时间步及多层网格方法,以加快格式的收敛速度。文中还对多层网格构造形式及其收敛效果作了讨论。通过对两种喷管模型流动的计算和与试验数据的对比,验证了该算法的精度。     

基于时间解析PIV的圆柱绕流尾迹特性研究

采用时间解析PIV（采样频率为1000Hz）在0.55m×0.4m声学风洞中测量了直径D=20mm圆柱后方7.5倍直径、圆柱两侧各3.3倍直径所围成范围内的绕流尾迹在雷诺数Re=2.74×10⁴下的非定常流场。针对PIV获得的速度场数据,进行流场和频谱特性分析,探讨了圆柱绕流尾迹中的平均流场和脉动流场特性,以及旋涡脱落的频率特性。提出了基于速度场之间相关性的相位平均分析方法,系统分析了圆柱上下两侧旋涡交替生成、脱落、发展并耗散的完整演化过程。结果表明：在圆柱后方存在一个低速回流区,其中心0.8D的位置附近是流动结构变化最剧烈的区域;圆柱后方1.9D位置附近是上/下两侧脱落旋涡交汇、耦合的区域,湍流脉动最强;圆柱绕流尾迹中,旋涡脱落频率对应的斯特劳哈尔数稳定在0.2左右;基于速度场之间相关性的相位平均分析方法简单有效,可以准确地识别绕流尾迹中旋涡交替脱落和发展的时空演化过程,在非定常流场测量方面具有普遍推广意义。     

格子玻尔兹曼和气体动理学通量算法及其应用进展

采用守恒律方程求解流体流动问题时,单元界面通量的计算尤为关键,该过程也被称为通量重构。由于离散控制方程的物理量定义在解点上,如何利用解点上的值来计算单元界面的通量,是计算流体力学最为关心的问题之一。针对该问题的研究已发展了各式各样的计算格式,例如完全基于数学重构的差分近似、基于部分物理重构的黎曼通量求解器以及近年发展起来的基于完全物理重构的气体动理学格式、格子Boltzmann通量算法和气体动理学通量算法。本文首先对几种典型的通量重构算法进展进行回顾和分析;然后着重介绍格子玻尔兹曼通量算法和气体动理学通量算法的研究进展及其相关应用;最后就该类算法存在的挑战和可能的研究方向进行展望。