Engineering Papers⌕ Search

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.

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

At least 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Collisional redistribution of radiation. II - The effects of degeneracy on the equations of motion for the density matrix. III - The equation of motion for the correlation function and the scattered spectrum

The effect of correlations between an absorber atom and perturbers in the binary-collision approximation are applied to degenerate atomic systems. A generalized absorption profile which specifies the final state of the atom after an absorption event is related to the total intensities of Rayleigh scattering and fluorescence from the atom. It is suggested that additional dynamical information to that obtainable from ordinary absorption experiments is required in order to describe redistributed atomic radiation. The scattering of monochromatic radiation by a degenerate atom is computed in a binary-collision approximation; an equation of motion is derived for the correlation function which is valid outside the quantum-regression regime. Solutions are given for the weak-field conditions in terms of generalized absorption and emission profiles that depend on the indices of the atomic multipoles.

Burnett, K.↗
DOE OSTI2022

Detailed analysis of excited-state systematics in a lattice QCD calculation of 𝑔 𝐴

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations ≳ 1 fm are commonly used to combat these issues. In this work, focusing on g A , we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2–1 fm), combined with a multistate analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while, for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Efficient discontinuous Galerkin implementations and preconditioners for implicit unsteady compressible flow simulations

This work presents and compares efficient implementations of high-order discontinuous Galerkin methods: a modal matrix-free discontinuous Galerkin (DG) method, a hybridizable discontinuous Galerkin (HDG) method, and a primal formulation of HDG, applied to the implicit solution of unsteady compressible flows. The matrix-free implementation allows for a reduction of the memory footprint of the solver when dealing with implicit time-accurate discretizations. HDG reduces the number of globally-coupled degrees of freedom relative to DG, at high order, by statically condensing element-interior degrees of freedom from the system in favor of face unknowns. The primal formulation further reduces the element-interior degrees of freedom by eliminating the gradient as a separate unknown. This paper introduces a p-multigrid preconditioner implementation for these discretizations and presents results for various flow problems. Benefits of the p-multigrid strategy relative to simpler, less expensive, preconditioners are observed for stiff systems, such as those arising from low-Mach number flows at high-order approximation. The p-multigrid preconditioner also shows excellent scalability for parallel computations. Finally, additional savings in both speed and memory occur with a matrix-free/reduced version of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Upscaling Reactive Transport and Clogging in Shale Microcracks by Deep Learning

Fracture networks in shales exhibit multiscale features. A rock system may contain a few main fractures and thousands of microcracks, whose length and aperture are orders of magnitude smaller than the former. It is computationally prohibitive to resolve all the fractures explicitly for such multiscale fracture networks. One traditional approach is to model the small-scale features (e.g., microcracks in shales) as an effective medium. Although this fracture-matrix conceptualization significantly reduces the problem complexity, there are classes of physical processes that cannot be accurately upscaled by effective medium approximations, e.g., microcrack clogging during mineral reactions. Here, we employ deep learning in place of effective medium theory to upscale physical processes in small-scale features. Specifically, we consider reactive transport in a fracture-microcrack network where microcracks can be clogged by precipitation. A deep learning multiscale algorithm is developed, in which the microcracks are upscaled as a wall boundary condition of the main fractures. The wall boundary condition is constructed by recurrent neural networks, which take concentration histories as input and predict the solute transport from main fractures to microcracks. The deep learning multiscale algorithm is firstly employed in specific scenarios, then a general model is developed which can work under various conditions. The new approach is validated against fully resolved simulations and an analytical solution, providing a reliable and efficient solution for problems that cannot be upscaled by effective medium models.

58 GEOSCIENCES↗
DOE OSTI2023

Data assimilation in operator algebras

We develop an algebraic framework for sequential data assimilation of partially observed dynamical systems. In this framework, Bayesian data assimilation is embedded in a nonabelian operator algebra, which provides a representation of observables by multiplication operators and probability densities by density operators (quantum states). In the algebraic approach, the forecast step of data assimilation is represented by a quantum operation induced by the Koopman operator of the dynamical system. Moreover, the analysis step is described by a quantum effect, which generalizes the Bayesian observational update rule. Projecting this formulation to finite-dimensional matrix algebras leads to computational schemes that are i) automatically positivity-preserving and ii) amenable to consistent data-driven approximation using kernel methods for machine learning. Moreover, these methods are natural candidates for implementation on quantum computers. Applications to the Lorenz 96 multiscale system and the El Niño Southern Oscillation in a climate model show promising results in terms of forecast skill and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

