【摘 要】
：
本文提出了一种梯度光滑法(Gradient Smoothing Method,GSM)与格子Boltzmann的耦合数值方法。基于宏观与介观两种流体计算方法的耦合,该方法能够实现多尺度数值模拟。根据计算需要,计算域可进行分区划分并建立耦合区域。在耦合区域边界上利用重构算子实现GSM中的宏观变量与LBM中的密度分布函数之间的流动信息传递,而且这两种变量在其对用的方法中均存储在网格节点上。
【机 构】
：
大连理工大学运载工程与力学学部船舶工程学院,大连 116024 大连理工大学运载工程与力学学部船舶
论文部分内容阅读
本文提出了一种梯度光滑法(Gradient Smoothing Method,GSM)与格子Boltzmann的耦合数值方法。基于宏观与介观两种流体计算方法的耦合,该方法能够实现多尺度数值模拟。根据计算需要,计算域可进行分区划分并建立耦合区域。在耦合区域边界上利用重构算子实现GSM中的宏观变量与LBM中的密度分布函数之间的流动信息传递,而且这两种变量在其对用的方法中均存储在网格节点上。
其他文献
对于越来越稀薄的流动,Navier-Stokes方程逐步失效部分归因于其单温度平衡近似.为了克服这个问题,基于Bhatnagar-Gross-Krook(BGK)方程和Landau-Teller-Jeans弛豫模型的连续多温度模型已经被提出(K.Xu et al,2008,并成功应用于二维非平衡流计算.在最近的研究中,已经成功发展出应用于平衡流动的两步四阶气体GKS(L.Pan et al,201
提出一种基于速度修正型浸入边界法(IBM)的PISO求解器，并采用开源CFD软件OpenFOAM对其进行验证.IBM通过欧拉网格中一系列的拉格朗日点来标识物体边界，而物体边界与流场通过回复力的方式相互影响，即在流体运动方程中加体积力项来实现.传统IBM在速度场的计算中事先求解回复力的值，因此其边界点难以完全满足无滑移边界条件.
本文以Hartmann的人工粘性激波捕捉方法为基础，提出了一种基于通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法，该方法形式简单，可向非均匀网格进行推广，编程易于实现。数值试验表明，对于等熵涡算例和高阶曲网格的圆柱算例，在大于二阶的情况下，该方法能有效保证解的精度，对光滑区的影响很小。
本文提出了一种模拟大密度比多相流的界面格子玻尔兹曼通量求解器(ILBFS)。不同于传统的有限体积多相流求解器,是利用宏观物理量来计算单元边界通量,ILBFS利用Chapman-Enskog展开重构宏观控制方程的对流项和扩散项,然后利用格子玻尔兹曼方程的局部解来计算单元边界通量。
多尺度气体流动在航空航天工程中广泛存在,例如,在微喷管或再入飞行器周围,流场尺度的分布跨越稀薄和连续流动区域。基于气体动理学理论的随机粒子方法,如DSMC,可以准确捕捉稀薄气体效应。然而,在分子碰撞占主导的连续流动区域,DSMC方法的计算效率变得非常低。为了克服DSMC方法的这一困难,可以在近连续区域采用模型方法简化分子碰撞过程,如BGK模型。
紧致格式比显式格式具有更好的谱特性，因此在湍流及气动噪声等多尺度流动模拟中备受青睐。为了进一步使紧致格式具有更好的色散特性和更低的耗散，本文提取紧致中心格式中两个自由度用于耗散色散谱特性优化，采用DRP优化思路为基于插值格式的有限差分方法设计了一系列优化紧致格式(L-Opts/7/9)。
在页岩气开采过程中，由于水的存在，流动通道内会出现气水两相同时流动的情况，而对于微纳尺度的流动而言，表面张力会在流动过程中占据主导地位，水的存在会堵塞流动通道，造成“水堵”现象。因此，气水两相界面的流动机理和力学性质对于分析气体在页岩中的流动状态而言至关重要。
本文利用CFD模型对活塞推动式层流涡环形成过程数值模拟,获取高空间分辨率和高时间分辨率的速度场,通过此速度场分别识别流场中的拉格朗日拟序结构。通过对类光测地线的有效计算的自动化识别程序,获取流场的椭圆型LCSs,其椭圆LCSs族的最外层称为拉格朗日相干涡结构,通过示踪识别不同时刻的椭圆LCSs,其在时空上具有相干性质,在一定时间内而不会演化为丝状结构；椭圆型LCSs非定常旋涡边界比基于涡量闭值的边
气体动理学格式(Gas-Kinetic Scheme GKS)是近年发展起来的一种基于介观气体动理论的新型数值方法.GKS利用BGK方程的近似解来构造网格界面上随时间演化的通量，自动耦合粘性和无粘通量，在连续流区的计算量和传统直接基于NS方程的CFD方法相当.由于BGK方程可以描述任意初始分布函数的演化过程，内涵自适应耗散机制，因而在多种流动，尤其是高马赫数粘性流动中能更好地兼顾精度与稳健性.
在多体分离、舱门开启、返回舱降落与螺旋桨滑流等非定常流动模拟中，常采用嵌套网格“挖洞”、找点和插值等技术。在不同网格之间通常采用插值方法解决网格之间的信息传递，实质就是根据背景网格的物理量分布在插值网格上重建物理量。嵌套网格间高精度插值是保证流场计算正确的重要基础。