Engineering PapersSearch

Engineering topics

Vargas, Mario

Publications and source records attributed to Vargas, Mario.

At least 19 records

Fragment Size Distribution for Ice Particle Impacts on a Glass Plate

This work presents the results of an experimental study of ice particle impacts on a flat glass plate. The experiment was conducted at the Ballistics Impact Laboratory of NASA Glenn Research Center. The main objective of the experiment was to gain understanding about the modifications needed to the experimental configuration for a future parametric study at a larger range of values for particle diameters and other parameters. This was achieved by studying the effect of the velocity of an impacting ice particle on the post-impact fragment size and distribution for a reduced range of impacting particle diameters. Pre-impact particle diameter and velocity data were captured with a high-speed side camera. Post-impact fragment data were captured in a single frame with a 29-megapixel camera located above and normal to the target. Repeat runs were conducted for ice particles with diameters ranging from 1.7 to 2.9 millimeters, impacting at velocities between 39 and 98 meters per second. The fragment areas were measured, and the corresponding equivalent diameters and histogram distributions were calculated. Analysis of the data showed that the average equivalent diameter for the fragments in a run was an order of magnitude smaller than the diameter of the impacting ice particle. The histograms for equivalent diameter distribution were nonnormal with long tails, with most of the fragments having equivalent diameters concentrated toward the minimum value of the fragment size that could be resolved. Factors affecting the accuracy of the data during the digital imaging analysis were identified. Needed modifications to the setup to handle small size ice particles and other testing conditions were also identified.

engine icing

Ice Roughness and Thickness Evolution on a Business Jet Airfoil

Experiments were performed in the Icing Research Tunnel at NASA Glenn Research Center (GRC) to investigate the ice roughness and thickness evolution on a 152.4-cm (60-in.) chord business jet airfoil exposed to both Appendix C and Appendix O (SLD - Super-cooled Large Droplet) icing conditions. The resulting measurements demonstrate that the average non-dimensional roughness and the stagnation point thickness scalings are similar to those demonstrated on symmetric wings. However, the surface variations of roughness and thickness exhibit significant differences from those observed on symmetric airfoils. The source of the roughness and thickness differences is the result of surface pressure, velocity and temperature distribution differences from the suction to the pressure sides of the airfoil. LEWICE (LEWis ICE accretion program - software developed at NASA Lewis Research Center - former name of the GRC) simulations are used to further investigate the influences of local collection efficiency and the local freezing fraction on the resulting ice roughness and thickness spatial variations.

icing codes

Ice Accretion Roughness Measurements and Modeling

Roughness on aircraft ice accretions is very important to the overall ice accretion process and to the resulting degradation in aircraft aerodynamic performance. Roughness enhances the local convection leading to more rapid ice accumulation rates, and roughness generates local flow perturbations that lead to higher skin friction. This paper presents 1) a review of the developments in ice shape three-dimensional laser scanning developed at NASA Glenn, 2) a review of the approach of McClain and Kreeger employed to characterize ice roughness evolution on an airfoil surface, and 3) a review of the experimental efforts that have been performed over the last five years to characterize, scale, and model ice roughness evolution physics.

SLD

Ice Roughness and Thickness Evolution on a Swept NACA 0012 Airfoil

Several recent studies have been performed in the Icing Research Tunnel (IRT) at NASA Glenn Research Center focusing on the evolution, spatial variations, and proper scaling of ice roughness on airfoils without sweep exposed to icing conditions employed in classical roughness studies. For this study, experiments were performed in the IRT to investigate the ice roughness and thickness evolution on a 91.44-cm (36-in.) chord NACA 0012 airfoil, swept at 30-deg with 0deg angle of attack, and exposed to both Appendix C and Appendix O (SLD) icing conditions. The ice accretion event times used in the study were less than the time required to form substantially three-dimensional structures, such as scallops, on the airfoil surface. Following each ice accretion event, the iced airfoils were scanned using a ROMER Absolute Arm laser-scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger to determine the spatial roughness variations along the surfaces of the iced airfoils. The resulting measurements demonstrate linearly increasing roughness and thickness parameters with ice accretion time. Further, when compared to dimensionless or scaled results from unswept airfoil investigations, the results of this investigation indicate that the mechanisms for early stage roughness and thickness formation on swept wings are similar to those for unswept wings.

glaze icing

Mechanism of Supercooled Water Droplet Breakup near the Leading Edge of an Airfoil

