• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.2
• /
• pp.209-218
• /
• 2017

Sea ice is the main factor affecting the safety of the Arctic engineering. However, traditional numerical methods derived from classical continuum mechanics have difficulties in resolving discontinuous problems like ice damage. In this paper, a non-local, meshfree numerical method called "peridynamics", which is based on integral form, was applied to simulate the interaction between level ice and a cylindrical, vertical, rigid structure at different velocities. Ice in the simulation was freshwater ice and simplified as elastic-brittle material with a linear elastic constitutive model and critical equivalent strain criterion for material failure in state-based peridynamics. The ice forces obtained from peridynamic simulation are in the same order as experimental data. Numerical visualization shows advantages of applying peridynamics on ice damage. To study the repetitive nature of ice force, damage zone lengths of crushing failure were computed and conclude that damage zone lengths are 0.15-0.2 times as ice thickness.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.2
• /
• pp.119-128
• /
• 2018

A numerical method is proposed to simulate level ice interaction with ship in transverse and longitudinal directions in time domain. A novel method is proposed to simulate non-symmetric transverse force in a stochastic way. On the basis of observations from the model tests, the simulation of longitudinal force combines the ice bending force acting on the waterline, submersion force below the waterline and ice friction forces caused by transverse force and ice floes rotation amidships. In the simulations the ship was fixed and towed through an intact ice sheet at a certain speed. The setup of the numerical simulation is similar to the ice tank setup as much as possible. The simulated results are compared with model tests data and the results show good agreement with the measurement.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.414-427
• /
• 2020

The finite element method is used to simulate the navigation of an ice-area bulk carrier in broken ice fields. The ice material is defined as elastic, and the simulations are accomplished at four model speeds and three ice concentrations. The movements of ice floes in the simulation are consistent with those in the model test, and the percentage deviation of the numerical ice resistance from the ice resistance in the model test can be controlled to be less than 15 %. The key characteristics of ice loads, including the average ice loads, extreme ice loads, and characteristic frequency, are analyzed thoroughly in a comprehensive manner. Moreover, the effects of sailing speed and ice concentration on the ice loads are analyzed. In particular, the stress distribution of ice floes is presented to help understand how model speed and concentration affect the ice loads. The "ice pressure" phenomenon is observed at 90 % ice concentration, and it is realistically reflected both in the time―and frequency―domain ice force curves.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.6
• /
• pp.815-823
• /
• 1999

A dynamic simulation program may be used for the development of effective control algorithms for the ice storage cooling system. Simplified effective dynamic models for an ice-on-coil type storage tank, a screw chiller, a water-to-air heat exchanger, three way valves, pipes, pumps, temperature sensors, and controllers were developed. And a dynamic simulation program for the ice storage cooling system was developed by using these dynamic models. Control algorithms for the full storage system were also selected and analyzed in order to show the effectiveness of these models. From the simulation results, it may be concluded that the simulation program developed in this study can be effectively used for the development of improved control algorithms for the ice storage cooling system.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.667-679
• /
• 2020

As the LNG carrier operates in ice covered waters, it is key to ensure the overall safety, which is related to the coupling effect of ice-breaking process and internal liquid sloshing. This paper focuses on the sloshing simulation of the ice-breaking LNG carrier, and the numerical method is proposed using Circumferential Crack Method (CCM) and Volume of Vluid (VOF) with two main key factors (velocity νx and force Fx). The ship motion analysis is carried out by CCM when the ship navigates in the ice-covered waters with a constant propulsion power. The velocity νx is gained, which is the initial excitation condition for the calculation of internal sloshing force Fx. Then, the ship motion is modified based on iterative computations under the union action of ice-breaking force and liquid sloshing load. The sloshing simulation under the LNG tank is studied with the modified ship motion. Moreover, an ice-breaking LNG ship with three-leaf type tank is used for case study. The internal LNG sloshing is simulated with three different liquid heights, including free surface shape and sloshing pressure distribution at a given moment, pressure curves at monitoring points on the bulkhead. This present method is effective to solve the sloshing simulation during ice-breaking process, which could be a good reference for the design of the polar ice-breaking LNG carrier.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.639-647
• /
• 2019

A realistic numerical simulation technology using a Lagrangian Fluid-Structure Interaction (FSI) model was combined with a fracture algorithm to predict the fluid-ice-structure interaction. The failure of ice was modeled as the tensile fracture of elastic material by applying a novel FSI model based on the Moving Particle Semi-implicit (MPS) method. To verify the developed fracture algorithm, a series of numerical simulations for 3-point bending tests with an ice beam were performed and compared with the experiments carried out in an ice room. For application of the developed FSI model, a dropping water droplet hitting a cantilever ice beam was simulated with and without the fracture algorithm. The simulation showed that the effects of fracture which can occur in the process of a FSI simulation can be studied.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.294-306
• /
• 2019

