Search NASA⌕ Search

SEARCH · Search NASA

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 37 records · Page 2

Computing Challenges in Coded Mask Imaging

This slide presaentation reviews the complications and challenges in developing computer systems for Coded Mask Imaging telescopes. The coded mask technique is used when there is no other way to create the telescope, (i.e., when there are wide fields of view, high energies for focusing or low energies for the Compton/Tracker Techniques and very good angular resolution.) The coded mask telescope is described, and the mask is reviewed. The coded Masks for the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) instruments are shown, and a chart showing the types of position sensitive detectors used for the coded mask telescopes is also reviewed. Slides describe the mechanism of recovering an image from the masked pattern. The correlation with the mask pattern is described. The Matrix approach is reviewed, and other approaches to image reconstruction are described. Included in the presentation is a review of the Energetic X-ray Imaging Survey Telescope (EXIST) / High Energy Telescope (HET), with information about the mission, the operation of the telescope, comparison of the EXIST/HET with the SWIFT/BAT and details of the design of the EXIST/HET.

Skinner, Gerald↗

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↗

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.↗

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.↗

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.↗

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.↗

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↗

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.↗

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.↗

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↗

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↗

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↗

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.↗

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↗

The Effect of Spaceflight on Bone Cell Cultures

Understanding the response of bone to mechanical loading (unloading) is extremely important in defining the means of adaptation of the body to a variety of environmental conditions such as during heightened physical activity or in extended explorations of space or the sea floor. The mechanisms of the adaptive response of bone are not well defined, but undoubtedly they involve changes occurring at the cellular level of bone structure. This proposal has intended to examine the hypothesis that the loading (unloading) response of bone is mediated by specific cells through modifications of their activity cytoskeletal elements, and/or elaboration of their extracellular matrices. For this purpose, this laboratory has utilized the results of a number of previous studies defining molecular biological, biochemical, morphological, and ultrastructural events of the reproducible mineralization of a primary bone cell (osteoblast) culture system under normal loading (1G gravity level). These data and the culture system then were examined following the use of the cultures in two NASA shuttle flights, STS-59 and STS-63. The cells collected from each of the flights were compared to respective synchronous ground (1G) control cells examined as the flight samples were simultaneously analyzed and to other control cells maintained at 1G until the time of shuttle launch, at which point they were terminated and studied (defined as basal cells). Each of the cell cultures was assayed in terms of metabolic markers- gene expression; synthesis and secretion of collagen and non-collagenous proteins, including certain cytoskeletal components; assembly of collagen into macrostructural arrays- formation of mineral; and interaction of collagen and mineral crystals during calcification of the cultures. The work has utilized a combination of biochemical techniques (radiolabeling, electrophoresis, fluorography, Western and Northern Blotting, and light microscopic immunofluorescence) and structural methods (conventional and high voltage electron microscopy, inununocytochemistry, stereomicroscopy, and 3D image reconstruction). The studies have provided new knowledge of aspects of bone cell development and structural regulation, extracellular matrix assembly, and mineralization during spaceflight and under normal gravity. The information has contributed to insights into the means in general by which cells respond and adapt to different conditions of gravity (loading). The data may as well have suggested an underlying basis for the observed loss of bone by vertebrates, including man, in microgravity; and these scientific results may have implications for understanding bone loss following fracture healing and extended periods of inactivity such as during long-term bedrest.

Landis, William J.↗

Implementation of a digital optical matrix-vector multiplier using a holographic look-up table and residue arithmetic

The design and implementation of a digital (numerical) optical matrix-vector multiplier are presented. The objective is to demonstrate the operation of an optical processor designed to minimize computation time in performing a practical computing application. This is done by using the large array of processing elements in a Hughes liquid crystal light valve, and relying on the residue arithmetic representation, a holographic optical memory, and position coded optical look-up tables. In the design, all operations are performed in effectively one light valve response time regardless of matrix size. The features of the design allowing fast computation include the residue arithmetic representation, the mapping approach to computation, and the holographic memory. In addition, other features of the work include a practical light valve configuration for efficient polarization control, a model for recording multiple exposures in silver halides with equal reconstruction efficiency, and using light from an optical fiber for a reference beam source in constructing the hologram. The design can be extended to implement larger matrix arrays without increasing computation time.

Habiby, Sarry F.↗

ATIC Experiment: Preliminary Results from the Flight in 2002

Abstract The Advanced Thin Ionization Calorimeter (ATIC) had successful Long Duration Balloon flights from McMurdo, Antarctica in both 2000 and 2002. The instrument consists of a Silicon matrix for charge measurement, a flared graphite target to induce nuclear interactions, scintillator strip hodoscopes for triggering and helping reconstruct trajectory, and a BGO calorimeter to measure the energy of incident particles. In this paper, we discuss the second flight, which lasted 20 days, starting on 12/29/02. Preliminary results from the on-going analysis of the data including the proton and helium spectra are reported.

Ahn, H. S.↗

The ATIC Long Duration Balloon Project

Long Duration Balloon (LDB) scientific experiments, launched to circumnavigate the south pole over Antarctica, have particular advantages compared to Shuttle or other Low Earth Orbit (LEO) missions in terms of cost, weight, scientific 'duty factor' and work force development. The Advanced Thin Ionization Calorimeter (ATIC) cosmic ray astrophysics experiment is a good example of a university-based project that takes full advantage of current LDB capability. The ATIC experiment is currently being prepared for its first LDB science flight that will investigate the charge composition and energy spectra of primary cosmic rays over the energy range from about 10(exp 10) to 10(exp 14) eV. The instrument is built around a fully active, Bismuth Germanate (BGO) ionization calorimeter to measure the energy deposited by the cascades formed by particles interacting in a thick carbon target. A highly segmented silicon matrix, located above the target, provides good incident charge resolution plus rejection of the 'backscattered' particles from the interaction. Trajectory reconstruction is based on the cascade profile in the BGO calorimeter, plus information from the three pairs of scintillator hodoscope layers in the target section above it. A full evaluation of the experiment was performed during a test flight occurring between 28 December 2000 and 13 January 2001 where ATIC was carried to an altitude of approx. 37 km above Antarctica by an approx. 850,000 cu m helium filled balloon for one circumnavigation of the continent. All systems behaved well, the detectors performed as expected, more than 43 gigabytes of engineering and cosmic ray event data was returned and these data are now undergoing preliminary data analysis. During the coming 2002-2003 Antarctica summer season, we are preparing for a ATIC science flight with approx. 15 to 30 days of continuous data collection in the near-space environment of LDB float altitudes.

Guzik, T. G.↗