Ice Shape Research Articles

Measuring Ice Pellets and Refrozen Wet Snow Using a Laser-Optical Disdrometer

Abstract This study aims to characterize the shapes and fall speeds of ice pellets formed in various atmospheric conditions and to investigate the possibility to use a laser-optical disdrometer to distinguish between ice pellets and other types of precipitation. To do so, four ice pellet events were documented using manual observations, macrophotography, and laser-optical disdrometer data. First, various ice pellet fall speeds and shapes, including spherical, bulged, fractured, and irregular particles, were associated with distinct atmospheric conditions. A higher fraction of bulged and fractured ice pellets was observed when solid precipitation was completely melted aloft while more irregular particles were observed during partial melting. These characteristics affected the diameter–fall speed relations measured. Second, the measurements of particles’ fall speed and diameter show that ice pellets could be differentiated from rain or freezing rain. Ice pellets larger than 1.5 mm tend to fall > 0.5 m s−1 slower than raindrops of the same size. In addition, the fall speed of a small fraction of ice pellets was < 2 m s−1 regardless of their size, as compared with a fall speed > 3 m s−1 for ice pellets with diameter > 1.5 mm. Video analysis suggests that these slower particles could be ice pellets passing through the laser-optical disdrometer after colliding with the head of the instrument. Overall, these findings contribute to a better understanding of the microphysics of ice pellets and their measurement using a laser-optical disdrometer. Significance Statement Ice pellets are challenging to forecast and to detect automatically. In this study, we documented the fall speed and physical characteristics of ice pellets during various atmospheric conditions using a combination of a laser-optical disdrometer, manual observations, and macrophotography images. Relationships were found between the shape and fall speed of ice pellets. These findings could be used to refine the parameterization of ice pellets in atmospheric models and, consequently, improve the forecast of impactful winter precipitation types such as freezing rain. Furthermore, they will also help to physically interpret laser-optical disdrometer data during ice pellets and freezing rain.

Aerodynamic prediction and roughness implementation toward ice accretion simulation

Ice accretion in aircraft leads to significant separation and surface roughness. Therefore, the simulation of icing is highly intricate and challenging. The aerodynamic prediction of an iced wing requires an accurate turbulence model. In this paper, a modified transition model for the Reynolds-averaged Navier–Stokes equation is developed. The model is implemented using a linear eddy viscosity transition model as a baseline. It is modified to accurately predict the flow separation and simulate the influence of roughness. Different cases are numerically tested based on the modified model, including two types of iced airfoils (glaze ice and rime ice) and two rough geometries. The results indicate that the model demonstrates engineering accuracy in terms of the flow separation and aerodynamic coefficient prediction. Moreover, this model could be adapted for a wide range of ice shapes. In addition, the present model is capable of accounting for the influence of surface roughness, leading to significant improvements in heat transfer coefficients that match the experimental data well.

Statistical Analysis of the Surface Roughness on Aircraft Icing

A statistical analysis of the surface roughness is performed on experimentally obtained ice shapes on an asymmetrical airfoil at Rec≈2⋅106. The ice shapes were generated in the Icing Wind Tunnel of the Technical University of Braunschweig under Appendix C and O conditions of the EASA airplane certification standards as part of the ICE GENESIS project. The photogrammetry method is used for the digitization of the experimental ice shapes, while statistical parameters such as the mean ice shape and the local root mean square (RMS) of the ice geometry are extracted using a traditional surface projection method, as well as a self-organizing maps approach. Results show the evolution of the statistical parameters over time and the influence of the freestream static temperature on these parameters. A comparison between the experimental values of the local RMS of the ice geometry and a correlation for roughness prediction is presented, showing a good match with the original formulation of the correlation for cases under Appendix C conditions while having a good match with Appendix O conditions when a temperature correction factor is applied to the formulation. Additionally, results show an almost linear growth of roughness over the whole accretion time.

Poly(l-alanine-co-l-threonine succinate) as a Biomimetic Cryoprotectant.