This work presents the results of an experimental study on supercooled droplet deformation and breakup near the leading edge of an airfoil. The results are compared to prior room temperature droplet deformation results to explore the effects of droplet supercooling. The experiments were conducted in the Adverse Environment Rotor Test Stand (AERTS) at The Pennsylvania State University. An airfoil model placed at the end of the rotor blades mounted onto the hub in the AERTS chamber was moved at speeds ranging between 50 and 80 m/sec. The temperature of the chamber was set at -20˚C. A monotonic droplet generator was used to produce droplets that fell from above, perpendicular to the path of the airfoil. The supercooled state of the droplets was determined by measurement of the temperature of the drops at various locations below the droplet generator exit. A temperature prediction code was also used to estimate the temperature of the droplets based on vertical velocity and the distance traveled by droplets from the droplet generator to the airfoil stagnation line. High speed imaging was employed to observe the interaction between the droplets and the airfoil. The high speed imaging provided droplet deformation information as the droplet approached the airfoil near the stagnation line. A tracking software program was used to measure the horizontal and vertical displacement of the droplet against time. It was demonstrated that to compare the effects of water supercooling on droplet deformation, the ratio of the slip velocity and the initial droplet velocity must be equal. A case with equal slip velocity to initial velocity ratios was selected for room temperature and supercooled droplet conditions. The airfoil velocity was 60 m/s and the slip velocity for both sets of data was 40 m/s. In these cases, the deformation of the weakly supercooled and warm droplets did not present different trends. The similar behavior for both environmental conditions indicates that water supercooling has no effect on particle deformation for the limited range of the weak supercooling of the droplets tested and the selected impact velocity. The assumption of a constant surface tension value was further supported by the equal trend of the Bond number obtained for supercooled and room temperature droplets.

in-flight icing

Computational Modeling Approaches to Multiscale Design of Icephobic Surfaces

To aid in the design of surfaces that prevent icing, a model and computational simulation of impact ice formation at the single droplet scale was implemented. The nucleation of a single supercooled droplet impacting on a substrate, in rime ice conditions, was simulated. Open source computational fluid dynamics (CFD) software was used for the simulation. To aid in the design of surfaces that prevent icing, a model of impact ice formation at the single droplet scale was proposed•No existing model simulates simultaneous impact and freezing of a single super-cooled water droplet•For the 10-week project, a low-fidelity feasibility study was the goal.

adhesion

Characterization of Ice Roughness Variations in Scaled Glaze Icing Conditions

Because of the significant influence of surface tension in governing the stability and breakdown of the liquid film in flooded stagnation regions of airfoils exposed to glaze icing conditions, the Weber number is expected to be a significant parameter governing the formation and evolution of ice roughness. To investigate the influence of the Weber number on roughness formation, 53.3-cm (21-in.) and 182.9-cm (72-in.) NACA 0012 airfoils were exposed to flow conditions with essentially the same Weber number and varying stagnation collection efficiency to illuminate similarities of the ice roughness created on the different airfoils. The airfoils were exposed to icing conditions in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. Following exposure to the icing event, the airfoils were then scanned using a ROMER Absolute Arm scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger (2013) to determine the spatial roughness variations along the surfaces of the iced airfoils. The roughness characteristics on each airfoil were then compared using the relative geometries of the airfoil. The results indicate that features of the ice shape and roughness such as glaze-ice plateau limits and maximum airfoil roughness were captured well by Weber number and collection efficiency scaling of glaze icing conditions. However, secondary ice roughness features relating the instability and waviness of the liquid film on the glaze-ice plateau surface are scaled based on physics that were not captured by the local collection efficiency variations.

Roughness Formation

Convection from Hemispherical and Conical Model Ice Roughness Elements in Stagnation Region Flows

To improve ice accretion prediction codes, more data regarding ice roughness and its effects on convective heat transfer are required. The Vertical Icing Studies Tunnel (VIST) at NASA Glenn Research was used to model realistic ice roughness in the stagnation region of a NACA 0012 airfoil. In the VIST, a test plate representing the leading 2% chord of the airfoil was subjected to flows of 7.62 m/s (25 ft/s), 12.19 m/s (40 ft/s), and 16.76 m/s (55 ft/s). The test plate was fitted with multiple surfaces or sets of roughness panels, each with a different representation of ice roughness. The sets of roughness panels were constructed using two element distribution patterns that were created based on a laser scan of an iced airfoil acquired in the Icing Research Tunnel at NASA Glenn. For both roughness patterns, surfaces were constructed using plastic hemispherical elements, plastic conical elements, and aluminum conical elements. Infrared surface thermometry data from tests run in the VIST were used to calculate area averaged heat transfer coefficient values. The values from the roughness surfaces were compared to the smooth control surface, showing convective enhancement as high as 400% in some cases. The data gathered during this study will ultimately be used to improve the physical modeling in LEWICE or other ice accretion codes and produce predictions of in-flight ice accretion on aircraft surfaces with greater confidence.

