浸没边界-格子玻尔兹曼方法的GPU并行加速

为了高效地模拟大规模三维流固耦合问题,提出了一种在GPU平台上并行加速浸没边界-格子玻尔兹曼(immersed boundary-lattice Boltzmann,IB-LB)耦合方法的策略。为了验证准确性和加速效果,计算了不同攻角下的三维平板绕流问题,得到了与实验结果吻合良好的升力和阻力系数随攻角的变化规律,并发现GPU内核函数的性能均受制于带宽,且其带宽利用率均达到70%。为了展示计算能力,模拟了中等雷诺数下三维拍动平板的绕流问题,证明所提方法和加速策略在模拟大规模复杂流固耦合问题上的有效性。

三维间断伽辽金玻尔兹曼方法的完全匹配层吸收边界条件研究

为解决计算气动声学数值模拟中远场边界反射问题,在三维间断伽辽金有限元玻尔兹曼方法框架下应用完全匹配层算法构造了吸收边界条件,并采用了两种不同的离散求解公式。利用三维高斯脉动源数值案例测试了两种公式的有效性和稳定性,发现仅有一种公式既能有效衰减反射波又具有良好的稳定性;接着研究了影响完全匹配层算法无反射性能的若干因素,其中无反射性能采用归一化误差范数来衡量。研究表明:衰减因子总是存在最优值,且最优值仅与吸收层的参数有关,而与高斯脉动源自身参数无关,故所定义的归一化误差范数具有一定的普适性;当衰减因子采用幂函数分布律时,幂指数取2时无反射性能更好,取4时无反射性能随衰减因子变化更平缓;吸收层厚度和区域半径越大,无反射性能越好,但相应的计算量也显著提高。研究工作可为构造更加实用有效的无反射边界条件提供参考。

计算气动声学中的伽辽金玻尔兹曼方法研究

为获得气动声学的高精度和低耗散特性的数值方法,发展了伽辽金玻尔兹曼方法和相应的无反射边界条件。首先,引入新粒子分布函数到格子玻尔兹曼BGK方程中并重构欧拉方程;然后,在空间上采用高精度的交点间断伽辽金有限元方法,在时间上采用显式五级四阶龙格库塔离散方法对解耦得到的对流步方程进行离散求解;最后,通过数值通量构造速度边界、声学硬壁面边界和无反射边界条件。采用包含声反射和多普勒效应的数值算例进行验证,可得模拟值与解析解吻合一致,从而证明了伽辽金玻尔兹曼方法和无反射边界条件用于气动声学计算的有效性和准确性。

一种稳定的浸没边界-格子波尔兹曼方法

基于强制边界条件的浸没边界-格子波尔兹曼方法(Immersed boundary-lattice Boltzmann method,IB-LBM)提出了一种稳定的IB-LBM。在该方法中,采用了具有良好稳定性的多松弛时间(Multi-relaxation time,MRT)模型代替原来的单松弛时间(Single-relaxation time,SRT)模型。因而,在保证无滑移边界条件能够精确满足的同时,该方法的稳定性和适用性都得到了提高。通过对中高雷诺数下二维圆柱绕流算例的模拟,计算所得的数值结果与其他文献中的结果吻合,从而验证了该方法的可行性。

沙纹床面振荡流边界层的三维格子玻尔兹曼模拟研究

将壁面自适应局部涡黏(WALE)模型引入三维格子玻尔兹曼模型,对沙纹床面上的振荡流边界层运动进行了模拟,并与Smagorinsky模型的模拟结果进行了比较.结果表明,在振荡流沙纹床面模拟中,WALE模型与Smagorinsky模型均获得与已有实验结果吻合较好的垂向流速剖面,但WALE模型可以更好地反映粗糙紊流的特性.引入WALE模型的LB模型模拟得到的沙纹床面形状摩阻系数和肤面摩阻系数结果与前人数模结果一致,床面总摩阻力中形状摩阻相比肤面摩阻占主要部分.该模型可进一步应用于振荡流作用下沙纹床面紊流边界层运动特性的系统研究.

一种改进的格子玻尔兹曼方法及其在流固耦合传热问题中的应用

