Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear algebra”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Optimizing Power–Frequency Droop Characteristics of Distributed Energy Resources

This paper outlines a procedure to design power-frequency droop slopes for distributed energy resources (DERs) installed in distribution networks to optimally participate in primary frequency response. In particular, the droop slopes are engineered such that DERs respond in proportion to their power ratings and they are not unfairly penalized in power provisioning based on their location in the distribution network. The main contribution of our approach is that a guaranteed level of frequency regulation can be guaranteed at the feeder head, while ensuring that the outputs of individual DERs conform to some well-defined notion of fairness. The approach we adopt leverages an optimization-based perspective and suitable linearizations of the power-flow equations to embed notions of fairness and information regarding the physics of the power flows within the distribution network into the droop slopes. Time-domain simulations from a differential algebraic equation model of the 39-bus New England test-case system augmented with three instances of the IEEE 37-node distribution-network with frequency-sensitive DERs are provided to validate our approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1978

Aircraft model prototypes which have specified handling-quality time histories

Several techniques for obtaining linear constant-coefficient airplane models from specified handling-quality time histories are discussed. The pseudodata method solves the basic problem, yields specified eigenvalues, and accommodates state-variable transfer-function zero suppression. The algebraic equations to be solved are bilinear, at worst. The disadvantages are reduced generality and no assurance that the resulting model will be airplane like in detail. The method is fully illustrated for a fourth-order stability-axis small motion model with three lateral handling quality time histories specified. The FORTRAN program which obtains and verifies the model is included and fully documented.

Johnson, S. H.↗
DOE OSTI2011

miniXyce

miniXyce uses a modified nodal analysis to approximate the governing differential algebraic equations (DAEs) for electronic devices and their coupling via a network and Kirchoff's conservation laws. The resulting nonlinear problem is solved iteratively using a fully coupled Newton method,which in turn results in a linear system that is solved using an iterative method.

Source record↗
DOE OSTI2015

An algebraic method for constructing stable and consistent autoregressive filters

In this paper, we introduce an algebraic method to construct stable and consistent univariate autoregressive (AR) models of low order for filtering and predicting nonlinear turbulent signals with memory depth. By stable, we refer to the classical stability condition for the AR model. By consistent, we refer to the classical consistency constraints of Adams–Bashforth methods of order-two. One attractive feature of this algebraic method is that the model parameters can be obtained without directly knowing any training data set as opposed to many standard, regression-based parameterization methods. It takes only long-time average statistics as inputs. The proposed method provides a discretization time step interval which guarantees the existence of stable and consistent AR model and simultaneously produces the parameters for the AR models. In our numerical examples with two chaotic time series with different characteristics of decaying time scales, we find that the proposed AR models produce significantly more accurate short-term predictive skill and comparable filtering skill relative to the linear regression-based AR models. These encouraging results are robust across wide ranges of discretization times, observation times, and observation noise variances. Finally, we also find that the proposed model produces an improved short-time prediction relative to the linear regression-based AR-models in forecasting a data set that characterizes the variability of the Madden–Julian Oscillation, a dominant tropical atmospheric wave pattern.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Exchange current density model for the contact-determined current-voltage behavior of solar cells

An analytic expression for the current–voltage [J(V)] behavior of a solar cell as limited by equilibrium exchange current densities of both carriers at both contacts is derived. The partial currents at both contacts to a generic semiconductor absorber are assumed to be linearly proportional to the excess carrier concentration at the interface with the contacts (e.g., as with Schottky-like contacts). The assumption that the quasi-Fermi levels in the absorber are approximately flat leads to an algebraic solution for the applied voltage as a function of current, which is inverted to obtain the analytic J(V) curve. The J(V) curve reveals distinct behavior associated with electrons and holes, separately, and allows for the determination of all critical performance parameters. In particular, it demonstrates how the characteristic features of the J(V) curve depend on the relative rate at which a particular carrier (electron or hole) is collected at one contact vs the other, rather than the relative rate of electron vs hole collection at a single contact. Furthermore, the model provides a unified explanation of how majority carrier extraction limitations cause nonideal J(V) behaviors such as S-shaped curves and dark/light crossover (i.e., failure of superposition). The efficacy and limitations of the model when applied to Schottky-type and doped semiconductor contacts are discussed. The work serves as a theoretical guide to scientists studying solar cells that are thought to be primarily limited by their contacts.

Solar cells↗
DOE OSTI2025

A Novel Approach for Computing Rigid Body Motion Using Linear Accelerations

