Simulations of collision of ice particles

The objective of this paper is to develop a realistic model for ice–structure interaction. To this end, the experiments made by Bridges et al. [Bridges FG, Hatzes A, Liu DNC. Structure, stability and evolution of Saturn’s rings. Nature 1984;309:333–5] in order to measure the coefficient of restitution for ice particles are thoroughly analyzed. One particularly troublesome aspect of the aforementioned experiments is fracture of the ice particles during a collision. In the present effort, the collisional properties of the ice particles are investigated using a Finite Element approach. It is found that a major challenge in modeling collision of the ice balls is the prediction of the onset of fracture and crack propagation in them. In simulations of a block of ice collision to a structure, it is crucial that fracture is determined correctly, as it will influence the collisional properties of the ice particles. The results of the simulation, considering fracture criterion implemented into the Finite Element Model [Zamankhan P, Bordbar M-H. Complex flow dynamics in dense granular flows. Part I: experimentation. J Appl Mech (T-ASME) 2006;73:648–57; Zamankhan P, Huang J. Complex flow dynamics in dense granular flows. Part II: simulations. J Appl Mech (T-ASME) 2007;74:691–702] together with a material model for the ice, imply that most of the kinetic energy dissipation occurs as a result of fracturing at the contact surface of the ice particles. The results obtained in the present study suggest that constitutive models such as those proposed by Brilliantov et al. [Brilliantov NV, Spahn F, Hertzsch JM, Poschel T. Model for collisions in granular gases. Phys Rev E;1996;53:5382–92] for collisions of ice particles are highly questionable.     

Class clustering destroys delay differentiation in priority queues

This paper considers a discrete-time priority queueing model with one server and two types (classes) of customers. Class-1 customers have absolute (service) priority over class-2 customers. New customer batches enter the system at the rate of one batch per slot, according to a general independent arrival process, i.e., the batch sizes (total numbers of arrivals) during consecutive time slots are i.i.d. random variables with arbitrary distribution. All customers entering the system during the same time slot (i.e., belonging to the same arrival batch) are of the same type, but customer types may change from slot to slot, i.e., from batch to batch. Specifically, the types of consecutive customer batches are correlated in a Markovian way, i.e., the probability that any batch of customers has type 1 or 2, respectively, depends on the type of the previous customer batch that has entered the system. Such an arrival model allows to vary not only the relative loads of both customer types in the arrival stream, but also the amount of correlation between the types of consecutive arrival batches. The results reveal that the amount of delay differentiation between the two customer classes that can be achieved by the priority mechanism strongly depends on the amount of such interclass correlation (or, class clustering) in the arrival stream. We believe that this phenomenon has been largely overlooked in the priority-scheduling literature.     

A moving boundary problem derived from heat and water transfer processes in frozen and thawed soils and its numerical simulation

The seasonal change in depths of the frozen and thawed soils within their active layer is reduced to a moving boundary problem,which describes the dynamics of the total ice content using an independent mass balance equation and treats the soil frost/thaw depths as moving(sharp)interfaces governed by some Stefan-type moving boundary conditions,and hence simultaneously describes the liquid water and solid ice states as well as the positions of the frost/thaw depths in soil.An adaptive mesh method for the moving boundary problem is adopted to solve the relevant equations and to determine frost/thaw depths,water content and temperature distribution.A series of sensitivity experiments by the numerical model under the periodic sinusoidal upper boundary condition for temperature are conducted to validate the model,and to investigate the effiects of the model soil thickness,ground surface temperature,annual amplitude of ground surface temperature and thermal conductivity on frost/thaw depths and soil temperature.The simulated frost/thaw depths by the model with a periodical change of the upper boundary condition have the same period as that of the upper boundary condition,which shows that it can simulate the frost/thaw depths reasonably for a periodical forcing.     

A doubly optimal stopping of a risk process (Poster Presentation)