含热源颗粒的方腔自然对流是一种常见的流固耦合热流动问题.本文提出一种使用分布函数修正的浸没边界格子玻尔兹曼方法来处理流固耦合和流场及温度场耦合问题,并把该方法运用到含移动热源颗粒的方腔自然对流的数值模拟中,分析了方腔内的速度,温度,热源颗粒的轨迹和受力,探究其热源颗粒对传热机理的改善,解决了含热源的自然对流情况下流固耦合传热问题.

基于改进格子玻尔兹曼方法的砰击涌浪模拟

针对船舶与落石入水砰击涌浪灾害问题,采用浸入式边界-格子玻尔兹曼方法对二维楔形块入水涌浪问题进行了研究。从格子玻尔兹曼方法、基于流体体积法(VOF)的自由液面模型以及浸入式边界法的原理出发,阐述了其可以用来模拟有自由液面的流固耦合问题的理论基础,并优化了界面格子处浸入式边界力的求解方法。采用此方法模拟了楔形块体入水抨击实验,通过对模拟的块体速度与位置和实验结果进行比较,两者吻合良好,论文此方法模拟有自由液面的流固耦合问题具有较好的可行性。

浸没边界格子Boltzmann方法的改进及转动圆柱绕流模拟

针对浸没边界格子Boltzmann方法计算效率不足的问题,提出一种改进浸没边界格子Boltzmann方法。通过使用作用力格子Boltzmann模型,并简化流固耦合作用力的计算方法,使计算流程得到简化,同时,对比原方法转动圆柱绕流计算结果,计算时间缩短一半以上,提高了数值计算效率。结合改进算法研究转动圆柱绕流的流场升阻力系数、压力系数、流场速度等流体参数场随转速比的关系,揭示了一定的转动圆柱绕流动力学规律。研究结果表明,改进方法准确可靠。

基于格子玻尔兹曼法的三维溃坝流数值模拟与验证

溃坝问题一直是验证数值模拟方法有效性的经典问题。利用基于格子玻尔兹曼离散方法研发的XFlow软件来模拟三维溃坝问题,以此验证其有效性并提出进一步的改进方向。为了对比前人所做的实验和模拟结果,计算了无障碍三维溃坝流相应监测点的压力、水位等要素,得到了水流演进过程,并提出了用增强固壁函数法处理固壁边界以增强计算的稳定性,另外还考虑了有矩形障碍物的三维溃坝问题。研究结果表明:该方法能捕捉到某一时间点的自由液面,相关数据计算结果准确,且无需划分网格,能减少工作量,但结果稳定性较差,合理运用边界处理方法能有效解决这一问题。

基于改进的浸没边界-格子Boltzmann方法的圆柱绕流仿真计算

基于Suzuki等人提出的浸没边界-格子Boltzmann方法,采用作用力格子Boltzmann模型,提出了一种改进的浸没格子Boltzmann算法。对比相同计算循环次数的圆柱绕流仿真结果发现,改进算法通过提高流场速度分布函数的计算精度来简化流固边界表面作用力计算方法,避免了使用迭代计算,简化了计算流程,使得在低雷诺数定常流动圆柱绕流仿真时,圆柱绕流流线的对称性得到了保证,计算更为稳定,且计算效率得到了显著提高。在定常流与非定常流的圆柱绕流仿真应用中,论文计算结果与其他文献中的结果吻合较好,表明了改进算法的有效性。

基于隐式扩散的直接力格式浸没边界格子Boltzmann方法

采用浸没边界格子Boltzmann (immersed boundary-lattice Boltzmann, IB-LB)模型执行动边界绕流数值模拟时,信息交互界面和边界力计算格式直接影响流动求解器的数值精度和计算效率.基于隐式扩散界面,一种改进的直接力格式IB-LB模型被提出.边界力表达式基于欧拉/拉格朗日变量同一性准则推导,转换矩阵描述的信息交互界面耦合了拉格朗日节点间的非同步运动.采用Richardson迭代数值求解关联边界力与无滑移速度约束的线性方程组,不仅克服了传统速度修正格式中矩阵求逆引起的计算效率问题,而且摆脱了算法稳定性对拉格朗日点分布的依赖.根据解析解已知的Taylor-Green涡流评估本文模型的数值模拟精度,结果表明改进的IB模型能够完整保留背景LB模型的二阶数值精度.静止圆柱和振荡圆柱绕流数值实验结果表明,当前模型在涉及复杂外形和运动界面的流动模拟中能够提供可靠的数值预测,满足力同一性的IB-LB模型能够有效抑制非定常流体力的伪物理震荡.波动翼型绕流模拟验证了当前模型的实用性,可在大变形柔性体流固耦合动力学问题中进一步推广.

