SEARCH · Engineering Papers

Results for “matrix approximation”

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 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Fe/Mn in olivine of carbonaceous meteorites

Olivines in primitive meteorites show a range of Fe/Mn both within one grain and among grains suggesting that they have recorded changing conditions during or after growth. Because olivine should be an early forming phase, Fe/Mn is used here to infer these earliest conditions. Initial Fe/Mn in cores of isolated, euhedral forsterite in both C2 and C3 meteorites ranges from 25 to 35 but differs at grain edge. Murchison (C2) forsterites show Fe/Mn approaching 1.0 at the grain edge while Ornans Fe/Mn is near 60 at grain edge. These values are lower than the matrix Fe/Mn for both meteorites and the distinct difference in zoning profile indicates different processes operating during and after grain growth. The Fe/Mn of bulk samples from a particular source such as the Moon is nearly constant. Individual samples show variation suggesting that there is some fractionation of Mn from Fe. Minerals have their individual ranges of Fe/Mn which has been used to recognize different types of olivine within one meteorite. Extreme values of Fe/Mn below 1.0 occur in forsterite from some IDP's, UOC matrix, and C1 meteorites. There are apparently no detailed studies of Fe/Mn variation within single olivine grains. Forsterite grains in C2 and C3 carbonaceous chondrites show complex zoning, and the nearly pure forsterites (Fo greater than approximately 99.5) have high levels of some minor elements including Ti, Al, V, and Sc. There is disagreement on the original source of these grains and both chondrule and vapor growth have been proposed. In addition, there is clear evidence that diffusion has affected the outer margins but in some cases the whole grain. Within the cores, the FeO range is limited, and if growing under constant conditions, the Fe/Mn should be near constant as there is little fractionation of Mn from Fe by forsterite. Additionally, there are apparently no co-crystallizing phases as evidenced by a lack of common inclusions in the forsterites. These observations are now followed by analyses of isolated olivine grains in C2 and C3 meteorites.

Steele, Ian M.↗

NASA NTRS2007

Atomic Data and Spectral Line Intensities for CA XVII

Electron impact collision strengths, energy levels, oscillator strengths and spontaneous radiative decay rates are calculated for Ca XVII. The configurations used are 2s(sup 2), 2s2p, 2p(sup 2), 2l3l', 214l' and 2s5l', with l = s,p and l' = s,p, d giving rise to 92 fine-structure levels in intermediate coupling. Collision strengths are calculated at seven incident energies (15, 30, 75, 112.5, 150, 187.5 and 225 Ry) for the transitions within the three lowest configurations corresponding to the 10 lowest energy levels, and five incident energies (75, 112.5, 150, 187.5 and 225 Ry) for transitions between the lowest five levels and the n = 3,4,5 configurations. Calculations have been carried out using the distorted wave approximation. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates of the present work, and R-Matrix results for the 2s2, 2s2p, 2p2 configurations available in the literature, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8)-10(exp 14)/cu cm at an electron temperature of log Te(K)=6.7, corresponding to the maximum abundance of Ca XVII. Spectral line intensities are calculated, and their diagnostic relevance L; discussed. This dataset will be made available in the next version of the CHIANTI database.

Bhatia, A.K.↗

DOE OSTI2023

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗

NASA NTRS1985

The accuracy of approximate solutions in the analysis of fracture of composites

This paper concerns the accuracy of three related mathematical models (developed by Hedgepeth, Eringen and Sendeckyj and Jones) used in the stress analysis and in fracture studies of continuous-fiber composites. These models have particular application in the investigation of fiber and matrix stresses in unidirectional composites in the region near a crack tip. The interest in such models is motivated by the desire to be able to simplify the equations of elasticity to the point that they can be solved in a relatively easy manner.

Goree, J. G.↗

NASA NTRS1985

The accuracy of approximate solutions in the analysis of fracture of composites

Goree, J. G.↗

NASA NTRS1977

