Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reduced order functions”

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 73 records · Page 4

Application of the method of integral relations to boundary layer flows over blunt bodies

Calculations of boundary layer flows past blunt bodies at angles of incidence are presented. Using the method of integral relations together with the method of lines, the full three-dimensional boundary layer equations are reduced to a system of first-order ordinary differential equations. The streamwise shear stress function and the cross-flow velocity component are represented as suitable functions of the streamwise velocity component. The role of the zone of dependence is automatically satisfied by the choice of differencing in the method of lines. Solutions correct to the second order are obtained in the positive shear region for flow over an ellipsoid at 30-deg incidence. The results are compared with corresponding finite difference solutions.

Modarress, D.↗

Linear stability analysis of the Vlasov-Poisson equations in high density plasmas in the presence of crossed fields and density gradients

The equations for the single-particle orbits in a nonneutral high density plasma in the presence of inhomogeneous crossed fields are obtained. Using these orbits, the linearized Vlasov equation is solved as an expansion in the orbital radii in the presence of inhomogeneities and density gradients. A model distribution function is introduced whose cold-fluid limit is exactly the same as that used in many previous studies of the cold-fluid equations. This model function is used to reduce the linearized Vlasov-Poisson equations to a second-order ordinary differential equation for the linearized electrostatic potential whose eigenvalue is the perturbation frequency.

Kaup, D. J.↗

Shape and Orientation of Mercury from Radar Ranging Data

If Mercury's radius is expanded in Legendre functions to the second degree and order, the systematic error in radar ranging can be reduced significantly. We interpret the expansion coefficients in terms of a best-fit ellipsoid displaced with respect to the center of mass in the equatorial plane. The ellipsoid's principal axes are rotated in the equatorial plane such that the long axis is aligned with cartographic longitude 15.3 deg +/- 2.9 deg (west). The pole location is consistent with the IAU pole, normal to Mercury's orbital plane. There is a significant equatorial ellipticity (a - b)/a = (540 +/- 54) x 10(exp -6). The center of figure is offset from the center of mass (C.F.-C.M.) by 640 +/- 78 m in the equatorial plane in the direction of cartographic longitude 319.5 deg +/- 6.9 deg. The magnitude of the equatorial center of figure offset implies an excess crustal thickness of 12 km or less, comparable to the Moon's excess. By comparing the equatorial ellipticity with the Mariner 10 gravity coefficient C(sub 22) and assuming Airy isostatic compensation, we conclude that Mercury's crustal thickness is in the range from 100 to 300 km.

Anderson, J. D.↗

Mercury's Global Topography from Radar Ranging Data

When Mercury's radius is expanded in Legendre functions to the second degree and order, the systematic error in radar ranging data is reduced substantially. Previously, data spanning an observing interval from 1966 to 1990 were used to infer an equatorial ellipticity (a - b)/a = (540 +/- 54) X 10(exp -6) and a center-of-figure minus center-of-mass offset of (640 +/- 78) m. The magnitude of this equatorial center of figure offset implies an excess crustal thickness of 12 km or less, comparable to the Moon's excess. By comparing the equatorial ellipticity with the Mariner 10 gravity field, and assuming Airy isostatic compensation, bounds on crustal thickness can be derived. Mercury's crustal thickness is in the range from 100 to 300 km. The Mercury radar ranging observing interval has been extended from 1966 to the present. In addition, improvements in data reduction techniques have resulted in a set of Mercury ranging data less affected by systematic error, in particular the biases introduced by local topographic variations. We use this new set of reduced ranging data to improve Mercury's global topography and center-of-figure minus center-of-mass offset. New results on crustal thickness are derived, and prospects for further improvement with Mercury Orbiter data are discussed.

Anderson, J. D.↗

The General Equations of Equilibrium of Rotationally Symmetric Membranes and some Static Solutions for Uniform Centrifugal Loading

