Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Detection of strength limiting defects in cellular glasses by dielectric measurements

The foamed glasses that are being considered for structural support in solar thermal energy system silvered glass concentrator mirrors are assessed for strength-limiting defects by a nondestructive technique which measures the capacitance of well defined regions. The feasibility of locating large defects is demonstrated for the case of one type of cellular mass. Observed capacitance variations are believed to be due to small voids or other irregularities in the material. If the defects are equal to or greater than approximately 20 times the glass matrix pore size, they can be resolved with sufficient accuracy to delineate their spatial extent.

Gammell, P. M.↗
NASA NTRS1984

A finite element method for nonlinear forced vibrations of beams

Techniques for defining a finite element model (FEM) for analysis of nonlinear vibrations in beam structures subjected to harmonic excitation are presented. The resulting model covers longitudinal deformation and inertial effects. The nonlinear oscillations of a beam element under forced excitation are modeled by a harmonic force matrix based on first order approximations of the Jacobian elliptic forcing function. Harmonic force and nonlinear stiffness matrices are derived and the nonlinear forced responses of beams are calculated under various boundary conditions. The results of FEM computations for simply-supported and clamped beams show that midplane stretching caused by large deflections increases the nonlinearity. Axially-restrained beams experience only hardening nonlinearity, while axially-free beams have reduced nonlinearity in deformation and inertia and an increase in linearity due to large deflection.

Mei, C.↗
NASA NTRS1989

Uncertainty estimates in geomagnetic field modeling

This paper presents an extension of the conventional uncertainty analysis which characterizes the sources of uncertainty in the coefficients of the geomagnetic field models. The new formalism accounts for the systematic errors introduced by the omission of such sources as the presence of crustal fields, the external fields, and the field from the truncated terms. The usefulness of this formalism depends on two critical conditions. The first of these is the knowledge of the statistical properties of the fields whose parameters are not solved in the analysis, i.e., the crustal field and the external field. The second critical point in the practical use of the method is the approximation used for the weight matrix.

Langel, R. A.↗
NASA NTRS1992

Rate coefficients for the excitation of infrared and ultraviolet lines in C II, N III, and O IV

New and improved collision strengths and Maxwellian-averaged rate coefficients for temperatures between 1000 and 40,000 K are presented for all the prominent transitions in the spectra of C II, N III, and O IV. The collision strengths show extensive autoionization structures that are delineated in detail and which enhance the rate coefficients for several transitions by a considerable amount. Maxwellian-averaged collision strengths are calculated for all possible fine-structure transitions along the states included in the eigenfunction expansion of the target ion. All calculations are carried out in the close-coupling approximation using the R-matrix method as adapted for the Opacity Project.

Blum, Robert D.↗
NASA NTRS1992

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

Work performed during the report period is summarized. The first numerical experiment which was performed on the North Carolina Supercomputer Center's CRAY-YMP machine during the second half of FY92 involved a 36 hour simulation of the CCOPE case study. This first coarse-mesh simulation employed the GMASS model with a 178 x 108 x 32 matrix of grid points spaced approximately 24 km apart. The initial data was comprised of the global 2.5 x 2.5 degree analyses as well as all available North American rawinsonde data valid at 0000 UTC 11 July 1981. Highly-smoothed LFM-derived terrain data were utilized so as to determine the mesoscale response of the three-dimensional atmosphere to weak terrain forcing prior to including the observed highly complex terrain of the northern Rocky Mountain region. It was felt that the model should be run with a spectrum of terrain geometries, ranging from observed complex terrain to no terrain at all, to determine how crucial the terrain was in forcing the mesoscale phenomena. Both convection and stratiform (stable) precipitation were not allowed in this simulation so that their relative importance could be determined by inclusion in forth-coming simulations. A full suite of planetary boundary layer forcing was allowed in the simulation, including surface sensible and latent heat fluxes employing the Blakadar PBL formulation. The details of this simulation, which in many ways could be considered the control simulation, including the important synoptic-scale, meso-alpha scale, and meso-beta scale circulations is described. These results are compared to the observations diagnosed by Koch and his colleagues as well as hypotheses set forth in the project proposal for terrain-influences upon the jet stream and their role in the generation of mesoscale wave phenomenon. The fundamental goal of the analyses being the discrimination among background geostrophic adjustment, terrain influences, and shearing instability in the initiation and maintainance of mesoscale internal wave phenomena. Based upon these findings, FY93 plans are discussed. A review of linear theory and theoretical modeling of a geostrophic zonal wind anomaly is included.

