Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearization”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Investigations of the linear and non-linear flow harmonics using the a multi-phase transport model

The multi-phase transport model (AMPT) is used to study the effects of the parton-scattering cross-sections (σ pp ) and hadronic re-scattering on the linear contributions to the flow harmonic v 4 , the non-linear response coefficients, and the correlations between different order flow symmetry planes in Au + Au collisions at 200 GeV. The model results, which agree with current experimental measurements, indicate that the higher-order flow harmonics are sensitive to the σ pp variations. However, the non-linear response coefficients and the correlations between different order flow symmetry planes are σ pp independent. Furthermore, these results suggest that further detailed experimental measurements which span a broad range of collision systems and beam energies could serve as an additional constraint for the theoretical models' calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Magnetic potential based formulation for linear and non-linear 3D RF sheath simulation

Abstract This paper reports a new numerical scheme to simulate the radio-frequency (RF) induced RF sheath, which is suitable for a large 3D simulation. In the RF sheath boundary model, the tangential component of the electric field ( E t ) is given by the gradient of a scalar electric field potential. We introduce two additional scalar potentials for the tangential components of the magnetic field, which effectively impose the normal electric displacement ( D n ) on the plasma sheath boundary condition via in-homogeneous Neumann boundary condition and constrain the tangential electric field on the surface as curl-free ( ∇ × E t = 0 ). In our approach, the non-linear sheath impedance is formulated as a natural extension of the large thickness (or asymptotic) sheath limit ( D n = 0 ), allowing for handling both asymptotic and non-linear regimes seamlessly. The new scheme is implemented using the Petra-M finite element method analysis framework and is verified with simulations in the literature. The significance of non-linearity is discussed in various plasma conditions. An application of this scheme to asymptotic RF sheath simulation on the WEST ICRF antenna side limiters is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Transformation matrices between non-linear and linear differential equations

In the linearization of systems of non-linear differential equations, those systems which can be exactly transformed into the second order linear differential equation Y"-AY'-BY=0 where Y, Y', and Y" are n x 1 vectors and A and B are constant n x n matrices of real numbers were considered. The 2n x 2n matrix was used to transform the above matrix equation into the first order matrix equation X' = MX. Specially the matrix M and the conditions which will diagonalize or triangularize M were studied. Transformation matrices P and P sub -1 were used to accomplish this diagonalization or triangularization to return to the solution of the second order matrix differential equation system from the first order system.

Sartain, R. L.↗
NASA NTRS1985

On the adequacy of meridional resolution of linear and quasi-linear barotropic models

The minimum meridional resolution needed for an adequate numerical simulation of the linear and 'quasi-linear' barotropic vorticity dynamics in the vicinity of a critical latitude (CL) is determined by using a semispectral nondivergent barotropic model on a sphere. The behavior of the analytic as well as the meridionally finite differenced form of the refractive index in the vicinity of a CL is discussed. An estimation of an upper bound on the latitudinal gridsize to be used in linear calculations is given, and constraint on the latitudinal grid size are considered. The dynamical model utilized is briefly described.

Nigam, S.↗
NASA NTRS1990

A linear quadratic regulator approach to the stabilization of uncertain linear systems

This paper presents a linear quadratic regulator approach to the stabilization of uncertain linear systems. The uncertain systems under consideration are described by state equations with the presence of time-varying unknown-but-bounded uncertainty matrices. The method is based on linear quadratic regulator (LQR) theory and Liapunov stability theory. The robust stabilizing control law for a given uncertain system can be easily constructed from the symmetric positive-definite solution of the associated augmented Riccati equation. The proposed approach can be applied to matched and/or mismatched systems with uncertainty matrices in which only their matrix norms are bounded by some prescribed values and/or their entries are bounded by some prescribed constraint sets. Several numerical examples are presented to illustrate the results.

Shieh, L. S.↗
NASA NTRS1991

Linear prediction filters for linear and nonlinear modeled geomagnetic activity

