共查询到20条相似文献，搜索用时 787 毫秒
1.
综合考虑欧拉坐标下有限差分法(FDM)在海冰动力学计算中的效率,以及拉格朗日坐标下光滑质点流体动力学方法(SPH)对海冰流变行为的精确模拟,本文发展了一种海冰动力学的混合拉格朗日-欧拉(HLE)数值方法。该方法首先在拉格朗日坐标下将海冰离散为若干个具有厚度、密集度的海冰质点,并由这些海冰质点通过Gauss函数对欧拉网格上的海冰参量进行积分插值;然后,在欧拉坐标下对海冰动量方程进行差分计算以确定各网格节点的海冰速度,并由此采用Gauss函数积分插值出拉格朗日坐标下各海冰质点的速度分布;最后,通过对海冰质点运动和分布的计算,确定出各海冰质点的位置、厚度和密集度等参量。采用该HLE方法对规则区域内的海冰堆积过程和涡动风场作用下的海冰动力演化趋势进行了数值试验;最后,采用该HLE方法对渤海海冰的动力过程进行了72h数值模拟,其计算结果与卫星遥感图像和现场观测资料吻合较好。以上计算结果均表明该HLE方法在海冰动力学数值模拟中具有较高的计算效率和模拟精度,可用于海冰动力过程的数值模拟。 相似文献
2.
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. 相似文献
3.
冰动力学的拉格郎日离散元模式 总被引：5，自引：0，他引：5
本文描述一种应用于河冰和海冰动力学研究的平滑质点流体力学的离散元模式。该离散元方法是一种自由拉格郎日方法,和传统的欧拉网格数值方法相比,它具有较小的数值扩散并且能够更加灵活地处理复杂的冰过程等许多优点。本文首先给出海冰动力学的离散元模式的基本方法,并与欧拉有限差分方法和半拉格郎日的质点一网格方法进行了对比,其次,本文还给出了离散元海冰模式在渤海的一个应用。 相似文献
4.
海冰动力过程的改进离散元模型及在渤海的应用 总被引：1，自引：0，他引：1
海冰的断裂、重叠和堆积等离散分布特性广泛地存在于极区和副极区的不同海域,并对海冰的生消、运移过程有着重要影响。针对海冰在不同尺度下的离散分布特点,发展海冰动力过程的离散元方法有助于完善海冰数值模式,提高海冰数值模拟的计算精度。为此,本文针对海冰生消运移过程中的非连续分布和形变特性,发展了适用于海冰动力过程的改进离散元模型(MDEM)。不同于传统离散元方法,该模型将海冰离散为具有一定厚度、尺寸和密集度的圆盘单元。海冰单元设为诸多浮冰块的集合体,其在运移和相互接触碰撞过程中,依照质量守恒发生单元尺寸、密集度和厚度的相应变化。基于海冰离散性和流变性的特点,该模型采用黏弹性接触本构模型计算单元间的作用力,并依据Mohr-Coulomb准则计算海冰法向作用下的塑性变形及切向摩擦力。为验证该模型的可靠性,本文对海冰在规则水域内的运移和堆积过程进行了分析,离散元计算结果与解析值相一致;此外,对旋转风场下海冰漂移规律的模拟进一步验证了本文方法的精确性。在此基础上,对渤海辽东湾的海冰动力过程进行了48h数值分析,计算结果与卫星遥感资料和油气作业区的海冰现场监测数据吻合良好。在下一步工作中将考虑海冰离散元模拟中的热力因素影响,发展具有冻结、断裂效应的海冰离散元模型,更精确地模拟海冰动力-热力耦合作用下的生消和运移过程。 相似文献
5.
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. 相似文献
6.
基于Mohr-Coulomb准则的黏弹-塑性海冰动力学本构模型 总被引：10，自引：1，他引：10
针对中小尺度下海冰动力作用过程中的漂移、重叠和堆积特征,在连续介质力学基础上建立了一个黏弹塑性海冰动力学本构模型.该模型主要包括四部分,即海冰在屈服前的Kelvin Vogit黏弹模型、海冰塑性屈服的Mohr Coulomb准则、屈服后相关联的正交流动法则和影响海冰强度的静水压力.采用黏弹塑性本构模型对规则海域内的海冰堆积过程进行了数值试验,其结果与传统黏塑性模型、改进的黏塑性模型和经典冰坝理论的解析解相符,可合理地模拟海冰堆积的动力过程.为进一步检验该黏弹塑性本构模型的适用性,对渤海海冰的动力过程进行了48h的数值模拟,计算结果与海冰卫星遥感现场观测资料一致.通过与传统的黏塑性和改进的黏塑性模型结果的对比分析,进一步验证了黏弹塑性模型在中小尺度海冰动力学模拟中的可靠性.在以上海冰堆积的数值试验和渤海海冰动力学模拟中均采用了光滑质点流体动力学(SPH)方法. 相似文献
7.
在寒冷地区,海冰在核电站取水结构物前的堆积会对取水口造成阻塞和损坏,进而影响核电设备的正常运作。以红沿河核电站的取水设施为例,综合考虑该工程海域冰区特点及风和流的作用等因素。建立了海冰的离散元模型,用于模拟海冰在结构物前的堆积和破坏过程。该离散元模型由规则排列的球体颗粒构成,颗粒间采用平行黏结模型进行黏结。考虑了海冰的平均尺寸、密集度及流速三个因素对海冰堆积过程的影响,对海冰堆积高度进行预判。其中,堆积高度随海冰的密集度和流速的增大而增大,而海冰的平均尺寸对堆积高度没有明显影响。分析结果表明,离散元数值模拟可用于评估和预测取水工程中海冰堆积造成堵塞的风险。 相似文献
8.
A viscoelastic-plastic constitutive model with Mohr-Coulomb yielding criterion for sea ice dynamics 总被引：2，自引：0，他引：2
A new viscoelastic-plastic （VEP） constitutive model for sea ice dynamics was developed based on continuum mechanics. This model consists of four components： Kelvin-Vogit viscoelastic model, Mohr-Coulomb yielding criterion, associated normality flow rule for plastic rehololgy, and hydrostatic pressure. The numerical simulations for ice motion in an idealized rectangular basin were made using smoothed particle hydrodynamics （SPH） method, and compared with the analytical solution as well as those based on the modified viscous plastic（VP） model and static ice jam theory. These simulations show that the new VEP model can simulate ice dynamics accurately. The new constitutive model was further applied to simulate ice dynamics of the Bohai Sea and compared with the traditional VP, and modified VP models. The results of the VEP model are compared better with the satellite remote images, and the simulated ice conditions in the JZ20-2 oil platform area were more reasonable. 相似文献
9.
10.
11.
The theoretical foundation of a wave–ice interaction model is reported in Part 1 of this study. The model incorporates attenuation of ocean surface waves by sea ice floes and the concomitant breaking of the floes by waves that determines the structure of the marginal ice zone (MIZ). A numerical implementation of the method is presented here. Convergence of the numerical method is demonstrated, as temporal and spatial grids are refined. A semi-analytical method, which does not require time-stepping, is also developed to validate the numerical results, when dispersion is neglected. The wave energy lost during ice breakage is parameterized, as part of the numerical method. Sensitivity studies are conducted in relation to the energy loss and also dispersive effects, the choice of the attenuation model, the properties of the wave field, and sea ice properties such as concentration, thickness and breaking strain. Example simulations intended to represent conditions in the Fram Strait in 2007, which exploit reanalyzed wave and ice model data, are shown to conclude the results section. These are compared to estimates of MIZ widths based on a concentration criteria, and obtained from remotely-sensed passive microwave images. 相似文献
12.
B. A. Kagan D. A. Romanenkov E. V. Sofina 《Izvestiya Atmospheric and Oceanic Physics》2007,43(6):766-784
The tidal ice drift is treated as an element of the three-dimensional tidal dynamics in a sea covered by ice. This dynamics is described by the QUODDY-4 finite-element model, and the tidal ice drift is described by a continuous viscous-elastic approximation. We present the results of modeling not only the tidal ice drift (M 2 wave) (its velocity, direction, and tidal variations in the concentration and pressure of ice compression) but also ice-induced changes in tidal dynamics and the residual tidal ice drift. The modeling results indicate that the maximum velocity of tidal ice drift, which is determined by a combination of various factors responsible for ice evolution and primarily by the horizontal gradient of the level and local tidal velocity, can be higher or lower than the velocity of the surface tidal current in the ice-free sea. This depends on the sign of deviations of tidal sea level elevations in the sea covered by ice from their values in the ice-free sea. In addition, it has been found that ice cover has a stronger effect on the energetics of tides than on their dynamics: the area-mean relative deviations constitute 1.5% for the density of the total tidal energy, 61.5% for the dissipation, 0.1% for the amplitudes of tidal sea level elevations, and 0.9% for the amplitudes of maximum barotropic tidal velocity. In this sense, the conclusion that the role of sea ice is insignificant in the formation of tides can be justified only partially. The main results of this paper are as follows: (1) the development of a module for tidal ice drift, (2) the inclusion of this module into the three-dimensional finite-element hydrothermodynamic model QUODDY-4 to extend its capabilities, and (3) the reproduction (on the basis of the modified model) of qualitative features of the practically important tidal ice drift and ice-induced changes in the tidal dynamics of marginal seas on the Siberian continental shelf. 相似文献
13.
Tuomo M. Saloranta 《地球，A辑:动力气象学与海洋学》2000,52(1):93-108
A numerical 1‐dimensional fine grid sea ice thermodynamic model is constructed accounting specially for: (1) slush formation via flooding and percolation of rain‐ and snow meltwater, (2) the consequent snow ice formation via slush freezing, and (3) the effects of snow compaction on heat diffusion in snow cover. The model simulations from ice winter period 1979–90 are viewed against corresponding observations at the Kemi fast ice station (65 °39.8' N, 24° 31.4' E). The 11‐year averaged model results show good overall consistency with corresponding total ice thickness observations. The model slightly overestimates the snow ice thickness and underestimates the snow thickness in February and March, which is mainly addressed to the model assumption of isostatic balance (i.e., slush formation via flooding), which was probably not fully satisfied at the coastal Kemi fast ice station. Supposing that this assumption is nevertheless generally valid away from the very coastal fast ice zone, an estimate for sea ice sensitivity to changes in winter precipitation rate is produced. Increased precipitation leads to an increase only in snow ice thickness with little change in total ice thickness, while a reduction in precipitation of more than {213}50% causes a significant increase in total ice thickness. The difference in modeled total ice thickness for the case of artificially neglecting snow ice physics is about 25%, which indicates the importance of including snow ice physics in a sea ice model dealing with the seasonal sea ice zone. 相似文献
14.
Qi Shu Fangli Qiao Jiping Liu Zhenya Song Zhiqiang Chen Jiechen Zhao Xunqiang Yin Yajuan Song 《海洋学报(英文版)》2021,40(10):65-75
To improve the Arctic sea ice forecast skill of the First Institute of Oceanography-Earth System Model (FIO-ESM) climate forecast system, satellite-derived sea ice concentration and sea ice thickness from the Pan-Arctic Ice-Ocean Modeling and Assimilation System (PIOMAS) are assimilated into this system, using the method of localized error subspace transform ensemble Kalman ?lter (LESTKF). Five-year (2014–2018) Arctic sea ice assimilation experiments and a 2-month near-real-time forecast in August 2018 were conducted to study the roles of ice data assimilation. Assimilation experiment results show that ice concentration assimilation can help to get better modeled ice concentration and ice extent. All the biases of ice concentration, ice cover, ice volume, and ice thickness can be reduced dramatically through ice concentration and thickness assimilation. The near-real-time forecast results indicate that ice data assimilation can improve the forecast skill significantly in the FIO-ESM climate forecast system. The forecasted Arctic integrated ice edge error is reduced by around 1/3 by sea ice data assimilation. Compared with the six near-real-time Arctic sea ice forecast results from the subseasonal-to-seasonal (S2S) Prediction Project, FIO-ESM climate forecast system with LESTKF ice data assimilation has relatively high Arctic sea ice forecast skill in 2018 summer sea ice forecast. Since sea ice thickness in the PIOMAS is updated in time, it is a good choice for data assimilation to improve sea ice prediction skills in the near-real-time Arctic sea ice seasonal prediction. 相似文献
15.
16.
17.
渤海悬浮颗粒物的三维输运模式 I.模式 总被引：18，自引：4，他引：14
为了长时间,大范围研究渤海中悬浮物的输运规律,研究应用并改造了备国汉堡大学的粒子追踪悬浮物输运的三维模式,模式中考虑了风、气压等气象要素以外海传入的潮波作用下的三维正压海流对悬浮颗粒物的输运．同时考虑了悬浮颗粒物的学降及再悬浮机制和底质中的细颗粒物的运动,将风浪的作用引入了悬浮物输运的计算,数值模式应用的粒五追踪方法,能较好地反映悬浮物浓度的迅速变化。 相似文献
18.
Numerical simulation for dynamical processes of sea ice 总被引：1，自引：0，他引：1
ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｆｏｒｄｙｎａｍｉｃａｌｐｒｏｃｅｓｓｅｓｏｆｓｅａｉｃｅＷｕＨｕｉｄｉｎｇ，ＢａｉＳｈａｎ，ＺｈａｎｇＺｈａｎｈａｉａｎｄＬｉＧｕｏｑｉｎｇ（ＲｅｃｅｉｖｅｄＭａｙ１６，１９９６；ａｃｃｅｐｔｅｄＪａｎｕａｒｙ... 相似文献
19.
The calculative method of heat transfer coefficient between ice cover and water is analyzed considering the heat balance at ice cover bottom firstly. The heat transfer coefficient is calculated with the meteorological, oceanographic data and sea ice conditions measured on the JZ20-2 Oil/Gas Platform in the Bohai Sea during the winter of 1997/1998. From the results, it is shown that the heat transfer coefficient is smaller in the freezing and melting periods, which is about 0.16× 10-3 and 0.04× 10-3 respectively. In the middle of ice season, the heat transfer coefficient has a larger value, which is about 0.5 × 10-3. Lastly, the influences of ice thickness and ice type on the heat transfer coefficient are discussed. With the heat transfer coefficient determined above, the oceanic heat flux in the winter of 1997～1998 is calculated, and its trend in the winter is analyzed. This study can be referenced in the sea ice numerical simulation and prediction in the Bohai Sea. 相似文献
20.
辽东湾JZ20—2海域海冰参数的概率分布 总被引：11，自引：1，他引：11
基于辽东湾JZ2 0 2海域 1996 2 0 0 0年 4个冬季的海冰定点观测资料和海冰数值模拟结果 ,对该海域的平整冰厚、冰速、冰向和压缩强度等海冰参数进行了概率分析 ,确定了各自的分布参数 ,并对冰速和冰向进行了联合概率分析。结果表明 :冰厚和冰速分别服从对数正态分布和瑞利分布 ,海冰压缩强度服从正态分布。计算结果可用于JZ2 0 2海域海洋结构可靠性设计和疲劳累积损伤分析的海冰参数 ,也可作为其邻近海域的参考资料 相似文献