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

来源 :吉林大学 | 被引量 : 2次 | 上传用户:cmcbst
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
格子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模型都给出了可以接受的结果。它可以用来研究螺旋波或涡卷的其他性质和其他的复偏微分方程。由于复变量的实虚部可以分别看作是两种物质的浓度,因此可以将此模型拓展到其他的双组分反应扩散系统的研究中。
其他文献
本论文采用高精度量子化学计算方法完全活化空间自洽场(CASSCF)和多组态二级微扰理论(CASPT2),对几种与星际和大气相关的重要的小分子自由基,及它们离子基态和激发态的性质进行了较为全面的理论研究。本文不仅为这几种小分子自由基的吸收光谱和光电离能谱提供理论指导,还有助于进一步开展理论和实验方面的工作。主要内容概括如下:1.本论文采用完全活化空间自洽场(CASSCF)和多组态二级微扰理论(CAS
天然抗性相关巨噬细胞蛋白2(Nramp2)是一个拥有重要生物活性功能的膜蛋白,广泛存在于人类、哺乳动物、细菌等生物体内,能够传输生命活动中所需要的二价金属阳离子。在Nramp2的十二个预测跨膜区中,第一跨膜区C端和与之相连的环区(exo-loop1)可能在该蛋白捕获金属离子方面起重要作用。我们采用核磁共振方法研究了Nramp2-exo-loop1片段及基于这个肽片段的系列突变体在水溶液和胶束溶液中
单质元素金属是物质世界的基本构成之一,其高压研究倍受关注。近期研究发现单质元素金属及其金属间化合物在高压下发生了复杂的结构相变,并呈现出一系列新奇的物理现象(如超导电性,超硬、电阻上升、熔点下降等等)。对这些体系的研究不仅能为天文和地球物理的研究提供知识储备,还可以揭示新的物理现象和物理规律。本文利用基于密度泛函理论的第一性原理方法系统的研究了Ca-Li(CaLi2)和Li-H(LiH2、LiH6
声波、电磁波的传播与散射现象是许多工程与科学领域中所关注的重要问题.然而,声波、电磁波的数值计算中却面临着许多挑战,特别是其中的无界域与高波数问题,仍有大量困难没有解决.本文针对这类问题中的数学模型,分别研究了一类近场散射问题的数值计算和无界表面散射问题的适定性及其有限元近似.论文分为三部分.第一部分为论文的第一章和第二章.我们首先阐述了声波、电磁波散射问题的研究背景,以及近场和无界表面散射问题的
原子里德堡(Rydberg)态能级的自然辐射寿命和朗德因子数据是人们了解原子的结构及电子动力学过程的重要资料,因此其实验和理论研究工作一直受到人们的重视。本论文利用时间分辨激光诱导荧光方法和塞曼量子拍技术,测量研究了锡原子5p7p组态9个能级,J = 1的5pnp (n = 10–13, 15–19),J = 2的5pnp (n = 10–13, 15–19, 27, 31, 32)和5pnf (
时标的概念最初是由Aulbach的博士生Hilger在1988年提出来的.起初引入这一概念的主要目的是为了建立一个研究连续和离散动力系统的统一框架.目前时标这一理论已经得到了非常广泛的发展.最近Potsche给出了时标上指数二分性概念.本文深入地研究了时标上指数二分性理论.本文共分为五章:第一章概述了本文研究的背景和预备知识.第二章首先给出时标上二次Lyapunov函数的定义,随后讨论了时标上指数
车辆的被动安全已经成为汽车工业的关键问题之一,而如何提高车辆在碰撞事故中的抗撞性(Crashworthiness)又成为汽车被动安全性分析中的一个关键环节,金属薄壁构件作为最传统、最有效的缓冲吸能装置,在车身吸能体系中得到了最广泛的应用。现今,在有限元模型中大多是利用壳单元模拟,也就是所谓的详细模型。使用壳单元可以保证对实际模型的几何和物理设计特点进行更真实的模拟,因此适合对产品进行详细的设计及评
统编小学语文教科书的"双线组元"编排,特别注重不同年段之间的纵向联系,充分体现了由浅到深的语文要素发展特点。教师应帮助学生通过学习《"诺曼底号"遇难记》,对这类文章产生更丰富的语感体验,使他们在生活中遇到类似语境时能借鉴课文中的表达方法。
利用基因工程操作技术克隆了4个分别含有表达产物为147个氨基酸的酵母转录激活因子核苷酸序列和(或)核定位信号(NLS)和(或)Ⅰ型人免疫缺陷病毒(HIV-1)穿膜肽(Tat)的嵌合核苷酸序列,并分别命名为GAL4/147、NLS-GAL4/147、Tat-GAL4/147和Tat-NLS-GAL4/147。将上述核苷酸序列克隆入原核表达载体pET28a(+),分别构建了原核表达质粒pET28-G(