We synthesized a series of [(l-Ala)x-co-(l-Thr succinate)y] (PATs), which are analogous to natural antifreezing glycoprotein with the structure of [l-Ala-l-Ala-l-Thr disaccharide]n, by varying the composition and degree of succinylation while fixing their molecular weight (Mn) and Ala/Thr ratio at approximately 10-12 kDa and 2:1, respectively. We investigated their ice recrystallization inhibition (IRI), ice nucleation inhibition (INI), dynamic ice shaping (DIS), thermal hysteresis (TH), and protein cryopreservation activities. Both IRI and INI activities were greater for PATs with higher l-Ala content (PATs-3 and PATs-4) than those with lower l-Ala content (PATs-1 and PATs-2). DIS activity with faceted crystal growth was clearly observed in PATs-2 and PATs-4 with a high degree of succinylation. TH was small with <0.1 °C for all PATs and slightly greater for PATs with a high l-Ala content. Except for PATs-1, the protein (lactate dehydrogenase, LDH) stabilization activity was excellent for all PATs studied, maintaining LDH activity as high as that of fresh LDH even after 15 freeze-thaw cycles. To conclude, the cryo-active biomimetic PATs were synthesized by controlling the l-Ala content and degree of succinylation. Our results showed that PATs with an l-Ala content of 65-70% and degree of succinylation of 12-19% exhibited the cryo-activities of IRI, INI, and DIS, and particularly promising properties for the cryoprotection of LDH protein.

Numerical Investigations of Synthetic Jet Control Effects on Iced Airfoils

Based on the numerical method, the effects of synthetic jet control (SJC) on the aerodynamic characteristics of iced airfoils are systematically analyzed for the first time. First, an analysis method for the airfoil flowfield under SJC in icing conditions is developed, which contains the CFD method, the icing prediction method, and the boundary conditions of the synthetic jet. Then, the variation in the aerodynamic characteristics of different iced airfoils under the SJC are analyzed, and the flow separation near the ice shape is discussed in detail. Finally, parameters, such as the jet position and angle, are quantified, and conclusions are obtained. The synthetic jet can reduce the separation area on the upper surface induced by the ice accretion. When the energy of the jet is sufficient, it interacts with the mainstream, reducing the separation vortex downstream of the jet orifice and greatly improving the aerodynamic characteristics of the iced airfoil. These results provide valuable insights for the application of synthetic jet technology in mitigating the adverse effects of ice accretion on the airfoil’s aerodynamic characteristics.

Experimental Study on Optimum Design of Aircraft Icing Detection Based on Large-Scale Icing Wind Tunnel

Icing detection is the premise and basis for the operation of aircraft icing protection system, and is the primary issue in flight safety assurance. At present, there is a lack of research methods and design reference for the layout optimization of ice detectors. Therefore, in order to simulate the real icing environment encountered by the aircraft more accurately, a large-scale icing wind tunnel was used to carry out experimental research on the icing characteristics of the sensor probes. A closed-loop experimental method including the typical condition selection, sensor array interference examination and ice shape repeatability verification was initially proposed. A stepwise optimization process and a sensitivity analysis on ambient conditions were combined to determine the optimal distribution for sensor installation. It is found that the water collection coefficient on the cylinder surface of the probe first increases and then decreases along the axial direction, reaching the extreme value at a certain position. The height of this extreme point will gradually increase with the development of the wall boundary layer, showing a variation range of 2~30 mm. Improper design may cause the sensor probe to fail to capture the point with the maximum icing thickness, affecting the sensitivity of icing detection. In addition, each probe position has different sensitivity to changes in flow parameters; the points with larger icing mass and lower sensitivity to changes in attack angle will have better detection effect. The measured data and analysis in the present work can provide a basis for the accurate design of icing sensor probes.

Computational Investigation of the Water Droplet Effects on Shapes of Ice on Airfoils

