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

基于SPH的流体仿真数值算法及工程应用研究

焦培刚  
【摘要】: 计算机数值仿真逐渐成为解决现代工程和科学问题的一条重要途径。数值仿真能为理论提供测试和检验,有助于对复杂的物理问题加深认识,甚至还能帮助解释和发现新现象。 基于网格的数值方法虽然已经有广泛的应用,但是在很多方面仍存在不足之处,比如在计算流体动力学中的大变形、运动物质交界面、自由表面等问题时,由于网格产生畸变导致计算误差过大或无法进行,从而使其在许多问题的应用上受到限制。 近年来,无网格法倍受关注,这种方法在许多应用中都优于传统的基于网格的有限元法、有限差分法以及有限体积法等数值方法。本文依托山东大学虚拟工程研究中心和新加坡国立大学ACES实验室,系统地研究了新一代无网格方法—光滑粒子流体动力学方法(SPH)在应用及扩展过程中的相关关键技术,实现了SPH方法的两相耦合应用和三维应用,在此基础上对三维管道气力输送过程进行了数值仿真,扩展了SPH方法的工程应用领域,为SPH方法真正辅助试验创造了条件,具有重要的理论和应用价值。 论文的主要内容如下: 系统总结了基于网格的数值方法,指出其在很多方面存在不足之处,由于网格畸变而使其在许多问题的应用上受到限制,在此基础上,详细分析了新一代无网格法—光滑粒子动力学方法的基本思想和求解问题的过程,由于问题域粒子之间没有任何连接,运算中使用的粒子取决于当前局部分布的粒子,得出了SPH方法的无网格属性和自适应性。系统分析了拉格朗日型的Navier-Stokes方程,应用SPH粒子近似法推导出了Navier-Stokes的密度方程、动量方程和能量方程的SPH表达式。 研究了流体动力学SPH数值计算中的关键技术,指出了为了使算法适合模拟不同的流体特性问题,必须对算法进行特别地处理。在SPH方法中有两种方式对密度进行展开:一种是密度求和法,另一种是连续密度法。虽然密度求和法所需要的计算量大,但因为密度求和法体现了SPH近似法的本质,所以使用得较为广泛。连续性密度法主要用于仿真具有强间断的问题。在SPH方法中,用核函数来近似狄拉克δ函数,分析了在实际应用中常用的两种核函数:高斯核函数和三次样条核函数。 在SPH的应用中,边界条件的处理既是该方法的优点,也是目前的薄弱环节。研究了处理固定边界条件的两种类型的虚粒子,第一种类型的虚粒子设置在固定边界上,第二种类型的虚粒子分布在边界的外部,通常在边界条件不断变化的场合下使用。第二种类型的虚粒子按以下的方式构造,即给定一个实粒子i,则在边界外与实粒子对称处分布一个虚粒子,这些虚粒子具有与相对应实粒子相同的压力和密度,但速度方向相反。为了避免数值震荡,提高算法的稳定性,防止粒子间相互接近时的非物理穿透,在SPH方法的动量方程中引入了人工粘度来进行修正。 本文使用跳蛙法进行时间积分,跳蛙法的优点是计算时所需要的存储量低,而且在每一次计算中只需要进行一次优化估值。给出了SPH算法的程序结构,SPH的基本方法和SPH中其它数值方面的伴随算法使得SPH程序具有一些特殊性质。这些特殊性质都包含在时间积分过程的主循环中。对腔内剪切流动问题、冲击管问题相关算例进行了研究,测试了SPH方法在不同的流体动力学问题中的实用性,对于以上问题,SPH方法均可得到令人满意的结果。 在分析SPH方法单相流基本原理和离散思想的基础上,提出了SPH的两相耦合流动仿真方法。由于SPH方法不需要背景网格,是一种纯无网格方法,在对问题进行粒子离散化的过程中,布置的粒子本身具有物理属性,因此,只要能够正确处理具有相同或者不同属性的粒子之间的相互作用,通过整体粒子的运动分布就能够辨析出不同物质之间的交界面,这样也就能够描述具有不同属性的两相流体的运动情况。在SPH两相耦合流动仿真中,可以通过将固体相或液体相离散为与气体相一样的SPH粒子,但其上携带的是固体相或液体相的物理属性来实现。 研究了保证SPH两相耦合流动数值仿真过程正常进行需要做的技术改进。将密度正则化引入SPH两相耦合流动数值仿真方程,解决了两相流体交界面附近由于密度差异过大导致的边界效应,提高了密度不同或不连续的两相流体交界面处的精度。研究了SPH两相耦合流动中的人工状态方程,通过在两相耦合流动中密度小的物质的状态方程中添加气体内聚力项,解决了密度较小的粒子的逃逸问题。研究了SPH两相耦合流动中的速度修正,引入XSPH中的速度修正方法,防止了粒子间的相互穿透,使得到的两相耦合流动交界面更加清晰和光滑。综合使用SPH两相耦合流动方法,实现了气泡上浮和溃坝两相流动的数值仿真,得到的结果验证了本文提出的SPH两相耦合流动方法的正确性和可行性。 实现了SPH方法在三维流体流动中的技术应用。分析了全配对搜索法和链表搜索法,将Bucket搜索优化算法和树形搜索优化算法引入到SPH方法中,并给出了程序实现,有效地提高了光滑粒子法的计算效率。深入研究了粒子对的相互作用,由于粒子间的相互作用是基于点对点方式的,因此可在SPH仿真分析中应用成对相互作用法来提高计算效率和节省存储空间。成对相互作用法是通过应用最近相邻粒子搜索法来实现的,并为之后进行的SPH求和法存储必须的数据。研究了提高计算结果精度的几种方法:在对核近似连续性理论分析的基础上,提出了初始粒子的配置方法,须尽量将粒子均匀分布,并且使所有粒子质量相同或者质量呈连续变化;采用一种简单高效的基于密度变换的可变光滑长度技术;对光滑长度进行优化和松弛更新,给出了优化和松弛过程的步骤和参数,有效解决了现有SPH方法在处理三维流体流动过程中计算误差大的问题。 应用SPH方法实现了三维管道气力输送过程的数值仿真,分析了影响管道气力输送过程的主要因素,推导了粒子速度和输送气体速度的关系,提出了拟合方法。研究了气力输送过程中管道边界的实现方法,给出了使用应变率及应力的迭代近似法求解粘性力的方法,并分析了管道中输送物料的分布及压力情况。提出了管道边界层流效应的仿真方法:通过设置边界虚粒子和实粒子的粘度参数,而在管道中部的粒子并不设置粘度参数,按无粘处理,能够模拟边界层流效应。输送过程中气流碰到刚性壁面后将产生壁面热量,研究了壁面热量能量转换,给出了考虑人工热量条件下的SPH能量方程。合理简化影响要素,综合运用前面的相关技术,实现了管道输送问题的SPH程序原型。 研究了管道气力输送SPH数值仿真方法的VC++实现技术,由于Fortran具有接近数学公式的自然描述,并且计算精度高,在计算机里具有很高的执行效率,是目前流行较广的一种面向过程的适用于科学计算的高级语言,因此整个SPH程序实现选用Fortran来编写,但是它在实现人机对话、交互处理方面不很理想,界面也不够美观,图形处理功能也比较欠缺。VC++是目前个人计算机上深受欢迎的面向对象的程序设计环境。VC++具有强大的功能,但在科学计算方面实现起来却比较复杂。因此,用VC++和Fortran进行混合编程可以编制出具有友好界面和良好计算功能的程序。本文通过综合运用SPH基本原理,SPH两相耦合流动方法,SPH三维仿真技术,实现了气力输送过程的SPH数值仿真,在保持计算效率的基础上,通过VC++调用Fortran生成的动态链接库(.dll)文件,实现了SPH程序的交互功能。 综上所述,本文提出的SPH两相耦合流动仿真方法,将SPH方法延伸到了两相耦合流动中,完善了现有的光滑粒子流体动力学应用理论。本文实现了SPH方法在三维仿真流动中的应用,将树形搜索优化算法引入到SPH方法中,有效地提高了光滑粒子法的计算效率,特别是在粒子数量较多的时候,计算效率优势更加明显。本文给出了提高计算结果精度的几种方法。首次将SPH方法应用在三维管道气力输送仿真中,扩展了SPH方法的工程应用领域,给出了考虑粘性和壁面热量转换的计算方法,分析了管道中输送物料的分布及压力情况,提出了管道边界层流效应的仿真方法。在保持计算效率的基础上,通过VC++调用Fortran生成的动态链接库文件,实现了SPH程序的交互功能,为SPH方法真正辅助试验创造了条件。这些工作具有重要的理论和工程应用价值。 本课题得到了国家建设高水平大学公派研究生项目(留金出[2007]3020)的支持。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 常满,郭秀梅;海面电磁波散射的数值仿真[J];地球物理学报;1986年02期
2 杨翊仁;;对称分段线性振子的次谐分叉分析[J];西南交通大学学报;1992年06期
3 李桦;贾地;范晓樯;李晓宇;;高超侧压进气道前/后掠的数值分析和比较[J];推进技术;2007年01期
4 周国宝;王林;;冲击载荷作用下海洋平台的数值仿真研究[J];中国海洋平台;2007年02期
5 王勇;王德荣;解东升;;道面浅埋爆炸效应数值模拟[J];解放军理工大学学报(自然科学版);2007年06期
6 邓荣兵;金先龙;陈向东;杜新光;沈建奇;李政;;爆炸冲击波作用下桥梁损伤效应的数值仿真[J];上海交通大学学报;2008年11期
7 高振国;;基于Flash建立的虚拟液压回路系统的研制[J];实验技术与管理;2008年08期
8 邓优林;;基于MatLab的一个新的混沌系统的仿真分析[J];现代计算机(专业版);2009年04期
9 焦培刚;周以齐;王喜仓;;自由表面流动的三维SPH数值仿真研究[J];山东大学学报(工学版);2009年06期
10 郭兰平;俞建宁;张建刚;张旭东;漆玉娟;;基于改进RBF神经网络对股价的演变预测[J];兰州交通大学学报;2010年01期
11 徐英;时家明;林志丹;;撞击物形状和速度对高速撞击结果的影响[J];弹箭与制导学报;2010年02期
12 陈鑫;余川;潘文武;曹晓阳;张生帅;;基于抛物面天线阵的空间功率合成技术[J];强激光与粒子束;2010年10期
13 林红;胡玉林;薛世峰;;ANSYS在工科材料力学教学中的应用[J];电脑学习;2011年01期
14 赵维加;黄健飞;杨斌;;运动弹性杆的一类数值仿真算法[J];系统仿真学报;2011年02期
15 李群宏;徐德贵;;一个类Lorenz系统的动力学分析[J];重庆理工大学学报(自然科学版);2011年02期
16 蔡虹;;基于耦合模理论均匀光纤光栅光谱仿真研究[J];太原师范学院学报(自然科学版);2011年01期
17 蔡虹;魏志峰;;基于传输矩阵理论线性啁啾光纤光栅光谱仿真研究[J];山西大学学报(自然科学版);2011年02期
18 曾亮;邹强;彭绍雄;;用LS-DYNA软件辅助“终点效应学”的教学[J];实验技术与管理;2011年04期
19 黄承绪;混沌动力系统的数值仿真研究[J];湖北工学院学报;2002年02期
20 蒋书敏;田立新;王学弟;;Arneodo混沌系统的控制[J];江苏大学学报(自然科学版);2005年06期
中国重要会议论文全文数据库 前10条
1 乐莉;王俊佳;闫军;;超高速撞击数值算法应用对比研究[A];’2004系统仿真技术及其应用学术交流会论文集[C];2004年
2 李中华;梁杰;;高空羽流流场数值仿真[A];第十届全国计算流体力学会议论文集[C];2000年
3 陈颖;邓小燕;;基于动脉系统旋动流原理的一种腔静脉滤器及其流场数值仿真[A];2008年全国生物流变学与生物力学学术会议论文摘要集[C];2008年
4 周凯;张洪军;殷兴景;苏中地;;旋进旋涡流量计结构特性数值模拟研究[A];中国仪器仪表学会第十一届青年学术会议论文集[C];2009年
5 詹杰民;李毓湘;;求解高阶Boussinesq方程数值方法的研究[A];第十六届全国水动力学研讨会文集[C];2002年
6 张厚美;;混合式盾构机在硬岩中的掘进性能数值仿真[A];中国土木工程学会第十二届年会暨隧道及地下工程分会第十四届年会论文集[C];2006年
7 王震;晏麓晖;曾首义;;平头弹侵彻金属靶板的SPH数值仿真[A];第16届全国结构工程学术会议论文集(第Ⅲ册)[C];2007年
8 袁化成;江海;刘栋梁;郭荣伟;;成发集团某系列TRT机组进排气蜗壳数值仿真研究[A];2010年全国炼铁生产技术会议暨炼铁学术年会文集(下)[C];2009年
9 陈众迎;龙连春;;三维编织复合材料的力学性能研究现状[A];北京力学会第15届学术年会论文摘要集[C];2009年
10 顾宝栋;陈忠明;;单元尺寸对冲击仿真精度影响分析[A];中国航空结构动力学专业组第十六届学术交流会论文集[C];2008年
中国博士学位论文全文数据库 前10条
1 焦培刚;基于SPH的流体仿真数值算法及工程应用研究[D];山东大学;2010年
2 杨启;大型船舶出坞数值仿真研究与试验[D];武汉理工大学;2003年
3 董曼玲;直流条件下厘米级间隙流注特性数值仿真及实验研究[D];华中科技大学;2012年
4 高慧婷;Fe/Zn复合材料及液压阻尼器抗震性能数值仿真研究[D];吉林大学;2004年
5 周超;膏体推进剂流变与输送特性研究[D];南京理工大学;2009年
6 蒋光彪;基于DLR型k-ε模型的AMG方法对渐扩管路内紊流数值仿真[D];湘潭大学;2009年
7 王晓煜;超磁致伸缩微位移执行器的系统建模与控制方法研究[D];大连理工大学;2007年
8 王芙昱;基于计算流体力学技术的个体化颅内动脉瘤的血流动力学模型与实验[D];中国人民解放军军医进修学院;2008年
9 颜世军;离心场中广义弹性体的动力学建模与数值分析[D];重庆大学;2011年
10 仇玉良;公路隧道复杂通风网络分析技术研究[D];长安大学;2005年
中国硕士学位论文全文数据库 前10条
1 林燕;铝电解槽槽膛内形在线监测的研究及应用[D];中南大学;2007年
2 刘绍华;微型燃气轮机燃烧室三维紊流数值仿真[D];重庆大学;2008年
3 赵金石;高压发电机并网问题的研究[D];哈尔滨理工大学;2009年
4 王真;单臂架起重机臂架结构变幅运动与动力数值仿真[D];武汉理工大学;2003年
5 王博;旋风分离器内气固两相运动的数值仿真研究[D];西安建筑科技大学;2003年
6 窦磊;火药破碎对压力异常现象影响的数值仿真[D];南京理工大学;2002年
7 王勇;重载货车中Huck系列连接器力学行为数值仿真研究[D];大连交通大学;2005年
8 朱孟巍;船舶与海洋平台碰撞的动力特性研究[D];武汉理工大学;2006年
9 张伯熹;航空发动机涡轮叶片的包容试验和数值仿真[D];浙江大学;2006年
10 侯荣国;磨料水射流切割性能和喷嘴内外流场的仿真研究[D];山东大学;2006年
中国重要报纸全文数据库 前10条
1 通讯员 徐瑞哲 张小果;三维仿真施工应用国内地下特大工程[N];中国建设报;2006年
2 记者  徐瑞哲;崇明隧道未动一土完成“彩排”[N];解放日报;2006年
3 沈建斌;兴澄特钢应用双套管气力输送技术[N];世界金属导报;2011年
4 记者 严风华;现代设计需建立CAE信息链[N];中国船舶报;2007年
5 王友明;压力容器计算精度取得大突破[N];中国质量报;2008年
6 河北粮食学校(050051) 李桂芹;粮食加工厂气力输送的常见故障与排除[N];河北科技报(农村版);2001年
7 北京信息与控制研究所 王正中;仿真学科的研究与思考[N];大众科技报;2007年
8 董仲平;气力输送及设备与选型[N];中国建材报;2006年
9 金敏 记者 陈铁;武汉“精细爆破”技术实现有害效应控制[N];中国质量报;2010年
10 ;仿真技术欲创学术新高[N];中国电脑教育报;2006年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978