3 He +𝛼 resonances in 7 Be

Resonances in 7 Be which decay into the 3 He+α exit channel have been measured with improved precision using preexisting data sets. The energy and width of the J π =7/2 - state have been extracted from an invariant-mass study of projectile-breakup products originating from interactions of an E/A=10.7-MeV 10 C beam on Be and C targets. The excitation energy of this state (from the pole of the S-matrix) is determined to be E*=4.545(6)~MeV, a factor of 8 improvement in precision as compared to the ENSDF value. This improvement is enabled by fine tuning the detector calibrations using calibration resonances in 6 Li, 7 Li, 6 Be, 9 B, and 12 C whose decay energies are known to high precision. The J π =7/2- resonance in 7 Be can now itself be used as a calibration resonance in invariant-mass experiments. This utility is particularly helpful for 3 He energy calibrations of CsI(Tl) detectors which are often used in detector arrays employed for measurements with fast beams. This utility is demonstrated with a data set associated with E/A=70 MeV 7 Be beams which are inelastically excited to the 7/2 - and 5/2 - 1 resonances. Again using the pole of the S-matrix as the definition of the resonance parameters, the fitted excitation energy of the 5/2 - 1 resonance is 6.376(17) MeV, approximately 300 keV lower than the ENSDF value. Finally, its width of 565(4) keV is roughly half of the ENSDF value.

energy levels↗
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 OSTI2022

Misanthropic entropy and renormalization as a communication channel

A central physical question is the extent to which infrared (IR) observations are sufficient to reconstruct a candidate ultraviolet (UV) completion. We recast this question as a problem of communication, with messages encoded in field configurations of the UV being transmitted to the IR degrees of freedom via a noisy channel specified by renormalization group (RG) flow, with noise generated by coarse graining/decimation. We present an explicit formulation of these considerations in terms of lattice field theory, where we show that the “misanthropic entropy” — the mutual information obtained from decimating neighbors — encodes the extent to which information is lost in marginalizing over/tracing out UV degrees of freedom. Our considerations apply both to statistical field theories as well as density matrix renormalization of quantum systems, where in the quantum case, the statistical field theory analysis amounts to a leading-order approximation. As a concrete example, we focus on the case of the 2D Ising model, where we show that the misanthropic entropy detects the onset of the phase transition at the Ising model critical point.

Physics↗
NASA NTRS1970

Pyrotechnic hazards classification and evaluation program. Phase 2, segment 1: Records and experience analysis

A comprehensive search review and analysis was made of various technical documents relating to both pyrotechnics and high explosives testing, handling, storage, manufacturing, physical and chemical characteristics and accidents and incidents. Of approximately 5000 technical abstracts reviewed, 300 applicable documents were analyzed in detail. These 300 documents were then converted to a subject matrix so that they may be readily referenced for application to the current programs. It was generally concluded that information in several important categories was lacking. Two of the more important categories were in pyrotechnics sensitivity testing and TNT equivalency testing. A general recommendation resulting from this study was that this activity continue and a comprehensive data bank be generated that would allow immediate access to a large volume of pertinent information in a relatively short period of time.

Source record↗
NASA NTRS1975

Simulation of aircraft crash and its validation

A joint FAA/NASA program is discussed which is aimed at developing a reliable technology for the design of crashworthy light aircraft. This program encompasses the development of analytical methods, the definition of a survivable crash envelope, and the design of improved seat and restraint systems. A facility for full-scale crash-simulation testing is described along with the test method and results of five full-scale crash tests of twin-engine light aircraft. The major goals of the analytical portion of the program are outlined, including the development and validation of the analytical technique using simplified structural specimens that approximate aircraft components, as well as the mathematical modeling of the complete airframe and its subsequent dynamic analysis by substructuring and matrix reduction techniques.

Alfaro-Bou, E.↗
NASA NTRS1977

Effects of mechanical and thermal cycling on composite and hybrid laminates with residual stresses

The effects of tensile load cycling and thermal cycling on residual stiffness and strength properties of the following composite and hybrid angle-ply laminates were studied: boron/epoxy, boron/polyimide, graphite/low-modulus epoxy, graphite/high-modulus epoxy, graphite/polyimide, S-glass/epoxy, graphite/Kevlar 49/epoxy, and graphite/S-glass/epoxy. Specimens of the first six types were mechanically cycled up to 90% of static strength. Those that survived 10 million cycles were tested statically to failure, and no significant changes in residual strength and modulus were noted. Specimens of all types were subjected to thermal cycling between room temperature and 411 K for the epoxy-matrix composites and 533 K for the polyimide-matrix composites. The residual strength and stiffness remained largely unchanged, except for the graphite/low-modulus epoxy, which showed reductions in both of approximately 35%. When low-temperature thermal cycling under tensile load was applied, there was a noticeable reduction in modulus and strength in the graphite/low-modulus epoxy and some strength reduction in the S-glass/epoxy.

