Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix reconstruction”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗
NASA NTRS1994

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗
DOE OSTI2026

Search for HH → bbτ⁺τ⁻ Using Run 3 Scouting Data Analyze b-tagging and tau-tagging Performance with Unified Particle Transformer

B-tagging and tau-tagging performances play an important role in the search for the rare event HH → bbτ⁺τ⁻. A transformer-based neural network, Unified Particle Transformer, is applied for both tagging tasks, and Run 3 proton–proton collision scouting data at center-of-mass energy of 13.6 TeV is used. The scouting data stream accepts events at a much higher rate compared to traditional triggers, but stores only the objects reconstructed in the trigger, no low-level detector information. Therefore, existing taggers trained for the offline event reconstruction cannot be used. Analysis of the SoftMax plots, ROC/AUC curves, confusion matrix, accuracy and losses are used to evaluate model performance. Specifically, the tagging efficiency of the signal and misidentification probability across multiple background processes are compared for varying working points. Different training samples with distinct distributions of jet flavors are utilized and related model performances are analyzed. Interpretability methods, such as Integrated Gradients, may further be applied to study the input features’ influence on the model’s decisions, providing insights into potential improvements.

Chen, Blair [Purdue U., West Lafayette; Fermilab]↗
NASA NTRS1982

Flying qualities criteria for GA single pilot IFR operations

The flying qualities criteria in general aviation (GA) to decrease accidents are discussed. The following in-flight research is discussed: (1) identification of key aerodynamic configurations; (2) implementation of an in-flight simulator; (3) mission matrix design; (4) experimental systems; (5) data reduction; (6) optimal flight path reconstruction. Some of the accomplished work is reported: an integrated flight testing and flight path reconstruction methodology was developd, high accuracy in trajectory estimation was achieved with an experimental setup, and a part of the flight test series was flown.

Bar-Gill, A.↗
NASA NTRS1981

Organic compounds in meteorites and their origins

The current investigation represents an extensively updated version of a review conducted by Anders et al. (1973). The investigation takes into account the literature through mid-1980. It is pointed out that Type 1 carbonaceous chondrites (C1) contain 6% of their cosmic complement of carbon, mainly in the form of organic matter. Most authors now agree that this material represents primitive prebiotic matter. The principal questions remaining are what abiotic processes formed the organic matter, and to what extent these processes took place in locales other than the solar nebula, such as interstellar clouds or meteorite parent bodes. The problem is approached in three stages. It is attempted to reconstruct the physical conditions during condensation from the clues contained in the inorganic matrix of the meteorite. The condensation behavior of carbon under these conditions is determined on the basis of thermodynamic calculations. Model experiments on the condensation of carbon are performed, and the synthesized compounds are compared with those actually found in meteorites.

Hayatsu, R.↗
NASA NTRS2000

Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absorption information. An algorithm, which combines a two dimensional (2D) matrix inversion with a one-dimensional (1D) Fourier transform inversion is used to obtain images of three dimensional hidden objects in turbid scattering media.

Alfano, Robert R.↗
NASA NTRS2018

An Investigation of Parallel Programming Techniques Applied to Monte Carlo Simulations for Post-Flight Reconstruction of Spacecraft Trajectory

Parallelizing software to execute on multi-core central processing units (CPUs) and graphics processing units (GPUs) can be challenging. For some fields outside of Computer Science, this transition comes with new issues. For example, memory limitations can require modifications to code not initially developed to run on GPUs. This work applies the Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) directive-based parallelization strategies on a Monte Carlo simulation approach for trajectory reconstruction enabling it to run on multi-core CPUs and GPUs. Large matrix operations are the most common use of GPUs, which are not present in this algorithm; however, the natural parallelism of independent trajectories in Monte Carlo simulations is exploited. Benchmarking data are presented comparing execution times of the software for single-thread CPUs, multi-thread CPUs with OpenMP, and multi-thread GPUs using OpenACC. These data were collected using nodes with Intel® Xeon® E5-2670 (Sandy Bridge) CPUs enhanced with NVIDIA® Tesla® K40 GPUs on the Pleiades Supercomputer cluster at the National Aeronautics and Space Administration (NASA) Ames Research Center (ARC) and a local Intel® Xeon Phi™ node at NASA Langley Research Center (LaRC).

Williams, R. Anthony↗
DOE OSTI2025

Structural response reconstruction using a system-equivalent singular vector basis

Here, this paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.

42 ENGINEERING↗
DOE OSTI2025

Quasiprobabilistic Readout Correction of Midcircuit Measurements for Adaptive Feedback via Measurement Randomized Compiling

