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 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Massively parallel quantum chemical density matrix renormalization group method

Here, we present, to the best of our knowledge, the first attempt to exploit the supercomputer platform for quantum chemical density matrix renormalization group (QCDMRG) calculations. We have developed the parallel scheme based on the in-house MPI global memory library, which combines operator and symmetry sector parallelisms, and tested its performance on three different molecules, all typical candidates for QC-DMRG calculations. In case of the largest calculation, which is the nitrogenase FeMo cofactor cluster with the active space comprising 113 electrons in 76 orbitals and bond dimension equal to 6000, our parallel approach scales up to approximately 2000 CPU cores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1981

Absorption and resonance Raman spectra of Pb2, Pb3, and Pb4 in xenon matrices

Matrix isolation techniques are used to investigate the spectra of lead molecules and, in particular, to obtain resonance Raman spectra of lead vapors isolated in solid xenon matrices. The presence of Pb2 is confirmed by the visible adsorption, and Raman spectra yield a vibrational frequency for the ground state of 108 per cm and a dissociation energy of 8200 per cm. A second resonance Raman progression indicates a Pb3 species of D3h symmetry. Finally, two additional Raman features at approximately 111 per cm spacing are evidence for a third species, tentatively identified as Pb4.

Stranz, D. D.↗

NASA NTRS2014

Effects of Nonsphericity on the Behavior of Lorenz-Mie Resonances in Scattering Characteristics of Liquid-Cloud Droplets

By using the results of highly accurate T-matrix computations for randomly oriented oblate and prolate spheroids and Chebyshev particles with varying degrees of asphericity, we analyze the effects of a deviation of water-droplet shapes from that of a perfect sphere on the behavior of Lorenz-Mie morphology-dependent resonances of various widths. We demonstrate that the positions and profiles of the resonances can change significantly with increasing asphericity. The absolute degree of asphericity required to suppress a Lorenz-Mie resonance is approximately proportional to the resonance width. Our results imply that numerical averaging of scattering characteristics of real cloud droplets over sizes may rely on a significantly coarser size-parameter resolution than that required for ideal, perfectly spherical particles.

matrices (mathematics)↗

NASA NTRS2017

A New Method of Absorption-Phase Nanotomography for 3D Observation of Mineral-Organic-Water Textiles and its Application to Pristine Carbonaceous Chondrites

Pristine carbonaceous chondrites contain fine-grained matrix, which is composed largely of amorphous silicates, sub-micron silicate and sulfide crystals, and organic materials. They are regarded as primitive dust in the early Solar System that have suffered minimal alteration in their parent bodies. The matrix generally has different lithologies; some of them are unaltered but some are more or less aqueously altered. Their textures have been examined in 2D usually by FE-SEM/EDS, TEM/EDS, nano-SIMS and micro-XRD. Observation of their complex fine textures, such as spatial relation between different lithologies in 3D, is important for understanding aggregation and alteration processes. Synchrotron radiation (SR)-based X-ray tomography reveals 3D structures nondestructively with high spatial resolution of approximately greater than 100 nm. We have developed a new technique using absorption contrasts called "dual-energy tomography" (DET) to obtain 3D distribution of minerals at SPring-8, SR facility in Japan, and applied successfully to Itokawa particles. Phase and absorption contrast images can be simultaneously obtained in 3D by using "scanning-imaging x-ray microscopy" (SIXM) at SPring-8, which can discriminate between void, water and organic materials. We applied this technique combined with FIB micro-sampling to carbonaceous chondrites to search for primitive liquid water. In this study, we combined the DET and SIXM to obtain three dimensional submicron-scale association between minerals, organic materials and water and applied this to pristine carbonaceous chondrites.

Tsuchiyama, A.↗

DOE OSTI2024

Denoising of imaginary time response functions with Hankel projections

Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost—in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies. Published by the American Physical Society 2024

Yu, Yang (ORCID:0000000186178878)↗

DOE OSTI2024