Here, a novel approach is presented for computing general rigid body motion based on a few known linear accelerations. This method utilizes linear acceleration data obtained from three distinct points on the body, all within a body-fixed reference frame. The only requirement is that the three chosen points must not be collinear. A system of differential-algebraic equations is derived, combining principles of rigid body kinematics with theory of the rotation group SO(3). These equations provide a framework for numerically computing various motion parameters, including angular velocity, angular acceleration, body orientation, velocity field, acceleration field, and displacement field. By numerically solving this system of equations, we can fully characterize rigid body motion in three-dimensional space. A numerical example is provided to demonstrate the practical implementation and efficacy of the proposed technique, illustrating its potential for accurate motion computation in various applications.

42 ENGINEERING↗
NASA NTRS1975

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, it is assumed that the total temperature is constant. Therefore, the energy equation can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed. It is shown that the system is strictly hyperbolic and well posed for the initial value problems. Boundary conditions are described such that the linearized system is well posed. The Hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗
NASA NTRS1989

Convergence of infinite dimensional sampled LQR problems - Theory and numerical results

A theory is developed for the convergence of the closed-loop solution to infinite-dimensional discrete-time linear-quadratic regulator (LQR) problems on the infinite time interval to the solution of a corresponding continuous-time LQR problem as the length of the sampling interval tends toward zero. Convergence of solutions to the operator algebraic Riccati equation and corresponding optimal feedback control gains is guaranteed under appropriate uniform stabilizability and detectability conditions and consistent sampling. Also presented are numerical results involving the optimal LQ control of a heat or diffusion equation, a hereditary or delay differential equation, and a hybrid system of ordinary and partial differential equations describing the transverse vibration of a cantilevered Voigt-Kelvin viscoelastic beam with tip mass.

Rosen, I. G.↗
NASA NTRS1992

Computation of the flow field in a centrifugal impeller with splitter blades

To support the design effort of the Space Transportation Main Engine (STME) Fuel Pump Stage, viscous flow calculations were performed in a centrifugal impeller with splitter blades. These calculations were carried out with a Navier-Stokes solver (MINT), which employs a linearized block-implicit Alternating Direction Implicit (ADI) procedure to iteratively solve a finite difference form of the system of conservation equations of mass, momentum, and energy in body-fitted coordinates. A computational grid was generated algebraically for the 'channel' between two main blades of the impeller and extended both upstream of the impeller inlet and downstream of the impeller exit so that the appropriate boundary conditions could be applied. The results of the calculations show that although the overall level of flow distortion near the impeller exit is not very large, there is a noticeable difference between the flow patterns in the two 'passages' (one passage between the pressure side of the splitter blade and the suction side of the next full blade).

Dejong, Frederik J.↗
NASA NTRS1979

Models of failure in linear systems

The concept of reliability in systems theory is discussed. Emphasis is placed on identifying failed or degraded systems and constructing a strategy to control the system in the failed mode. Systems are assumed to be linear and transformations into the failed mode occur smoothly. The family into which a nominal system is embedded is examined along with the relations within the family. The class of systems into which a given system can be algebraically transformed is investigated.

Martin, C. F.↗
NASA NTRS1981

Finite-analytic numerical solution of heat transfer in two-dimensional cavity flow

Heat transfer in cavity flow is numerically analyzed by a new numerical method called the finite-analytic method. The basic idea of the finite-analytic method is the incorporation of local analytic solutions in the numerical solutions of linear or nonlinear partial differential equations. In the present investigation, the local analytic solutions for temperature, stream function, and vorticity distributions are derived. When the local analytic solution is evaluated at a given nodal point, it gives an algebraic relationship between a nodal value in a subregion and its neighboring nodal points. A system of algebraic equations is solved to provide the numerical solution of the problem. The finite-analytic method is used to solve heat transfer in the cavity flow at high Reynolds number (1000) for Prandtl numbers of 0.1, 1, and 10.

Chen, C.-J.↗
DOE OSTI2015

Application of Adomian Decomposition for Multi-Machine Power System Simulation

In this paper, the merits of Adomian Decomposition Method (ADM) have been investigated for the time domain simulation of multi-machine power systems. ADM is an analytical approximation method for the solution of linear, nonlinear, deterministic and stochastic operator equations. It approximates the nonlinearities in the equations using Adomian polynomial series and that can form a rapidly convergent sequence of analytic functions. We apply the ADM approach for simulation of differential algebraic equations representing the detailed multimachine power system. The ADM approach has been applied on the widely used IEEE 3 generator 9 bus system and IEEE 10 generator 39 bus system for 3 fault simulations. In our simulations we found that the ADM approach is faster than the standard trapezoidal time integration method for the comparable accuracy.