Quantum measurements are a fundamental component of quantum computing. However, on present-day quantum computers, measurements can be more error prone than quantum gates and are susceptible to nonunital errors as well as nonlocal correlations due to measurement crosstalk. While readout errors can be mitigated in postprocessing, this is inefficient in the number of qubits due to a combinatorially large number of possible states that need to be characterized. In this work, we show that measurement errors can be tailored into a simple stochastic error model using randomized compiling, enabling the efficient mitigation of readout errors via quasiprobability distributions reconstructed from the measurement of a single preparation state in an exponentially large confusion matrix. We demonstrate the scalability and power of this approach by correcting readout errors without matrix inversion on a large number of different preparation states applied to a register of eight superconducting transmon qubits. Moreover, we show that this method can be extended to midcircuit measurements used for active feedback via quasiprobabilistic error cancellation, and we demonstrate the correction of measurement errors on an ancilla qubit used to detect and actively correct bit-flip errors on an entangled memory qubit. Our approach enables the correction of readout errors on large numbers of qubits and offers a strategy for correcting readout errors in adaptive circuits in which the results of midcircuit measurements are used to perform conditional operations on nonlocal qubits in real time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Organic compounds in meteorites

The problem of whether organic compounds originated in meteorites as a primary condensate from a solar gas or whether they were introduced as a secondary product into the meteorite during its residence in a parent body is examined by initially attempting to reconstruct the physical conditions during condensation (temperature, pressure, time) from clues in the inorganic matrix of the meteorite. The condensation behavior of carbon under these conditions is then analyzed on the basis of thermodynamic calculations, and compounds synthesized in model experiments on the condensation of carbon are compared with those actually found in meteorites. Organic compounds in meteorites seem to have formed by catalytic reactions of carbon monoxide, hydrogen, and ammonia in the solar nebula at 360 to 400 K temperature and about 3 to 7.6 microtorr pressure. The onset of these reactions was triggered by the formation of suitable catalysts (magnetite, hydrated silicates) at these temperatures.

Anders, E.↗
NASA NTRS1999

Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: ##EQU1## wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absoption information. An algorithm, which combines a two dimensional (2D) matrix inversion with a one-dimensional (1D) Fourier transform inversion is used to obtain images of three dimensional hidden objects in turbid scattering media.

Alfano, Robert R.↗
NASA NTRSDate not supplied

Multi-Scale Thermo-Mechanical Modeling of Porous 3D Woven TPS Materials

This work summarizes the process to compute and analyze the thermal conductivity and mechanical properties of 3D woven TPS materials such as 3MDCP (3D Mid-Density Carbon-Phenolic). The PuMA [1,2] software, developed at NASA Ames, was used to characterize different 3MDCP samples from their constituents' data, averaging their thermal conductivity and elastic properties in the three main directions to obtain their effective orthotropic thermo-mechanical properties. This was performed in multiple steps: firstly, TPS samples were digitally reconstructed using micro computed tomography (µCT) and their constituents were segmented; then, the porous matrix phase was analyzed at the micro-scale and these results were used, along with the fibers’ constituents information, to model the tows at the meso-scale; finally, results from the constituents were homogenized and used to model the thermo-mechanical behavior of the unit cell at the macro-scale. Additionally, to gain a better understanding of the tows’ morphology and distribution of the carbon-phenolic blended fibers, microscopy images of the tows’ cross-section were segmented using deep learning techniques and analyzed with PuMA. This workflow can be applied to any TPS material to computationally obtain its thermo-mechanical properties, enabling more informed TPS design and manufacturing choices.

Multi-Scale↗
NASA NTRSDate not supplied

Multi-Scale Thermo-Mechanical Modeling of Porous 3D Woven TPS Materials

This work summarizes the process to compute and analyze the thermal conductivity and mechanical properties of 3D woven TPS materials such as 3MDCP (3D Mid-Density Carbon-Phenolic). The PuMA software, developed at NASA Ames, was used to characterize different 3MDCP samples from their constituents' data, averaging their thermal conductivity and elastic properties in the three main directions to obtain their effective orthotropic thermo-mechanical properties. This was performed in multiple steps: firstly, TPS samples were digitally reconstructed using micro computed tomography (µCT) and their constituents were segmented; then, the porous matrix phase was analyzed at the micro-scale and these results were used, along with the fibers’ constituents information, to model the tows at the meso-scale; finally, results from the constituents were homogenized and used to model the thermo-mechanical behavior of the unit cell at the macro-scale. Additionally, to gain a better understanding of the tows’ morphology and distribution of the carbon-phenolic blended fibers, microscopy images of the tows’ cross-section were segmented using deep learning techniques and analyzed with PuMA. This workflow can be applied to any TPS material to computationally obtain its thermo-mechanical properties, enabling more informed TPS design and manufacturing choices.

Multi-Scale↗
NASA NTRS2024

Constraining the Thermal and Textural Evolution of Evolved Achondrite DOM 10100: Implications for Volatile Loss

Identification of evolved achrondritic clasts suggests that some early planetesimals experienced extensive differentiation and formed chemically evolved crusts. Many of these samples contain volatile-bearing phases (i.e., apatite, merrillite, nominally anhydrous minerals), which can be used to constrain volatile budgets and abundancies via in-situ measurements. These ancient crustal rocks are subject to high heat, moderate pressures, potentially hot mineral-rich fluids, and degassing that has likely modified many of these volatiles, and their abundancies do not represent the composition of the primary magma from which these achondrites formed. In this work we attempt to deconvolve the thermal evolution of these samples to provide petrologic information that can be used for making inferences regarding their volatile history, and ultimately constraining volatile budgets on planetary bodies in the early Solar System. Here, we combine in-situ major element geochemistry, thermodynamic modeling, and detailed textural and microstructural observations to reconstruct the thermal and textural evolution of howardite DOM 10100, which contains a dacite within its matrix. Previous analyses of several sections of DOM 10100 suggest that the dacitic clast is evidence for partial melting of Vesta’s differentiated basaltic crust, making it an ideal candidate for evolved crustal studies. While its igneous petrogenesis has been well constrained, there is still opportunity to better characterize its post-crystallization thermal history, including possible thermal metamorphism and cooling.

