收藏本站
《吉林大学》 2011年
收藏 | 手机打开
二维码
手机客户端打开本文

格子Boltzmann模拟动边界流动问题

郑海成  
【摘要】:格子Boltzmann方法(Lattice Boltzmann Method, LBM)是一种新兴的介观数值方法,在1988年由McNamara等首先提出。LBM起源于格子气自动机(Lattice Gas Automata, LGA),用实数f代替布尔量来表示系统平均后的局部粒子分布函数,以及用Boltzmann方程代替格子气自动机的演化方程,是LBM诞生的标志。LBM的粒子分布函数只与粒子相邻的节点信息相关,因此适合并行运算。粒子分布函数的演化方程是线性的,因此容易编程。LBM的边界处理也很简单,适合一些复杂边界问题。 本文采用LBM二维最常用的D2Q9模型,相关的演化方程为其中ρ是宏观密度,ωα是加权系数, c s是格子声速,D2Q9模型中具体取值如下宏观压强P =ρcs2,宏观密度和速度由下面的计算式决定为了计算方便,一个循环过程分为碰撞和传递两个过程,对应的方程形式为碰撞过程: 由表达式可以看出,碰撞过程只发生在节点上,因此又被称作局部碰撞。迁移过程是指碰撞后的粒子分布函数按照速度方向传递到相邻的节点上,在固体边界处,一般是将通过各种方法得到的虚拟分布函数传递到流体节点上。本文采用的处理方法是Lallemand等[100]提出的格式。 从LBM微观演化方程出发,应用Chapman-Enskog展开,可以恢复到宏观控制方程。时间偏导数和分布函数的Chapman-Enskog展开形式为将LBM演化方程左端第一项Taylor展开,按照ε各阶前的系数列出方程组,遵循质量和动量守恒,结合一些合理的假设,可以恢复到宏观控制方程。 在数值模拟中,初始条件和边界条件是影响运算稳定性和收敛性的重要因素。LBM数值模拟中,边界条件处理尤为重要。一些简单的边界处理格式包括反弹格式和镜面反射格式、周期性边界处理格式、充分发展边界处理格式;已知边界处物理条件的动力学边界处理格式;将边界处分布函数分成平衡态和非平衡态两部分,非平衡态部分的演化形式与平衡态的部分相似;另外还有复杂边界的虚拟平衡态分布函数方法,更是包含了运动边界的影响。 数值模拟的第一部分是并列旋转双圆柱绕流问题,其工程背景是多钝体绕流。作为一种消涡减阻方法,两个并列圆柱的逆向旋转可以有效地抑制旋涡的脱落。旋涡的消减效果与旋转速度和两个圆柱间的间距有关,临界旋转速度是这类问题的研究重点。本文选取的旋转速度比为0.5、1.0、1.5和2.0,间距比的范围是1.2~3.0。通过对数值结果的分析,得到以下几点结论: 1.在小旋转速度比(静止圆柱除外)下,小间距比的流型一般为单钝体流型,大间距比的流型一般为反相对称流型。 2.在较大旋转速度比下,小间距比的流型为单钝体流型或者旋涡脱落完全受到抑制,大间距比的流型为反相对称流型。 3.两个圆柱间流体的流量与间距比和旋转速度比都有关,它随着间距比增加而增加,随着旋转速度比增加而减少。 4.流型与两个圆柱间流体流量有关,流量大时容易形成旋涡脱落,流型为偏流流型或对称流型,流量小时一般为单钝体流型。 数值模拟的第二部分是运动平板附近圆柱绕流问题,其工程背景是赛车等问题。雷诺数Re和间隙率G/ D是两个重要参数,判定圆柱下方旋涡完全受到抑制时的临界间隙率是这类问题的重点。本文计算的雷诺数Re = 600,间隙率0 .1≤G /D≤1.0。间隙率不同时,间隙处的流向平均速度分布、圆柱后的旋涡脱落形态、升阻力系数、斯特劳哈尔数等都会表现出不同的特性。本文分别对这些流场特性作了详尽的分析,得到一些有用的结论,总结出一般性的判定临界间隙率的几种方法: 1.通过间隙处平均流向速度的剖面寻找。一般地,间隙率小于临界值时,速度剖面只有一个极大值;间隙率大于临界值时,速度剖面会有两个极大值;在临界间隙率时,速度剖面的两个极大值相等。因此,可以通过平均流向速度得到临界间隙率。 2.根据旋涡脱落形态判定临界间隙率。间隙率小于临界值时,圆柱只有远离平板的一侧有旋涡脱落;间隙率大于临界值时,圆柱两侧都有旋涡脱落;临界间隙率时,圆柱靠近平板一侧刚刚出现旋涡脱落现象。因此,可以通过观察圆柱靠近平板一侧是否有旋涡脱落现象来判定临界间隙率。 3.通过阻力系数判定临界间隙率。阻力系数随时间变化曲线中,间隙率大于临界值时,阻力系数有两个极大值;间隙率小于临界值时,阻力系数只有一个极大值;临界间隙率时,阻力系数曲线中会出现一小段平台。另外,间隙率达到临界值时,平均阻力系数最大。因此,阻力系数也可以作为一种判定临界间隙率的方法。 4.通过斯特劳哈尔数判定间隙率。随着间隙率由小到大增加,斯特劳哈尔数也逐渐增加。在间隙率达到临界值时,斯特劳哈尔数也进入平缓阶段,临界值是一个转折点。 数值模拟的第三部分是周期性变形圆柱的绕流问题,这可以看作是生物游动时自适应现象的初步研究。研究这个问题的目的是看圆柱变形是否能够达到消涡减阻的目的。圆柱变形主要是指变形幅度和变形频率,本文选取变形幅度相对圆柱半径的变形幅度比为0.05、0.075、0.1和0.15,变形频率相对孤立圆柱旋涡自然脱落频率的变形频率比的范围是0.8-2.0。通过对数值结果的分析,得到如下规律: 1.较小的变形频率比和变形幅度比对阻力的影响都不大; 2.当变形频率比较小时,变形幅度比值越大,阻力越小;当变形频率比较大时,变形幅度比值越大,阻力越大,并且阻力随着变形频率比增加而增加; 3.小变形频率比和变形幅度比下,旋涡脱落模态以2S为主;大变形幅度下,旋涡脱落以P+S和2P为主; 4.当ar =0.15、fr =1.2时,变形圆柱所受阻力小于孤立圆柱所受的阻力,圆柱的周期性变形确实有消涡减阻的作用。
【学位授予单位】:吉林大学
【学位级别】:博士
【学位授予年份】:2011
【分类号】:O357.1

