Search NASASearch

Engineering topics

David L Rigby

Publications and source records attributed to David L Rigby.

Ice Shape Analysis of an eVTOL Propeller in Forward Flight at the NASA Glenn Icing Research Tunnel

Advanced Air Mobility (AAM) introduces many novel electric vertical takeoff and landing (eVTOL) aircraft configurations for which the effect of icing is not well understood. While icing computational tools have often aided in the design and certification of conventional aircraft, experimental data is needed to support the development and validation of such tools for eVTOL applications. To investigate the icing phenomenon relevant to eVTOL aircraft, NASA Glenn Research Center developed a general-purpose propeller test stand for conducting fundamental icing research on electrically driven propellers in the Icing Research Tunnel (IRT). A 10-day test entry in the IRT was completed in March, 2023, with the goal of generating ice accretions on a non-proprietary propeller geometry under well-characterized conditions. Three carbon fiber propellers of diameters 0.610, 0.711, and 0.914 m (24, 28, and 36 in) were used during the test. Various parametric sweeps of cloud and operating conditions were performed to evaluate the sensitivity of the parameters on the resulting ice accretion characteristics such as ice thickness, mass, location (impingement limits), and type (glaze vs rime). Ice shapes were reviewed for their repeatability, correctness of trends, and sensitivity to parameter variations. The data from the test is being used to support the development of the ice accretion solver, GlennICE.

Icing

Three Dimensional Surface Redefinition Method for Computational Ice Accretion Solvers

Computational tools have been increasing in maturity and are thus commonly used in the engineering design process. Advancements in NASA’s current state of the art computational ice accretion tool are required to tackle the icing challenges of tomorrow. GlennICE, a next generation ice accretion solver, is under development at NASA to tackle these challenges. One of the elements of this migration is transitioning from a quasi-three dimensional strip theory based ice accretion methodology to a fully three dimensional methodology. This requires construction of a method to redefine or extrude a discretized or tessellated surface geometry based on the predicted volumetric ice growth for each tessellated surface triangle. This paper describes the methodology that is employed in the GlennICE software, and assesses the performance of the method in replicating a well defined analytical test case.

Computational

Numerical Investigation of Heat Transfer on the NASA SIDRM Model

NASA has developed a Simulated Inter-Duct Research Model (SIDRM) to study ice crystal icing in aircraft engines. A numerical investigation has been completed using the NASA in-house code, GlennHT. For each flow condition, many different wall thermal boundary conditions were imposed. Six uniform wall temperature cases were generated, with wall temperature relative to the freestream total temperature of 0.85, 0.90, 0.95, 1.05, 1.10, and 1.15. In addition, five cases with uniform heat flux were generated. Heat fluxes were chosen to generate average wall temperatures in the range explored by the uniform wall temperature cases. Included in the uniform heat flux cases is the adiabatic case where heat flux equals zero. The results of these simulations are then used to quantify the degree to which the heat transfer coefficient depends on wall thermal boundary condition. It is seen that the heat transfer coefficient is fairly independent of wall thermal boundary condition, as expected. However, it would appear that, at least in some locations, it may be prudent to attempt to correlate a variation of heat transfer coefficient with wall thermal boundary condition.

Ice Crystal Icing

An Automated Refinement Process for Particle Trajectory Methods in GlennICE

Computational methods for ice accretion can simulate the impact of water drops and ice crystals on an aircraft surface in a Lagrangian reference frame or in the Eulerian reference frame. In the Eulerian reference frame, particles are considered a continuous fluid while in the Lagrangian frame individual particle trajectories are calculated. Methods that use the Eulerian reference frame are typically easier to develop as established modules used for continuum mechanics can be leveraged. The Eulerian systems can also be faster since the user does not have to simulate millions of particles in order to achieve good results. It is imperative therefore that a Lagrangian method optimize the release points of trajectories such that accurate solutions can be obtained while minimizing as much as possible the number of trajectories computed. This paper will present a methodology for this refinement process and demonstrate its effectiveness on sample three dimensional test cases.

William B Wright

GlennICE 2.2 Capabilities and Results

GlennICE (Glenn Icing Computational Environment) is a computational tool designed to calculate ice growth on complex three-dimensional geometries using the input from a user-supplied computational fluid dynamics (CFD) solution for the geometry of interest. The NASA John H. Glenn Research Center at Lewis Field is developing this tool to aid those evaluating, designing and certifying aircraft, engines, and aircraft components for flight in icing conditions. This domestically available software is being developed to enable the introduction of new icing physics into a computational environment in a manner that is open for evaluation and eventual use by industry, academia, and other government organizations. This paper will document the current capabilities for version 2.1 of this software and provide example cases with comparison to available experimental data.

Icing

Utilization of Streamtubes to Analyze the Physical Interaction of a Dispersed Cloud with the CRM65 Hybrid Midspan Model

Recently there have been numerous efforts to identify the relevant icing physics related to the formation of complex three dimensional features, sometimes referred to as ’scallops’, on swept wings. However, much of the physics is still not well understood. This paper computationally investigates the interaction of the icing cloud with the 65 percent Common Research Model (CRM65). Both the interaction with an uniced model, and the interaction with a representative simulated three dimensional ice accretion are analyzed. Preliminary results suggest that the liquid water content increases near the aerodynamic body. For small droplets on the uniced geometry, the particle velocity vector becomes nearly parallel with the aerodynamic model resulting in a low value of collection efficiency. When a three dimensional feature with a length scale much smaller than the leading edge of the airfoil is introduced into the flow, the impingement of small particles can become significantly more perpendicular to the particle velocity vector in the region of these features. Since the liquid water content near the body has increased due to the interaction with the larger scale features of the aerodynamic model, i.e. the leading edge of the swept wing, extremely high collection efficiency is observed. These results suggest that three dimensional features are likely a significant physical driver that should be modeled in some capacity when simulating the impingement of a cloud on an aerodynamic model.

Icing

GlennICE Simulation of 24, 28 and 36 Inch Diameter eVTOL Propellers in Forward Flight

Flow solutions and water collection results are presented for three different propeller sizes. Flow solutions were generated using the NASA FUN3D software, and water collection was simulated using NASA GlennICE software. The geometries considered represent experimental articles that have been tested in the NASA Icing Research Tunnel. Each configuration has four blades. For the experiment, the same spinner centerbody is used for each propeller size. The different sized blades are similar for radii larger than about 40% of total span. Operating conditions were chosen to match operating points where experimental data was collected. Results are presented for both single bin simulations, as well as 7-bin Langmuir D distributions. For 80 microns MVD the Langmuir D distribution showed small differences compared to the single bin. For 15 microns MVD the Langmuir D distribution had a larger impact, especially on the spinner. In the experiment, the centerbody was unchanged as propeller size was changed. In this paper, additional simulations were completed with the 36 inch case perfectly scaled down to 24 and 28 inch. The perfect geometric scaling allows for an assessment of the effect of constant centerbody size on the results. While slight differences in pressure distribution on the propeller are visible between scaled and non-scaled centerbody, the effect seems quite minimal regarding water collection. All simulations were run very cold to ensure rime ice, since handling runback water in the non-inertial frame is beyond the scope of the present work. Ice shapes were generated using a single time step growth. The total mass on the blade and spinner was generally underpredicted for 15 microns MVD, and overpredicted for 80 microns MVD. When comparing ice shapes at operating conditions the numerical results produce very similar ice shapes.

Aircraft Icing