Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear projection”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

A study of the use of linear programming techniques to improve the performance in design optimization problems

This project has two objectives. The first is to determine whether linear programming techniques can improve performance when handling design optimization problems with a large number of design variables and constraints relative to the feasible directions algorithm. The second purpose is to determine whether using the Kreisselmeier-Steinhauser (KS) function to replace the constraints with one constraint will reduce the cost of total optimization. Comparisons are made using solutions obtained with linear and non-linear methods. The results indicate that there is no cost saving using the linear method or in using the KS function to replace constraints.

Young, Katherine C.↗
NASA NTRS1993

A homotopy algorithm for digital optimal projection control GASD-HADOC

The linear-quadratic-gaussian (LQG) compensator was developed to facilitate the design of control laws for multi-input, multi-output (MIMO) systems. The compensator is computed by solving two algebraic equations for which standard closed-loop solutions exist. Unfortunately, the minimal dimension of an LQG compensator is almost always equal to the dimension of the plant and can thus often violate practical implementation constraints on controller order. This deficiency is especially highlighted when considering control-design for high-order systems such as flexible space structures. This deficiency motivated the development of techniques that enable the design of optimal controllers whose dimension is less than that of the design plant. A homotopy approach based on the optimal projection equations that characterize the necessary conditions for optimal reduced-order control. Homotopy algorithms have global convergence properties and hence do not require that the initializing reduced-order controller be close to the optimal reduced-order controller to guarantee convergence. However, the homotopy algorithm previously developed for solving the optimal projection equations has sublinear convergence properties and the convergence slows at higher authority levels and may fail. A new homotopy algorithm for synthesizing optimal reduced-order controllers for discrete-time systems is described. Unlike the previous homotopy approach, the new algorithm is a gradient-based, parameter optimization formulation and was implemented in MATLAB. The results reported may offer the foundation for a reliable approach to optimal, reduced-order controller design.

Collins, Emmanuel G., Jr.↗
NASA NTRS1987

Space power demonstrator engine, phase 1

The design, analysis, and preliminary test results for a 25 kWe Free-Piston Stirling engine with integral linear alternators are described. The project is conducted by Mechanical Technology under the direction of LeRC as part of the SP-100 Nuclear Space Power Systems Program. The engine/alternator system is designed to demonstrate the following performance: (1) 25 kWe output at a specific weight less than 8 kg/kW; (2) 25 percent efficiency at a temperature ratio of 2.0; (3) low vibration (amplitude less than .003 in); (4) internal gas bearings (no wear, no external pump); and (5) heater temperature/cooler temperature from 630 to 315 K. The design approach to minimize vibration is a two-module engine (12.5 kWe per module) in a linearly-opposed configuration with a common expansion space. The low specific weight is obtained at high helium pressure (150 bar) and high frequency (105 Hz) and by using high magnetic strength (samarium cobalt) alternator magnets. Engine tests began in June 1985; 16 months following initiation of engine and test cell design. Hydrotest and consequent engine testing to date has been intentionally limited to half pressure, and electrical power output is within 15 to 20 percent of design predictions.

Source record↗
NASA NTRS1993

Conjugate gradient based projection - A new explicit methodology for frictional contact

With special attention towards the applicability to parallel computation or vectorization, a new and effective explicit approach for linear complementary formulations involving a conjugate gradient based projection methodology is proposed in this study for contact problems with Coulomb friction. The overall objectives are focussed towards providing an explicit methodology of computation for the complete contact problem with friction. In this regard, the primary idea for solving the linear complementary formulations stems from an established search direction which is projected to a feasible region determined by the non-negative constraint condition; this direction is then applied to the Fletcher-Reeves conjugate gradient method resulting in a powerful explicit methodology which possesses high accuracy, excellent convergence characteristics, fast computational speed and is relatively simple to implement for contact problems involving Coulomb friction.

Tamma, Kumar K.↗
NASA NTRS2005

Optimal Prediction of Clocks from Finite Data

This talk is about optimal linear prediction of processes with stationary dth increments, which serve as a class of models for random clock disturbances. The predictor is obtained by orthogonal projection on the affine space of estimators whose errors are invariant to additive polynomials of degree < d. The projection conditions give a system of linear equations thatcan be solved straightforwardly for the regression coefficients. If the data are equally spaced, then the predictor can be obtained by an extension of Levinson's algorithm.