On a Spectral Method for β -particle Bound Excitation Collisions in Kilonovae

The interaction of β-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova (KN) transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti following the seminal formulation of Bethe, in a spectral solver of the Vlasov–Maxwell–Boltzmann equation. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon–matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric KN ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for β-particles in the KN is a bulk velocity parameter u = 0, a thermal velocity parameter α = 0.5c, and a 9 × 9 × 9 mode velocity basis set (Hermite orders of 0–8 in each dimension). For interior-ejecta sample zones, we estimate that the ratio of thermalization from large-angle (≳2fdg5) bound excitation scattering to total thermalization is ~0.002–0.003.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2011

Spatially Resolved, In Situ Carbon Isotope Analysis of Archean Organic Matter

Spatiotemporal variability in the carbon isotope composition of sedimentary organic matter (OM) preserves information about the evolution of the biosphere and of the exogenic carbon cycle as a whole. Primary compositions, and imprints of the post-depositional processes that obscure them, exist at the scale of individual sedimentary grains (mm to micron). Secondary ion mass spectrometry (SIMS) (1) enables analysis at these scales and in petrographic context, (2) permits morphological and compositional characterization of the analyte and associated minerals prior to isotopic analysis, and (3) reveals patterns of variability homogenized by bulk techniques. Here we present new methods for in situ organic carbon isotope analysis with sub-permil precision and spatial resolution to 1 micron using SIMS, as well as new data acquired from a suite of Archean rocks. Three analytical protocols were developed for the CAMECA ims1280 at WiscSIMS to analyze domains of varying size and carbon concentration. Average reproducibility (at 2SD) using a 6 micron spot size with two Faraday cup detectors was 0.4 %, and 0.8 % for analyses using 1 micron and 3 micron spot sizes with a Faraday cup (for C-12) and an electron multiplier (for C-13). Eight coals, two ambers, a shungite, and a graphite were evaluated for micron-scale isotopic heterogeneity, and LCNN anthracite (delta C-13 = -23.56 +/- 0.1 %, 2SD) was chosen as the working standard. Correlation between instrumental bias and H/C was observed and calibrated for each analytical session using organic materials with H/C between 0.1 and 1.5 (atomic), allowing a correction based upon a C-13H/C-13 measurement included in every analysis. Matrix effects of variable C/SiO2 were evaluated by measuring mm to sub-micron graphite domains in quartzite from Bogala mine, Sri Lanka. Apparent instrumental bias and C-12 count rate are correlated in this case, but this may be related to a crystal orientation effect in graphite. Analyses of amorphous Archean OM suggest that instrumental bias is consistent for 12C count rates as low as 10% relative to anthracite. Samples from the ABDP-9 (n=3; Mount McRae Shale, approximately 2.5 Ga), RHDH2a (n=2; Carrawine Dolomite and Jeerinah Fm, approximately 2.6 Ga), WRL1 (n=3; Wittenoom Fm, Marra Mamba Iron Formation, and Jeerinah Fm, approximately 2.6 Ga), and SV1 (n=1; Tumbiana Fm, approximately 2.7 Ga) drill cores, each previously analyzed for bulk organic carbon isotope composition, yielded 100 new, in situ data from Neoarchean sedimentary OM. In these samples, delta C-13 varies between -53.1 and -28.3 % and offsets between in situ and bulk compositions range from -8.3 to 18.8%. In some cases, isotopic composition and mode of occurrence (e.g. morphology and mineral associations) are statistically correlated, enabling the identification of distinct reservoirs of OM. Our results support previous evidence for gradients of oxidation with depth in Neoarchean environments driven by photosynthesis and methane metabolism. The relevance of these findings to questions of bio- and syngenicity as well as the alteration history of previously reported Archean OM will be discussed.

Williford, Kenneth H.↗

NASA NTRS2017

Assessment of the Gaussian Covariance Approximation over an Earth-Asteroid Encounter Period

