ELCIRC河口模型优点之一是水平法向速度的全导数采用逆向追踪计算。这种欧拉-拉格朗日模拟必须采用精确的追踪算法。就ELCIRC现有模型中插值方法的双线性映射误差进行了一定的讨论,给出另外一种较为简化的迭代算法。通过弯道水槽算例分析可知,迭代算法在计算精度和计算有效性上与原有的程序基本一致。
The numerical algorithm of ELCIRC (Eulerian Lagrangian CIRCulation) has an important advantage: the total derivatives of the horizontal normal velocity are discretized using Lagrangian back tracking. The use of accurate tracking algorithms is strongly recommended for this Eulerian Lagrangian simulation. Tracking errors in inverse bilinear mapping for quadrangles is discussed and a more simple iterative method is recommended. A calculation example of bend flume shows that this iterative method has the same accuracy and computational efficiency with the original method.
[1]OLIVEIRA A,BAPTISTA A M. On the Role of Tracking on Eulerian Lagrangian Solutions of the Transport Equation[J]. Advances in Water Resources,1998,21(7) 539-554.
[2]ZHANG Ying long,BAPTISTA A M,MYERS E P. A Cross scale Model for 3 D Baroclinic Circulation in Estuary Plume Shelf Systems I. Formulation and Skill Assessment [J]. Continental Shelf Research,2005,25(7/8) 935-972.
[3] 杨飞,杨文俊,杨森. ELCIRC源程序代码分析[J]. 长江科学院院报,2013,30(5) 97-102.(YANG Fei,YANG Wen jun,YANG Sen. Algorithm of ELCIRC Source Code[J]. Journal of Yangtze River Scientific Research Institute,2013,30(5)97-102.(in Chinese)
[4]李健.香溪河水华数值模拟研究[D].北京清华大学,2012.(LI Jian. Numerical Simulation on Xiangxi River Algal Bloom[D]. Beijing Tsinghua University,2012.(in Chinese)