用格子Boltzmann方法模拟方柱绕流

基于格子Boltzmann方法的原理和对边界条件处理的便捷性,分别对单方柱和并列双方柱绕流进行了研究.对于单方柱绕流,分别研究了不同雷诺数和阻塞比对它的影响,给出了不同雷诺数下单方柱绕流的流线图和涡量图,以及平均阻力系数和最大升力系数随雷诺数的变化曲线,得到单方柱产生卡门涡街的临界雷诺数,也给出了不同阻塞比下单方柱绕流...     

平板间方柱绕流的格子Boltzmann方法模拟

为研究平板间方柱绕流上下平板对置于其中的方柱绕流所产生的影响,采用格子Boltzmann方法对二维平板间低雷诺数(Re=100)方柱绕流问题进行了数值模拟研究.分析了3种不同阻塞比下,平板边壁对方柱的升、阻力系数、Strouhal数和尾涡流场的影响.结果表明:平板对方柱绕流特性有明显的影响,随着阻塞比的增加,阻力系数和Strouhal数均增大,与无边壁相比阻力系数可增加达30%,而升力系数却随之减小.计算结果与相关实验数据相吻合,验证了格子Boltzmann方法对钝体绕流非定常问题模拟的有效性。     

基于动态Smagorinsky-MRT-LBM的圆柱绕流湍流研究

采用多松弛格子玻尔兹曼与大涡模拟的动态Smagorinsky亚格子涡粘模型相结合的方法对圆柱绕流湍流进行模拟,研究300至10 000不同雷诺数下的单圆柱绕流湍流涡的周期性变化。研究结果表明:计算所得的斯考特数值和阻力系数值与文献中实验结果以及有限元仿真结果基本一致,说明该方法对于不同雷诺数的湍流模拟是合理的。     

来流湍流度对圆柱绕流特性的影响

采用大涡模拟方法,在雷诺数为1.4?05时,研究来流湍流度对圆柱绕流特性的影响规律。研究表明：湍流度会使得自由剪切层间的相互作用更加剧烈,使得流场的三维特性更为强烈;同时加强了尾流区湍流能量的相互传递,减少能量的耗散,从而延长了回流区长度。随着回流长度的增加,圆柱背风面的平均风压系数得以增大,最终导致圆柱阻力系数随之减小。在尾流湍流得到增强的同时,尾流区会形成多尺度的湍流,从而促使漩涡形成区形成大涡量的漩涡,导致尾流区的漩涡呈现多频率的脱落。     

用格子Boltzmann方法模拟圆柱绕流的振荡边界层

提出了一种用于模拟圆柱绕流的振荡边界层流动的新方法——旋转系统平均法。该系综由一些相同的系统组成,每个系统除使用不同的网格外,具有相同的流动。经模拟圆柱绕流的振荡边界层发现,流体质点的流动具有平行于振荡方向离开圆柱、竖直于振荡方向流向圆柱的特点,同时也发现在紧邻圆柱出现的被称为二次流的4个涡旋流动现象。数值模拟结果符合经典的理论结果。该方法也可用于其他圆柱绕流问题。     

用格子Boltzmann方法模拟圆柱绕流的振荡边界层

提出了一种用于模拟圆柱绕流的振荡边界层流动的新方法--旋转系统平均法.该系综由一些相同的系统组成,每个系统除使用不同的网格外,具有相同的流动.经模拟圆柱绕流的振荡边界层发现,流体质点的流动具有平行于振荡方向离开圆柱、竖直于振荡方向流向圆柱的特点,同时也发现在紧邻圆柱出现的被称为二次流的4个涡旋流动现象.数值模拟结果符合经典的理论结果.该方法也可用于其他圆柱绕流问题.     

基于动态亚格子模型的方柱绕流大涡模拟

目的 对高雷诺数下的方柱绕流进行数值模拟,得到方柱下流场时间平均结果以及湍流流场的流场特征,用于验证亚格子模型对大涡模拟(LES)计算精确性的影响.方法 基于有限体积法(FVM)的流体计算开源软件OpenFOAM,采用两种不同动态亚格子模型下的大涡模拟方法,模拟高雷诺数(Re =2.2 ×104)下的方柱绕流.结果 通过数值计算得到流场完整的平均和脉动结构,并得到阻力系数、升力系数、斯托拉哈数等重要的气动力参数,经对比与试验实际情况吻合.结论 高雷诺数复杂流场下亚格子应力的封闭形式对大涡模拟的计算结果有重要的影响,采用动态模拟技术一方程模型计算结果优于常规的Smagorinsky模型,在高雷诺数复杂流场下的计算是高效的,为结构抗风的数值研究提供了参考.     

用多松弛格子Boltzmann方法模拟三维水击波

旨在探索用多松弛格子Boltzmann方法模拟多维水击波的可行性.采用三维方法模拟了一维水电站简单管水击现象,得到准确的直接水击和间接水击过程,证明其能准确模拟压力波传播与反射特性;模拟了二维拐角衍射,得到与其他方法吻合的结果;详细分析了断波撞击三维台阶产生的反射和衍射现象,获得合理的压力分布.研究表明:多松弛格子Boltzmann方法稳定性好,精度较高,能够准确模拟水击波的传播、反射和衍射,可推广到实际问题研究中.     

