Search NASA⌕ Search

Engineering topics

Sankar, Lakshmi N.

Publications and source records attributed to Sankar, Lakshmi N..

At least 19 records

Application of Extended Messinger Models to Complex Geometries

Since, ice accretion can significantly degrade the performance and the stability of an airborne vehicle, it is imperative to be able to model it accurately. While ice accretion studies have been performed on airplane wings and helicopter blades in abundance, there are few that attempt to model the process on more complex geometries such as fuselages. This paper proposes a methodology that extends an existing in-house Extended Messinger solver to complex geometries by introducing the capability to work with unstructured grids and carry out spatial surface streamwise marching. For the work presented here commercial solvers such as STAR-CCM+ and ANSYS Fluent are used for the flow field and droplet dispersed phase computations. The ice accretion is carried out using an in-house icing solver called GT-ICE. The predictions by GT-ICE are compared to available experimental data, or to predictions by other solvers such as LEWICE and STAR-CCM+. Three different cases with varying levels of complexity are presented. The first case considered is a commercial transport airfoil, followed by a three-dimensional MS(1)-317 swept wing. Finally, ice accretion calculations performed on a Robin fuselage have been discussed. Good agreement with experimental data, where applicable, is observed. Differences between the ice accretion predictions by different solvers have been discussed.

Gupta, Avani↗

Extraction of Dynamic Inflow Models for Coaxial and Tandem Rotors from CFD Simulations

The dynamic inflow coupling with rotor/body dynamics is crucial in the analysis of stability and control law design for helicopters. Over the past several decades, finite-state inflow models for single rotor configurations in hover, forward flight, and maneuver have developed (Ref.1-3). By capturing the interference effects between rotors, the extension of pressure potential finite state inflow model has promising result for coaxial rotor configuration (Ref.4-6). Recently, the focus of the dynamic inflow modeling has shifted to tandem rotor configurations (Ref.7, 8). The development of the dynamic inflow models for tandem rotor configuration still have some limitations due to the lack of knowledge of rotor-to-rotor interference, and rotor-wake interference. Experimental methods, and computational fluid dynamics methods are commonly used to understand the rotor performance and rotor airload variations, and measure or predict inflow velocity distributions at the rotor desk. The inflow distributions are subsequently used to improve the dynamic inflow models. Tandem rotor configurations have been studied experimentally and computationally for several decades (Ref.9-12). Sweet (Ref.10) observed that a tandem rotor with 76-percent-radius overlap required 14% more induced power at hovering condition, relative to an isolated rotor of equivalent disk area. Sweet also found that, above a shaft-to-shaft distance of 1.03 diameter, the performance of the tandem rotor was nearly the same as two isolated rotors. The objective of the present study is to apply computational fluid dynamics simulations of tandem rotors for the extraction of dynamic inflow models. The extended methodology is first validated by comparing the computed induced power against test data. Subsequently inflow distributions and wake structures are analyzed.

Extraction of Dynamic Inflow Models↗

Development and Validation of Physics Based Models for Ice Shedding

Calculations for ice accretion and shedding are presented for a model scale rotor in hover. The aerodynamic characteristics of the rotor are first computed using a combined blade element-momentum theory. The effective angles of attack, and the local flow velocity are used within the NASA Glenn solver LEWICE to estimate the collection efficiency. The computed convection efficiency and the surface pressure distribution from a panel method within LEWICE are used to estimate the ice accretion over the rotor blades over an elapsed time interval. Finally, a force balance approach is used to establish shedding events where the centrifugal force over the ice mass exceeds the adhesive forces at the rotor surface and the cohesive forces between adjacent masses of ice. Preliminary comparisons with test data acquired at the Pennsylvania State Icing Research Tunnel are presented. Sensitivity of the ice shedding events to surface roughness, adhesive strength, cohesive strength, and ambient conditions is discussed.

Gupta, Avani↗

Ice Accretion Modeling using an Eulerian Approach for Droplet Impingement

A three-dimensional Eulerian analysis has been developed for modeling droplet impingement on lifting bodes. The Eulerian model solves the conservation equations of mass and momentum to obtain the droplet flow field properties on the same mesh used in CFD simulations. For complex configurations such as a full rotorcraft, the Eulerian approach is more efficient because the Lagrangian approach would require a significant amount of seeding for accurate estimates of collection efficiency. Simulations are done for various benchmark cases such as NACA0012 airfoil, MS317 airfoil and oscillating SC2110 airfoil to illustrate its use. The present results are compared with results from the Lagrangian approach used in an industry standard analysis called LEWICE.

