• Journal of the Society of Naval Architects of Korea
• /
• v.52 no.3
• /
• pp.180-186
• /
• 2015

This study proposed a new method of preparing granular ethylene glycol/aliphatic detergent (EG/AD) model ice that has both strength and uniform thickness. Various sheets of granular model ice prepared in ice tanks are surveyed and their preparation procedures are analyzed. We not only made a new granular model ice using the EG/AD solution but also measured its thickness, strength, and density. In addition, we found that the strength of the model ice could be controlled by varying the time and air temperature in the consolidation phase. Based on the results of this study, we verified that granular EG/AD model ice can be prepared more uniformly and effectively than columnar EG/AD model ice. This study is intended to contribute to reducing the time required for the ice model test and the operation of the ice model basin.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.19 no.8
• /
• pp.593-598
• /
• 2007

The present study has dealt with the measuring method of ice slurry viscosity using falling sphere viscometer. The experimental apparatus was composed by test section and high-speed video system. And the spheres used in this study were alumina and glass. The main parameters were ice packing factor (IPF) and falling velocity of sphere so the acquired results were discussed for these parameters. The viscosity of ice slurry was calculated by using measured falling velocity and moving distance at instantaneous time and the Stokes hypothesis was used for this calculation. It was clarified that possible measuring range was $IPF\;=\;0.06{\sim}0.14$ of this type of measuring device and measuring method. In addition, it was clarified that the viscosity of ice slurry increased to increase of ice packing factor (IPF) of ice slurry.

• 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.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.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.11
• /
• pp.43-55
• /
• 2018

AGF (Artificial Ground Freezing) method is a temporary ground improvement method which can apply to all types of soil with the purpose of high stiffness and low hydraulic conductivity. However, the groundwater flow and the heterogeneity of the ground increase the uncertainty of the ice-column formation which hinders the reliability of this method. The effects of groundwater flow and layered heterogeneity on ice-wall integrity by AGF method were analyzed using finite element analysis program for a coupled thermo-hydro phenomena in the freezing ground. Groundwater flow changes circular ice-column into elliptical shapes and increases the time required for the formation of ice walls. The previous theoretical formula overestimated the completion time of the ice wall and the critical groundwater velocity by neglecting the thermal interaction between adjacent ice-columns. Numerical results presented the corrected formula and verified the proposed equation for the dimensionless ice-wall completion time. In the layered heterogeneous ground, the thickness of the layer with higher hydraulic conductivity and its relative magnitude were found to be important factors in the ice-wall completion time and critical velocity.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.5
• /
• pp.330-338
• /
• 2021

This paper focuses on the improvement of the influence coefficient matrix method for estimation of local ice load on the icebreaking research vessel ARAON. The influence coefficient matrix relates ice pressure on the hull plate to the measured/calculated hull strain/stress. Conventionally von Mises equivalent stresses representing hull stresses and ice pressure acting on the hull plate are utilized to assemble the influence coefficient matrix. Because of the three dimensional features of the ship-ice collision process, an enhanced method to assemble the influence coefficient matrix is derived considering ice loads in the X, Y, and Z direction simultaneously. Furthermore the location of ice loads acting on hull-plate may fall outside the measuring sensor area, and the enhanced influence coefficient matrix is modified to reduce the difference between the actual and the estimated ice loads by expanding the domain outward from the sensor area. The developed method for enhanced influence coefficient matrix is applied to IBRV ARAON during the 2019 Antarctic ice field test and the local ice loads in three directions are efficiently calculated compared to those by a conventional method.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.4
• /
• pp.451-462
• /
• 1996

In this study, which focuses on ice storage, a fundamental study in cooling and solidification was performed, including the interesting phenomena of density inversion, supercooling and dendritic ice. A numerical study was performed for natural convection and ice formation considering existence of subcooling and dendritic ice were analyzed numerically by using finite difference method and boundary fixing method. In the mesh, the solid fraction was introduced with adding as a term to the energy conservation equation. A flow in the dendrite was modelled as a flow in a porous medium, and the momentum conservation equation was modified to incorporate resistance forces involved in flows through porous media. A numerical solution of the time dependencies of dendrite area and dense ice front was successfully obtained, and the numerical results were good agreement with experimental results. Based on this methodology, a discussion was made of phenomena and characteristics of cooling and freezing processes under various conditions.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.2
• /
• pp.117-127
• /
• 2019

Vessel's ice speed performances will be verified in ice sea trial but environmental conditions of ice fields are changeable according to the weather condition of ice trial area. Speed performance has to correct in the no wind, wave and current etc. after sea trial. Especially finding ice fields which is exact the same as owner's ice thickness and strength requirements is not easy. Therefore speed correction according to environment condition has to be done after sea trial measurements. Correction methods for ice thickness, ice strength, wave, wind and ship draft, trim, ice drift etc. are checked in ice sea trial based on literature review such as ISO standard, ITTC recommendation, journal papers and proceedings of conferences. Possibility of application for current and ice drift correction in ice field are discussed and measuring schemes and procedures of correction methods are described in this paper. All of correction schemes are calculated for 'Araon' which is ice breaking research vessel with Arctic and Antarctic ice field test results. Analyzed results shows that Araon is satisfied with her official ice speed performance of 3 knots with 10MW power at 1m ice thickness, 570kPa ice flexural strength.

• 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.5 no.1
• /
• pp.116-131
• /
• 2013

The resistance performance of an icebreaking cargo vessel in pack ice conditions was investigated numerically and experimentally using a recently developed finite element (FE) model and model tests. A comparison between numerical analysis and experimental results with synthetic ice in a standard towing tank was carried out. The comparison extended to results with refrigerated ice to examine the feasibility of using synthetic ice. Two experiments using two different ice materials gave a reasonable agreement. Ship-ice interaction loads are numerically calculated based on the fluid structure interaction (FSI) method using the commercial FE package LS-DYNA. Test results from model testing with synthetic ice at the Pusan National University towing tank, and with refrigerated ice at the National Research Council's (NRC) ice tank, are used to validate and benchmark the numerical simulations. The designed ice-going cargo vessel is used as a target ship for three concentrations (90%, 80%, and 60%) of pack ice conditions. Ice was modeled as a rigid body but the ice density was the same as that in the experiments. The numerical challenge is to evaluate hydrodynamic loads on the ship's hull; this is difficult because LS-DYNA is an explicit FE solver and the FSI value is calculated using a penalty method. Comparisons between numerical and experimental results are shown, and our main conclusions are given.