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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-E-J-02: 4D Digital Tomosynthesis Based On Algebraic Image Reconstruction and Total-Variation Minimization for the Improvement of Image Quality

Purpose: In this paper, we implemented the four-dimensional (4D) digital tomosynthesis (DTS) imaging based on algebraic image reconstruction technique and total-variation minimization method in order to compensate the undersampled projection data and improve the image quality. Methods: The projection data were acquired as supposed the cone-beam computed tomography system in linear accelerator by the Monte Carlo simulation and the in-house 4D digital phantom generation program. We performed 4D DTS based upon simultaneous algebraic reconstruction technique (SART) among the iterative image reconstruction technique and total-variation minimization method (TVMM). To verify the effectiveness of this reconstruction algorithm, we performed systematic simulation studies to investigate the imaging performance. Results: The 4D DTS algorithm based upon the SART and TVMM seems to give better results than that based upon the existing method, or filtered-backprojection. Conclusion: The advanced image reconstruction algorithm for the 4D DTS would be useful to validate each intra-fraction motion during radiation therapy. In addition, it will be possible to give advantage to real-time imaging for the adaptive radiation therapy. This research was supported by Leading Foreign Research Institute Recruitment Program (Grant No.2009-00420) and Basic Atomic Energy Research Institute (BAERI); (Grant No. 2009-0078390) through the National Research Foundation of Korea(NRF) funded by the Ministry of Science, ICT and Future Planning (MSIP)

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A low-rank solver for the stochastic unsteady Navier–Stokes problem

Here we study a low-rank iterative solver for the unsteady Navier–Stokes equations for incompressible flows with a stochastic viscosity. The equations are discretized using the stochastic Galerkin method, and we consider an all-at-once formulation where the algebraic systems at all the time steps are collected and solved simultaneously. The problem is linearized with Picard’s method. To efficiently solve the linear systems at each step, we use low-rank tensor representations within the Krylov subspace method, which leads to significant reductions in storage requirements and computational costs. Combined with effective mean-based preconditioners and the idea of inexact solve, we show that only a small number of linear iterations are needed at each Picard step. The proposed algorithm is tested with a model of flow in a two-dimensional symmetric step domain with different settings to demonstrate the computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Managing uncertainty in data-driven simulation-based optimization

Optimization using data from complex simulations has become an attractive decision-making option, due to ability to embed high-fidelity, non-linear understanding of processes within the search for optimal values. Due to lack of tractable algebraic equations, the link between simulations and optimization is oftentimes a surrogate metamodel. However, several forms of uncertainty exist within the cycle that links simulation data, to metamodels, to optimization. Uncertainty may originate from parameters of the simulation, or the form and fitted parameters of the metamodel. This paper reviews different literatures that are relevant to surrogate-based optimization and proposes different strategies for handling uncertainty, by combining machine learning with stochastic programming, robust optimization, and discrepancy modeling. We show that incorporating uncertainty management within simulation-based optimization leads to more robust solutions, which protect the decision-maker from infeasible solutions. Here, we present the results of our proposed approaches through a case study for direct-air capture through temperature swing adsorption.

42 ENGINEERING↗
DOE OSTI2015

The mimetic finite difference method for elliptic and parabolic problems with a staggered discretization of diffusion coefficient

Numerical schemes for nonlinear parabolic equations based on the harmonic averaging of cell-centered diffusion coefficients break down when some of these coefficients go to zero or their ratio grows. To tackle this problem, we propose new mimetic finite difference schemes that use a staggered discretization of the diffusion coefficient. The primary mimetic operator approximates div (k•); the derived (dual) mimetic operator approximates - ∇(•). The new mimetic schemes preserve symmetry and positive-definiteness of the continuum problem which allows us to use algebraic solvers with optimal complexity. We perform detailed numerical analysis of the new schemes for linear elliptic problems and a specially designed linear parabolic problem that has solution dynamics typical for nonlinear problems. We also show that the new schemes are competitive with the state-of-the-art schemes for steady-state problems but provide much more accurate solution dynamics for the transient problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Oscillator strengths, first-order properties, and nuclear gradients for local ADC(2)

We describe theory and implementation of oscillator strengths, orbital-relaxed first-order properties, and nuclear gradients for the local algebraic diagrammatic construction scheme through second order. The formalism is derived via time-dependent linear response theory based on a second-order unitary coupled cluster model. The implementation presented here is a modification of our previously developed algorithms for Laplace transform based local time-dependent coupled cluster linear response (CC2LR); the local approximations thus are state specific and adaptive. The symmetry of the Jacobian leads to considerable simplifications relative to the local CC2LR method; as a result, a gradient evaluation is about four times less expensive. Test calculations show that in geometry optimizations, usually very similar geometries are obtained as with the local CC2LR method (provided that a second-order method is applicable). As an exemplary application, we performed geometry optimizations on the low-lying singlet states of chlorophyllide a.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Scalar-fluid theories: cosmological perturbations and large-scale structure

