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

用于复Ginzburg-Landau方程的格子Boltzmann方法

张建影  
【摘要】:格子Boltzmann方法是上世纪末兴起的一种用于流体力学模拟及复杂系统计算的数值方法。它具有算法简单有效、计算并行、复杂边界容易处理等优点,因此受到了广泛关注。近几年,格子Boltzmann方法广泛应用于线性及非线性偏微分方程的求解及模拟。本文用此方法研究复Ginzburg-Landau方程。方程如下:其中β,a,d都是复常数,(?)2是Laplace算子。这是一个具有复变量的振幅方程,描述的是Hopf分叉附近的系统的行为。它在超导、超流、化学系统等领域中,都有广泛的应用。本文重点研究螺旋波及三维涡卷的动力学演化。 复Ginzburg-Landau方程具有复变量和复方程参数,因此本文构建直接用于复变量的格子Boltzmann模型。 在源项的二阶假设下,使用多尺度技术,通过对具有源项及复分布函数的格子Boltzmann方程做Chapman分析,得到一系列关于复平衡态分布函数的偏微分方程。采用二维FHP格子和三维主立方体格子,给出了用于二维及三维复Ginzburg-Landau方程的格子Boltzmann模型。此模型具有复平衡态分布、复附加分布、复格子Boltzmann方程及平衡态分布的复高阶矩形式。同时,针对附加分布函数,我们引入了一个可调参数,以便对模型做出适当调整来改善结果。模型具有实数的单松弛时间因子,因此模型的表达与用于实偏微分方程的模型形式一致。这使得模型保留了用于实方程的模型的简单性。为便于数值模拟,采用Neumann边界条件。边界上的平衡态分布函数由宏观边界的表达式给出。假定在边界处的分布函数为平衡态分布。因为二维及三维的高阶离散速度矩表达式不同,本文分别给出其误差表达式。由于二阶假设导致误差反弹现象,使得模型具有一阶精度。 在本文中,首先研究螺旋波的动力学演化。在对稳定螺旋波的研究中,发现Im(a)有两个值分别标志着其旋转和振动性质的改变。另外,螺旋端点是一种拓扑缺陷点,称为量子涡。在这些点处,理论上|A|=0,但格子Boltzmann模型的结果在2%到8%之间。穿过量子涡剖面的|A|在缺陷点附近线性增加,这与理论预测一致。第二,研究了不同条件下不稳定螺旋波的失稳性质。在特定参数条件下,形成的单螺旋波会破碎成众多的不规则小螺旋,系统进入螺旋湍流态。在螺旋破碎之前,有低振幅波向外传播,在计算区域中心和低振幅波之间存在着振幅的靶波。 与二维问题相比,三维问题的计算消耗明显增加。因此,为研究涡卷性质,将网格调整为50×50×50。为测试三维模型的有效性,先用其模拟准二维空间中的螺旋波。结果表明,当前的网格尺寸可以用来研究三维涡卷的某些性质。然后,在不同的初始条件下,分别给出了各种稳定涡卷的演化,如涡环、涡线为螺旋线的涡卷、不规则多个涡卷和其他类型的涡卷。涡环运动的结果与经典结果一致。螺旋状涡线的涡卷与涡环的运动相似,但最后会形成直涡线的稳定涡卷。对于随机给定的初始条件形成的涡卷,其涡线形成时一般均为弯曲状,经过相互作用,涡线将因拉伸而运动出计算区域以外。因此,在演化中涡线数目逐渐减少,与二维情形不同,最终涡线将全部消失。另外,对于其他两种直线型涡卷,当涡线两端位于计算区域两相邻侧面时,涡线将弯成圆弧,然后沿着半径方向运动至消失;而当涡线两端位于两相对侧面时,将演化成直线型涡卷。因此,涡卷的演化结果为或者随着演化消失,全场变成均匀态,或者经过拉伸变形后形成与坐标轴平行的直涡线的涡卷,这符合正涡线拉伸理论。 对二维螺旋波及三维涡卷的模拟和研究,格子Boltzmann模型都给出了可以接受的结果。它可以用来研究螺旋波或涡卷的其他性质和其他的复偏微分方程。由于复变量的实虚部可以分别看作是两种物质的浓度,因此可以将此模型拓展到其他的双组分反应扩散系统的研究中。
【学位授予单位】:吉林大学
【学位级别】:博士
【学位授予年份】:2011
【分类号】:O35

