Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

Airfoil Design and Optimization by the One-Shot Method

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Lagrange multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Kuruvila, G.↗

Integration of an Autopilot for a Micro Air Vehicle

Two autopilots providing autonomous flight capabilities are presented herein. The first is the Pico-Pilot, demonstrated for the 12-inch size class of micro air vehicles. The second is the MicroPilot MP2028(sup g), where its integration into a 36-inch Zagi airframe (tailless, elevons only configuration) is investigated and is the main focus of the report. Analytical methods, which include the use of the Advanced Aircraft Analysis software from DARCorp, were used to determine the stability and control derivatives, which were then validated through wind tunnel experiments. From the aerodynamic data, the linear, perturbed equations of motion from steady-state flight conditions may be cast in terms of these derivatives. Using these linear equations, transfer functions for the control and navigation systems were developed and feedback control laws based on Proportional, Integral, and Derivative (PID) control design were developed to control the aircraft. The PID gains may then be programmed into the autopilot software and uploaded to the microprocessor of the autopilot. The Pico-Pilot system was flight tested and shown to be successful in navigating a 12-inch MAV through a course defined by a number of waypoints with a high degree of accuracy, and in 20 mph winds. The system, though, showed problems with control authority in the roll and pitch motion of the aircraft: causing oscillations in these directions, but the aircraft maintained its heading while following the prescribed course. Flight tests were performed in remote control mode to evaluate handling, adjust trim, and test data logging for the Zagi with integrated MP2028(sup g). Ground testing was performed to test GPS acquisition, data logging, and control response in autonomous mode. Technical difficulties and integration limitations with the autopilot prevented fully autonomous flight from taking place, but the integration methodologies developed for this autopilot are, in general, applicable for unmanned air vehicles within the 36-inch size class or larger that use a PID control based autopilot.

Platanitis, George↗

Quantification of the Land Surface and Brown Ocean Influence on Tropical Cyclone Intensification Over Land

An investigation of TC Kelvin in February 2018 over northeast Australia was conducted to understand the mechanisms of the Brown Ocean effect (BOE), and to develop a comprehensive analysis framework for land falling tropical cyclones (TCs) in the process. NASA’s Land Information System (LIS) coupled to the NASA Unified WRF (NU-WRF) system was employed as the numerical model framework for 12 land/soil moisture perturbation experiments. Impacts of soil moisture and surface enthalpy flux conditions on TC Kelvin were investigated by closely evaluating simulated track and intensity, midlevel atmospheric thermodynamic properties, vertical wind shear, total precipitable water (TPW), and surface moisture flux. The results suggest that there were recognized differentiations among the sensitivity simulations as a result of land surface (e.g., soil moisture and texture) conditions. However, the intensification of TC Kelvin over land was more strongly related to atmospheric moisture advection and the diurnal cycle of solar radiation (i.e., radiative cooling) than to overall soil moisture conditions or surface fluxes. The analysis framework employed here for TC Kelvin can serve as a foundation to specifically quantify the factors governing the BOE. It also demonstrates that the BOE is not a binary influence (i.e. all or nothing), but instead operates in a continuum from largely to minimally influential such that it could be utilized to help improve prediction of inland effects for all land falling TCs.

Jinwoong Yoo↗

Stochastic stability and instability of model ecosystems

In this work, we initiate a stability study of multispecies communities in stochastic environment by using Ito's differential equations as community models. By applying the direct method of Liapunov, we obtain sufficient conditions for stability and instability in the mean of the equilibrium populations. The conditions are expressed in terms of the dominant diagonal property of community matrices, which is a suitable mechanism for resolving the central problem of 'complexity vs stability' in model ecosystems. As a by-product of this analysis we exhibit important structural properties of the stochastic density-dependent models, and establish tolerance of community stability to a broad class of nonlinear time-varying perturbations.

Ladde, G. S.↗

Design sensitivity analysis with Applicon IFAD using the adjoint variable method

A numerical method is presented to implement structural design sensitivity analysis using the versatility and convenience of existing finite element structural analysis program and the theoretical foundation in structural design sensitivity analysis. Conventional design variables, such as thickness and cross-sectional areas, are considered. Structural performance functionals considered include compliance, displacement, and stress. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. That is, design sensitivity analysis software does not have to be imbedded in an existing finite element code. The finite element structural analysis program used in the implementation presented is IFAD. Feasibility of the method is shown through analysis of several problems, including built-up structures. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of a finite difference perturbation.

