海冰动力学的混合拉格朗日-欧拉数值方法

李海; 季顺迎; 沈洪道; 王瑞学; 岳前进

海冰动力学的混合拉格朗日-欧拉数值方法

1.
国家海洋环境预报中心, 北京, 100081
2.
大连理工大学工业装备结构分析国家重点实验室, 辽宁, 大连, 116023
3.
美国Clarkson, 大学土木与环境工程系, Potsdam, 纽约, 13699-5710
4.
大连理工大学工业装备结构分析国家重点实验室, 辽宁, 大连, 116023;美国Clarkson, 大学土木与环境工程系, Potsdam, 纽约, 13699-5710

基金项目: 国家自然科学基金资助项目(40206004;10710101004)

计量
- 文章访问数: 821
- HTML全文浏览量: 6
- PDF下载量: 1220
- 被引次数: 0
出版历程
- 收稿日期: 2007-03-12
- 修回日期: 2007-07-15

A hybrid Lagrangian-Eulerian numerical model for sea ice dynamics

1.
National Marine Environment Forecast Center, Beijing 100081, China
2.
State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China
3.
Department of Civil and Environmental Engineering, Clarkson University, New York 13699-5710, USA
4.
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

摘要

摘要: 综合考虑欧拉坐标下有限差分法(FDM)在海冰动力学计算中的效率,以及拉格朗日坐标下光滑质点流体动力学方法(SPH)对海冰流变行为的精确模拟,本文发展了一种海冰动力学的混合拉格朗日-欧拉(HLE)数值方法。该方法首先在拉格朗日坐标下将海冰离散为若干个具有厚度、密集度的海冰质点,并由这些海冰质点通过Gauss函数对欧拉网格上的海冰参量进行积分插值;然后,在欧拉坐标下对海冰动量方程进行差分计算以确定各网格节点的海冰速度,并由此采用Gauss函数积分插值出拉格朗日坐标下各海冰质点的速度分布;最后,通过对海冰质点运动和分布的计算,确定出各海冰质点的位置、厚度和密集度等参量。采用该HLE方法对规则区域内的海冰堆积过程和涡动风场作用下的海冰动力演化趋势进行了数值试验;最后,采用该HLE方法对渤海海冰的动力过程进行了72h数值模拟,其计算结果与卫星遥感图像和现场观测资料吻合较好。以上计算结果均表明该HLE方法在海冰动力学数值模拟中具有较高的计算效率和模拟精度,可用于海冰动力过程的数值模拟。
- 混合拉格朗日-欧拉法 /
- 海冰动力学 /
- 数值方法 /
- 黏塑性本构模型
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.
- hybrid Lagrangian-Eulerian model /
- sea ice dynamics /
- numerical method /
- viscous plastic constitutive model

HTML全文

参考文献(30)

郭智昌,赵进平.北极海冰数值模拟研究述评[J].海洋与湖沼,1998,29(2):219-228.

MOR1TZ R E,BITZ C M,STEIG E.Dynamics of recent climate change in the Arctic[J].Science,2002,297:1497-1502.

刘钦政,黄嘉佑,白珊,等.全球冰-海洋耦合模式的海冰模拟[J].地学前缘,2000,7(增刊):219-230.

丁德文.工程海冰学概论[M].北京:海洋出版社,1999.

岳前进,季顺迎,于学兵.局地海冰数值预测的冰激平台结构响应计算[J].海洋工程,2003,21(2):32-37.

吴辉碇,杨国金,张方俭,等.渤海海冰设计作业条件[M].北京:海洋出版社,2001.

白珊,吴辉碇.渤海的海冰数值预报[J].气象学报,1998,56(2):139-151.

HIBLER W D.A dynamic and thermodynamic sea ice model[J].Journal Physical Oceanography,1979,9:815-846.

ZHANG Z H.On modeling ice dynamics of semi-enclosed seasonally ice-covered seas[R]//Report Series in Geophysics,Helsinki,2000.

王仁树.渤海海冰的数值试验[J].海洋学报,1984,6(4):572-580.

吴辉碇,白珊,张占海.海冰动力学过程的数值模拟[J].海洋学报,1998,20(2):1-13.

苏洁,吴辉碇,白珊.等.渤海冰-海洋耦合模式:Ⅱ.个例研究[J].海洋学报,2005,27(2):18-28.

ZHANG J,ROTHROCK D.Modeling Arctic sea ice with an efficient plastic solution[J].Journal of Geophysical Research,2000,105 (C2):3325-3338.

HUNKE E C,DUKOWICZ J K.An Elastic-Viscous-Plastic Model for Sea Ice Dynamics[J].Journal of Physical Oceanography,1997,27:1849-1867.

HOPKINS M A,FRANKENSEIN S,THORNDIKE A S.Formation of an aggregate scale in Arctic sea ice[J].Journal of Geophysical Research,2004,109(C01032):1-10.

HUANG Z J,SAVAGE S 13.Particle-in-cell and finite difference approaches for the study of marginal ice zone problems[J].Cold Regions Science and Technology,1998,28:1-28.

OVERLAND J E,MCUTT S L,SALO S,et al.Arctie sea ice as a granular plastic[J].Journal oI Geophysical Research,1998,103(C10):21845-21868.

SCHULSON E M.Compressive shear faults within arctic sea ice:Fracture on scales large and small[J].Journal of Geophysical Research,2004,109(C07016):1-23.

DEMPSEY J P.Research trends in ice mechanics[J].International Journal of Solids and Structures,2000,37,131-153.

HIBLER W D.Sea ice fracturing on the large scale[J].Engineering Fracture Mechanics,2001,68:2013-2043.

GUTFRAIND R,SAVAGE S B.Smoothed particle hydrodynamics for the simulation of broken-ice fields,Mohr-Coulomb-Type rheology and frictional boundary conditions[J].Journal of Computational Physics,1997,134:203-215.

SHEN H T,SU Jun-shan,LIU Lian-wu.SPH simulation of river ice dynamics[J].Journal of Computational Physics,2000,165:752-770.

WANG Z I.A Coastal Sea Ice Model with D/screte Parcel Method[R]//Chrkson University Internal Report.Potsdam:Clarkson University,2000.

SHEN H H,ACKLEY S F,YANG Y.Limiting diameter of pancake ice[J].Journal of Geophysical Research,2004,109(C12-35):1-10.

PRITCHARD R,MUELLER A,Hanzlick D,et al.Forecasting Bering sea ice edge behavior[J].Journal of Geophysical Research,1990,95:775-788.

RHEEM C K,YAMAGUCH1 H,KATO H.Distributed mass/discrete floe model for pack ice rheology computation[J].Journal of Marine Science and Technology,1997,2(2):101-121.

FLATO G M.A particle-in-cell sea-ice model[J].Atmosphere-Ocean,1993,31(3):339-358.

刘煜,吴辉碇,张占海,等.基于质点-网格模式的海冰厚度变化过程数值模拟[J].海洋学报,2006,28(2):14-21.

季顺迎,岳前进,赵凯.渤海海冰动力学的质点网格法数值模拟[J].水动力学研究与进展,2003,18(6):748-760.

SHEN H T,CHEN Y C,WAKE A,et al.Lagrangian discrete parcel simulation of river ice dynamics[J].International Journal of Offshore and Polar Engineering,1993,3(4):328-332.

施引文献

资源附件(0)

访问统计