Lin, Yuh-Lang↗
NASA NTRS1993

Excitation rate coefficients and line ratios for the optical and ultraviolet transitions in S II

New calculations are reported for electron excitation collision strengths, rate coefficients, transition probabilities, and line ratios for the astrophysically important optical and UV lines in S II. The collision strengths are calculated in the close coupling approximation using the R-matrix method. The present calculations are more extensive than previous ones, including all transitions among the 12 lowest LS terms and the corresponding 28 fine-structure levels in the collisional-radiative model for S II. While the present rate coefficients for electron impact excitation are within 10-30 percent of the previous values for the low-lying optical transitions employed as density diagnostics of H II regions and nebulae, the excitation rates for the UV transitions 4S super 0 sub 3/2 - 4Psub 1/2,3/2,5/2 differ significantly from earlier calculations, by up to factor of 2. We describe temperature and density sensitive flux ratios for a number of UV lines. The present UV results are likely to be of interest in a more accurate interpretation of S II emission from the Io plasma torus in the magnetosphere of Jupiter, as well as other UV sources observed from the IUE, ASTRO 1, and the HST.

Cai, Wei↗
NASA NTRS1995

Heat Treatment Effects on the Tensile Properties and Microstructures of a SiC/RBSN Composite in Nitrogen

The room-temperature tensile properties and constituent microstructures of a unidirectionally reinforced SiC/reaction bonded silicon nitride (RBSN) composite have been investigated after heat treatments at 1400, 1600, or 1800 C in nitrogen for up to 100 hr. The composite consisted of approximately 24 vol% of aligned 140 micron diameter, continuous length, chemically vapor deposited SiC fibers in an approximately 40% porous silicon nitride matrix. The composites heat treated at 1400 C for up to 100 hr showed elastic modulus, first matrix cracking strength, and ultimate tensile strength values similar to those of the as-fabricated composites, but those heat treated for 1 hr beyond this temperature displayed losses in all three properties. Recrystallization of the SiC fibers, reaction between the carbon-rich interface coating on the fibers and the RBSN matrix, and dissociation of the RBSN matrix are the reasons for the loss of mechanical properties.

Bhatt, R. T.↗
NASA NTRS1999

An Analysis of the Macroscopic Tensile Behavior of a Nonlinear Nylon Reinforced Elastomeric Composite System Using MAC/GMC

A special class of composite laminates composed of soft rubbery matrices and stiff reinforcements made of steel wires or synthetic fibers is examined, where each constituent behaves in a nonlinear fashion even in the small strain domain. Composite laminates made of piles stacked at alternating small orientation angles with respect to the applied axial strain are primarily dominated by the nonlinear behavior of the reinforcing fibers. However; composites with large ply orientations or those perpendicular to the loading axis, will approximate the behavior of the matrix phase and respond in even a more complex fashion for arbitrarily stacked piles. The geometric nonlinearity due to small cord rotations during loading was deemed here to have a second order effect and consequently dropped from any consideration. The user subroutine USRMAT within the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC), was utilized to introduce the constituent material nonlinear behavior. Stress-strain behavior at the macro level was experimentally generated for single and multi ply composites comprised of continuous Nylon-66 reinforcements embedded in a carbon black loaded rubbery matrix. Comparisons between the predicted macro composite behavior and experimental results are excellent when material nonlinearity is included in the analysis. In this paper, a brief review of GMC is provided, along with a description of the nonlinear behavior of the constituents and associated constituent constitutive relations, and the improved macro (or composite) behavior predictions are documented and illustrated.

Assaad, Mahmoud↗
NASA NTRS2000