Frederick, Marjorie C.↗

Implementation of design sensitivity analysis with existing finite element codes

A numerical method is presented to implement structural design sensitivity analysis theory, using the versatility and convenience of existing finite element structural analysis programs. Design variables such as thickness and cross-sectional areas of components of individual members and built-up structures are considered. Structural performance functionals considered include displacement and stress. The method is also applicable for eigenvalue problem design sensitivity analysis. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. Thus, design sensitivity analysis software does not have to be embedded in an existing finite element code. Feasability of the method is shown through analysis of several problems, including a built-up structure. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of finite difference perturbations.

Choi, Kyung K.↗

A simplified model of the modified Hadley circulation

In this paper, a simplified model of modified Hadley circulation that is amenable to analysis is developed to find the low-latitude upper-tropospheric circulation which would exist in the absence of transient eddies. The model is developed as a perturbation of the ideal Hadley circulation solution due to zonally assymmetric thermal forcing and is solved semianalytically. The modified Hadley circulation, found using forcing derived from observed outgoing long-wave radiation, exhibits some similarities to observed features of the tropical and subtropical time-mean upper-tropospheric flow, particularly with regard to the velocity potential.

Schneider, Edwin K.↗

Singular perturbations and time scales (SPaTS) in discrete control systems-An overview

Recent developments in the theory of singular perturbations and time scales (SPaTS) in discrete control systems are reviewed. Sources of discrete models and the effect of the discretizing interval on the model are examined. The analysis of two-time scale systems is presented to bring out typical characteristic features of SPaTS. In the control of the two-time scale systems, the important issue of multirate sampling is addressed.

Naidu, D. S.↗

Propagation of large amplitude Alfven waves in the solar wind neutral sheet

Analysis of solar wind fluctuation data show that the correlation between velocity and magnetic field fluctuations decreases when going farther away from the Sun. This decorrelation can be attributed either to the time evolution of the fluctuations, carried away by the solar wind, or to the interaction between the solar wind neutral sheet and Alfven waves. To check this second hypothesis we have numerically studied the propagation of Alfven waves in the solar wind neutral sheet. The initial conditions have been set up in order to guarantee B(exp 2) = const, so that the following numerical evolution is only due to the inhomogeneity in the background magnetic field. The analysis of the results shows that compressive structures are formed, mainly in the neutral sheet where they have been identified as pressure balanced structures, i.e., tangential discontinuities. Fast perturbations, which are also produced, have a tendency to leave the simulation domain, propagating also perpendicularly to the mean magnetic field. For this reason the level of fast perturbations is always smaller with respect to the previously cited plasma balanced structures, which are slow mode perturbations. A comparison between the numerical results and some particular observational issues is also presented.

Malara, F.↗

Physics of Boundaries and their Interactions in Space Plasmas

This report describes the work done by SciberNet, Inc. during the month of August. We have resolved the issues associated with the implementation of the dipole field in our large scale hybrid simulations of the magnetopause. We have setup several runs and will spend the next several months analyzing the data. The results will be presented at the Fall AGU. We are also continuing our analysis of the 3-D simulations of thin current sheets at the magnetopause, paying special attention to the conditions under which Kelvin-Helmholtz would lead to sizable perturbations of the magnetopause. In a related study, we are in the process of developing a new kinetic linear code that would for the first time enable us to examine the linear properties of the Kelvin-Helmholtz instability in the fully kinetic regime. Finally, we are continuing our code development to include inflow-outflow boundary conditions in our 2-D and 3-D hybrid codes. We are also comparing the different methods of code parallelization in order to extend the limits of our calculations.

Omidi, Nojan↗

Development of Advanced Methods of Structural and Trajectory Analysis for Transport Aircraft

This paper develops a near-optimal guidance law for generating minimum fuel, time, or cost fixed-range trajectories for supersonic transport aircraft. The approach uses a choice of new state variables along with singular perturbation techniques to time-scale decouple the dynamic equations into multiple equations of single order (second order for the fast dynamics). Application of the maximum principle to each of the decoupled equations, as opposed to application to the original coupled equations, avoids the two point boundary value problem and transforms the problem from one of a functional optimization to one of multiple function optimizations. It is shown that such an approach produces well known aircraft performance results such as minimizing the Brequet factor for minimum fuel consumption and the energy climb path. Furthermore, the new state variables produce a consistent calculation of flight path angle along the trajectory, eliminating one of the deficiencies in the traditional energy state approximation. In addition, jumps in the energy climb path are smoothed out by integration of the original dynamic equations at constant load factor. Numerical results performed for a supersonic transport design show that a pushover dive followed by a pullout at nominal load factors are sufficient maneuvers to smooth the jump.