The general equations of equilibrium of a rotationally symmetrical body undergoing large arbitrary deflections but small strains are derived by means of the principle of virtual work. All deflections are referred to a reference surface which, in general, may be taken to coincide with any intermediate state of the membrane. Static solutions of the general equations for steadily rotating membranes are considered for the cases in which the final or the initial configuration of the membrane is specified. For problems in which the final configuration is specified, conditions are given under which such a configuration is obtainable from a stress-free initial state. Such conditions being satisfied, it is shown that the stresses and radial displacements are determinate while axial displacements may be obtained by quadrature. For problems in which the initial configuration is given, the determination of the stresses and deflections is reduced to the solution of a second-order nonlinear differential equation for a stress function. An approximate method of solution based on boundary-layer considerations is presented for the case of a spinning sphere of constant thickness. Specific examples are given for both classes of problems.

Simmonds, James G.↗

Wavefront Sensing for WFIRST with a Linear Optical Model

In this paper we develop methods to use a linear optical model to capture the field dependence of wavefront aberrations in a nonlinear optimization-based phase retrieval algorithm for image-based wavefront sensing. The linear optical model is generated from a ray trace model of the system and allows the system state to be described in terms of mechanical alignment parameters rather than wavefront coefficients. This approach allows joint optimization over images taken at different field points and does not require separate convergence of phase retrieval at individual field points. Because the algorithm exploits field diversity, multiple defocused images per field point are not required for robustness. Furthermore, because it is possible to simultaneously fit images of many stars over the field, it is not necessary to use a fixed defocus to achieve adequate signal-to-noise ratio despite having images with high dynamic range. This allows high performance wavefront sensing using in-focus science data. We applied this technique in a simulation model based on the Wide Field Infrared Survey Telescope (WFIRST) Intermediate Design Reference Mission (IDRM) imager using a linear optical model with 25 field points. We demonstrate sub-thousandth-wave wavefront sensing accuracy in the presence of noise and moderate undersampling for both monochromatic and polychromatic images using 25 high-SNR target stars. Using these high-quality wavefront sensing results, we are able to generate upsampled point-spread functions (PSFs) and use them to determine PSF ellipticity to high accuracy in order to reduce the systematic impact of aberrations on the accuracy of galactic ellipticity determination for weak-lensing science.

Jurling, Alden S.↗

Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals

Recently, significant progress has been made in the handling of singular and nearly-singular potential integrals that commonly arise in the Boundary Element Method (BEM). To facilitate object-oriented programming and handling of higher order basis functions, cancellation techniques are favored over techniques involving singularity subtraction. However, gradients of the Newton-type potentials, which produce hypersingular kernels, are also frequently required in BEM formulations. As is the case with the potentials, treatment of the near-hypersingular integrals has proven more challenging than treating the limiting case in which the observation point approaches the surface. Historically, numerical evaluation of these near-hypersingularities has often involved a two-step procedure: a singularity subtraction to reduce the order of the singularity, followed by a boundary contour integral evaluation of the extracted part. Since this evaluation necessarily links basis function, Green s function, and the integration domain (element shape), the approach ill fits object-oriented programming concepts. Thus, there is a need for cancellation-type techniques for efficient numerical evaluation of the gradient of the potential. Progress in the development of efficient cancellation-type procedures for the gradient potentials was recently presented. To the extent possible, a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. However, since the gradient kernel involves singularities of different orders, we also require that the transformation leaves remaining terms that are analytic. The terms "normal" and "tangential" are used herein with reference to the source element. Also, since computational formulations often involve the numerical evaluation of both potentials and their gradients, it is highly desirable that a single integration procedure efficiently handles both.

Fink, Patrick W.↗

Quickbird Satellite in-orbit Modulation Transfer Function (MTF) Measurement Using Edge, Pulse and Impulse Methods for Summer 2003

