• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.3
• /
• pp.217-230
• /
• 2019

In this study, a fatigue damage estimation for an ice-going vessel navigating through broken ice fields was carried out. A numerical model to simulate the interaction between ice and structure developed using the finite element method was introduced. Time series of stresses calculated by the proposed model and the corresponding fatigue analysis results are presented. The numerical model enables the long time analysis through an efficient interaction model, the application of the periodic media analysis and the convolution integral, and it allows the stress time history to be extracted directly using the finite element method. To describe the probability distribution of stress amplitudes, the 2-parameter Weibull model was applied to the calculated stress time history, and the fatigue damage was calculated using the Palmgren-Miner rule. Finally, the fatigue damage considering the ice conditions of the Baltic Sea was calculated using the proposed method and LR method, and the results were compared to each other.

• Ocean Systems Engineering
• /
• v.6 no.4
• /
• pp.345-362
• /
• 2016

This study aims to develop the numerical method to estimate level ice impact load and investigate the dynamic interaction between an arctic Spar with sloped surface and drifting level ice. When the level ice approaches the downward sloped structure, the interaction can be decomposed into three sequential phases: the breaking phase, when ice contacts the structure and is bent by bending moment; the rotating phase, when the broken ice is submerged and rotated underneath the structure; and the sliding phase, when the submerged broken ice becomes parallel to the sloping surface causing buoyancy-induced fictional forces. In each phase, the analytical formulas are constructed to account for the relevant physics and the results are compared to other existing methods or standards. The time-dependent ice load is coupled with hull-riser-mooring coupled dynamic analysis program. Then, the fully coupled program is applied to a moored arctic Spar with sloped surface with drifting level ice. The occurrence of dynamic resonance between ice load and spar motion causing large mooring tension is demonstrated.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.6
• /
• pp.556-563
• /
• 2019

The ice keel gouge and seabed interaction is one of the major considerations in the design of an Arctic pipeline system. Ice keel and seabed interaction engineering models based on experimental data, which give an explicit equation for estimating the ice gouging depth, have been suggested. The suggested equations usually overestimate the ice keel gouging depth. In addition, various types of numerical analyses have been carried out to verify the suggested engineering model equations in comparison to the experimental data. However, most of numerical analysis results were also overestimated compared with the laboratory experimental data. In this study, a numerical analysis considering the contact condition and geostatic stress was carried out to predict the ice keel gouging depth and compared with the previous studies. Considering the previously mentioned conditions, more accurate results were produced compared with the laboratory experiment results and the error rate was reduced compared to previous numerical analysis studies.

• Tunnel and Underground Space
• /
• v.16 no.1 s.60
• /
• pp.38-49
• /
• 2006

LNG storage in lined rock cavern demands various techniques concerned with rock mechanics, thermo-mechanics and hydrogeology in design, construction and maintenance stage. LNG pilot cavern was constructed in Daejeon in order to verify these techniques. In this paper, evaluation of drainage system and ice ring formation was studied by numerical simulation. By Modflow analysis in the viewpoint of aquifer and Seep/W analysis in the viewpoint of flow system, it was verified that the drainage system in the pilot cavern was efficiently operated. Since ice ring formation can be simulated by interactive relation between heat transfer and water flow, coupled analysis of those was performed. In this analysis, the position of ice ring was presumed and it was demonstrated that the formation is affected by velocity and direction of groundwater flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.456-459
• /
• 2010

Icing is one of the most serious hazards for aircraft. The amount and rate of icing depend on a number of meteorogical and aerodynamic factors. Of primary importance are amount of liquid water content of droplets, their size, the temperature of aircraft surfaces, the collection efficiency, and the extent of supercooled droplets. In this study, in-flight icing analysis of low reynolds number high aspect ratio wing is carried out by using FENSAP-ICE. Each liquid water contents with altitude is obtained from FAR 25 Appendix-C. And the collectoin efficiency is calculated to check out the ice accretion position of wing with two angles of attack. The degradation of aerodynamic characteristics of aircraft are figured out by investigating the accretion of rime and glaze ice.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.411-418
• /
• 2018

