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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗
DOE OSTI2024

Multiparticle integral and differential correlation functions

This paper formalizes the use of integral and differential cumulants for measurements of multiparticle event-by-event transverse momentum fluctuations, rapidity fluctuations, as well as net-charge fluctuations. This enables the introduction of multiparticle balance functions, defined based on differential correlation functions (factorial cumulants), that suppress two- and three-prong resonance decays effects and enable measurements of underlying long-range correlations obeying quantum number conservation constraints. These multiparticle balance functions satisfy simple sum rules determined by quantum number conservation. It is additionally shown that these multiparticle balance functions arise as an intrinsic component of high-order net-charge cumulants. This implies that the magnitude of these cumulants, measured in a specific experimental acceptance, is strictly constrained by charge conservation and primarily determined by the rapidity and momentum width of these balance functions. The paper also presents techniques to reduce the computation time of differential correlation functions up to order n = 10 based on the methods of moments. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Adaptive Methods for Radial Basis Functions

Radial basis functions (RBFs) are a powerful tool for constructing high-order accurate reduced representations of scattered data in arbitrary dimension and on manifolds. We present a method of constructing data approximations in which we utilize a functional tail to capture a global background profile and a RBF neural network (NN) to capture the smaller-scale features. In the RBF NN the RBF centers, matrix shape parameters were selected adaptively for each RBF. We also utilized a geodesic notion of distance on the manifold on which the data lies, e.g., the spherical geodesic for data on the sphere. Although each of these ideas have been been investigated separately in previous works, their combination into a single algorithm is novel. We defined a machine learning problem in which these properties are learned to minimize the data reduction error. We demonstrate the algorithm for applications of scattered data reduction in the plane and on the sphere.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2021

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗
NASA NTRS1992

Surface transport vehicles and supporting technology requirements

Requirements have been identified for surface transport vehicles which allow remote scientific exploration on the moon, as well as lunar resource recovery and emplacement of a permanent base on the lunar surface. Attention is given to the results of a design study which developed configurational concepts for lunar surface transport vehicles and inferred technology-development requirements, with a view to a phased program of implementation. Distinct benefits are noted for the design of simple vehicle platforms with high commonality, in order to reduce logistical-support requirements and maximize functional flexibility. Two generic vehicle classed are defined.

Matijevic, J. R.↗
NASA NTRS2006

Coating MCPs with AlN and GaN

A development effort underway at the time of reporting the information for this article is devoted to increasing the sensitivity of microchannel plates (MCPs) as detectors of photons and ions by coating the MCPs with nitrides of elements in period III of the periodic table. Conventional MCPs are relatively insensitive to slowly moving, large-mass ions for example, ions of biomolecules under analysis in mass spectrometers. The idea underlying this development is to coat an MCP to reduce its work function (decrease its electron affinity) in order to increase both (1) the emission of electrons in response to impingement of low-energy, large-mass ions and (2) the multiplying effect of secondary electron emission. Of particular interest as coating materials having appropriately low or even negative electron affinities are gallium nitride, aluminum nitride, and ternary alloys of general composition Al(x)Ga(1-x)N (where 0<x<1). These materials exhibit attractively high degrees of chemical, mechanical, and thermal stability plus acceptably high resistance to sputtering. The electron-excitation cross sections of these materials are expected to exceed those of several other materials (including diamond) that are, variously, in use or under development for the same purpose. Moreover, by doping these materials with silicon, one can render them partly electrically conductive, thereby suppressing the undesired accumulation of electric charge that could otherwise occur during bombardment by ions. For experiments, thin films of AlN and GaN both undoped and doped with Si were deposited on commercial MCPs by radio-frequency molecular-beam epitaxy (also known as plasma-assisted molecular-beam epitaxy) at temperatures <200 C. This deposition technique is particularly suitable because (1) MCPs cannot withstand the higher deposition-substrate temperatures used to decompose constituent compounds in some other deposition techniques and (2) in this technique, the constituent Al, Ga, and N are supplied in elemental form, so that there is no need for thermal decomposition at the substrate surface. The nitride films thus formed were, variously, amorphous or polycrystalline. The nitride films were coated with surface layers of gold <100 thick. The MCPs were tested in a standard configuration in which the output stage of a first MCP was coupled to the input stage of a second MCP. Each pair of MCPs was mounted in a standard holder that included front and back contact rings and an anode for collecting the output electrons of the second MCP. The MCP pairs were biased at potentials between 1.7 and 1.9 kV, and count rates measured after preamplification and discrimination. To enable a direct comparison, in one pair, the second MCP was uncoated while the first MCP was coated over half its surface. The coated and uncoated sides of the half-coated MCP were exposed to fluxes of argon ions at kinetic energies of 1.0 and 0.5 keV. At 1.0 keV, the count rate for the coated side was about 2.3 times greater than that for the uncoated side; at 0.5 keV, the count rate for the coated side was about 1.8 times greater than that for the uncoated side.

