SEARCH · Search NASA
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.
We describe and test the mathematical background for using checksum methods to validate results returned by a numerical subroutine operating in a fault- prone environment that causes unpredictable errors in data.
Optical hardware has been developed to measure the depth of defects in the Space Shuttle Orbiter's windows. In this hardware, a mirror is translated such that its position corresponds to the defect's depth, so the depth measurement problem is transferred to a mirror-position measurement problem. This is preferable because the mirror is internal to the optical system and thus accessible. Based on requirements supplied by the window inspectors, the depth of the defects needs to be measured over a range of 200 microns with a resolution of about 100 nm and an accuracy of about 400 nm. These same requirements then apply to measuring the position of the mirror, and in addition, since this is a scanning system, a response time of about 10 ms is needed. A market search was conducted and no sensor that met these requirements that also fit into the available housing volume (less than one cubic inch) was found, so a novel sensor configuration was constructed to meet the requirements. This new sensor generates a nearly linearly varying magnetic field over a small region of space, which can easily be sampled, resulting in a voltage proportional to position. Experiments were done with a range of inductor values, drive voltages, drive frequencies, and inductor shapes. A rough mathematical model was developed for the device that, in most aspects, describes how it operates and what electrical parameters should be chosen for best performance. The final configuration met all the requirements, yielding a small rugged sensor that was easy to use and had nanometer resolution over more than the 200-micron range required. The inductive position sensor is a compact device (potentially as small as 2 cubic centimeters), which offers nanometer-position resolution over a demonstrated range of nearly 1 mm. One of its advantages is the simplicity of its electrical design. Also, the sensor resolution is nearly uniform across its operational range, which is in contrast to eddy current and capacitive sensors whose sensitivity is dependent upon position.
This paper will summarize the modal analysis and model refinement results for the International Space Station (ISS) Flight-2A mated configuration. To maintain structural integrity of the ISS, structural loading distributions have been rigorously analyzed through numerical simulation and taken into account during the design of the structure and mission operations. The accuracy of the analysis results is directly affected by the precision of mathematical models and estimated input forces. A mathematical model of an ISS configuration is composed of individual component math models. Each component model is required to be correlated with ground test data.
The evaluation of multilayer ceramic capacitors (MLCCs) with base-metal electrodes (BMEs) for potential NASA space project applications requires an in-depth understanding of their reliability. The reliability of an MLCC is defined as the ability of the dielectric material to retain its insulating properties under stated environmental and operational conditions for a specified period of time t. In this presentation, a general mathematic expression of a reliability model for a BME MLCC is developed and discussed. The reliability model consists of three parts: (1) a statistical distribution that describes the individual variation of properties in a test group of samples (Weibull, log normal, normal, etc.), (2) an acceleration function that describes how a capacitors reliability responds to external stresses such as applied voltage and temperature (All units in the test group should follow the same acceleration function if they share the same failure mode, independent of individual units), and (3) the effect and contribution of the structural and constructional characteristics of a multilayer capacitor device, such as the number of dielectric layers N, dielectric thickness d, average grain size r, and capacitor chip size S. In general, a two-parameter Weibull statistical distribution model is used in the description of a BME capacitors reliability as a function of time. The acceleration function that relates a capacitors reliability to external stresses is dependent on the failure mode. Two failure modes have been identified in BME MLCCs: catastrophic and slow degradation. A catastrophic failure is characterized by a time-accelerating increase in leakage current that is mainly due to existing processing defects (voids, cracks, delamination, etc.), or the extrinsic defects. A slow degradation failure is characterized by a near-linear increase in leakage current against the stress time; this is caused by the electromigration of oxygen vacancies (intrinsic defects). The two identified failure modes follow different acceleration functions. Catastrophic failures follow the traditional power-law relationship to the applied voltage. Slow degradation failures fit well to an exponential law relationship to the applied electrical field. Finally, the impact of capacitor structure on the reliability of BME capacitors is discussed with respect to the number of dielectric layers in an MLCC unit, the number of BaTiO3 grains per dielectric layer, and the chip size of the capacitor device.
The Public Entry Risk Assessment (PERA) program addresses risk to the public from shuttle or other spacecraft re-entry trajectories. Managing public risk to acceptable levels is a major component of safe spacecraft operation. PERA is given scenario inputs of vehicle trajectory, probability of failure along that trajectory, the resulting debris characteristics, and field size and distribution, and returns risk metrics that quantify the individual and collective risk posed by that scenario. Due to the large volume of data required to perform such a risk analysis, PERA was designed to streamline the analysis process by using innovative mathematical analysis of the risk assessment equations. Real-time analysis in the event of a shuttle contingency operation, such as damage to the Orbiter, is possible because PERA allows for a change to the probability of failure models, therefore providing a much quicker estimation of public risk. PERA also provides the ability to generate movie files showing how the entry risk changes as the entry develops. PERA was designed to streamline the computation of the enormous amounts of data needed for this type of risk assessment by using an average distribution of debris on the ground, rather than pinpointing the impact point of every piece of debris. This has reduced the amount of computational time significantly without reducing the accuracy of the results. PERA was written in MATLAB; a compiled version can run from a DOS or UNIX prompt.
The force-reflecting teleoperator breadboard is described. It is the first system among available Research and Development systems with the following combined capabilities: (1) The master input device is not a replica of the slave arm. It is a general purpose device which can be applied to the control of different robot arms through proper mathematical transformations. (2) Force reflection generated in the master hand controller is referenced to forces and moments measured by a six DOF force-moment sensor at the base of the robot hand. (3) The system permits a smooth spectrum of operations between full manual, shared manual and automatic, and full automatic (called traded) control. (4) The system can be operated with variable compliance or stiffness in force-reflecting control. Some of the key points of the system are the data handling and computing architecture, the communication method, and the handling of mathematical transformations. The architecture is a fully synchronized pipeline. The communication method achieves optimal use of a parallel communication channel between the local and remote computing nodes. A time delay box is also implemented in this communication channel permitting experiments with up to 8 sec time delay. The mathematical transformations are computed faster than 1 msec so that control at each node can be operated at 1 kHz servo rate without interpolation. This results in an overall force-reflecting loop rate of 200 Hz.
Over the last decades flight vehicles such as aircraft and helicopters entering service and requiring increased operational effectiveness have with few exceptions experienced prolonged flight test development to achieve full certification. In many cases the original requirements had later to be reduced to enable release to service. The impact on the customer, and manufacturer has been considerable leading to increased costs and or reduced operational capabilities. These costly experiences are largely a result of the flight vehicle not behaving as modelled and designed. The evaluation of flight test data can be used as a tool for validating windtunnel results and mathematical models describing the flight dynamical behaviour. In this sense the uncertainty of important aerodynamic stability and control parameters can be reduced and the confidence of aircraft mathematical models improved. An additional important factor comes from the implementation of active control systems offering the promise of significantly increased flight vehicle performance and operational capability. This approach extends the traditional trade-offs between aerodynamics, structures and propulsion systems to include full- time, full-authority fly-by-wire/light systems. It is imperative that the aerodynamic stability and control parameters of such integrated flight and propulsion control systems have to turn out inflight as predicted, since inherent stability margins will be lower and the flight control system must correct these deficiencies to provide flight critical redundancy and safety. With the methodology of system identification from flight tests it is possible to sense the control inputs and the flight vehicle reactions Such as accelerations, rates and attitudes. The mathematical model, e.g. the model structure and parameters, has to be determined from the relationship of the measured control inputs and the system's responses. The aim of this symposium was to review the present state of the art of flight vehicle system and parameter identification techniques, and to provide a critical appraisal of current methods developed and applied to flight test data in a number of NATO nations. Particular emphasis was placed on practical aspects and lessons learned in order to generate information useful to the flight test community in industry and government agencies. The technical papers share invaluable experience and emphasize the advances of flight vehicle system identification over the last years to the point where confidence and robustness level is now reasonably high. The symposium covered overviews of identification methodologies, flight test techniques, recent aircraft and helicopter application programs, and a session of short papers covering up-to-the-minute flight test results. A final discussion included prepared comments from experts and concluded with key issues learned in the application of system identification and future research needs. The essential benefits to NATO nations can be condensed as follows: More accurate mathematical models for high bandwidth flight control systems, Improved assessment and evaluation of flying qualities, High fidelity mathematical models for flight vehicle development and mission training simulators, and generally, Reduced flight test time and costs.
Two forms of high-performance torque control for motor drives have been described in the literature: field orientation control and direct torque control. Field orientation control has been the method of choice for previous NASA electromechanical actuator research efforts with induction motors. Direct torque control has the potential to offer some advantages over field orientation, including ease of implementation and faster response. However, the most common form of direct torque control is not suitable for the highspeed, low-stator-flux linkage induction machines designed for electromechanical actuators with the presently available sample rates of digital control systems (higher sample rates are required). In addition, this form of direct torque control is not suitable for the addition of a high-frequency carrier signal necessary for the "self-sensing" (sensorless) position estimation technique. This technique enables low- and zero-speed position sensorless operation of the machine. Sensorless operation is desirable to reduce the number of necessary feedback signals and transducers, thus improving the reliability and reducing the mass and volume of the system. This research was directed at developing an alternative form of direct torque control known as a "deadbeat," or inverse model, solution. This form uses pulse-width modulation of the voltage applied to the machine, thus reducing the necessary sample and switching frequency for the high-speed NASA motor. In addition, the structure of the deadbeat form allows the addition of the high-frequency carrier signal so that low- and zero-speed sensorless operation is possible. The new deadbeat solution is based on using the stator and rotor flux as state variables. This choice of state variables leads to a simple graphical representation of the solution as the intersection of a constant torque line with a constant stator flux circle. Previous solutions have been expressed only in complex mathematical terms without a method to clearly visualize the solution. The graphical technique allows a more insightful understanding of the operation of the machine under various conditions.
Various Symbolic Manipulation Programs (SMP) were tested to check the functioning of their commands and suitability under various operating systems. Support systems for SMP were found to be relatively better than the one for MACSYMA. The graphics facilities for MACSYMA do not work as expected under the UNIX operating system. Not all commands for MACSYMA function as described in the manuals. Shape representation is a central issue in computer graphics and computer-aided design. Aside from appearance, there are other application dependent, desirable properties like continuity to certain order, symmetry, axis-independence, and variation-diminishing properties. Several shape representations are studied, which include the Osculatory Method, a Piecewise Cubic Polynomial Method using two different slope estimates, Piecewise Cubic Hermite Form, a method by Harry McLaughlin, and a Piecewise Bezier Method. They are applied to collected physical and chemical data. Relative merits and demerits of these methods are examined. Kinematics of a single link, non-dissipative robot arm is studied using MACSYMA. Lagranian is set-up and Lagrange's equations are derived. From there, Hamiltonian equations of motion are obtained. Equations suggest that bifurcation of solutions can occur, depending upon the value of a single parameter. Using the characteristic function W, the Hamilton-Jacobi equation is derived. It is shown that the H-J equation can be solved in closed form. Analytical solutions to the H-J equation are obtained.
Nonlinear observers were originally developed based on the ideas of variable structure control, and for the purpose of detecting disturbances in complex systems. In this anomaly detection application, these observers were designed for estimating the distributed state of fluid flow in a pipe described by a class of advection equations. The observer algorithm uses collected data in a piping system to estimate the distributed system state (pressure and velocity along a pipe containing liquid gas propellant flow) using only boundary measurements. These estimates are then used to further estimate and localize possible anomalies such as leaks or foreign objects, and instrumentation metering problems such as incorrect flow meter orifice plate size. The observer algorithm has the following parts: a mathematical model of the fluid flow, observer control algorithm, and an anomaly identification algorithm. The main functional operation of the algorithm is in creating the sliding mode in the observer system implemented as software. Once the sliding mode starts in the system, the equivalent value of the discontinuous function in sliding mode can be obtained by filtering out the high-frequency chattering component. In control theory, "observers" are dynamic algorithms for the online estimation of the current state of a dynamic system by measurements of an output of the system. Classical linear observers can provide optimal estimates of a system state in case of uncertainty modeled by white noise. For nonlinear cases, the theory of nonlinear observers has been developed and its success is mainly due to the sliding mode approach. Using the mathematical theory of variable structure systems with sliding modes, the observer algorithm is designed in such a way that it steers the output of the model to the output of the system obtained via a variety of sensors, in spite of possible mismatches between the assumed model and actual system. The unique properties of sliding mode control allow not only control of the model internal states to the states of the real-life system, but also identification of the disturbance or anomaly that may occur.
A mixing chamber used in rocket engine testing at the NASA Stennis Space Center is modelled by a system of two nonlinear ordinary differential equations. The mixer is used to condition the thermodynamic properties of cryogenic liquid propellant by controlled injection of the same substance in the gaseous phase. The three inputs of the mixer are the positions of the valves regulating the liquid and gas flows at the inlets, and the position of the exit valve regulating the flow of conditioned propellant. Mixer operation during a test requires the regulation of its internal pressure, exit mass flow, and exit temperature. A mathematical model is developed to facilitate subsequent controller designs. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as uniqueness of the equilibrium point, feedback linearizability and local stability are shown to hold under conditions having direct physical interpretation. The existence of fixed valve positions that attain a desired operating condition is also shown. Validation of the model against real data is likewise provided.
In 2015, a series of systems analysis studies were conducted on John F. Kennedy Airport in New York (NY) in a collaborative effort between NASA and the Port Authority of New York and New Jersey (PANYNJ). This work was performed to build a deeper understanding of NY airspace and operations to determine the improvements possible through operational changes with tools currently available, and where new technology is required for additional improvement. The analysis was conducted using tool-based mathematical analyses, video inspection and evaluation using recorded arrival/departure/surface traffic captured by the Aerobahn tool (used by Kennedy Airport for surface metering), and aural data archives available publically through the web to inform the video segments. A discussion of impacts of trajectory and operational choices on capacity is presented, including runway configuration and usage (parallel, converging, crossing, shared, independent, staggered), arrival and departure route characteristics (fix sharing, merges, splits), and how compression of traffic is staged. The authorization in March of 2015 for New York to use reduced spacing under the Federal Aviation Administration (FAA) Wake Turbulence Recategorization (RECAT) also offers significant capacity benefit for New York airports when fully transitioned to the new spacing requirements, and the impact of these changes for New York is discussed. Arrival and departure capacity results are presented for each of the current day Kennedy Airport configurations. While the tools allow many variations of user-selected conditions, the analysis for these studies used arrival-priority, no-winds, additional safety buffer of 5% to the required minimum spacing, and a mix of traffic typical for Kennedy. Two additional "novel" configurations were evaluated. These configurations are of interest to Port Authority and to their airline customers, and are believed to offer near-term capacity benefit with minimal operational and equipage changes. One of these is the addition of an Optimized Profile Descent (OPD) route to runways 22L and 22R, and the other is the simultaneous use of 4 runways, which is not currently done at Kennedy. The background and configuration for each of these is described, and the capacity results are presented along with a discussion of drawbacks and enablers for each.
The primary experiment on the Geodynamics Experimental Ocean Satellite-3 (GEOS-3) is the radar altimeter. This experiment's major objective is to demonstrate the utility of measuring the geometry of the ocean surface; i.e., the geoid. Results obtained from this experiment so far indicate that the planned objectives of measuring the topography of the ocean surface with an absolute accuracy of + or - 5 m can be met and perhaps exceeded. The GEOS-3 satellite altimeter measurements have an instrument precision in the range of + or - 25 cm to + or - 50 cm when the altimeter is operating in the 'short pulse' mode. After one year's operations of the altimeter, data from over 5000 altimeter passes have been collected. With the mathematical models developed and the altimeter data presently available, mapping of local areas of ocean topography has been realized to the planned accuracy levels and better. This paper presents the basic data processing methods employed and some interesting results achieved with the early data. Plots of mean sea surface heights as inferred by the altimeter measurements are compared with a detailed 1 by 1 deg gravimetric geoid.
The overall objective of this program was the development of control algorithms that allow the concurrent operation of slewing, pointing, vibration, and shape control subsystems. This objective is important for near-term space surveillance missions that require the rapid retargeting and precise pointing of large flexible satellites. The success of these missions requires the design and concurrent operation of the various interacting control subsystems. There were two phases conducted: phase 1 was mathematical model development, and phase 2 was control development. The program is detailed and major conclusions given.
Mathematical models for the photovoltaic conversion of laser power were developed. These models simulate the operation of planar and vertical junction photovoltaic converters and are described in detail.
The spatial resolution that can be attained using scanning linear arrays (consisting of discrete IR solid-state detectors) for image acquisition is considered, and a filtered local projection (FLP) method is described which efficiently combines all available scan information into one rectangular grid without the need for explicit interpolation. Mathematically, the FLP algorithm consists of a localized summation followed by an inverse-filter operation, and it has application to nonlinear restoring techniques. The present method is applied, using a linear array, to simulated data for staggered parallel scans and to multiple scan directions. Noise effects and limitations of the technique are also considered.
A technique is described which can be used for design of feedback controllers for high-performance aircraft operating in flight conditions in which nonlinearities significantly affect performance. Designs are performed on a mathematical model of the longitudinal dynamics of a hypothetical aircraft similar to proposed supermaneuverable flight test vehicles. Nonlinear controller designs are performed using truncated solutions of the Hamilton-Jacobi-Bellman equation. Preliminary results show that the method yields promising results.
A three-dimensional, linearized, Euler analysis is being developed to provide a comprehensive and efficient unsteady aerodynamic analysis for predicting the aeroacoustic and aeroelastic responses of axial-flow turbomachinery blading. The mathematical models needed to describe nonlinear and linearized, inviscid, unsteady flows through a blade row operating within a cylindrical annular duct are presented in this report. A numerical model for linearized inviscid unsteady flows, which couples a near-field, implicit, wave-split, finite volume analysis to far-field eigen analyses, is also described. The linearized aerodynamic and numerical models have been implemented into the three-dimensional unsteady flow code, LINFLUX. This code is applied herein to predict unsteady subsonic flows driven by wake or vortical excitations. The intent is to validate the LINFLUX analysis via numerical results for simple benchmark unsteady flows and to demonstrate this analysis via application to a realistic wake/blade-row interaction. Detailed numerical results for a three-dimensional version of the 10th Standard Cascade and a fan exit guide vane indicate that LINFLUX is becoming a reliable and useful unsteady aerodynamic prediction capability that can be applied, in the future, to assess the three-dimensional flow physics important to blade-row, aeroacoustic and aeroelastic responses.