stationary increments↗
NASA NTRS1987

Highly Non-Linear Optical (NLO) organic crystals

This research project involves the synthesis and characterization of organic materials having powerful nonlinear optical (NLO) properties and the growth of highly ordered crystals and monomolecular films of these materials. Research in four areas is discussed: theoretical design of new materials, characterization of NLO materials, synthesis of new materials and development of coupling procedures for forming layered films, and improvement of the techniques for vapor phase and solution phase growth of high quality organic crystals. Knowledge gained from these experiments will form the basis for experiments in the growth of these crystals.

Harris, J. Milton↗
NASA NTRS1993

Summation by parts, projections, and stability

We have derived stability results for high-order finite difference approximations of mixed hyperbolic-parabolic initial-boundary value problems (IBVP). The results are obtained using summation by parts and a new way of representing general linear boundary conditions as an orthogonal projection. By slightly rearranging the analytic equations, we can prove strict stability for hyperbolic-parabolic IBVP. Furthermore, we generalize our technique so as to yield strict stability on curvilinear non-smooth domains in two space dimensions. Finally, we show how to incorporate inhomogeneous boundary data while retaining strict stability. Using the same procedure one can prove strict stability in higher dimensions as well.

Olsson, Pelle↗
NASA NTRS1987

Image reconstruction from multiple 1-D scans using filtered localized projection

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.

Frieden, B. Roy↗
DOE OSTI2026

Local reduced-order modeling for electrostatic plasmas by physics-informed solution manifold decomposition

Despite advancements in high-performance computing and modern numerical algorithms, computational cost remains prohibitive for multi-query kinetic plasma simulations. Here, in this work, we develop data-driven reduced-order models (ROMs) for collisionless electrostatic plasma dynamics, based on the kinetic Vlasov-Poisson equation. Our ROM approach projects the equation onto a linear subspace defined by the proper orthogonal decomposition (POD) modes. We introduce an efficient tensorial method to update the nonlinear term using a precomputed third-order tensor. We capture multiscale behavior with a minimal number of POD modes by decomposing the solution manifold into multiple time windows and creating temporally local ROMs. We consider two strategies for decomposition: one based on the physical time and the other based on the electric field energy. Applied to the 1D1V Vlasov–Poisson simulations, that is, prescribed E-field, Landau damping, and two-stream instability, we demonstrate that our ROMs accurately capture the total energy of the system both for parametric and time extrapolation cases. The temporally local ROMs are more efficient and accurate than the single ROM. In addition, in the two-stream instability case, we show that the energy-windowing reduced-order model (EW-ROM) is more efficient and accurate than the time-windowing reduced-order model (TW-ROM). With the tensorial approach, EW-ROM solves the equation approximately 90 times faster than Eulerian simulations while maintaining a maximum relative error of 7.5% for the training data and 11% for the testing data.

Electrostatic plasmas↗
DOE OSTI2020

Fermilab Cryogenic Test Facilities

Fermi National Accelerator Laboratory (Fermilab) has multiple cryogenic test facilities, undertaking testing of superconducting magnets, Superconducting Radio Frequency (SRF) cavities, SRF cryomodules and other helium cryogenic components. The test areas within Fermilab include: Meson Cryogenic Test facility (Meson), Industrial Building 1 (IB-1), Heavy Assembly Building/Illinois Accelerator Research Center (HAB/IARC), New Muon Lab, (NML) and the Cryomodule Test Facility (CMTF). Meson and HAB/IARC utilize repurposed Tevatron era reciprocating engine based cryogenic refrigerator systems, whereas CMTF and IB-1 utilize gas bearing turbine based coldboxes. Each of these test areas support the various Fermilab projects and collaborations including the Linear Coherent Light Source II (LCLS-II), Proton Improvement Plan II (PIP-II), High Luminosity Large Hadron Collider (HL-LHC), and Mu2e. This paper outlines the diverse and extensive cryogenic test capabilities within Fermilab.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Strong Turbulent Alfvén Wave Interactions in a Laboratory Plasma (Final Scientific/Technical Report)