Problems which lead to an optimal stopping of a risk process are considered. Let an insurance company be endowed with an initial capital a > 0, receive insurance premiums and pay out successive claims. The losses occur according to renewal process. At any moment the company may broaden or narrow down the offer, what entails the change of the parameters. These changes concern the rate of income, the intensity of renewal process and the distribution of claims. The model of the risk process with two types of claims stream is considered. After the change the management wants to know the moment of the maximal value of the capital assets. Our goal is to find two optimal stopping times: the best moment of change the parameters and the moment of maximal value of the capital assets. A dynamic programming method to calculate the expected capital at that times is used. Based on the model which combine two types of risk the model of reinsurance with two firms is formulated. In this case the aim is to find for the firms the equilibrium strategy. The equilibrium is constructed in class of strategies driven by their risks. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The choice of constitutive relations for an ice cover

The constitutive relations between the internal stresses and the deformation parameters of a sea ice cover, which are used in the AIDJEX elastoplastic model and Hibler's non-linearly viscous model, are investigated. It is shown that the structural instability of the ice cover with respect to plastic shear deformations is a consequence of the associated flow rule used in these models. The use of constitutive relations which violate the associated flow rule, but which are in good agreement with the physical properties of granular media, is suggested. An ice cover damage parameter and an empirical equation which describes the change in this parameter are introduced into the treatment. Energy relations are investigated.     

Water wave scattering by two sharp discontinuities in the surface boundary conditions

^** Corresponding author. Email: biren{at}isical.ac.in The problem of water wave scattering by two sharp discontinuitiesin the surface boundary conditions involving infinitely deepwater is examined here by reducing it to two coupled Carleman-typesingular integral equations. The discontinuities arise due tothe presence of two types of non-interacting materials floatingon the surface, one type being in the form of an infinite stripof finite width sandwiched between another type. The non-interactingmaterials form an inertial surface which is a mass-loading modelof floating ice and is regarded as a material of uniform surfacedensity having no elastic property. The two integral equationsare solved approximately by assuming the two discontinuitiesto be widely separated, and approximate analytical expressionsfor the reflection and transmission coefficients are also obtained.This problem has applications in wave propagation through stripsof frazil or pancake ice modelled as floating inertial surfaces.Numerical results for the reflection coefficient are depictedgraphically against the wave number for different values ofthe surface densities of the two types of floating materials.The main feature of the graphs is the oscillatory nature ofthe reflection coefficient and occurrence of zero reflectionfor an increasing sequence of discrete values of the wave number.A direct analytical treatment to solve the integral equationsnumerically, when the separation length between the two discontinuitiesis arbitrary, is also indicated. For the case of more than twodiscontinuities the solution methodology of the correspondingscattering problem is described briefly.     

黄河小型清沟封冻过程现场模拟实验与实测数据分析

利用黄河封冻阶段的冰层,开16m2敞流水面模拟清沟,通过雷达测冰厚的方法,观测敞流水面因冰花封冻过程,并测量敞流水面封冻后冰盖厚度分布,并用费氏台观察冰晶体结构.通过将实测冰厚数据绘成等高线图以及冰厚剖面图,证实封河过程中,由周边向中间封冻的清沟除了与冰冻度日有关外,还与清沟下水流速度和方向以及冰花含量有关.水流携带冰花流至清沟时,冰花会漂浮在清沟中.由于水的粘性作用,水流使清沟中的水起动,然后冲击下游冰壁引起涡旋.涡旋场在清沟上游形成较平静区域,可以积累大量冰花;而在下游冰壁处水体运动活跃,冰花无法堆积.经过长时间冻结,上游处冰花被冻结成为较厚的冰花冰,而下游逐渐冻结成较薄的冰花冰.之后,柱状冰在冰花冰下开始生长.最后经过1个月的冻结,形成了敞流水面下游处冰较薄,上游处冰较厚的非对称冻结厚度分布,并在与水流垂直的方向表现为对称冻结厚度分布.     

The diffraction of plane gravitational waves by the edge of an ice cover