The spatial characteristics of an imaging system cannot be expressed by a single number or simple statement. However, the Modulation Transfer Function (MTF) is one approach to measure the spatial quality of an imaging system. Basically, MTF is the normalized spatial frequency response of an imaging system. The frequency response of the system can be evaluated by applying an impulse input. The resulting impulse response is termed the Point Spread function (PSF). This function is a measure of the amount of blurring present in the imaging system and is itself a useful measure of spatial quality. An underlying assumption is that the imaging system is linear and shift-independent. The Fourier transform of the PSF is called the Optical Transfer Function (OTF) and the normalized magnitude of the OTF is the MTF. In addition to using an impulse input, a knife-edge in technique has also been used in this project. The sharp edge exercises an imaging system at all spatial frequencies. The profile of an edge response from an imaging system is called an Edge Spread Function (ESF). Differentiation of the ESF results in a one-dimensional version of the Point Spread Function (PSF). Finally, MTF can be calculated through use of Fourier transform of the PSF as stated previously. Every image includes noise in some degree which makes MTF of PSF estimation more difficult. To avoid the noise effects, many MTF estimation approaches use smooth numerical models. Historically, Gaussian models and Fermi functions were applied to reduce the random noise in the output profiles. The pulse-input method was used to measure the MTF of the Landsat Thematic Mapper (TM) using 8th order even functions over the San Mateo Bridge in San Francisco, California. Because the bridge width was smaller than the 30-meter ground sample distance (GSD) of the TM, the Nyquist frequency was located before the first zero-crossing point of the sinc function from the Fourier transformation of the bridge pulse. To avoid the zero-crossing points in the frequency domain from a pulse, the pulse width should be less than the width of two pixels (or 2 GSD's), but the short extent of the pulse results in a poor signal-to-noise ratio. Similarly, for a high-resolution satellite imaging system such as Quickbird, the input pulse width was critical because of the zero crossing points and noise present in the background area. It is important, therefore, that the width of the input pulse be appropriately sized. Finally, the MTF was calculated by taking ratio between Fourier transform of output and Fourier transform of input. Regardless of whether the edge, pulse and impulse target method is used, the orientation of the targets is critical in order to obtain uniformly spaced sub-pixel data points. When the orientation is incorrect, sample data points tend to be located in clusters that result in poor reconstruction of the edge or pulse profiles. Thus, a compromise orientation must be selected so that all spectral bands can be accommodated. This report continues by outlining the objectives in Section 2, procedures followed in Section 3, descriptions of the field campaigns in Section 4, results in Section 5, and a brief summary in Section 6.

Helder, Dennis↗

A Three-Dimensional Thermodynamically Based Function for the Progressive Failure of Unidirectional Composites

A thermodynamically-based work potential theory for modeling progressive damage for laminated, unidirectional composites assuming plane stress (2D Schaperystheory) is extended to three dimensioanl (3D). An internal state variable, S, is defined to account for the dissipated energy due to damage evolution in the form of microstructure changes in the matrix. With the stationary of the total work potential with respect to the internal state variable, a thermodynamically-consistent set of evolution equationsis derived. The internal state variable is related to the transverse and shear modulithrough microdamage functions. In the first part of this work, coupon specimensare prepared to conduct experiments to characterize the relations between the internalstate variable and the transverse modulus as well as shear modulus. The informationis subsequent used for the prediction of three point bending test. In the second partof this work, objectivity is studied. Three separate methods utilizing different definitions of a reduced internal state variable or of the order of the polynomials are used to represent the matrix micro-damage functions are employed. The three methods are implemented in a user defined subroutine within a commercial finite element method software package. Results from numerical simulations of a center-notched composites panel are compared. The agreement in the maximum stress predictions among the three methods indicates that objectivity, with respect to the functional form of themicrodamage functions, is satisfied.

carbon fiber reinforced polymer composites↗

Digital robust active control law synthesis for large order systems using constrained optimization

This paper presents a direct digital control law synthesis procedure for a large order, sampled data, linear feedback system using constrained optimization techniques to meet multiple design requirements. A linear quadratic Gaussian type cost function is minimized while satisfying a set of constraints on the design loads and responses. General expressions for gradients of the cost function and constraints, with respect to the digital control law design variables are derived analytically and computed by solving a set of discrete Liapunov equations. The designer can choose the structure of the control law and the design variables, hence a stable classical control law as well as an estimator-based full or reduced order control law can be used as an initial starting point. Selected design responses can be treated as constraints instead of lumping them into the cost function. This feature can be used to modify a control law, to meet individual root mean square response limitations as well as minimum single value restrictions. Low order, robust digital control laws were synthesized for gust load alleviation of a flexible remotely piloted drone aircraft.

Mukhopadhyay, Vivek↗

Quality of Candidate Flights and Submission Prediction in Collaborative Digital Departure Reroute

Collaborative Digital Departure Reroute (CDDR) enables the reroute of flights using a flight operator proposed set of alternative route options, referred to as Trajectory Option Set (TOS), in order to reduce delay on the airport's surface and in the Metroplex environment. The reroute functionality is enabled through NASA's Digital Information Platform (DIP). TOS candidate flights are defined as flights with an alternative route with delay savings greater than the flight operator defined relative trajectory cost. This paper analyzes the TOS candidate flights at Dallas/Fort Worth International Airport (KDFW) in the North Texas Metroplex to gain insight into which candidate flights are higher quality through a scoring method. This insight will inform refinements to help CDDR focus on high quality reroute opportunities. Binary classification models for predicting the flight operator's submission of candidate flights are also explored in this paper.

Machine Learning↗

Quality of Candidate Flights and Submission Prediction in Collaborative Digital Departure Reroute

Collaborative Digital Departure Reroute (CDDR) enables the reroute of flights using a flight operator proposed set of alternative route options, referred to as Trajectory Option Set (TOS), in order to reduce delay on the airport's surface and in the Metroplex environment. The reroute functionality is enabled through NASA's Digital Information Platform (DIP). TOS candidate flights are defined as flights with an alternative route with delay savings greater than the flight operator defined relative trajectory cost. This paper analyzes the TOS candidate flights at Dallas/Fort Worth International Airport (KDFW) in the North Texas Metroplex to gain insight into which candidate flights are higher quality through a scoring method. This insight will inform refinements to help CDDR focus on high quality reroute opportunities. Binary classification models for predicting the flight operator's submission of candidate flights are also explored in this paper.

Machine Learning↗

In-Flight Evaluation of the Traffic Aware Planner on the NASA HU-25A Guardian Aircraft

NASA’s Traffic Aware Planner (TAP) software is a research-prototype decision support tool that provides pilots with time- and fuel-saving route recommendations that optimize their current trajectory. The software runs on a first-of-a-kind system architecture onboard three aircraft in revenue service conducting operational evaluations with a major domestic airline. Therefore, significant NASA-internal testing is required prior to releasing the software to the partner airline. This paper describes a flight test plan that exercises the functionality of the TAP software in a representative operational environment, describes the system architecture developed and implemented for the NASA Langley HU-25A Guardian aircraft to support the test objectives, presents outcomes of the flight test campaign, and discusses use cases that demonstrate the value of flight testing for this activity.Research into flight path optimization of transport aircraft conducted by the National Aeronautics and SpaceAdministration (NASA) has produced an operational concept known as Traffic Aware Strategic Aircrew Requests(TASAR) [1, 2]. This near-term concept [3] provides the aircrew with a flight deck decision support tool known asthe Traffic Aware Planner (TAP). The TAP software leverages a growing number of information sources on the flightdeck to make time- and fuel-saving route optimization recommendations to the aircrew while en route. The aircrewcan then use the suggestions provided by the tool to make route change requests with a greater likelihood of acceptanceby air traffic control (ATC). Since TASAR is a concept intended for the current operational environment, it isintentionally designed to have no safety-critical impact or require any changes to current Federal AviationAdministration (FAA) rules and procedures [4, 5].The research prototype TAP system [6–8], explained further in Section III.C, continually incorporates up-to-dateaircraft state data from onboard avionics, as well as the latest position of surrounding traffic, the most recent windforecast, and the most recent convective weather forecast, in order to calculate candidate trajectory modifications thatimprove upon the current active route. These trajectories account for user-selectable objective functions [3] of reducedfuel burn, reduced flight time, or an airline-derived combination of factors known as trip cost. Previous analyses andsimulations have estimated substantial savings for airlines employing this technique within the U.S. National AirspaceSystem (NAS) [9–11]. Operational evaluations with Alaska Airlines seek to validate these projected benefits usingmeasured data while simultaneously providing benefits to the airline [12, 13].The TAP software has undergone a number of human-in-the-loop simulations [14] and flight test activities[15–17] in order to validate the operational concept, evaluate human factors considerations (e.g., workload, usability,distraction, etc.), and to assess the ability of the software to function in a representative operational environment (e.g.,connected to live avionics data, using in-flight internet connectivity, etc.). However, these simulations and flight testcampaigns did not account for the hardware architecture implemented on the three aircraft for Alaska Airlines’operational evaluations of the TAP software. Therefore, a need was identified to thoroughly test the functionality ofthe software in a similar hardware architecture to that of the partner airline’s aircraft. Information regarding testapparatus and environments used to evaluate TAP prior to testing on the HU-25A can be found in reference [18].A campaign of flight trials on a NASA aircraft, the HU-25A Guardian, was conducted to ensure that the researchprototype TAP system functions well in a configuration similar to the Alaska Airlines aircraft prior to deployment.This airborne, networked environment enables an assessment of the operational factors unique to the flight environment. Additionally, this activity evaluated the effectiveness and benefit of new TAP functionality andoperation in a relevant flight environment while allowing the rapid prototyping of new concepts and features.This paper is organized as follows: Section II discusses the details of the flight test plan, flight profiles, and theduties of personnel involved with conducting flight operations. Section III describes the test platform, avionicsequipage, and system architecture. Section IV presents a discussion of results, and Section V contains concludingremarks.

Underwood, Matthew C.↗

Order reduction of z-transfer functions via multipoint Jordan continued-fraction expansion

The order reduction problem of z-transfer functions is solved by using the multipoint Jordan continued-fraction expansion (MJCFE) technique. An efficient algorithm that does not require the use of complex algebra is presented for obtaining an MJCFE from a stable z-transfer function with expansion points selected from the unit circle and/or the positive real axis of the z-plane. The reduced-order models are exactly the multipoint Pade approximants of the original system and, therefore, they match the (weighted) time-moments of the impulse response and preserve the frequency responses of the system at some characteristic frequencies, such as gain crossover frequency, phase crossover frequency, bandwidth, etc.

Lee, Ying-Chin↗

Optimal model reduction and frequency-weighted extension

In this paper the quadratically optimal model reduction problem for single-input, single-output systems is considered. The reduced order model is determined by minimizing the integral of the magnitude-squared of the transfer function error. It is shown that the numerator coefficients of the optimal approximant satisfy a weighted least squares problem and, on this basis, a two-step iterative algorithm is developed combining a least squares solver with a gradient minimizer. The existence of globally optimal stable solutions to the optimization problem is established, and convergence of the algorithm to stationary values of the cost function is proved. The formulation is extended to handle the frequency-weighted optimal model reduction problem. Three examples demonstrate the optimization algorithm.

Spanos, J. T.↗

A new algorithm for L2 optimal model reduction

In this paper the quadratically optimal model reduction problem for single-input, single-output systems is considered. The reduced order model is determined by minimizing the integral of the magnitude-squared of the transfer function error. It is shown that the numerator coefficients of the optimal approximant satisfy a weighted least squares problem and, on this basis, a two-step iterative algorithm is developed combining a least squares solver with a gradient minimizer. Convergence of the proposed algorithm to stationary values of the quadratic cost function is proved. The formulation is extended to handle the frequency-weighted optimal model reduction problem. Three examples demonstrate the optimization algorithm.

Spanos, J. T.↗

Reduced-Order Modeling for Flutter/LCO Using Recurrent Artificial Neural Network

The present study demonstrates the efficacy of a recurrent artificial neural network to provide a high fidelity time-dependent nonlinear reduced-order model (ROM) for flutter/limit-cycle oscillation (LCO) modeling. An artificial neural network is a relatively straightforward nonlinear method for modeling an input-output relationship from a set of known data, for which we use the radial basis function (RBF) with its parameters determined through a training process. The resulting RBF neural network, however, is only static and is not yet adequate for an application to problems of dynamic nature. The recurrent neural network method [1] is applied to construct a reduced order model resulting from a series of high-fidelity time-dependent data of aero-elastic simulations. Once the RBF neural network ROM is constructed properly, an accurate approximate solution can be obtained at a fraction of the cost of a full-order computation. The method derived during the study has been validated for predicting nonlinear aerodynamic forces in transonic flow and is capable of accurate flutter/LCO simulations. The obtained results indicate that the present recurrent RBF neural network is accurate and efficient for nonlinear aero-elastic system analysis

Yao, Weigang↗

A technique for increasing the accuracy of the numerical inversion of the Laplace transform with applications

A technique is introduced which extends the range of useful approximation of numerical inversion techniques to many cycles of an oscillatory function without requiring either the evaluation of the image function for many values of s or the computation of higher-order terms. The technique consists in reducing a given initial value problem defined over some interval into a sequence of initial value problems defined over a set of subintervals. Several numerical examples demonstrate the utility of the method.

Berger, B. S.↗