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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1985

On the factorization of block-tridiagonals without storage constraints

In many programs solving difference equations, problem size is restricted by the number of available memory cells. A strategy has been developed to permit trade-offs between the number of floating point operations required and storage requirements for the solution of certain problems such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement will roughly double the number of computations. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5 sq m + 1 temporary storage cells. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.

Merriam, M. L.↗

DOE OSTI2018

Technical Note: An efficient method for accelerating the spin-up process for process-based biogeochemistry models

To better understand the role of terrestrial ecosystems in the global carbon cycle and their feedbacks to the global climate system,process-based biogeochemistry models need to be improved with respect tomodel parameterization and model structure. To achieve these improvements,the spin-up time for those differential equation-based models needs to be shortened. Here, an algorithm for a fast spin-up was developed by finding the exact solution of a linearized system representing the cyclo-stationary state of a model and implemented in a biogeochemistry model, the Terrestrial Ecosystem Model (TEM). With the new spin-up algorithm, we showed that the model reacheda steady state in less than 10 years of computing time, while the original method requires more than 200 years on average of model run. For the test sites with five different plant functional types, the new method saves over 90 % of the original spin-up time in site-level simulations. In North American simulations, average spin-up time savings for all grid cells is 85 % for either the daily or monthly version of TEM. Thus, the developed spin-up method shall be used for future quantification of carbon dynamics at fine spatial and temporal scales.

58 GEOSCIENCES↗

NASA NTRS1974

A Computing Method for Sound Propagation Through a Nonuniform Jet Stream

Understanding the principles of jet noise propagation is an essential ingredient of systematic noise reduction research. High speed computer methods offer a unique potential for dealing with complex real life physical systems whereas analytical solutions are restricted to sophisticated idealized models. The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions and a more suitable approach was needed. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Padula, S. L.↗

NASA NTRS1974

A computing method for sound propagation through a nonuniform jet stream

The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Padula, S. L.↗

DOE OSTI2018

Diagonal and circulant or skew-circulant splitting preconditioners for spatial fractional diffusion equations

We propose three new preconditioners: diagonal and optimal-circulant splitting preconditioner, diagonal and skew-circulant splitting preconditioner, and diagonal and optimal-skew-circulant splitting preconditioner for solving the diagonal-plus-Toeplitz linear system discretized from the spatial fractional diffusion equations. Theoretical analysis shows that these three preconditioners can make the eigenvalues of the preconditioned matrices be clustered around 1, especially when the grids of the discretizations are refined. These results coincide with the one about the diagonal and circulant splitting preconditioner constructed recently by Bai et al. (Numer Linear Algebra Appl 24:e2093, 2017). Numerical experiments exhibit that the proposed preconditioners can significantly improve the convergence of the Krylov subspace iteration methods like GMRES and BiCGSTAB, and they outperform the diagonal and circulant splitting preconditioner as well.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Frequency analysis via the method of moment functionals

Several variants are presented of a linear-in-parameters least squares formulation for determining the transfer function of a stable linear system at specified frequencies given a finite set of Fourier series coefficients calculated from transient nonstationary input-output data. The basis of the technique is Shinbrot's classical method of moment functionals using complex Fourier based modulating functions to convert a differential equation model on a finite time interval into an algebraic equation which depends linearly on frequency-related parameters.

Pearson, A. E.↗

NASA NTRS1991

Frequency analysis via the method of moment functionals

Pearson, A. E.↗

NASA NTRS1972

The optimal control of merging aircraft - Implementation of the hybrid air traffic controller.

The control of merging aircraft is formulated as a finite-time, quadratic optimal control problem of a linear system with state and control constraints. The purpose of this paper is to demonstrate that the Hybrid Air Traffic Controller (HAC), which has been previously developed as a solution to this problem, may be easily implemented. Use is made of both the properties of the algebraic solution to the matrix Riccati equation and the structure of the linear model. This approach results in a real-time synthesis procedure for the HAC which does not rely on iterative numerical integration techniques.

Schatz, J. G.↗

NASA NTRS1992

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗

NASA NTRS2001

State-Dependent Pseudo-Linear Filter for Spacecraft Attitude and Rate Estimation

