基于Mohr-Coulomb准则的黏弹-塑性海冰动力学本构模型

针对中小尺度下海冰动力作用过程中的漂移、重叠和堆积特征,在连续介质力学基础上建立了一个黏弹塑性海冰动力学本构模型.该模型主要包括四部分,即海冰在屈服前的Kelvin Vogit黏弹模型、海冰塑性屈服的Mohr Coulomb准则、屈服后相关联的正交流动法则和影响海冰强度的静水压力.采用黏弹塑性本构模型对规则海域内的海冰堆积过程进行了数值试验,其结果与传统黏塑性模型、改进的黏塑性模型和经典冰坝理论的解析解相符,可合理地模拟海冰堆积的动力过程.为进一步检验该黏弹塑性本构模型的适用性,对渤海海冰的动力过程进行了48h的数值模拟,计算结果与海冰卫星遥感现场观测资料一致.通过与传统的黏塑性和改进的黏塑性模型结果的对比分析,进一步验证了黏弹塑性模型在中小尺度海冰动力学模拟中的可靠性.在以上海冰堆积的数值试验和渤海海冰动力学模拟中均采用了光滑质点流体动力学(SPH)方法.     

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

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

Numerical simulation for dynamical processes of sea ice

ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｆｏｒｄｙｎａｍｉｃａｌｐｒｏｃｅｓｓｅｓｏｆｓｅａｉｃｅＷｕＨｕｉｄｉｎｇ，ＢａｉＳｈａｎ，ＺｈａｎｇＺｈａｎｈａｉａｎｄＬｉＧｕｏｑｉｎｇ（ＲｅｃｅｉｖｅｄＭａｙ１６，１９９６；ａｃｃｅｐｔｅｄＪａｎｕａｒｙ...     

海冰模式和数值预报

根据渤海冰情,在海冰动力学和热力学研究基础上,提出一种海冰动力一热力模式。该模式包括平整冰、堆积冰和开阔水三要素,采用粘一塑性本构关系计算冰内应力。利用参数化方法处理变形函数和热力增长率。该冰模式与潮流模式耦合研究渤海冰－潮相互作用,并与ＥＣＯＭ－Ｓｉ海洋模式耦合进行预报试验。该冰模式与大气模式联接,或应用气象部门数值天气预报产品进行海冰业务数值预报。本文还给出了预报结果。     

A hybrid Lagrangian-Eulerian numerical model for sea-ice dynamics

A hybrid Lagrangian-Eulerian(HLE) method is developed for sea ice dynamics,which combines the high computational efficiency of finite difference method(FDM) with the high numerical accuracy of smoothed particle hydrodynamics(SPH).In this 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 Eularian gird nodes using the Gaussian interpolation function.The FDM is used to determine the ice velocities at Eulerian grid nodes,and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also.The thicknesses and concentrations of ice particles are determined based on their new locations.With the HLE numerical model,the ice ridging process in a rectangular basin is simulated,and the simulated results are validated with the analytical solution.This method is also applied to the simulation of sea ice dynamics in a vortex wind field.At last,this HLE model is applied to the Bohai Sea,and the simulated concentration,thickness and velocity match the satellite images and the field observed data well.     

海冰动力过程的改进离散元模型及在渤海的应用

海冰的断裂、重叠和堆积等离散分布特性广泛地存在于极区和副极区的不同海域,并对海冰的生消、运移过程有着重要影响。针对海冰在不同尺度下的离散分布特点,发展海冰动力过程的离散元方法有助于完善海冰数值模式,提高海冰数值模拟的计算精度。为此,本文针对海冰生消运移过程中的非连续分布和形变特性,发展了适用于海冰动力过程的改进离散元模型(MDEM)。不同于传统离散元方法,该模型将海冰离散为具有一定厚度、尺寸和密集度的圆盘单元。海冰单元设为诸多浮冰块的集合体,其在运移和相互接触碰撞过程中,依照质量守恒发生单元尺寸、密集度和厚度的相应变化。基于海冰离散性和流变性的特点,该模型采用黏弹性接触本构模型计算单元间的作用力,并依据Mohr-Coulomb准则计算海冰法向作用下的塑性变形及切向摩擦力。为验证该模型的可靠性,本文对海冰在规则水域内的运移和堆积过程进行了分析,离散元计算结果与解析值相一致;此外,对旋转风场下海冰漂移规律的模拟进一步验证了本文方法的精确性。在此基础上,对渤海辽东湾的海冰动力过程进行了48h数值分析,计算结果与卫星遥感资料和油气作业区的海冰现场监测数据吻合良好。在下一步工作中将考虑海冰离散元模拟中的热力因素影响,发展具有冻结、断裂效应的海冰离散元模型,更精确地模拟海冰动力-热力耦合作用下的生消和运移过程。     

Surface wave simulation during winter with sea ice in the Bohai Sea

