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

30 2 Vol1 30, No1 2 2008 3 ACTA OCEANOLOGICA SINICA March 2008 - 1 2* 3 2, 3 2 李海 , 季顺迎 , 沈洪道 , 王瑞学 , 岳前进 ( 1. , 100081; 2. , 116023; 3. Clarkso , Potsdam, 13699- 5710) : 综合考虑欧拉坐标下有限差分法( FDM) 在海冰动 学计算中的效率, 以及拉格朗日坐标下 光滑质点流体动 学方法( SPH) 对海冰流变行为的精确模拟, 本文发展了一种海冰动 学的混合 拉格朗日- 欧拉( HLE) 数值方法该方法首先在拉格朗日坐标下将海冰离散为若干 个具有厚度密 集度的海冰质点, 并由这些海冰质点通过Gauss 函数对欧拉网格上的海冰参量进行积分插值; 然 后, 在欧拉坐标下对海冰动量方程进行差分计算以确定各网格节点的海冰速度, 并由此采用Gauss 函数积分插值出拉格朗日坐标下各海冰质点的速度分布; 最后, 通过对海冰质点运动和分布的计 算, 确定出各海冰质点的位置厚度和密集度等参量采用该HLE 方法对规则区域内的海冰堆积 过程和涡动风场作用下的海冰动 演化趋势进行了数值试验; 最后, 采用该HLE 方法对渤海海冰 的动 过程进行了72 h 数值模拟, 其计算结果与卫星遥感图像和现场观测资料吻合较好以上计 算结果均表明该HLE 方法在海冰动 学数值模拟中具有较高的计算效率和模拟精度, 可用于海 冰动 过程的数值模拟 : 混合拉格朗日- 欧拉法; 海冰动 学; 数值方法; 黏塑性本构模型 : P731. 15 :A :0253-4 193( 2008) 02-0001- 11 1 , [ 8) 9] , [ 1) 2] ; , ( FDM) Baltic Okhotsk , [ 7, 9) 12] , [ 3) 4] , ( LSOR) ( ADI) , , , [ 5) 7] [ 13) 14] , , , : , , , , 20 70 , ; , , : 2007-03-12; : 2007- 07-15 : ( 10710101004) : ( 1971) ) , , , , , E-mail: lihai@ mefc. gov. c * : , E-mail: jisy@ dlut. edu. c 2 30 , , , [ 14 ) 16] 2 , 2. 1