Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix decomposition”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Dynamics of Molecular Hydrogen in Hypersaline Microbial Mars

Early Earth microbial communities that centered around the anaerobic decomposition of organic molecular hydrogen as a carrier of electrons, regulator of energy metabolism, and facilitator of syntroph'c microbial interactions. The advent of oxygenic photosynthetic organisms added a highly dynamic and potentially dominant term to the hydrogen economy of these communities. We have examined the daily variations of hydrogen concentrations in cyanobacteria-dominated microbial mats from hypersaline ponds in Baja California Sur, Mexico. These mats bring together phototrophic and anaerobic bacteria (along with virtually all other trophic groups) in a spatially ordered and chemically dynamic matrix that provides a good analog for early Earth microbial ecosystems. Hydrogen concentrations in the photic zone of the mat can be three orders of magnitude or more higher than in the photic zone, which are, in turn, an order of magnitude higher than in the unconsolidated sediments underlying the mat community. Within the photic zone, hydrogen concentrations can fluctuate dramatically during the diel (24 hour day-night) cycle, ranging from less than 0.001% during the day to nearly 10% at night. The resultant nighttime flux of hydrogen from the mat to the environment was up to 17% of the daytime oxygen flux. The daily pattern observed is highly dependent on cyanobacterial species composition within the mat, with Lyngbya-dominated systems having a much greater dynamic range than those dominated by Microcoleus; this may relate largely to differing degrees of nitrogen-fixing and fermentative activity in the two mats. The greatest H2 concentrations and fluxes were observed in the absence of oxygen, suggesting an important potential feedback control in the context of the evolution of atmospheric composition. The impact of adding this highly dynamic photosynthetic term to the hydrogen economy of early microbial ecosystems must have been substantial. From an evolutionary standpoint, the H2 generated in mats could have represented a very important new source of electrons and energy - but one that could not be harnessed without substantial adaptation to the highly variable chemistry of the mat surface. In addition, the emergent chemistry of anaerobic communities is often highly dependent on ambient hydrogen concentrations, so that incorporation of these communities into photosynthetic mats could have significantly affected the composition and flux of reduced "biosignature' gases to the environment.

Hoehler, Tori M.↗
NASA NTRS2013

Conducting nanotubes or nanostructures based composites, method of making them and applications

An electromagnetic interference (EMI) shielding material includes a matrix of a dielectric or partially conducting polymer, such as foamed polystyrene, with carbon nanotubes or other nanostructures dispersed therein in sufficient concentration to make the material electrically conducting. The composite is formed by dispersing the nanotube material in a solvent in which the dielectric or partially conducting polymer is soluble and mixing the resulting suspension with the dielectric or partially conducting polymer. A foaming agent can be added to produce a lightweight foamed material. An organometallic compound can be added to enhance the conductivity further by decomposition into a metal phase.

Gupta, Mool C.↗
NASA NTRSDate not supplied

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

Time-resolved, high-speed self-aligned focusing schlieren images were acquired in the NASA Glenn Research Center Transonic Turbine Blade Cascade facility to help understand the aerodynamic behavior of high-pressure, thick trailing edge turbine blades. The trailing edge thickness of 9% of axial chord tested represents simulated ceramic matrix composite fabrication constraints, which was verified previously to possess a high-loss flow regime at high inlet turbulence conditions over a narrow range of Reynolds numbers and at a fixed design exit Mach number of 0.74. Our high-speed images, which were acquired at 10 distinct Reynolds numbers, show a significant increase in energy from flow oscillations due to transonic vortex shedding at Reynolds numbers corresponding to the high loss conditions. For those conditions, strong acoustic waves turn into shock waves. Spectral Proper Orthogonal Decomposition of the high-speed images shows acoustic waves from trailing edge vortex shedding at all conditions, with increased spectral energy at the high-loss conditions and slightly increasing frequency (about 6%) as a function of Reynolds number. Analysis of potential feedback timing is performed using velocity fields from a previous LES simulation, considering different feedback mechanisms. Most noteworthy is the acoustic/shock-boundary layer interaction mechanism on the suction surface at the blade geometric throat, which likely plays an important role in realistic curved blade passages.