The paper presents the results of studying the effects of droplet diameters on the NACA0012 airfoil ice accretion, which have been obtained in the 3D numerical simulation of icing. To simulate the motion of water droplets as a multiphase medium, the Eulerian approach is used, which assumes that water droplets have spherical shapes, do not undergo deformation and breakup, do not interact with each other, and that coalescence/fragmentation of droplets does not take place. Both monodisperse (of the same size) and polydisperse (of various sizes) droplets are considered; they are represented by the spectral Langmuir distributions. These spectral distributions take into account the polydisperse nature of droplets and provide a higher efficiency in predicting ice shapes. The obtained ice shapes on an airfoil are compared with the available experimental and calculated data. It should be noted according to the simulation results that the use of the standard size of droplet diameter equal to 20 μm does not allow for obtaining correct shapes of ice on the leading edge of the wing profile not at all temperature regimes. For temperatures from −20 °C to −10 °C, there is a noticeable difference compared to the experimental data. At the same time, for this temperature range, the use of the Langmuir spectral distribution of droplet diameters relative to 15 μm provides a better agreement of the formed ice forms with the experiment.

Where to live? Landfast sea ice shapes emperor penguin habitat around Antarctica.

Predicting species survival in the face of climate change requires understanding the drivers that influence their distribution. Emperor penguins (Aptenodytes forsteri) incubate and rear chicks on landfast sea ice, whose extent, dynamics, and quality are expected to vary substantially due to climate change. Until recently, this species' continent-wide observations were scarce, and knowledge on their distribution and habitat limited. Advances in satellite imagery now allow their observation and characterization of habitats across Antarctica at high resolution. Using circumpolar high-resolution satellite images, unique fast ice metrics, and geographic and biological factors, we identified diverse penguin habitats across the continent, with no significant difference between areas with penguins or not. There is a clear geographic partitioning of colonies with respect to their defining habitat characteristics, indicating possible behavioral plasticity among different metapopulations. This coincides with geographic structures found in previous genetic studies. Given projections of quasi-extinction for this species in 2100, this study provides essential information for conservation measures.

A numerical simulation of the distribution and the variation law of the liquid water content in icing wind tunnel

To solve the aircraft icing problem, simulation of the high-altitude clouds as the icing environment becomes necessary. Clouds always contain many micro-scale subcooled droplets. The accurate control of the icing environment simulation is challenging. To investigate the distribution and the variation law of the liquid water content in the icing wind tunnel, a simulation model of the icing wind tunnel was developed with the reprogramming software FLUENT. The model was validated with the parameters of ice shape and maximum ice thickness on the model of SC0710 airfoil with a constant cross-section at zero angle of attack. The errors of the maximum ice thickness, the maximum and mean deviation are 5.06%, 0.22 mm and 0.13 mm, respectively. Then, the influences of the parameters of water droplets at the nozzle exit and airflow were simulated and analyzed. The results indicate that increasing the initial velocity or temperature of water droplets at nozzle exit will promote the uniformity of the liquid water content (LWC) distribution and reduce the average LWC in the test section. The average LWC of the test section is significantly reduced due to a large amount of water droplets hitting the wall under the action of inertia once the medium volume diameter (MVD) exceeds 40 μm.

Numerical Investigation on the influence of meteorological conditions on wind turbine icing

Keeping the blade structure and shape in the design condition is the premise and guarantee for efficient operation of wing turbine, but wet and cold climate will lead to icing blade, then the blade shapes were changed and the output power of the wind turbine is reduced. A numerical frame for wind turbine icing was developed in present work, which consists of a numerical method for impingement of water droplets based on the multiple frames of reference and Euler method and a 3D heat and mass transferring method for wind turbine icing, then the effects of several meteorological parameters, including the temperature, Liquid Water Content (LWC), Mean Volume Diameter (MVD) and the inflow velocity, on ice shape of the blade were investigated in the numerical frame. The present method can provide supports for prediction of icing process and aerodynamic characteristic after icing of wind turbine blades.

A dynamic collection efficiency formulation implemented in the icing code DICEPS

