一种基于二维光滑粒子法的流体仿真方法

张海超; 郑丹晨; 边茂松; 韩敏

doi:10.7498/aps.65.244701

搜索

期刊名称:

年:

起始页:

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种基于二维光滑粒子法的流体仿真方法

张海超 , 郑丹晨 , 边茂松 , 韩敏

一种基于二维光滑粒子法的流体仿真方法

张海超, 郑丹晨, 边茂松, 韩敏

物理学报. 2016, 65(24): 244701.

doi: 10.7498/aps.65.244701.

摘要

针对大场景下流体仿真计算复杂度高的问题，本文以浅水方程为基础，提出一种改进的二维光滑粒子方法.该方法中使用光滑粒子法离散二维浅水方程，将水深作为粒子的属性，把计算复杂度降到二维的程度；同时为了提高邻域粒子的搜索效率，提出一种基于动态网格的邻近粒子搜索方法；并使用虚粒子和惩罚力相结合的方法处理边界条件以高效率应对复杂边界；渲染时，首先将粒子映射并插值到规则网格内得到流体表面，避免三维流体表面重构复杂度高的问题，最后利用OpenGL着色语言实现加速渲染，从而达到大场景下流体实时仿真.

关键词:

作者及机构信息

1.
大连理工大学电子信息与电气工程学部, 大连 116023

通信作者: 韩敏, minhan@dlut.edu.cn

基金项目: 国家自然科学基金（批准号：61374154）资助的课题.

参考文献

[1]	Lucy L B 1977 Astron. J. 82 1013
[2]	Prakash M, Rothauge K, Cleary P W 2014 Appl. Math. Model. 38 1534
[3]	Kipfer P, Westermann R 2006 Proceedings of Graphics Interface 2006 Quebec City, Canada, June 7-9, 2006 p41
[4]	Ata R, Soulaïmani A 2005 Int. J. Numer. Meth. Fl. 47 139
[5]	de Leffe M, Le Touzé D, Alessandrini B 2010 J. Hydraul. Res. 48 118
[6]	Lee H, Han S 2010 Visual Comput. 26 865
[7]	Solenthaler B, Bucher P, Chentanez N, Mller M, Gross M 2011 Proceedings of Workshop in Virtual Reality Interactions and Physical Simulations Lyon, France, December 5-6, 2011 p39
[8]	He X W, Liu N, Wang G P, Zhang F J, Li S, Shao S D, Wang H A 2012 ACM T. Graphic. 31 439
[9]	Cornelis J, Ihmsen M, Teschner M 2015 Comput. Graph-UK 52 72
[10]	He J, Chen X, Wang Z Y, Cao C, Yan H, Peng Q S 2010 Visual Comput. 26 243
[11]	Liu H, Qiang H F, Chen F Z, Han Y W, Fan S J 2015 Acta Phys. Sin. 64 094701 (in Chinese)[刘虎, 强洪夫, 陈福振, 韩亚伟, 范树佳2015物理学报 64 094701]
[12]	Han Y W, Qiang H F, Zhao J L, Gao W R 2013 Acta Phys. Sin. 62 044702 (in Chinese)[韩亚伟, 强洪夫, 赵玖玲, 高巍然2013物理学报 62 044702]
[13]	Hu D A, Long T, Xiao Y H, Han X, Gu Y T 2014 Comput. Method. Appl. M. 276 266
[14]	Xia X L, Liang Q H 2015 Environ. Modell Softw. 75 28
[15]	Goswami P, Schlegel P, Solenthaler B, Pajarola R 2010 Proceedings of the 2010 ACM SIGGRAPH/Eurographics Symposium on Computer Animation Madrid, Spain, July 2-4, 2010 p55
[16]	Yang L P, Li S, Hao A, Qin H 2012 Comput. Graph. Forum. 31 2037
[17]	Rodriguez-Paz M, Bonet J 2005 Comput. Struct. 83 1396
[18]	Liu M B, Liu G R, Lam K Y 2003 Comput. Appl. Math. 155 263
[19]	Ihmsen M, Orthmann J, Solenthaler B, Kolb A, Teschner M 2014 Proceedings of Eurographics 2014 State of the Art Reports Strasbourg, April 7-11, 2014 p21

施引文献

