SEARCH · Engineering Papers

Results for “linear operator theory”

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

DOE OSTI2020

Zero-pairing and zero-temperature limits of finite-temperature Hartree-Fock-Bogoliubov theory

Recently, variational Hartree-Fock-Bogoliubov (HFB) mean-field equations were shown to possess a mathematically well-defined zero-pairing limit, independently of the closed- or open-shell character of the system under consideration. This limit is non-trivial for open-shell systems such that HFB theory does {\it not} reduce to the Hartree-Fock (HF) formalism in all cases. The present work extends this analysis to finite-temperature HFB (FTHFB) theory by investigating the behavior of this more general formalism in the combined zero-temperature and zero-pairing limits. The zero-pairing and zero-temperature limits of the FTHFB statistical density operator constrained to carry an arbitrary (integer) number of particles A on average is worked out analytically and realized numerically using a two-nucleon interaction. While the FTHFB density operator reduces to the projector corresponding to a pure HF Slater determinant for closed-shell nuclei, the FTHFB formalism does not reduce to the HF theory in all cases in the zero-temperature and zero-pairing limits, i.e. for open-shell nuclei. However, the fact that a nucleus can be of open-shell character in these joint limits is necessarily the result of some symmetry restrictions. Whenever it is the case, the non-trivial description obtained for open-shell systems is shown to depend on the order with which both limits are taken, i.e. the two limits do not commute for these systems. When the zero-temperature limit is performed first, the FTHFB density operator is demoted to a projector corresponding to a pure state made out of a linear combination of a finite number of Slater determinants with different (even) numbers of particles. When the zero-pairing limit is performed first, the FTHFB density operator remains a statistical mixture of a finite number of Slater determinants with both even and odd particle numbers. While the entropy (pairing density) is zero in the first (second) case, it does not vanish in the second (first) case in spite of the temperature (pairing) tending towards zero. The difference between both limits can have striking consequences for the (thermal) expectation values of observables. For instance, the particle-number variance does not vanish in either case and has limiting values that differ by a factor of two in both cases. In conclusion, while in the textbook situation associated with closed-shell nuclei Hartree-Fock-Bogoliubov (finite-temperature Hartree-Fock) theory reduces to Hartree-Fock theory in the zero-pairing (zero-temperature) limit, the present analysis demonstrates that a non trivial and unexpected limit is obtained for this formalism in open-shell systems. This result sheds a new light on certain aspects of this otherwise very well-studied many-body formalism.

20 ≤ A ≤ 38↗

DOE OSTI2022

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Quantum simulation of second-quantized Hamiltonians in compact encoding

We describe methods for simulating general second-quantized Hamiltonians using the compact encoding, in which qubit states encode only the occupied modes in physical occupation number basis states. These methods apply to second-quantized Hamiltonians composed of a constant number of interactions, i.e., linear combinations of ladder operator monomials of fixed form. Compact encoding leads to qubit requirements that are optimal up to logarithmic factors. Here, we show how to use sparse Hamiltonian simulation methods for second-quantized Hamiltonians in compact encoding, give explicit implementations for the required oracles, and analyze the methods. We also describe several example applications including the free boson and fermion theories, the φ 4 -theory, and the massive Yukawa model, all in both equal-time and light-front quantization. Our methods provide a general-purpose tool for simulating second-quantized Hamiltonians, with optimal or near-optimal scaling with error and model parameters.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1989

An algorithm for the solution of dynamic linear programs

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation scheme.

Psiaki, Mark L.↗

NASA NTRS1988

Crossed orbit interferometry - Theory and experimental results from SIR-B

Crossed orbit interferometry, which can perform measurements with only one antenna making two images of a scene during two separate passes and can operate even if the orbits are not parallel, is discussed and tested using SIR-B data. It is found that a Doppler refocusing of the SAR azimuth correlation, involving a resampling of one of the imgages in the cross-track direction, is necessary to remove the linear shift of the scene. The refocusing process also involves a terrain dependent resampling in the azimuth direction. A method for finding tie points to guide the resampling is discussed and a coarse altitude map derived only from the tie points is presented. Spatial heterodyned interferograms that contain the effects of the crossed orbit geometry are presented and a theoretical model is developed to explain them. The model is extended to calculate altitudes from the interferograms, and a final altitude map is presented.

Gabriel, Andrew K.↗

DOE OSTI2024

Parton distributions from boosted fields in the Coulomb gauge

We propose a new method to calculate parton distribution functions (PDFs) from lattice correlations of boosted quarks and gluons in the Coulomb gauge. Compared to the widely used gauge-invariant Wilson-line operators, these correlations greatly simplify the renormalization thanks to the absence of linear power divergence. Additionally, they enable access to larger off-axis momenta under preserved 3D rotational symmetry, as well as enhanced long-range precision that facilitates the Fourier transform. We verify the factorization formula that relates this new observable to the quark PDF at one-loop order in perturbation theory. Moreover, through a lattice calculation of the pion valence quark PDF, we demonstrate the aforementioned advantage and features of the Coulomb gauge correlation and show that it yields consistent results with the gauge-invariant method. This opens the door to a more efficient way to calculate parton physics on the lattice.