Trailing Edge↗
NASA NTRS1973

Decomposition-aggregation stability analysis

This report presents the development and description of the decomposition aggregation approach to stability investigations of high dimension mathematical models of dynamic systems. The high dimension vector differential equation describing a large dynamic system is decomposed into a number of lower dimension vector differential equations which represent interconnected subsystems. Then a method is described by which the stability properties of each subsystem are aggregated into a single vector Liapunov function, representing the aggregate system model, consisting of subsystem Liapunov functions as components. A linear vector differential inequality is then formed in terms of the vector Liapunov function. The matrix of the model, which reflects the stability properties of the subsystems and the nature of their interconnections, is analyzed to conclude over-all system stability characteristics. The technique is applied in detail to investigate the stability characteristics of a dynamic model of a hypothetical spinning Skylab.

Siljak, D. D.↗
NASA NTRSDate not supplied

Hot Hydrogen Testing of W-coated dUN Kernels for Nuclear Thermal Propulsion

Tungsten-coated kernels of uranium nitride (UN) are a possible reactor fuel component for nuclear thermal propulsion. Hot hydrogen testing of such kernels (approx. 250μm diameter, coating thickness 5μm) has been performed at temperatures from 1800°C to 2300°C for 30min each, using depleted uranium nitride (dUN). In addition to direct observations and mass loss measurements, the samples were analyzed by XRD and SEM/EDS after each run. Decomposition of dUN and the formation of molten uranium happened at all temperatures tested with reaction rates increasing with temperature, despite the tungsten coating. At 2300°C, the kernels disintegrated and molten uranium seeped through the wall and bottom of the tungsten crucible employed. These hot hydrogen tests allowed to test the integrity of the tungsten coating and addressed the operating limits of the loose tungsten-coated dUN kernels. Improved results are expected upon consolidation of the tungsten-coated dUN kernels in a Mo-W alloy matrix.

nuclear fuels↗
NASA NTRS1991

Accuracy and speed in computing the Chebyshev collocation derivative

We studied several algorithms for computing the Chebyshev spectral derivative and compare their roundoff error. For a large number of collocation points, the elements of the Chebyshev differentiation matrix, if constructed in the usual way, are not computed accurately. A subtle cause is is found to account for the poor accuracy when computing the derivative by the matrix-vector multiplication method. Methods for accurately computing the elements of the matrix are presented, and we find that if the entities of the matrix are computed accurately, the roundoff error of the matrix-vector multiplication is as small as that of the transform-recursion algorithm. Results of CPU time usage are shown for several different algorithms for computing the derivative by the Chebyshev collocation method for a wide variety of two-dimensional grid sizes on both an IBM and a Cray 2 computer. We found that which algorithm is fastest on a particular machine depends not only on the grid size, but also on small details of the computer hardware as well. For most practical grid sizes used in computation, the even-odd decomposition algorithm is found to be faster than the transform-recursion method.

Don, Wai-Sun↗
NASA NTRS1998

An Infrared Spectroscopic and Density Functional Theoretical Investigation of the Reaction Products of Laser-Ablated Scandium and Titanium Atoms with Nitric Oxide

Laser-ablated Sc and Ti atoms have been reacted with NO during condensation in excess argon. Matrix infrared spectra show that the major products are the side bonded Sc[NO] species and the inserted NScO and NTiO molecules based on isotopic substitution (15NI6O and 15N18O) and DIFT calculations of isotopic frequencies, which provide a match for two modes in three isotopic modifications for each molecule. The NScO and NTiO molecules are nitride/oxides with M-O stretching modes only 46-88/cm below the diatomic metal oxides but M-N stretching modes 314-442/cm lower than the diatomic metal nitride molecules. The ScN, ScO, TiN, and TiO molecules are observed as decomposition products. Evidence is also presented for the nitrosyls ScNO and TiNO, the Sc[NO]+ cation, and the NTiO- anion.

