• Atmosphere
• /
• v.30 no.1
• /
• pp.75-89
• /
• 2020

Current in-situ airborne probes that measure the sizes of ice crystals smaller than 50 ㎛ are based on the concept that the measured intensity of light scattered by a particle in the forward and/or backward direction can be converted to particle size. The relationship between particle size and scattered light used in forward scattering probes is based on Mie theory, which assumes the refractive index of particle is known and all particles are spherical. Not only are small crystals not spherical, but also there are a wide variety of non-spherical shapes. Although it is well known that the scattering properties of non-spherical ice crystals differ from those of spherical shapes, the impacts of non-sphericity on derived in-situ particle size distributions are unknown. Thus, precise relationships between the intensity of scattered light and particle size and shape are required, as based on accurate calculations of scattering properties of ice crystals. In this study, single-scattering properties of ice crystals smaller than 50 ㎛ are calculated at a wavelength of 0.55 ㎛ using a numerically exact method (i.e., discrete dipole approximation). For these calculations, hexagonal ice crystals with varying aspect ratios are used to represent the shapes of natural small ice crystals to determine the errors caused by non-spherical ice crystals measured by forward scattering probes. It is shown that the calculated errors in sizing nonspherical ice crystals are at least 13% and 26% in forward (4~12°) and backward (168~176°) directions, respectively, and maximum errors are up to 120% and 132%.

• 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.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.2
• /
• pp.796-808
• /
• 2019

In this paper, a Smoothed Particle Hydrodynamics (SPH) method is extended to simulate the ice failure process and ice-ship interactions. The softening elastoplastic model integrating Drucker-Prager yield criterion is embedded into the SPH method to simulate the failure progress of ice. To verify the accuracy of the proposed SPH method, two benchmarks are presented, which include the elastic vibration of a cantilever beam and three-point bending failure of the ice beam. The good agreement between the obtained numerical results and experimental data indicates that the presented SPH method can give the reliable and accurate results for simulating the ice failure progress. On this basis, the extended SPH method is employed to simulate level ice interacting with sloping structure and three-dimensional ice-ship interaction in level ice, and the numerical data is validated through comparing with experimental results of a 1:20 scaled Araon icebreaker model. It is shown the proposed SPH model can satisfactorily predict the ice breaking process and ice breaking resistance on ships in ice-ship interaction.

• Korean Journal of Remote Sensing
• /
• v.34 no.6_2
• /
• pp.1299-1310
• /
• 2018

In this study, we analyzed distribution and movement trends using in-situ observations and particle tracking methods to understand the movement of the drift ice in the Arctic Ocean. The in-situ movement data of the drift ice in the Arctic Ocean used ITP (Ice-Tethered Profiler) provided by NOAA (National Oceanic and Atmospheric Administration) from 2009 to 2018, which was analyzed with the location and speed for each year. Particle tracking simulates the movement of the drift ice using daily current and wind data provided by HYCOM (Hybrid Coordinate Ocean Model) and ECMWF (European Centre for Medium-Range Weather Forecasts, 2009-2017). In order to simulate the movement of the drift ice throughout the Arctic Ocean, ITP data, a field observation data, were used as input to calculate the relationship between the current and wind and follow up the Lagrangian particle tracking. Particle tracking simulations were conducted with two experiments taking into account the effects of current and the combined effects of current and wind, most of which were reproduced in the same way as in-situ observations, given the effects of currents and winds. The movement of the drift ice in the Arctic Ocean was reproduced using a wind-imposed equation, which analyzed the movement of the drift ice in a particular year. In 2010, the Arctic Ocean Index (AOI) was a negative year, with particles clearly moving along the Beaufort Gyre, resulting in relatively large movements in Beaufort Sea. On the other hand, in 2017 AOI was a positive year, with most particles not affected by Gyre, resulting in relatively low speed and distance. Around the pole, the speed of the drift ice is lower in 2017 than 2010. From seasonal characteristics in 2010 and 2017, the movement of the drift ice increase in winter 2010 (0.22 m/s) and decrease to spring 2010 (0.16 m/s). In the case of 2017, the movement is increased in summer (0.22 m/s) and decreased to spring time (0.13 m/s). As a result, the particle tracking method will be appropriate to understand long-term drift ice movement trends by linking them with satellite data in place of limited field observations.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.4
• /
• pp.475-482
• /
• 1998

