收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

用格子Boltzmann方法模拟电磁力作用下的流动问题

李秀文  
【摘要】:流体与钝体发生相对运动时,会在钝体表面形成粘性边界层,当雷诺数较大时会在钝体表面某处产生脱体旋涡,使钝体运动受阻增大,使流场的能量减少,从而导致钝体振荡,产生噪音。从工程角度讲,人们希望抑制这一现象,即有效的控制钝体绕流。 近十年来,人们研究了电磁力控制流动边界层分离的方法。其原理是利用电磁力改变流体在边界层的流动,控制流体绕流的形态,进而达到抑制边界层分离、消涡和减阻的目的。该电磁控制方法不同于传统的磁流体力学(MHD),对于传统的磁流体力学,导电的液体的导电率一般为σ~10 6S /m,因此,由该电流电磁感应得到的电磁力足够大,能够影响和改变流体的流动。但对于海水和其他的弱电解质溶液,电流的导电率σ只有10 S /m ,从而由该电流感生的电磁力太小不能够影响流体的流动,所以这需要外加一个电场与磁场耦合产生足够大的电磁力去改变流体的流动。这种特殊的磁流体力学通常称为电磁流体力学(EMHD)。 对于EMHD控制问题,本文分析了电磁场和流体的基本方程。将由电极条和磁极条交叉排列组成的电磁激活板包覆在钝体表面,将其置于弱电解质溶液中,会产生电磁力,即作用于流体上的一种体积力,如果磁条为永久磁铁,电极为恒定电压,则满足Maxwell方程,其中E是电场强度,B是磁感应强度,J是电流密度, u为流动速度矢量,μ0是磁导率,σ是电解质导电率。σ( u×B )是由磁场诱导出来的电流密度。由于感生电流太小,对外加磁场几乎无影响,于是可以忽略不记,则第一个方程可写为j =σE ,因此,我们考虑只有外加电场E和外加磁场B产生的电磁力情况下的电磁流体力学,其表达式为: Lorentz力作为一种体积力出现在流体的控制方程中,因此,基于特征长度和特征速度,弱电解质溶液流动的无量纲形式的流体控制方程为式中,u是流体的速度,p为压力,雷诺数Re = u∞d /γ,γ为流体运动学粘性系数,u∞为来流速度,F是单位质量流体受到的电磁力,电磁力在钝体表面呈法向指数衰减。 对于上述电磁力控制流动的Navier-Stokes方程,通常利用有限体积法、有限差分法、有限元法等数值方法求解。但在模拟电磁流体时会遇到边界处理复杂、并行效率低等困难。所以,我们从介观层次上考虑,利用格子Boltzmann方法计算电磁控制问题,研究边界层流动的特性。与宏观方法不同,该方法的特点是边界处理简单、并行效率高,适合于电磁流体的微观模拟。基于以上认识,本文利用格子Boltzmann方法对电磁流体中一些基本的问题进行研究,通过对这些问题的微观模拟来加深人们对它们的认识。 格子Boltzmann方法是从微观动力学角度出发,将流体的宏观运动看作是大量微观粒子运动的统计平均结果,宏观的物理量由微观粒子的统计平均得到,建立离散的速度模型,在满足质量、动量和能量守恒的条件下,得出粒子分布函数,然后对粒子分布函数进行统计计算。在模拟流体流动和相关的输运现象方面,它具有算法简单、可靠,适合并行计算,边界处理简单等优越性。格子Boltzmann方程的离散化速度模型为:其中, eα是粒子沿第α个方向运动的速度,是沿eα方向运动的粒子速度分布函数和平衡态分布函数,δt是时间步长,τ是无量纲松弛时间, x i是空间的位置,t是时间, gα( xi , t)是与电磁力有关,被表示为:其中F是电磁体力。格子Boltzmann方程中的平衡态分布函数可由Maxwell-Boltzmann分布函数对速度u的二阶Taylor展开获得,可表示为:其中,δx是格子长度, x流体的密度和速度可以由下面两方程计算得到αρu =∑e +δF. 对于钝体表面的曲线边界,我们将采用反弹边界和浸没边界法去处理,利用动量转化法和浸没边界法中的修正速度项去计算钝体表面边界受力。 首先我们利用格子Boltzmann方法数值计算圆柱绕流的电磁控制,研究流向电磁力对边界层结构的影响,利用二阶精度的曲线边界处理技术处理边界,利用动量转化法计算曲线边界的受力。通过计算得到了在不同电磁力作用下的圆柱绕流的流线、等涡线以及阻力系数、升力系数。通过分析各个量的变化,发现在弱电解质溶液中,电磁力能够改变圆柱绕流的边界层结构,使圆柱表面的分离点后移,从而有效抑制旋涡脱落。同时,电磁力与流体流动同向时,电磁力能够减阻、减震;电磁力与流体流动反向时,电磁力有较好的制动控制作用。 其次利用格子Boltzmann方法数值模拟电磁力作用下的椭圆柱绕流,研究电磁力对流场的影响,同样利用二阶精度的反弹边界处理技术处理边界,利用动量转化法计算曲线边界的受力,计算结果表明,电磁力可以改变椭圆柱表面边界层的结构,抑制边界层分离,减少阻力。雷诺数的不同、横轴与纵轴之比的不同,椭圆柱体攻角大小的不同以及电磁力大小的不同等等,这些物理量对边界层的控制效应有很大的影响。对于攻角为零的椭圆柱,在雷诺数相同时,没有电磁力,椭圆柱横轴与纵轴之比越大阻力系数越小;对流场施加相同的电磁力,随着横轴与纵轴之比值的增加,阻力系数下降的幅度越来越小。对于横轴与纵轴之比为定值的椭圆柱,在雷诺数相同时,无电磁力,椭圆柱体攻角越大,尾迹旋涡摆动越大,对流场施加相同大小的电磁力,椭圆柱攻角越大减阻效果越差;在攻角相同时,没有电磁力,雷诺数越大边界层越容易分离,对流场施加相同大小的电磁力,雷诺数越小边界层的分离越容易抑制。 对较大雷诺数Re = 2.4×105下翼型体绕流的电磁控制进行了数值研究,将大涡模拟(LES)引入到LBM中,得到滤波后的LBM演化方程。利用二阶精度的反弹边界处理边界,利用动量转化法计算曲线边界的受力,计算较大雷诺数下的翼型绕流的电磁控制。无电磁力作用时,流体绕过翼型,随着攻角的增大,流体边界层发生分离,这会导致翼型表面阻力增加,同时,计算出该翼型绕流的失速迎角是12。对流场施加电磁力,可以有效控制翼型绕流形态。计算结果显示,电磁力可以控制边界层流体的流动,使边界层的分离点后移,抑制旋涡脱落,有较好的减阻效应。随着电磁力增大,翼型绕流的前缘分离逐渐被抑制,可使翼型的升力逐渐提高,阻力进一步减小,电磁力愈大,升力提升和阻力减小的幅度愈大。 最后我们利用IB-LBM方法数值模拟电磁力作用下的翼型绕流,将浸没边界法和格子Boltzmann方法有效的结合在一起,综合两种方法的优点去计算流体问题,使算法更加简单,边界处理更容易。对于小雷诺数Re = 500,考察了一定的攻角条件下,不同大小的电磁力对翼型表面的流体速度、边界层分离以及升阻力变化的影响。结果表明:作用于翼型上表面与流动方向相同的电磁体积力,可以有效的改变流体边界层结构;较大的电磁力更有利于抑制翼型的前缘涡和后缘涡以及边界层的分离;在足够大电磁力的作用下,可以完全消除前缘涡和后缘涡,抑制边界层分离,并且较大程度地提高翼型表面的升力和减少阻力。本文还拟合出阻力系数C D与电磁作用参数N之间的定量关系式,线性关系C D=? c1 N+ c2,升力系数C L与作用参数N的定量关系式,指数关系


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 王振林,李振亚;自旋S=1座稀释蜂窝格子的伊辛模型[J];物理学报;1990年08期
2 吴坚,梁珍璇,李荫藩;拉格朗日坐标下的三角形网格重分技术[J];数值计算与计算机应用;1994年04期
3 赵立宽;岳晓鹏;许春明;;关于斯坦因豪斯棋盘问题的又一证明[J];中国科技信息;2008年24期
4 ;思维擂台赛[J];发明与创新(综合版);2009年06期
5 子愚;谈格子[J];中国青年科技;1997年11期
6 苏先基;云纹方法——介绍一个有效的应变场测量方法[J];力学进展;1973年06期
7 ;相变与临界观象(Ⅴ)——模型理论综述[J];物理学进展;1981年04期
8 ;游戏中心[J];新知客;2007年11期
9 胡祖熾;数值流体动力学[J];力学进展;1973年05期
10 讴歌;;格子间,一条虫[J];科学大观园;2008年09期
11 ;晶体的结构与结合力和热現象的关系是怎样的?[J];物理;1955年01期
12 张景贵;文建国;;光学格子中孤子的相互作用[J];光电子技术;2007年03期
13 ;游戏中心[J];新知客;2007年10期
14 章文一;;网上卖格子[J];发明与创新(综合版);2008年03期
15 柯大观,张宏,童勤业;格子复杂性和符号序列的细粒化[J];物理学报;2005年02期
16 索洪敏;;超级模数博弈的存在性[J];应用数学学报;2006年02期
17 章睿齐;;成功要马上动手[J];青年科学;2007年12期
18 D.H.克龙 ,孟广娥;用多重二次方程进行重力地形改正[J];物探与化探;1979年02期
19 ;相变与临界现象(Ⅱ)——Ising模型[J];物理学进展;1981年03期
20 李荫藩,曹亦明;任意三角形(或四边形)网格的格子中流体方法[J];计算数学;1981年04期
中国重要会议论文全文数据库 前10条
1 方科;陈锋;苏雁军;;一种环境水质的格子Boltzmann建模研究[A];'2010系统仿真技术及其应用学术会议论文集[C];2010年
2 赫万恒;钱跃竑;;浅水波方程的格子Boltzmann模拟[A];力学与工程应用(第十三卷)[C];2010年
3 董宇红;;基于格子Boltzmann方程的大涡模拟及湍流时空关联性的研究[A];第十二届现代数学和力学会议论文集[C];2010年
4 ;The Gap between Linearized Poisson-Boltzmann Theory and DLVO Theory[A];第七届全国液体和软物质物理学术会议程序册及论文摘要集[C];2010年
5 ;A Fluctuating Lattice-Boltzmann Model for Direct Numerical Simulation of Particle's Brownian Motion[A];中国颗粒学会第六届学术年会暨海峡两岸颗粒技术研讨会论文集(下)[C];2008年
6 李志辉;方明;唐少强;;气体运动论统一算法与DSMC方法描述Boltzmann方程的一致性研究[A];北京力学会第18届学术年会论文集[C];2012年
7 戴传山;刘学章;;格子Boltzmann方法用于多孔介质与自由流体开口腔体内自然对流的数值模拟研究[A];中国地球物理学会第二十七届年会论文集[C];2011年
8 梁功有;曾忠;张永祥;张良奇;姚丽萍;邱周华;;封闭方腔内自然对流的格子Boltzmann方法模拟[A];重庆力学学会2009年学术年会论文集[C];2009年
9 戴传山;李琪;;微纳米颗粒群外绕恒壁温圆管沉积特性[A];中国化工学会2011年年会暨第四届全国石油和化工行业节能节水减排技术论坛论文集[C];2011年
10 闻斌;;区组大小为3和4的最优格子带洞填充设计[A];苏州市自然科学优秀学术论文汇编(2008-2009)[C];2010年
中国博士学位论文全文数据库 前10条
1 李秀文;用格子Boltzmann方法模拟电磁力作用下的流动问题[D];吉林大学;2011年
2 史秀波;用于波动方程的格子Boltzmann方法及数值模拟研究[D];吉林大学;2010年
3 刘芳;格子Boltzmann方法求解偏微分方程的相关研究[D];吉林大学;2011年
4 张建影;用于复Ginzburg-Landau方程的格子Boltzmann方法[D];吉林大学;2011年
5 王勇;格子Boltzmann方法在热声领域的应用及热声谐振管可视化实验研究[D];西安交通大学;2009年
6 武芳芳;用格子Boltzmann方法求解一类变系数偏微分方程[D];吉林大学;2012年
7 李俏杰;格子Boltzmann方法在几类典型偏分方程及高速压制成形中的应用[D];中南大学;2012年
8 郑林;微尺度流动与传热传质的格子Boltzmann方法[D];华中科技大学;2010年
9 闫铂;反应扩散方程的格子Boltzmann方法及数值模拟[D];吉林大学;2011年
10 谭玲燕;用格子Boltzmann方法模拟圆柱的搅动流动及减阻[D];吉林大学;2011年
中国硕士学位论文全文数据库 前10条
1 叶丽娜;基于格子Boltzmann模型的Kuramoto-Sivashinsky方程的数值模拟研究[D];吉林大学;2011年
2 沈峰;均匀各向同性湍流的格子Boltzmann模拟[D];华中科技大学;2009年
3 陈林婕;几类非线性偏微分方程的格子Boltzmann方法[D];福建师范大学;2010年
4 邓玉瑾;热声谐振管中声振荡的格子Boltzmann模拟[D];华中科技大学;2009年
5 王英英;基于格子Boltzmann方法的海流能水轮机大雷诺数模拟研究[D];中国海洋大学;2011年
6 赵胤智;基于格子Boltzmann方法的轨道动边界模拟[D];华中科技大学;2011年
7 皮兴才;格子Boltzmann方法及其在微流动中的应用研究[D];国防科学技术大学;2011年
8 梁义强;基于格子Boltzmann方法的混相驱替过程模拟研究[D];大连理工大学;2012年
9 周楠;基于Boltzmann原理烟雾蔓延的仿真[D];西安建筑科技大学;2010年
10 王辉;基于格子Boltzmann方法的建筑流场仿真研究[D];西安建筑科技大学;2010年
中国重要报纸全文数据库 前10条
1 尚德荣;格子裙给冬天点颜色[N];中国质量报;2000年
2 本报记者 顾一琼;“格子铺”赢得年轻人青睐[N];文汇报;2008年
3 YMG记者 王波 凌云鹏 通讯员 潘鲁宁;“爱心格子店”培育爱心[N];烟台日报;2010年
4 王炳坤陈光明;格子店“试水”鞍山,投资百元可当老板[N];新华每日电讯;2008年
5 罗闻;格子装尽显柔情[N];中国商报;2001年
6 廖晓生;花90元租个“格子”当老板[N];南昌日报;2008年
7 胡彦殊;“格子铺”成长要迈三道坎[N];四川日报;2008年
8 本报记者 傅盛宁实习生 吴翀;“格子店”:低成本创业[N];深圳商报;2008年
9 本报实习记者 王培根;用百万格子赚百万美元![N];电脑报;2005年
10 林志芳;武汉“格子铺”无人问津[N];市场报;2007年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978