Influence of snow cover on the parameters flexural-gravity waves in ice cover

This paper presents a mathematical model for the analysis of the stress-strain state of ice cover in the presence of a snow layer under dynamic loads. Numerical calculations were made using the model of indestructible ice, and the results were compared with experimental data.     

Modelling the rheology of sea ice as a collection of diamond-shaped floes

In polar oceans, seawater freezes to form a layer of sea ice of several metres thickness that can cover up to 8% of the Earth’s surface. The modelled sea ice cover state is described by thickness and orientational distribution of interlocking, anisotropic diamond-shaped ice floes delineated by slip lines, as supported by observation. The purpose of this study is to develop a set of equations describing the mean-field sea ice stresses that result from interactions between the ice floes and the evolution of the ice floe orientation, which are simple enough to be incorporated into a climate model. The sea ice stress caused by a deformation of the ice cover is determined by employing an existing kinematic model of ice floe motion, which enables us to calculate the forces acting on the ice floes due to crushing into and sliding past each other, and then by averaging over all possible floe orientations. We describe the orientational floe distribution with a structure tensor and propose an evolution equation for this tensor that accounts for rigid body rotation of the floes, their apparent re-orientation due to new slip line formation, and change of shape of the floes due to freezing and melting. The form of the evolution equation proposed is motivated by laboratory observations of sea ice failure under controlled conditions. Finally, we present simulations of the evolution of sea ice stress and floe orientation for several imposed flow types. Although evidence to test the simulations against is lacking, the simulations seem physically reasonable.     

Modeling ocean wave propagation under sea ice covers

Operational ocean wave models need to work globally, yet current ocean wave models can only treat ice covered regions crudely. The purpose of this paper is to provide a brief overview of ice effects on wave propagation and different research methodology used in studying these effects. Based on its proximity to land or sea, sea ice can be classified as: landfast ice zone, shear zone, and the marginal ice zone. All ice covers attenuate wave energy. Only long swells can penetrate deep into an ice cover. Being closest to open water, wave propagation in the marginal ice zone is the most complex to model. The physical appearance of sea ice in the marginal ice zone varies. Grease ice, pancake ice,brash ice, floe aggregates, and continuous ice sheet may be found in this zone at different times and locations. These types of ice are formed under different thermal-mechanical forcing. There are three classic models that describe wave propagation through an idealized ice cover: mass loading,thin elastic plate, and viscous layer models. From physical arguments we may conjecture that mass loading model is suitable for disjoint aggregates of ice floes much smaller than the wavelength, thin elastic plate model is suitable for a continuous ice sheet, and the viscous layer model is suitable for grease ice. For different sea ice types we may need different wave ice interaction models. A recently proposed viscoelastic model is able to synthesize all three classic models into one. Under suitable limiting conditions it converges to the three previous models. The complete theoretical framework for evaluating wave propagation through various ice covers need to be implemented in the operational ocean wave models. In this review, we introduce the sea ice types, previous wave ice interaction models, wave attenuation mechanisms,the methods to calculate wave reflection and transmission between different ice covers, and the effect of ice floe breaking on shaping the sea ice morphology. Laboratory experiments,field measurements and numerical simulations supporting the fundamental research in wave-ice interaction models are discussed. We conclude with some outlook of future research needs in this field.     

Stress-strain state of ice cover during aircraft takeoff and landing

We consider the linear unsteady motion of an IL-76TD aircraft on ice. Water is treated as an ideal incompressible liquid, and the liquid motion is considered potential. Ice cover is modeled by an initially unstressed uniform isotropic elastic plate, and the load exerted by the aircraft on the ice cover with consideration of the wing lift is modeled by regions of distributed pressure of variable intensity, arranged under the aircraft landing gear. The effect of the thickness and elastic modulus of the ice plate, takeoff and landing regimes on stress-strain state of the ice cover used as a runway.     

On the characterization of the pressure-sinkage relationship of snow covers containing an ice layer

Based on an analysis performed by Meyerhof and on field data obtained in two snow covers, a method for estimating the bearing capacity of an ice layer in a snow cover is proposed. A procedure for characterizing the pressure-sinkage relationship of a snow cover containing an ice layer is suggested. The objective is to provide an improved characterization of the behaviour of snow-covered terrain for the prediction of vehicle mobility.     

海冰动力学数值模拟中改进的PIC方法

为了准确地模拟海冰的动力过程,需要建立精确有效的数值方法。本文结合质点网格法(PIC)和光滑质点流体动力学方法(SPH)发展了一种改进的PIC方法。该方法在欧拉坐标下对海冰动量方程进行差分计算,在拉格朗日坐标下进行海冰质点位移、厚度和密集度计算,并采用Gauss函数进行欧拉网格点与拉格朗日质点间海冰参数的交互插值。采用改进的PIC方法对规则区域内的海冰堆积过程进行了数值试验,对渤海海冰的动力过程进行了72小时数值模拟。计算结果均表明改进的PIC方法具有计算量小,计算结果平稳精确的优点,可很好地适用于海冰动力作用过程的数值模拟。     