圆柱绕流涡致振动的平面湍流数值模拟

为了研究工程中存在的湍流涡致振动问题,应用平面湍流k-ε模型和重整化群k-ε模型,结合分块耦合方法及任意拉格朗日欧拉法（ALE）数值模拟了圆柱绕流涡致振动,求解基于非交错网格系统, 其中包含满足连续性方程的压力泊松方程解法.计算了湍动能、湍流耗散率、湍流黏性系数及柱面涡量的变化情况.结果表明,当雷诺数为5 000～2×105时计算所得的阻力系数与实验数据吻合良好.研究发现,在涡致振动情况下的湍流耗散率比固定圆柱涡脱落情况下小得多,湍流黏性系数与雷诺数呈对数律增长.     

粘性流体圆柱绕流的有限元模拟

对粘性流体圆柱绕流进行了数值模拟。根据三阶Taylor展开弱解的概念,在二阶分裂步Taylor-Galerkin有限元法基础上,采用时间推进和张量分析的方法推导了对流扩散通用微分方程的ETG有限元离散格式。对Reynolds数为32、102和250的粘性流体圆柱红流流场进行了数值仿真,数值模拟结果与实验结果吻合良好。表明隐含流线迎风耗散作用并具有三阶部分展开的Taylor-Galerkin有限元法稳定性好,精度较高,可用于计算包含对流扩散耦合传递动量的绕流问题及类似的流动问题。     

通用的插值补充格子Boltzmann方法应用于计算气动声学

提出将通用的插值补充格子Boltzmann方法（GILBM）应用于非均匀网格进行计算气动声学研究. 通过顶盖驱动方腔流、低雷诺数圆柱绕流算例验证GILBM数值模拟方法的正确性. 在此基础上,将该方法应用于高斯脉冲传播、周期性声源传播、二维圆柱绕流的气动噪声计算. 研究结果表明,GILBM方法可以在非均匀网格上较好地模拟高斯脉冲及周期性声源声波的传播过程,计算结果与解析解较吻合. GILBM方法能够模拟非均匀贴体网格下的圆柱曲面边界由于涡脱落造成的气动噪声的产生和传播过程,可以较好地捕捉到近场及远场的声压传播. 圆柱绕流声学特性呈现出偶极子现象,计算结果与参考文献较吻合,证明采用GILBM方法在非均匀网格中模拟声传播问题的正确性及求解气动声学问题的可行性.     

Numerical simulation of flow around square cylinder using different low-Reynolds number turbulence models

ABE-KONDOH-NAGANO, ABID, YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field. Meanwhile, drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed. The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results, and they are significantly better than those of ABE-KONDOH-NAGANO, ABID and LAUNDER-SHARMR models. The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model, while the correction of the turbulence near the wall has minor influence on the simulation results. For ABE-KONDOH-NAGANO, ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy. With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.     

基于插值的Lattice Boltzmann方法非均匀网格算法

提出了一种基于插值的latticeBoltzmann方法的非均匀网格算法 ,给出了计算步骤 ;对典型的二维恒定方腔流、恒定和非恒定圆柱绕流问题进行了模拟分析 ,所得结果的正确性证明了该算法是有效、可行和实用的     

多孔介质内流体流动的大涡格子Boltzmann方法研究

为研究多孔介质内流动随Re数变化的特点,采用结合Smagorinsky亚格子模型的格子Boltzmann方法(LES-LBM)对多孔介质内流动进行了数值模拟.结果表明：多孔介质内的单相流动在高Re数时会表现出复杂的非线性现象;LES-LBM克服了传统LBGK方法模拟高Re数流动时容易产生数值不稳定的缺点,能清晰地描述出多孔介质内流动存在的3个区域,即低速时的线性达西区、过渡区和高速时的非线性二次区;不同Re下的流线图还说明微观的惯性作用最终导致了多孔介质宏观上的非线性现象,多孔介质流动呈现明显的多尺度特征.进一步分析计算结果可以证明：LES-LBM方法能准确地验证Darcy-Forchhimer阻力方程,Darcy-Forchhimer总阻力随Re数增加而增加,随孔隙率增加而减小,并且小孔隙率下的Forchhimer阻力占总阻力比例小于大孔隙率时的比例.     

血液在分叉管栓塞的晶格玻尔兹曼方法研究

基于对血栓形成机理问题,建立了血液在分叉血管中流动的晶格玻尔兹曼模型,研究在不同压差、压积下血栓的形成机理。模型用刚性小球代替红细胞,为了使模拟更符合实际,在管壁与小球之间增加了摩擦力。结果表明,压积越大越容易发生堵塞,压差越大越不容易发生堵塞,在入口以及分叉处最易发生堵塞。     

锥形管中血液栓塞的晶格玻尔兹曼方法

针对血液循环中血液栓塞问题,基于晶格玻尔兹曼方法建立二维锥形管模型,以不可变形的多个刚性粒子代表红细胞模拟血液在锥形管中的栓塞,考察粒子在流场中的平均速率,对粒子在管道中的运动进行跟踪,统计粒子在管道中的位置,求出粒子在管道不同水平位置的容积率,并记录栓塞瞬间整个管道粒子的分布情况。研究结果表明,压积越大越易发生栓塞,压差越大越不易发生栓塞,锥形管的出口处容易发生栓塞。