SEARCH · Search NASA

Results for “Model order reduction”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1993

The analytical representation of viscoelastic material properties using optimization techniques

This report presents a technique to model viscoelastic material properties with a function of the form of the Prony series. Generally, the method employed to determine the function constants requires assuming values for the exponential constants of the function and then resolving the remaining constants through linear least-squares techniques. The technique presented here allows all the constants to be analytically determined through optimization techniques. This technique is employed in a computer program named PRONY and makes use of commercially available optimization tool developed by VMA Engineering, Inc. The PRONY program was utilized to compare the technique against previously determined models for solid rocket motor TP-H1148 propellant and V747-75 Viton fluoroelastomer. In both cases, the optimization technique generated functions that modeled the test data with at least an order of magnitude better correlation. This technique has demonstrated the capability to use small or large data sets and to use data sets that have uniformly or nonuniformly spaced data pairs. The reduction of experimental data to accurate mathematical models is a vital part of most scientific and engineering research. This technique of regression through optimization can be applied to other mathematical models that are difficult to fit to experimental data through traditional regression techniques.

Hill, S. A.↗

NASA NTRS2002

Acoustic and Laser Doppler Anemometer Results for Confluent and 12-Lobed E(exp 3) Mixer Exhaust Systems for Subsonic Jet Noise Reduction

The research described in this report has been funded by NASA Glenn Research Center as part of the Advanced Subsonic Technologies (AST) initiative. The program operates under the Large Engine Technologies (LET) as Task Order #3 1. Task Order 31 is a three year research program divided into three subtasks. Subtask A develops the experimental acoustic and aerodynamic subsonic mixed flow exhaust system databases. Subtask B seeks to develop and assess CFD-based aero-acoustic methods for subsonic mixed flow exhaust systems. Subtask B relies on the data obtained from Subtask A to direct and calibrate the aero-acoustic methods development. Subtask C then seeks to utilize both the aero-acoustic data bases developed in Subtask A and the analytical methods developed in Subtask B to define improved subsonic mixed-flow exhaust systems. The mixed flow systems defined in Subtask C will be experimentally demonstrated for improved noise reduction in a scale model aero-acoustic test conducted similarly to the test performed in Subtask A. The overall object of this Task Order is to develop and demonstrate the technology to define a -3EPNdB exhaust system relative to 1992 exhaust system technology.

Salikuddin, M.↗

NASA NTRS1998

Flammability of Heterogeneously Combusting Metals

Most engineering materials, including some metals, most notably aluminum, burn in homogeneous combustion. 'Homogeneous' refers to both the fuel and the oxidizer being in the same phase, which is usually gaseous. The fuel and oxidizer are well mixed in the combustion reaction zone, and heat is released according to some relation like q(sub c) = delta H(sub c)c[((rho/rho(sub 0))]exp a)(exp -E(sub c)/RT), Eq. (1) where the pressure exponent a is usually close to unity. As long as there is enough heat released, combustion is sustained. It is useful to conceive of a threshold pressure beyond which there is sufficient heat to keep the temperature high enough to sustain combustion, and beneath which the heat is so low that temperature drains away and the combustion is extinguished. Some materials burn in heterogeneous combustion, in which the fuel and oxidizer are in different phases. These include iron and nickel based alloys, which burn in the liquid phase with gaseous oxygen. Heterogeneous combustion takes place on the surface of the material (fuel). Products of combustion may appear as a solid slag (oxide) which progressively covers the fuel. Propagation of the combustion melts and exposes fresh fuel. Heterogeneous combustion heat release also follows the general form of Eq.(1), except that the pressure exponent a tends to be much less than 1. Therefore, the increase in heat release with increasing pressure is not as dramatic as it is in homogeneous combustion. Although the concept of a threshold pressure still holds in heterogeneous combustion, the threshold is more difficult to identify experimentally, and pressure itself becomes less important relative to the heat transfer paths extant in any specific application. However, the constants C, a, and E(sub c) may still be identified by suitable data reduction from heterogeneous combustion experiments, and may be applied in a heat transfer model to judge the flammability of a material in any particular actual-use situation. In order to support the above assertions, two investigations are undertaken: 1) PCT data are examined in detail to discover the pressure dependence of heterogeneous combustion experiment results; and 2) heterogeneous combustion in a PCT situation is described by a heat transfer model, which is solved first in simplified form for a simple actual-use situation, and then extended to apply to PCT data reduction (combustion constant identification).

Jones, Peter D.↗

NASA NTRS1992

An unsymmetric Lanczos algorithm for damped structural dynamics systems