Kim, Joe Woong↗

Computational Evaluation of the Steady and Pulsed Jet Effects on the Performance of a Circulation Control Wing Section

Circulation Control technology is a very effective way of achieving high lift forces required by aircraft during take-off and landing. This technology can also directly control the flow field over the wing. Compared to a conventional high-lift system, a Circulation Control Wing (CCW) can generate comparable or higher lift forces during take-off/landing with fewer or no moving parts and much less complexity. In this work, an unsteady three-dimensional Navier-Stokes analysis procedure has been developed and applied to Circulation Control Wing configurations. The effects of 2-D steady jets and 2-D pulsed jets on the aerodynamic performance of CCW airfoils have been investigated. It is found that a steady jet can generate very high lift at zero angle of attack without stall, and that a small amount of blowing can eliminate vortex shedding at the trailing edge, a potential noise source. It is also found that a pulsed jet can achieve the same high lift as a steady jet at lower mass flow rates, especially at a high frequency, and that the Strouhal number has a more dominant effect on the pulsed jet performance than just the frequency or the free-stream velocity.

Liu, Yi↗

Numerical Simulations of the Steady and Unsteady Aerodynamic Characteristics of a Circulation Control Wing Airfoil

The aerodynamic characteristics of a Circulation Control Wing (CCW) airfoil have been numerically investigated, and comparisons with experimental data have been made. The configuration chosen was a supercritical airfoil with a 30 degree dual-radius CCW flap. Steady and pulsed jet calculations were performed. It was found that the use of steady jets, even at very small mass flow rates, yielded a lift coefficient that is comparable or superior to conventional high-lift systems. The attached flow over the flap also gave rise to lower drag coefficients, and high L/D ratios. Pulsed jets with a 50% duty cycle were also studied. It was found that they were effective in generating lift at lower reduced mass flow rates compared to a steady jet, provided the pulse frequency was sufficiently high. This benefit was attributable to the fact that the momentum coefficient of the pulsed jet, during the portions of the cycle when the jet was on, was typically twice as much as that of a steady jet.

Liu, Yi↗

Computational Evaluation of the Steady and Pulsed Jet Effects on the Performance of a Circulation Control Wing Section

Circulation Control Wing (CCW) technology is a very effective way of achieving very high lift coefficients needed by aircraft during take-off and landing. This technology can also be used to directly control the flow field over the wing. Compared to a conventional high-lift system, a Circulation Control Wing (CCW) can generate the required values of lift coefficient C(sub L,max) during take-off/landing with fewer or no moving parts and much less complexity. Earlier designs of CCW configurations used airfoils with a large radius rounded trailing edge to maximize the lift benefit. However, these designs also produced very high drag. These high drag levels associated with the blunt, large radius trailing edge can be prohibitive under cruise conditions when Circulation Control is no longer necessary. To overcome this difficulty, an advanced CCW section, i.e., a circulation hinged flap was developed to replace the original rounded trailing edge CC airfoil. This concept developed by Englar is shown. The upper surface of the CCW flap is a large-radius arc surface, but the lower surface of the flap is flat. The flap could be deflected from 0 degrees to 90 degrees. When an aircraft takes-off or lands, the flap is deflected as in a conventional high lift system. Then this large radius on the upper surface produces a large jet turning angle, leading to high lift. When the aircraft is in cruise, the flap is retracted and a conventional sharp trailing edge shape results, greatly reducing the drag. This kind of flap does have some moving elements that increase the weight and complexity over an earlier CCW design. But overall, the hinged flap design still maintains most of the Circulation Control high lift advantages, while greatly reducing the drag in cruising condition associated with the rounded trailing edge CCW design. In the present work, an unsteady three-dimensional Navier-Stokes analysis procedure has been developed and applied to this advanced CCW configuration. The solver can be used in both a 2-D and a 3-D mode, and can thus model airfoils as well as finite wings. The jet slot location, slot height, and the flap angle can all be varied easily and individually in the grid generator and the flow solver. Steady jets, pulsed jets, the leading edge and trailing edge blowing can all be studied with this solver.

Liu, Yi↗

Dynamic Stall Characteristics of Drooped Leading Edge Airfoils