Turbulence is a fundamental and widespread aspect of plasma behavior, occurring in nearly every natural plasma system that is observed, from clusters of distant galaxies to plasmas within the solar system. Plasma continuously streaming from the sun, known as the solar wind, has been particularly fruitful in understanding plasma turbulence due to the availability of spacecraft data. An ensuring mystery in solar wind turbulence is the unexpected observation of residual energy – there is excess energy in the magnetic fluctuations compared to the velocity fluctuations even though an Alfvén wave in Magnetohydrodynamics has equal amounts of energy in each. The project studied strong non-linear interaction between counter-propagating Alfvén waves in the lab for the first time. Dr. Mel Abler, who was hired as a postdoctoral researcher under this grant, obtained the first observation of residual energy in a counter-propagating Alfvén wave interaction. Dr. Abler showed experimentally that the interaction of two counter-propagating Alfvén waves produces daughter modes that are not not normal modes of the system and only exist in the presence of the nonlinear drive. These Alfvén quasimodes contain either negative residual energy (more energy in magnetic fluctuations) or positive residual energy (more energy in velocity fluctuations).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Component mode synthesis for model order reduction of nonclassically damped systems

Based on a component mode synthesis approach, a model-order-reduction method for linear structures with arbitrary linear damping has been developed. Projection matrices are introduced to make the method applicable to systems having rigid-body freedom. To test the method, eigenvalues of a reduced-order model of a free-free beam with nonproportional damping were compared to exact eigenvalues and to eigenvalues obtained using two other model-reduction strategies. The present model-reduction strategy proved to be decidedly superior.

Craig, Roy R., Jr.↗
NASA NTRS1992

Software For Generating Graphs And Charts

Common Graphics Library (CGL) computer program designed to enable users to generate graphs and charts of quality sufficient for publication or optical projection. Quickly and easily generates linear, logarithmic, bar, pie, and composite charts. Features of linear charts include automatic scaling, increasing or decreasing numerical axes, and character axes. Features of pie charts include segment labels, exploded segments, and chart keys. One of general features, horizontal or vertical orientation of figures. Written in ANSI FORTRAN 77.

Hammond, Dana P.↗
NASA NTRS1991

Integrated adaptive filtering and design for control experiments of flexible structures

A novel method is presented of identifying a state space model and a state estimator for linear stochastic systems from input and output data. The method is primarily based on the relations between the state space model and the finite difference model for linear stochastic systems derived through projection filters. It is proven that least squares identification of a finite difference model converges to the model derived from the projection filters. System pulse response samples are computed from the coefficients of the finite difference model. In estimating the corresponding state estimator gain, a z-domain method is used. First the deterministic component of the output is subtracted out, and then the state estimator gain is obtained by whitening the remaining signal. Experimental example is used to illustrate the feasibility of the method.

Huang, Jen-Kuang↗
NASA NTRS1992

Neural processing of gravity information

The goal of this project was to use the linear acceleration capabilities of the NASA Vestibular Research Facility (VRF) at Ames Research Center to directly examine encoding of linear accelerations in the vestibular system of the cat. Most previous studies, including my own, have utilized tilt stimuli, which at very low frequencies (e.g., 'static tilt') can be considered a reasonably pure linear acceleration (e.g., 'down'); however, higher frequencies of tilt, necessary for understanding the dynamic processing of linear acceleration information, necessarily involves rotations which can stimulate the semicircular canals. The VRF, particularly the Long Linear Sled, has promise to provide controlled pure linear accelerations at a variety of stimulus frequencies, with no confounding angular motion.

Schor, Robert H.↗
NASA NTRS1999

Use of the MATRIXx Integrated Toolkit on the Microwave Anisotropy Probe Attitude Control System

Recent advances in analytical software tools allow the analysis, simulation, flight code, and documentation of an algorithm to be generated from a single source, all within one integrated analytical design package. NASA's Microwave Anisotropy Probe project has used one such package, Integrated Systems' MATRIXx suite, in the design of the spacecraft's Attitude Control System. The project's experience with the linear analysis, simulation, code generation, and documentation tools will be presented and compared with more traditional development tools. In particular, the quality of the flight software generated will be examined in detail. Finally, lessons learned on each of the tools will be shared.

Ward, David K.↗