Bensaoula, Abdelhakim↗
NASA NTRSDate not supplied

MODIS Reflective Solar Band Calibration Improvements using Pseudo-Invariant Desert Targets

To provide the best science data quality, an accurate characterization of the response versus scan angle (RVS) function is critical for the MODIS reflective solar bands (RSB) on-orbit calibration. In every MODIS operational scan, the Earth’s surface, referred to here as Earth view (EV), the space view (SV) port, and the onboard calibrators are viewed via a two-sided scan mirror. The RVS is defined as the sensor’s relative response as a function the angle of incidence (AOI) to the scan mirror. Many different approaches have been developed to derive the time-dependent RVS and its look-up table (LUT) applied to MODIS Level 1B (L1B) products since calibration Collection 4. For most MODIS RSB, the on-board calibrators can reasonably track the RVS change with time. In practice, their RVS is derived using data from on-board calibrators and the EV mirror side ratio (for mirror side 2). For Terra bands 1-4, 8-10 and Aqua bands 1-4, 8-9, an enhancement has been employed in Collections 6 and 6.1 (C6/C6.1) by using Earth scene response trending from pseudo-invariant desert sites in addition to the onboard calibrators. The current C6/C6/1 RVS algorithm is focused on fitting the EV data at each AOI over time and then deriving the relative change at different AOI. The EV response trending is currently fitted with multiple segments over time. Alternatively, the EV responses can be fit first as a function of AOI before fitting temporally in order to reduce the dependence on the stability of the desert site. These pre-treatment methods on the EV data provide improvement in the derived calibration coefficients. However, evidence of insufficient calibration is still observed in the MODIS L1B reflectance data, especially in the form of differences between the mirror sides. In this paper, we review the current methodologies that utilize the EV response trends from the pseudo-invariant Libyan desert targets to supplement the gain derived from the onboard calibrators. An improvement is then proposed and investigated such that a sliding window average (SWA) is used to pre-process the raw EV data. The SWA parameters are carefully selected using trade-off studies to accurately track the Earth scene response trending in multiple cases to overcome the reflectance differences between two mirror sides. Calibration results show improvements for both Aqua and Terra MODIS RSB L1B data products. This new adjustment has been included in the recently delivered Collection 7 LUT that will be evident in the L1B products expected to be released in late 2021.

MODIS↗
DOE OSTI2021

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Using dynamic mode decomposition to predict the dynamics of a two-time non-equilibrium Green’s function

Computing the numerical solution of the Kadanoff–Baym equations, a set of nonlinear integral differential equations satisfied by the two-time Green's functions derived from many-body perturbation theory for a quantum many-body system away from equilibrium, is a challenging task. Recently, we have successfully applied dynamic mode decomposition (DMD) to construct a data driven reduced order model that can be used to extrapolate the time-diagonal of a two-time Green's function from numerical solutions of the KBE within a small time window. In this paper, we extend the previous work and use DMD to predict off-diagonal elements of the two-time Green's function. We partition the two-time Green's function into a number of one-time functions along the diagonal and subdiagonals of the two-time window as well as in horizontal and vertical directions. We use DMD to construct separate reduced order models to predict the dynamics of these one-time functions in a two-step procedure. We extrapolate along diagonal and several subdiagonals within a subdiagonal band of a two-time window in the first step. In the second step, we use DMD to extrapolate the Green's function outside of the sub-diagonal band. In conclusion, we demonstrate the efficiency and accuracy of this approach by applying it to a two-band Hubbard model problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Reduced-Order Models of Static Power Grids based on Spectral Clustering