The diffraction of plane surface gravitational waves by the edge of an ice cover lying on the surface of an incompressible fluid of infinitely great depth is considered. The ice cover is simulated by a thin elastic plate. The wave reflection and transmission coefficients are determined when it interacts with the ice cover. A wave field is constructed in the fluid under the conditions that a periodic lumped force and a lumped moment act on the edge of the ice cover. It is shown that as the incident wavelength increases the reflection coefficient tends to zero and the transmission coefficient tends to unity.     

Oscillations of a Stratified Liquid Partially Covered with Ice (General Case)

We study the problem of small motions of an ideal stratified liquid whose free surface consists of three regions: liquid surface without ice, a region of elastic ice, and a region of crumbled ice. The elastic ice is modeled by an elastic plate. The crumbled ice is understood as weighty particles of some matter floating on the free surface. Using the method of orthogonal projection of boundary conditions on a moving surface and the introduction of auxiliary problems, we reduce the original initial boundary value problem to an equivalent Cauchy problem for a second-order differential equation in a Hilbert space. We obtain conditions under which there exists a strong (with respect to time) solution of the initial boundary value problem describing the evolution of the hydrodynamic system under consideration.

     

Hybrid Grid-Characteristic Schemes for Arctic Seismic Problems

Doklady Mathematics - The process of conducting seismic exploration of gas fields with a source and receivers installed on the surface of an ice island is considered. A model that includes an ice...     

Exact and approximate expressions for resistance coefficients of a colloidal sphere approaching a solid surface at intermediate Reynolds numbers

A mathematical model has been developed to describe the force of liquid flow acting on a colloidal spherical particle as it approaches a solid surface at intermediate-Reynolds-number-flow regime. The model has incorporated bispherical coordinates to determine a stream function for the flow disturbed by the sphere. The stream function was then used to derive the flow force on the particle as a function of the inter-surface separation distance. The force equation was related to the modified Stokes equation to obtain an exact analytical expression for the correction factor to the Stokes law. Finally, a rational approximation is presented, which is in good agreement with the exact numerical result, and can be readily applied to more general particle–surface interactions involving short-range hydrodynamics associated with colloidal particles in the near vicinity of a large solid collector surface at intermediate Reynolds number of the supporting flow.     

Resonance interactions of waves in an ice channel

The properties of low-amplitude surface waves propagating in an ice channel are investigated in the shallow-water approximation. The ice cover is modelled either by a rigid cap or by a thin elastic plate floating on a liquid surface. It is shown that an ice channel is a waveguide for surface waves. The dispersive properties of the natural oscillations of the liquid in the channel are investigated. The resonance velocities of the motion of the load on the channel surface, at which the amplitude of the forced oscillations of the liquid increases without limit in time, are determined. The decay instability of the natural oscillations of high harmonics with respect to waves of the first mode is demonstrated. The process is described by the standard equations for non-linear three-wave interaction. The investigations lead to the conclusion that critical modes of motion of a boat are realizable in an ice channel.     

海冰弯曲破坏下的破坏准则数值模拟方法研究

海冰的弯曲破坏进程的研究已经称为冰工程界的一项重要课题.通过理论分析结合数值计算对海冰弯曲破坏数值模拟方法进行探索:探讨适合于海冰弯曲应力状态下的本构模型和破坏准则;在LS-DYNA中模拟圆环形冰排弯曲破坏和冰排与锥体结构相互作用的进程,提出了一种对海冰破坏准则验证的技术思路;通过所提出的技术思路对基于三轴压缩试验的Derradji破坏准则进行了修正,使其能够适用于海冰的弯曲破坏.     

Modelling of macroscopical transient frost processes and microscopical driving forces