Ardema, Mark D.↗

Analysis of Stellar Occultation and Imaging Data for the Rings of Saturn and Uranus

The principal goals of this research were as follows: 1. To carry out photometric modelling of Saturn's B Ring; 2. To re-examine dynamical perturbations of the Uranian c ring, 'While delivering relevant occultation data to the PDS Rings Node; and 3. Prepare data from the Saturn ring-plane crossing observations for delivery to the Rings Node.

Nicholson, Philip D.↗

Solar Dynamo Driven by Periodic Flow Oscillation

We have proposed that the periodicity of the solar magnetic cycle is determined by wave mean flow interactions analogous to those driving the Quasi Biennial Oscillation in the Earth's atmosphere. Upward propagating gravity waves would produce oscillating flows near the top of the radiation zone that in turn would drive a kinematic dynamo to generate the 22-year solar magnetic cycle. The dynamo we propose is built on a given time independent magnetic field B, which allows us to estimate the time dependent, oscillating components of the magnetic field, (Delta)B. The toroidal magnetic field (Delta)B(sub phi) is directly driven by zonal flow and is relatively large in the source region, (Delta)(sub phi)/B(sub Theta) much greater than 1. Consistent with observations, this field peaks at low latitudes and has opposite polarities in both hemispheres. The oscillating poloidal magnetic field component, (Delta)B(sub Theta), is driven by the meridional circulation, which is difficult to assess without a numerical model that properly accounts for the solar atmosphere dynamics. Scale-analysis suggests that (Delta)B(sub Theta) is small compared to B(sub Theta) in the dynamo region. Relative to B(sub Theta), however, the oscillating magnetic field perturbations are expected to be transported more rapidly upwards in the convection zone to the solar surface. As a result, (Delta)B(sub Theta) (and (Delta)B(sub phi)) should grow relative to B(sub Theta), so that the magnetic fields reverse at the surface as observed. Since the meridional and zonai flow oscillations are out of phase, the poloidal magnetic field peaks during times when the toroidal field reverses direction, which is observed. With the proposed wave driven flow oscillation, the magnitude of the oscillating poloidal magnetic field increases with the mean rotation rate of the fluid. This is consistent with the Bode-Blackett empirical scaling law, which reveals that in massive astrophysical bodies the magnetic moment tends to increase with the angular momentum of the fluid.

Mayr, Hans G.↗

Gravity Estimation from a Simulated GRACE Mission: Short vs. Long Arcs

We present simulations of gravity estimation from a GRACE-like satellite mission: low-low intersatellite tracking with a precision of order 1 micron s(exp -1) yielding gravity fields of degree and order 120. We employ a unique parameterization of the intersatellite baseline vector which allows the gravity estimation to be performed (relatively) independently of the GPS (global positioning system) tracking data once sufficiently accurate orbits are obtained. This considerably simplifies data processing during the gravity estimation. During that process only certain components of the baseline parameterization need be adjusted; other components are uncorrelated with gravity and may be adopted unchanged from the initial GPS orbits. The technique is also amenable to very short arcs of data. We present comparisons of gravity estimation from 30 days of observations with arcs of length 15 minutes vs. arcs of one day. Our 'truth' field is the EGM96 (Earth Gravitational Model) model; our prior field is a degree-70 clone of EGM96, perturbed from it by amounts comparable to the standard errors of EGM96 (and identically zero for degrees 71-120). For a high inclination orbit, the short-arc analysis recovers low order gravity coefficients remarkably well, although higher order terms, especially sectorial terms, are understandably less accurate. The simulations suggest that either long or short arcs of GRACE data are likely to improve parts of the geopotential spectrum by several orders of magnitude. This is especially so for low order coefficients, which are markedly improved for all degrees through 120.

Rowlands, David D.↗

Anomalistic Disturbance Torques during the Entry Phase of the Mars Exploration Rover Missions: A Telemetry and Mars-Surface Investigation