手机知网App
【相似文献】
中国期刊全文数据库 前10条
1 ;智慧城堡[J];大科技(科学之谜);2011年06期
2 李明;唐涛;陈鼎汉;;V型三能级原子双模光场系统中光场压缩性质[J];物理学报;2011年07期
3 周小莉;;杨图法确定同科电子原子态[J];台州学院学报;2011年03期
4 董海玲;;GI/G/1排队系统队长的强大数定律和中心极限定理[J];数学理论与应用;2011年02期
5 ;[J];;年期
6 ;[J];;年期
7 ;[J];;年期
8 ;[J];;年期
9 ;[J];;年期
10 ;[J];;年期
中国重要会议论文全文数据库 前10条
1 方科;陈锋;苏雁军;;一种环境水质的格子Boltzmann建模研究[A];'2010系统仿真技术及其应用学术会议论文集[C];2010年
2 赫万恒;钱跃竑;;浅水波方程的格子Boltzmann模拟[A];力学与工程应用(第十三卷)[C];2010年
3 董宇红;;基于格子Boltzmann方程的大涡模拟及湍流时空关联性的研究[A];第十二届现代数学和力学会议论文集[C];2010年
4 戴传山;刘学章;;格子Boltzmann方法用于多孔介质与自由流体开口腔体内自然对流的数值模拟研究[A];中国地球物理学会第二十七届年会论文集[C];2011年
5 ;The Gap between Linearized Poisson-Boltzmann Theory and DLVO Theory[A];第七届全国液体和软物质物理学术会议程序册及论文摘要集[C];2010年
6 ;A Fluctuating Lattice-Boltzmann Model for Direct Numerical Simulation of Particle's Brownian Motion[A];中国颗粒学会第六届学术年会暨海峡两岸颗粒技术研讨会论文集(下)[C];2008年
7 梁功有;曾忠;张永祥;张良奇;姚丽萍;邱周华;;封闭方腔内自然对流的格子Boltzmann方法模拟[A];重庆力学学会2009年学术年会论文集[C];2009年
8 李志辉;方明;唐少强;;气体运动论统一算法与DSMC方法描述Boltzmann方程的一致性研究[A];北京力学会第18届学术年会论文集[C];2012年
9 闻斌;;区组大小为3和4的最优格子带洞填充设计[A];苏州市自然科学优秀学术论文汇编(2008-2009)[C];2010年
10 陈瑜;蔡庆东;;基于四叉树网格的格子玻尔兹曼方法[A];北京力学会第14届学术年会论文集[C];2008年
中国重要报纸全文数据库 前10条
1 记者 晁会芹;“甩手老板的格子铺”[N];菏泽日报;2010年
2 常工贸 凃贤平;“创业格子间”首现常州[N];常州日报;2011年
3 本报记者 胡敏;格子铺在困境中突围[N];东方烟草报;2011年
4 本报记者 董超标 本报通讯员 蔡重武;听格子铺掌门说“试水要领”[N];新华日报;2011年
5 通讯员 张生军;密云万名协管员认领“格子田”[N];北京日报;2011年
6 本报记者 杨文;新时尚:租个“格子”搞经营[N];山西日报;2011年
7 尚德荣;格子裙给冬天点颜色[N];中国质量报;2000年
8 实习记者 常青;贵阳悄然流行“格子铺”[N];贵阳日报;2008年
9 本报记者 顾一琼;“格子铺”赢得年轻人青睐[N];文汇报;2008年
10 YMG记者 王波 凌云鹏 通讯员 潘鲁宁;“爱心格子店”培育爱心[N];烟台日报;2010年
中国博士学位论文全文数据库 前10条
1 赖惠林;几类非线性波方程的格子玻尔兹曼模型的数值分析及模拟[D];福建师范大学;2012年
2 刘芳;格子Boltzmann方法求解偏微分方程的相关研究[D];吉林大学;2011年
3 郑林;微尺度流动与传热传质的格子Boltzmann方法[D];华中科技大学;2010年
4 王勇;格子Boltzmann方法在热声领域的应用及热声谐振管可视化实验研究[D];西安交通大学;2009年
5 张建影;用于复Ginzburg-Landau方程的格子Boltzmann方法[D];吉林大学;2011年
6 谭玲燕;用格子Boltzmann方法模拟圆柱的搅动流动及减阻[D];吉林大学;2011年
7 李俏杰;格子Boltzmann方法在几类典型偏分方程及高速压制成形中的应用[D];中南大学;2012年
8 史秀波;用于波动方程的格子Boltzmann方法及数值模拟研究[D];吉林大学;2010年
9 闫铂;反应扩散方程的格子Boltzmann方法及数值模拟[D];吉林大学;2011年
10 李秀文;用格子Boltzmann方法模拟电磁力作用下的流动问题[D];吉林大学;2011年
中国硕士学位论文全文数据库 前10条
1 王辉;基于格子Boltzmann方法的建筑流场仿真研究[D];西安建筑科技大学;2010年
2 张李超;微尺度流体流动和混合的LBM模拟[D];清华大学;2010年
3 李娟;基于格子Boltzmann方法的单通道内两相流的研究[D];大连理工大学;2010年
4 周文宁;复杂微通道内流体流动的格子Boltzmann模拟[D];大连理工大学;2010年
5 邓玉瑾;热声谐振管中声振荡的格子Boltzmann模拟[D];华中科技大学;2009年
6 党亚兵;基于格子Boltzmann的等离子射流数值模拟[D];河南理工大学;2011年
7 陈林婕;几类非线性偏微分方程的格子Boltzmann方法[D];福建师范大学;2010年
8 梁义强;基于格子Boltzmann方法的混相驱替过程模拟研究[D];大连理工大学;2012年
9 沈峰;均匀各向同性湍流的格子Boltzmann模拟[D];华中科技大学;2009年
10 徐超;基于格子Boltzmann方法的海流能水轮机翼型叶片水动力特性研究[D];中国海洋大学;2010年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026