A one-sided, unsymmetric block Lanczos algorithm is proposed for the model reduction of structural dynamics systems with unsymmetric damping and/or stiffness matrices. The algorithm is a three-term iteration scheme, which transforms the system matrix into an almost skew-symmetric, block-tridiagonal form. The Lanczos reduced-order model is guaranteed to be stable if the full-order system is stable. For unstable systems, a shifting method is available. Also, the algorithm offers flexibility in the choice of starting vectors and thus can yield more accurate reduced-order models. A linear system example and a plane truss structure example are used to show the efficacy of the proposed method.

Su, Tzu-Jeng↗

NASA NTRS2011

Waterhammer Transient Simulation and Model Anchoring for the Robotic Lunar Lander Propulsion System

Waterhammer transients have the potential to adversely impact propulsion system design if not properly addressed. Waterhammer can potentially lead to system plumbing, and component damage. Multi-thruster propulsion systems also develop constructive/destructive wave interference which becomes difficult to predict without detailed models. Therefore, it is important to sufficiently characterize propulsion system waterhammer in order to develop a robust design with minimal impact to other systems. A risk reduction activity was performed at Marshall Space Flight Center to develop a tool for estimating waterhammer through the use of anchored simulation for the Robotic Lunar Lander (RLL) propulsion system design. Testing was performed to simulate waterhammer surges due to rapid valve closure and consisted of twenty-two series of waterhammer tests, resulting in more than 300 valve actuations. These tests were performed using different valve actuation schemes and three system pressures. Data from the valve characterization tests were used to anchor the models that employed MSCSoftware.EASY5 v.2010 to model transient fluid phenomena by using transient forms of mass and energy conservation. The anchoring process was performed by comparing initial model results to experimental data and then iterating the model input to match the simulation results with the experimental data. The models provide good correlation with experimental results, supporting the use of EASY5 as a tool to model fluid transients and provide a baseline for future RLL system modeling. This paper addresses tasks performed during the waterhammer risk reduction activity for the RLL propulsion system. The problem of waterhammer simulation anchoring as applied to the RLL system is discussed with results from the corresponding experimental valve tests. Important factors for waterhammer mitigation are discussed along with potential design impacts to the RLL propulsion system.

Stein, William B.↗

NASA NTRS1984

Stellar fibril magnetic systems. I - Reduced energy state

The remarkable fibril structure of the magnetic fields at the surface of the sun (with fibrils compressed to 1,000-2,000 gauss) lies outside existing statistical theories of magnetohydrodynamic turbulence. The total energy of the fibril field is enhanced by a factor of more than 100 above the energy for the mean field in a continuum state. The magnetic energy density within a fibril is of the order of 100 times the local kinetic energy density, so that no simple application of equipartition principles is possible. It is pointed out that the total energy of the atmosphere (thermal + gravitational + magnetic) is reduced by the fibril state of the field by avoiding the magnetic inhibition of the convective overturning, suggesting that the formation of the observed intense fibril state may be in response to the associated energy reduction. Calculation of the minimum total energy of a polytropic atmosphere permeated by magnetic fibrils yields theoretical fibril fields of the order of 1-5 kilogauss when characteristics appropriate to the solar convective zone are introduced, in rough agreement with the actual fields of 1-2 kilogauss. The polytrope model, although crude, establishes that a large reduction in total energy is made possible by the fibril state.

Parker, E. N.↗

NASA NTRS1993

Higher-Order Ionospheric Effects on the GPS Observables and Means of Modeling Them

Based on realistic modeling of the electron density of the ionosphere and using a dipole moment approximation for the earth magnetic field, we are able to estimate the effect of the ionosphere on the GPS signal for a ground user. Our study implies reduction of the second-order ionospheric error to less than 2 mm on the average and, therefore, potential improvement of orbit determination and baseline solutions.

ionosphere↗

DOE OSTI2026

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗

DOE OSTI2025

An evaluation of multi-fidelity methods for quantifying uncertainty in projections of ice-sheet mass change

Abstract. This study investigated the computational benefits of using multi-fidelity statistical estimation (MFSE) algorithms to quantify uncertainty in the mass change of Humboldt Glacier, Greenland, between 2007 and 2100 using a single climate change scenario. The goal of this study was to determine whether MFSE can use multiple models of varying cost and accuracy to reduce the computational cost of estimating the mean and variance of the projected mass change of a glacier. The problem size and complexity were chosen to reflect the challenges posed by future continental-scale studies while still facilitating a computationally feasible investigation of MFSE methods. When quantifying uncertainty introduced by a high-dimensional parameterization of the basal friction field, MFSE was able to reduce the mean-squared error in the estimates of the statistics by well over an order of magnitude when compared to a single-fidelity approach that only used the highest-fidelity model. This significant reduction in computational cost was achieved despite the low-fidelity models used being incapable of capturing the local features of the ice-flow fields predicted by the high-fidelity model. The MFSE algorithms were able to effectively leverage the high correlation between each model's predictions of mass change, which all responded similarly to perturbations in the model inputs. Consequently, our results suggest that MFSE could be highly useful for reducing the cost of computing continental-scale probabilistic projections of sea-level rise due to ice-sheet mass change.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1987