This paper presents the development and performance of a special algorithm for estimating the attitude and angular rate of a spacecraft. The algorithm is a pseudo-linear Kalman filter, which is an ordinary linear Kalman filter that operates on a linear model whose matrices are current state estimate dependent. The nonlinear rotational dynamics equation of the spacecraft is presented in the state space as a state-dependent linear system. Two types of measurements are considered. One type is a measurement of the quaternion of rotation, which is obtained from a newly introduced star tracker based apparatus. The other type of measurement is that of vectors, which permits the use of a variety of vector measuring sensors like sun sensors and magnetometers. While quaternion measurements are related linearly to the state vector, vector measurements constitute a nonlinear function of the state vector. Therefore, in this paper, a state-dependent linear measurement equation is developed for the vector measurement case. The state-dependent pseudo linear filter is applied to simulated spacecraft rotations and adequate estimates of the spacecraft attitude and rate are obtained for the case of quaternion measurements as well as of vector measurements.

Bar-Itzhack, Itzhack Y.↗

NASA NTRS1986

Stability with large step sizes for multistep discretizations of stiff ordinary differential equations

One-leg and multistep discretizations of variable-coefficient linear systems of ODEs having both slow and fast time scales are investigated analytically. The stability properties of these discretizations are obtained independent of ODE stiffness and compared. The results of numerical computations are presented in tables, and it is shown that for large step sizes the stability of one-leg methods is better than that of the corresponding linear multistep methods.

Majda, George↗

DOE OSTI2016

A method to the impact assessment of the returning grazing land to grassland project on regional eco-environmental vulnerability

The Chinese government has conducted the Returning Grazing Land to Grassland Project (RGLGP) across large portions of grasslands from western China since 2003. In order to explore and understand the impact in the grassland's eco-environment during the RGLGP, we utilized Projection Pursuit Model (PPM) and Geographic Information System (GIS) to develop a spatial assessment model to examine the ecological vulnerability of the grassland. Our results include five indications: (1) it is practical to apply the spatial PPM on ecological vulnerability assessment for the grassland. This methodology avoids creating an artificial hypothesis, thereby providing objective results that successfully execute a multi-index assessment process and analysis under non-linear systems in eco-environments; (2) the spatial PPM is not only capable of evaluating regional eco-environmental vulnerability in a quantitative way, but also can quantitatively demonstrate the degree of effect in each evaluation index for regional eco-environmental vulnerability; (3) the eco-environment of the Xianshui River Basin falls into the medium range level. The normalized difference vegetation index (NDVI) and land use cover and change (LUCC) crucially influence the Xianshui River Basin's eco-environmental vulnerability. Generally, in the Xianshui River Basin, regional eco-environmental conditions improved during 2000 and 2010. The RGLGP positively affected NDVI and LUCC structure, thereby promoting the enhancement of the regional eco-environment; (4) the Xianshui River Basin divides its ecological vulnerability across different levels; therefore our study investigates three ecological regions and proposes specific suggestions for each in order to assist in eco-environmental protection and rehabilitation; and lastly that (5) the spatial PPM established by this study has the potential to be applied on all types of grassland eco-environmental vulnerability assessments under the RGLGP and under the similar conditions in the Returning Agriculture Land to Forest Project (RALFP). However, when establishing an eco-environmental vulnerability assessment model, it is necessary to choose suitable evaluation indexes in accordance with regional eco-environmental characteristics. - Highlights: • We present a method for regional eco-environmental vulnerability assessment. • The method combines Projection Pursuit Model with Geographic Information System. • The Returning Grazing Land to Grassland Project is crucial to environment recovery. • The method is more objective to assess regional eco-environmental vulnerability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Preliminary Development of a Drift-Flux Model Based Pin-Resolved Core Thermal-Hydraulic Code for Dedicated Use in Efficient High-Fidelity Neutronics Coupled Calculations