On the basis of the Computational Fluid Dynamics technique (CFD) combined with the overlap grid method, this paper establishes a numerical simulation method to study the problem of ice-propeller interaction in viscous flow and carries out a simulation forecast of the hydrodynamic performance of an ice-class propeller and flow characteristics when in the proximity of milling-shape ice (i.e., an ice block with a groove cut by a high-speed revolving propeller). We use a trimmed mesh in the entire calculation domain and use the overlap grid method to transfer information between the domains of propeller rotation calculation and ice-surface computing. The grid is refined in the narrow gap between the ice and propeller to ensure the accuracy of the flow field. Comparison with the results of the experiment reveals that the error of the hydrodynamic performance is within 5%. This confirms the feasibility of the calculation method. In this paper, we calculate the exciting force of the propeller, analyze the time domain of the exciting force, and obtain the curve of the frequency domain using a Fourier transform of the time-domain curve of the exciting force. The existence of milling-shape ice before the propeller can greatly disturb the wake flow field. Unlike in open water, the propeller bearing capacity shows a downward trend in three stages, and fluctuating pressure is more disordered near the ice.

• Proceedings of the KSRS Conference
• /
• v.2
• /
• pp.892-895
• /
• 2006

The atmospheric effects on the retrieval of sea ice concentration from passive microwave sensors are examined using simulated data typical for the Arctic summer. The simulation includes atmospheric contributions of cloud liquid water and water vapor and surface wind on surface emissivity on the microwave signatures. A plane parallel radiative transfer model is used to compute brightness temperatures at SSM/I frequencies over surfaces that contain open water, first-year (FY) ice and multi-year (MY) ice and their combinations. Synthetic retrievals in this study use the NASA Team (NT) algorithm for the estimation of sea ice concentrations. This study shows that if the satellite sensor’s field of view is filled with only FY ice the retrieval is not much affected by the atmospheric conditions due to the high contrast between emission signals from FY ice surface and the signals from the atmosphere. Pure MY ice concentration is generally underestimated due to the low MY ice surface emissivity that results in the enhancement of emission signals from the atmospheric parameters. Simulation results in marginal ice areas also show that the atmospheric and surface effects tend to degrade the accuracy at low sea ice concentration. FY ice concentration is overestimated and MY ice concentration is underestimated in the presence of atmospheric water and surface wind at low ice concentration. In particular, our results suggest that strong surface wind is more important than atmospheric water in contributing to the retrieval errors of total ice concentrations over marginal ice zones.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.59-69
• /
• 2019

This paper focuses on the numerical simulation of ice abrasion induced by unbreakable ice floe. Under the assumption that unbreakable floes behave as rigid body, the Discrete Element Method (DEM) was applied to simulate the interaction between a fixed structure and ice floes. DEM is a numerical technique which is eligible for computing the motion and effect of a large number of particles. In DEM simulation, individual ice floe was treated as single rigid element which interacts with each other following the given interaction rules. Interactions between the ice floes and structure were defined by soft contact and viscous Coulomb friction laws. To derive the details of the interactions in terms of interaction parameters, the Finite Element Method (FEM) was employed. An abrasion process between a structure and an ice floe was simulated by FEM, and the parameters in DEM such as contact stiffness, contact damping coefficient, etc. were calibrated based on the FEM result. Resultantly, contact length and contact path length, which are the most important factors in ice abrasion prediction, were calculated from both DEM and FEM and compared with each other. The results showed good correspondence between the two results, providing superior numerical efficiency of DEM.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.394-405
• /
• 2020

To assess the station-keeping performance of floating structures in the Arctic region, the ice load should be considered along with other environmental loads induced by waves, wind, and currents. However, present methods for performance evaluation in the time domain are not effective in terms of time and cost. An ice load generation module is proposed based on the experimental data measured at the KRISO ice model basin. The developed module was applied to a time domain simulation. Using the results of a captive model test conducted in multiple directions, the statistical characteristics of ice loads were analyzed and processed so that an ice load corresponding to an arbitrary angle of the structure could be generated. The developed module is connected to commercial dynamic analysis software (OrcaFlex) as an external force input. Station-keeping simulation in the time domain was conducted for the same floating structure used in the model test. The mooring system was modeled and included to reflect the designed operation scenario. Simulation results show the effectiveness of the proposed ice generation module and its application to station-keeping performance evaluation. Considering the generated ice load, the designed structure can maintain a heading angle relative to ice up to 4°. Station-keeping performance is enhanced as the heading angle conforms to the drift direction. It is expected that the developed module will be used as a platform to verify station-keeping algorithms for Arctic floating structures with a dynamic positioning system.