流速逆变张量隐式求解方法及其在航道港池流场计算中的应用
An implicit method using contravariant velocity components and its application to the calculations in a harbour-channel area
-
摘要: 运用高分辨率的边界适应网格进行流体动力学数值计算时,如何提高计算稳定性和减少计算量成为数值求解的关键性问题.在非正交的边界适应坐标系中,每个动量方程中同时出现了两个交叉方向的水位偏导数项,给隐式求解带来困难,而显式格式下的时间步长由于受与空间步长有关的Courant-Friedrichs-Lewy条件限制,计算量成倍增加.本文从广义曲线坐标系下浅海动力学方程组出发,导出了流速的逆变张量所满足的动量方程组,使方程中的水位偏导数项变成了沿某一协变基向量方向占优的形式,方便地采用了交替方向隐式差分格式,从而提高了计算稳定性并减小了计算量.本文通过对澳门海域航道和港池中流场的计算,证实了该模式是一种进行高分辩率数值计算的有效方法.Abstract: The key problem in the computation of fluid dynamics using fine boundary-fitted grids is how to improve the numerical stability and decrease the calculating quantity.To solve this problem,implicit schemes should be adopted since explicit schemes may bring about a great increase in computation quantity according to the Courant-Friedrichs-Lewy condition.Whereas the adoption of implicit schemes is difficult to be realized because of the existence of two partial derivatives of surface elevations with respect to variables of alternative direction coordinates in each momentum equation in non-rectangular coordinates.With an aim to design an implicit scheme in non-rectangular coordinate,new momentum equations with the contravariant components of velocity vector are derived based on the shallow water dynamic equations in generalized curvilinear coordinates.In each equation the coefficients before the two derivatives of surface elevations have different orders ofmagnitude,i.e.,the derivative with the larger coefficient may play a more important role than that with the smaller one.With this advantage,the ADI scheme can then be easily employed to improve the numerical stability and decrease the calculating quantity.The calculation in a harbour and a channel in Macao nearshore area shows that the implicit model is effective in calculating current fields in small size areas.
-
Key words:
- Numerical model /
- contravariant component of velocity vector /
- implicit scheme /
-
Sheng Y P.Modelling coastal and estuarine processes using boundary-fitted grids.In;Proceedings 3rd Int'l Symposium on River Sedimentation,Garland.New York,1986,1416-1 442 Haase W,Misegrades K.Naar M.Adaptive grids in numerical fluid dynamics.Int.J.Numer.Methods Fluids.1985,5.515-528 刘卓,曾庆存.自适应网格在大气海洋问题中的初步应用.大气科学,1994,18(6):641-648 Shi Fengyan.Sun Wenxin.A variable boundary model of storm surge flooding in generalized curvilinear grids.Int.J.Numer.Methods Fluids,1995,21,641-651 史峰岩,孙文心.岸界适应干湿网格变边界模型及风暴潮漫滩模拟.高科技研究中的数值计算(蒋伯成,张树春主编),长沙:国防科技大学出版社,1995.134-141 Brackbill J U,Saltaman J S.Adaptive zoning for singular problems in two dimensions.Journal of Computational Physics,1982,46.342-368
点击查看大图
计量
- 文章访问数: 836
- HTML全文浏览量: 9
- PDF下载量: 856
- 被引次数: 0