Vibrational and rotational transitions in low-energy electron-diatomic-molecule collisions. I - Close-coupling theory in the moving body-fixed frame. II - Hybrid theory and close-coupling theory: An /l subscript z-prime/-conserving close-coupling approximation

A detailed vibrational-rotational (V-R) close-coupling formulation of electron-diatomic-molecule scattering is developed in which the target molecular axis is chosen to be the z-axis and the resulting coupled differential equation is solved in the moving body-fixed frame throughout the entire interaction region. The coupled differential equation and asymptotic boundary conditions in the body-fixed frame are given for each parity, and procedures are outlined for evaluating V-R transition cross sections on the basis of the body-fixed transition and reactance matrix elements. Conditions are discussed for obtaining identical results from the space-fixed and body-fixed formulations in the case where a finite truncated basis set is used. The hybrid theory of Chandra and Temkin (1976) is then reformulated, relevant expressions and formulas for the simultaneous V-R transitions of the hybrid theory are obtained in the same forms as those of the V-R close-coupling theory, and distorted-wave Born-approximation expressions for the cross sections of the hybrid theory are presented. A close-coupling approximation that conserves the internuclear axis component of the incident electronic angular momentum (l subscript z-prime) is derived from the V-R close-coupling formulation in the moving body-fixed frame.

Choi, B. H.↗

DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Recent Advances toward Efficient Calculation of Higher Nuclear Derivatives in Quantum Chemistry

In this article, we provide an overview of state-of-the-art techniques that are being developed for efficient calculation of second and higher nuclear derivatives of quantum mechanical (QM) energy. Calculations of nuclear Hessians and anharmonic terms incur high costs and memory and scale poorly with system size. Three emerging classes of methods—machine learning (ML), automatic differentiation (AD), and matrix completion (MC)—have demonstrated promise in overcoming these challenges. We illustrate studies that employ unsupervised ML methods to reduce the need for multiple Hessian calculations in dynamics simulations and those that utilize supervised ML to construct approximate potential energy surfaces and estimate Hessians and anharmonic terms at reduced cost. By extension, if electronic structure operations could be written in a manner similar to functions underlying ML methods, rapid differentiation or AD routines can be employed to inexpensively calculate higher arbitrary-order derivatives. While ML approaches are typically black-box, we describe methods such as compressed sensing (CS) and MC, which explicitly leverage problem-specific mathematical properties of higher derivatives such as sparsity and low-rank, to complete higher derivative information using only a small, incomplete sample. The three classes of methods facilitate reliable predictions of observables ranging from infrared spectra to thermal conductivity and constitute a promising way forward in accurately capturing otherwise intractable higher-order responses of QM energy to nuclear perturbations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2024

FEED Study of CarbonCapture Inc DAC and CarbonCure Utilization Technologies Using United States Steel’s Gary Works Plant Waste Heat (Final Report)

