基于格子Boltzmann方法的气—液—固耦合动力学数值研究

来源 :哈尔滨工程大学 | 被引量 : 0次 | 上传用户:xcnyy_007
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
气-液-固耦合是一个多相流(气相、液相和固相)研究问题,是水利机械、液压机械、以及船舶与海洋工程等领域中的重要课题,直接与设备的防护和优化设计相关。流体流经水下结构物时,会发生绕流现象,造成结构上载荷的变化。载荷的波动,轻者会引起结构振动,重者甚至会造成结构的损毁。特别是在军事装备领域,即使是结构的振动,也将严重影响水下航行器(如潜艇)的深隐身性能,危及航行器和舰载人员的安全。受流场环境的影响,设备与流体间的相互作用是一个非常复杂的过程,常伴随气泡的产生,如螺旋桨附近的空化气泡、管道和阀门流中的气泡等。因此,对气-液-固耦合动力学的研究不但具有基础的学术意义,而且能够为实际工程应用提供一定的技术保障。气-液-固耦合动力学研究的关键是耦合问题,即对流-固界面、气-液界面,以及界面间的相互作用进行研究。界面的产生、运动、破裂和融合过程是一个表现在宏观层面,却决定于微观层面的非常复杂的物理行为。目前,虽然主要依赖于实验研究,但实验研究的代价昂贵,而且受限于实验条件和观察手段;理论研究是另一种主要研究手段,但其求解复杂,往往需要大量假设,受限于一般或者通解的情况。数值研究是区别于前两种的第三种主要研究手段,随着现代高性能计算机和数值计算方法的发展,其越来越被认为是研究多相流动力学过程和相关机理的重要手段。本文从水下界面运动特性研究的角度出发,采用基于分子动理论的格子Boltzmann方法(LBM)对流-固界面、气-液界面和两者间相互作用的多相耦合问题进行数值模拟研究。格子Boltzmann方法常被视为介于宏观层面和微观层面之间的一种介观无网格粒子算法,它利用密度分布函数来确定流域内的粒子分布,并通过多尺度展开,自下而上建立与宏观控制方程,如Navier-Stokes方程的联系。不同于传统的有限差分、有限元和有限体积等数值方法,其具有微观动理论特性,有着传统数值方法无法比拟的优点,如清晰的物理背景、天然的可并行性和灵活的参数可控性等。针对气-液-固耦合动力学研究的重要性,依托格子Boltzmann方法在计算流体领域内的优势,本文通过发展不同的数值模型,系统地开展了对流-固界面、气-液界面以及两种界面共存时运动特性的研究工作。主要的研究内容和取得的成果有:(1)提出了一种适用于处理静止流-固界面的格子Boltzmann方法。在该方法中,根据Half-way反弹模型在实际物理边界与毗邻的流体节点位置之间建立一层虚拟边界,并结合有限差分方法获得虚拟边界上的宏观变量,以虚拟边界代替实际边界的壁面反弹过程,计入了对实际边界相对于网格线偏移量的考虑。该方法虽然适用于静止流-固界面的处理,但也能拓展到对特定的运动流-固界面(即在运动过程中不涉及到流体节点和固体节点间的转变)进行处理。基于一些经典算例进行了验证,包括平板间的Poiseuille流、平板间的Couette流、静止圆柱绕流和旋转Couette流。(2)提出了一种适用于处理运动流-固界面的格子Boltzmann方法。在该方法中,引入了浸没边界方法,基于欧拉-拉格朗日坐标体系对流场中复杂运动界面进行描述。采用隐式速度纠正算法,保证了流-固界面处的非渗透、非滑移的边界条件。同时,采用多松弛系数矩阵代替单一的松弛系数提高方法的数值稳定性,最终形成了一种基于隐式速度纠正的多松弛浸没边界-格子Boltzmann方法。通过一些经典算例进行了方法验证,并以此为基础,研究了来流与旋转平板间的耦合特性。(3)提出了 一种适用于处理气-液界面的格子Boltzmann方法及其并行计算模型。该方法基于自由能多相流原理,利用两组格子Boltzmann方程对整个两相流系统进行描述,其中一组用来求解计入了流场粘性和表面张力影响的Navier-Stokes方程,另一组用来求解Cahn-Hilliard方程。通过采用分步操作,提高了方法的数值稳定性,并利用OpenMP多核共享内存式并行计算技术,构建了处理气-液界面的三维并行模型。通过Laplace定律、气泡的形态等验证了方法的可靠性,并以此为基础,对单、双气泡的运动特性进行了研究。(4)研究了狭域中气泡运动特性,包括管道内气泡和壁面驱动剪切流中的气泡。通过管道内气泡运动特性的研究,一方面验证了本文提出的界面处理方法在气-液-固耦合系统中应用的可靠性,另一方面探究了阻塞率和偏心率对气泡的拓扑变形和运动轨迹的影响。通过对壁面驱动剪切流中气泡运动特性的研究,一方面验证了 Capillary数和雷诺数对气泡变形的预报,另一方面对比了壁面驱动剪切流和自由场剪切流中气泡的运动特性,探究了壁面的影响。(5)研究了敞域中气泡与平、曲壁面的耦合特性,通过改变气泡与壁面结构的特征尺寸比和流体属性,观察自由场中气泡与两种类型壁面耦合的过程。采用多松弛格子Boltzmann模型,来提高数值稳定性和可实现的流场粘度范围,建立一维和二维管道流内的声波模拟模型,研究了水下气泡脉动或其他扰动产生的波传播过程。
其他文献
体育报道因其自身写作手法的灵活性等特点,已率先迈向后现代的领地。开始了对传统价值取向的根本性颠覆,并以消解“元叙事”、有意增加文本的空白等手段调动受众参与其问,加入解
从近年来的上海信息化建设来看,创新的IT理念和技术,赢得了越来越多的认可,获得了更为广泛的应用,对于上海社会、经济、文化等领域的发展起到了积极的推动作用。而在众多作出突出
《现代医院管型》杂志总国家级、卫生管理类专业学术期刊。由国家卫生健康委员会主管,中国医院协会和中国医科大学上小。中国医科大学附属盛京医院承办。国内统-刊号CN21-149
摩擦力实验难观察,教师难直观演示,使摩擦力演示实验可视程度低。牛顿三定律可实现摩擦力实验的转换设计,从物理学发展来看牛顿三定律与摩擦力有着千丝万缕的联系,了解两者间
上海市无管局党委日前召开世博保障工作党员大会,动员各党支部充分发挥政治优势、组织优势和群众工作优势,动员广大共产党员充分发挥先锋模范作用,以最好的精神状态、最大的工作
ue*M#’#dkB4##8#”专利申请号:00109“7公开号:1278062申请日:00.06.23公开日:00.12.27申请人地址:(100084川C京市海淀区清华园申请人:清华大学发明人:隋森芳文摘:本发明属于生物技
为探究皂化棉籽油脚与改性膨润土复合的最佳工艺,通过单因素试验和正交试验研究分析了膨润土添加量、皂化棉籽油脚添加量、棉籽油添加量、搅拌时间和搅拌温度对复合工艺的影
2007年是中国消费电子渠道的转折之年,也是最明显的裂变之年。而转折和裂变的背后,一场列强厮杀,在所难免。
In English, as in most other languages, we can say the same thing in a number of ways.For example, to advice. Giving Advice is a Language function. There are va