Icing Physics

Convective Enhancement of Icing Roughness Elements in Stagnation Region Flows

To improve existing ice accretion simulation codes, more data regarding ice roughness and its effects on convective heat transfer are required. To build on existing research on this topic, this study used the Vertical Icing Studies Tunnel (VIST) at NASA Glenn Research to model realistic ice roughness in the stagnation region of a NACA 0012 airfoil. Using the VIST, a test plate representing the leading 2% chord of the airfoil was subjected to flows of 7.62 m/s (25 ft/s), 12.19 m/s (40 ft/s), and 16.76 m/s (55 ft/s). The test plate was fitted with 3 surfaces, each with a different representation of ice roughness: 1) a control surface with no ice roughness, 2) a surface with ice roughness with element height scaled by 10x and streamwise rough zone width from the stagnation point scaled by 10x, and 3) a surface with ice roughness with element height scaled by 10x and streamwise rough zone width from the stagnation point scaled by 25x. Temperature data from the tests were recorded using an infrared camera and thermocouples imbedded in the test plate. From the temperature data, a convective heat transfer coefficient map was created for each case. Additional testing was also performed to validate the VIST's flow quality. These tests included five-hole probe and hot-wire probe velocity traces to provide flow visualization and to study boundary layer formation on the various test surfaces. The knowledge gained during the experiments will help improve ice accretion codes by providing heat transfer coefficient validation data and by providing flow visualization data helping understand current and future experiments performed in the VIST.

Icing

Ice Particle Impacts on a Flat Plate

An experimental study was conducted at the Ballistic Laboratory of NASA Glenn Research Center to study the impact of ice particles on a stationary flat surface target set at 45 degrees with respect to the direction of motion of the impinging particle (Figure 1). The experiment is part of NASA efforts to study the physics involved in engine power-loss events due to ice-crystal ingestion and ice accretion formation inside engines. These events can occur when aircraft encounter high-altitude convective weather.

ice particle

Ice Particle Impacts on a Moving Wedge

This work presents the results of an experimental study of ice particle impacts on a moving wedge. The experiment was conducted in the Adverse Environment Rotor Test Stand (AERTS) facility located at Penn State University. The wedge was placed at the tip of a rotating blade. Ice particles shot from a pressure gun intercepted the moving wedge and impacted it at a location along its circular path. The upward velocity of the ice particles varied from 7 to 12 meters per second. Wedge velocities were varied from 0 to 120 meters per second. Wedge angles tested were 0, 30, 45, and 60. High speed imaging combined with backlighting captured the impact allowing observation of the effect of velocity and wedge angle on the impact and the post-impact fragment behavior. It was found that the pressure gun and the rotating wedge could be synchronized to consistently obtain ice particle impacts on the target wedge. It was observed that the number of fragments increase with the normal component of the impact velocity. Particle fragments ejected immediately after impact showed velocities higher than the impact velocity. The results followed the major qualitative features observed by other researchers for hailstone impacts, even though the reduced scale size of the particles used in the present experiment as compared to hailstones was 4:1.

Impacts

Ice Particle Impacts on a Moving Wedge

This work presents the results of an experimental study of ice particle impacts on a moving wedge. The experiment was conducted in the Adverse Environment Rotor Test Stand (AERTS) facility located at Penn State University. The wedge was placed at the tip of a rotating blade. Ice particles shot from a pressure gun intercepted the moving wedge and impacted it at a location along its circular path. The upward velocity of the ice particles varied from 7 to 12 meters per second. Wedge velocities were varied from 0 to 120 meters per second. Wedge angles tested were 0 deg, 30 deg, 45 deg, and 60 deg. High speed imaging combined with backlighting captured the impact allowing observation of the effect of velocity and wedge angle on the impact and the post-impact fragment behavior. It was found that the pressure gun and the rotating wedge could be synchronized to consistently obtain ice particle impacts on the target wedge. It was observed that the number of fragments increase with the normal component of the impact velocity. Particle fragments ejected immediately after impact showed velocities higher than the impact velocity. The results followed the major qualitative features observed by other researchers for hailstone impacts, even though the reduced scale size of the particles used in the present experiment as compared to hailstones was 4:1.

Roughness