Sea ice can attenuate wave energy significantly when waves propagate through ice covers.In this study,a third-generation wave model called simulating wave nearshore(SWAN)was advanced to include damping of wave energy due to friction in the boundary layer below the ice.With the addition of an eddy viscosity wave-ice model,the resulting new SWAN model was applied to simulate wave height in the Bohai Sea during the freezing winter.Its performance was validated with available buoy data near the ice edge,and the new model showed an improvement in accuracy because it considered the ice effect on waves.We then performed a wave hindcast for the Bohai Sea during a freezing period in the winter of 2016 that had the severest ice conditions in recent years and found that the mean significant wave height changed by approximately 16.52%.In the Liaodong Bay,where sea ice concentration is highest,the change reached 32.57%,compared with the most recent SWAN model version.The average influence of sea ice on wave height simulation was also evaluated over a five-year(2013-2017)hindcast during January and February.We found that the wave height decrease was more significant in storm conditions even the eddy viscosity wave-ice model itself showed no advantage on damping stronger waves.     

Observations and modeling of the ice-ocean conditions in the coastal Chukchi and Beaufort Seas

The Chukchi and Beaufort Seas include several important hydrological features： inflow of the Pacific water, Alaska coast current （ ACC ）, the seasonal to perennial sea ice cover, and landfast ice ＇along the Alaskan coast. The dynamics of this coupled ice-ocean system is important for both regional scale oceanography and large-scale global climate change research. A mumber of moorings were deployed in the area by JAMSTEC since 1992, and the data revealed highly variable characteristics of the hydrological environment. A regional high-resolution coupled ice-ocean model of the Chukchi and Beaufort Seas was established to simulate the ice-ocean environment and unique seasonal landfast ice in the coastal Beaufort Sea. The model results reproduced the Beaufort gyre and the ACC. The depthaveraged annual mean ocean currents along the Beaufort Sea coast and shelf hreak compared well with data from four moored ADCPs, but the simulated velocity had smaller standard deviations, which indicate small-scale eddies were frequent in the region. The model resuits captured the sea,real variations of sea ice area as compared with remote sensing data, and the simulated sea ice velocity showed an ahnost stationary area along the Beaufort Sea coast that was similar to the observed landfast ice extent. It is the combined effects of the weak oceanic current near the coast, a prevailing wind with an onshore component, the opposite direction of the ocean current, and the blocking hy the coastline that make the Beaufort Sea coastal areas prone to the formation of landfast ice.     

海冰动力学过程的数值模拟

讨论了海冰动力学性质并阐述决定海冰漂移的动量平衡,冰脊和水道形成及确定冰应力与形变、强度之间关系的海冰流变学.提出了模拟海冰动力学过程的数值模式,模式中冰厚分布由开阔水、平整冰和堆积冰3种要素表示.在这3要素的预报方程中引入形变函数,采用一种参数化方法模拟冰脊和水道.为了表示冰内相互作用,将海冰作为一种非线性粘性可压缩物质,采用粘-塑性本构关系.本文还概述和讨论了模式中所采用的数值方法,应用此模式模拟了渤海、波罗的海的波的尼亚湾和拉布拉多海的冰漂移.渤海冰漂移模拟结果明显地显示出潮周期变化,还模拟了渤海的冰脊和水道,进行了海冰流变学参数的敏感性试验.并将此冰模式与大气模式和边界层模式联接,给出渤海海冰预报结果.     

An elastic–viscous–plastic sea ice model formulated on Arakawa B and C grids

The elastic–viscous–plastic (EVP) sea ice rheology has been introduced in the large-scale Louvain-la-Neuve sea-Ice Model, version 2 (LIM2), and its performance has been evaluated. Centred difference versions of the rheology have been implemented on both an Arakawa B grid and a C grid, and their performance have been intercompared in coupled simulations with the Nucleus for European Modelling of the Ocean (NEMO) model. Integrations with both implementations lead to fairly similar results which compare well with observations and with previous LIM simulations. The C grid version, however, offers a number of advantages: (a) easier ice coupling with NEMO, which is itself defined on a C grid; (b) possibility of representing ice transport across one-cell-wide straits and passages; (c) better representation of inertial-plastic compressive waves. For these reasons, we recommend the use of the C grid EVP formulation of the ice dynamics in future LIM applications.     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.     

Improvement of short-term sea ice forecast in the southern Okhotsk Sea

In this study, a numerical model of 7-day forecast of sea ice produced by the Japan Meteorological Agency was improved by the following approaches. First, a new ice dynamic model was introduced: the distributed mass/discrete floe model. The model takes account of discrete characteristics of ice floes and well simulates the ice edge location at low computational cost. Secondly, the grid size was reduced to 5 × 5 km for the future high resolution forecasts. Next, the sea surface current data was examined because it significantly influences sea ice movement. We applied two new datasets of HINO and Okhotsk Ocean General Circulation Model (Okhotsk OGCM), which are estimated by numerical simulations, for the 7-day forecast of sea ice. Ice southward speed in January and the whorl formations in February and March were well reproduced with Okhotsk OGCM datasets. Finally, the ocean heat flux at the ice-ocean interface was refined. As a result, we achieved an ice edge error reduction from 30.8 km to 23.5 km.     