The Technical University of Braunschweig developed the 2D icing code DICEPS, used for the simulation of ice accretion on aircraft components exposed to icing conditions. In its formulations, DICEPS presents an improvement in the traditional formulation which icing codes use for the calculation of the collection efficiency, which works together with a multi-sub-step algorithm used for the surface deformation during an icing-step. This new formulation continuously updates the impingement mass flow based on the change of the wall-normal as the surface deforms during the multi-sub-step routine in an icing-step, resulting in a dynamic collection efficiency. In the current work, DICEPS and its dynamic collection efficiency formulation (DCEF) are presented and validated through several test cases under three categories: prediction of collection efficiency, rime ice shapes and glaze ice shapes. DICEPS managed to successfully replicate distributions of collection efficiency along with an implemented splashing model. A substantial improvement in the prediction of rime ice shapes via single-step simulations was achieved using the DCEF, while the traditional formulation requires several icing-steps to obtain a converged ice shape. This is also seen when rime ice results are compared to those of FENSAP-ICE, which require as well several icing-steps. In the case of glaze ice, the DCEF presents a time limit since prominent ice features will heavily influence the real collection efficiency, making the DCEF inaccurate after this limit. Nevertheless, the effects that geometrical changes have over parameters such as the heat transfer imply that a multi-step approach is eventually required. By analyzing the results, the DCEF shows that it effectively resembles the physical behavior that geometrical changes impose on the collection efficiency for particles with large inertia, while smaller particles present some inaccuracies, suggesting a proper extension of DCEF for these cases.

Ice melting in salty water: layering and non-monotonic dependence on the mean salinity

The presence of salt in seawater strongly affects the melt rate and the shape evolution of ice, both of utmost relevance in ice–ocean interactions and thus for the climate. To get a better quantitative understanding of the physical mechanics at play in ice melting in salty water, we numerically investigate the lateral melting of an ice block in stably stratified saline water. The developing ice shape from our numerical results shows good agreement with the experiments and theory from Huppert &amp; Turner (J. Fluid Mech., vol. 100, 1980, pp. 367–384). Furthermore, we find that the melt rate of ice depends non-monotonically on the mean ambient salinity: it first decreases for increasing salt concentration until a local minimum is attained, and then increases again. This non-monotonic behaviour of the ice melt rate is due to the competition among salinity-driven buoyancy, temperature-driven buoyancy and salinity-induced stratification. We develop a theoretical model based on the force balance which gives a prediction of the salt concentration for which the melt rate is minimal, and is consistent with our data. Our findings give insight into the interplay between phase transitions and double-diffusive convective flows.

A rotorcraft in-flight ice detection framework using computational aeroacoustics and Bayesian neural networks

This work develops a novel ice detection framework specifically suitable for rotorcraft using computational aeroacoustics and Bayesian neural networks. In an offline phase of the work, the acoustic signature of glaze and rime ice shapes on an oscillating wing are computed. In addition, the aerodynamic performance indicators corresponding to the ice shapes are also monitored. These performance indicators include the lift, drag, and moment coefficients. A Bayesian neural network is subsequently trained using projected Stein variational gradient descent to create a mapping from the acoustic signature generated by the iced wings to predict their performance indicators along with quantified uncertainty that is highly important for time- and safety-critical decision-making scenarios. While the training is carried out fully offline, usage of the Bayesian neural network to make predictions can be conducted rapidly online allowing for an ice detection system that can be used in real time and in-flight.

Self-assembly of tamarind seed polysaccharide via enzymatic depolymerization and degalactosylation enhanced ice recrystallization inhibition activity

