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

用守恒型跟踪算法数值模拟流体中的介质界面

Mohammed Aman Ullah(阿曼)  
【摘要】:现今界面的不稳定性在流体力学的研究中是一个很重要的问题,从物理实验到数值模拟等方面有关界面不稳定性的研究都引起了人们的众多关注。一定的物理机制可使得初始状态时流体界面小振幅的扰动演变成大振幅的扰动,这就是所谓的界面不稳定性。引发界面不稳定性的物理机制可以是激波的作用,它造成Richtmyer-Meshkov不稳定性,可以是引力,它造成Rayleigh-Taylor不稳定性,或者速度切变,它造成Kelvin-Helmholtz不稳定性。对界面不稳定性的研究在科学和工程方面有着广泛的应用,例如水下爆炸、雾滴沉积、惯性约束聚变(ICF),及超新星爆炸等自然现象。 对于描述界面不稳定性的流体方程组,如Euler方程组或Navier-Stokes方程组,无论是可压的还是不可压缩情况,其精确解都很难找到,因此通过理论研究所得到的结论是非常有限的。有关界面问题的物理实验通常十分困难,且由于实验条件的限制,物理实验所得到的结果通常不够准确,尤其是在流动的后期阶段关于界面的细节方面不够准确。因此,现今关于界面问题的数值研究引起了广泛关注,人们希望通过数值模拟给出这类问题更清晰的解释。 目前有两类关于界面不稳定的问题引起了人们在数值模拟及实验上的关注,这两个问题分别是激波-气泡作用的问题和单模Richtmyer-Meshkov不稳定性问题。 在激波-气泡问题中,一个圆柱形或者球形的充满了氦气,R22或者其他任何气体的气泡被放置在一个充满空气的方形激波管中,见[28,39,38],一个在空气中运动的激波从右侧撞击气泡并和其作用。它们之间的作用在开始阶段产生了不同的激波,并且在激波相互作用后气泡严重变形。这类问题有两个重要的东西需要研究:1)开始阶段产生的激波以及它们的传播速度和跨过他们的压力变化,2)后期阶段气泡发生严重形变后的形状。有关这个问题的详细讨论将在第四章给出。 在单模Richtmyer-Meshkov不稳定问题中,一个激波管中有两种不同的流体被一个界面分开,且给界面一个正弦波的扰动,然后一激波与这个界面作用。激波的作用使得界面变得不稳定,并且使得界面的小扰动变成具有气泡和尖状的非线性构造[4]。在单模Richtmyer-Meshkov不稳定问题研究中,扰动的振幅和振幅的增长率的理论分析结果,物理实验以及数值模拟的结果是否吻合,是一个非常重要的问题,见[33,32,31,27]。有关这个问题的更详细的讨论将在第五章给出。 对两类问题的数值模拟具有双重的意义。首先,它们是研究流体力学的重要的基准问题,通过对比数值计算结果和在早期通过物理实验所得到的结果,人们可以验证他们的算法的物理上的合理性和数学上的有效性。第二,人们希望通过数值模拟来研究不稳定性的后期阶段,如在激波-气泡问题中后期的界面结构严重变形,以及在单模Richtmyer-Meshkov不稳定问题中的振幅扰动和振幅增长率,以期得到更好的结果。 模拟这两类问题所使用过的数值方法有很多,其中包括基于Euler网格的方法[67,57,77,42],Lagrange或者ALE方法[51,21],level set方法[63,81]VOF方法[74],和界面追踪法[84,24,52]。一般来说,这些方法在开始阶段得到的结果是令人满意的,和物理实验所得到的结论是相吻合的,而且比物理试验更清晰。但是对于后期阶段的数值模拟大多不令人满意,主要困难有复杂界面的构造,网格的大形变,以及数值耗散等,这些困难使人们不能得到很好的结果。因此,对这两类问题的数值模拟对现今的计算流体来说仍是一个巨大的挑战。 在过去的二十年里,对二维的Euler方程组茅德康发展了一种跟踪法,见[52,53]。关于二维Euler方程组可以见下一章的(2.2.1)。该方法和所有的跟踪法一样,几乎没有数值耗散。而且相对于其它跟踪法,有如下两个方面的特点:1)间断曲线由Euler方程组的守恒律性质跟踪,而不是通过曲线上Rienmann问题的解所得到的传播速度跟踪。2)间断曲线的运动由从Euler方程组所导出的一维偏微分方程所刻画,跟踪是通过笛卡儿子网格上局部离散这些一维偏微分方程来实现的。这样的设计使得茅德康的跟踪法相对其它的跟踪法更简单,它在笛卡儿网格上运行,并且没有用到自适应网格。更为重要的是它是守恒的。这个方法已经被成功地运用于各种数值模拟实验,并且在一维和二维的情况下都表现出其高效性,见[52,53,54,55,56,92,93]。 在本文中作者用茅德康的跟踪法对Haas-Sturtevant两个激波-气泡问题和Benjamin的和Meshkov的单模Richtmyer-Meshkov不稳定问题进行数值模拟。对两个激波-气泡相互作用的问题,对空气-氦气的情形数值运算进行到后期的983微秒时刻,对于空气-R22情形运算进行到后期的1020微秒时刻。关于定性方面的比较,作者用Matlab绘制了密度的彩色图像,并和试验获得的阴影图像进行比较[26]。关于定量比较,作者绘制了早期的激波和气泡界面上一些特征点的x—t运动轨迹图,用线性最小二乘法计算了他们的速度,并且和实验及其他方法所得到的结果进行了比较。此外也给出了压力的分布,并和实验的所得到的关于压力的轨迹作比较。不管是定性比较还是定量比较都能看出,用我们的数值方法所得到的结果和实验结果吻合的更好。此外,我们的数值模拟给出了在这两种情况下后期气泡发生严重形变后的图形,这些图形能够显示出气泡在那段时间的拓扑构造。最后讨论了涡旋的产生。 对两类单模RM-不稳定问题,我们计算了振幅和振幅的增长率。在这两种情形下,计算得到的振幅和振幅的增长率和Holmes等人的数值结果吻合的很好[32],尽管我们的结果比他们的更大一些。正如Holmes等人所给出的解释[32],界面的增长速度将会下降,这种界面的减速是由于反射激波和折射激波相互作用产生的次波所导致,它是一种非线性压缩现象。他们给出了一些特定时刻压力的图形,从这些图形上可以捕捉到这种非线性压缩现象。为了证明这个观点,我们用数值方法算出了这些特殊时刻的压力图,从图上也可以捕捉到这种非线性压缩现象。 本文的结构如下:在第二章中,我们简单的描述了控制方程,可压流体力学的Euler方程组,以及介质界面。第三章中,我们简单的描述了茅德康的守恒型追踪法。第四章中我们对Haas-Sturtevant的两个Shock-bubble相互作用问题的进行了数值模拟。我们在4.1节描述了这两个问题;我们在4.2节介绍了关于这两个问题已有的数值计算方面的工作;4.3节描述了数值模拟的设置,其中包括一网格加密的研究,以确保我们的数值模拟在所用网格上是可靠的;4.4节中给出了空气-氦气情况下的的数值结果;4.5节给出了空气-R-22情况下的数值结果。在每一部分中都包含了两部分,对数值模拟的定性和定量的评估以及数值模拟中的所得到的新发现以及理解。在第五章中给出了两个单模RM-不稳定性的基准问题。给出了振幅扰动和振幅增长率的定义,我们对问题进行了数值模拟,并且和Holmes得到的数值结果作比较。此外,对于RM-不稳定性问题,我们还给出了后期的渐进解。我们提出了RM不稳定性的长时间渐近解。计算了长时间的钉状和泡状速度。泡状结构速度与势能流模型的比较,以及钉状结构速度与最佳拟合曲线的比较,最后给出了结论。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 黄唐,毛国良,姜贵庆,周伟江;二维流场、热、结构一体化数值模拟[J];空气动力学学报;2000年01期
2 翟雨阳;刘钦甫;张立德;唐勇;刘振兴;;数值模拟技术在海坨子油田开发调整研究中的应用[J];上海地质;2008年02期
3 陶建华;;海岸、港口工程中波、流数值模拟的关键技术和应用[J];天津科技;2009年02期
4 颜怡霞;黎维芬;郝志明;;平板切削过程的数值模拟和实验研究[J];固体力学学报;2013年S1期
5 远航,于定勇;潮流与泥沙数值模拟回顾及进展[J];海洋科学进展;2004年01期
6 文岑;蒋友祥;赵海燕;;溃坝问题数值模拟研究综述[J];中国科技信息;2010年21期
7 彭道贵;;精细数值模拟技术在海上方案实施跟踪中的应用[J];内江科技;2012年02期
8 陈耀松;创新与构思——力学小议之二[J];力学与实践;2002年03期
9 戴世强;;认识水波驾驭水波——评介陶建华的专著《水波的数值模拟》[J];水动力学研究与进展A辑;2007年02期
10 姜利;王嘉松;肖进;余斌;叶国强;;隧道风塔废气扩散的数值模拟研究[J];水动力学研究与进展A辑;2012年01期
11 马丽娟;徐丰;胡非;张德良;;侧加热腔体内重力波演化过程的数值模拟[J];力学与实践;2006年05期
12 陈锋,戴蓬勃,顾钧扬,黄凤来;BMC固化过程的数值模拟[J];中国塑料;1999年07期
13 都桂梅;胡文斌;杨昌智;;架空方式对某住宅小区室外风环境影响的数值模拟[J];建筑科学;2008年12期
14 朱春东;戴太亮;魏浡;王华君;;CFD数值模拟课程的教学研究与实践[J];教育教学论坛;2012年18期
15 张先锋;黄正祥;;数值模拟应用于弹丸实验教学的实践[J];实验室研究与探索;2007年02期
16 易卫国;杨谦;李群松;;稀薄颗粒流体对弯管冲蚀的数值模拟[J];湖南师范大学自然科学学报;2012年05期
17 张治畴;姚军;;PADE′逼近在FEL数值模拟中的应用[J];计算物理;1996年01期
18 马成刚;邓建军;谢敏;;Rod-pinch二极管理论及数值模拟[J];强激光与粒子束;2007年02期
19 朱桂萍;王健;;混沌摆系统的动力学分析和数值模拟[J];扬州大学学报(自然科学版);2008年03期
20 吴文娟;师永民;王小军;刘洪涛;秦小双;王磊;柴智;李晓敏;;超低渗油气藏非对称压裂数值模拟理论及应用[J];北京大学学报(自然科学版);2012年06期
中国重要会议论文全文数据库 前10条
1 刘慈群;郭白奇;;有弥散和吸附的径向渗流的数值模拟[A];第五届全国水动力学学术会议暨第十五届全国水动力学研讨会文集[C];2001年
2 袁帅;张汉昭;古成钢;;爆炸密封的初步数值模拟研究[A];第五届全国强动载效应及防护学术会议暨复杂介质/结构的动态力学行为创新研究群体学术研讨会论文集[C];2013年
3 孙远翔;徐琳;;混凝土和钢筋混凝土冲击试验的数值模拟[A];第22届全国结构工程学术会议论文集第Ⅲ册[C];2013年
4 陈刚;陈忠富;;45钢Taylor撞击变形的数值模拟[A];中国力学学会学术大会'2009论文摘要集[C];2009年
5 沈成康;吴建成;翁智远;赵永言;;混凝土板在弹体冲击下的响应试验与数值模拟[A];第八届全国结构工程学术会议论文集(第Ⅲ卷)[C];1999年
6 刘仁培;董祖珏;魏艳红;;焊接凝固裂纹数值模拟前处理系统[A];第十次全国焊接会议论文集(第1册)[C];2001年
7 余春祥;;弹丸穿钢靶姿态变化的数值模拟[A];中国工程物理研究院科技年报(2005)[C];2005年
8 买买提明·艾尼;;现代数值模拟方法与工程实际应用[A];第22届全国结构工程学术会议论文集第Ⅰ册[C];2013年
9 刘瑞祥;陈立亮;周建新;廖敦明;;铸造数值模拟技术在重机行业的应用[A];首届中国大型铸锻件制造技术发展论坛论文集[C];2006年
10 张勋;徐丰;;数值模拟倒三角腔内分层流体的瞬态对流[A];北京力学会第20届学术年会论文集[C];2014年
中国博士学位论文全文数据库 前8条
1 Mohammed Aman Ullah(阿曼);用守恒型跟踪算法数值模拟流体中的介质界面[D];上海大学;2011年
2 杨坚;缝洞型碳酸盐岩油藏缝洞单元数值模拟理论与方法研究[D];中国石油大学;2007年
3 张立文;数值模拟技术在金属材料固态加工中的应用[D];大连理工大学;2004年
4 史锋;野外风沙跃移运动的观测实验与数值模拟[D];兰州大学;2010年
5 王志东;三维自由面湍流场数值模拟及其在水利工程中的应用[D];河海大学;2004年
6 曹玉春;流化床垃圾焚烧炉内流动和燃烧污染物生成数值模拟研究[D];浙江大学;2005年
7 张冬娟;板料冲压成形回弹理论及有限元数值模拟研究[D];上海交通大学;2007年
8 Peter Robert Williams;星系的数值模拟[D];中国科学院研究生院(上海天文台);2004年
中国硕士学位论文全文数据库 前10条
1 孙光明;圆球粒子沉降的数值模拟及实验研究[D];浙江大学;2006年
2 齐殿军;大庆长垣高台子油田葡萄花油层精细油藏描述及数值模拟一体化研究[D];吉林大学;2012年
3 屈春叶;离心泵的内部三维流场数值模拟研究[D];太原理工大学;2013年
4 邹星;结构物出水过程的数值模拟[D];大连理工大学;2013年
5 刘春元;热毛细对流的研究及其数值模拟[D];河北工业大学;2003年
6 陈陶;激光与物质热作用过程的数值模拟研究[D];南京理工大学;2004年
7 钱声华;井下它激自增流量理论及数值模拟[D];西南石油学院;2004年
8 严世强;基于GIS的海上泄漏化学品行为数值模拟研究[D];大连海事大学;2004年
9 李建龙;小型离心风扇数值模拟分析[D];南京理工大学;2012年
10 黄全军;涡轮叶栅三维湍流流动的数值模拟[D];哈尔滨工程大学;2005年
中国重要报纸全文数据库 前5条
1 本报实习记者 蒋文雯 通讯员 罗耀华;宝钢数值模拟研究和应用领跑钢铁行业[N];中国冶金报;2012年
2 通讯员 朱江 冯明农;风能项目组研讨数值模拟产品共享[N];中国气象报;2010年
3 杨雄飞;炼钢与环境技术的研究方法[N];世界金属导报;2007年
4 上海交通大学 潘健生;材料热处理过程的数值模拟[N];科技日报;2002年
5 张厚淼 张世明;数值模拟技术“照亮”黑暗油藏[N];中国石化报;2013年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978