For large-scale interconnected power systems that cover large geographical areas, certain electrical studies are required so that appropriate decisions ensure system reliability and low cost. For such studies, it is often neither practical nor necessary to model in detail the entire power system, which is increasingly complex due to a more diverse range of grid assets to choose from in both short and long-term planning. The goal of this paper is to present a methodology to reduce the order of large-scale power networks based on spectral graph theory given that current methods for static network reduction are not scalable. A brief analysis of some spectral clustering properties to determine which graph Laplacian matrix should be used and why is included. The analysis shows that the utilization of the normalized graph Laplacian is more advantageous for clustering purposes. Techniques are proposed to approximate cost functions for the aggregated generators. This is done via linear regression. The reduced-order model obtained with the proposed methodology has an accuracy above 94% and solves the scalability issue commonly present in other reduction methods. If the utilization of the reduced-order model is either constrained to load levels above mid-peak demand, or cost functions of aggregated units are approximated via a piecewise quadratic approach, then the error distribution is in the order of 10^-3. .

Baquedano-Aguilar, Mario D.↗
DOE OSTI2022

Cost function for low-dimensional manifold topology assessment

Abstract In reduced-order modeling, complex systems that exhibit high state-space dimensionality are described and evolved using a small number of parameters. These parameters can be obtained in a data-driven way, where a high-dimensional dataset is projected onto a lower-dimensional basis. A complex system is then restricted to states on a low-dimensional manifold where it can be efficiently modeled. While this approach brings computational benefits, obtaining a good quality of the manifold topology becomes a crucial aspect when models, such as nonlinear regression, are built on top of the manifold. Here, we present a quantitative metric for characterizing manifold topologies. Our metric pays attention to non-uniqueness and spatial gradients in physical quantities of interest, and can be applied to manifolds of arbitrary dimensionality. Using the metric as a cost function in optimization algorithms, we show that optimized low-dimensional projections can be found. We delineate a few applications of the cost function to datasets representing argon plasma, reacting flows and atmospheric pollutant dispersion. We demonstrate how the cost function can assess various dimensionality reduction and manifold learning techniques as well as data preprocessing strategies in their capacity to yield quality low-dimensional projections. We show that improved manifold topologies can facilitate building nonlinear regression models.

42 ENGINEERING↗
NASA NTRS1975

Multidimensional Hermite-Gaussian quadrature formulae and their application to nonlinear estimation

A simplified technique is proposed for calculating multidimensional Hermite-Gaussian quadratures that involves taking the square root of a matrix by the Cholesky algorithm rather than computation of the eigenvectors of the matrix. Ways of reducing the dimension, number, and order of the quadratures are set forth. If the function f(x) under the integral sign is not well approximated by a low-order algebraic expression, the order of the quadrature may be reduced by factoring f(x) into an expression that is nearly algebraic and one that is Gaussian.

Mcreynolds, S. R.↗
NASA NTRS1989

Optimization of the structural and control system for LSS with reduced-order model

The objective is the simultaneous design of the structural and control system for space structures. The minimum weight of the structure is the objective function, and the constraints are placed on the closed loop distribution of the frequencies and the damping parameters. The controls approach used is linear quadratic regulator with constant feedback. A reduced-order control system is used. The effect of uncontrolled modes is taken into consideration by the model error sensitivity suppression (MESS) technique which modified the weighting parameters for the control forces. For illustration, an ACOSS-FOUR structure is designed for a different number of controlled modes with specified values for the closed loop damping parameters and frequencies. The dynamic response of the optimum designs for an initial disturbance is compared.

Khot, N. S.↗
NASA NTRSDate not supplied

Work Function Matching Passive Lunar Dust Mitigation Coating Preparation for Lunar Flight Opportunity

With the Artemis program, NASA isplanning longer stays on the surface, with more activities that have the potential to put the astronauts and equipment in contact with greater quantities of lunar dust. The success of these missions will depend on our understanding of material interactions with lunar dust and the development of ways to mitigate dust effects. This is particularly truein cases where exposure to dust will lead to failure of components, unacceptable loss of power or thermal control, unacceptable loss of visibility, or health issues.Optically transparent, sputter deposited,work function matching coatings are being developed at NASA Glenn Research Center to reduce adhesion of dust to windows, lenses and display panels by matching the minimum energy to remove an electron from the surface to that oflunar dust in order to reduce adhesion due to charge transfer. One of these work function matching coatings will be tested as part of Aegis Aerospace Inc.’s Regolith Adherence Characterization experiment going to the lunar surface on a Commercial Lunar Payload Services (CLPS) lander in 2023. Preparation of the work function matching coating and initial characterization prior to delivery for flight integration will be discussed.

Lunar dust, passive mitigation, CLPS demonstration↗
NASA NTRS1959

Trajectory Control for Vehicles Entering the Earth's Atmosphere at Small Flight Path Angles