Ice Roughness in Short Duration SLD Icing Events

Ice accretion codes depend on models of roughness parameters to account for the enhanced heat transfer during the ice accretion process. While mitigating supercooled large droplet (SLD or Appendix O) icing is a significant concern for manufacturers seeking future vehicle certification due to the pending regulation, historical ice roughness studies have been performed using Appendix C icing clouds which exhibit mean volumetric diameters (MVD) much smaller than SLD clouds. Further, the historical studies of roughness focused on extracting parametric representations of ice roughness using multiple images of roughness elements. In this study, the ice roughness developed on a 21-in. NACA 0012 at 0deg angle of attack exposed to short duration SLD icing events was measured in the Icing Research Tunnel at the NASA Glenn Research Center. The MVD's used in the study ranged from 100 micrometer to 200 micrometers, in a 67 m/s flow, with liquid water contents of either 0.6 gm/cubic meters or 0.75 gm/cubic meters. The ice surfaces were measured using a Romer Absolute Arm laser scanning system. The roughness associated with each surface point cloud was measured using the two-dimensional self-organizing map approach developed by McClain and Kreeger (2013) resulting in statistical descriptions of the ice roughness.

Icing

Drag Coefficient of Water Droplets Approaching the Leading Edge of an Airfoil

This work presents results of an experimental study on droplet deformation and breakup near the leading edge of an airfoil. The experiment was conducted in the rotating rig test cell at the Instituto Nacional de Tecnica Aeroespacial (INTA) in Madrid, Spain. An airfoil model was placed at the end of the rotating arm and a monosize droplet generator produced droplets that fell from above, perpendicular to the path of the airfoil. The interaction between the droplets and the airfoil was captured with high speed imaging and allowed observation of droplet deformation and breakup as the droplet approached the airfoil near the stagnation line. Image processing software was used to measure the position of the droplet centroid, equivalent diameter, perimeter, area, and the major and minor axes of an ellipse superimposed over the deforming droplet. The horizontal and vertical displacement of each droplet against time was also measured, and the velocity, acceleration, Weber number, Bond number, Reynolds number, and the drag coefficients were calculated along the path of the droplet to the beginning of breakup. Results are presented and discussed for drag coefficients of droplets with diameters in the range of 300 to 1800 micrometers, and airfoil velocities of 50, 70 and 90 meters/second. The effect of droplet oscillation on the drag coefficient is discussed.

drag

Mechanism of Water Droplet Breakup Near the Leading Edge of an Airfoil

This work presents results of an experimental study on droplet deformation and breakup near the leading edge of an airfoil. The experiment was conducted in the rotating rig test cell at the Instituto Nacional de Tecnica Aeroespacial (INTA) in Madrid, Spain. The airfoil model was placed at the end of the rotating arm and a monosize droplet generator produced droplets that fell from above, perpendicular to the path of the airfoil. The interaction between the droplets and the airfoil was captured with high speed imaging and allowed observation of droplet deformation and breakup as the droplet approached the airfoil near the stagnation line. Image processing software was used to measure the position of the droplet centroid, equivalent diameter, perimeter, area, and the major and minor axes of an ellipse superimposed over the deforming droplet. The horizontal and vertical displacement of each droplet against time was also measured, and the velocity, acceleration, Weber number, Bond number, Reynolds number, and the drag coefficients were calculated along the path of the droplet to the beginning of breakup. Droplet deformation is defined and studied against main parameters. The high speed imaging allowed observation of the actual mechanism of breakup and identification of the sequence of configurations from the initiation of the breakup to the disintegration of the droplet. Results and comparisons are presented for droplets of diameters in the range of 500 to 1800 microns, and airfoil velocities of 70 and 90 m/sec.

Vargas, Mario

Droplet Deformation Prediction With the Droplet Deformation and Breakup Model (DDB)

The Droplet Deformation and Breakup Model was used to predict deformation of droplets approaching the leading edge stagnation line of an airfoil. The quasi-steady model was solved for each position along the droplet path. A program was developed to solve the non-linear, second order, ordinary differential equation that governs the model. A fourth order Runge-Kutta method was used to solve the equation. Experimental slip velocities from droplet breakup studies were used as input to the model which required slip velocity along the particle path. The center of mass displacement predictions were compared to the experimental measurements from the droplet breakup studies for droplets with radii in the range of 200 to 700 mm approaching the airfoil at 50 and 90 m/sec. The model predictions were good for the displacement of the center of mass for small and medium sized droplets. For larger droplets the model predictions did not agree with the experimental results.

Vargas, Mario