This paper introduces a new method of ice load generation in the time domain for the station-keeping performance evaluation of Arctic offshore structures. This method is based on the ice load spectrum and mean ice load. Recently, there has been increasing interest in Arctic offshore technology for the exploration and exploitation of the Arctic region because of the better accessibility to the Arctic ocean provided by the global warming effect. It is essential to consider the ice load during the development of an Arctic offshore structure. In particular, when designing a station-keeping system for an Arctic offshore structure, a consideration of the ice load acting on the vessel in the time domain is essential to ensure its safety and security. Several methods have been developed to consider the ice load in the time domain. However, most of the developed methods are computationally heavy because they consider every ice floe in the sea ice field to calculate the ice load acting on the vessel. In this study, a new approach to generate the ice load in the time domain with computational efficiency was suggested, and its feasibility was examined. The ice load spectrum and mean ice load were acquired from a numerical analysis with GPU-event mechanics (GEM) software, and the ice load with the varying heading of a vessel was reconstructed to show the feasibility of the proposed method.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.3
• /
• pp.208-215
• /
• 2011

A full-scale field trial in ice-covered sea is one of the most important tasks in the design of icebreaking ships. The first Korean icebreaking research vessel 'ARAON', after her delivery in late 2009, had a sea ice field trial in the Arctic Sea during July-August, 2010. This paper describes the test procedures and data analysis on the icebreaking performance of the IBRV ARAON. The data gathered from the icebreaking performance test in the Chukchi Sea and the Beaufort Sea during the Arctic voyage of ARAON includes the speed and engine power of the ship as well as sea ice thickness and strength data. The air temperature, wind speed and heading of the ship were also measured during each sea ice trial. The ARAON was designed to break 1 m thick level ice with a flexural strength of 630kPa at a continuous speed of 3knots. She is registered as a KR POLAR 10 class ship. The principal dimensions of ARAON are 110 m, 19 m and 6.8 m in length, breadth and draft respectively. She is equipped with four 3,500kW diesel-electric main engines and two Azipod type propulsion motors. Four sea ice trials were carried out to understand the relationship between the engine power and the ship speed, given the Arctic ice condition. The analysis shows that the ARAON was able to operate at 1.5knots in a 2.5m thick medium ice floe condition with the engine power of 5MW, and the speed reached 3.1 knots at the same ice floe condition when the power increased to 6.6MW. She showed a good performance of speed in medium ice floe compared to the speed performance in level ice. More detailed analysis is summarized in this paper.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.5
• /
• pp.562-568
• /
• 2008

In this study, the target FSICR class is 1A whose target thickness of the brash ice is 46 mm in model scale. Normally ice floes for brash ice do not exceed 2 m in full scale, so the model ice sheet was cut by about 10 cm by 10 cm using hand saws. Since the target thickness of brash ice is 46 mm, 46 mm ice sheet makes one layer brash ice. For 23 mm thickness ice sheet, two layers should be accumulated to reach 46mm brash ice thickness. For 15mm thickness ice sheet, three layers need to be accumulated as the same as those in 23 mm ice sheet. New methodology to produce a brash ice was proposed. The results showed that it would be important to use multi-layer rather than single layer possibly because of significant thrust deduction from the propeller-ice interaction in the present ice condition (FSICR 1A).

• 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.57 no.6
• /
• pp.353-364
• /
• 2020

This paper presents the semi-analytical ice load calculation methods that are useful to simulate the ice-breaking process. Since the semi-analytical methods rely on the previously developed closed form equations or numerical analysis results, the user's exact understanding for the equations must be supported in order to use the methods properly. In this study, various failure modes of ice such as local crushing, in-plane splitting failure, out-of-plane bending failure and radial or circumferential cracking with rotation of the broken ice floe are considered. Based on the presented methods, the fracture modes were evaluated according to the size and thickness of ice. In addition, time series analysis for the ice-breaking process was performed on several ice conditions and the results were analyzed.