Recently a new Lagrangian framework was introduced to describe interactions between scalar fields and relativistic perfect fluids. This allows two consistent generalizations of coupled quintessence models: non-vanishing pressures and a new type of derivative interaction. The implications of these to the formation of cosmological large-scale structure are uncovered here at the linear order. The full perturbation equations in the two cases are derived in a unified formalism and their Newtonian, quasi-static limit is studied analytically. Requiring the absence of an effective sound speed term in the coupled dark matter fluid restricts the Lagrangian to be a linear function of the matter number density. This leaves new potentially viable classes of both algebraically and derivatively interacting models wherein the coupling may impact the background expansion dynamics and imprint new signatures into the large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1980

Alternatives for jet engine control

Alternatives to linear quadratic regulator theory in the linear case are examined along with nonlinear modelling and optimization approaches for global control. Context for the studies has been set by the DYNGEN digital simulator and by models generated for various phases of the F100 Multivariable Control Synthesis Program. With respect to the linear alternatives, the multivariable frequency domain is stressed. Progress is reported in both the direct algebraic approach to exact model matching, by means of stimulating work on the basic computational issues, and in the indirect generalized Nyquist approach. With respect to nonlinear modelling and optimization, the emphasis is twofold: the development of analytical nonlinear models of the jet engine and the use of these models in conjunction with techniques of mathematical programming in order to study global control over nonincremental portions of the flight envelope. The possibility of using tensor methods is explored.

Sain, M. K.↗
NASA NTRS1985

Accretion disk oscillations - A local analysis in a disk of finite thickness

Two types of oscillations are observed to occur in dwarf novae: 'coherent' and 'quasi-periodic' oscillations. These may be associated with the pulsation of the white dwarf or the accretion disk components of the dwarf nova. Here a local (short-wavelength) analysis is utilized to study the oscillation of a self-consistent, two-dimensional model of an accretion disk. The linearized equations describing adiabatic, inviscid, nonaxisymmetric oscillations are used to derive a fifth-order algebraic equation for the (complex) pulsation frequency of the disk. The solutions of this equation for various values of the wavevector k reveal that the disk is capable of supporting (1) a pair of high-frequency acoustic modes (p-modes); (2) a pair of intermediate-frequency modes which may share the characteristics of internal gravity waves (g-modes) and inertial waves; and (3) a mode associated with a dynamical instability (purely imaginary frequency). The role played by the shear in determining the stability or instability of these modes is also considered. Finally, the global oscillation frequencies of the disk are discussed.

Carroll, B. W.↗
NASA NTRS1986

Time scale analysis of a closed-loop discrete optimal control system

A two-time scale discrete control system is considered. The closed-loop optimal linear quadratic (LQ) regulator for the system requires the solution of a full-order algebraic matrix Riccati equation. Alternatively, the original system is decomposed into reduced-order slow and fast subsystems. The closed-loop optimal control of the subsystems requires the solution of two algebraic matrix Riccati equations of order lower than that required for the full-order system. A composite, closed-loop suboptimal control is created from the sum of the slow and fast feedback optimal controls. Numerical results obtained for an aircraft model show a very close agreement between the exact (optimal) solutions and computationally simpler composite (suboptimal) solutions. The main advantage of the method is the considerable reduction in the overall computational requirements for the closed-loop optimal control of digital flight systems.

Naidu, D. S.↗
NASA NTRS1988

Developing integrated parametric planning models for budgeting and managing complex projects

The applicability of integrated parametric models for the budgeting and management of complex projects is investigated. Methods for building a very flexible, interactive prototype for a project planning system, and software resources available for this purpose, are discussed and evaluated. The prototype is required to be sensitive to changing objectives, changing target dates, changing costs relationships, and changing budget constraints. To achieve the integration of costs and project and task durations, parametric cost functions are defined by a process of trapezoidal segmentation, where the total cost for the project is the sum of the various project cost segments, and each project cost segment is the integral of a linearly segmented cost loading function over a specific interval. The cost can thus be expressed algebraically. The prototype was designed using Lotus-123 as the primary software tool. This prototype implements a methodology for interactive project scheduling that provides a model of a system that meets most of the goals for the first phase of the study and some of the goals for the second phase.

Etnyre, Vance A.↗
NASA NTRS1989

Applications of algebraic image operators to model-based vision

A highly structured and compact algebraic representation of grey-level images is expanded. Addition and multiplication are defined for the set of all grey-level images, which can then be described as polynomials of two variables. Utilizing this new algebraic structure, an innovative, efficient edge-detection scheme is devised. A robust method for linear feature extraction is developed by combining the techniques of a Hough transform and a line follower with this new edge detection scheme. The major advantage of this feature extractor is its general, object-independent nature. Target attributes, such as line segment lengths, intersections, angles of intersection, and endpoints are derived by the feature extraction algorithm and employed during model matching. The feature extractor and model matcher are being incorporated into a distributed robot-control system.