Polysaccharides are becoming potential candidates for developing food-grade cryoprotectants due to their extensive accessibility and health-promoting effects. However, unremarkable ice recrystallization inhibition (IRI) activity and high viscosity limit their practical applications in some systems. Our previous study found a galactoxyloglucan polysaccharide from tamarind seed (TSP) showing moderate IRI activity. Herein, the enhancement of the IRI performance of TSP via enzymatic depolymerization and degalactosylation-induced self-assembly was reported. TSP was depolymerized and subsequently removed ∼40 % Gal, which induced the formation of supramolecular rod-like fiber self-assembles and exhibited a severalfold enhancement of IRI. Ice shaping assay did not show obvious faceting of ice crystals, indicating that both depolymerized and self-assembled TSP showed very weak binding to ice. Molecular dynamics simulation confirmed the absence of molecular complementarity with ice. Further, it highlighted that degalactosylation did not cause significant changes in local hydration properties of TSP from the view of a single oligomer. The inconsistency between molecular simulation and macroscopic IRI effect proposed that the formation of unique supramolecular self-assemblies may be a key requirement for enhancing IRI activity. The findings of this study provided a new opportunity to enhance the applied potential of natural polysaccharides in food cryoprotection.

A Survey Method for Drift Ice Characteristics of the Yellow River Based on Shore-Based Oblique Images

Acquisition of continuous drift ice characteristic parameters such as ice size, shape, concentration, and drift velocity are of great importance for understanding river freezing and thawing processes. This study acquired hourly oblique images captured by a shore-based camera in the winter of 2021–2022 on the Yellow River, China. The pixel point scale method for correcting oblique images is provided. The 61 lines were measured at the calibration site and the absolute error between the measured value and the calculated value was in the range of 0.009–0.850 m, with a mean error of 0.236 m. After the correction of oblique images, the true equivalent diameter of drift ice during the freezing period ranged from 0.52–13.10 m with a mean size of 3.36 m, which was larger than that of 2.30 m during the thawing period which ranged from 0.20–12.54 m. It was found that the size of drift ice increased with time during the freezing period and decreased with time during the thawing period. The fractal dimension and roundness were used to represent drift ice shape. The fractal dimension ranged from 1.0–1.3 and the roundness ranged from 0.1–1.0. A Gaussian distribution was used to estimate drift ice size and shape distributions. There is a nonlinear relationship between ice concentration and drift velocity, which can be well expressed by the logistic function. In the future, drift ice parameters for more years and hydrometeorological data for the same time need to be accumulated, which helps to analyze the freezing and thawing patterns of river ice.

Centrifugation forced ice-concentrate separation to enhance suspension freezing desalination performance

Suspension freeze desalination is a promising route for treatment of high salinity wastewater, but the poor ice-concentrate separation often generates decrease of desalination efficiency. In this work, it is demonstrated that centrifugation forced dewatering can assist in ice purification and enhance suspension freeze desalination performance. The positive correlation and good linear fitting relationship between ice purity and dewatering rate were identified at various centrifugation speeds, centrifugation time, and ice layer thickness. Changes of ice shape were recorded by low-temperature optical microscope observation system. The free water removal and ice floc destruction induced by centrifugal field have an important contribution to ice purification. Increasing centrifugation speed from 100 r/min to 500 r/min, the corresponding relative centrifugal force increased 27-fold, creating sharp increases in ice purity (rate constant k = 0.1605) and dewatering rate (rate constant D = 0.0563). Due to the strong binding characteristics and small portion of bound water, ice desalination efficiency per unit centrifugal force dropped when centrifugation speed increased from 1000 r/min to 3000 r/min. Besides, the multi-step dewatering and purification mechanism involved in ice centrifugation was depicted. The suspension freeze technology hybrid with centrifugation not only enhances desalination efficiency but also is beneficial to obtain high-purity ice.

Scallop ice shape characteristics of swept wing based on large-scale icing wind tunnel experiment

Scallop ice is a special phenomenon that occurs during swept wing aircraft passing through icing clouds. It poses a great challenge for the icing safety assessment that the complex scallop ice shape feature and its mechanism are still unclear. In this work, a large-scale icing wind tunnel experiment of swept wing designed by NACA0012 airfoil is conducted in the Icing Wind Tunnel of China Aerodynamics Research and Development Center. The detailed three-dimensional ice shapes under 0°, 15°, 30° and 45° swept angles are obtained by laser scanning technology. The experimental results show that with the swept angle increasing from 0° to 45°, the 2D double ice horn structures show certain spanwise variation, and finally transform into complete scallop ice with ice thickness greatly enhanced in the stagnation line region. The empirical mode decomposition of the spanwise ice curve captures the high-frequency fluctuation on the scallop ice caused by the small-scale roughness element, while the trend with low frequency is not obvious. Based on the experimental data, a new complete scallop ice geometric model, named 5Points-5Lines-2Arcs (5P-5L-2A) model, is proposed, which can provide important basis for the quantitative description of complex scallop ice shape.