In assessing the risk an asteroid may pose to the Earth, the asteroids state is often predicted for many years, often decades. Only by accounting for the asteroids initial state uncertainty can a measure of the risk be calculated. With the asteroids state uncertainty growing as a function of the initial velocity uncertainty, orbit velocity at the last state update, and the time from the last update to the epoch of interest, the asteroids position uncertainties can grow to many times the size of the Earth when propagated to the encounter risk corridor. This paper examines the merits of propagating the asteroids state covariance as an analytical matrix. The results of this study help to bound the efficacy of applying different metrics for assessing the risk an asteroid poses to the Earth. Additionally, this work identifies a criterion for when different covariance propagation methods are needed to continue predictions after an Earth-encounter period.

planetary defense↗

NASA NTRS2017

Assessment of the Gaussian Covariance Approximation over an Earth-Asteroid Encounter Period

uncertainty characterization↗

DOE OSTI2022

Developments in Performance and Portability for MadGraph5_aMC@NLO

Event generators simulate particle interactions using Monte Carlo techniques, providing the primary connection between experiment and theory in experimental high energy physics. These software packages, which are the first step in the simulation worflow of collider experiments, represent approximately 5 to 20% of the annual WLCG usage for the ATLAS and CMS experiments. With computing architectures becoming more heterogeneous, it is important to ensure that these key software frameworks can be run on future systems, large and small. In this contribution, recent progress on porting and speeding up the Madgraph5_aMC@NLO event generator on hybrid architectures, i.e. CPU with GPU accelerators, is discussed. The main focus of this work has been in the calculation of scattering amplitudes and "matrix elements", which is the computational bottleneck of an event generation application. For physics processes limited to QCD leading order, the code generation toolkit has been expanded to produce matrix element calculations using C++ vector instructions on CPUs and using CUDA for NVidia GPUs, as well as using Alpaka, Kokkos and SYCL for multiple CPU and GPU architectures. Performance is reported in terms of matrix element calculations per time on NVidia, Intel, and AMD devices. The status and outlook for the integration of this work into a production release usable by the LHC experiments, with the same functionalities and very similar user interfaces as the current Fortran version, is also described.

Valassi, Andrea↗

NASA NTRS1973

Stress wave calculations in composite plates using the fast Fourier transform.

The protection of composite turbine fan blades against impact forces has prompted the study of dynamic stresses in composites due to transient loads. The mathematical model treats the laminated plate as an equivalent anisotropic material. The use of Mindlin's approximate theory of crystal plates results in five two-dimensional stress waves. Three of the waves are flexural and two involve in-plane extensional strains. The initial value problem due to a transient distributed transverse force on the plate is solved using Laplace and Fourier transforms. A fast computer program for inverting the two-dimensional Fourier transform is used. Stress contours for various stresses and times after application of load are obtained for a graphite fiber-epoxy matrix composite plate. Results indicate that the points of maximum stress travel along the fiber directions.

Moon, F. C.↗

NASA NTRS1977

Implicit approximate-factorization schemes for the efficient solution of steady transonic flow problems

Implicit approximate-factorization algorithms (AF) are developed for the solution of steady-state transonic flow problems. The performance of the AF solution method is evaluated relative to that of the standard solution method for transonic flow problems, successive line over-relaxation (SLOR). Both methods are applied to the solution of the nonlinear, two-dimensional transonic small-disturbance equation. Results indicate that the AF method requires substantially less computer time than SLOR to solve the nonlinear finite-difference matrix equation for a transonic flow field. This increase in computational efficiency is achieved with no appreciable increase in computer storage or coding complexity.

Ballhaus, W. F.↗

NASA NTRS1992

LDEF materials: An overview of the interim findings