Methods of controlling the trajectories of high-drag-low-lift vehicles entering the earth's atmosphere at angles of attack near 90 deg and at initial entry angles up to 3 deg are studied. The trajectories are calculated for vehicles whose angle of attack can be held constant at some specified value or can be perfectly controlled as a function of some measured quantity along the trajectory. The results might be applied in the design of automatic control systems or in the design of instruments which will give the human pilot sufficient information to control his trajectory properly during an atmospheric entry. Trajectory data are compared on the basis of the deceleration, range, angle of attack, and, in some cases, the rate of descent. The aerodynamic heat-transfer rate and skin temperature of a vehicle with a simple heat-sink type of structure are calculated for trajectories made with several types of control functions. For the range of entry angles considered, it is found that the angle of attack can be controlled to restrict the deceleration down to an arbitrarily chosen level of 3g. All the control functions tried are successful in reducing the maximum deceleration to the desired level. However, in order to avoid a tendency for the deceleration to reach an initial peak decrease, and then reach a second peak, some anticipation is required in the control function so that the change in angle of attack will lead the change in deceleration. When the angle of attack is controlled in the aforementioned manner, the maximum rate of aerodynamic heat transfer to the skin is reduced, the maximum skin temperature of the vehicle is virtually unaffected, and the total heat absorbed is slightly increased. The increase in total heat can be minimized, however, by maintaining the maximum desired deceleration for as much of the trajectory as possible. From an initial angle of attack of 90 deg, the angle-of-attack requirements necessary to maintain constant values of deceleration (1g to 4g) and constant values of rate of descent (450 to 1,130 ft/sec) as long as it is aerodynamically practical are calculated and are found to be moderate in both magnitude and rate. Entry trajectories made with these types of control are presented and discussed.

Eggleston, John M.↗
NASA NTRS1989

Autonomous frequency domain identification: Theory and experiment

The analysis, design, and on-orbit tuning of robust controllers require more information about the plant than simply a nominal estimate of the plant transfer function. Information is also required concerning the uncertainty in the nominal estimate, or more generally, the identification of a model set within which the true plant is known to lie. The identification methodology that was developed and experimentally demonstrated makes use of a simple but useful characterization of the model uncertainty based on the output error. This is a characterization of the additive uncertainty in the plant model, which has found considerable use in many robust control analysis and synthesis techniques. The identification process is initiated by a stochastic input u which is applied to the plant p giving rise to the output. Spectral estimation (h = P sub uy/P sub uu) is used as an estimate of p and the model order is estimated using the produce moment matrix (PMM) method. A parametric model unit direction vector p is then determined by curve fitting the spectral estimate to a rational transfer function. The additive uncertainty delta sub m = p - unit direction vector p is then estimated by the cross spectral estimate delta = P sub ue/P sub uu where e = y - unit direction vectory y is the output error, and unit direction vector y = unit direction vector pu is the computed output of the parametric model subjected to the actual input u. The experimental results demonstrate the curve fitting algorithm produces the reduced-order plant model which minimizes the additive uncertainty. The nominal transfer function estimate unit direction vector p and the estimate delta of the additive uncertainty delta sub m are subsequently available to be used for optimization of robust controller performance and stability.

Yam, Yeung↗
NASA NTRS2005

Reduced-Order Modeling of a Heaving Airfoil

A reduced-order model of a flapping airfoil is developed using Proper Orthogonal Decomposition (POD). The proper basis functions, developed from snapshots of full Navier-Stokes simulations, are used for a Galerkin projection of the governing equations. The resulting coupled, nonlinear ordinary di.erential equations have a low dimension because the first few basis members capture most of the energy of the flow. The reduced-order model is used to simulate heaving motions that are both similar to and different from the motion(s) used to generate the basis functions, and the errors in the model are quantified. Several methods are used to generate mode sets that can be used over a range of heaving parameters, including snapshots from one, two, and multiple Navier-Stokes simulations. As snapshots from additional simulations are added to the decomposition, the mode sets become richer and can simulate a wider range of parameter space, at some computational cost. Whereas the POD method is fully applicable in three dimensions, the simulation technique based on a body-fixed and body-fitted grid suffers large overhead when extended to three dimensions. To reduce the overhead, an embedding technique is discussed which embeds the solid wing into a fixed Cartesian grid. The wing, which can now have multiple pieces and also be flexible, is represented by a distribution of body forces. This distribution is determined to give exactly the flow around a flapping wing.

Haj-Hariri, H.↗