Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equation systems”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Bounded Linear Stability Analysis - A Time Delay Margin Estimation Approach for Adaptive Control

This paper presents a method for estimating time delay margin for model-reference adaptive control of systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent the conventional model-reference adaptive law by a locally bounded linear approximation within a small time window using the comparison lemma. The locally bounded linear approximation of the combined adaptive system is cast in a form of an input-time-delay differential equation over a small time window. The time delay margin of this system represents a local stability measure and is computed analytically by a matrix measure method, which provides a simple analytical technique for estimating an upper bound of time delay margin. Based on simulation results for a scalar model-reference adaptive control system, both the bounded linear stability method and the matrix measure method are seen to provide a reasonably accurate and yet not too conservative time delay margin estimation.

Nguyen, Nhan T.↗
NASA NTRS1991

Physical lumping methods for developing linear reduced models for high speed propulsion systems

In gasdynamic systems, information travels in one direction for supersonic flow and in both directions for subsonic flow. A shock occurs at the transition from supersonic to subsonic flow. Thus, to simulate these systems, any simulation method implemented for the quasi-one-dimensional Euler equations must have the ability to capture the shock. In this paper, a technique combining both backward and central differencing is presented. The equations are subsequently linearized about an operating point and formulated into a linear state space model. After proper implementation of the boundary conditions, the model order is reduced from 123 to less than 10 using the Schur method of balancing. Simulations comparing frequency and step response of the reduced order model and the original system models are presented.

Immel, S. M.↗
DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Accuracy and convergence of a finite element algorithm for turbulent boundary layer flow

The Galerkin-Weighted Residuals formulation is employed to derive an implicit finite element solution algorithm for the nonlinear parabolic partial differential equation system governing turbulent boundary layer flow. Solution accuracy and convergence with discretization refinement are quantized in several error norms using linear and quadratic basis functions. Richardson extrapolation is used to isolate integration truncation error in all norms, and Newton iteration is employed for all equation solutions performed in double-precision. The mathematical theory supporting accuracy and convergence concepts for linear elliptic equations appears extensible to the nonlinear equations characteristic of turbulent boundary layer flow.

Soliman, M. O.↗
DOE OSTI2020

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

17 WIND ENERGY↗
NASA NTRS1993

An efficient model for coupling structural vibrations with acoustic radiation

The scattering of an incident wave by a flexible panel is studied. The panel vibration is governed by the nonlinear plate equations while the loading on the panel, which is the pressure difference across the panel, depends on the reflected and transmitted waves. Two models are used to calculate this structural-acoustic interaction problem. One solves the three dimensional nonlinear Euler equations for the flow-field coupled with the plate equations (the fully coupled model). The second uses the linear wave equation for the acoustic field and expresses the load as a double integral involving the panel oscillation (the decoupled model). The panel oscillation governed by a system of integro-differential equations is solved numerically and the acoustic field is then defined by an explicit formula. Numerical results are obtained using the two models for linear and nonlinear panel vibrations. The predictions given by these two models are in good agreement but the computational time needed for the 'fully coupled model' is 60 times longer than that for 'the decoupled model'.

Frendi, Abdelkader↗
NASA NTRS1990

Simulated dynamic response of a servovalve controlled hydraulic actuator

A general purpose math model of a servovalve controlled hydraulic actuator system is derived. The system consists of a linear actuator with unequal piston areas, a single stage servovalve, a gas charged hydraulic accumulator, and the interconnecting piping. The state equations are integrated using the Advanced Continuous Simulation Language (ACSL) for determining the system's dynamic response characteristics. Using this generalized hydraulic actuator system model, response characteristics were determined for various servovalve commands.

Babcock, Dale A.↗
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↗
DOE OSTI2021

Linear embedding of nonlinear dynamical systems and prospects for efficient quantum algorithms

The simulation of large nonlinear dynamical systems, including systems generated by discretization of hyperbolic partial differential equations, can be computationally demanding. Such systems are important in both fluid and kinetic computational plasma physics. This motivates exploring whether a future error-corrected quantum computer could perform these simulations more efficiently than any classical computer. In this work, we describe a method for mapping any finite nonlinear dynamical system to an infinite linear dynamical system (embedding) and detail three specific cases of this method that correspond to previously studied mappings. Then we explore an approach for approximating the resulting infinite linear system with finite linear systems (truncation). Using a number of qubits only logarithmic in the number of variables of the nonlinear system, a quantum computer could simulate truncated systems to approximate output quantities if the nonlinearity is sufficiently weak. Other aspects of the computational efficiency of the three detailed embedding strategies are also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Pseudospin-triplet pairing in iron-chalcogenide superconductors