The effect of a gravity wave drag parameterization scheme on GLA fourth order GCM forecasts

Ten-day forecast experiments have been performed to determined whether the introduction of a simple orographic gravity wave drag scheme into a fourth-order GCM would reduce the climate drift of the fine resolution model and improve the model's medium range predictive skill. Error reduction due to the gravity waves is found in stratospheric predictions, where the improvement is confined mainly to the zonal mean component. Improvements are noted in the Northern Hemisphere climatology, where low level westerlies are weakened and shifted poleward, and in the Southern Hemisphere, where the roaring forties and fifties are better simulated.

Helfand, H. M.↗

NASA NTRS2000

Gravity and Nonconservative Force Model Tuning for the GEOSAT Follow-On Spacecraft

The US Navy's GEOSAT Follow-On spacecraft was launched on February 10, 1998 and the primary objective of the mission was to map the oceans using a radar altimeter. Three radar altimeter calibration campaigns have been conducted in 1999 and 2000. The spacecraft is tracked by satellite laser ranging (SLR) and Doppler beacons and a limited amount of data have been obtained from the Global Positioning Receiver (GPS) on board the satellite. Even with EGM96, the predicted radial orbit error due to gravity field mismodelling (to 70x70) remains high at 2.61 cm (compared to 0.88 cm for TOPEX). We report on the preliminary gravity model tuning for GFO using SLR, and altimeter crossover data. Preliminary solutions using SLR and GFO/GFO crossover data from CalVal campaigns I and II in June-August 1999, and January-February 2000 have reduced the predicted radial orbit error to 1.9 cm and further reduction will be possible when additional data are added to the solutions. The gravity model tuning has improved principally the low order m-daily terms and has reduced significantly the geographically correlated error present in this satellite orbit. In addition to gravity field mismodelling, the largest contributor to the orbit error is the non-conservative force mismodelling. We report on further nonconservative force model tuning results using available data from over one cycle in beta prime.

Lemoine, Frank G.↗

DOE OSTI2024

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗

NASA NTRS1982

Pole placement and order reduction in two-time-scale control systems through Riccati iteration

A transformation of variables taken from singular perturbations may be applied to two-time-scale linear systems in state space form to reduce the system to block-diagonal form with slow and fast modes decoupled. The transformation is easily computed by applying the new Riccati iteration. The iteration yields a solution to the nonsymmetric algebraic Riccati equation obtained by partitioning the original system matrix A. The numerical procedure is initiated with the trivial iterate L(0) = 0, and is globally convergent to the desired unique time scale decoupling solution. After transformation, the decoupled system may be used in controller design to achieve exact closed loop pole placement in the slow subsystem without altering the poles of the fast subsystem. The decoupled form may also be used to reduce system order by wetting a small parameter to zero. Provided the fast subsystem is stable, the order reduction can be expected to yield a good approximation to the original system. These methods are demonstrated using the 16th order linear model of a turbofan engine.

Anderson, L. R.↗

NASA NTRSDate not supplied

Microstructure Controlled Multi-Layer Porous Material Liner Tested On The Advanced Noise Control Fan

A novel concept of multilayer micro structurally controlled noise mitigation porous material was developed and experimentally validated in the laboratory by normal incidence impedance test. Extensive parametric studies were made to understand various parameters that influence acoustic absorption and transmission loss based on the test specimen size, particle size, porosity, pore size, and its distribution through the thickness. The results were reported in many publications. This work was focused on applying the concept and assessing the validity of the technology to aircraft turbomachine engines. The NASA Glenn Research Center’s Advanced Noise Control Fan (ANCF) low-speed test rig was selected. The liner insert required a 4-ft internal diameter cylinder with a 2-inch thick and 9-inch axial length. The liner was fabricated in eight segments and assembled in an aluminum test rig. The liner was tested for both tonal and broadband noise at the Notre Dame Turbomachinery Lab White Field test facility, where the ANCF is located. The study concludes that: Microstructurally controlled porous material is very effective in dissipating the sound energy; the liner is effective for tonal noise (2 to 13 dB reduction in noise levels) that constitutes the major part of the overall fan noise, especially the blade passing frequency which is generally the dominating component; Also, the liner is effective for broadband noise with reductions in sound pressure level of about 3 dB for the shaft order greater than 20. The proposed two-layer liner model demonstrated the validity of the concept. More studies are required to optimize the design.