The flight and retrieval of the National Aeronautics and Space Administration's Long Duration Exposure Facility (LDEF) provided an opportunity for the study of the low-Earth orbit (LEO) environment and long-duration space environmental effect (SEE) on materials that are unparalleled in the history of the U.S. space program. The 5.8-year flight of LDEF greatly enhanced the potential value of materials data from LDEF to the international SEE community, compared to that of the original 1-year flight plan. The remarkable flight attitude stability of LDEF enables specific analyses of various individual and combined effects of LEO environmental parameters on identical materials of the same space vehicle. NASA recognized the potential by forming the LDEF Space Environmental Effects on Materials Special Investigation Group (MSIG) to address the greatly expanded materials and LEO space environment parameter analysis opportunities available in the LDEF structure, experiment trays, and corollary measurements, so that the combined value of all LDEF materials data to current and future space missions will be assessed and documented. This paper provides an overview of the interim LDEF materials findings of the Principal Investigators and the Materials Special Investigation Group. These revelations are based on observations of LEO environmental effects on materials made in-space during LDEF retrieval and during LDEF tray deintegration at the Kennedy Space Center, and on findings of approximately 1.5 years of laboratory analyses of LDEF materials by the LDEF materials scientists. These findings were extensively reviewed and discussed at the MSIG-sponsored LDEF Materials Workshop '91. The results are presented in a format which categorizes the revelations as 'clear findings' or 'confusing/unexplained findings' and resultant needs for new space materials developments and ground simulation testing/analytical modeling in seven categories: environmental parameters and data bases; LDEF contamination; thermal control coatings and protective treatments; polymers and films; polymer-matrix composites; metals, ceramics, and optical materials; and systems-related materials. General outlines of findings of the other LDEF Special Investigation Groups (Ionizing Radiation, Meteoroid and Debris, and Systems) are also included. The utilization of LDEF materials data for future low-earth orbit missions is also discussed, concentrating on Space Station Freedom. Some directions for continuing studies of LDEF materials are outlined. In general, the LDEF data is remarkable consistent; LDEF will provide a 'benchmark' for materials design data bases for satellites in low-Earth orbit. Some materials were identified to be encouragingly resistant to LEO SEE for 5.8-years; other 'space qualified' materials displayed significant environmental degradation. Molecular contamination was widespread; LDEF offers an unprecedented opportunity to provide a unified perspective of unmanned LEO spacecraft contamination mechanisms. New material development requirements for long-term LEO missions have been identified and current ground simulation testing methods/data for new, durable materials concepts can be validated with LDEF results. LDEF findings are already being integrated into the design of Space Station Freedom.

Stein, Bland A.↗

NASA NTRS1993

Accuracy of the Generalized Self-Consistent Method in Modelling the Elastic Behaviour of Periodic Composites

Local stress and strain fields in the unit cell of an infinite, two-dimensional, periodic fibrous lattice have been determined by an integral equation approach. The effect of the fibres is assimilated to an infinite two-dimensional array of fictitious body forces in the matrix constituent phase of the unit cell. By subtracting a volume averaged strain polarization term from the integral equation we effectively embed a finite number of unit cells in a homogenized medium in which the overall stress and strain correspond to the volume averaged stress and strain of the constrained unit cell. This paper demonstrates that the zeroth term in the governing integral equation expansion, which embeds one unit cell in the homogenized medium, corresponds to the generalized self-consistent approximation. By comparing the zeroth term approximation with higher order approximations to the integral equation summation, both the accuracy of the generalized self-consistent composite model and the rate of convergence of the integral summation can be assessed. Two example composites are studied. For a tungsten/copper elastic fibrous composite the generalized self-consistent model is shown to provide accurate, effective, elastic moduli and local field representations. The local elastic transverse stress field within the representative volume element of the generalized self-consistent method is shown to be in error by much larger amounts for a composite with periodically distributed voids, but homogenization leads to a cancelling of errors, and the effective transverse Young's modulus of the voided composite is shown to be in error by only 23% at a void volume fraction of 75%.

Walker, Kevin P.↗

NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

NASA NTRS2020

Ultrasonic Lamb Wave Inspection of Composite Defects

