Artificial resistance of Eulerian-Lagrangian method in 3D numerical model
-
摘要: 在三维数模的实践中,发现使用欧拉-拉格朗日方法(ELM)离散对流项会引发计算水位随时间步长减小而增大的异常现象,本文结合数值实验和数学分析对其进行了阐明和研究。数学推导分析表明,这些异常是由于ELM中的线性插值与水平流速的非线性垂线分布之间的不适应所致,插值误差随着插值点离单元侧面距离的减小而增大。在此基础上,对异常现象的本质进行了解释。此外,通过使用垂向σ坐标网格上的三维水动力学模型进行计算,得出ELM的时间阻力现象具有普遍性,但对边界适应能力强的网格系统可以减轻其影响。Abstract: The artificial resistance of the Eulerian-Lagrangian Method (ELM) in a 3D hydrodynamic model is observed and characterized by a rising free-surface elevation when the time step decreases. The phenomenon is illustrated and investigated by combining the numerical experiments and the formal analysis. The analysis indicates that the abnormality results from the linear interpolation of the ELM which does not match the nonlinear distribution of the velocity. Based on the analysis, the phenomena are explained. Moreover, the computations by a 3D hydrodynamic model using the sigma-coordinate indicate that the artificial resistance of the ELM is universal, but the grid which follows the topography and free surface better will alleviate its negative influence.
-
Key words:
- hydraulics /
- numerical methods /
- three-dimensional /
- Eulerian-Lagrangian method /
- artificial resistance
-
[1] BAPTISTA A M.Solution of advection-dominated transport by Eulerian-Lagrangian methods using the backwards method of characteristics[D].Cambridge:Massachusetts Institute of Technology,1987. [2] 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):539-554. [3] DIMOU K.3-D hybrid Eulerian-Lagrangian particle tracking model for simula ting mass transport in coastal water bodies[R].Cambridge:Massachusetts Institute of Technology,1992. [4] JANKOWSKI J A.A non-hydrostatic model for free surface flows[D].Hannover:Hannover University,1999. [5] ZHANG Y L,BAPTISTA A M,MYERS E P.A cross-scale model for 3D baroclinic circulation in estuary-plume-shelf systems:I.Formulation and skill assessment[J].Continental Shelf Research,2004,24(18):2187-2214. [6] CASULLI V,ZANOLLI P.Semi-implicit numerical modeling of nonhydrostatic free-surface flows for environmental problems[J].Mathematical and Computer Modeling,2002,36:1131-1149. [7] 王志力,陆永军,左利钦.强潮河口三维无结构网格盐度数学模型[J].海洋工程,2008,26(2):43-53.(WANG Zhi-li,LU Yong-jun,ZUO Li-qin.Unstructured 3D baroclinic model of current and salt for strong tidal estuary[J].Ocean Engineering,2008,26(2):43-53.(in Chinese)) [8] OLIVEIRA A,FORTUNATO A B.Towards an oscillation-free,mass conservative,Eulerian–Lagrangian transport model[J].Journal of Computational Physics,2002,183:142-164. [9] WARNER J C,SHERWOOD C R,ARANGO H G,et al.Performance of four turbulence closure models implemented using a generic length scale method[J].Ocean Modeling.2005,8(1/2):81-113. [10] 胡德超.三维水沙运动及河床变形数学模型研究[D].北京:清华大学,2009.(HU De-chao.Studies on the three-dimensional numerical for free-surface flow and suspended-sediment transpont[D].Beijing:Tsinghua University,2009.(in Chinese)) -

计量
- 文章访问数: 12
- HTML全文浏览量: 2
- PDF下载量: 582
- 被引次数: 0