Photoionization and Recombination

Theoretically self-consistent calculations for photoionization and (e + ion) recombination are described. The same eigenfunction expansion for the ion is employed in coupled channel calculations for both processes, thus ensuring consistency between cross sections and rates. The theoretical treatment of (e + ion) recombination subsumes both the non-resonant recombination ("radiative recombination"), and the resonant recombination ("di-electronic recombination") processes in a unified scheme. In addition to the total, unified recombination rates, level-specific recombination rates and photoionization cross sections are obtained for a large number of atomic levels. Both relativistic Breit-Pauli, and non-relativistic LS coupling, calculations are carried out in the close coupling approximation using the R-matrix method. Although the calculations are computationally intensive, they yield nearly all photoionization and recombination parameters needed for astrophysical photoionization models with higher precision than hitherto possible, estimated at about 10-20% from comparison with experimentally available data (including experimentally derived DR rates). Results are electronically available for over 40 atoms and ions. Photoionization and recombination of He-, and Li-like C and Fe are described for X-ray modeling. The unified method yields total and complete (e+ion) recombination rate coefficients, that can not otherwise be obtained theoretically or experimentally.

Nahar, Sultana N.↗
NASA NTRS1997

Recent Progress in Parallel Schur Complement Preconditioning for Computational Fluid

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a nonoverlapping Schur complement procedure for arbitrary triangulated domains. The triangulation is first partitioned using the METIS multi-level $k$-way partitioning code. This partitioning of the triangulation induces a natural 2x2 partitioning of the demoralization matrix. By considering various inverse approximations of the 2x2 system we have developed a family of robust preconditioning techniques. The performance of these approximations will be discussed and numerous examples shown to illustrate the efficiency of the technique.

Barth, Tim↗
NASA NTRS1997

Progress in Parallel Schur Complement Preconditioning for Computational Fluid Dynamics

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a non-overlapping Schur complement procedure for arbitrary triangulated domains. The ultimate goal of this research is to develop scalable preconditioning algorithms for fluid flow discretizations on parallel computing architectures. In our implementation of the Schur complement preconditioning technique, the triangulation is first partitioned into a number of subdomains using the METIS multi-level k-way partitioning code. This partitioning induces a natural 2X2 partitioning of the p.d.e. discretization matrix. By considering various inverse approximations of the 2X2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 and the SGI Power Challenge array using MPI message passing protocol. A number of example CFD calculations will be presented to illustrate and assess various Schur complement approximations.

Barth, Timothy J.↗
NASA NTRS2005

Optical and Thermo-optical Properties of Polyimide-Single-Walled Carbon Nanotube Films: Experimental Results and Empirical Equations

The incorporation of single-walled carbon nanotubes (SWNTs) into the bulk of space environmentally durable polymers at loading levels greater than or equal to 0.05 wt % has afforded thin films with surface and volume resistivities sufficient for electrostatic charge mitigation. However, the optical transparency at 500 nm decreased and the thermo-optical properties (solar absorptivity and thermal emissivity) increased with increaed SWNT loading. These properties were also dependent on film thickness. The absorbance characteristics of the films as a function of SWNT loading and film thickness were measured and determined to follow the classical Beer-Lambert law. Based on these results, an empirical relationship was derived and molar absorptivities determined for both the SWNTs and polymer matrix to provide a predictive approximation of these properties. The molar absorptivity determined for SWNTs dispersed in the polymer was comparable to reported solution determined values for HiPco SWNTs.

Smith, Joseph G., Jr.↗
NASA NTRS2012

A Newton-Raphson Method Approach to Adjusting Multi-Source Solar Simulators