It is shown that the Faraday loop analog model of geomagnetic activity exhibits both the directly driven and loading-unloading magnetospheric responses to solar wind input. It is further shown that the directly driven component is a linear response to loading while the loading-unloading response is nonlinear. Linear prediction filters which relate model input to output are discussed. By either allowing or suppressing the loading-unloading model response filters that relate to nonlinear or linear dynamics, respectively, have been computed. Filters that described the directly driven response are finite ranged; they asymptote to zero with increasing lag on a time scale that is fixed by the dissipation rate of the model. Filters that describe the nonlinear total response are infinite ranged; they asymptote with increasing lag to large amplitude periodic oscillations. Some implications of these infinite ranged filters are discussed.

Klimas, A. J.↗
NASA NTRS1992

Linear and non-linear failure analysis of composite laminates with transverse shear

A finite element computational procedure has been developed to find linear and nonlinear (von Karman) first-ply failure loads of composite laminates subjected to in-plane and transverse loads. The FEM is based on first-order shear deformation theory and several phenomenological failure criteria. It is found that the failure loads predicted by various failure criteria differ from one another by a maximum of 35 percent in the case of linear loads and 50 percent in the case of nonlinear loads; the failure locations differ from each other in a random way. It is demonstrated that the difference between the linear and nonlinear failure loads is large in the case of transversely loaded simply supported laminates and thin plates. The difference is found to be quite small in the case of in-plane (tensile) loading and thick laminates.

Reddy, Y. S. N.↗
NASA NTRS1992

Efficient yet accurate solution of the linear transport equation in the presence of internal sources - The exponential-linear-in-depth approximation

The present solutions to the linear transport equation pertain to monoenergetic particles' interaction with a multiple scattering/absorbing layered medium with a general anisotropic internal source term. Attention is given to a novel exponential-linear approximation to the internal source, as a function of scattering depth, which furnishes an at-once efficient and accurate solution to the linear transport equation through its reduction of the spatial mesh size. The great superiority of the proposed method is demonstrated by the numerical results obtained in the illustrative cases of (1) an embedded thermal source and (2) a rapidly varying beam pseudosource.

Kylling, Arve↗
NASA NTRS1998

Trellises and Trellis-Based Decoding Algorithms for Linear Block Codes: An Iterative Decoding Algorithm for Linear Block Codes Based on a Low-Weight Trellis Search - Part 3

For long linear block codes, maximum likelihood decoding based on full code trellises would be very hard to implement if not impossible. In this case, we may wish to trade error performance for the reduction in decoding complexity. Sub-optimum soft-decision decoding of a linear block code based on a low-weight sub-trellis can be devised to provide an effective trade-off between error performance and decoding complexity. This chapter presents such a suboptimal decoding algorithm for linear block codes. This decoding algorithm is iterative in nature and based on an optimality test. It has the following important features: (1) a simple method to generate a sequence of candidate code-words, one at a time, for test; (2) a sufficient condition for testing a candidate code-word for optimality; and (3) a low-weight sub-trellis search for finding the most likely (ML) code-word.

Lin, Shu↗
DOE OSTI2017

Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory

We present two efficient iterative algorithms for solving the linear response eigen- value problem arising from the time dependent density functional theory. Although the matrix to be diagonalized is nonsymmetric, it has a special structure that can be exploited to save both memory and floating point operations. In particular, the nonsymmetric eigenvalue problem can be transformed into a product eigenvalue problem that is self-adjoint with respect to a K-inner product. This product eigenvalue problem can be solved efficiently by a modified Davidson algorithm and a modified locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm that make use of the K-inner product. The solution of the product eigenvalue problem yields one component of the eigenvector associated with the original eigenvalue problem. However, the other component of the eigenvector can be easily recovered in a postprocessing procedure. Therefore, the algorithms we present here are more efficient than existing algorithms that try to approximate both components of the eigenvectors simultaneously. The efficiency of the new algorithms is demonstrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

PB-AM: An open-source, fully analytical linear poisson-boltzmann solver