手机知网App
【引证文献】
中国重要会议论文全文数据库 前1条
1 郑海成;拾兵;张靖;;带有自由端的可变形膜绕流数值研究[A];第十三届全国水动力学学术会议暨第二十六届全国水动力学研讨会论文集——B水动力学基础[C];2014年
中国博士学位论文全文数据库 前1条
1 付博;基于格子Boltzmann方法的Rayleigh对流及其对界面传质影响的研究[D];天津大学;2012年
【参考文献】
中国期刊全文数据库 前10条
1 施鎏鎏;余俊;万津津;刘应征;;近壁方柱绕流非定常特性的TR-PIV测量[J];工程热物理学报;2009年08期
2 沈智军,袁光伟,沈隆钧;格子Boltzmann方法求解Burgers方程(英文)[J];计算物理;2000年Z1期
3 刘曰武,周富信,闫广武;Lattice Boltzmann方法模拟多孔介质Klinkenberg效应(英文)[J];计算物理;2003年02期
4 聂德明,林建忠;格子Boltzmann方法中的边界条件[J];计算物理;2004年01期
5 杜睿;施保昌;;格子Boltzmann方法中的曲边界处理[J];计算物理;2006年04期
6 施卫平,胡守信,阎广武;用改进的格子Boltzmann方法模拟浅水长波问题[J];计算物理;1997年Z1期
7 熊盛武,李元香,康立山,陈炬桦,阮剑;热动力学格子Boltzmann模型[J];计算物理;1998年04期
8 贾复 ,魏启鲲;流体力学年鉴1985年第17卷[J];力学进展;1985年04期
9 邵传平;;钝体尾流控制机理及方法研究进展[J];力学进展;2008年03期
10 祖迎庆,施卫平;用格子Boltzmann方法模拟流场中可变形膜的运动[J];力学学报;2005年02期
中国博士学位论文全文数据库 前1条
1 丁丽霞;用于模拟粘性流体流动的格子Boltzmann方法[D];吉林大学;2009年
【共引文献】
中国期刊全文数据库 前10条
1 张德兴;不可压缩流的晶格玻耳兹曼模型的讨论[J];安徽师范大学学报(自然科学版);2003年04期
2 闫广武;;A LAGRANGIAN LATTICE BOLTZMANN METHOD FOR EULER EQUATIONS[J];Acta Mechanica Sinica;1998年02期
3 ;A novel incompressible finite-difference lattice Boltzmann equation for particle-laden flow[J];Acta Mechanica Sinica;2005年06期
4 伊厚会;李红梅;;颗粒在流体中相互作用的晶格玻尔兹曼方法模拟[J];滨州学院学报;2010年06期
5 ;Highly Efficient Lattice Boltzmann Model for Compressible Fluids:Two-Dimensional Case[J];Communications in Theoretical Physics;2009年10期
6 陈锋;许爱国;张广财;李英骏;;Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows[J];Communications in Theoretical Physics;2010年12期
7 常茹;;Numerical simulation of wind flow fields around buildings using the lattice Boltzmann method[J];Journal of Chongqing University(English Edition);2011年02期
8 李维仲;董波;宋永臣;;水平通道内黏性指进现象LBM模拟[J];大连理工大学学报;2012年03期
9 朱益华;陶果;方伟;王胜奎;;低渗气藏中气体渗流Klinkenberg效应研究进展[J];地球物理学进展;2007年05期
10 田智威;;煤层气渗流理论及其研究进展[J];地质科技情报;2010年01期
中国重要会议论文全文数据库 前10条
1 陈瑜;蔡庆东;;格子玻尔兹曼多相流模拟中无滑移边条件选取[A];北京力学会第18届学术年会论文集[C];2012年
2 毕继红;任洪鹏;丁代伟;于会超;;串列双圆柱静止绕流的二维数值仿真分析[A];第20届全国结构工程学术会议论文集(第Ⅰ册)[C];2011年
3 刘天成;刘高;葛耀君;吴宏波;曹丰产;;基于拓展LB方法的桥梁结构数值风洞软件研发及应用[A];第十四届全国结构风工程学术会议论文集(下册)[C];2009年
4 王兴勇;索丽生;程永光;;用Lattice Boltzmann方法模拟二维管道连接处的流场[A];第十六届全国水动力学研讨会文集[C];2002年
5 王嘉松;;附加分离盘控制隔水管涡激振动的高分辨率数值模拟研究[A];第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会论文集[C];2009年
6 王兴勇;郭军;刘树坤;谢省宗;;Lattice Boltzmann方法的分块-耦合算法[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
7 刘建军;裴桂红;薛强;;微流边界层对孔隙介质渗流的影响[A];第一届中国水利水电岩土力学与工程学术讨论会论文集(上册)[C];2006年
8 张华;;液滴冲击液面的仿真研究[A];第三届全国水力学与水利信息学大会论文集[C];2007年
9 李维仲;张晓红;董波;孙涛;;单汽泡沿倾斜表面运动特性研究[A];高等学校工程热物理第十九届全国学术会议论文集[C];2013年
10 郑海成;拾兵;张靖;;带有自由端的可变形膜绕流数值研究[A];第十三届全国水动力学学术会议暨第二十六届全国水动力学研讨会论文集——B水动力学基础[C];2014年
中国博士学位论文全文数据库 前10条
1 吴锋;基于GPU并行计算的数值模拟与燃煤锅炉系统的优化研究[D];浙江大学;2010年
2 张邵波;纳米流体强化传热的实验和数值模拟研究[D];浙江大学;2009年
3 闫铂;反应扩散方程的格子Boltzmann方法及数值模拟[D];吉林大学;2011年
4 张建影;用于复Ginzburg-Landau方程的格子Boltzmann方法[D];吉林大学;2011年
5 刘芳;格子Boltzmann方法求解偏微分方程的相关研究[D];吉林大学;2011年
6 荣伏梅;基于格子Boltzmann方法的多孔介质REV尺度的流动与强化传热研究[D];华中科技大学;2011年
7 周陆军;磁流体流动及能量传递特性的多尺度研究[D];南京理工大学;2010年
8 赵凯;基于孔隙尺度的多孔介质流动与传热机理研究[D];南京理工大学;2010年
9 栗怀广;钝体空气动力学的格子玻尔兹曼方法[D];西南交通大学;2010年
10 王勇;格子Boltzmann方法在热声领域的应用及热声谐振管可视化实验研究[D];西安交通大学;2009年
中国硕士学位论文全文数据库 前10条
1 陈海明;深水立管束涡激振动分析和数值模拟[D];中国海洋大学;2010年
2 李士伟;类客车体气动特性及减阻措施的试验研究和数值模拟[D];中国空气动力研究与发展中心;2010年
3 陈林婕;几类非线性偏微分方程的格子Boltzmann方法[D];福建师范大学;2010年
4 叶丽娜;基于格子Boltzmann模型的Kuramoto-Sivashinsky方程的数值模拟研究[D];吉林大学;2011年
5 党亚兵;基于格子Boltzmann的等离子射流数值模拟[D];河南理工大学;2011年
6 仇正中;波浪对海底管线作用的数值模拟[D];大连理工大学;2011年
7 李宁;格子Boltzmann方法研究蒸汽泡在过冷水中的传热传质规律[D];大连理工大学;2011年
8 陈海峰;基于Hopf-Cole变换Burgers方程的有限元方法研究[D];内蒙古大学;2011年
9 王锐;虚拟环境中海洋场景建模与绘制研究及实现[D];安徽大学;2011年
10 吴小霞;GPU高性能计算技术在晶格玻尔兹曼方法模拟中的应用[D];广西师范大学;2011年
【同被引文献】
中国期刊全文数据库 前8条
1 王勇,张泽廷;气液传质相界面湍动现象的实验及分析[J];北京化工大学学报(自然科学版);2002年02期
2 王勇,张泽廷;吸收过程中界面湍动对传质的影响[J];北京化工大学学报(自然科学版);2002年06期
3 沙勇;李樟云;林芬芬;吐松;肖宗源;叶李艺;;气液传质界面湍动现象投影观察[J];化工学报;2010年04期
4 凌锦龙;徐敏虹;俞丽丽;;乙酸乙酯-异丙醇二元系的黏度行为和表面性质[J];化工学报;2012年01期
5 梁博;袁希钢;;界面Rayleigh对流现象模拟的格子Boltzmann方法[J];化学工程;2008年06期
6 祖迎庆,施卫平;用格子Boltzmann方法模拟流场中可变形膜的运动[J];力学学报;2005年02期
7 王亮;陈宗芳;付强;苗仁德;王明;;鱼类自主游动水动力学特性的数值模拟[J];力学学报;2012年01期
8 沙勇 ,成弘 ,于艺红;The Numerical Analysis of the Gas-Liquid Absorption Process Accompanied by Rayleigh Convection[J];Chinese Journal of Chemical Engineering;2002年05期
中国博士学位论文全文数据库 前4条
1 王勇;格子Boltzmann方法在热声领域的应用及热声谐振管可视化实验研究[D];西安交通大学;2009年
2 陈炜;气液界面Rayleigh-Bénard-Marangoni对流现象实验测量及传质研究[D];天津大学;2010年
3 张云;多孔介质中流动的格子Boltzmann模拟[D];中国石油大学;2011年
4 曾建邦;基于气泡生长及多相流动的格子Boltzmann模型及应用研究[D];重庆大学;2011年
【二级参考文献】
中国期刊全文数据库 前10条
1 ;The Mechanism of Drag Reduction around Bodies of Revolution Using Bionic Non-Smooth Surfaces[J];Journal of Bionics Engineering;2007年02期
2 施卫平,胡守信,阎广武;计算可压缩流体流动的LB多速度模型[J];吉林大学自然科学学报;1996年02期
3 胡守信;直管道中绕平板平面流动的格子气仿真[J];计算物理;1989年04期
4 施卫平;;用Lattice Boltzmann方法模拟1/4圆腔内的定常层流运动[J];计算物理;2008年01期
5 阎广武,胡守信,施卫平;用Latice Boltzmann方法确定多孔介质的渗透率[J];计算物理;1997年01期
6 李元香,康立山;线性和拟线性抛物型方程的格点模型[J];计算物理;1997年Z1期
7 邵传平;;钝体尾流控制机理及方法研究进展[J];力学进展;2008年03期
8 祖迎庆,施卫平;用格子Boltzmann方法模拟流场中可变形膜的运动[J];力学学报;2005年02期
9 邵传平;王建明;;较高Re数圆柱尾流的控制[J];力学学报;2006年02期
10 程鸣,凌国平,庄业高;均匀流绕旋转圆柱分离流动的数值模拟[J];水动力学研究与进展(A辑);1990年01期
【相似文献】
中国期刊全文数据库 前1条
1 谷黎;姜水生;;振动Couette流稀薄效应的直接模拟Monte Carlo方法[J];南昌大学学报(理科版);2008年04期
中国重要会议论文全文数据库 前10条
1 ;The Gap between Linearized Poisson-Boltzmann Theory and DLVO Theory[A];第七届全国液体和软物质物理学术会议程序册及论文摘要集[C];2010年
2 ;A Fluctuating Lattice-Boltzmann Model for Direct Numerical Simulation of Particle's Brownian Motion[A];中国颗粒学会第六届学术年会暨海峡两岸颗粒技术研讨会论文集(下)[C];2008年
3 方科;陈锋;苏雁军;;一种环境水质的格子Boltzmann建模研究[A];'2010系统仿真技术及其应用学术会议论文集[C];2010年
4 赫万恒;钱跃竑;;浅水波方程的格子Boltzmann模拟[A];力学与工程应用(第十三卷)[C];2010年
5 李志辉;方明;唐少强;;气体运动论统一算法与DSMC方法描述Boltzmann方程的一致性研究[A];北京力学会第18届学术年会论文集[C];2012年
6 董宇红;;基于格子Boltzmann方程的大涡模拟及湍流时空关联性的研究[A];第十二届现代数学和力学会议论文集[C];2010年
7 戴传山;刘学章;;格子Boltzmann方法用于多孔介质与自由流体开口腔体内自然对流的数值模拟研究[A];中国地球物理学会第二十七届年会论文集[C];2011年
8 梁功有;曾忠;张永祥;张良奇;姚丽萍;邱周华;;封闭方腔内自然对流的格子Boltzmann方法模拟[A];重庆力学学会2009年学术年会论文集[C];2009年
9 吴南健;;采用耦合单电子的量子器件和量子电路(英文)[A];信息科学与微电子技术:中国科协第三届青年学术年会论文集[C];1998年
10 戴传山;李琪;;微纳米颗粒群外绕恒壁温圆管沉积特性[A];中国化工学会2011年年会暨第四届全国石油和化工行业节能节水减排技术论坛论文集[C];2011年
中国重要报纸全文数据库 前1条
1 记者 王进;华中科大船海学院六项目获国家自然基金资助[N];中国船舶报;2010年
中国博士学位论文全文数据库 前10条
1 武芳芳;用格子Boltzmann方法求解一类变系数偏微分方程[D];吉林大学;2012年
2 王颖;一些关于Boltzmann方程的解的存在与稳定性问题[D];中国科学技术大学;2012年
3 刘艳红;电磁波的格子Boltzmann方法研究[D];吉林大学;2013年
4 李婷婷;弹性动力学变形的格子Boltzmann方法研究[D];吉林大学;2013年
5 史秀波;用于波动方程的格子Boltzmann方法及数值模拟研究[D];吉林大学;2010年
6 闫铂;反应扩散方程的格子Boltzmann方法及数值模拟[D];吉林大学;2011年
7 刘芳;格子Boltzmann方法求解偏微分方程的相关研究[D];吉林大学;2011年
8 张建影;用于复Ginzburg-Landau方程的格子Boltzmann方法[D];吉林大学;2011年
9 王勇;格子Boltzmann方法在热声领域的应用及热声谐振管可视化实验研究[D];西安交通大学;2009年
10 李俏杰;格子Boltzmann方法在几类典型偏分方程及高速压制成形中的应用[D];中南大学;2012年
中国硕士学位论文全文数据库 前10条
1 周楠;基于Boltzmann原理烟雾蔓延的仿真[D];西安建筑科技大学;2010年
2 杨士梅;基于Boltzmann方法绕流的数值计算[D];大庆石油学院;2010年
3 光明;由线性化的Boltzmann方程的解得到二粒子Boltzmann方程的色散关系[D];内蒙古大学;2011年
4 皮兴才;格子Boltzmann方法及其在微流动中的应用研究[D];国防科学技术大学;2011年
5 魏泓睿;关于Boltzmann方程的反向解[D];清华大学;2012年
6 李霞;用Lattice Boltzmann方法处理传热传质中的固液相变问题[D];浙江师范大学;2011年
7 宾超华;用格子Boltzmann方法研究中医推拿和分叉管中的血液流[D];广西师范大学;2014年
8 许鹤林;格子Boltzmann方法理论及其在流体动力学中的应用研究[D];复旦大学;2010年
9 叶丽娜;基于格子Boltzmann模型的Kuramoto-Sivashinsky方程的数值模拟研究[D];吉林大学;2011年
10 沈峰;均匀各向同性湍流的格子Boltzmann模拟[D];华中科技大学;2009年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026