Shortly after landing on Mars, post-flight analysis of the "Spirit" entry data suggested that the vehicle experienced large, anomalistic oscillations in angle-of-attack starting at about M=6. Similar analysis for "Opportunity " found even larger oscillations starting immediately after maximum dynamic pressure at M=14. Where angles-of-attack of 1-2 degrees were expected from maximum dynamic pressure to drogue deployment, the reconstructions suggested 4 to 9 degrees. The next Mars lander, 2007 Phoenix project, was concerned enough to recommend further exploration of the anomalies. Detailed analysis of "Opportunity" data found significant anomalies in the hypersonic aerodynamic torques. The analysis showed that these torques were essentially fixed in the spinning vehicle. Nearly a year after landing, the "Oportunity" rover took pictures of its aeroshell on the surface, which showed that portions of the aeroshell thermal blanket assembly still remained. This blanket assembly was supposed to burn off very early in the entry. An analysis of the aeroshell photographs led to an estimate of the aerodynamic torques that the remnants could have produced. A comparison of two estimates of the aerodynamic torque perturbations (one extracted from telemetry data and the other from Mars surface photographs) showed exceptional agreement. Trajectory simulations using a simple data derived torque perturbation model provided rigid body motions similar to that observed during the "Opportunity" entry. Therefore, the case of the anomalistic attitude behavior for the "Opportunity" EDL is now considered closed and a suggestion is put forth that a similar event occurred for the "Spirit" entry as well.

Tolson, Robert H.↗

Rotor-Fuselage Interaction: Analysis and Validation with Experiment

The problem of rotor-fuselage aerodynamic interaction has to be considered in industry applications from various aspects. First, in order to increase helicopter speed and reduce operational costs, rotorcraft tend to be more and more compact, with a main rotor closer to the fuselage surface. This creates significant perturbations both on the main rotor and on the fuselage, including steady and unsteady effects due to blade and wake passage and perturbed inflow at the rotor disk. Furthermore,the main rotor wake affects the tail boom, empennage and anti-torque system. This has important consequences for helicopter control and vibrations at low speeds and also on tail rotor acoustics (main rotor wake-tail rotor interactions). This report describes the US Army-France MOD cooperative work on this problem from both the theoretical and experimental aspects. Using experimental 3D velocity field and fuselage surface pressure measurements, three codes that model the interactions of a helicopter rotor with a fuselage are compared. These comparisons demonstrate some of the strengths and weaknesses of current models for the combined rotor-fuselage analysis.

Berry, John D.↗

An Aircraft Separation Algorithm with Feedback and Perturbation

A separation algorithm is a set of rules that tell aircraft how to maneuver in order to maintain a minimum distance between them. This paper investigates demonstrating that separation algorithms satisfy the FAA requirement for the occurrence of incidents by means of simulation. Any demonstration that a separation algorithm, or any other aspect of flight, satisfies the FAA requirement is a challenge because of the stringent nature of the requirement and the complexity of airspace operations. The paper begins with a probability and statistical analysis of both the FAA requirement and demonstrating meeting it by a Monte Carlo approach. It considers the geometry of maintaining separation when one plane must change its flight path. It then develops a simple feedback control law that guides the planes on their paths. The presence of feedback control permits the introduction of perturbations, and the stochastic nature of the chosen perturbation is examined. The simulation program is described. This paper is an early effort in the realistic demonstration of a stringent requirement. Much remains to be done.

White, Allan L.↗

The symmetric turbulent plane wake downstream of a sharp trailing edge

The analysis and modeling of the symmetric turbulent plane wake downstream of a sharp trailing edge is addressed. A compact description of the flow near the trailing edge is formulated using the results of a previous asymptotic analysis. The new description retains the two-layered structure identified in the previous work, and it clarifies the principal dynamics of the flow in the near-wake outer layer, away from the wake centerline. For zero-pressure-gradient flow, the near-wake outer layer is shown to be represented to leading order by the similarity solution that governs the outer region of the surface boundary layer. The leading perturbation in the outer layer due to the developing near-wake inner-layer flow is identified, and this is shown to be asymptotically smaller than undetermined higher-order terms associated with the surface boundary-layer flow. Results of the new near-wake analysis are used to formulate an algebraic eddy viscosity model for wake flow predictions at arbitrary distances from the trailing edge. The model is used in a numerical solution of the boundary layer equations, and computed velocity and Reynolds stress profiles are shown to compare well with experimental data.

Bogucz, E. A.↗