模拟MKDV方程的格子Boltzmann方法

目前,格子Boltzmann方法已被广泛应用于模拟各种非线性物理方程。文中用D1Q2模型给出MKDV方程的带修正项的BGK型格子Boltzmann方法。通过对演化方程的泰勒展开并应用多尺度技术恢复了宏观方程。数值模拟表明文中所建立的模型是有效的。     

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

基于二维紧凑型四阶格子模型的激波管流动模拟

为了提高格子的稳定性，使用Hermite展开方法，构建了新的二维四阶紧凑型格子模型，即D2Q37A。比较了D2Q37A和与Philippi给出的紧凑型格子模型(D2Q37B)的稳定性。在相同的碰撞频率下，与D2Q37B相比，D2Q37A可以模拟初始密度比更高的一维激波管流动。这表明D2Q37A与现有格子模型相比，具有更好的稳定性。详细给出了适用于高阶格子模型的边界条件实现方式。此边界条件实现方式保留了体现LBM（lattice Boltmann method）粒子特性的迁移 碰撞机制。用以上给出的格子模型和边界条件处理方式模拟激波管流动，得到的模拟结果和解析解吻合得很好。这表明所给出的边界处理方式是可行的。此边界格式同样可以用于其他类型的流动和边界。      

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

广义守恒相场简化多相流格子Boltzmann方法

针对不可压缩、非混溶的复杂多相流问题，提出一种广义守恒相场简化多相流格子Boltzmann方法。此方法运用早前发展的简化多相流格子Boltzmann方法（simplified multiphase lattice Boltzmann method,SMLBM），通过采用带有拉格朗日算子的广义守恒相场方程来控制界面的演化并确保每个相的体积和总质量守恒。此外，在单松弛格子Boltzmann方法框架内，SMLBM是通过预测-校正策略来模拟流体系统和跟踪界面演化，其计算过程中仅需要考虑平衡态分布函数的演化，并且平衡态分布函数可直接从宏观量计算得到，因而具有良好的稳定性、高计算效率和边界条件易于实施的优点。本方法继承了SMLBM的优势，能够解决由不同流体组分之间的大密度比和大黏度比引起的界面处大压力梯度问题。为了验证本方法的稳定性和准确性，模拟了包括拉普拉斯定律、液滴透镜、三相泊肃叶流以及复合液滴铺展在内的四个多相流算例。结果表明，本方法能有效地模拟密度比达1 200和黏度比达500的复杂界面算例。     

高超声速湍流高效模拟算法

采用了直接过滤的Navier-Stokes(N-S)方程组对高速可压缩湍流进行研究。针对高超声速湍流的非线性流动特性,对N-S方程直接过滤推导了大尺度湍流流场的控制方程,更精确地反映高速湍流的可压缩性,建立了可压缩湍流的大涡模拟TDM模型。使用传统的Smagorinsky模型的非线性推广,采用基于非Favre过滤的超声速可压缩湍流的大涡模拟模型,应用Caylay-Hamilton定理,建立可压缩湍流大涡模拟的非线性亚格子模型,并发展为动力学模式,模型中的两个常数通过当地流场动态的确定,消除了可调经验常数的影响。针对构造的高超声速湍流大涡模型开发相应的高效并行算法。     

关于格子Boltzmann方法的几个问题(Ⅰ)

格子气方法及后续发展的格子Boltzmann方法,是一种完全不同于传统流场计算的基于介观物理的新方法.由于其内在的优势,尽管该方法问世不过几十年的时间,就在传统计算流体力学难于处理的研究领域得到了广泛的应用,被誉为是现代流体力学的一场变革.文章主要论述了格子Boltzmann方法有关学术思想的源头,论述在该论题发展过程中的一些关键性工作.     

烧蚀外形转捩流动的数值模拟

本文采用雷诺平均Navier-Stokes方程和修正Baldwin-Lomax代数湍流模型,以Roe二阶流通量差分分裂格式进行离散,数值模拟了典型烧蚀外形的高超声速流场,重点讨论了流动转捩的判定及其与复杂流态的关联.结果表明:数值模拟出的球钝锥外形转捩点位置与工程估算结果符合较好;对工程估算中存在困难的奶头状凹陷烧蚀外形,数值模拟也可给出符合流场结构的转捩点位置.     

有限体积格子Boltzmann方法用于近空间连续流区绕流模拟

