Numerical investigation of the effects of icing on fixed and rotary wing aircraft

A 2-D multi-element airfoil code was modified to study the effects of icing on the aerodynamic characteristics of high lift systems. In each zone of the flow field, the solver numerically integrates the 2-D compressible Navier-Stokes equations using a time marching scheme. The surface pressure distribution is generated over a GAW 130 airfoil/flap combination for a flap setting of 25 degrees, and an angle of attack equal to 5 degrees, at a freestream Mach number equal to 0.3. A series of calculations were performed to determine the effects of small scale ice build up on the high lift characteristics of this arifoil/flap combination. The appendix summarizes this progress. Joint studies on correlation of a 3-D iced wing code with experimental data reviewed new measured laser Doppler velocimeter data in the separated region behind the leading edge ice shape. A version of the iced wing analysis using the Roe scheme was developed to evaluate the poor correlation between the computed and measured velocities in the separated region. Work on the extension of the wing-alone analysis to wing body configuration began with modifications to the 3-D iced wing analyses to accept externally generated grids and multi-block grids

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

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

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

During the past two decades, there has been significant progress in the field of numerical simulation of unsteady compressible viscous flows. At present, a variety of solution techniques exist such as the transonic small disturbance analyses (TSD), transonic full potential equation-based methods, unsteady Euler solvers, and unsteady Navier-Stokes solvers. These advances have been made possible by developments in three areas: (1) improved numerical algorithms; (2) automation of body-fitted grid generation schemes; and (3) advanced computer architectures with vector processing and massively parallel processing features. In this work, the GMRES scheme has been considered as a candidate for acceleration of a Newton iteration time marching scheme for unsteady 2-D and 3-D compressible viscous flow calculation; from preliminary calculations, this will provide up to a 65 percent reduction in the computer time requirements over the existing class of explicit and implicit time marching schemes. The proposed method has ben tested on structured grids, but is flexible enough for extension to unstructured grids. The described scheme has been tested only on the current generation of vector processor architecture of the Cray Y/MP class, but should be suitable for adaptation to massively parallel machines

Simulation of unsteady rotational flow over propfan configuration

During the past decade, aircraft engine manufacturers and scientists at NASA have worked on extending the high propulsive efficiency of a classical propeller to higher cruise Mach numbers. The resulting configurations use highly swept twisted and very thin blades to delay the drag divergence Mach number. Unfortunately, these blades are also susceptible to aeroelastic instabilities. This was observed for some advanced propeller configurations in wind tunnel tests at NASA Lewis Research Center, where the blades fluttered at cruise speeds. To address this problem and to understand the flow phenomena and the solid fluid interaction involved, a research effort was initiated at Georgia Institute of Technology in 1986, under the support of the Structural Dynamics Branch of the NASA Lewis Research Center. The objectives of this study are: (1) the development of solution procedures and computer codes capable of predicting the aeroelastic characteristics of modern single and counter-rotation propellers; and (2) the use of these solution procedures to understand physical phenomena such as stall flutter, transonic flutter, and divergence

Studies of unsteady viscous flows using a two-equation model of turbulence

A two equation model of turbulence, based on the turbulent kinetic energy and energy dissipation, suitable for prediction of unsteady viscous flows, was developed. Also, the performance of the two equation model was compared with simpler algebraic models such as the Baldwin-Lomax two layer eddy viscosity model, and a model by Johnson and King which accounts for upstream history of the turbulent kinetic energy. A brief discussion of this study is given

On the Analysis of Trajectories of Gradient Descent in the Optimization of Deep Neural Networks

Theoretical analysis of the error landscape of deep neural networks has garnered significant interest in recent years. In this work, we theoretically study the importance of noise in the trajectories of gradient descent towards optimal solutions in multi-layer neural networks. We show that adding noise (in different ways) to a neural network while training increases the rank of the product of weight matrices of a multi-layer linear neural network. We thus study how adding noise can assist reaching a global optimum when the product matrix is full-rank (under certain conditions). We establish theoretical foundations between the noise induced into the neural network - either to the gradient, to the architecture, or to the input/output to a neural network - and the rank of product of weight matrices. We corroborate our theoretical findings with empirical results.Comment: 4 pages + 1 figure (main, excluding references), 5 pages + 4 figures (appendix

ADINE: An Adaptive Momentum Method for Stochastic Gradient Descent

Two major momentum-based techniques that have achieved tremendous success in optimization are Polyak's heavy ball method and Nesterov's accelerated gradient. A crucial step in all momentum-based methods is the choice of the momentum parameter $m$ which is always suggested to be set to less than $1$. Although the choice of $m < 1$ is justified only under very strong theoretical assumptions, it works well in practice even when the assumptions do not necessarily hold. In this paper, we propose a new momentum based method $\textit{ADINE}$, which relaxes the constraint of $m < 1$ and allows the learning algorithm to use adaptive higher momentum. We motivate our hypothesis on $m$ by experimentally verifying that a higher momentum ($\ge 1$) can help escape saddles much faster. Using this motivation, we propose our method $\textit{ADINE}$ that helps weigh the previous updates more (by setting the momentum parameter $> 1$), evaluate our proposed algorithm on deep neural networks and show that $\textit{ADINE}$ helps the learning algorithm to converge much faster without compromising on the generalization error.Comment: 8 + 1 pages, 12 figures, accepted at CoDS-COMAD 201

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

Application of Navier-Stokes analysis to stall flutter

A solution procedure was developed to investigate the two-dimensional, one- or two-dimensional flutter characteristics of arbitrary airfoils. This procedure requires a simultaneous integration in time of the solid and fluid equations of motion. The fluid equations of motion are the unsteady compressible Navier-Stokes equations, solved in a body-fitted moving coordinate system using an approximate factorization scheme. The solid equations of motion are integrated in time using an Euler implicit scheme. Flutter is said to occur if small disturbances imposed on the airfoil attitude lead to divergent oscillatory motions at subsequent times. The flutter characteristics of airfoils in subsonic speed at high angles of attack and airfoils in high subsonic and transonic speeds at low angles of attack are investigated. The stall flutter characteristics are also predicted using the same procedure