Daniel, I. M.↗
NASA NTRS1984

On the stability of the decelerating laminar boundary layer

The stability of a decelerating boundary-layer flow is investigated experimentally and numerically. Experimentally, a flat plate having a Blasius boundary layer is decelerated in an 18 m towing tank. The boundary layer becomes unstable to two-dimensional waves, which break down into three-dimensional patterns, hairpin vortices, and finally turbulent bursts when the vortices lift off the wall. The unsteady boundary-layer equations are solved numerically to generate instantaneous velocity profiles for a range of boundary and initial conditions. A quasi-steady approximation is invoked and the stability of local velocity profiles is determined by solving the Orr-Sommerfeld equation using Chebyshev matrix methods. Comparisons are made between the numerical predictions and the experimentally observed instabilities.

Gad-El-hak, M.↗
NASA NTRS1982

Scout trajectory error propagation computer program

Since 1969, flight experience has been used as the basis for predicting Scout orbital accuracy. The data used for calculating the accuracy consists of errors in the trajectory parameters (altitude, velocity, etc.) at stage burnout as observed on Scout flights. Approximately 50 sets of errors are used in Monte Carlo analysis to generate error statistics in the trajectory parameters. A covariance matrix is formed which may be propagated in time. The mechanization of this process resulted in computer program Scout Trajectory Error Propagation (STEP) and is described herein. Computer program STEP may be used in conjunction with the Statistical Orbital Analysis Routine to generate accuracy in the orbit parameters (apogee, perigee, inclination, etc.) based upon flight experience.

Myler, T. R.↗
NASA NTRS1988

PIA update: Correlation analyses of mass spectra

The PIA instrument aboard the Giotto spacecraft (a time of flight spectrometer) has been presented elsewhere. The mass spectra used in this analysis were decoded and mass numbers assigned according to the presence of carbon and silver, using the global values for these elements in their spectral absence. The results presented here were obtained using a frequency of occurrence based on analysis which correlated how often mass numbers appear in the mass spectra and which mass numbers tend to occur together in the same spectra; no amplitude information is utilized. The data are presented as plots of mass vs coincident mass for different subsets of the PIA data set, with both axes having units of atomic mass. Frequency contours are plotted at approximately five percent contour intervals, relative to the maximum AMU occurrence in that plot. The plots presented are symmetrical about the matrix diagonal, i.e., every mass is coincident with itself in a given spectra.

Mason, L. W.↗
NASA NTRS1990

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗
NASA NTRS1989

Low-order design and high-order simulation of active closed-loop control for aerospace structures under construction

Partially constructed/assembled structures in space are complicated enough but their dynamics will also be operating in closed-loop with feedback controllers. The dynamics of such structures are modeled by large-scale finite element models. The model dimension L is extremely large (approximately 10,000) while the numbers of actuators (M) and sensors (P) are small. The model parameters M(sub m) mass matrix, D(sub o) damping matrix, and K(sub o) stiffness matrix, are all symmetric and sparse (banded). Thus simulation of open-loop structure models of very large dimension can be accomplished by special integration techniques for sparse matrices. The problem of simulation of closed-loop control of such structures is complicated by the addition of controllers. Simulation of closed-loop controlled structures is an essential part of the controller design and evaluation process. Current research in the following areas is presented: high-order simulation of actively controlled aerospace structures; low-order controller design and SCI compensation for unmodeled dynamics; prediction of closed-loop stability using asymptotic eigenvalue series; and flexible robot manipulator control experiment.

Balas, Mark J.↗
NASA NTRS1995

Krylov methods for compressible flows

We investigate the application of Krylov methods to compressible flows, and the effect of implicit boundary conditions on the implicit solution of nonlinear problems. Two defect-correction procedures, namely, approximate factorization (AF) for structured grids and ILU/GMRES for general grids, are considered. Also considered here are Newton-Krylov matrix-free methods that we combined with the use of mixed discretization schemes in the implicitly defined Jacobian and its preconditioner. Numerical experiments that show the performance of our approaches are then presented.

Tidriri, M. D.↗