Helicopters in high-speed forward flight usually experience large regions of dynamic stall over the retreating side of the rotor disk. The rapid variations in the lift and pitching moments associated with the stall process can result in vibratory loads, and can cause fatigue and failure of pitch links. In some instances, the large time lag between the aerodynamic forces and the blade motion can trigger stall flutter. A number of techniques for the alleviation of dynamic stall have been proposed and studied by researchers. Passive and active control techniques have both been explored. Passive techniques include the use of high solidity rotors that reduce the lift coefficients of individual blades, leading edge slots and leading edge slats. Active control techniques include steady and unsteady blowing, and dynamically deformable leading edge (DDLE) airfoils. Considerable amount of experimental and numerical data has been collected on the effectiveness of these concepts. One concept that has not received as much attention is the drooped-leading edge airfoil idea. It has been observed in wind tunnel studies and flight tests that drooped leading edge airfoils can have a milder dynamic stall, with a significantly milder load hysteresis. Drooped leading edge airfoils may not, however, be suitable at other conditions, e.g. in hover, or in transonic flow. Work needs to be done on the analysis and design of drooped leading edge airfoils for efficient operation in a variety of flight regimes (hover, dynamic stall, and transonic flow). One concept that is worthy of investigation is the dynamically drooping airfoil, where the leading edge shape is changed roughly once-per-rev to mitigate the dynamic stall.

Sankar, Lakshmi N.↗

Blade vortex interaction noise reduction techniques for a rotorcraft

An active control device for reducing blade-vortex interactions (BVI) noise generated by a rotorcraft, such as a helicopter, comprises a trailing edge flap located near the tip of each of the rotorcraft's rotor blades. The flap may be actuated in any conventional way, and is scheduled to be actuated to a deflected position during rotation of the rotor blade through predetermined regions of the rotor azimuth, and is further scheduled to be actuated to a retracted position through the remaining regions of the rotor azimuth. Through the careful azimuth-dependent deployment and retraction of the flap over the rotor disk, blade tip vortices which are the primary source for BVI noise are (a) made weaker and (b) pushed farther away from the rotor disk (that is, larger blade-vortex separation distances are achieved).

Charles, Bruce D.↗

Blade vortex interaction noise reduction techniques for a rotorcraft

An active control device for reducing blade-vortex interactions (BVI) noise generated by a rotorcraft, such as a helicopter, comprises a trailing edge flap located near the tip of each of the rotorcraft's rotor blades. The flap may be actuated in any conventional way, and is scheduled to be actuated to a deflected position during rotation of the rotor blade through predetermined regions of the rotor azimuth, and is further scheduled to be actuated to a retracted position through the remaining regions of the rotor azimuth. Through the careful azimuth-dependent deployment and retraction of the flap over the rotor disk, blade tip vortices which are the primary source for BVI noise are (a) made weaker and (b) pushed farther away from the rotor disk (that is, larger blade-vortex separation distances are achieved).

Charles, Bruce D.↗

Computational strategies for three-dimensional flow simulations on distributed computer systems

This research effort is directed towards an examination of issues involved in porting large computational fluid dynamics codes in use within the industry to a distributed computing environment. This effort addresses strategies for implementing the distributed computing in a device independent fashion and load balancing. A flow solver called TEAM presently in use at Lockheed Aeronautical Systems Company was acquired to start this effort. The following tasks were completed: (1) The TEAM code was ported to a number of distributed computing platforms including a cluster of HP workstations located in the School of Aerospace Engineering at Georgia Tech; a cluster of DEC Alpha Workstations in the Graphics visualization lab located at Georgia Tech; a cluster of SGI workstations located at NASA Ames Research Center; and an IBM SP-2 system located at NASA ARC. (2) A number of communication strategies were implemented. Specifically, the manager-worker strategy and the worker-worker strategy were tested. (3) A variety of load balancing strategies were investigated. Specifically, the static load balancing, task queue balancing and the Crutchfield algorithm were coded and evaluated. (4) The classical explicit Runge-Kutta scheme in the TEAM solver was replaced with an LU implicit scheme. And (5) the implicit TEAM-PVM solver was extensively validated through studies of unsteady transonic flow over an F-5 wing, undergoing combined bending and torsional motion. These investigations are documented in extensive detail in the dissertation, 'Computational Strategies for Three-Dimensional Flow Simulations on Distributed Computing Systems', enclosed as an appendix.

Sankar, Lakshmi N.↗

Development of iterative techniques for the solution of unsteady compressible viscous flows