Gurrala, Gurunath↗
DOE OSTI2016

Efficient Adjoint Computation of Hybrid Systems of Differential Algebraic Equations with Applications in Power Systems

Sensitivity analysis is an important tool to describe power system dynamic behavior in response to parameter variations. It is a central component in preventive and corrective control applications. The existing approaches for sensitivity calculations, namely, finite-difference and forward sensitivity analysis, require a computational effort that increases linearly with the number of sensitivity parameters. In this work, we investigate, implement, and test a discrete adjoint sensitivity approach whose computational effort is effectively independent of the number of sensitivity parameters. The proposed approach is highly efficient for calculating trajectory sensitivities of larger systems and is consistent, within machine precision, with the function whose sensitivity we are seeking. This is an essential feature for use in optimization applications. Moreover, our approach includes a consistent treatment of systems with switching, such as DC exciters, by deriving and implementing the adjoint jump conditions that arise from state and time-dependent discontinuities. The accuracy and the computational efficiency of the proposed approach are demonstrated in comparison with the forward sensitivity analysis approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗
DOE OSTI2011

Quantum fields in toroidal topology

The standard representation of c*-algebra is used to describe fields in compactified space-time dimensions characterized by topologies of the type {Gamma}{sub D}{sup d}=(S{sup 1}){sup d}xM{sup D-d}. The modular operator is generalized to introduce representations of isometry groups. The Poincare symmetry is analyzed and then we construct the modular representation by using linear transformations in the field modes, similar to the Bogoliubov transformation. This provides a mechanism for compactification of the Minkowski space-time, which follows as a generalization of the Fourier integral representation of the propagator at finite temperature. An important result is that the 2x2 representation of the real-time formalism is not needed. The end result on calculating observables is described as a condensate in the ground state. We initially analyze the free Klein-Gordon and Dirac fields, and then formulate non-abelian gauge theories in {Gamma}{sub D}{sup d}. Using the S-matrix, the decay of particles is calculated in order to show the effect of the compactification. - Highlights: > C*-algebra is used to describe fields in compactified space-time dimensions. > The space-time is characterized by toroidal topologies. > Representations of the Poincare group are studied by using the modular operator. > We derive non-abelian gauge theories in compactified regions of space-time. > We show the compactification effect in the decay of particles using the S-matrix.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Conformal geometrodynamics regained: Gravity from duality

There exist several ways of constructing general relativity from ‘first principles’: Einstein’s original derivation, Lovelock’s results concerning the exceptional nature of the Einstein tensor from a mathematical perspective, and Hojman–Kuchař-Teitelboim’s derivation of the Hamiltonian form of the theory from the symmetries of space–time, to name a few. Here I propose a different set of first principles to obtain general relativity in the canonical Hamiltonian framework without presupposing space–time in any way. I first require consistent propagation of scalar spatially covariant constraints (in the Dirac picture of constrained systems). I find that up to a certain order in derivatives (four spatial and two temporal), there are large families of such consistently propagated constraints. Then I look for pairs of such constraints that can gauge-fix each other and form a theory with two dynamical degrees of freedom per space point. This demand singles out the ADM Hamiltonian either in (i) CMC gauge, with arbitrary (finite, non-zero) speed of light, and an extra term linear in York time, or (ii) a gauge where the Hubble parameter is conformally harmonic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Feedback Control Approaches for Restoration of Power Grids from Blackouts

The automated restoration of power systems with variable energy resources is a timely problem to tackle. Automated restoration advice can support operators in deciding on strategic actions to restore power grids from a blackout with a mix of conventional and renewable generation resources. To this end, this paper frames the restoration process of power grids with solar resources as a nonlinear dynamic model with algebraic constraints in discrete time which is steered by feedback control loops. We discuss two feedback-control strategies based on greedy and reinforcement learning algorithms, and contrast their performance with restoration plans generated by a mixed-integer linear program. We found that the reinforcement learning algorithm infers restoration actions faster than the greedy one. However, the tuning process of the reinforcement learning parameters is slower than for the greedy one.

machine learning↗
DOE OSTI2014

On linear groups of degree 2 over a finite commutative ring

Let p > 3 be a prime number and F{sub p} be a field of p elements. Let K be a commutative and associative ring obtained by adjoining to F{sub p} an element α such that α satisfies a polynomial over F{sub p} and a polynomial of the least degree whose root is α can be written as a product of distinct polynomials irreducible over F{sub p}. We prove that the special linear group SL{sub 2}(K) is generated by two elementary transvections ( (table) ), ( (table) )

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