• Journal of the Society of Naval Architects of Korea
• /
• v.57 no.1
• /
• pp.23-34
• /
• 2020

This paper covers the development of prediction techniques for ice load on ice-breakers operating in continuous ice-breaking under level ice conditions using particle-based continuum mechanics. Ice is assumed to be a linear elastic material until the fracture occurs. The maximum normal stress theory is used for the criterion of fracture. The location of the crack can be expressed using a local scalar function consisting of the gradient of the first principal stress and the corresponding eigen-vector. This expression is used to determine the relative position of particle pair to the new crack. The Hertz contact model is introduced to consider the collisions between ice fragments and the collisions between hull and ice fragments. In order to verify the developed technique, the simulation results for the three-point bending problems of ice-specimen and the continuous ice-breaking problem around a wedge-type model ship with bow angle of 20° are compared with the experimental results carrying out at Korea Research Institute of Ships and Ocean Engineering (KRISO).

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

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

• Journal of Nutrition and Health
• /
• v.46 no.4
• /
• pp.382-390
• /
• 2013

This study was designed to investigate the amount of free sugar according to each beverage category in coffee and beverage stores. The groups were categorized as 15 groups based on the kind of beverage material. The beverage groups contributing to total sugar per 100 mL were milk + syrup or powder, hot (12.9 g), ade (12.6 g), milk + syrup or powder + crushed ice (11.9 g), and espresso shot + milk + syrup + crushed ice (11.4 g). The beverage groups contributing to free sugar per 100 mL were ade (12.6 g), milk + syrup or powder + crushed ice (10.8 g), espresso shot + milk + syrup + crushed ice (10.3 g), and milk + syrup or powder, hot (9.7 g). The beverage groups contributing to total sugar (energy) per portion size were milk + syrup or powder + crushed ice 56.6 g (332.3 kcal), espresso shot + milk + syrup + crushed ice 49.3 g (333.4 kcal), milk + syrup or powder, hot 46.3 g (372.1 kcal), and milk + syrup or powder, ice 38.1 g (325.9 kcal). The beverage groups contributing to free sugar per portion size were milk + syrup or powder + crushed ice 51.2 g, espresso shot + milk + syrup + crushed ice 44.9 g, ade 37.1 g, milk + syrup or powder, hot 34.6 g, and milk + syrup or powder, ice 30.1 g. The percent of average free sugar per portion size of the WHO recommendation (free sugars <10% of total energy; <50 g/2,000 kcal) was milk + syrup or powder + crushed ice 102.4%, espresso shot + milk + syrup + crushed ice 89.8%, ade 74.1%, and milk + syrup or powder, hot 69.2%. The proportion of beverage in excess of WHO recommendation per portion size was 14.6% in espresso shot + milk + syrup + crushed ice, 22.7% in ade, and 10.9% in milk + syrup or powder, hot. Therefore, in coffee and beverage stores, menu development with reduced sugar content is needed, and nutrition information should be provided through sugar nutrition labeling.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.16-22
• /
• 2012

In this study, we setup a ground-based scatterometer using an antenna of which the center frequency is 9.5 GHz (X-band), and measured radar backscatterings from air/ice and ice/water interfaces to extract ice thickness. Both of air/ice and ice/water interfaces make strong radar backscatterings and so we can clearly identify two peaks in measured data by scatterometer. By using the distance of two peaks and refractive index of ice, we confirmed that it is possible to measure ice thickness. Field survey was performed at the downstream of Jiam River flowing into Chuncheon Lake. We measured radar backscattering from river ice along a survey path and extracted ice thickness. The ice thickness map of the downstream of Jiam River was produced by using kriging which is one of well known interpolation methods. The ice thickness was about 50 cm along the mainstream while ice was thin as 30 ~ 40 cm at a fast-flowing meander. Ice thickness was particularly thinner at some locations than that of surrounding areas even in the mainstream region of constant flow. This was because of impurities in ice or artificially formed refrozen holes after fishing. We expect that this study helps to expand utilization field of X-band SAR and airborne scatterometer system.

• Atmosphere
• /
• v.31 no.1
• /
• pp.113-141
• /
• 2021

