基于物理的波浪实时模拟方法

J4 ›› 2014, Vol. 36 ›› Issue (03): 508-512.

基于物理的波浪实时模拟方法

张志毅，吕之华

（西北农林科技大学信息工程学院，陕西杨凌 712100）

收稿日期:2012-09-20 修回日期:2012-12-19 出版日期:2014-03-25 发布日期:2014-03-25
基金资助:
教育部留学回国人员科研启动基金资助项目(K314020901)；中央高校基本科研业务费专项基金资助项目(Z109021004)

A method of real time wave simulation based on physical

ZHANG Zhiyi,L Zhihua

(College of Information Engineering,Northwest A&F University,Yangling 712100,China)

Received:2012-09-20 Revised:2012-12-19 Online:2014-03-25 Published:2014-03-25

摘要/Abstract

摘要：

为解决广域范围内波浪状态的实时计算和可视化问题，结合无粘滞流体力学的物理模型，提出了一种以三角形为控制单元的有限体积简化算法。该算法的优点在于：以不规则边界区域的非规则三角网格为基础，通过简化通量向量分裂方法获得三角形控制单元的边界数值通量，能快速逼近二维浅水方程的解进而模拟非规则边界浅水的实时流动。实验结果显示，所提方法能在符合现实世界物理规律的前提下较好地实现大规模波浪的实时可视化模拟。

关键词: 浅水方程, 限定Delaunay三角剖分, 有限体积法, 通量向量分裂法

Abstract:

In order to solve the realtime computation and visualization of the wave state in the wide area, by combining the physical model of nonviscous fluid dynamics, a simplified finite volume algorithm is proposed, which uses triangle as control unit. The advantage of the proposed algorithm is that: Based on irregular triangular mesh of irregular boundary area, the boundaries' numerical flux of triangles can be obtained by simplifying fluxvector splitting method. Then shallow water flow with irregular boundary can be realtime simulated by fast approaching the solution of the twodimensional shallow water equations. Experimental results show that the realtime visual simulation of the wave can be better achieved by using this method in the physical laws of the real world under the premise.

Key words: shallow water equations;constrained Delaunay triangulation;finite volume method;flux-vector splitting method

张志毅，吕之华. 基于物理的波浪实时模拟方法[J]. J4, 2014, 36(03): 508-512.

ZHANG Zhiyi,L Zhihua. A method of real time wave simulation based on physical [J]. J4, 2014, 36(03): 508-512.

参考文献

［1］Tan Weiyan. Computational shallow water hydrodynamics—the application of the finite volume method［M］. Beijing:Tsinghua University Press, 1998.(in Chinese)
［2］Causon D M, Mingham C G. Introductory finite difference methods for PDEs［EB/OL］.［20100516］. http://bookboon.methodsforpdesebook.
［3］Li Renxian.The basis of finite volume method［M］. Beijing:National Defence Industry Press, 2008.(in Chinese)
［4］Versyeeg H K,Malalasekera W.An introduction to computational fluid dynamics［M］. England:Longman Group Ltd, 1995.
［5］Bowyer A. Computing dirichlet tessellations［J］. The Computer Journal, 1980,24(2):162166.
［6］Chippada S, Dawson C N，Martinet M L, et al. Compute methods appL［J］. Mech Eng，1998,151(12):105129.
［7］Rogers B M, Fujihara M, Borthwick A G L. Adaptive Qtree Godunovtype scheme for shadow water equation［J］. International Journal for Numerical Methods in Fluids, 2001, 35(3):247280.
［8］Pan C, Dai S, Chen S. Numerical simulation for 2D shallow water equations by using Godunovtype scheme with unstructured mesh［J］. Journal of Hydrodynamics, 2005,18(4):475480.
［9］Yang Qin. Constrained Delaunay triangulation mesh generation technology［M］. Beijing:Publishing House of Electronics Industry, 2005.(in Chinese)
［10］Du Min. Unstructured triangulation grid generation and its application in 2D hydrodynamics model［D］. Tianjin:Tianjin University, 2005.(in Chinese)
［11］Zhang Zhiyi, Konno K, Tokuyama Y. Curve mesh reconstruction based on mountain contours ［J］. Journal of the Institute of Image Information and Television Engineers, 2006, 60(11):18031810.
［12］Zhao Dihua, Yao Qi, Jiang Yan, et al. FVS scheme in 2D depthaveraged flowpollutants modeling［J］. Advances In Water Science, 2002, 13(6):701705.(in Chinese)
［13］Zhang Liqiong,Cui Guangbai,Yang Jue.Flux vector splitting and finite volume method in the application of hydraulic modeling［J］. Water Resources and Hydropower Engineering, 2001,8(32):811.(in Chinese)
［14］Roshandel A,Hedayat N,Kiamanesh H.Simulation of dam break using finite volume method［J］. World Academy of Science, Engineering and Technology, 2010,48(4):112115.
［15］Spekerijse S P. Second order accurate upwind solutions of the 2D steady Euler equations by the use of a defect correction method［C］∥Proc of the 2nd European Conference on Multigrid Method, 1986:285300.
［16］Steger J L, Warming R. Flux vector splitting of the inviscid gasdynamic equations with application to finite difference methods［J］. Journal of Computational Physics, 1981, 40(2):263293.
［17］Pang Mingyong. The real time dynamic simulation method of two dimensional water waves based on the discrete model［J］. Shuili Xuebao, 2007,38(11):13581363.(in Chinese)
［18］Gao Zhiyi, Yu Fujiang, Xu Fuxiang. Preparing 3D animation of sea wave forecasting［J］. Marine Science Bulletin, 2011,30(2):173178.(in Chinese)
［19］Li Yunfei,Cheng Tiantian,He Wei.Simulation method for lake surface wave［J］. Journal of System Simulation,2009,21(23):75077510.(in Chinese)
附中文参考文献：
［1］谭维炎. 计算浅水动力学—有限体积法的应用［M］．北京：清华大学出版社, 1998.
［3］李人宪. 有限体积法基础［M］. 北京：国防工业出版社, 2008.
［9］杨钦. 限定Delaunay三角网格剖分技术［M］. 北京：电子工业出版社,2005.
［10］杜敏. 非结构三角网格生成及其在二维水动力学模型中的应用［D］. 天津：天津大学, 2005.
［12］赵棣华, 姚琪, 蒋艳, 等. 通量向量分裂格式的二维水流水质模拟［J］. 水科学进展, 2002, 13(6):701705.
［13］张丽琼, 崔广柏, 杨珏. 通量向量分裂格式及有限体积法在水流模拟中的应用［J］. 水利水电技术,2001,8(32):811.
［17］庞明勇. 基于离散模型的二维水波实时动态模拟方法［J］. 水利学报, 2007,38(11):13581363.
［18］高志一，于福江，许富祥. 海浪预报三维动画计算原理与制作方法［J］. 海洋通报, 2011,30(2):173178.
［19］李云飞，程甜甜，何伟. 一种湖面波浪模拟的方法［J］. 系统仿真学报, 2009,21(23):75077510.

[1]	王军华欧阳洁. 系数快速振荡的抛物型方程求解的新方法[J]. J4, 2007, 29(7): 67-70.
[2]	杨继业欧阳洁. 多重网格技术对SIMPLER算法的加速性能研究[J]. J4, 2006, 28(8): 81-85.