In this paper, an experimental study was conducted to investigate the performance of the spherical ice particle production system which uses the technology of water spray in a vacuum chamber for increasing the heat transfer area. As a result, following conclusions were obtained. The diffusion-controlled evaporation model agreed relatively well with experiments. The spray flow rate influences the performance of the system rather than any other factors, for example, the spray nozzle position, the nozzle number. As the spray rate increases, the system efficiency increases. It is due to the entrainment of small droplets among the spray with the exhausted vapor. Thus the system should be designed and operated to prevent the entrainment. Assuming the compressor efficiency to be 70%, the COP of the system reaches highly up to 6 at a maximum spray rate. Under the conditions, however, the rigid ice layer is obtained because ice particles bond together with increase of the spray rate. Therefore, the spray rate should be controlled properly to make the spherical ice particles.

• Polymer(Korea)
• /
• v.29 no.3
• /
• pp.282-287
• /
• 2005

A novel ice particle leaching method for fabrication of porous and biodegradable PLGA scaffold has been proposed for the application to tissue engineering. After uniform mixing of poly(L-lactide-co-glycolide) (PLGA) and bovine serum albumin-fluorescein isothiocyanate (FITC-BSA), the FITC-BSA loaded scaffold was fabricated by adding various ratio of ice particle. The release profiles of FITC-BSA were examined using pH 7.4 PBS for 28 days at $37^{circ}$. The release amount was determined by fluorescence intensity by using the fluorescence spectrophotometer and the morphological change of the scaffolds was observed by scanning electron microscope. The release initial burst of BSA containing scaffolds was lower than that of simple dipping scaffolds resulting in constant release aspect. Although the BSA concentration increased. the initial burst was not increased. As a result of this study, it can be suggested that ice particle leaching method for the tissue engineered scaffold miff be very useful and it is possible to impregnate with water soluble factors like cytokine. We suggest that ice particle leaching method may be useful to tissue engineered organ regeneration.

• Atmosphere
• /
• v.22 no.1
• /
• pp.73-85
• /
• 2012

The collision efficiency data for collision between graupel or hail particles and cloud drops that take into account the differences of particle density are applied to the Takahashi cloud model. The original setting assumes that graupel or hail collision efficiency is the same as that of the cloud drops of the same volume. The Takahashi cloud model is run with the new collision efficiency data and the results are compared with those with the original. As an initial condition, a thermodynamic profile that can initiate strong convection is provided. Three different CCN concentration values and therefore three initial cloud drop spectra are prescribed that represent maritime (CCN concentration = 300 $cm^{-3}$), continental (1000 $cm^{-3}$) and extreme continental (5000 $cm^{-3}$) air masses to examine the aerosol effects on cloud and precipitation development. Increase of CCN concentration causes cloud drop sizes to decrease and cloud drop concentrations to increase. However, the concentration of ice particles decreases with the increase of CCN concentration because small drops are difficult to freeze. These general trends are well captured by both model runs (one with the new collision efficiency data and the other with the original) but there are significant differences: with the new data, the development of cloud and raindrop formation are delayed by (1) decrease of ice collision efficiency, (2) decrease of latent heat from riming process and (3) decrease of ice crystals generated by ice multiplication. These results indicate that the model run with the original collision efficiency data overestimates precipitation rates.

• Proceedings of the SAREK Conference
• /
• 2008.06a
• /
• pp.1105-1110
• /
• 2008

An experimental study was performed to control Ice Packing Factor (IPF) of ice slurry in a pipe in a real time. This paper presented the concept that IPF can be adjusted by the amount of the solution contained to ice slurry. Based on this concept, we designed IPF controller consisting of the outlet tube providing ice slurry and the upper tube discharging only a solution through holes, and investigated the technical validity and efficiency of the controller experimentally. As a result, the original proposed IPF controller could not control IPF of ice slurry in a pipe. This is because an ice of ice slurry was drained out into not only the outlet but also the upper of the controller due to the size of the holes relatively large compared to the ice particle. Therefore, we changed the hole size of IPF controller surface using fine meshes and then, observed that IPF in a pipe was increased by $4{\sim}7$ percent when the hole size was $80{\mu}m$ and less.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.21 no.9
• /
• pp.496-504
• /
• 2009

In order to investigate the feasibility of a direct ice slurry transporting system for the purpose of district cooling, a case study of field application is performed. The research aims include the field measurement of ice packing factor, the performance of coldness delivery, and the branching characteristics of ice slurry. Two representative types of pipe branch are dealt with in this work. For the slurry flow with ice volume fraction of 0.16 or less, the pipe blocking due to aggregation is not observed. Based on the time-wise variation of temperature in the storage tank, a calculating method of ice packing factor is newly developed, which seems to be useful when the brine concentration is unknown. It is confirmed that the mass flow rate of ice slurry per unit cooling load is markedly reduced with increasing the ice content. The pumping power also decreases, but remains unchanged for high ice fractions. The distribution of ice particle before and after branching shows a good uniformity within the range of 5% difference, but yields a unique trend depending on the flow rate.