The impacts of ice clouds on the energy budget of the Earth and their representation in climate models have been identified as important and unsolved problems. Ice clouds consist almost exclusively of non-spherical ice crystals with various shapes and sizes. To determine the influences of ice clouds on solar and infrared radiation as required for remote sensing retrievals and numerical models, knowledge of scattering and microphysical properties of ice crystals is required. A conventional method for representing the radiative properties of ice clouds in satellite retrieval algorithms and numerical models is to combine measured microphysical properties of ice crystals from field campaigns and pre-calculated single-scattering libraries of different shapes and sizes of ice crystals, which depend heavily on microphysical and scattering properties of ice crystals. However, large discrepancies between theoretical calculations and observations of the radiative properties of ice clouds have been reported. Electron microscopy images of ice crystals grown in laboratories and captured by balloons show varying degrees of complex morphologies in sub-micron (e.g., surface roughness) and super-micron (e.g., inhomogeneous internal and external structures) scales that may cause these discrepancies. In this study, the current idealized models representing morphologies of ice crystals and the corresponding numerical methods (e.g., geometric optics, discrete dipole approximation, T-matrix, etc.) to calculate the single-scattering properties of ice crystals are reviewed. Current problems and difficulties in the calculations of the single-scattering properties of atmospheric ice crystals are addressed in terms of cloud microphysics. Future directions to develop physically consistent ice-crystal models are also discussed.

• Atmosphere
• /
• v.23 no.1
• /
• pp.39-45
• /
• 2013

A parameterization for the scattering of longwave radiation by ice clouds has been developed based on spectral scattering property calculations with shapes and sizes of ice crystals. For this parameterization, the size distribution data by Fu (1996) and by Michell and Arnott (1994) are used. The shapes of ice crystal considered in this study are plate, solid column, hollow column, bullet-rosette, droxtal, aggregate, and spheroid. The properties of longwave scattering by ice crystals are presented as a function of the extinction coefficient, single-scattering albedo, and asymmetry factor. The heating rate and flux by the radiative parameterization model are calculated for wide range of ice crystal sizes, shapes, and optical thickness. The results are compared with the calculated results using a six-stream discrete ordinate scattering algorithm and Chou's method. The new method (with various habits and size distributions) provides a good simulation of the scattering properties and cooling rate in optically thin clouds (optical thickness < 5). Depending on the inclusion of scattering by ice clouds, the errors in the calculation of the cooling rates are significantly different.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.1
• /
• pp.43-50
• /
• 2018

Recently, various efforts have been made to develop oil and gas in the Arctic Ocean. It is very important to consider the load caused by ice in designing floating structures in the area. The magnitude of the ice load and its impact on a structure should be considered. In this paper, we analyze ice parameters affecting the DP performance of FPSO with a DP-assisted mooring system. Several ice characteristics are selected, and the resulting ice load is calculated using GEM software. Numerous simulations are conducted while changing the values of the parameters, and DP capability plots are generated to visualize the effects of changing these parameters. It is shown that the ice drift speed and thickness are the major properties to be considered in DP system design. The limitations of the analysis and future work are discussed in the conclusion.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.2 no.2
• /
• pp.59-66
• /
• 1990

Various types of ice distribution under low temperature greatly influence the environment of the Arctic and Antarctic Oceans. To understand fundamentals of ice properties such as Polar glaciers, icebergs and sea ice, this study focuses on the material behaviors and failure mechanisms of polycrystalline ice. Utilizing the continuum damage theory, a three-dimensional constitutive model to describe creep deformation characteristics in the glacial flow is developed in consideration of micro-cracking as the major physical process of ice deformation. The numerical model is compared with the published experimental data especially in uniaxial constant stress creep tests. The model can simulate primary and secondary creeps as well as tertiary creep characteristics due to the microcrack accumulation.

• Korean Journal of Remote Sensing
• /
• v.29 no.1
• /
• pp.1-9
• /
• 2013

To classify the landfast ice in the north of the Greenland, observation data, multifrequency Synthetic Aperture Radar (SAR) images and texture images were used. The total four types of sea ice are first year ice, highly deformed ice, ridge and moderately deformed ice. The texture images that were processed by K-means algorithm showed higher accuracy than the ones that were processed by SAR images; however, overall accuracy of maximum likelihood algorithm using texture images did not show the highest accuracy all the time. It turned out that when using K-means algorithm, the accuracy of the multi SAR images were higher than the single SAR image. When using the maximum likelihood algorithm, the results of single and multi SAR images are differ from each other, therefore, maximum likelihood algorithm method should be used properly.