Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Computationally efficient modeling and simulation of large scale systems

A method of simulating operation of a VLSI interconnect structure having capacitive and inductive coupling between nodes thereof. A matrix X and a matrix Y containing different combinations of passive circuit element values for the interconnect structure are obtained where the element values for each matrix include inductance L and inverse capacitance P. An adjacency matrix A associated with the interconnect structure is obtained. Numerical integration is used to solve first and second equations, each including as a factor the product of the inverse matrix X.sup.1 and at least one other matrix, with first equation including X.sup.1Y, X.sup.1A, and X.sup.1P, and the second equation including X.sup.1A and X.sup.1P.

Jain, Jitesh↗
NASA NTRS2012

Computationally efficient modeling and simulation of large scale systems

A method of simulating operation of a VLSI interconnect structure having capacitive and inductive coupling between nodes thereof. A matrix X and a matrix Y containing different combinations of passive circuit element values for the interconnect structure are obtained where the element values for each matrix include inductance L and inverse capacitance P. An adjacency matrix A associated with the interconnect structure is obtained. Numerical integration is used to solve first and second equations, each including as a factor the product of the inverse matrix X.sup.-1 and at least one other matrix, with first equation including X.sup.-1Y, X.sup.-1A, and X.sup.-1P, and the second equation including X.sup.-1A and X.sup.-1P.

Jain, Jitesh↗
NASA NTRS2004

Laboratory Evidence for a Key Intermediate in the Venus Atmosphere: Peroxychloroformyl Radical

For two decades, the peroxychloroformyl radical, ClC(O)OO, has played a central role in models of the chemical stability of the Venus atmosphere. No confirmation, however, has been possible in the absence of laboratory measurements for ClC(O)OO. We report the isolation of ClC(O)OO in a cryogenic matrix and its infrared and ultraviolet spectral signatures. These experiments show that ClC(O)OO is thermally and photolytically stable in the Venus atmosphere. These experimental discoveries validate the existence of ClC(O)OO, confirm several longstanding model assumptions, and provide a basis for the astronomical search for this important radical species.

photochemistry↗
NASA NTRS1994

Conformational Characteristics of Poly(tetrafluoroethylene) (PTFE) Based Upon Ab Initio Electronic Structure Calculations on Model Molecules

Conformational energy contours of perfluoroalkanes, determined from ab initio calculations, confirm the well-known spitting of trans states into two minima at plus or minus 17 degrees but also show that the gauche states split as well, with minima at plus or minus 124 degrees and plus or minus 84 in order to relieve steric crowding. The directions of such split distortions from the perfectly staggered states are strongly coupled for adjacent pairs of bonds in a manner identical to the intradyad pair for poly (isobutylene) chains. These conformational characteristics are fully represented by a six-state rotational isomeric state (RIS) model for PTFE comprised of t(+), t(-), g(sup +)+, g(sup +)-, g(sup -) + and g(sup -)-states, located at the split energy minima. The resultant 6 x 6 statistical weight matrix is described by first-order interaction parameters for the g+(+) (ca. 0.6 kcal/mol) and g+- (ca. 2.0 kcal/mol) states, and second order parameters for the g(sup +)+g(sup +)+ (ca 0.6 kcal/mol) and g(sup +)+g(sup -)+ (ca. 1.0 kcal/mol) states. This six-state RIS model, without adjustment of the geometric or energy parameters as determined from the ab initio calculations, predicts the unperturbed chain dimensions and the fraction of gauche bonds as a function of temperature for PTFE in good agreement with available experimental values.

Smith, Grant D.↗
NASA NTRS2017

X-Ray Spectro-Polarimetry with Photoelectric Polarimeters