Lerner, Bao-Ting↗
NASA NTRS1990

On the numerical simulation of spatial disturbances in blunt-nose flat plate flow

The feasibility of the simulation of the evolution of spatially amplified disturbances in realistic airfoil flows using a standard Navier-Stokes airfoil code is considered for a blunt-nose flat plate flow. A basic stationary flow and its linear-stability-theory characteristics are analyzed. A computational C-grid is generated using an algebraic grid-generation technique. It is concluded that further improvements to the computational grid and to the code are necessary to make the proposed simulations economical.

Laurien, Eckart↗
DOE OSTI2012

ML Enhancements via the Calculation of Rigid Body Modes (RBMs) for Mechanics Problems Implemented within the Albany Code Base

The algebraic multigrid approach known as smoothed aggregation is very efficient at solving systems that arise from elasticity problems. In order to construct an efficient algebraic multilevel method, a multigrid solver should be provided with a small set of vectors that represent the error components that are difficult to resolve. It is well-known that for linear elasticity problems, these components correspond to the the so-called rigid body modes (RBMs). The present document summarizes some new development within the Albany code base that has enabled the application of algebraic multigrid preconditioners from the ML package of Trilinos to mechanics problems implemented within Albany via a new function that calculates the RBMs using information about the problem's underlying mesh. The performance of these preconditioners is evaluated on four problems: a 3D static elasticity problem, a 3D non-linear elasticity problem, a 3D thermo-elasticity problem, and a 3D thermo-poro-plasticity problem. The tests reveal the superiority of the ML preconditioners over ILU preconditioners from the Trilinos Ifpack package for mechanics problems in Albany.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

On the scalability of the Albany/FELIX first-order Stokes approximation ice sheet solver for large-scale simulations of the Greenland and Antarctic ice sheets

We examine the scalability of the recently developed Albany/FELIX finite-element based code for the first-order Stokes momentum balance equations for ice flow. We focus our analysis on the performance of two possible preconditioners for the iterative solution of the sparse linear systems that arise from the discretization of the governing equations: (1) a preconditioner based on the incomplete LU (ILU) factorization, and (2) a recently-developed algebraic multigrid (AMG) preconditioner, constructed using the idea of semi-coarsening. A strong scalability study on a realistic, high resolution Greenland ice sheet problem reveals that, for a given number of processor cores, the AMG preconditioner results in faster linear solve times but the ILU preconditioner exhibits better scalability. A weak scalability study is performed on a realistic, moderate resolution Antarctic ice sheet problem, a substantial fraction of which contains floating ice shelves, making it fundamentally different from the Greenland ice sheet problem. Here, we show that as the problem size increases, the performance of the ILU preconditioner deteriorates whereas the AMG preconditioner maintains scalability. This is because the linear systems are extremely ill-conditioned in the presence of floating ice shelves, and the ill-conditioning has a greater negative effect on the ILU preconditioner than on the AMG preconditioner.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Relative control effectiveness technique with application to airplane control coordination

A method to select optimal combinations of the control variables of a linear system is reported. The combinations are chosen so that the control channels have their principal influences on selected fundamental modes of the system. A series of algebraic maximization problems is used to maximize the effects of the control channels on selected modes while simultaneously minimizing the effects on the remaining modes. The method is applied to the lateral and directional control of a linearized airplane model having ailerons, a rudder, and differential tail surfaces. Integration of these control eliminates oscillations present in the roll rate for a step lateral-control input and improves the sideslip response with reduced rolling motions for a step directional-control input. Inclusion of thrust-vectoring engine nozzles improves the roll rate capability of the airplane.

Lallman, F. J.↗
DOE OSTI2014

Radiative observables for linearized gravity on asymptotically flat spacetimes and their boundary induced states

We discuss the quantization of linearized gravity on globally hyperbolic, asymptotically flat, vacuum spacetimes, and the construction of distinguished states which are both of Hadamard form and invariant under the action of all bulk isometries. The procedure, we follow, consists of looking for a realization of the observables of the theory as a sub-algebra of an auxiliary, non-dynamical algebra constructed on future null infinity ℱ⁺. The applicability of this scheme is tantamount to proving that a solution of the equations of motion for linearized gravity can be extended smoothly to ℱ⁺. This has been claimed to be possible provided that a suitable gauge fixing condition, first written by Geroch and Xanthopoulos [“Asymptotic simplicity is stable,” J. Math. Phys. 19, 714 (1978)], is imposed. We review its definition critically, showing that there exists a previously unnoticed obstruction in its implementation leading us to introducing the concept of radiative observables. These constitute an algebra for which a Hadamard state induced from null infinity and invariant under the action of all spacetime isometries exists and it is explicitly constructed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