J S Gorce↗
NASA NTRS2002

Charge Resolution of the Silicon Matrix of the ATIC Experiment

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CRT particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. The silicon matrix of the ATIC spectrometer is designed to resolve individual elements from proton to iron. To provide this resolution careful calibration of each pixel of the silicon matrix is required. Firstly, for each electronic channel of the matrix the pedestal value was subtracted taking into account its drift during the flight. The muon calibration made before the flight was used then to convert electric signals (in ADC channel number) to energy deposits in each pixel. However, the preflight muon calibration was not accurate enough for the purpose, because of lack of statistics in each pixel. To improve charge resolution the correction was done for the position of Helium peak in each pixel during the flight . The other way to set electric signals in electronics channels of the Si-matrix to one scale was correction for electric channel gains accurately measured in laboratory. In these measurements it was found that small different nonlinearities for different channels are present in the region of charge Z > 20. The correction for these non-linearities was not done yet. In linear approximation the method provides practically the same resolution as muon calibration plus He-peak correction. For searching a pixel with the signal of primary particle an indication from the cascade in the calorimeter was used. For this purpose a trajectory was reconstructed using weight centers of energy deposits in BGO layers. The point of intersection of this trajectory with Si-matrix and its RMS was determined. The pixel with maximal signal in 3sigma region was taken as sought. The signal in this pixel was corrected by trajectory zenith angle. The preliminary results on charge resolution of the Si-matrix in the range from protons to iron are presented.

Zatsepin, V. I.↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
DOE OSTI2026

Edge Reconstruction in a Quantum Spin Hall Insulator

We study interaction-driven edge reconstruction in a quantum spin Hall insulator described by the Bernevig-Hughes-Zhang model with Kanamori-Hubbard interactions using the real-space density matrix renormalization group method in both the grand-canonical and canonical ensembles. For a two-dimensional cylinder with a smooth edge, we identify discrete particle-number transitions that lead to a spin-polarized edge state stabilized by an emergent ferromagnetic exchange interaction. The reconstruction is orbital-selective, occurring predominantly in the 𝑠-orbital channel. Our results reveal a microscopic mechanism for emergent fluctuating moments at the edge that could compromise the topological protection of helical edge states by time-reversal symmetry.

Soni, Rahul [ORNL] (ORCID:0000000317714299)↗
NASA NTRS2003

Structural Anomalies Detected in Ceramic Matrix Composites Using Combined Nondestructive Evaluation and Finite Element Analysis (NDE and FEA)

Most reverse engineering approaches involve imaging or digitizing an object and then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. The rapid prototyping technique builds high-quality physical prototypes directly from computer-aided design files. This fundamental technique for interpreting and interacting with large data sets is being used here via Velocity2 (an integrated image-processing software, ref. 1) using computed tomography (CT) data to produce a prototype three-dimensional test specimen model for analyses. A study at the NASA Glenn Research Center proposes to use these capabilities to conduct a combined nondestructive evaluation (NDE) and finite element analysis (FEA) to screen pretest and posttest structural anomalies in structural components. A tensile specimen made of silicon nitrite (Si3N4) ceramic matrix composite was considered to evaluate structural durability and deformity. Ceramic matrix composites are being sought as candidate materials to replace nickel-base superalloys for turbine engine applications. They have the unique characteristics of being able to withstand higher operating temperatures and harsh combustion environments. In addition, their low densities relative to metals help reduce component mass (ref. 2). Detailed three-dimensional volume rendering of the tensile test specimen was successfully carried out with Velocity2 (ref. 1) using two-dimensional images that were generated via computed tomography. Subsequent, three-dimensional finite element analyses were performed, and the results obtained were compared with those predicted by NDE-based calculations and experimental tests. It was shown that Velocity2 software can be used to render a three-dimensional object from a series of CT scan images with a minimum level of complexity. The analytical results (ref. 3) show that the high-stress regions correlated well with the damage sites identified by the CT scans and the experimental data. Furthermore, modeling of the voids collected via NDE offered an analytical advantage that resulted in more accurate assessments of the material s structural strength. The top figure shows a CT scan image of the specimen test section illustrating various hidden structural entities in the material and an optical image of the test specimen considered in this study. The bottom figure represents the stress response predicted from the finite element analyses (ref .3 ) for a selected CT slice where it clearly illustrates the correspondence of the high stress risers due to voids in the material with those predicted by the NDE. This study is continuing, and efforts are concentrated on improving the modeling capabilities to imitate the structural anomalies as detected.

Abdul-Aziz, Ali↗