We derive a generalization of forward fitting for X-ray spectroscopy to include linear polarization of X-ray sources, appropriate for the anticipated next generation of space-based photoelectric polarimeters. We show that the inclusion of polarization sensitivity requires joint fitting to three observed spectra, one for each of the Stokes parameters, I(E), U(E), and Q(E). The equations for StokesI (E) (the total intensity spectrum) are identical to the familiar case with no polarization sensitivity, and for which the model-predicted spectrum is obtained by a convolution of the source spectrum, F (E), with the familiar energy response function,(E) R(E,E), where (E) and R(E,E) are the effective area and energy redistribution matrix, respectively. In addition to the energy spectrum, the two new relations for U(E) and Q(E) include the source polarization fraction and position angle versus energy, a(E), and 0(E), respectively, and the model-predicted spectra for these relations are obtained by a convolution with the modulated energy response function, (E)(E) R(E,E), where(E) is the energy-dependent modulation fraction that quantifies a polarimeters angular response to 100 polarized radiation. We present results of simulations with response parameters appropriate for the proposed PRAXyS Small Explorer observatory to illustrate the procedures and methods, and we discuss some aspects of photoelectric polarimeters with relevance to understanding their calibration and operation.

polarization↗
NASA NTRS1991

An improved error assessment for the GEM-T1 gravitational model

Several tests have been designed to estimate the correct error variances for the GEM-T1 gravitational solution that was derived exclusively from satellite tracking data. The basic method uses both independent and dependent subset data solutions and produces a coefficient-by-coefficient estimate of the model uncertainties. The GEM-T1 errors have been further analyzed using a method based on eigenvalue-eigenvector analysis, which calibrates the entire covariance matrix. Dependent satellite data sets and independent altimetric, resonant satellite, and surface gravity data sets all confirm essentially the same error assessment. The calibration test results yield very stable calibration factors, which vary only by approximately 10 percent over the range of tests performed. Based on these calibrated error estimates, GEM-T1 is a significantly improved solution, which to degree and order 8 is twice as accurate as earlier satellite derived models like GEM-L2. Also, by being complete to degree and order 36, GEM-T1 is more complete and has significantly reduced aliasing effects that were present in previous models.

Lerch, F. J.↗
NASA NTRS1977

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

An accurate mathematical model for sound propagation in axisymmetric aircraft engine ducts with compressible mean flow is reported. The model is based on the usual perturbation of the basic fluid mechanics equations for small motions. Mean flow parameters are derived in the absence of fluctuating quantities and are then substituted into the equations for the acoustic quantities which were linearized by eliminating higher order terms. Mean swirl is assumed to be zero from the restriction of axisymmetry. A linear rectangular serendipity element is formulated from these equations using a Galerkin procedure and assembled in a special purpose computer program in which the matrix map for a rectangular mesh was specifically coded. Representations of the fluctuating quantities, mean quantities and coordinate transformations are isoparametric. The global matrix is solved by foreward and back substitution following an L-U decomposition with pivoting restricted internally to the blocks. Results from the model were compared with results from several alternative analyses and yielded satisfactory agreement.

Abrahamson, A. L.↗
NASA NTRS2017

Environmental Barrier Coating Fracture, Fatigue and High-Heat-Flux Durability Modeling and Stochastic Progressive Damage Simulation

Advanced environmental barrier coatings will play an increasingly important role in future gas turbine engines because of their ability to protect emerging light-weight SiC/SiC ceramic matrix composite (CMC) engine components, further raising engine operating temperatures and performance. Because the environmental barrier coating systems are critical to the performance, reliability and durability of these hot-section ceramic engine components, a prime-reliant coating system along with established life design methodology are required for the hot-section ceramic component insertion into engine service. In this paper, we have first summarized some observations of high temperature, high-heat-flux environmental degradation and failure mechanisms of environmental barrier coating systems in laboratory simulated engine environment tests. In particular, the coating surface cracking morphologies and associated subsequent delamination mechanisms under the engine level high-heat-flux, combustion steam, and mechanical creep and fatigue loading conditions will be discussed. The EBC compostion and archtechture improvements based on advanced high heat flux environmental testing, and the modeling advances based on the integrated Finite Element Analysis Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program will also be highlighted. The stochastic progressive damage simulation successfully predicts mud flat damage pattern in EBCs on coated 3-D specimens, and a 2-D model of through-the-thickness cross-section. A 2-parameter Weibull distribution was assumed in characterizing the coating layer stochastic strength response and the formation of damage was therefore modeled. The damage initiation and coalescence into progressively smaller mudflat crack cells was demonstrated. A coating life prediction framework may be realized by examining the surface crack initiation and delamination propagation in conjunction with environmental degradation under high-heat-flux and environment load test conditions.

Mechanical creep and fatigue loading condition↗
NASA NTRS1996

Optimization of Residual Stresses in MMC's through Process Parameter Control and the use of Heterogeneous Compensating/Compliant Interfacial Layers. OPTCOMP2 User's Guide

A user's guide for the computer program OPTCOMP2 is presented in this report. This program provides a capability to optimize the fabrication or service-induced residual stresses in unidirectional metal matrix composites subjected to combined thermomechanical axisymmetric loading by altering the processing history, as well as through the microstructural design of interfacial fiber coatings. The user specifies the initial architecture of the composite and the load history, with the constituent materials being elastic, plastic, viscoplastic, or as defined by the 'user-defined' constitutive model, in addition to the objective function and constraints, through a user-friendly data input interface. The optimization procedure is based on an efficient solution methodology for the inelastic response of a fiber/interface layer(s)/matrix concentric cylinder model where the interface layers can be either homogeneous or heterogeneous. The response of heterogeneous layers is modeled using Aboudi's three-dimensional method of cells micromechanics model. The commercial optimization package DOT is used for the nonlinear optimization problem. The solution methodology for the arbitrarily layered cylinder is based on the local-global stiffness matrix formulation and Mendelson's iterative technique of successive elastic solutions developed for elastoplastic boundary-value problems. The optimization algorithm employed in DOT is based on the method of feasible directions.

Pindera, Marek-Jerzy↗
NASA NTRS1975

High Reynolds number tests of a C-141A aircraft semispan model to investigate shock-induced separation

Results from a high Reynolds number transonic wind tunnel investigation are presented. Wing chordwise pressure distributions were measured over a matrix of Mach numbers and angles-of-attack for which shock-induced separations are known to exist. The range of Reynolds number covered by these data nearly spanned the gap between previously available wind tunnel and flight test data. The results are compared with both flight and low Reynolds number data, and show that use of the semispan test technique produced good correlation with the prior data at both ends of the Reynolds number range, but indicated strong sensitivity to details of the test setup.

Blackerby, W. T.↗
NASA NTRS2011

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Manchiraju, Sivom↗
NASA NTRS1989

Orthogonal canonical forms for second-order systems

The authors prove that a linear second-order system with arbitrary damping cannot be reduced to Hessenberg-triangular form by means of orthogonal transformations, while this reduction is always possible for the modal damping commonly assumed for models of flexible structures. The type of canonical form obtainable by means of orthogonal transformations acting on a second-order system is heavily dependent on the type of damping considered. If the damping matrix is merely positive semi-definite symmetric, it is generally not possible to obtain a reduction to Hessenberg-triangular form, while this reduction is trivial for zero or Rayleigh damping. If damping is modal, however, as is commonly assumed in structural models, the reduction exists and is nontrivial. Furthermore, reduction to triangular second-order Schur form is always possible for such damping: this canonical form appears likely to have applications to second-order system theory.

Williams, Trevor↗
NASA NTRSDate not supplied

Turbulence Effects on Shaped Booms: Finite Impulse Response Filter Development

Numerical simulations of propagation through turbulent atmospheres can quantify effects on ground waveforms, but such simulations are computationally expensive. To enable quick turnaround analyses as required by NASA’s Quesst Mission, updating the N-wave filtering approach developed by researchers at The Pennsylvania State University to include shaped booms is proposed as an alternative method for estimating turbulence effects on acoustic metrics more quickly. Beginning with a nearfield pressure cylinder modeled after the on-design X-59 configuration, a database of propagation results at 45 turbulence conditions was compiled using nonlinear turbulence propagation modeling code (KZKFourier) and used as input to a process for generating finite impulse response (FIR) filters. Ground waveforms distorted by turbulence were selected to represent mean and mean±standard deviation levels for six metrics, and corresponding FIR filters were generated through a matrix deconvolution process. In order to evaluate how well the FIR filters perform, additional KZKFourier verification cases were devised with different input conditions, and results used as a benchmark. Convolution of shaped boom waveforms modeled using nonturbulent propagation simulations with the new FIR filters showed better agreement on average with KZKFourier statistical results than the N-wave-based FIR filters.

turbulence↗
NASA NTRS1975

Geological model for Boulder 1 at Station 2, South Massif, Valley of Taurus-Littrow

A possible geological model for the origin and history of the materials that make up Boulder 1 is proposed on the basis of firm and probable regional, local, and boulder geological constraints. These constraints are described in detail, unresolved questions are considered, and a model is presented which appears to satisfy all the firm constraints and most of the probable constraints. According to this model, the crystallization of plagioclase and other ANT-suite phases now present in the boulder as clasts and matrix materials took place during the melted-shell stage of lunar history; the original rocks were greatly modified during the cratered-highland stage; and the events that determined the major characteristics of the boulder occurred during the large-basin stage.

Schmitt, H. H.↗
NASA NTRS1983

On identifiability of flexible structure parameters

This report investigates the identifiability of modal parameters of flexible structures. Expressions are derived for Cramer-Rao lower bounds for the modal parameters, that is, frequencies, damping ratios and mode shapes or slopes. The optimal initial state, which maximizes the trace of the Fisher information matrix in the absence of persistent input, is obtained. The concepts discussed are applied to a finite-element model of the 122 meter hoop/column antenna. The numerical results show that the identifiability of the structural frequencies is excellent, followed by that of the damping ratios and the mode-slopes.

Joshi, S. M.↗
NASA NTRS1983

Robot arm dynamic model reduction for control

General methods are described by which the mathematical complexities of explicit and exact state equations of robot arms can be reduced to a simplified and compact state equation representation without introducing significant errors into the robot arm dynamic model. The model reduction methods are based on homogeneous coordinates and on the Langrangian algorithm for robot arm dynamics, and utilize matrix, vector and numeric analysis techniques. The derivation of differential vector representation of centripetal and Coriolis forces which has not yet been established in the literature is presented.

Bejczy, A. K.↗
NASA NTRS1987

A creep-rupture model of filament-wound spherical pressure vessels

The creep-rupture model is that of a quasi-isotropic filament-wound spherical pressure vessel, subjected to internal pressurization. The matrix material of the composite system is assumed to be linearly viscoelastic. Internal damage resulting from the relaxation of the matrix and the corresponding increase in microcracks is represented by a functional relationship between circumferential strain and transverse modulus. The numerical solution to this nonlinear problem is an iterative technique, whereby the elastic-viscoelastic correspondence principle is employed. In the Laplace domain, the associated elastic solution is obtained and this solution is inverted by the multidata method to yield the time-dependent solution.

Dozier, Jan D.↗
NASA NTRS1988

A general-purpose approach to computer-aided dynamic analysis of a flexible helicopter

A general purpose mathematical formulation is described for dynamic analysis of a helicopter consisting of flexible and/or rigid bodies that undergo large translations and rotations. Rigid body and elastic sets of generalized coordinates are used. The rigid body coordinates define the location and the orientation of a body coordinate frame (global frame) with respect to an inertial frame. The elastic coordinates are introduced using a finite element approach in order to model flexible components. The compatibility conditions between two adjacent elements in a flexible body are imposed using a Boolean matrix, whereas the compatibility conditions between two adjacent bodies are imposed using the Lagrange multiplier approach. Since the form of the constraint equations depends upon the type of kinematic joint and involves only the generalized coordinates of the two participating elements, then a library of constraint elements can be developed to impose the kinematic constraint in an automated fashion. For the body constraints, the Lagrange multipliers yield the reaction forces and torques of the bodies at the joints. The virtual work approach is used to derive the equations of motion, which are a system of differential and algebraic equations that are highly nonlinear. The formulation presented is general and is compared with hard-wired formulations commonly used in helicopter analysis.

Agrawal, Om P.↗