Composites have been extensively used in aerospace engineering due to their advantages of light weight, high strength, and engineering design flexibility. Manufacturing defects such as wrinkle and porosity can affect the performance of the composites and may lead to failure in the end, while damage such as delamination, fiber fraction, and matrix cracking can directly cause failure of the composites. In this paper, a Lamb wave based nonintrusive nondestructive evaluation system, which employs piezoelectric transducer for actuation and scanning laser Doppler vibrometer for wavefield sensing, is presented for typical composite defect and damage inspection and evaluation. Two composite panels with different geometry (flat or curved) and with various embedded defects (wrinkle and delamination) are inspected using the nonintrusive Lamb wave system. Both the wrinkles and delamination are detected from the wavefield and approximately quantified through wavefield imaging methods.

Wrinkles↗

DOE OSTI2025

Laser powder bed fusion of ODS Fe–Cr–Al (0.3Zr, 0.3Y 2 O 3 ): Unveiling processing-microstructure- mechanical property relationships

Here, this study investigates the fabrication of oxide dispersion strengthened (ODS) Fe-Cr-Al alloys via laser powder bed fusion (LPBF) with strategic additions of 0.3 wt% Zr and 0.3 wt% Y 2 O 3 for enhanced mechanical performance in nuclear applications. Systematic processing parameter optimization yielded three distinct conditions: one low-density product with significant defects and two near-full-density materials with improved consolidation. Comprehensive characterization confirmed single-phase α-ferrite matrix formation with successful incorporation of Y-, Zr-, O-, and C-rich precipitates characteristic of ODS alloys. However, precipitate density remained low (∼10 7 cm −3 ), resulting in sink strength values substantially below optimal levels for radiation resistance. Microhardness values (mid-200s HV) correlated inversely with grain size following the Hall-Petch relationship, indicating grain boundary strengthening as the dominant mechanism rather than precipitation strengthening. The optimized processing conditions achieved excellent mechanical properties with room temperature yield strength of approximately 500 MPa and 30 % elongation, demonstrating superior strength-ductility synergy compared to other additively manufactured ODS materials and performance consistent with literature values for LPBF-processed ODS-FeCrAl alloys. This investigation reveals both the potential and limitations of LPBF processing for ODS Fe-Cr-Al alloys. While successful defect-free fabrication was achieved, results highlight the critical need for systematic optimization of processing parameters and post-processing heat treatments to enhance precipitate density for effective dispersion strengthening and radiation resistance while maintaining additive manufacturing advantages.

Additive manufacturing↗

NASA NTRS1990

Effect of atmospherics on beamforming accuracy

Two mathematical representations of noise due to atmospheric turbulence are presented. These representations are derived and used in computer simulations of the Bartlett Estimate implementation of beamforming. Beamforming is an array processing technique employing an array of acoustic sensors used to determine the bearing of an acoustic source. Atmospheric wind conditions introduce noise into the beamformer output. Consequently, the accuracy of the process is degraded and the bearing of the acoustic source is falsely indicated or impossible to determine. The two representations of noise presented here are intended to quantify the effects of mean wind passing over the array of sensors and to correct for these effects. The first noise model is an idealized case. The effect of the mean wind is incorporated as a change in the propagation velocity of the acoustic wave. This yields an effective phase shift applied to each term of the spatial correlation matrix in the Bartlett Estimate. The resultant error caused by this model can be corrected in closed form in the beamforming algorithm. The second noise model acts to change the true direction of propagation at the beginning of the beamforming process. A closed form correction for this model is not available. Efforts to derive effective means to reduce the contributions of the noise have not been successful. In either case, the maximum error introduced by the wind is a beam shift of approximately three degrees. That is, the bearing of the acoustic source is indicated at a point a few degrees from the true bearing location. These effects are not quite as pronounced as those seen in experimental results. Sidelobes are false indications of acoustic sources in the beamformer output away from the true bearing angle. The sidelobes that are observed in experimental results are not caused by these noise models. The effects of mean wind passing over the sensor array as modeled here do not alter the beamformer output as significantly as expected.

Alexander, Richard M.↗