格子Boltzmann方法多松弛模型局部网格加密算法研究

发布时间：2017-12-15 18:35

  本文关键词：格子Boltzmann方法多松弛模型局部网格加密算法研究

  更多相关文章： 格子Boltzmann方法 多松弛时间模型 网格加密 方腔流 Couette流

【摘要】：格子 Boltzmann 方法(Lattice Boltzmann Method,LBM)是一种新的高效的流体力学计算方法,具有许多独特的优势,已经被广泛应用到流体力学问题研究中。LBM把流体假想成可以沿规则格子移动并仅在格点处相互碰撞的粒子,通过求解碰撞-迁移方程得到一组速度分布函数。压力,速度等宏观参数可以通过求解分布函数的矩方程得到。在流体中,动量、能量、质量的传递往往是不同步的,多松弛(Multiple-Relaxation-Time,MRT)模型在矩空间中有多个松弛参数,把粒子的松弛过程分解成几个独立的松弛过程。通过调整碰撞矩阵的参数将动量,能量和质量的传递区分开来,不仅提高了计算稳定性也提高了计算精度;另外构造线性变换矩阵使得碰撞-迁移后得到的量具有物理意义。在矩空间完成碰撞以后,做逆变换,变换到速度空间,然后流动从一个格点沿离散速度方向迁移到下一个格点。流场中经常有一些区域的物理量变化剧烈,空间和时间梯度较大,特别是在凸角、边缘等处,为了得到精确的物理量的数据,如作用在粒子上的力和力矩等,需要对计算区域的网格进行加密,局部网格加密方法可使物理量变化平稳,提高计算效率。在远离边界和流-固交界面时,流动比较平缓,粗网格可以满足精度要求,可以在流场的大部分区域应用。为此,把计算区域分为粗网格区,过渡区和加密区,不同区域的分布函数通过过渡区的边界进行传递。由于不同区域的网格密度不同,过渡区边界上某些点的数据不能直接由传递得到,这些点称为插值点,为尽量减少误差,使用三次样条插值方法计算插值点上的数据。不同区域的时间步长不同,需要在时间上进行插值,采用三点的Lagrange插值计算。区域界面上的物理量如密度、速度、应力等应保持连续,应用Chapman-Enskog分析推导宏观方程,建立界面上分布函数的转换关系,通过选取合适的松弛参数来提高计算的精度和数值稳定性。边界条件的处理在LBM中占有重要地位,直接关系到计算的精度和效率。流动在特定方向具有周期性时,将具有周期性的区域取出作为计算区域,在边界施加周期性边界条件。对于复杂的曲面边界,采用插值格式对其进行修正;对于一般的平直边界应用反弹格式对边界上碰撞后的分布函数进行修正。Lid-driven cavity flow(方腔流)是检验数值方法的一个经典算例,对雷诺数为1000和400的方腔流进行模拟以验证加密方法的有效性。方腔中位于移动盖板下的左右两个角为奇点,物理量变化剧烈,为得到奇点附近的物理量,对方腔的左上角和右上角进行网格加密。流动初始为静止状态,在经过充分长的时间后达到稳定状态。对比分析不同网格结构得到的速度、压力、涡量、应力,结果表明加密网格处理后的数据更加精确,在奇点附近有非常明显的降噪效果,得到的速度曲线能够更好地与经典算例吻合,降低了压力的振荡幅度,应力曲线的振荡幅度也明显减缓,在应力变化剧烈的区域能够捕捉到应力的变化细节,提高了模拟的质量。对含粒子的Couette流动,应用常规网格、全场加密、局部二重加密和局部三重加密等四种网格结构进行分析计算,程序采用先碰撞后迁移的结构,在两个算例中对计算结果进行了比较。在固定粒子算例中,三种网格结构所得的曳力和举力基本相同,其模拟结果可以作为基准来与运动粒子算例进行比较。在运动粒子算例中,过渡区和加密区与粒子协同运动。使用四种网格结构对比作用在粒子上的曳力,举力和力矩。结果表明三重网格局部加密结构中力和力矩的波动幅度明显小于其他网格结构中的波动幅度,计算准确度得到了提高,并且所需的计算时间少于全场二重加密的计算时间,提高了计算效率。应力的轮廓在加密区与过渡区的交界面都是连续的,证实了局部网格加密方法的正确性。
【学位授予单位】：山东大学
【学位级别】：硕士
【学位授予年份】：2017
【分类号】：O35;TP309.7

【相似文献】

中国期刊全文数据库 前10条

1 ;Lattice Boltzmann simulation for the spiral waves in the excitable medium[J];Communications in Nonlinear Science & Numerical Simulation;2000年04期

2 孙成海,王保国,沈孟育;Adaptive Lattice Boltzmann Model for Compressible Flows[J];Tsinghua Science and Technology;2000年01期

3 张英魁 ,卢旭光;Boltzmann Equations with Quantum Effects (1):Long Time Behavior of Spatial Decay Solutions[J];Tsinghua Science and Technology;2002年03期

4 张英魁 ,卢旭光;Boltzmann Equations with Quantum Effects (2): Entropy Identity, Existence and Uniqueness of Spatial Decay Solutions[J];Tsinghua Science and Technology;2002年03期