Baryons↗

NASA NTRS1951

Integrals and integral equations in linearized wing theory

The formulas of subsonic and supersonic wing theory for source, doublet, and vortex distributions are reviewed and a systematic presentation is provided which relates these distributions to the pressure and to the vertical induced velocity in the plane of the wing. It is shown that care must be used in treating the singularities involved in the analysis and that the order of integration is not always reversible. Concepts suggested by the irreversibility of order of integration are shown to be useful in the inversion of singular integral equations when operational techniques are used. A number of examples are given to illustrate the methods presented, attention being directed to supersonic flight speed.

Lomax, Harvard↗

NASA NTRS2008

Rotorcraft Technology for HALE Aeroelastic Analysis

Much of technology needed for analysis of HALE nonlinear aeroelastic problems is available from rotorcraft methodologies. Consequence of similarities in operating environment and aerodynamic surface configuration. Technology available - theory developed, validated by comparison with test data, incorporated into rotorcraft codes. High subsonic to transonic rotor speed, low to moderate Reynolds number. Structural and aerodynamic models for high aspect-ratio wings and propeller blades. Dynamic and aerodynamic interaction of wing/airframe and propellers. Large deflections, arbitrary planform. Steady state flight, maneuvers and response to turbulence. Linearized state space models. This technology has not been extensively applied to HALE configurations. Correlation with measured HALE performance and behavior required before can rely on tools.

Young, Larry↗

DOE OSTI2021

Stochastic Continuous-time Flexibility Scheduling and Pricing in Wholesale Electricity Markets

Large-scale integration of intermittent renewable energy sources (RES) is calling for additional flexibility resources as well as more advanced modeling and optimization techniques to account for the increasing uncertainty and variability in power systems operation. As the RES integration gains momentum, the magnitude and frequency of their variations increase, which may trigger ramping scarcity events in real-time power systems operation. This necessitates revisiting the present definition of power systems flexibility and reserve services to reflect their robustness and adequacy towards sub-interval variations of the load and RES, as well as adjusting the operation models to accommodate the new reserve services. This project took a fundamental approach and aimed at developing continuous-time scheduling and pricing model that accurately models the continuous-time variations of load and RES and efficiently deploys the ramping capability of flexible resources to compensate the sources of variability and uncertainty in the market. In this regard, this project pursued the following goals: Developing stochastic multi-fidelity continuous-time optimization models for scheduling of energy storage (ES) systems and flexible loads in wholesale energy markets; Developing the theory and practices of continuous-time locational marginal pricing for valuating energy storage systems and flexible loads in wholesale energy markets; Developing function space solution approach to convert the proposed stochastic multi-fidelity continuous-time optimization models into tractable mixed-integer linear optimization models; and Defining flexibility reserve as a new type of reserve in markets that would enable ultimate participation of energy storage devices in provision of services to compensate the variability and uncertainty of RES in electricity markets. This project successfully completed all five major tasks defined in the SOPO, and produced 8 high-impact journal papers, 6 conference papers, 3 published U.S. patents, and one web-based software for continuous-time operation optimization of power systems. The application of the proposed flexibility reserve and the stochastic multi-fidelity continuous-time operation scheduling models would modify the forward commitment and schedule of generating units, ES devices and flexible loads, and would line up the resources in such a way that the composition of available resources is better prepared to respond to the sub-hourly variations of the load and renewable resources in real-time operation. Therefore, this project paves the way to sustainable, reliable, and economic integration of renewable energy resources in power system, supporting the progress towards reaching the national targets on energy independence. Even if the proposed models offers a radically different point of view as compared to existing models, it does not alter fundamentally the architecture of power systems operations, nor the complexity of the scheduling problem, so the integration of this project in power systems is extremely practical.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Novel Relaxation Time Approximation to the Relativistic Boltzmann Equation

In this work, we show that the widely used relaxation time approximation to the relativistic Boltzmann equation contains basic flaws, being incompatible with micro- and macroscopic conservation laws if the relaxation time depends on energy or general matching conditions are applied. We propose a new approximation that fixes such fundamental issues and maintains the basic properties of the linearized Boltzmann collision operator. We show how this correction affects transport coefficients, such as the bulk viscosity and particle diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1973

Hydrogen evolution in nickel-water heat pipes.

A study was made of the evolution of hydrogen gas in nickel-water heat pipes for the purpose of investigating methods of accelerated life testing. The data were analyzed in terms of a phenomenological corrosion model of heat pipe degradation which incorporates corrosion and oxidation theory and contains parameters which can be determined by experiment. The gas was evolved with a linear time dependence and an exponential temperature dependence with an activation energy of 1.03 x 10 to the minus 19th joules. A flow-rate dependence of the gas evolution was found in the form of a threshold. The results were used to predict usable lifetimes of heat pipes operated at normal operating conditions from results taken under accelerated operating conditions.