Action of periodic surface pressure on an ice cover in the vicinity of a vertical wall

This paper presents the solution of the linear hydroelastic problem for steady forced vibrations of a semi-infinite ice cover under the effect of localized external load. The ice cover is simulated by a viscoelastic thin plate, the thickness of the fluid layer is assumed to be small, and the shallow water theory is used. The fluid is limited by a solid vertical wall, and the rectilinear edge of the elastic plate adjacent to the wall can be both free and clamped. The solution is obtained with the help of the Fourier integral transform. The behavior of the ice cover is studied depending on the frequency of the external load and boundary conditions on the edge of the plate. It is shown that, in the case of a free edge of the plate, there are considerable deflections on the edge, which could be comparable with deflections at the center of the pressure impact region. It is established that, due to the existence of wave movements of the type of edge waves, the external load energy is transferred to larger distances along the free edge, and there are significant bending moments on the edge of the clamped plate, which can lead to fracture of the ice cover with sufficiently great intensity of the external load.     

Solitary wave packets beneath a compressed ice cover

A family of plane solitary wave packets of a small (but finite) amplitude on the surface of an ideal incompressible fluid of finite depth beneath an ice cover is described. The solitary wave trains correspond to solutions of the two-dimensional system of Euler’s equations of an ideal incompressible fluid of the type of a traveling wave which decreases at infinity and has identical phase and group velocities. The ice cover is simulated by an elastic Kirchhoff-Love plate freely floating on the fluid surface in the compressed state.     

Waves generated on an ice cover by a source pulsating in fluid

A point source of variable intensity located at rest in a plane infinitely deep fluid layer under an ice cover is considered. The general expression for perturbations of the fluid-ice interface is obtained. In the case of the long operation in the pulsating regime the wave established on the ice cover is found.     

Oscillations of a cylindrical body submerged in a fluid with ice cover

The linear plane problem of oscillations of an elliptic cylinder in an ideal incompressible fluid of finite depth in the presence of an ice cover of finite length is solved. The ice cover is modeled by an elastic plate of constant thickness. The hydrodynamic loads acting on the body are determined as functions of the oscillation frequency and the positions of the cylinder and plate.     

Behavior of a semi-infinite ice cover under periodic dynamic impact

Oscillations of a semi-infinite ice cover in an ideal incompressible liquid of finite depth under local time-periodic axisymmetric load are considered. The ice cover is simulated by a thin elastic plate of constant thickness. An analytical solution of the problem is obtained using the Wiener–Hopf method. The asymptotic behavior of the amplitudes of oscillations of the plate and the liquid in the far field is studied. It is shown that the propagation of waves in the far field is uneven: in some directions, the waves propagate with a significantly greater amplitude.     

Numerical modelling of ice jam resistance to main channel flow

In northern countries, subfreezing temperatures during the winter season result in the formation of ice covers on most rivers. Towards the end of the winter season, during the spring break-up period, stationary ice covers become weak in strength and break up. The resulting broken ice pieces or ice floes are significantly larger in thickness and have a rougher undersurface relative to sheet ice and impose higher hydraulic resistance. The downstream movement of the ice floes may be arrested under conditions such as an intact ice cover, bridge piers or channel constrictions, among others, thereby initiating a break-up ice jam. These ice jams most often have been observed to cause very high water stages. Detrimental effects caused by these high water levels encompass those of operational and design-related problems such as the flooding of communities due to ice-jam-induced backwater, flood risk assessments, altering of the open water flow regime, bed scour and flooding of bridges. The ability to predict the influence of an ice jam on the main flow is of considerable importance in river engineering and can be viewed upon by its effects on the variation in the water surface levels. All other information is dependent on the foregoing. The ice jam influence on the main flow can be regarded with respect to local and global standpoints. The primary objective of this study is to formulate the influence of the ice jam on the main channel flow. The formulation is then coupled with a two-dimensional numerical model for the simulation of the water flow regime. The data from different laboratory experiments on ice jams are reproduced numerically. Various simulations are then carried out to compute the water surface levels and velocities in channels under ice jam conditions. The numerical results are then compared with the laboratory data. Results show that the mathematical formulation developed to predict the water surface levels and velocities along the ice jam length as well as upstream and downstream of its leading and trailing edges respectively gives satisfactory predictions.     

Interaction of surface and flexural-gravity waves in ice cover with a vertical wall

The problem of the interaction of surface and flexural-gravity waves with a vertical barrier is solved in a two-dimensional formulation. It is assumed that the fluid is ideal and incompressible, has infinite depth, and is partially covered with ice. The ice cover is modeled by an elastic plate of constant thickness. The eigenfrequencies and eigenmodes of oscillation of the floating elastic ice plate, the deflection and deformation of ice, and the forces acting on the wall are determined.     