浸没边界–简化热格子Boltzmann方法研究及其应用

针对流固耦合传热问题,本文提出了一种基于浸没边界–简化热格子玻尔兹曼方法 (immersed boundary method–simplified thermal lattice Boltzmann method,IB-STLBM)的耦合模型.不同于传统的格子玻尔兹曼方法使用分布函数演化流场和温度场,简化热格子玻尔兹曼方法 (simplified thermal lattice Boltzmann method,STLBM)的演化过程不需要依赖分布函数,只涉及平衡态分布函数和非平衡态分布函数,能够直接演化宏观量,极大减小了计算过程中所占用的虚拟内存,简化了边界条件的实现方式,同时具有较高的稳定性.传统的浸没边界法对流场的计算采用欧拉网格,对固体边界采用拉格朗日网格,认为固体边界是对流场产生某种体积力.在应用浸没边界法时,汲取介观的思想,把固体的介入看作是对流场的干扰,打破了固体附近流体介观微团颗粒原始的平衡状态,这种干扰可以看作是在耦合边界上产生的一个非平衡项,可用非平衡态分布函数来表示.基于此,在模型中浸没边界法与简化热格子玻尔兹曼方法更紧密联系在一起,更大程度发挥二者的优点,整个计算过程更加简单直观,符合物理特性.通过对热圆柱绕流和内含热颗粒的封闭方腔自然对流问题的模拟以及对其结果的分析,验证了该算法在求解流固耦合传热问题的有效性和可行性.

基于曲线坐标系浸没边界方法的折流燃烧室模拟

航空发动机折流燃烧室几何结构复杂,其高保真数值模拟需要高效的网格与边界条件处理方法。采用曲线坐标系隐式浸没边界方法结合大涡模拟-概率密度函数输运方程湍流燃烧模型开发自研软件,并实现WP11中折流燃烧室的高保真模拟。流动模拟中准确解析了该燃烧室中的三股主要气流,且三股气流分别约占进口流量的75%,12.5%和12.5%。两相燃烧模拟中针对拉格朗日框架下的液滴运动和欧拉框架下的湍流燃烧采用不同网格标记,模拟得到的出口径向温度分布规律与实验一致,平均相对误差为17.95%,表明基于本方法开发的自研软件能准确模拟折流燃烧室中的两相湍流燃烧现象。

利用高精度浸没边界-格子Boltzmann流固耦合格式模拟2维刚体自由沉降运动

自由刚体的流固耦合问题是计算流体动力学领域的热点问题,也是实际工程应用中的常见问题。为了利用基于迭代外力的浸没边界-格子Boltzmann(immersed boundary-lattice Boltzmann,IB-LB)耦合格式计算自由刚体的流固耦合问题,首先在原控制方程中添加自由刚体的运动方程组,然后利用IB-LB耦合格式离散并数值求解该数学模型。以两小球在水箱中的自由沉降运动为例,比较提出方法与其他流固耦合算法的计算结果。最后模拟富有挑战性的矩形平板在水箱中的飘动过程,捕捉到与实验数据相吻合的刚体运动轨迹和流态结构,通过对比提出方法和现有同类算法的单步平均耗时与平板平均速度,证明推广后的算法在复杂流固耦合问题模拟上有较高的计算效率和精度。

基于浸没边界法的水库变动水面模拟及验证

水库水面波动具有复杂的流动边界。对复杂边界条件的处理,将直接影响水资源的开发利用效果。为模拟真实水库自由液面的大幅变动,利用浸没边界法来处理水库自由液面的变动问题,比较水库水温分布、水库下泄水温和坝前垂向水温分布,并通过简单算例予以验证。结果表明:该方法具有良好的可行性和准确性,改进了水库变动水面预测方法,可用于研究水库水面的时空变化规律。