The artificial saturation phenomenon due to freeze–thaw cycles is described by a multi–phase and multiscale model [1,2,3] formulated within the Theory of Porous Media, [4]. It represents partially saturated concrete as a mixture of 5 interacting constituents φ^α, namely the solid skeleton φ^s, the bulk water φ^l, the pore volume occupied by vapour φ^v, the ice φⁱ and the gel water phase φ^p. Most relevant for the model is the distinction between two length scales and their characteristic time scales. The boundary is marked where macroscopic bulk conditions change to surface physics and chemistry. Surface physics and chemistry acting on the nano–scale affect fundamental properties of concrete and consequently the durability of concrete against freeze–thaw. At the macroscopic scale the model describes transient conditions (i.e. water–uptake, heat transport, volume dilatation of 9%, phase change of first order considering hysteresis) which are characterized by a relatively long time period to reach equilibrium in contrast to the processes modelled on the microstructure. At the microscopic scale the model represents the nanoscopic CSH–gel system consisting of solid CSH and water as a linked system of both components basing on the concept of the “Solid–Liquid Gel System” [5]. In the constribution the numerical results of the model are presented with focus on the evaluation of the process zone during the penetration of the melting front into the matrix. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Ice Accretion on Cables of Various Cross-sections

Poots & Rodgers (1976) gave analytical and numerical resultson the icing of a cable having circular cross-section. The windcarrying freezing fog or supercooled cloud droplets impingeson the cable and a spike of rime ice grows outward into thewind. A linear Stokes law of resistance was assumed in calculatingthe drag force on individual droplets. Consequently the theoreticalpredictions on ice accretion are shown to be restricted to smallcable radii and high wind speeds (i.e. for large Stokes number). The purpose of this paper is to give detailed numerical solutionsfor ice accretion on cables of various cross-sections. Empiricaldrag formulae for droplets in motion, with local Reynolds numbersup to 10², are now employed in the numerical solution of droplettrajectories and impaction velocities on the ice accretion surface.Cables having circular, elliptical, square, and equilateraltriangular cross-section are investigated; also included areresults for a cylinder of elliptic cross-section the major axisof which is aligned at an angle to the upstream wind direction.For fixed wind speed and cable size and cross-section the massof ice accreted is given as a function of the time. These resultsare relevant to the evolution of ice accretion on structures.     

K-means聚类算法在SIR传染病模型中的应用研究

基于SIR传染病模型,建立了具有K-means聚类算法的SIR元胞自动机模拟模型.通过对分别服从高斯分布和随机均匀分布的两类初始感染源的分析与模拟,给出了疾病感染半径与隔离半径对疾病传播的影响.结果显示:在两种不同类型的初试分布下,感染者的最大值分别与疾病感染传播半径和隔离半径呈正相关与负相关关系,感染者数量随时间的变...     

Mathematical Model of Ice Melting on Transmission Lines

During ice storms, ice forms on high voltage electrical lines. This ice formation often results in downed lines and has been responsible for considerable damage to life and property as was evidenced in the catastrophic ice storm of Quebec recently. There are two main aspects, viz., the formation of ice and its timely mitigation. In this paper, we mathematically model the melting of ice due to a higher current applied to the transmission wire. The two dimensional cross-section contains four layers consisting of the transmission wire, water due to melting of ice, ice, and the atmosphere. The model includes heat equations for the various regions with suitable boundary conditions. Heat propagation and ice melting are expressed as a Stefan-like problem for the moving boundary between the layers of ice and water. The model takes into account gravity which leads to downward motion of ice and to forced convection of heat in the water layer. In this paper, the results are applied to the case when the cross-sections are concentric circles to yield melting times for ice dependent on the increase in intensity of the electrical flow in the line. This research has been supported in part by Manitoba Hydro and NSERC.     

流面上的Navier-Stokes方程及其维数分裂方法

该文首先提出了流面和流层的概念,然后推导出了半测地坐标系下流层内的三维NS (Navier-Stokes)方程,以及流面上的二维NS方程.通过引入流面上的流函数,得到了流函数方程的非线性初边值问题,并讨论了方程解的存在性和唯一性.基于以上讨论,提出了求解三维NS方程的维数分裂方法, 并给出了算例.