A research effort was initiated at Georgia Tech in February 1991 on the development of efficient techniques for the computation of 2-D and 3-D unsteady compressible flow problems. It was found that in 2-D unsteady viscous flow applications, the generalized minimal residual (GMRES) scheme was able to significantly improve the accuracy and stability characteristics of an existing 2-D ADI (Alternating Direction Implicit) time marching scheme. That is, the GMRES/ADI combination allowed 10 to 20 times larger time steps compared to an ADI scheme. Because the GMRES algorithm requires 5 to 10 times the CPU work compared to the ADI scheme, the combined GMRES/ADI scheme yields a net factor of 2 savings in CPU cost. During the past year, we also experimented with GMRES/multigrid/ADI combination. The purpose of this combination was to compute the low frequency components of the change in the flow properties from one time step to the next on a coarse grid. This strategy reduces the memory requirements of the GMRES method roughly by a factor of 4-8 for steady flow problems.

Sankar, Lakshmi N.↗

Development of iterative techniques for the solution of unsteady compressible viscous flows

The development of efficient iterative solution methods for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations is discussed. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. In this work, another approach based on the classical conjugate gradient method, known as the Generalized Minimum Residual (GMRES) algorithm is investigated. The GMRES algorithm has been used in the past by a number of researchers for solving steady viscous and inviscid flow problems. Here, we investigate the suitability of this algorithm for solving the system of non-linear equations that arise in unsteady Navier-Stokes solvers at each time step.

Sankar, Lakshmi N.↗

Numerical investigation of performance degradation of wings and rotors due to icing

The aerodynamic load characteristics and the performance degradation of moderate aspect ratio wings and rotors with simulated glaze leading-edge ice have been studied using a three-dimensional, compressible Navier-Stokes solver. The effect of a splitter plate at the wing root on both clean and iced wing configurations has been studied and the results are compared with the experiment. A significant difference has been observed with and without splitter plates in the magnitude of flow separation and aerodynamic loading at the inboard stations for the iced wing at 8-deg angle of attack. Inviscid calculations were performed and compared with viscous calculations to investigate whether the performance of iced swept wings can be inexpensively predicted using Euler methods. It is shown that inviscid calculations predict higher aerodynamic loading than viscous calculations, and cannot model separation effects. A typical nonlifting helicopter rotor in forward flight condition is also studied, and the penalty due to the leading-edge ice formation on the required torque is numerically demonstrated.

Kwon, Oh J.↗

Development of iterative techniques for the solution of unsteady compressible viscous flows

Efficient iterative solution methods are being developed for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. Thus, the extra work required by iterative schemes can also be designed to perform efficiently on current and future generation scalable, missively parallel machines. An obvious candidate for iteratively solving the system of coupled nonlinear algebraic equations arising in CFD applications is the Newton method. Newton's method was implemented in existing finite difference and finite volume methods. Depending on the complexity of the problem, the number of Newton iterations needed per step to solve the discretized system of equations can, however, vary dramatically from a few to several hundred. Another popular approach based on the classical conjugate gradient method, known as the GMRES (Generalized Minimum Residual) algorithm is investigated. The GMRES algorithm was used in the past by a number of researchers for solving steady viscous and inviscid flow problems with considerable success. Here, the suitability of this algorithm is investigated for solving the system of nonlinear equations that arise in unsteady Navier-Stokes solvers at each time step. Unlike the Newton method which attempts to drive the error in the solution at each and every node down to zero, the GMRES algorithm only seeks to minimize the L2 norm of the error. In the GMRES algorithm the changes in the flow properties from one time step to the next are assumed to be the sum of a set of orthogonal vectors. By choosing the number of vectors to a reasonably small value N (between 5 and 20) the work required for advancing the solution from one time step to the next may be kept to (N+1) times that of a noniterative scheme. Many of the operations required by the GMRES algorithm such as matrix-vector multiplies, matrix additions and subtractions can all be vectorized and parallelized efficiently.

Sankar, Lakshmi N.↗

An efficient hybrid scheme for the analysis of counter rotating propellers

An efficient solution procedure has been developed for analyzing inviscid unsteady flow past counter rotating propellers. This scheme is first order accurate in time and second order in space, and has been extended to fourth order accuracy in the axial direction. The solution procedure has been applied to a 2-bladed SR-7 single rotation propeller and to a GE F7/A7 counter rotation propeller. The pressure coefficients and the global quantities, power and thrust, show good comparison with experimental measurements.

Srivastava, R.↗