5 ;New Boundary Treatment Methods for Lattice Boltzmann Method[J];Wuhan University Journal of Natural Sciences;2003年01期

6 ;Gas kinetic algorithm for flows in Poiseuille-like microchannels using Boltzmann model equation[J];Science in China(Series G:Physics,Mechanics & Astronomy);2005年04期

7 田红晓;刘全生;阿其拉图;;二粒子Boltzmann方程组的初始层解[J];内蒙古大学学报(自然科学版);2013年04期

8 黄祖洽,丁鄂江;Boltzmann方程的理论和应用[J];物理学进展;1986年03期

9 刘清荣;;线性Boltzmann方程解的存在性和界[J];纯粹数学与应用数学;1987年00期

10 陈建宁;Boltzmann方程的两个解的平均作为二粒子Boltzmann方程系的解[J];数学物理学报;1990年03期

中国重要会议论文全文数据库 前10条

1 Yu CHEN;Shulong TENG;Hirotada OHASHI;;On the Lattice Boltzmann Modeling of Multi-Phase Flows[A];Genetic Algorithm and Its Application to Physics, Life Science and Engineering--Proceedings of CCAST (World Laboratory) Workshop[C];1999年

2 冯士德;;多种粒子相耦合的格子Boltzmann模型[A];“力学2000”学术大会论文集[C];2000年

3 冯士德;;多组份流体的格子模型Boltzmann[A];第十四届全国水动力学研讨会文集[C];2000年

4 梁宏;施保昌;;基于相场理论的不可压格子Boltzmann两相流模型及其应用[A];第十六届全国流体力学数值方法研讨会2013论文集[C];2013年

5 丁丽霞;施卫平;郑海成;;格子Boltzmann方法模拟圆腔流[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集（下）[C];2007年

6 李志辉;方明;唐少强;;气体运动论统一算法与DSMC方法描述Boltzmann方程的一致性研究[A];北京力学会第18届学术年会论文集[C];2012年

7 谭玲燕;李秀文;李伟杰;;用Lattice Boltzmann方法模拟腔体内搅动流体[A];第十六届全国流体力学数值方法研讨会2013论文集[C];2013年

8 李志辉;彭傲平;张涵信;;求解Boltzmann模型方程的高性能并行计算方法研究[A];中国力学大会——2013论文摘要集[C];2013年

9 冯士德;毛江玉;任荣彩;;格子Boltzmann数值流体力学[A];钱学森技术科学思想与力学论文集[C];2001年

10 彭傲平;李志辉;吴俊林;;求解Boltzmann模型方程的气体运动论隐式格式及应用研究[A];中国力学大会——2013论文摘要集[C];2013年

中国博士学位论文全文数据库 前10条

1 李元;格子Boltzmann方法的应用研究[D];中国科学技术大学;2009年

2 胡安杰;多相流动格子Boltzmann方法研究[D];重庆大学;2015年

3 刘克同;新型格子Boltzmann流固耦合算法研究及其在桥梁风工程中的应用[D];哈尔滨工业大学;2015年

4 毛煜东;微纳尺度传热问题的理论分析和格子Boltzmann数值模拟[D];山东大学;2015年

5 刘智翔;改进的格子Boltzmann方法研究及大规模并行计算[D];上海大学;2014年

6 谢海琼;基于格子Boltzmann方法模拟热毛细对流[D];重庆大学;2015年

7 尹显利;弹性波方程的格子Boltzmann方法研究[D];吉林大学;2016年

8 罗忠贤;联合格子Boltzmann方法和X-ray CT成像技术研究开级配沥青混合料中水的输运规律[D];西南交通大学;2015年

9 樊迎哲;扰动框架下Vlasov-Maxwell-Boltzmann方程组的整体适定性及相关问题的研究[D];武汉大学;2017年

10 李钱焕;耦合对流扩散方程的格子Boltzmann模型及应用研究[D];华中科技大学;2016年

中国硕士学位论文全文数据库 前10条

1 叶丽娜;基于格子Boltzmann模型的Kuramoto-Sivashinsky方程的数值模拟研究[D];吉林大学;2011年

2 郭宇隆;基于格子Boltzmann方法的气液混合流体模拟[D];华南理工大学;2015年

3 伍祥超;格子Boltzmann方法在流体流动及传热中的应用研究[D];昆明理工大学;2015年

4 郭槛菲;格子Boltzmann方法在气固两相流中的应用研究[D];山东大学;2015年

5 韩文骥;适合飞溅流动的格子Boltzmann方法数值研究[D];中国舰船研究院;2015年

6 周训;单气泡沿过热曲面运动的格子Boltzmann方法模拟[D];大连理工大学;2015年

7 薛福;基于格子Boltzmann方法的饱和土壤渗流与传热数值研究[D];大连理工大学;2015年

8 罗丰;基于Boltzmann方程的一致气体动力学格式对连续稀薄气体流动的数值模拟的研究[D];湘潭大学;2015年

9 冯笑含;基于Boltzmann的煤层微孔道气体流动规律研究[D];东北石油大学;2015年

10 王朝冰;基于格子Boltzmann方法的毛细管节流特性研究[D];大连海事大学;2016年

，

本文编号：1293016