Bauschlicher, Charles W., Jr.↗
NASA NTRS2009

Decomposition Algorithm for Global Reachability Analysis on a Time-Varying Graph with an Application to Planetary Exploration

Hot air (Montgolfiere) balloons represent a promising vehicle system for possible future exploration of planets and moons with thick atmospheres such as Venus and Titan. To go to a desired location, this vehicle can primarily use the horizontal wind that varies with altitude, with a small help of its own actuation. A main challenge is how to plan such trajectory in a highly nonlinear and time-varying wind field. This paper poses this trajectory planning as a graph search on the space-time grid and addresses its computational aspects. When capturing various time scales involved in the wind field over the duration of long exploration mission, the size of the graph becomes excessively large. We show that the adjacency matrix of the graph is block-triangular, and by exploiting this structure, we decompose the large planning problem into several smaller subproblems, whose memory requirement stays almost constant as the problem size grows. The approach is demonstrated on a global reachability analysis of a possible Titan mission scenario.

Dijkstra's algorithm↗
NASA NTRS1995

Implementation of a multiblock sensitivity analysis method in numerical aerodynamic shape optimization

A multiblock sensitivity analysis method is applied in a numerical aerodynamic shape optimization technique. The Sensitivity Analysis Domain Decomposition (SADD) scheme which is implemented in this study was developed to reduce the computer memory requirements resulting from the aerodynamic sensitivity analysis equations. Discrete sensitivity analysis offers the ability to compute quasi-analytical derivatives in a more efficient manner than traditional finite-difference methods, which tend to be computationally expensive and prone to inaccuracies. The direct optimization procedure couples CFD analysis based on the two-dimensional thin-layer Navier-Stokes equations with a gradient-based numerical optimization technique. The linking mechanism is the sensitivity equation derived from the CFD discretized flow equations, recast in adjoint form, and solved using direct matrix inversion techniques. This investigation is performed to demonstrate an aerodynamic shape optimization technique on a multiblock domain and its applicability to complex geometries. The objectives are accomplished by shape optimizing two aerodynamic configurations. First, the shape optimization of a transonic airfoil is performed to investigate the behavior of the method in highly nonlinear flows and the effect of different grid blocking strategies on the procedure. Secondly, shape optimization of a two-element configuration in subsonic flow is completed. Cases are presented for this configuration to demonstrate the effect of simultaneously reshaping interfering elements. The aerodynamic shape optimization is shown to produce supercritical type airfoils in the transonic flow from an initially symmetric airfoil. Multiblocking effects the path of optimization while providing similar results at the conclusion. Simultaneous reshaping of elements is shown to be more effective than individual element reshaping due to the inclusion of mutual interference effects.

Lacasse, James M.↗
NASA NTRS2017

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗
NASA NTRS1991

Application of powder densification models to the consolidation processing of composites

Unidirectional fiber reinforced metal matrix composite tapes (containing a single layer of parallel fibers) can now be produced by plasma deposition. These tapes can be stacked and subjected to a thermomechanical treatment that results in a fully dense near net shape component. The mechanisms by which this consolidation step occurs are explored, and models to predict the effect of different thermomechanical conditions (during consolidation) upon the kinetics of densification are developed. The approach is based upon a methodology developed by Ashby and others for the simpler problem of HIP of spherical powders. The complex problem is devided into six, much simpler, subproblems, and then their predicted contributions are added to densification. The initial problem decomposition is to treat the two extreme geometries encountered (contact deformation occurring between foils and shrinkage of isolated, internal pores). Deformation of these two geometries is modelled for plastic, power law creep and diffusional flow. The results are reported in the form of a densification map.

Wadley, H. N. G.↗
NASA NTRS1980

Burning characteristics and fiber retention of graphite/resin matrix composites