NASA Glenn Research Center has been using an in house designed X25 based multi-source solar simulator since 2003. The simulator is set up for triple junction solar cells prior to measurements b y adjusting the three sources to produce the correct short circuit current, lsc, in each of three AM0 calibrated sub-cells. The past practice has been to adjust one source on one sub-cell at a time, iterating until all the sub-cells have the calibrated Isc. The new approach is to create a matrix of measured lsc for small source changes on each sub-cell. A matrix, A, is produced. This is normalized to unit changes in the sources so that Ax(delta)s = (delta)isc. This matrix can now be inverted and used with the known Isc differences from the AM0 calibrated values to indicate changes in the source settings, (delta)s = A ·'x.(delta)isc This approach is still an iterative one, but all sources are changed during each iteration step. It typically takes four to six steps to converge on the calibrated lsc values. Even though the source lamps may degrade over time, the initial matrix evaluation i s not performed each time, since measurement matrix needs to be only approximate. Because an iterative approach is used the method will still continue to be valid. This method may become more important as state-of-the-art solar cell junction responses overlap the sources of the simulator. Also, as the number of cell junctions and sources increase, this method should remain applicable.

Snyder, David B.↗
NASA NTRS2013

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties↗
DOE OSTI2023

Exact Gaussian processes for massive datasets via non-stationary sparsity-discovering kernels

Abstract A Gaussian Process (GP) is a prominent mathematical framework for stochastic function approximation in science and engineering applications. Its success is largely attributed to the GP’s analytical tractability, robustness, and natural inclusion of uncertainty quantification. Unfortunately, the use of exact GPs is prohibitively expensive for large datasets due to their unfavorable numerical complexity of $$O(N^3)$$ O ( N 3 ) in computation and $$O(N^2)$$ O ( N 2 ) in storage. All existing methods addressing this issue utilize some form of approximation—usually considering subsets of the full dataset or finding representative pseudo-points that render the covariance matrix well-structured and sparse. These approximate methods can lead to inaccuracies in function approximations and often limit the user’s flexibility in designing expressive kernels. Instead of inducing sparsity via data-point geometry and structure, we propose to take advantage of naturally-occurring sparsity by allowing the kernel to discover—instead of induce—sparse structure. The premise of this paper is that the data sets and physical processes modeled by GPs often exhibit natural or implicit sparsities, but commonly-used kernels do not allow us to exploit such sparsity. The core concept of exact, and at the same time sparse GPs relies on kernel definitions that provide enough flexibility to learn and encode not only non-zero but also zero covariances. This principle of ultra-flexible, compactly-supported, and non-stationary kernels, combined with HPC and constrained optimization, lets us scale exact GPs well beyond 5 million data points.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Efficient solution of parabolic equations by Krylov approximation methods

Numerical techniques for solving parabolic equations by the method of lines is addressed. The main motivation for the proposed approach is the possibility of exploiting a high degree of parallelism in a simple manner. The basic idea of the method is to approximate the action of the evolution operator on a given state vector by means of a projection process onto a Krylov subspace. Thus, the resulting approximation consists of applying an evolution operator of a very small dimension to a known vector which is, in turn, computed accurately by exploiting well-known rational approximations to the exponential. Because the rational approximation is only applied to a small matrix, the only operations required with the original large matrix are matrix-by-vector multiplications, and as a result the algorithm can easily be parallelized and vectorized. Some relevant approximation and stability issues are discussed. We present some numerical experiments with the method and compare its performance with a few explicit and implicit algorithms.

Gallopoulos, E.↗
NASA NTRS1985

Interdiffusional effects between tungsten fibers and an iron-nickel-base alloy

Tungsten fibers in the INCOLOY 903 alloy were annealed for over 100 hours at 1038 C and 1200 C. It was found that interdiffusion results in the formation of a reaction zone. SEM-EDS probe analysis showed that the chemistries across this zone were constant, suggesting the zone was a compound phase. The composition of the compound was estimated to be that of a mu-type phase. The local chemistry (in atomic percent) at the reaction zone/alloy matrix interface was found to be approximately 8 pct W, 1.2 pct Nb, 40 pct Fe, 14 pct Co, and 36 pct Ni. In addition, recrystallization was observed in both the remaining tungsten fiber and the nearby INCOLOY 903 matrix after annealing at 1200 C, but not at 1038 C. The results of this study suggest that reaction zone growth kinetics can be minimized by the reduction of Co and Fe and the increase of W in the matrix alloy.

Caulfield, T.↗
NASA NTRS1995

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