Abstract Understanding the pairing symmetry is a crucial theoretical aspect in the study of unconventional superconductivity for interpreting experimental results. Here we study superconductivity of electron systems with both spin and pseudospin-1/2 degrees of freedom. By solving linearized gap equations, we derive a weak coupling criterion for the even-parity spin-singlet pseudospin-triplet pairing. It can generally mix with the on-site s -wave pairing since both of them belong to the same symmetry representation ( A 1 g ) and their mixture could naturally give rise to anisotropic intra-band pairing gap functions with or without nodes. This may directly explain why some of the iron-chalcogenide superconductors are fully gapped (e.g. FeSe thin film) and some have nodes (e.g. LaFePO and LiFeP). We also find that the anisotropy of gap functions can be enhanced when the principal rotation symmetry is spontaneously broken in the normal state such as nematicity, and the energetic stabilization of pseudospin-triplet pairings indicates the coexistence of nematicity and superconductivity. This could be potentially applied to bulk FeSe, where gap anisotropy has been experimentally observed.

Physics↗
DOE OSTI2024

Modeling In-Reactor Temperature Gradients in Lithium Aluminate Ceramic Pellets

In reactor thermal performance of lithium aluminate ceramic pellets are calculated to determine operating temperatures if internal components shift during irradiation. A 2D conduction model was setup to model this scenario with additional physics, including thermal expansion of components to reactor temperatures, creep down of the cladding due to differential pressure between the inside of the rod and reactor coolant system, and net radiation exchange between components. The non-linear coupling of the governing equations presented convergence challenges due to the very different effects from a small temperature change on the radiation vs radial and circumferential conduction equations. Initial temperature guesses for the multidimensional solution were provided from conservative 1D models simulating bounding scenarios. Convergence was greatly improved by taking small steps towards the root solutions and averaging with initial guesses. The final method robustly converged all 132 points at an axial plane after a few dozen iterations. The temperature results show that the additional physics have a secondary effect as the dominant heat transfer mechanism is radial conduction across the gas gap followed by circumferential conduction around components. Circumferential conduction in the lithium aluminate pellets is enhanced by the relatively high thermal conductivity of this ceramic.

Carstens, Nathan A.↗
NASA NTRS1972

Control of linear uncertain systems utilizing mismatched state observers

The control of linear continuous dynamical systems is investigated as a problem of limited state feedback control. The equations which describe the structure of an observer are developed constrained to time-invarient systems. The optimal control problem is formulated, accounting for the uncertainty in the design parameters. Expressions for bounds on closed loop stability are also developed. The results indicate that very little uncertainty may be tolerated before divergence occurs in the recursive computation algorithms, and the derived stability bound yields extremely conservative estimates of regions of allowable parameter variations.

Goldstein, B.↗
NASA NTRS1972

On the stability of numerical integration routines.

Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.

Glover, K.↗
NASA NTRS1973

On the stability of numerical integration routines for ordinary differential equations.

Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.

Glover, K.↗
NASA NTRS1976

The motion and stability of spin stabilized spacecraft with hinged appendages

The dynamics and stability of a spin stabilized spacecraft with a hinged appendage system are treated analytically and numerically. The hinged system consists of a central hub with masses attached to massless booms of fixed length whose orientation relative to the main part can change. The equations of motion for the hinged system, with viscous damping at both hinge points, are linearized about the nominal equilibrium position where the booms are orthogonal to the nominal spin axis. Analytic stability criteria indicate that hinge damping must be present and that for limiting cases, where the spin axis is an axis of symmetry, certain conditions relating the hinge point offset coordinates to the moment of inertia ratio and end masses must be satisfied. The hinge damping is always required for the nominal deployment of hinge members.

Sellappan, R.↗
NASA NTRS1983

Direct and implicit optical matrix-vector algorithms

New direct and implicit algorithms for optical matrix-vector and systolic array processors are considered. Direct rather than indirect algorithms to solve linear systems and implicit rather than explicit solutions to solve second-order partial differential equations are discussed. In many cases, such approaches more properly utilize the advantageous features of optical systolic array processors. The matrix-decomposition operation (rather than solution of the simplified matrix-vector equation that results) is recognized as the computationally burdensome aspect of such problems that should be computed on an optical system. The Householder QR matrix-decomposition algorithm is considered as a specific example of a direct solution. Extensions to eigenvalue computation and formation of matrices of special structure are also noted.

Casasent, D.↗