长江科学院院报 ›› 2013, Vol. 30 ›› Issue (5): 97-102.DOI: 10.3969/j.issn.1001-5485.2013.05.021

• 科技简报 • 上一篇    

ELCIRC源程序代码分析

杨飞1,杨文俊1,杨森2   

  1. 1.长江科学院 水利部江湖治理与防洪重点实验室,武汉 430010;2.重庆市武隆县水务局,重庆 408500
  • 收稿日期:2012-04-25 修回日期:2013-04-28 出版日期:2013-04-28 发布日期:2013-04-28
  • 作者简介:杨 飞(1985-),男,山东泰安人,硕士研究生,从事水力学及河流动力学研究,(电话)13810609717(电子信箱)yangfeihaoyun@163.com。
  • 基金资助:
    国家自然科学基金面上项目(51079008)

Algorithm of ELCIRC Source Code

YANG Fei1,YANG Wen-jun1,YANG Sen2   

  1. 1.Key Laboratory of River Regulation and Flood Control of MWR, Yangtze River Scientific Research Institute, Wuhan 430010, China; 2.Water Affairs Bureau of Wulong County, Chongqing City,Chongqing 408500,China
  • Received:2012-04-25 Revised:2013-04-28 Online:2013-04-28 Published:2013-04-28

摘要: 开放源程序代码ELCIRC是采用基于水平向无结构网格、垂向z坐标体系和半隐格式的欧拉-拉格朗日有限体积/有限差分方法来解浅水方程。研究了ELCIRC中除控制方程组求解以外的部分,分析了源程序在插值计算、拓扑结构、定解条件、分层信息和干湿法等技术上的具体操作,尤其是对欧拉拉格朗日模型特有的逆向追踪算法做了详细阐述。同时讨论了一些细节问题,在不恰当的地方给出一些参考和建议,避免小数做分母、相近数相减引起的较大误差。然而,这些经过实际应用的算法原理,可以为我们开发水动力学模型提供有价值的参考。

关键词: ELCIRC , 逆向追踪 , 源代码 , 三维水动力学模型 , 算法

Abstract: Open-source code ELCIRC(Eulerian-Lagrangian CIRCulation) solves the shallow water equations using a semi-implicit Eulerian-Lagrangian finite volume/finite difference method with horizontally unstructured grids and vertically unstretched z-coordinates. ELCIRC source code aside from the governing equation solution is analyzed in this paper. The operation of interpolation,topological structure,definite condition,hierarchy information,and wetting and drying method are described in detail. Backtracking,as the key technology of ELCIRC,is expounded comprehensively. Moreover,some detailed problems are discussed,and a few references and suggestions are given to avoid the errors caused by similar number subtraction and by employing small number as the denominator. ELCIRC has been applied in practice and could serve as a valuable reference for developing hydrodynamic models.

Key words: ELCIRC , backtrack , source code , 3-D hydrodynamic model , algorithm

中图分类号: