Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix equations”

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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Registering and resampling images in STSDAS

Registering different images can be difficult, especially if the images to be registered are images at different wavelengths, where features in one image may look entirely different or be absent from the second image. Using two new packages soon to be added to the STSDAS package, REGISTER and RESAMPLE, this job is done automatically. The REGISTER package allows the user to determine the amount of translation, rotation, and/or magnification needed to make two images, spectra, or time series congruent. The methods implemented to compute the registration parameters use: a set of the pixel coordinates of the same features identified in two files, or the FITS coordinate transformation parameters in the headers of two data files, or a single feature identified as the peak of a cross-correlation between two vectors. The coefficients describing the registration are defined by the equations (for a two-dimensional image): x = a + bx + cy, y = d + ex + fy, where (X,Y) are the pixel coordinates of a feature in the reference image, (x,y) are the pixel coordinates of a feature in the secondary image, and the computed coefficients are a, b, c, d, e, and f. Results may be produced by linking the output of REGISTER to RESAMPLE in a command language procedure. The output from REGISTER and the input to RESAMPLE consists of a matrix of coefficients (a through f above) fully specifying the registration. The RESAMPLE package resamples simple vector or image data for a given amount of translation, rotation (images only), and magnification, or reflection of the science data. Specific options included are: image rotation about the FITS reference pixel; scale changes, i.e. magnification or demagnification (for images, independently on both axes); simple translation; reflection (for images, about one or both axes); and resampling and registration to a reference data set. Output from the RESAMPLE task is the resampled image which may then be displayed and compared with the reference image.

Williamson, R. L., II↗
NASA NTRS2018

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗
NASA NTRS1980

Linear state feedback, quadratic weights, and closed loop eigenstructures

Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalue and the directional derivatives of closed loop eigenvectors. An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties. An algorithm is presented that can be used to select a feedback gain matrix for the linear state feedback problem which produces a specified asymptotic eigenstructure. Another algorithm is given to compute the asymptotic eigenstructure properties inherent in a given set of quadratic weights. Finally, it is shown that optimal root loci for nongeneric problems can be approximated by generic ones in the nonasymptotic region.

Thompson, P. M.↗
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 NTRS1983

A non-linear viscoelastic characterisation method for matrix resin composites

Lifetime prediction of matrix resin composites is an important problem because of their viscoelastic character. If the Time-Temperature-Stress-Superposition principle is valid, results from short time tests at high stress levels can give the necessary information for the prediction of long time behavior. At such stress levels the viscoelastic response is nonlinear. In this contribution the basis of a computer routine is presented starting from a thermodynamical formulation given by Schapery. This method gives in an accurate way the different material parameters in that nonlinear constitutive equation for different stresses and temperatures.

Hiel, C. C.↗
NASA NTRS1978

Modal content of noise generated by a coaxial jet in a pipe

Noise generated by air flow through a coaxial obstruction in a long, straight pipe was investigated with concentration on the modal characteristics of the noise field inside the pipe and downstream of the restriction. Two measurement techniques were developed for separation of the noise into the acoustic duct modes. The instantaneous mode separation technique uses four microphones, equally spaced in the circumferential direction, at the same axial location. The time-averaged mode separation technique uses three microphones mounted at the same axial location. A matrix operation on time-averaged data produces the modal pressure levels. This technique requires the restrictive assumption that the acoustic modes are uncorrelated with each other. The measured modal pressure spectra were converted to modal power spectra and integrated over the frequency range 200-6000 Hz. The acoustic efficiency levels (acoustic power normalized by jet kinetic energy flow), when plotted vs. jet Mach number, showed a strong dependence on the ratio of restriction diameter to pipe diameter. The acoustic energy flow analyses based on the thermodynamic energy equation and on the results of Mohring both resulted in orthogonality properties for the eigenfunctions of the radial mode shape equation. These orthogonality relationships involve the eigenvalues and derivatives of the radial mode shape functions.

Kerschen, E. J.↗
NASA NTRS1999

Inversion of Gravity Fields From the Spacecraft Orbital Data Using an Adjoint Operator Approach

In perturbation approximation, the forward problem of orbital dynamics (equations with initial conditions) is linear with respect to variations of coordinates and/or velocities of the spacecraft and to corresponding variations of the gravity field in the models used. The linear operator adjoint to the linear operator of such forward problem turns out to be instrumental in inversion of differences between observed and predicted coordinates/velocities in terms of the updates of harmonics in the initial gravity field model. Based on this approach, the solution of resulting adjoint problem of orbital dynamics can be used to directly evaluate the matrix of partial derivatives of observable differences with respect to the gravity field harmonics. General discussion of the adjoint problem of orbital dynamics is given and an example of a mathematical formalism for the practical retrieval algorithm is presented.

Ustinov, E. A.↗
NASA NTRS1980

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1)the elimination of large matrix storage requirements, and (2)shorter solution times under most conditions.

Baumeister, K. J.↗
NASA NTRS1980

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing/acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1) the elimination of large matrix storage requirements, and (2) shorter solution times under most conditions.

Baumeister, K. J.↗
NASA NTRS2012

Scattering and Absorption Properties of Polydisperse Wavelength-sized Particles Covered with Much Smaller Grains

Using the results of direct, numerically exact computer solutions of the Maxwell equations, we analyze scattering and absorption characteristics of polydisperse compound particles in the form of wavelength-sized spheres covered with a large number of much smaller spherical grains.The results pertain to the complex refractive indices1.55 + i0.0003,1.55 + i0.3, and 3 + i0.1. We show that the optical effects of dusting wavelength-sized hosts by microscopic grains can vary depending on the number and size of the grains as well as on the complex refractive index. Our computations also demonstrate the high efficiency of the new superposition T-matrix code developed for use on distributed memory computer clusters.

Dlugach, Jana M.↗
NASA NTRS1991

Parallel variable-band Choleski solvers for computational structural analysis applications on vector multiprocessor supercomputers

A Choleski method used to solve linear systems of equations that arise in large scale structural analyses is described. The method uses a novel variable-band storage scheme and is structured to exploit fast local memory caches while minimizing data access delays between main memory and vector registers. Several parallel implementations of this method are described for the CRAY-2 and CRAY Y-MP computers demonstrating the use of microtasking and autotasking directives. A portable parallel language, FORCE, is also used for two different parallel implementations, demonstrating the use of CRAY macrotasking. Results are presented comparing the matrix factorization times for three representative structural analysis problems from runs made in both dedicated and multi-user modes on both the CRAY-2 and CRAY Y-MP computers. CPU and wall clock timings are given for the various parallel methods and are compared to single processor timings of the same algorithm. Computation rates over 1 GIGAFLOP (1 billion floating point operations per second) on a four processor CRAY-2 and over 2 GIGAFLOPS on an eight processor CRAY Y-MP are demonstrated as measured by wall clock time in a dedicated environment. Reduced wall clock times for the parallel methods relative to the single processor implementation of the same Choleski algorithm are also demonstrated for runs made in multi-user mode.

Poole, E. L.↗
NASA NTRS2000

Fatigue, Creep-Fatigue, and Thermomechanical Fatigue Life Testing of Alloys

The fatigue crack initiation resistance of an alloy is determined by conducting a series of tests over a range of values of stress amplitude or strain range. The observed number of cycles to failure is plotted against the stress amplitude or strain range to obtain a fatigue curve. The fatigue properties quoted for an alloy are typically the constants used in the equation(s) that describe the fatigue curve. Fatigue lives of interest may be as low as 10(exp 2) or higher than 10(exp 9) cycles. Because of the enormous scatter associated with fatigue, dozens of tests may be needed to confidently establish a fatigue curve, and the cost may run into several thousands of dollars. To further establish the effects on fatigue life of the test temperature, environment, alloy condition, mean stress effects, creep-fatigue effects, thermomechanical cycling, etc. requires an extraordinarily large and usually very costly test matrix. The total effort required to establish the fatigue resistance of an alloy should not be taken lightly. Fatigue crack initiation tests are conducted on relatively small and presumed to be initially crack-free, samples of an alloy that are intended to be representative of the alloy's metallurgical and physical condition. Generally, samples are smooth and have uniformly polished surfaces within the test section. Some may have intentionally machined notches of well-controlled geometry, but the surface at the root of the notch is usually not polished. The purpose of polishing is to attain a reproducible surface finish. This is to eliminate surface finish as an uncontrolled variable. Representative test specimen geometries will be discussed later. Test specimens are cyclically loaded until macroscopically observable cracks initiate and eventually grow to failure. Normally, the fatigue failure life of a specimen is defined as the number of cycles to separation of the specimen into two pieces. Alternative definitions are becoming more common, particularly for low-cycle fatigue testing, wherein some prescribed indication of impending failure due to cracking is adopted. Specific criteria will be described later. As a rule, cracks that develop during testing are not measured nor are the test parameters intentionally altered owing to the presence of cracking.

