Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear system”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗
DOE OSTI2019

Quantum algorithm for simulating the wave equation

We present a quantum algorithm for simulating the wave equation under Dirichlet and Neumann boundary conditions. The algorithm uses Hamiltonian simulation and quantum linear system algorithms as subroutines. It relies on factorizations of discretized Laplacian operators to allow for polynomially improved scaling in truncation errors and improved scaling for state preparation relative to general purpose quantum algorithms for solving linear differential equations. Relative to classical algorithms for simulating the D-dimensional wave equation, our quantum algorithm achieves exponential space savings and achieves a speedup which is polynomial for fixed D and exponential in D. Furthermore, we also consider using Hamiltonian simulation for Klein-Gordon equations and Maxwell's equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Centralized, decentralized, and independent control of a flexible manipulator on a flexible base

The dynamics and control of a flexible manipulator arm with payload mass on a flexible base in space are considered. The controllers are provided by one torquer at the center of the base and one torquer at the connection joint of the robot and the base. The nonlinear dynamics of the system is modeled by applying the finite element method and Lagrangian formula. Three control strategies are considered and compared, i.e., centralized control, decentralized control, and independent control. All these control designs are based on the linear quadratic regulator theory. A mathematical decomposition is used in the decentralization process so that the coupling between the subsystems is weak, while a physical decomposition is used in the independent control design process. For both the decentralized and the independent controls, the stability of the overall linear system is checked before a numerical simulations is initiated. Two numerical examples show that the response of the independent control system are close to those of the centralized control system, while the responses of the decentralized control system are not.

Li, Feiyue↗
DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗
NASA NTRS1976

Regions of constrained maximum likelihood parameter identifiability

This short paper considers the parameter-identification problem of general discrete-time, nonlinear, multiple input-multiple output dynamic systems with Gaussian white distributed measurement errors. Knowledge of the system parameterization is assumed to be available. Regions of constrained maximum likelihood (CML) parameter identifiability are established. A computation procedure employing interval arithmetic is proposed for finding explicit regions of parameter identifiability for the case of linear systems.

Lee, C.-H.↗
NASA NTRS1975

Regions of constrained maximum likelihood parameter identifiability

This paper considers the parameter identification problem of general discrete-time, nonlinear, multiple-input/multiple-output dynamic systems with Gaussian-white distributed measurement errors. Knowledge of the system parameterization is assumed to be known. Regions of constrained maximum likelihood (CML) parameter identifiability are established. A computation procedure employing interval arithmetic is proposed for finding explicit regions of parameter identifiability for the case of linear systems. It is shown that if the vector of true parameters is locally CML identifiable, then with probability one, the vector of true parameters is a unique maximal point of the maximum likelihood function in the region of parameter identifiability and the CML estimation sequence will converge to the true parameters.

Lee, C.-H.↗
NASA NTRS1984

Refinement and evaluation of helicopter real-time self-adaptive active vibration controller algorithms

A Real-Time Self-Adaptive (RTSA) active vibration controller was used as the framework in developing a computer program for a generic controller that can be used to alleviate helicopter vibration. Based upon on-line identification of system parameters, the generic controller minimizes vibration in the fuselage by closed-loop implementation of higher harmonic control in the main rotor system. The new generic controller incorporates a set of improved algorithms that gives the capability to readily define many different configurations by selecting one of three different controller types (deterministic, cautious, and dual), one of two linear system models (local and global), and one or more of several methods of applying limits on control inputs (external and/or internal limits on higher harmonic pitch amplitude and rate). A helicopter rotor simulation analysis was used to evaluate the algorithms associated with the alternative controller types as applied to the four-bladed H-34 rotor mounted on the NASA Ames Rotor Test Apparatus (RTA) which represents the fuselage. After proper tuning all three controllers provide more effective vibration reduction and converge more quickly and smoothly with smaller control inputs than the initial RTSA controller (deterministic with external pitch-rate limiting). It is demonstrated that internal limiting of the control inputs a significantly improves the overall performance of the deterministic controller.

Davis, M. W.↗
NASA NTRS1985

Nested subcritical flows within supercritical systems

In supercritical systems the design inlet and outlet pressures are maintained above the thermaodynamic critical pressure P sub C. Designers rely on this simple rule of thumb to circumvent problems associated with a subcritical pressure regime nested within the supercritical pressure system along with the uncertainties in heat transfer, fluid mechanics, and thermophysical property variations. The simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines and linear systems, where nested two-phase regions can exist. Examples for a free-jet expansion with backpressure greater than P sub C and a rotor (bearing) with ambient pressure greater than P sub C illustrate the existence of subcritical pressure regimes nested within supercritical systems.

Hendricks, R. C.↗
NASA NTRS1975

System balance analysis for vector computers

The availability of vector processors capable of sustaining computing rates of 10 to the 8th power arithmetic results pers second raised the question of whether peripheral storage devices representing current technology can keep such processors supplied with data. By examining the solution of a large banded linear system on these computers, it was found that even under ideal conditions, the processors will frequently be waiting for problem data.

Knight, J. C.↗
NASA NTRS2010

Low-Friction, Low-Profile, High-Moment Two-Axis Joint