Dynamic ice loads on conical structures

Two series of model tests were performed to observe the dynamic ice loads on conical structures. The variable testing parameters include the water line diameter of the model cone and ice parameters. During small water line diameter tests, two-time breaking is found to be the typical failure of ice on steep conical structure, and also be controlled by other factors, such as ice speed and the cone angle. During big water line diameter tests, the ice sheet failed nonsimultaneously around the cone. Several independent zones of bending were found in the nonsimultaneous failure process of ice. With the increase of the ratio of D/h and the number of independent zones, the total ice force was found being gradually reduced.     

Water wave propagation over multiple porous barriers with variable porosity in the presence of an ice cover

Meccanica - A semi-analytical method is applied to investigate the propagation of flexural wave over multiple bottom-standing porous barriers with variable porosity beneath an ice cover under the...     

Hydrodynamic loads acting on an oscillating cylinder submerged in a stratified fluid with ice cover

The two-dimensional problem of steady oscillations of a horizontal cylinder submerged in a linearly stratified fluid layer whose upper boundary is ice cover is considered in a linear treatment using the Boussinesq approximation. The method of mass sources distributed along the body contour is used for the internal wave generation regime, and the integral equation for the disturbed pressure in the fluid is used for the regime of no internal waves. The hydrodynamic load acting on the body was calculated as a function of the oscillation frequency for the case of a continuous ice cover and for special cases (broken ice, free surface, and rigid lid).     

水下爆炸冲击波载荷作用下冰层破碎特性及其影响因素

水下爆炸破冰是复杂的爆破工程,为了研究冰层在水下爆炸冲击波载荷作用下的破碎特性及规律,利用几何动力分析软件LS-DYNA对水下爆炸破冰的过程进行数值模拟,并将计算结果与实验结果进行对比,误差在8%以内,验证了数值模型的有效性。根据本文中的建模方法及建立的模型,计算不同的实验工况:实验场地环境不变,调整爆距分析不同爆距下冰层破碎特性;调整药量、爆距和冰厚,通过正交设计方法设计9组实验方案,应用灰色系统理论对3种因素进行分析,建立了各个因素与破冰半径之间的灰色关联度系数及灰色关联度。分析结果表明:药量为100 g,冰厚为29 cm,水深为2.9 m,爆距范围为0.3~1.5 m破冰的半径范围为0~1.1 m,最佳爆距范围在0.3~0.45 m之间;根据以上9种工况的分析可知,药量(100、200、300 g)、爆距(0.3、0.6、0.9 m)和冰厚(24、28、32 cm)对破冰半径的影响的主次关系依次为爆距、药量、冰厚。     

Experimental investigation of ice slurry flow pressure drop in horizontal tubes

Pressure drop behaviour of ice slurry based on ethanol–water mixture in circular horizontal tubes has been experimentally investigated. The secondary fluid was prepared by mixing ethyl alcohol and water to obtain initial alcohol concentration of 10.3% (initial freezing temperature ?4.4 °C). The pressure drop tests were conducted to cover laminar and slightly turbulent flow with ice mass fraction varying from 0% to 30% depending on test conditions. Results from flow tests reveal much higher pressure drop for higher ice concentrations and higher velocities in comparison to the single phase flow. However for ice concentrations of 15% and higher, certain velocity exists at which ice slurry pressure drop is same or even lower than for single phase flow. It seems that higher ice concentration delay flow pattern transition moment (from laminar to turbulent) toward higher velocities. In addition experimental results for pressure drop were compared to the analytical results, based on Poiseulle and Buckingham–Reiner models for laminar flow, Blasius, Darby and Melson, Dodge and Metzner, Steffe and Tomita for turbulent region and general correlation of Kitanovski which is valid for both flow regimes. For laminar flow and low buoyancy numbers Buckingham–Reiner method gives good agreement with experimental results while for turbulent flow best fit is provided with Dodge–Metzner and Tomita methods.Furthermore, for transport purposes it has been shown that ice mass fraction of 20% offers best ratio of ice slurry transport capability and required pumping power.     

Increasing the efficiency of the resonance method for breaking an ice cover with simultaneous movement of two air cushion vehicles

This paper describes a study on the possibility of increasing the efficiency of the resonance method for breaking an ice cover due to the flexural-gravitational wave interference arising during simultaneous movement of several air cushion vehicles.     

Numerical modeling of ice behavior under high velocity impacts

In this work a constitutive relation for ice at high strain rates and an algorithm for its numerical integration are developed. This model is based on the Drucker–Prager plasticity criteria, which allows a different behavior in tension and in compression. In addition a failure criteria, based on pressure cut-offs, is implemented to describe the ice damage. In order to validate the constitutive model, numerical simulations were compared with experimental results, in which ice cylinders were impacted against a steel plate, allowing the measurement of the contact load. Three different numerical solvers are used in order to analyze its performance to appropriately modeling the ice behavior.