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

基于大涡模拟的移动粒子半隐式法研究及其应用

潘徐杰  
【摘要】: 自由表面是船舶与海洋工程领域的一个重要研究课题,对自由表面问题的模拟是船舶与海洋工程领域,计算流体力学的难点之一。虽然一些基于欧拉网格的自由表面处理方法在模拟带自由表面问题时能够取得不错的结果,但当遇到大变形自由表面、液面破碎、飞溅等问题时,基于欧拉网格的数值方法将使计算量增加并会遇到一定困难。移动粒子半隐式法(Moving-Particle Semi-Implicit Method, MPS)是最近兴起的一种基于拉格朗日(Lagrange)粒子的无网格法, MPS法中粒子的运动性使得该方法在自由表面尤其是大变形自由表面的模拟中有着出色的表现。研究表明,MPS法在船舶与海洋工程领域有着广阔的应用前景,但是鉴于该方法提出不久,算法自身仍有一些缺陷制约着MPS法的发展。因此,本文的目的在于研究MPS法理论基础,针对该方法目前比较突出的几点不足之处,讨论缓解或解决这些缺点的方法,进一步完善MPS法的理论,使其在船舶与海洋工程领域有着更好的适应能力,最后再结合大涡模拟,将MPS法应用至具有大变形自由表面的液体晃荡问题中。本文主要有以下几点内容和创新点: 一、使用了一种防止粒子碰撞或穿越的人工修正步。MPS法中的拉格朗日粒子有着很强的运动随机性,使得在模拟过程中粒子之间相互碰撞以及粒子穿越边界条件的情况频繁发生,通过使用该人工修正步以避免上述情况的发生。数值实验表明,使用了该人工修正步后,模拟可以顺利进行,合适的修正距离可以保证粒子碰撞或穿越现象在整个模拟中出现的概率极低。 二、讨论了不同核函数在不同类型模拟中的适用性。核函数是MPS法中的重要参数,也是MPS法中为数不多的可调参数之一,以往的研究中,往往通过调整不同的核函数来获取更好的模拟结果,但是没有针对核函数的特性进行系统的研究。通过讨论不同核函数在不同类型模拟中的适用性,可以指导在不同模拟中核函数的选取情况。 三、研究了压力振荡现象,提出了以面积-时间平均为扩展的缓解方法,讨论了该方法对压力振荡现象的缓解效果。长期以来,压力振荡现象一直是制约MPS法广泛应用的最大原因,在早期MPS法的研究中,大多的注意力放在自由表面的表现上,压力问题因振荡现象而被忽视。通过使用面积-时间平均法,对不同的面积平均方式和不同的时间平均段进行讨论,选取合适的面积-时间平均方法,可将MPS法的压力结果缓解至工程实践可接受的范围之内,这为MPS法的进一步应用做好了准备,并且基于面积-时间的压力平均法也符合MPS法理论的实际意义。 四、提出一种全新的基于邻居搜索的自由表面混合判别方法。通过粒子数密度来判别自由表面是MPS法在自由表面处理上的一大特色,该自由表面判别方法不但物理意义明确,而且基本不消耗计算资源,有着极高的效率。但是由于MPS法的算法流程问题,该方法会不可避免的产生自由表面的误判情况,这样不但会给模拟结果带来误差,还会加重模拟中的压力振荡现象,实际模拟表明这种自由表面误判的情况较为严重。新的自由表面判别方法由于是基于邻居搜索的,既而可以避免误判情况的发生,考虑到使用粒子数密度来判别自由表面的方法有着很高的效率,所以将基于邻居搜索的自由表面判别方法同基于数密度的自由表面判别方法结合使用,只在自由表面误判较为频繁的粒子数密度段使用邻居搜索,这样在不影响模拟效率的前提下,很大程度上避免了非自由表面粒子被误判的情况发生。 五、基于大涡模拟的MPS法在具有大变形自由表面晃荡问题中的应用。工程中大部分问题都是湍流问题,而以往绝大多数的MPS法研究中均未考虑湍流,所以实现MPS的大涡模拟有着积极的意义。论文通过Smagorinsky涡粘模型,实现了MPS法的大涡模拟,讨论了大涡模拟中微分模型的合适作用距离。最后结合之前提出的方法,成功的模拟了3种不同情况下,具有大变形自由表面的晃荡问题,包括在二维水平方向上晃荡的共振方形液箱,二维平面内横摇的模型化的共振LNG液舱,以及在二维平面内3方向耦合运动的全尺寸共振LNG棱柱形液舱。上述不同方法在模拟中同时应用,均有很好的表现,模拟能够得到满意的结果。 在具体的使用中,人工修正步能够很好的避免模拟中的碰撞和穿越问题;不同核函数的适用性能够准确的指导不同类型模拟的核函数选取;基于邻居搜索的混合自由表面判别方法能够准确的识别自由表面,并且对模拟的效率整体影响不大;而基于面积-时间平均的压力振荡缓解方法可以得到很好的压力结果,将原MPS法完全不能使用,幅度极大的振荡结果缓解至工程上可接受范围之内,并且压力结果同实验结果、其他数值模拟结果都极为接近,在晃荡的拍击压力上也有不错的表现。 实践证明,提出的新技术或方法能够适应于MPS法的应用,并且很好的弥补了原MPS法在相关方面的不足。另外经过进一步完善MPS法大涡在大变形自由表面晃荡问题中有良好的表现,能够模拟出各种物理现象,这对工程实践有着重要的指导意义。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 王兵,张会强,王希麟,郭印诚,林文漪;不同亚格子模式在后台阶湍流流动大涡模拟中的应用[J];工程热物理学报;2003年01期
2 杨武兵,张会强,王希麟,郭印诚,林文漪;空间发展平板混和层流动的大涡模拟-3D与2D模拟的比较[J];工程热物理学报;2004年06期
3 刘有军;平面突扩流动非稳定性的大涡模拟[J];计算力学学报;1998年02期
4 王雪梅;大涡模拟模式研究对流层湍流[J];热带气象学报;2001年03期
5 樊洪明,何钟怡,王小华;弯曲管段内流动的大涡模拟[J];水动力学研究与进展A辑;2001年01期
6 武生智;赵玉军;;沙粒在沙波纹表面上运动的大涡模拟[J];甘肃科学学报;2010年03期
7 苏铭德;直方管内充分发展湍流的大涡模拟第一部分[J];空气动力学学报;1995年01期
8 魏英杰,何钟怡;槽道中方形障碍物绕流的大涡模拟[J];水动力学研究与进展A辑;2003年04期
9 额日其太,邹正平,王强;湍流边界层拟序结构的大涡模拟研究[J];工程热物理学报;2004年04期
10 曾卓知;王有远;刘文光;;后台阶流动的大涡模拟[J];南昌航空工业学院学报(自然科学版);2006年01期
11 肖隐利;宋文艳;杨顺华;;大涡模拟在超音速燃烧领域的应用研究[J];世界科技研究与发展;2008年03期
12 童兵,祝兵,周本宽;方柱绕流的数值模拟[J];力学季刊;2002年01期
13 金晗辉,罗坤,郑友取,樊建人,岑可法;大涡模拟二维平面尾迹湍流拟序结构[J];自然科学进展;2002年12期
14 蒋平,郭印诚,张会强,王希麟,林文漪;矩形射流中的流向涡分布特性及作用[J];清华大学学报(自然科学版);2005年08期
15 李斌;吴颂平;;基于非Favre过滤的大涡模拟[J];北京航空航天大学学报;2007年10期
16 应展烽;陈志华;徐晓磊;董刚;范宝春;;可压缩混合层流动的三维大涡模拟[J];南京理工大学学报(自然科学版);2008年03期
17 苏铭德,R.Friedrich;沿纵向弯曲壁面湍流边界层的大涡模拟[J];力学学报;1984年04期
18 杨建明,吴玉林,曹树良;流体机械中高雷诺数流动的大涡模拟[J];工程热物理学报;1998年02期
19 蔡树群,王文质;大涡模拟及其在海洋湍流数值模拟中的应用[J];海洋通报;1999年05期
20 吴海燕;王昱;孙明波;汪洪波;周进;;颗粒两相流自由剪切层大涡模拟[J];航空计算技术;2008年03期
中国重要会议论文全文数据库 前10条
1 苗世光;蒋维楣;;城市边界层大涡模拟初步研究[A];新世纪气象科技创新与大气科学发展——中国气象学会2003年年会“城市气象与科技奥运”分会论文集[C];2003年
2 张凯杨;许春晓;崔桂香;张兆顺;;平面射流及标量扩散的大涡模拟[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
3 周磊;徐春龙;柳贺;韩君;;4气门柴油机进气及喷雾过程的大涡模拟研究[A];第十七届全国大功率柴油机学术年会论文集[C];2011年
4 陆夕云;庄礼贤;;热分层剪切湍流大涡模拟动力学亚格子尺度模型[A];全国流体力学青年研讨会论文集[C];2001年
5 陈建国;符松;;综合体育馆假想火灾的数值模拟[A];计算流体力学研究进展——第十一届全国计算流体力学会议论文集[C];2002年
6 杨庆山;李启;;亚格子模型在钝体绕流大涡模拟中的比较[A];第十四届全国结构风工程学术会议论文集(下册)[C];2009年
7 刘奕;郭印诚;张会强;王希麟;林文漪;;不同速度比下平面射流拟序结构发展过程的大涡模拟[A];自然、工业与流动——第六届全国流体力学学术会议论文集[C];2001年
8 王永康;石磊;刘向军;;循环流化床内气相流动的大涡模拟研究[A];冶金企业自动化、信息化与创新——全国冶金自动化信息网建网30周年论文集[C];2007年
9 李廷秋;;大涡模拟风浪中动物体的流固耦合问题[A];第二十三届全国水动力学研讨会暨第十届全国水动力学学术会议文集[C];2011年
10 钱炜祺;;可压缩翼型绕流大涡模拟初步研究[A];第十届全国分离流、旋涡和流动控制会议论文集[C];2004年
中国博士学位论文全文数据库 前10条
1 张斌;大涡模拟滤波网格分析及网格自适应控制研究与应用[D];上海交通大学;2011年
2 朱伟亮;基于大涡模拟的CFD入口条件及脉动风压模拟研究[D];北京交通大学;2011年
3 关晖;湍流射流的大涡模拟及其动力学特性研究[D];河海大学;2005年
4 史万里;涡轮流场大涡模拟与拟序结构研究[D];南京航空航天大学;2012年
5 李雪松;基于可压缩流方法的大涡模拟及其工程化应用[D];中国科学院研究生院(工程热物理研究所);2006年
6 周磊;大涡模拟在燃油喷雾过程及多孔介质发动机中应用的研究[D];大连理工大学;2010年
7 刘戈;基于代表性互动小火焰(RIF)模型的内燃机喷雾燃烧大涡模拟研究[D];大连理工大学;2012年
8 邹高万;船舶机舱火灾热流场特性研究[D];哈尔滨工程大学;2005年
9 陈巨辉;基于大涡模拟—颗粒二阶矩的两相流动与反应数值模拟[D];哈尔滨工业大学;2013年
10 白志刚;三维波浪大涡模拟及其对结构的动态效应研究[D];天津大学;2004年
中国硕士学位论文全文数据库 前10条
1 周阳;建筑物脉动风压场大涡模拟入流边界条件的生成研究[D];广东工业大学;2012年
2 鲍龙;隧道内高速列车外部绕流大涡模拟初探[D];兰州交通大学;2013年
3 胡应明;采用大涡模拟技术发展NAPA软件的分离流动模拟能力[D];南京航空航天大学;2009年
4 孙学尧;基于缓解压力振荡MPS法的数值水池研究[D];上海交通大学;2011年
5 朱凤琦;低雷诺数条件下可压缩槽道湍流的大涡模拟[D];哈尔滨工业大学;2012年
6 江川;基于大涡模拟评估GRAPES模式对对流边界层的模拟性能[D];中国气象科学研究院;2013年
7 梁旭;平壁湍流边界层拟序结构的大涡模拟[D];西安建筑科技大学;2004年
8 肖红林;槽道湍流的大涡模拟[D];天津大学;2004年
9 陈建平;低亚音速和跨音速矩形柱绕流大涡模拟研究[D];清华大学;2012年
10 郭大为;带壁面抽吸的槽道湍流大涡模拟[D];天津大学;2005年
中国重要报纸全文数据库 前1条
1 本报记者 徐玢;汗洒科研主阵地[N];科技日报;2008年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978