The University of Illinois at Urbana-Champaign (UIUC) led this project to produce a front-end engineering design (FEED) study of an advanced Direct Air Capture and Utilization System (DACUS) system that can remove a minimum of 5,000 tonnes/yr net of carbon dioxide from air (based on cradle-to-gate LCA) and utilizing the CO 2 to produce low carbon intensity ready mix concrete. The designed system, if built, would be larger than any currently existing Direct Air Capture (DAC) collector in the U.S. Such carbon capture technologies are critical to meeting the goals of the DOE’s program to accelerate climate-critical technology. In addition to the power sector, industrial facilities for the manufacture of steel and cement/concrete are among the major sources of anthropogenic CO 2 . DAC is a promising new technology for reducing CO 2 , a potent greenhouse gas, in the atmosphere but is expensive, in part due to the energy required to adsorb and desorb captured CO 2 during cycles. By integrating CarbonCapture Inc. (CCI) DAC modules at United States Steel's Gary Works (USS) and utilizing the site's waste heat, energy, and location this project evaluates the feasibility of utilizing the captured CO 2 and the logistics of transportation. CarbonCapture Inc. has developed an innovative DAC system using novel adsorbents to cost-effectively capture CO 2 . The captured, liquified gas will be trucked to ready-mix concrete plants within the region, the closest of which is approximately 3.5 miles away, where CarbonCure will inject it into concrete during the mixing process at the facilities. The carbon dioxide reacts with concrete, mineralizing into calcium carbonate (CaCO 3 ), permanently locking the greenhouse gas into the matrix of the building material. This FEED study demonstrated a full CO 2 value chain for DACUS from industrial facilities. It also provided a means for Visage Energy Corp. (Visage) to assess the impact of this holistic approach on job creation, regional economic impact, and environmental justice. The project team also included Sargent & Lundy (S&L) to provide the constructability review and costing of the integration of the DAC with the steel plant. Ecotek Engineering USA, LLC designed the outside battery limit (OSBL) infrastructure to connect the DAC and the plant. Activities performed during the project included: (1) Project Management Plan; (2) Technology Maturation Plan (TMP); (3) Initial Workforce Readiness Plan; (4) Workforce Readiness Plan; (5) Front-End Engineering Design (FEED) Study; (6) Project Design Basis; (7) Hazards and Operability (HAZOP) Study; (8) Constructability Review; (9) Project Cost Assessment; (10) Logistics Analysis of CO 2 Transportation to the Utilization Site; (11) Business Case Analysis; (12) Life Cycle Analysis (LCA); (13) Environmental Health and Safety (EH&S) Analysis; (14) Environmental Justice Analysis; and (15) Economic Revitalization and Job Creation Outcomes Analysis. This report provides a summary of the outcomes and results of the project, which was performed between Oct. 1, 2022, through Sept. 30, 2024.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1984

A finite element method for nonlinear forced vibrations of rectangular plates

The finite element method has been extended to determine the response of large amplitude forced vibrations of thin plates. A harmonic force matrix of a rectangular element under uniform harmonic excitation is developed for nonlinear forced vibration analysis. Inplane deformation and inertia are both considered in the formulation. Results obtained are compared with simple elliptic response, perturbation and other approximation solutions.

Mei, C.↗

NASA NTRS1986

Implicit solution methods in computational fluid dynamics

A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented. The methods are demonstrated in fully vectorized codes for a CRAY type architecture. The emphasis will be on the Euler equations in generalized coordinates.

Pulliam, Thomas H.↗

NASA NTRS1987

Self-lubricating coatings for high-temperature applications

Some present-day aeropropulsion systems impose severe demands on the thermal and oxidative stability of lubricant, bearing, and seal materials. These demands will be much more severe for operational systems around the turn of the century. Solid lubricants with maximum temperature capabilities of about 1100 C are known. Unfortunately, none of the solid lubricants with the highest temperature capabilities are effective below approximately 400 C. However, research shows that silver and stable fluorides, such as calcium and barium fluoride act synergistically to provide lubrication from below room temperature to approximately 900 C. Plasma-sprayed, self-lubricating composite coatings that were developed at Lewis are described. Background information is given on coatings, designed as PS100 and PS101, that contain the solid lubricants in a Nichrome matrix. These coatings have low friction coefficients over a wide temperature range, but they have inadequate wear resistance for some long-duration applications. Wear resistance was dramatically improved in a recently developed coating PS200, by replacing the Nichrome matrix material with metal-bonded chromium carbide containing dispersed silver and calcium fluoride/barium fluoride eutectic (CaF2/BaF2). The lubricants control friction and the carbide matrix provides excellent wear resistance. Successful tests of these coatings are discussed.

Sliney, Harold E.↗

NASA NTRS1988

Thermal stresses in composite tubes using complementary virtual work

This paper addresses the computation of thermally induced stresses in layered, fiber-reinforced composite tubes subjected to a circumferential gradient. The paper focuses on using the principle of complementary virtual work, in conjunction with a Ritz approximation to the stress field, to study the influence on the predicted stresses of including temperature-dependent material properties. Results indicate that the computed values of stress are sensitive to the temperature dependence of the matrix-direction compliance and matrix-direction thermal expansion in the plane of the lamina. There is less sensitivity to the temperature dependence of the other material properties.