Graphite fiber reinforced resin matrix composites were subjected to controlled burning conditions to determine their burning characteristics and fiber retention properties. Small samples were burned with a natural gas fired torch to study the effects of fiber orientation and structural flaws such as holes and slits that were machined into the laminates. Larger laminate samples were burned in a modified heat release rate calorimeter. Unidirectional epoxy/graphite and polyimide/graphite composites and boron powder filled samples of each of the two composite systems were burn tested. The composites were exposed to a thermal radiation of 5.3 Btu/sq ft-sec in air. Samples of each of the unfilled composite were decomposed anaerobically in the calorimeter. Weight loss data were recorded for burning and decomposition times up to thirty-five minutes. The effects of fiber orientation, flaws, and boron filler additives to the resins were evaluated. A high char forming polyimide resin was no more effective in retaining graphite fibers than a low char forming epoxy resin when burned in air. Boron powder additions to both the polyimide and the epoxy resins stabilized the chars and effectively controlled the fiber release.

Bowles, K. J.↗
NASA NTRS1998

Parallel 3D Multi-Stage Simulation of a Turbofan Engine

A 3D multistage simulation of each component of a modern GE Turbofan engine has been made. An axisymmetric view of this engine is presented in the document. This includes a fan, booster rig, high pressure compressor rig, high pressure turbine rig and a low pressure turbine rig. In the near future, all components will be run in a single calculation for a solution of 49 blade rows. The simulation exploits the use of parallel computations by using two levels of parallelism. Each blade row is run in parallel and each blade row grid is decomposed into several domains and run in parallel. 20 processors are used for the 4 blade row analysis. The average passage approach developed by John Adamczyk at NASA Lewis Research Center has been further developed and parallelized. This is APNASA Version A. It is a Navier-Stokes solver using a 4-stage explicit Runge-Kutta time marching scheme with variable time steps and residual smoothing for convergence acceleration. It has an implicit K-E turbulence model which uses an ADI solver to factor the matrix. Between 50 and 100 explicit time steps are solved before a blade row body force is calculated and exchanged with the other blade rows. This outer iteration has been coined a "flip." Efforts have been made to make the solver linearly scaleable with the number of blade rows. Enough flips are run (between 50 and 200) so the solution in the entire machine is not changing. The K-E equations are generally solved every other explicit time step. One of the key requirements in the development of the parallel code was to make the parallel solution exactly (bit for bit) match the serial solution. This has helped isolate many small parallel bugs and guarantee the parallelization was done correctly. The domain decomposition is done only in the axial direction since the number of points axially is much larger than the other two directions. This code uses MPI for message passing. The parallel speed up of the solver portion (no 1/0 or body force calculation) for a grid which has 227 points axially.

Turner, Mark G.↗
NASA NTRS2004

Lunar and Planetary Science XXXV: Astrobiology

The presentations in this session are: 1. A Prototype Life Detection Chip 2. The Geology of Atlantis Basin, Mars, and Its Astrobiological Interest 3. Collecting Bacteria Together with Aerosols in the Martian Atmosphere by the FOELDIX Experimental Instrument Developed with a Nutrient Detector Pattern: Model Measurements of Effectivity 4. 2D and 3D X-ray Imaging of Microorganisms in Meteorites Using Complexity Analysis to Distinguish Field Images of Stromatoloids from Surrounding Rock Matrix in 3.45 Ga Strelley Pool Chert, Western Australia 4. Characterization of Two Isolates from Andean Lakes in Bolivia Short Time Scale Evolution of Microbiolites in Rapidly Receding Altiplanic Lakes: Learning How to Recognize Changing Signatures of Life 5. The Effect of Salts on Electrospray Ionization of Amino Acids in the Negative Mode 6. Determination of Aromatic Ring Number Using Multi-Channel Deep UV Native Fluorescence 7. Microbial D/H Fractionation in Extraterrestrial Materials: Application to Micrometeorites and Mars 8. Carbon Isotope Characteristics of Spring-fed Iron-precipitating Microbial Mats 9. Amino Acid Survival Under Ambient Martian Surface UV Lighting Extraction of Organic Molecules from Terrestrial Material: Quantitative Yields from Heat and Water Extractions 10. Laboratory Detection and Analysis of Organic Compounds in Rocks Using HPLC and XRD Methods 11. Thermal Decomposition of Siderite-Pyrite Assemblages: Implications for Sulfide Mineralogy in Martian Meteorite ALH84001 Carbonate Globules 12. Determination of the Three-Dimensional Morphology of ALH84001 and Biogenic MV-1 Magnetite: Comparison of Results from Electron Tomography and Classical Transmission Electron Microscopy 13. On the Possibility of a Crypto-Biotic Crust on Mars Based on Northern and Southern Ringed Polar Dune Spots 14. Comparative Planetology of the Terrestrial Inner Planets: Implications for Astrobiology 15. A Possible Europa Exobiology 16. A Possible Biogeochemical Model for Titan