Halford, Gary R.↗
NASA NTRS2010

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.↗
NASA NTRS1972

Determination of stability derivatives from flight data using a Newton-Raphson minimization technique

A modified Newton-Raphson or quasilinearization minimization technique for determining stability derivatives from flight data was developed and compared with simple-equations, analog-matching, least-squares, and Shinbrot methods of analysis. For the data analyzed, the solutions computed by using the estimates obtained from the Newton-Raphson technique fit the data and determined coefficients adequately. A further modification to include a priori information was found to be useful. A model statistically similar to the flight data was analyzed using the same methods (excluding analog matching), and the Newton-Raphson technique was found to yield superior estimates. An approximate Cramer-Rao bound was compared with the error covariance matrix of the model and was found to provide information about the reliability of the individual estimates obtained. The technique was successfully applied to data obtained from a light airplane, a large supersonic airplane, and a lifting body vehicle. It was shown that the reliability of the estimates of a given coefficient obtained from these vehicles depends upon the data analyzed.

Iliff, K. W.↗
NASA NTRS1974

A generalized-capacity-matrix technique for computing aerodynamic flows

A numerical generalized-capacity-matrix technique is developed for application to aerodynamic flow computations. This technique allows the very fast direct (noniterative) numerical elliptic solvers to be used in problems with arbitrary internal boundaries and with a wide class of boundary conditions, including numerical application of the Kutta condition on an airfoil without iteration. Accuracy, speed, and usefulness of the technique are demonstrated with linear problems for potential flows over airfoil shapes. The method's main advantages, however, can be exploited within iterative procedures for a variety of complex flow problems governed by systems of equations not necessarily elliptic or linear.

Martin, E. D.↗
NASA NTRS1993

Nested Krylov methods and preserving the orthogonality

Recently the GMRESR inner-outer iteraction scheme for the solution of linear systems of equations was proposed by Van der Vorst and Vuik. Similar methods have been proposed by Axelsson and Vassilevski and Saad (FGMRES). The outer iteration is GCR, which minimizes the residual over a given set of direction vectors. The inner iteration is GMRES, which at each step computes a new direction vector by approximately solving the residual equation. However, the optimality of the approximation over the space of outer search directions is ignored in the inner GMRES iteration. This leads to suboptimal corrections to the solution in the outer iteration, as components of the outer iteration directions may reenter in the inner iteration process. Therefore we propose to preserve the orthogonality relations of GCR in the inner GMRES iteration. This gives optimal corrections; however, it involves working with a singular, non-symmetric operator. We will discuss some important properties, and we will show by experiments that, in terms of matrix vector products, this modification (almost) always leads to better convergence. However, because we do more orthogonalizations, it does not always give an improved performance in CPU-time. Furthermore, we will discuss efficient implementations as well as the truncation possibilities of the outer GCR process. The experimental results indicate that for such methods it is advantageous to preserve the orthogonality in the inner iteration. Of course we can also use iteration schemes other than GMRES as the inner method; methods with short recurrences like GICGSTAB are of interest.

Desturler, Eric↗
NASA NTRS1990

A study of void effects on the interlaminar shear strength of unidirectional graphite fiber reinforced composites

A study was conducted to evaluate the effect of voids on the interlaminar shear strength (ILSS) of a polyimide matrix composite system. The graphite/PRM-15 composite was chosen for study because of the extensive amount of experience that has been amassed in the processing of this material. Composite densities and fiber contents of more than thirty different laminates were measured along with ILSS. Void contents were calculated and the void geometry and distribution were noted using microscopic techniques such as those used in metallography. It was found that there was a good empirical correlation between ILSS and composite density. The most acceptable relationship between the ILSS and density was found to be a power equation which closely resembles theoretically derived expressions. An increase in scatter in the strength data was observed as the void content increased. In laminates with low void content, the void appears to be more segregated in one area of the laminate. It was found that void free composites could be processed in matched metal die molds at pressures greater than 1.4 and less than 6.9 MPa.

Bowles, Kenneth J.↗
NASA NTRS2011

Numerical Implementation of a Multiple-ISV Thermodynamically-Based Work Potential Theory for Modeling Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