Anderson, W. T.↗

NASA NTRS2003

Transient Modeling of Hybrid Rocket Low Frequency Instabilities

A comprehensive dynamic model of a hybrid rocket has been developed in order to understand and predict the transient behavior including instabilities. A linearized version of the transient model predicted the low-frequency chamber pressure oscillations that are commonly observed in hybrids. The source of the instabilities is based on a complex coupling of thermal transients in the solid fuel, wall heat transfer blocking due to fuel regression rate and the transients in the boundary layer that forms on the fuel surface. The oscillation frequencies predicted by the linearized theory are in very good agreement with 43 motor test results obtained from the hybrid propulsion literature. The motor test results used in the comparison cover a very wide spectrum of parameters including: 1) four separate research and development programs, 2) three different oxidizers (LOX, GOX, N2O), 3) a wide range of motor dimensions (i.e. from 5 inch diameter to 72 inch diameter) and operating conditions and 4) several fuel formulations. A simple universal scaling formula for the frequency of the primary oscillation mode is suggested.

Karabeyoglu, M. Arif↗

NASA NTRS1972

Application of a stochastic inverse to the geophysical inverse problem

The inverse problem for gross earth data can be reduced to an undertermined linear system of integral equations of the first kind. A theory is discussed for computing particular solutions to this linear system based on the stochastic inverse theory presented by Franklin. The stochastic inverse is derived and related to the generalized inverse of Penrose and Moore. A Backus-Gilbert type tradeoff curve is constructed for the problem of estimating the solution to the linear system in the presence of noise. It is shown that the stochastic inverse represents an optimal point on this tradeoff curve. A useful form of the solution autocorrelation operator as a member of a one-parameter family of smoothing operators is derived.

Jordan, T. H.↗

NASA NTRS1982

Simultaneous stabilization and simultaneous pole placement by nonswitching dynamic compensation

The 'simultaneous stabilization problem' is defined and theorems are proposed for its solution. The problem consists in answering the question: given an r-tuple G sub 1(s), G sub r(s) of p x m proper transfer functions, does there exist a compensator K(s) such that the closed loop systems G sub 1(s) (I+K(s)G sub 1(s)) (-1), G sub r(s) (I+K(s) G sub r(s)) (-1) are (internally) stable. This question arises in reliability theory, where G sub 2(s), G sub r(s) represents a plant G sub 1(s) operating in various modes of failure and K(s) is a nonswitching stabilizing compensator. It is important in the stability analysis and design of a plant which can be switched into various operating modes. The simultaneous stabilization problem can also apply to the stabilization of a nonlinear system which is linearized at several equilibria. Conditions are defined for pole placement and the generalized Sylvestor matrix is discussed.

Ghosh, B. K.↗

NASA NTRS1983

Simultaneous stabilization and simultaneous pole-placement by nonswitching dynamic compensation

Ghosh, B. K.↗

NASA NTRS1980

Nonlinear state estimation and feedback control of nonlinear and bilinear distributed parameter systems

This paper presents a theory of nonlinear state observers for nonlinear and bilinear distributed parameter systems. Convergence results are proved for these observers. Linear feedback control derived from such state observers is applied to the distributed parameter system and conditions are presented for closed-loop stability. The emphasis is on finite dimensional state observers and controllers (which can be implemented with on-line computers) and conditions for their successful operation with infinite dimensional distributed parameter systems.

Balas, M. J.↗

NASA NTRS1992

Adaptive control of Space Station with control moment gyros

An adaptive approach to Space Station attitude control is investigated. The main components of the controller are the parameter identification scheme, the control gain calculation, and the control law. The control law is a full-state feedback space station baseline control law. The control gain calculation is based on linear-quadratic regulator theory with eigenvalues placement in a vertical strip. The parameter identification scheme is a recursive extended Kalman filter that estimates the inertias and also provides an estimate of the unmodeled disturbances due to the aerodynamic torques and to the nonlinear effects. An analysis of the inertia estimation problem suggests that it is possible to estimate Space Station inertias accurately during nominal control moment gyro operations. The closed-loop adaptive control law is shown to be capable of stabilizing the Space Station after large inertia changes. Results are presented for the pitch axis.

Bishop, Robert H.↗

NASA NTRS1988

A new approach to the oxidational theory of mild wear

The wear and friction of 316 stainless steel was investigated in relation to oxidation at the real area of contact. An approach describing the oxidation characteristics of mild wear in terms of linear oxidation kinetics is presented. Linear oxidation kinetics more closely approaches the physical conditions involved where oxide growth on asperities is continuously interrupted by spalling and removal. The mild wear rate may now be described mathematically by two empirical equations in terms of the operational variables and the oxidational constants without using the critical oxide thickness, which is dependent on asperity size, shape and temperature, and is difficult to measure. The resulting activation energy for oxidational wear was significantly lower than that for static oxidation governed by diffusion controlled kinetics.

Hong, H.↗