[1]	Lucy L B 1977 Astron. J. 82 1013
[2]	Prakash M, Rothauge K, Cleary P W 2014 Appl. Math. Model. 38 1534
[3]	Kipfer P, Westermann R 2006 Proceedings of Graphics Interface 2006 Quebec City, Canada, June 7-9, 2006 p41
[4]	Ata R, Soulaïmani A 2005 Int. J. Numer. Meth. Fl. 47 139
[5]	de Leffe M, Le Touzé D, Alessandrini B 2010 J. Hydraul. Res. 48 118
[6]	Lee H, Han S 2010 Visual Comput. 26 865
[7]	Solenthaler B, Bucher P, Chentanez N, Mller M, Gross M 2011 Proceedings of Workshop in Virtual Reality Interactions and Physical Simulations Lyon, France, December 5-6, 2011 p39
[8]	He X W, Liu N, Wang G P, Zhang F J, Li S, Shao S D, Wang H A 2012 ACM T. Graphic. 31 439
[9]	Cornelis J, Ihmsen M, Teschner M 2015 Comput. Graph-UK 52 72
[10]	He J, Chen X, Wang Z Y, Cao C, Yan H, Peng Q S 2010 Visual Comput. 26 243
[11]	Liu H, Qiang H F, Chen F Z, Han Y W, Fan S J 2015 Acta Phys. Sin. 64 094701 (in Chinese)[刘虎, 强洪夫, 陈福振, 韩亚伟, 范树佳2015物理学报 64 094701]
[12]	Han Y W, Qiang H F, Zhao J L, Gao W R 2013 Acta Phys. Sin. 62 044702 (in Chinese)[韩亚伟, 强洪夫, 赵玖玲, 高巍然2013物理学报 62 044702]
[13]	Hu D A, Long T, Xiao Y H, Han X, Gu Y T 2014 Comput. Method. Appl. M. 276 266
[14]	Xia X L, Liang Q H 2015 Environ. Modell Softw. 75 28
[15]	Goswami P, Schlegel P, Solenthaler B, Pajarola R 2010 Proceedings of the 2010 ACM SIGGRAPH/Eurographics Symposium on Computer Animation Madrid, Spain, July 2-4, 2010 p55
[16]	Yang L P, Li S, Hao A, Qin H 2012 Comput. Graph. Forum. 31 2037
[17]	Rodriguez-Paz M, Bonet J 2005 Comput. Struct. 83 1396
[18]	Liu M B, Liu G R, Lam K Y 2003 Comput. Appl. Math. 155 263
[19]	Ihmsen M, Orthmann J, Solenthaler B, Kolb A, Teschner M 2014 Proceedings of Eurographics 2014 State of the Art Reports Strasbourg, April 7-11, 2014 p21

引用本文:

Citation:

计量

文章访问数: 1515
PDF下载量: 251
被引次数: 0

物理学报. 2016, 65(24): 244701.

doi: 10.7498/aps.65.244701.

一种基于二维光滑粒子法的流体仿真方法

1. 大连理工大学电子信息与电气工程学部, 大连 116023

通信作者: 韩敏, minhan@dlut.edu.cn

基金项目:

国家自然科学基金（批准号：61374154）资助的课题.

收稿日期: 2016-04-27
修回日期: 2016-07-22
刊出日期: 2016-12-05

摘要: 针对大场景下流体仿真计算复杂度高的问题，本文以浅水方程为基础，提出一种改进的二维光滑粒子方法.该方法中使用光滑粒子法离散二维浅水方程，将水深作为粒子的属性，把计算复杂度降到二维的程度；同时为了提高邻域粒子的搜索效率，提出一种基于动态网格的邻近粒子搜索方法；并使用虚粒子和惩罚力相结合的方法处理边界条件以高效率应对复杂边界；渲染时，首先将粒子映射并插值到规则网格内得到流体表面，避免三维流体表面重构复杂度高的问题，最后利用OpenGL着色语言实现加速渲染，从而达到大场景下流体实时仿真.

关键词:

English Abstract

A fluid simulation method based on two-dimensional smoothed particle hydrodynamics

1.
Faculty of Electronic Information and Electrical Engineering, Dalian University of Technology, Dalian 116023, China

Fund Project:
Project supported by the National Natural Science Foundation of China (Grant No. 61374154).

Received Date: 27 April 2016
Accepted Date: 22 July 2016
Published Online: 05 December 2016

Abstract: Smoothed particle hydrodynamics (SPH) method is a kind of meshless method, which is used to solve the problem of fluid simulation without complex operations of the grids. To reduce the computational complexity, SPH method based on the two-dimensional shallow water equations is employed to establish a fluid model. In large scale scenes, taking into account the high computational complexity and the serious distortion problems, in this paper we introduce an improved two-dimensional SPH algorithm according to the shallow water equations. The proposed method with two-dimensional complexity is obtained by discretizing the two-dimensional shallow water equations with SPH, and the depth of water is introduced as the particle's property. The problem of increased amount of calculation cannot be well solved by using traditional neighboring particle search method based on tree structure. To improve the efficiency of search and simplify the search operation of neighborhood particles, in this paper we introduce a point-in-box search algorithm and put forward a neighboring particles searching method on the basis of dynamic grid. Besides, for large scale scenes, by considering that the virtual particle method provides slow computation speed with complex boundary condition, the type-one virtual particles are utilized to ensure that the borders can be calculated precisely by combining the punish force to prevent the phenomenon of particle boundary penetrating. Therefore, a method is further obtained to handle boundary condition efficiently by combining the virtual particles with punish force in this paper. In the process of rendering, the fluid surface is first determined by mapping and interpolating particles into regular grids without the complex reconstruction of surface in three-dimensional. Then, we utilize the bilinear interpolation method to deal with the problem of missing values, and the surface grids are further densified. With OpenSceneGraph three-dimensional render engine, OpenGL Shading Language is adopted to speed up the rendering speed, and in this way, the real-time fluid simulation of large scale scenes can be further achieved. With the basic KD tree searching method employed in the simulations, the comparative experiments are provided to verify effectiveness of the proposed searching method based on dynamic grid. Given the data set obtained from random points, experimental results demonstrate that the method in this paper can be used to solve the problem of neighboring particles searching in large scale scenes. To show the effectiveness of the proposed method on the basis of the virtual particles and the punish force, another experiment based on the collapsing of a water column is further provided. Besides, in this paper we conduct an experiment on a certain actual reservoir terrain to prove that the proposed method can be applied to fluid simulation of large scale scenes.

Keywords:

全文HTML

参考文献 (19)

下载: 全尺寸图片幻灯片

返回文章

搜索

留言板