Hyer, M. W.↗

NASA NTRS1990

Finite element mesh refinement criteria for stress analysis

This paper discusses procedures for finite-element mesh selection and refinement. The objective is to improve accuracy. The procedures are based on (1) the minimization of the stiffness matrix race (optimizing node location); (2) the use of h-version refinement (rezoning, element size reduction, and increasing the number of elements); and (3) the use of p-version refinement (increasing the order of polynomial approximation of the elements). A step-by-step procedure of mesh selection, improvement, and refinement is presented. The criteria for 'goodness' of a mesh are based on strain energy, displacement, and stress values at selected critical points of a structure. An analysis of an aircraft lug problem is presented as an example.

Kittur, Madan G.↗

NASA NTRS1998

Accelerated Training for Large Feedforward Neural Networks

In this paper we introduce a new training algorithm, the scaled variable metric (SVM) method. Our approach attempts to increase the convergence rate of the modified variable metric method. It is also combined with the RBackprop algorithm, which computes the product of the matrix of second derivatives (Hessian) with an arbitrary vector. The RBackprop method allows us to avoid computationally expensive, direct line searches. In addition, it can be utilized in the new, 'predictive' updating technique of the inverse Hessian approximation. We have used directional slope testing to adjust the step size and found that this strategy works exceptionally well in conjunction with the Rbackprop algorithm. Some supplementary, but nevertheless important enhancements to the basic training scheme such as improved setting of a scaling factor for the variable metric update and computationally more efficient procedure for updating the inverse Hessian approximation are presented as well. We summarize by comparing the SVM method with four first- and second- order optimization algorithms including a very effective implementation of the Levenberg-Marquardt method. Our tests indicate promising computational speed gains of the new training technique, particularly for large feedforward networks, i.e., for problems where the training process may be the most laborious.

Stepniewski, Slawomir W.↗

NASA NTRS2014

Fuzzy Model-based Pitch Stabilization and Wing Vibration Suppression of Flexible Wing Aircraft.

This paper presents a fuzzy nonlinear controller to regulate the longitudinal dynamics of an aircraft and suppress the bending and torsional vibrations of its flexible wings. The fuzzy controller utilizes full-state feedback with input constraint. First, the Takagi-Sugeno fuzzy linear model is developed which approximates the coupled aeroelastic aircraft model. Then, based on the fuzzy linear model, a fuzzy controller is developed to utilize a full-state feedback and stabilize the system while it satisfies the control input constraint. Linear matrix inequality (LMI) techniques are employed to solve the fuzzy control problem. Finally, the performance of the proposed controller is demonstrated on the NASA Generic Transport Model (GTM).

Stability and Control↗

NASA NTRS2014

Processing and Characterization of Carbon Nanotube Composites

Recent advances in the synthesis of large-scale quantities of carbon nanotubes (CNT) have provided the opportunity to study the mechanical properties of polymer matrix composites using these novel materials as reinforcement. Nanocomp Technologies, Inc. currently supplies large sheets with dimensions up to 122 cm x 244 cm containing both single-wall and few-wall CNTs. The tubes are approximately 1 mm in length with diameters ranging from 8 to 12 nm. In the present study being conducted at NASA Langley Research Center (LaRC), single and multiple layers of CNT sheets were infused or coated with various polymer solutions that included commercial toughened-epoxies and bismaleimides, as well as a LaRC developed polyimide. The resulting CNT composites were tested in tension using a modified version of ASTM D882-12 to determine their strength and modulus values. The effects of solvent treatment and mechanical elongation/alignment of the CNT sheets on the tensile performance of the composite were determined. Thin composites (around 50 wt% CNT) fabricated from acetone condensed and elongated CNT sheets with either a BMI or polyimide resin solution exhibited specific tensile moduli approaching that of toughened epoxy/ IM7 carbon fiber unidirectional composites.

Can, Roberto J.↗

DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