Amodified discrete element model for sea ice dynamics

Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.     

基于广义接触模型的离散元方法及其对船体冰载荷的分析

冰载荷是影响船舶结构安全的主要环境载荷,因而确定破冰船的冰载荷特征可用于船体结构设计及船舶安全预警分析。海冰作为典型的脆性材料,可采用单元间具有黏结失效效应的离散元方法模拟其破碎行为。为有效地模拟海冰破碎的渐进过程,采用由局部黏结点组成黏结单元的广义接触模型。为分析该模型中直接决定黏结单元失效的强度参数的影响,采用海冰单轴压缩和三点弯曲试验数值模拟进行分析。通过模型试验确定了失效模型中法向黏结强度、切向黏结强度和Mohr-Coulomb准则中的接触摩擦系数的影响。将该数值模型应用于极地科学考察船"雪龙号"的平整冰区冰载荷分析中,并与船舶整体冰阻力经验公式进行了对比验证。     

预测细粒酒精模型冰物理和力学参数的实验关系

依据模型冰物理性质控制其力学性质的原则 ,利用能量平衡观点建立了冰内未冻结液体含量同制冰过程中的冰池温度、降温及升温时间、模型冰厚度 (喷雾时间×流量 )的理论关系 ,以其量纲表达式作为预测模型冰物理和力学参数的综合控制指标 ;通过系统测量酒精溶液冰点温度、冰池大厅气温、模型冰温度、冰内未冻液含量、冰密度、冰压缩强度、冰弯曲强度和冰变形模量 ,获得了不同冰点温度模型冰的物理和力学参数同该综合控制指标的实验关系 ;实现了用人为可控制因子定量控制模型冰物理和力学参数的目的     

Improved sea-ice radiative processes in a global coupled climate model

1Introduction Seaiceplaysanimportantroleinmoderating heatandmoistureexchangesbetweentheatmosphere andtheoceanathighlatitudes.Seaicealsointeracts withthebroaderclimatesystembythepositiveice albedofeedback(Curryetal.,1995),whichamplifies projectedclimatewarmingatthehighlatitudes,andby theoceanicfeedbackinvolvingicegrowthandmelt, whichinfluencesglobalthermohalinecirculation(i.e., theNorthAtlanticDeepWaterandtheAntarcticBot- tomWater)(Walsh,1983;Barryetal.,1993). Recently,theimplementationofas…     

Modeling regional sediment transport and shoreline response in the vicinity of tidal inlets on the Long Island coast,United States

A new numerical model was developed to simulate regional sediment transport and shoreline response in the vicinity of tidal inlets based on the one-line theory combined with the reservoir analogy approach for volumetric evolution of inlet shoals. Sand bypassing onshore and sheltering effects on wave action from the inlet bar and shoals were taken into account. The model was applied to unique field data from the south coast of Long Island, United States, including inlet opening and closure. The simulation area extended from Montauk Point to Fire Island Inlet, including Shinnecock and Moriches Inlets. A 20-year long time series of hindcast wave data at three stations along the coast were used as input data to the model. The capacity of the inlet shoals and bars to store sand was estimated based on measured cross-sectional areas of the inlets as well as on comprehensive bathymetric surveys of the areas around the inlet. Several types of sediment sources and sinks were represented, including beach fills, groin systems, jetty blocking, inlet bypassing, and flood shoal and ebb shoal feeding. The model simulations were validated against annual net longshore transport rates reported in the literature, measured shorelines, and recorded sediment volumes in the flood and ebb shoal complexes. Overall, the model simulations were in good agreement with the measured data.     

普里兹湾海冰季节性变化的高分辨率数值模拟

普里兹湾海冰以一年冰为主,海冰覆盖存在较大的季节性变化.海冰的分布及其季节性变化主要受当地大气环流及海流的影响.基于一个海洋-海冰耦和模式,模拟研究了该海区海冰的季节性变化特征.海洋模式基于MIT环流模式(MITgcm),海冰动力学模式参考Hibler类型的VP模型,热力学过程取自Winton三层模型.模式区域覆盖整个...     

细粒酒精模型冰的弯曲强度试验研究

基于不同制冰条件和工的细粒酒精模型冰层原位悬臂弯曲强度试验,给出细料酒精模型冰弯曲强度;分析了考虑和不考虑不浮力效应对弯曲强度的差异;利用向睡向下弯曲强度结果给出该模型冰各向同性的证据;得到该冰弯曲强度与其密度、冰内未冻液含量和制冷过程可控制因子的关系,并作为定量控制模型冰质量的基础。     

粘塑性海冰流变模型的稳定性分析

对Hibler 1977提出的具有粘塑性本构关系的海冰模型进行了修改和稳定性分析.在已有的工作中,应用这样一个粘塑性本构关系进行数值积分时大多都必须引入人工扩散项.这里我们用数值逼近技巧提出计算相邻浮冰之间相互作用力的新方法,并且证明修改后的模型无论是海冰的收敛还是发散流动都是稳定的,模式积分无需引入人工扩散项.以二维...