Here, we present the open source distributed software package Poisson-Boltzmann Analytical Method (PB-AM), a fully analytical solution to the linearized PB equation, for molecules represented as non-overlapping spherical cavities. The PB-AM software package includes the generation of outputs files appropriate for visualization using visual molecular dynamics, a Brownian dynamics scheme that uses periodic boundary conditions to simulate dynamics, the ability to specify docking criteria, and offers two different kinetics schemes to evaluate biomolecular association rate constants. Given that PB-AM defines mutual polarization completely and accurately, it can be refactored as a many-body expansion to explore 2- and 3-body polarization. Additionally, the software has been integrated into the Adaptive Poisson-Boltzmann Solver (APBS) software package to make it more accessible to a larger group of scientists, educators, and students that are more familiar with the APBS framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Analysis of Monte Carlo accelerated iterative methods for sparse linear systems: Analysis of Monte Carlo accelerated iterative methods for sparse linear systems

Here, we consider hybrid deterministic-stochastic iterative algorithms for the solution of large, sparse linear systems. Starting from a convergent splitting of the coefficient matrix, we analyze various types of Monte Carlo acceleration schemes applied to the original preconditioned Richardson (stationary) iteration. We expect that these methods will have considerable potential for resiliency to faults when implemented on massively parallel machines. We also establish sufficient conditions for the convergence of the hybrid schemes, and we investigate different types of preconditioners including sparse approximate inverses. Numerical experiments on linear systems arising from the discretization of partial differential equations are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Kinetic-j : A computational kernel for solving the linearized Vlasov equation applied to calculations of the kinetic, configuration space plasma current for time harmonic wave electric fields

We present the Kinetic-j code, a computational kernel for evaluating the linearized Vlasov equation withapplication to calculating the kinetic plasma response (current) to an applied time harmonic wave electricfield. This code addresses the need for a configuration space evaluation of the plasma current to enablekinetic full-wave solvers for waves in hot plasmas to move beyond the limitations of the traditional Fourierspectral methods. We benchmark the kernel via comparison with the standard k-space forms of the hotplasma conductivity tensor.

Green, David L.↗
DOE OSTI2017

Linear permeability evolution of expanding conduits due to feedback between flow and fast phase change: Linear Permeability Evolution of Conduit

Conduits are ubiquitous and critical pathways for many fluids relevant for geophysical processes such as magma, water, and gases. Predicting flow through conduits is challenging when the conduit geometry coevolves with the flow. We theoretically show that the permeability (k) of a conduit whose walls are eroding due to fast phase change increases linearly with time because of a self-reinforcing mechanism. This simple result is surprising given complex feedbacks between flow, transport, and phase change. The theory is congruent with previous experimental observations of fracture dissolution in calcite. Supporting computational fracture dissolution experiments showed that k only slightly increases until the dissolution front reaches the narrowest conduit constriction, after which the linear evolution of k manifests. Finally, the theory holds across multiple scales and a broad range of Peclet and Damkohler numbers and thus advances the prediction of dynamic mass fluxes through expanding conduits in various geologic and environmental settings

58 GEOSCIENCES↗
DOE OSTI2017

GPU implementation of the linear scaling three dimensional fragment method for large scale electronic structure calculations

LS3DF, namely linear scaling three-dimensional fragment method, is an efficient linear scaling ab initio total energy electronic structure calculation code based on a divide-and-conquer strategy. In this paper, we present our GPU implementation of the LS3DF code. Our test results show that the GPU code can calculate systems with about ten thousand atoms fully self-consistently in the order of 10 min using thousands of computing nodes. This makes the electronic structure calculations of 10,000-atom nanosystems routine work. This speed is 4.5–6 times faster than the CPU calculations using the same number of nodes on the Titan machine in the Oak Ridge leadership computing facility (OLCF). Such speedup is achieved by (a) carefully re-designing of the computationally heavy kernels; (b) redesign of the communication pattern for heterogeneous supercomputers.

Computer Science↗