Experimental and Numerical Analysis of Incompressible Flow over an Iced Airfoil

Determining the aerodynamic characteristics of iced airfoil is an important step in aircraft design. The goal of this work is to study experimentally and numerically an iced airfoil to assess the aerodynamic penalties associated with presence of ice on the airfoil surface. Three iced shapes were tested on NACA 0012 straight wing at zero and non-zero angles of attack, at Reynolds No. equal to (3.36*105). The 2-D steady state continuity and momentum equations have been solved utilizing finite volume method to analyze the turbulent flow over a clean and iced airfoil. The results show that the ice shapes affected the aerodynamic characteristics due to the change in airfoil shape. The experimental results show that the horn iced airfoil consumes more power than the other shapes of ice, its value was (44.4W). The horn iced shape has the worst effect on the airfoil than the other shapes. The present results are compared with previously reported results; it is found there is a very good agreement between them. A comparison between the experimental and computational results of the presented work were pursing the same behavior.

FastICENet: A real-time and accurate semantic segmentation model for aerial remote sensing river ice image

River ice semantic segmentation is a crucial task, which can provide us with information for river monitoring, disaster forecasting, and transportation management. Previous works mainly focus on higher accuracy acquirement, while efficiency is also important for reality usage. In this paper, a real-time and accurate river ice semantic segmentation network is proposed, named FastICENet. The general architecture consists of two branches, i.e., a shallow high-resolution spatial branch and a deep context semantic branch, which are carefully designed for the scale diversity and irregular shape of river ice in remote sensing images. Then, a novel Downsampling module and a dense connection block based on a lightweight Ghost module are adopted in the context branch to reduce the computation cost. Furthermore, a learnable upsampling strategy DUpsampling is utilized to replace the commonly used bilinear interpolation to improve the segmentation accuracy. We deploy detailed experiments on three publicly available datasets, named NWPU_YRCC_EX, NWPU_YRCC2, and Alberta River Ice Segmentation Dataset. The experimental results demonstrate that our method achieves state-of-the-art performance with competing methods, on the NWPU_YRCC_EX dataset, we can achieve the segmentation speed as 90.84FPS and the segmentation accuracy as 90.770% mIoU, which also illustrates the good leverage between accuracy and speed. Our code is available at https://github.com/nwpulab113/FastICENet

Ice Accretion: Image Post-Processing Measurement Techniques for 2D Ice Shapes

Ice accretion poses substantial safety hazards for the manned and unmanned aviation industries. Its study is essential for icing events risk assessment and for the development of efficient ice protection systems. The existing ice accretion measurement techniques—casting, molding, and laser-scanning—are time-consuming, sometimes cumbersome to use, and highly expensive, while hand tracing is inexpensive, but has lower accuracy and time-consuming post-processing. This work presents two low-cost, fast, and easy-to-use measurement techniques for 2D ice accretion profiles. Both employ algorithms of automatic ice shape detection, one based on unmediated image-processing, another based on the processing of manual ice tracings. The techniques are applied to ice accretion experiments conducted in an icing wind tunnel at low Reynolds numbers, and their results are validated against ice thickness caliper measurements. A comparison of the results shows that both techniques accurately measure the leading-edge ice thickness and the 2D shape of the ice accretion profiles. One technique is faster, with higher measurement accuracy, but produces interrupted-line 2D ice profiles and requires good lighting conditions, while the other generates continuous-line 2D profiles and has no application restriction, but it is slower, with lower accuracy. A discussion is conducted, aiming to help one determine the best applications for each ice accretion measurement technique presented.