| 海冰动力学的混合拉格朗日-欧拉数值方法 |
| |
| 引用本文： | 李海,季顺迎,沈洪道,王瑞学,岳前进.海冰动力学的混合拉格朗日-欧拉数值方法[J].海洋学报,2008,30(2):1-11. |
| |
| 作者姓名： | 李海 季顺迎 沈洪道 王瑞学 岳前进 |
| |
| 作者单位： | 1.国家海洋环境预报中心, 北京, 100081 |
| |
| 摘 要： | 综合考虑欧拉坐标下有限差分法(FDM)在海冰动力学计算中的效率,以及拉格朗日坐标下光滑质点流体动力学方法(SPH)对海冰流变行为的精确模拟,本文发展了一种海冰动力学的混合拉格朗日-欧拉(HLE)数值方法。该方法首先在拉格朗日坐标下将海冰离散为若干个具有厚度、密集度的海冰质点,并由这些海冰质点通过Gauss函数对欧拉网格上的海冰参量进行积分插值;然后,在欧拉坐标下对海冰动量方程进行差分计算以确定各网格节点的海冰速度,并由此采用Gauss函数积分插值出拉格朗日坐标下各海冰质点的速度分布;最后,通过对海冰质点运动和分布的计算,确定出各海冰质点的位置、厚度和密集度等参量。采用该HLE方法对规则区域内的海冰堆积过程和涡动风场作用下的海冰动力演化趋势进行了数值试验;最后,采用该HLE方法对渤海海冰的动力过程进行了72h数值模拟,其计算结果与卫星遥感图像和现场观测资料吻合较好。以上计算结果均表明该HLE方法在海冰动力学数值模拟中具有较高的计算效率和模拟精度,可用于海冰动力过程的数值模拟。
|
| 关 键 词： | 混合拉格朗日-欧拉法 海冰动力学 数值方法 黏塑性本构模型 |
| 文章编号： | 0253-4193(2008)02-0001-11 |
| 收稿时间： | 2007/3/12 0:00:00 |
| 修稿时间： | 2007年3月12日 |
A hybrid Lagrangian-Eulerian numerical model for sea ice dynamics |
| |
| LI Hai,JI Shun-ying,SHEN Hung Tao,WANG Rui-xue and YUE Qian-jin.A hybrid Lagrangian-Eulerian numerical model for sea ice dynamics[J].Acta Oceanologica Sinica (in Chinese),2008,30(2):1-11. |
| |
| Authors: | LI Hai JI Shun-ying SHEN Hung Tao WANG Rui-xue and YUE Qian-jin |
| |
| Institution: | National Marine Environment Forecast Center, Beijing 100081, China;State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China;Department of Civil and Environmental Engineering, Clarkson University, New York 13699-5710, USA;State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China;Department of Civil and Environmental Engineering, Clarkson University, New York 13699-5710, USA;State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China |
| |
| Abstract: | In this study,a hybrid Lagrangian-Eulerian(HLE)method is developed for sea ice dynamics.The method adopts the advantages of the high computational efficiency of Particle-in-Cell method(PIC)and the high numerical accuracy of Lagragian Smoothed Particle Hydrodynamics(SPH).In the HLE model,the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations.These ice variables are interpolated to the Eulerian grid nodes using the Gaussian int erpolation function.Thereafter,the FDM is used to determine the ice velocities on the Eulerian grid nodes,then the velocities of Lagrangian ice particles are interpolated from the grid velocity with the Gaussian function.The Lagrangian ice particles are displaced with the interpolated velocities,and the thicknesses and concentrations of ice particles are then determined based on their new locations.With the HLE numerical model,the ice ridging process in a rectangular pool is simulated,and the simulated results are validated with the analytical solution.The method is also applied to the simulation of sea ice dynamics under vortex wind field.At last,the HLE model is used to simulate the sea ice dynamics of Bohai Sea,and the simulated concentration,thickness and velocity match the satellite images and the field observed data wel. |
| |
| Keywords: | hybrid Lagrangian-Eulerian model sea ice dynamics numerical method viscous plastic constitutive model |
| 本文献已被 维普 万方数据 等数据库收录！ |
| 点击此处可从《海洋学报》浏览原始摘要信息 |
| 点击此处可从《海洋学报》下载免费的PDF全文 |
|