The two-axis joint is a mechanical device that provides two-degrees-of-freedom motion between connected components. A compact, moment-resistant, two-axis joint is used to connect an electromechanical actuator to its driven structural members. Due to the requirements of the overall mechanism, the joint has a low profile to fit within the allowable space, low friction, and high moment-reacting capability. The mechanical arrangement of this joint can withstand high moments when loads are applied. These features allow the joint to be used in tight spaces where a high load capability is required, as well as in applications where penetrating the mounting surface is not an option or where surface mounting is required. The joint consists of one base, one clevis, one cap, two needle bearings, and a circular shim. The base of the joint is the housing (the base and the cap together), and is connected to the grounding structure via fasteners and a bolt pattern. Captive within the housing, between the base and the cap, are the rotating clevis and the needle bearings. The clevis is attached to the mechanical system (linear actuator) via a pin. This pin, and the rotational movement of the clevis with respect to the housing, provides two rotational degrees of freedom. The larger diameter flange of the clevis is sandwiched between a pair of needle bearings, one on each side of the flange. During the assembly of the two-axis joint, the circular shims are used to adjust the amount of preload that is applied to the needle bearings. The above arrangement enables the joint to handle high moments with minimal friction. To achieve the high-moment capability within a low-profile joint, the use of depth of engagement (like that of a conventional rotating shaft) to react moment is replaced with planar engagement parallel to the mounting surface. The needle bearings with the clevis flange provide the surface area to react the clevis loads/moments into the joint housing while providing minimal friction during rotation. The diameter of the flange and the bearings can be increased to react higher loads and still maintain a compact surface mounting capability. This type of joint can be used in a wide variety of mechanisms and mechanical systems. It is especially effective where precise, smooth, continuous motion is required. For example, the joint can be used at the end of a linear actuator that is required to extend and rotate simultaneously. The current design application is for use in a spacecraft docking-system capture mechanism. Other applications might include industrial robotic or assembly line apparatuses, positioning systems, or in the motion-based simulator industry that employs complex, multi-axis manipulators for various types of motions.

Lewis, James L.↗
DOE OSTI2024

Parallel interior-point solver for block-structured nonlinear programs on SIMD/GPU architectures

Here, we investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses SIMD/GPU accelerators locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations. We demonstrate our method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the reduction method is 50x faster than the sparse linear solver HSL MA57 running in serial on the CPU, and 6x faster than Pardiso running in parallel on CPU on the same number of processes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Two-phase flows within systems with ambient pressure

In systems where the design inlet and outlet pressures are maintained above the thermodynamic critical pressure, it is often assumed that two phase flows within the system cannot occur. Designers rely on this simple rule of thumb to circumvent problems associated with a highly compressible two phase flow occurring within the supercritical pressure system along with the uncertainties in rotordynamics, load capacity, heat transfer, fluid mechanics, and thermophysical property variations. The simple rule of thumb is adequate in many low power designs but is inadequate for high performance turbomachines and linear systems, where two phase regions can exist even though outlet pressure is greater than critical pressure. Rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two phase zone can differ significantly from those for a single-phase zone. Using the Reynolds equation the angular velocity, eccentricity, geometry, and ambient conditions are varied to determine the point of two phase flow incipience.

Hendricks, R. C.↗
NASA NTRS1989

Structural modeling for control design (articulated multibody component representation)

High gain, high frequency flexible responses in gimbaled multibody systems are discussed. Their origin and physical significance are described in terms of detailed mass and stiffness modeling at actuator/sensor interfaces. Guyan Reduction, Generalized Dynamic Reduction, inadequate mass modeling detail, as well as system mode truncation, are shown to suppress the high gain high frequency response and thereby lose system flexibility important for stability and performance predictions. Model validation by modal survey testing is shown to risk similar loss of accuracy. Difficulties caused by high frequency responses in component mode simulations, such as DISCOS, and also linearized system mode simulations, are described, and approaches for handling these difficulties are discussed.

Haugse, E. D.↗
DOE OSTI2023

ETG turbulence in a tokamak pedestal

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking at positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗
NASA NTRS1976

Information and display requirements for independent landing monitors

The ways an Independent Landing Monitor (ILM) may be used to complement the automatic landing function were studied. In particular, a systematic procedure was devised to establish the information and display requirements of an ILM during the landing phase of the flight. Functionally, the ILM system is designed to aid the crew in assessing whether the total system (e.g., avionics, aircraft, ground navigation aids, external disturbances) performance is acceptable, and, in case of anomaly, to provide adequate information to the crew to select the least unsafe of the available alternatives. Economically, this concept raises the possibility of reducing the primary autoland system redundancy and associated equipment and maintenance costs. The required level of safety for the overall system would in these cases be maintained by upgrading the backup manual system capability via the ILM. A safety budget analysis was used to establish the reliability requirements for the ILM. These requirements were used as constraints in devising the fault detection scheme. Covariance propagation methods were used with a linearized system model to establish the time required to correct manually perturbed states due to the fault. Time-to-detect and time-to-correct requirements were combined to devise appropriate altitudes and strategies for fault recovery.

Karmarker, J. S.↗
NASA NTRS1979

Decoupling control of a long flexible beam in orbit

The paper presents a method of control for large flexible systems using state variable feedback, with a long flexible beam given as an example. These feedback gains are selected: (1) based on the decoupling of the original coordinates and to obtain proper damping and (2) by applying the linear regulator problem to the individual modal coordinates separately. It is shown that the linear control law thus obtained are then evaluated by numerical integration of the non-linear system equations. Also included are results showing the effects (control spillover) on the uncontrolled modes when the number of controllers is less than the number of modes, and the effects of inaccurate knowledge of the control influence coefficients which lead to errors in the calculated feedback gains.

Reddy, A. S. S. R.↗
NASA NTRS1981

Accelerated convergence of structured banded systems using constrained corrections

An efficient iterative method for solving a structured banded system of equations is described. The method was developed for a full potential flow program and uses a basic interation step, a dynamic relation step, and a multigrid concept of constraining iterative corrections. The solution of a large linear system of equations is examined. Efficient iterative methods have become attractive for large problems. In the nonlinear cases, these iterations may be effectively merged to improve convergence rates.

Kneile, K.↗