Liner↗

NASA NTRS1995

Intermittent turbulence in solar wind from the south polar hole

The magnetic fields measured by the Ulysses spacecraft are used to study solar wind turbulence in the fast solar wind from the south polar hole. The spacecraft was at about 46 deg south latitude and 3.9 AU. For a magnetic field with a Gaussian distribution the power spectrum (second-order structure function) is sufficient to completely characterize the turbulence. However, the actual distribution is non-Gaussian so that the effects of intermittency must be taken into account. The observed spectral exponents include effects of intermittency and cannot be directly compared with the standard second-order spectral theories such as the Kolmogorov and Kraichnan theories. To permit a better comparison of the observations with the theoretical models, we study the structure characteristics of the data. We find the exponents of the second-order structure functions (power spectra) and the higher-order normalized structure functions for the components of the magnetic fields. We show that these sets of exponents can be approximately described by two basic numbers: the spectral exponent and the intermittency exponent. The intermittency exponent characterizes correlation properties of the energy cascade from large to small scales. Before comparing the observations to the theoretically expected values, a reduction must be made to the observed spectral exponent. The amount of the reduction depends on both the intermittency exponent and the model of the energy cascade assumed in the turbulence theory. We reduce the measured spectral indices according to a simple model for Alfven turbulence that is described here. We then compare our reduced spectral indices with second-order spectral theory. The reduced spectral indices for the period range of 1 min to about a half hour are remarkably constant and in good agreement with the value of 3/2. Thus our treatment is self-consistent. Our tentative conclusion is that the high-frequency turbulence appears to agree with the model of random-phased Alfven waves. This tentative conclusion must be tested by further theoretical and observational work.

Ruzmaikin, A. A.↗

DOE OSTI2024

Exponential Improvements in the Simulation of Lattice Gauge Theories Using Near-Optimal Techniques

We report a first-of-its-kind analysis on post-Trotter simulation of U(1), SU(2), and SU(3) lattice gauge theories including fermions in arbitrary spatial dimension. We provide explicit circuit constructions as well as T-gate counts and logical qubit counts for Hamiltonian simulation. We find a reduction of up to 25 orders of magnitude in space-time volume over Trotter methods for simulations of non-Abelian lattice gauge theories relevant to the standard model. This improvement results from our algorithm having polynomial scaling with the number of colors in the gauge theory, achieved by utilizing oracle constructions relying on the sparsity of physical operators, in contrast to the exponential scaling seen in state-of-the-art Trotter methods, which employ explicit mappings onto Pauli operators. Our work demonstrates that the use of advanced algorithmic techniques leads to dramatic reductions in the cost of simulating fundamental interactions, bringing it in step with resources required for first-principles quantum simulation of chemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

(2 + δ)-dimensional theory of the electromechanics of lipid membranes: Electrostatics

The coupling of electric fields to the mechanics of lipid membranes gives rise to intriguing electromechanical behavior, as, for example, evidenced by the deformation of lipid vesicles in external electric fields. Electromechanical effects are relevant for many biological processes, such as the propagation of action potentials in axons and the activation of mechanically gated ion channels. Currently, a theoretical framework describing the electromechanical behavior of arbitrarily curved and deforming lipid membranes does not exist. Purely mechanical models commonly treat lipid membranes as two-dimensional surfaces, ignoring their finite thickness. While holding analytical and numerical merit, this approach cannot describe the coupling of lipid membranes to electric fields and is thus unsuitable for electromechanical models. In a sequence of articles, we derive an effective surface theory of the electromechanics of lipid membranes, called the (2 + δ)-dimensional theory, which has the advantages of surface descriptions while accounting for finite thickness effects. The present article proposes a generic dimension reduction procedure relying on low-order spectral expansions. This procedure is applied to the electrostatics of lipid membranes to obtain the (2 + δ)-dimensional theory that captures potential differences across and electric fields within lipid membranes. Finally, this model is tested on different geometries relevant for lipid membranes, showing good agreement with the corresponding three-dimensional electrostatics theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1983

Model deformation system

The development of a system to measure model deflections encountered in the National Transonic Facility is discussed. The goal is to be able to measure peak deflections of up to 3 in. with accuracies to within 0.0025 in. over an area 1 m square as the model pitches through an included angle of 30 deg. Stereophotogrammetric techniques are being implemented, with the initial system being an extension of standard techniques. A second system, which will be all electronic, is under development. Both techniques require targets to be strategically placed on the model. Active targets are being developed for location in the model in order to maximize the signal-to-noise ratio and to approximate a point source. Image processing techniques and stereophotogrammetric data reduction programs are being implemented to perform the data reduction tasks.

Holmes, H. K.↗