曲壁面两相流的简单扩散界面浸没边界法

针对在欧拉-拉格朗日网格框架下曲边界上的两相流润湿边界条件实施问题,本文提出了一种简单扩散界面浸没边界法。在此方法的框架下,静止和移动曲边界的湿润边界条件都可以采用平板边界方式处理。在该方法中,我们在实际边界内外一个拉格朗日网格间距处添加了两层虚拟贴体拉格朗日网格,并通过插值得到三层拉格朗日点上的序参数,从而得到实际的曲边界润湿条件。该方案实施思路简单,不仅适用于静止曲边界,还适用于复杂的运动曲边界。为了验证该方法的有效性,本文在广义守恒相场格子玻尔兹曼方法(LBM)框架下,模拟了静态液滴在圆柱上的铺展、动态液滴撞击圆柱、圆柱在水池中沉降三个问题。结果表明该方法能够有效地模拟曲面上的两相流动,具有良好的精度、鲁棒性和质量守恒性。

格子玻尔兹曼方法中的特征无反射边界条件研究

采用多尺度展开技术从耦合双分布函数方程中恢复出守恒型的连续、动量和能量方程。在边界上取分布函数的平衡态形式恢复出欧拉方程,方便地引入了特征无反射边界条件。在此基础上采用二阶非平衡态外推格式来考虑黏性和传热对声传播的影响。通过对四种经典的数值算例的考查,特征无反射边界条件均能避免边界上伪反射波的产生。格子玻尔兹曼方法与特征无反射边界条件的结合能够准确地预测噪声的传播。

基于时驱硬球算法与格子玻尔兹曼方法的颗粒流体系统直接数值模拟

实现了一种直接数值模拟颗粒流体系统的耦合算法,颗粒间相互作用由时驱硬球算法描述,而流体的控制方程采用格子玻尔兹曼方法求解,流固耦合用浸入运动边界法实现.该方法使用欧拉网格求解流场,拉格朗日网格跟踪颗粒,避免了非结构化贴体网格方法需要重新划分网格的问题.通过模拟两个圆形颗粒在黏性流体中的沉降过程,成功地复现了经典的Drafting-Kissing-Tumbling(DKT)过程,验证了耦合算法的有效性.

粗糙床面振荡紊流边界层运动的三维格子玻尔兹曼模拟

该文采用格子玻尔兹曼(LB)方法建立了紧密排列颗粒粗糙床面振荡紊流边界层模拟的三维数学模型。收集粗糙床面波浪水槽实验结果,对水槽实验条件下的振荡流边界层进行模拟,模拟得到的平均流速、紊动特性、等效粗糙高度、摩阻系数和边界层厚度等与实验结果进行了全面对比,并对三维涡旋结构进行了分析。结果表明,所建立的LB数学模型可以合理模拟粗糙床面振荡流边界层流动,能够反映三维涡旋特性,可进一步应用于粗糙床面振荡流紊动特性以及床面粗糙参数的研究。

三角形布置圆柱体群绕流特性与流动机理研究

基于浸没边界-多松弛-格子玻尔兹曼方法,对Re=100时等边三角形布置三圆柱体群绕流问题进行数值模拟,分析间距比(Kd)与来流角度(α)两个关键参数对流体力系数及流场特性的影响,揭示其流动机理。计算结果表明:不同流场下尾流模式分为单漩涡、单漩涡体向双漩涡过渡、不规则的双漩涡、规则的双漩涡、双漩涡向三漩涡过渡、不规则的三漩涡和规则的三漩涡模式,在α=30°时流场下游涡街更加规则,尾流模式转变更快。三圆柱所受的时均流体力系数受间隙流流速影响较大,各来流角度下随间距比(除小间距比外)的变化趋势基本一致,尾流模式转变会导致三圆柱群的流体力系数均方根值发生显著变化,在Kd≥3.5后上游圆柱流体力系数均方根值有较大起伏。除α=0°工况外,在Kd≥2.5后中游圆柱的涡脱落频率受α和Kd影响较小,数值相对稳定。