针对航天器再入解体形成残骸碎片的近空间绕流计算问题,拓展了格子Boltzmann方法在可压缩流动模拟的能力。引入有限体积隐式格式求解耦合双分布函数格子Boltzmann模型方程及采用圆函数为基础构造的D2Q13离散速度模型;引入IMEX-RK格式进行了时间项离散解决模型方程的源项刚性问题;对Riemann问题、平板双马赫反射问题、RAE2822翼型跨声速绕流等近空间连续流区、可压缩典型案例进行了数值模拟。通过比较分析,初步验证了耦合双分布函数有限体积格子Boltzmann方法对连续流区、可压缩流动的模拟能力;进一步开展了方柱形解体残骸的超声速绕流模拟,得到与N-S方程计算结果吻合一致的绕流结果,证实经改进的耦合双分布函数有限体积格子Boltzmann方法对解体残骸碎片绕流问题具有较好的模拟能力。     

黏性彻体力建模方法及其应用

研究了将Navier-Stokes方程中黏性项简化为黏性彻体力的建模方法,并将该模型加入已有的计算流体力学(CFD)程序中,实现了用三维欧拉方程加黏性彻体力模型快速模拟叶轮机三维非定常流动的方法.使用该方法对进气周向总压畸变条件下的某双级压气机非定常流场进行了模拟,并与雷诺平均的Navier-Stokes(RANS)程...     

用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

结合响应面法与网格法的低速翼型优化设计

本文首先了介绍响应面法和网格法的基本原理,翼型绕流场分析使用有限体积法空间离散和RUNGE-KUTTA时间迭代格式进行N-S方程计算,由翼型表面压力和摩擦力数值积分得到升力和阻力.然后,以NACA0012和GAw-1为算例,结合响应面法和网格法进行低速翼型的减阻优化设计,算例表明:在应用响应面法进行低速翼型优化时,采用...     

轴流压气机叶尖射流扩稳数值计算模型

基于体积力方法将压气机叶片对气流的作用力和做功简化为源项,对环形无叶片通道求解带源项的Navier-Stokes(N-S)方程,实现了压气机全通道流场的计算,通过在转子上游机匣施加射流建立了轴流压气机叶尖射流扩稳数值模型.模型中的源项根据叶片排进口气流参数进行变化,通过判断施加的人工扰动变化确定压气机失稳边界.采用该模型对某跨声速单级压气机在有、无射流下的性能和稳定性进行了计算分析,无射流下计算的该压气机性能和稳定边界与试验结果都能较好地一致,其中稳定边界流量仅相差0.7%;有射流下射流对压气机具有明显的扩稳效果,将计算结果与国外公开文献的结果进行了对比,验证了模型的可靠性.      

Prediction of Boundary Layer Transition Based on Modeling of Laminar Fluctuations Using RANS Approach

This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transport equations, separately, to compute laminar kinetic energy, turbulent kinetic energy, and dissipation rate in a flow field. It needs neither correlations of intermittency factors nor knowledge of the transition onset. Two transition tests are carried out: flat plate boundary layer under zero ...     

应用等离子体进行翼型前缘分离控制的数值模拟研究

将等离子体对中性气体的作用模型化为彻体力矢量,求解带源项的Navier-Stokes方程,数值模拟了等离子体激励器在NACA0015翼型大迎角下的分离控制效果,彻体力为净电荷在外加电场作用下的电场力.解拉普拉斯方程得到外加电场分布,等离子体中的净电荷分布由泊松方程给出.为了较好地模拟分离涡的发展,采用了雷诺平均与大涡模拟相结合的脱体涡模拟(Detached Eddy Simulation)方法.通过与实验结果对比,发现该模型能较好地模拟等离子体激励器的控制作用.     

应用等离子体进行翼型升力控制的数值模拟研究

将等离子体对中性气体的作用模型化为彻体力矢量,求解带源项的Navier-Stokes方程,数值模拟了在NLF(1)-0213翼型上表面60%弦长处安装等离子体激励器对升力的控制效果.彻体力为净电荷在外加电场作用下的电场力.解拉普拉斯方程得到外加电场分布,等离子体中的净电荷分布由泊松方程给出.升力线计算结果与实验值吻合,激励器工作时,升力线向上平移,控制效果与襟翼类似.     

格构式塔架三维动力风荷载的风洞试验研究

通过底座天平测力风洞试验,测量了三种典型的格构式塔架(输电塔、通讯塔和电视塔)的顺风向、横风向和扭转向动力风荷载。由于模型本身的质量不够轻、刚度不够大,一阶共振力的影响不可忽视,本文作者推导了消除模型一阶共振影响的实用公式,进而得到格构式塔架顺、横、扭三维广义荷载谱和均方根广义力系数。最后,本文分析了这三种格构式塔架的动力风荷载特征,得到了几点重要的结论。