As the high-fidelity multiphysics coupled core calculation becomes more practical with the advent of high-performance computing technology, efficient coupling of a neutronics code with a thermal-hydraulic (T/H) code at the pin level poses a new challenge. The COBRA-TF (CTF) code which is one of the well-known subchannel analysis codes based on the two-fluid, three-field model was coupled with many neutron transport codes such as MPACT, TORTTD, and nTRACER due to the easiness of coupling by its well-established coupling interface. Nevertheless, the calculation efficiency of such coupled codes is not satisfactory because of the heavy computational burden of CTF. The MATRA subchannel analysis code of KAERI which is based on the homogeneous equilibrium model (HEM) was coupled with neutronics codes as well. The choice of the HEM and the marching scheme in MATRA reduces the computational burden significantly. However, the applicability of the HEM for the two-phase flow conditions is limited to normal operating conditions of pressurized water reactors (PWRs) and therefore a better model which would allow wider applicability than HEM is necessary in order to deal with the core states involving nontrivial voiding. The drift-flux model (DFM) which uses 5 equations for two-fluid analyses by adopting the mixture velocity instead of two phasic velocities can be a good compromise between the two-fluid model and the HEM. With the DFM, the computational cost can be reduced and moreover the modeling of the interfacial area constitutive relations can be avoided. The phase separation for the mixture velocity can be incorporated by using the drift-flux parameters, C0 and Vgj, which can be obtained by using empirical formulas. Several subchannel codes such as WOSUB and ASSERT were developed based on the DFM. It was also adopted in some system analysis codes such as RETRAN3D. In the work here, the development of a DFM based pin-resolved T/H analysis code, ESCOT (Efficient Simulator of Core Thermal-hydraulics), is introduced for the use in efficient coupled calculations with a high-fidelity neutronics code. The field equations and the numerical methods are briefly described. Then the efficiency tests of various linear system solvers including the Krylov subspace methods (KSMs) are presented. A test problem is solved by ESCOT, MATRA and CTF for performance comparison and validation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Simulations of plasmas and fluids using anti-symmetric models

ALMA (Anti-symmetric, Large-Moment, Accelerated) is a fast, flexible, and scalable toolkit designed to solve hyperbolic conservation law systems in hybrid supercomputers. Here this manuscript describes the theoretical background and implementation of ALMA, which uses the anti-symmetric formulation of fluids to obtain simple, robust, and easily paralellizable code. Practical GPU acceleration is realized on entire applications with an overall gain factor of 2 to 4. ALMA also provides a parallel, GPU accelerated sparse solver based on geometric multigrid, capable of diagonalizing linear systems with 239 unknowns. Here we demonstrate ALMA's scaling and performance in petascale supercomputers and use standard fluid models to verify the overall approach with canonical benchmark problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2012

An efficient implementation of the Chebyshev Rational Approximation Method (CRAM) for solving the burnup equations

The Chebyshev Rational Approximation Method (CRAM) has been recently introduced by the authors for solving the burnup equations with excellent results. This method has been shown to be capable of simultaneously solving an entire burnup system with thousands of nuclides both accurately and efficiently. The method was prompted by an analysis of the spectral properties of burnup matrices and it can be characterized as the best rational approximation on the negative real axis. The coefficients of the rational approximation are fixed and have been reported for various approximation orders. In addition to these coefficients, implementing the method only requires a linear solver. This paper describes an efficient method for solving the linear systems associated with the CRAM approximation. The introduced direct method is based on sparse Gaussian elimination where the sparsity pattern of the resulting upper triangular matrix is determined before the numerical elimination phase. The stability of the proposed Gaussian elimination method is discussed based on considering the numerical properties of burnup matrices. Suitable algorithms are presented for computing the symbolic factorization and numerical elimination in order to facilitate the implementation of CRAM and its adoption into routine use. The accuracy and efficiency of the described technique are demonstrated by computing the CRAM approximations for a large test case with over 1600 nuclides. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1985

Characteristics pertaining to a stiffness cross-coupled Jeffcott model

Rotordynamic studies of complex systems utilizing multiple degree-of-freedom analysis have been performed to understand response, loads, and stability. In order to understand the fundamental nature of rotordynamic response, the Jeffcott rotor model has received wide attention. The purpose of this paper is to provide a generic rotordynamic analysis of a stiffness cross-coupled Jeffcott rotor model to illustrate characteristics of a second order stiffness-coupled linear system. The particular characteristics investigated were forced response, force vector diagrams, response orbits, and stability. Numerical results were achieved through a fourth order Runge-Kutta method for solving differential equations and the Routh Hurwitz stability criterion. The numerical results were verified to an exact mathematical solution for the steady state response.

Spanyer, K. L.↗

DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1982

On the inversion of block tridiagonals without storage constraints

A strategy was developed to permit trade-offs between the number of floating point operations required and the storage requirements for the solution of certain difference problems, such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement roughly doubles the number of computations. Reducing the storage more than this tends to make the number of computations prohibitively large. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5m(2) + 1 temporary storage cells. In many cases m is a constant and quite small. For example, in solving a factored form of the three-dimensional Navier-Stokes equations, the size m of the block tridiagonals is 5. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.

Merriam, M. L.↗