• Journal of the Korea Computer Graphics Society
• /
• v.26 no.5
• /
• pp.1-13
• /
• 2020

In this paper, we propose a new method to stable simulation the directional ice shape by coupling of freezing solver and viscous water flow. The proposed ice modeling framework considers viscous fluid flow in the direction of ice growth, which is important in freezing simulation. The water simulation solution uses the method of applying a new viscous technique to the IISPH(Implicit incompressible SPH) simulation, and the ice direction and the glaze effect use the proposed anisotropic freezing solution. The condition in which water particles change state to ice particles is calculated as a function of humidity and new energy with water flow. Humidity approximates a virtual water film on the surface of the object, and fluid flow is incorporated into our anisotropic freezing solution to guide the growth direction of ice. As a result, the results of the glaze and directional freezing simulations are shown stably according to the flow direction of viscous water.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.8
• /
• pp.747-757
• /
• 2010

The ice accreted on the airfoil is one of the critical drivers that causes the degradation of aerodynamic performance as well as aircraft accidents. Hence, an efficient numerical code to predict the accreted ice shape is crucial for the successful design of de-icing and anti-icing devices. To this end, a numerical code has been developed for the prediction of glaze ice accretion shape on 2D airfoil. Constant Source-Doublet method is used for the purpose of computational efficiency and heat transfer in the icing process is accounted for by Messinger model. The computational results are thoroughly compared against available experiments and other computation codes such as LEWICE and TRAJICE. The direction and thickness of ice horn are shown to yield similar results compared to the experiments and other codes. In addition, the effects of various parameters - temperature, free-stream velocity, liquid water contents, and droplet diameter - on the ice shape are systematically analyzed through parametric studies.

• 한국전산유체공학회:학술대회논문집
• /
• 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 the Korean Society for Aeronautical & Space Sciences
• /
• v.46 no.8
• /
• pp.611-620
• /
• 2018

For the past decades, the analytic model for distributed surface roughness has been developed to improve the accuracy of the icing simulation code. However, it remains limitations to validate the developed model and determine the empirical parameters due to the absence of the quantitative experimental data which were focused on the surface state. To this end, the experimental study conducted to analyze the ice covered surface state from a micro-perspective. Above all, the tendency of the smooth zone width which occurs near the stagnation point has been quantitatively analyzed. It is observed that the smooth zone width is increased as growing the ambient temperature and freestream velocity. Next, the characteristics of the ice covered surface under rime and glaze ice have been analyzed. For rime ice conditions, ice elements are developed as the opaque circular corn in the opposite direction of freestream. The height and interval of each circular corn are increased as rising the ambient temperature. For glaze ice conditions, numerous lumps of translucent ice can be observed. This is because the beads formed by gravity concentrate and froze on the lower surface.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.11
• /
• pp.1105-1114
• /
• 2010

Ice shapes accumulated on the aircraft surfaces are categorized into rime and glaze ice, which are highly dependent on various parameters such as ambient temperature, liquid water contents (LWC), mean volumetric droplet diameter and freestream velocity. In this study, quantitative analyses on the ice accretion have been attempted in a systematical manner and the key findings are as follows. First, the increase of freestream velocity can cause tremendous change in the ice accumulation such as the growth of ice accretion area, ice heading direction and maximum thickness of ice horn. Second, LWC is found to be linearly proportional to the ice accretion area. Third, the effects of ambient temperature on incoming water mass seem to be relatively small in comparison with LWC and freestream velocity. Finally, it was shown that MVD has only a little influence on ice shapes. However, it may increase the ice accretion area by increasing the droplet impacting range.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.2
• /
• pp.352-364
• /
• 2017

A novel 3D mathematical model for water film runback and icing on a rotating surface is established in this work, where both inertial forces caused by the rotation and shear forces due to the air flow are taken into account. The mathematical model of the water film runback and energy conservation of phase transition process is established, with a cyclical average method applied to simulate the unsteady parameters variation at angles of attack. Ice accretion on a conical spinner surface is simulated and the results are compared with the experimental data to validate the presented model. Then Ice accretion on a cowling surface is numerically investigated. Results show that a higher temperature would correspond to a larger runback ice area and thinner ice layer for glaze ice. Rotation would enhance the icing process, while it would not significantly affect the droplet collection efficiency for an axi-symmetric surface. In the case at angle of attack, the effect of rotation on ice shape is appreciable, ice would present a symmetric shape, while in a stationary case the shape is asymmetric.

• Korean Journal of Dental Materials
• /
• v.44 no.3
• /
• pp.197-206
• /
• 2017

The hardening effect by ice-quenching after oxidation of a Pd-Ag-Sn-Au metal-ceramic alloy during porcelain firing simulation was investigated by means of hardness test, field emission scanning electron microscopic observations, and X-ray diffraction analysis. The hardness decreased by ice-quenching after oxidation, which was induced by the homogenization of the ice-quenched specimen. The decreased hardness by ice-quenching after oxidation was recovered from the wash stage which was the first stage of the remaining firing process for bonding porcelain. After wash stage, the hardness of the ice-quenched specimens decreased during the subsequent porcelain firing process. But the final hardness of the ice-quenched specimens after oxidation was higher than that of the specimens cooled at stage 0 after oxidation. The increase in hardness of the specimens during the first firing process was caused by the lattice strains generated at the interface between the face-centered cubic Pd-Ag-rich matrix and the face-centered tetragonal Pd3(Sn, Ga, In) precipitate. The decrease in hardness of the specimens during the remaining firing process was caused by the microstructural coarsening.