Source record↗
NASA NTRS1998

3D Navier-Stokes Time Accurate Solutions Using Multipartitioning Parallel Computation Methodology

A parallel CFD code solving 3D time accurate Navier-Stokes equations with multipartitioning parallel Methodology is being developed in collaboration with Ohio State University within the Air Vehicle Directorate, at Wright Patterson Air Force Base. The advantage of the multipartitioning parallel method is that the domain decomposition will not introduce domain boundaries for the implicit operators. A ring structure data communication is employed so that the implicit time accurate method can be implemented for multi-processors with the same accuracy as for the single processor. No sub-iteration is needed at the domain boundaries. The code has been validated for some typical unsteady flows, which include Coutte Flow, flow passing a cylinder. The code now is being employed for a large scale time accurate wall jet transient flow computation. 'ne preliminary results are promising. The mesh has been refined to capture more details of the flow field. The mesh refinement computation is in progress and would be difficult to successfully implement without the parallel computation techniques used. A modified version of the code with more efficient inversion of the diagonalized block matrix is currently being tested.

Zha, Ge-Cheng↗
NASA NTRS2001

Algorithms and Application of Sparse Matrix Assembly and Equation Solvers for Aeroacoustics

An algorithm for symmetric sparse equation solutions on an unstructured grid is described. Efficient, sequential sparse algorithms for degree-of-freedom reordering, supernodes, symbolic/numerical factorization, and forward backward solution phases are reviewed. Three sparse algorithms for the generation and assembly of symmetric systems of matrix equations are presented. The accuracy and numerical performance of the sequential version of the sparse algorithms are evaluated over the frequency range of interest in a three-dimensional aeroacoustics application. Results show that the solver solutions are accurate using a discretization of 12 points per wavelength. Results also show that the first assembly algorithm is impractical for high-frequency noise calculations. The second and third assembly algorithms have nearly equal performance at low values of source frequencies, but at higher values of source frequencies the third algorithm saves CPU time and RAM. The CPU time and the RAM required by the second and third assembly algorithms are two orders of magnitude smaller than that required by the sparse equation solver. A sequential version of these sparse algorithms can, therefore, be conveniently incorporated into a substructuring for domain decomposition formulation to achieve parallel computation, where different substructures are handles by different parallel processors.

Watson, W. R.↗
NASA NTRS2008

Tungsten and Barium Transport in the Internal Plasma of Hollow Cathodes

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the flow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushedback to the emitter surface by the electric field and drag from the xenon ion flow. Thisbarium ion flux is sufficient to maintain a barium surface coverage at the downstream endgreater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length,so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollowcathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

Hollow cathodes↗
NASA NTRS2009

Barium Depletion in Hollow Cathode Emitters

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the ow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushed back to the emitter surface by the electric field and drag from the xenon ion flow. This barium ion flux is sufficient to maintain a barium surface coverage at the downstream end greater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length, so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollow cathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

plasma discharges↗