SEARCH · Engineering Papers

Results for “Derived Matrices”

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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS2011

TPSAS-NF1676L-12839-DND

Hyper-spectral remote sensors, such as Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer (IASI), provides top of atmospheric radiances with high information content on atmospheric and surface properties. In order to analyze these data in real time, fast, and accurate forward and inverse models are needed. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The method has been successfully applied to AIRS, IASI, and NAST-I data. By applying the same method to the current and future hyperspectral sounders, the derived products will have less errors due to biases introduced by different retrieval methods. We have compared the retrieved products to radiosondes, aircraft measurements and other validation dataset. Simulations have done to quantify retrieval errors associated with the retrievals. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a subset of the channels. The PCRTM forward model has also been used to performance end-to-end sensor performance simulations of the Climate Absolute Radiance and Refractivity Observatory (CLARREO).

Xu Liu↗

DOE OSTI2024

The geometry of the modular bootstrap

Abstract We explore the geometry behind the modular bootstrap and its image in the space of Taylor coefficients of the torus partition function. In the first part, we identify the geometry as an intersection of planes with the convex hull of moment curves onR + ⊗ℤ, with boundaries characterized by the total positivity of generalized Hankel matrices. We phrase the Hankel constraints as a semi-definite program, which has several advantages, such as the validity of bounds irrespective of spin truncation. We derive bounds on the gap, twist-gap, and the space of Taylor coefficients themselves. We find that if the gap is above$$ {\Delta }_{\textrm{gap}}^{\ast } $$ ∆ gap ∗ , where$$ \frac{c-1}{12}<{\Delta}_{\textrm{gap}}^{\ast }<\frac{c}{12} $$ c − 1 12 < Δ gap ∗ < c 12 , all coefficients become bounded on both sides and kinks develop in the space. In the second part, we propose an analytic method of imposing the integrality condition for the degeneracy number in the spinless bootstrap, which leads to a non-convex geometry. We find that even at very low derivative order this condition rules out regions otherwise allowed by bootstraps at high derivative order.

Physics↗

NASA NTRS1986

Sensitivity analysis and approximation methods for general eigenvalue problems

Optimization of dynamic systems involving complex non-hermitian matrices is often computationally expensive. Major contributors to the computational expense are the sensitivity analysis and reanalysis of a modified design. The present work seeks to alleviate this computational burden by identifying efficient sensitivity analysis and approximate reanalysis methods. For the algebraic eigenvalue problem involving non-hermitian matrices, algorithms for sensitivity analysis and approximate reanalysis are classified, compared and evaluated for efficiency and accuracy. Proper eigenvector normalization is discussed. An improved method for calculating derivatives of eigenvectors is proposed based on a more rational normalization condition and taking advantage of matrix sparsity. Important numerical aspects of this method are also discussed. To alleviate the problem of reanalysis, various approximation methods for eigenvalues are proposed and evaluated. Linear and quadratic approximations are based directly on the Taylor series. Several approximation methods are developed based on the generalized Rayleigh quotient for the eigenvalue problem. Approximation methods based on trace theorem give high accuracy without needing any derivatives. Operation counts for the computation of the approximations are given. General recommendations are made for the selection of appropriate approximation technique as a function of the matrix size, number of design variables, number of eigenvalues of interest and the number of design points at which approximation is sought.

Murthy, D. V.↗

NASA NTRS1988

Prediction of the microyield strength of polymer matrix composites

A theory for the prediction of the creep and the microyield strength of composite laminates with viscoelastic matrices has been developed and tested. Precise creep curves from three laminates with arbitrary orientations provide the necessary input data. The time-dependent lamina properties are derived for subsequent use with standard laminate codes. Both the level and the shape of creep and strain recovery curves can be predicted. Extensions of the theory to include the effects of microcracking are outlined.

Wolff, Ernest G.↗

NASA NTRS1988

System identification using modal residues and roots

This presentation will address an approach for using modal residues and roots, typically derived from vibration testing, to solve for the lumped parameters of a finite element model. The uniqueness of the approach is that the mass, stiffness, or damping matrices of a structure do not need to be known or estimated a-priori. Instead just the connectivity of the structure needs to be estimated. In addition, the presentation will demonstrate how models with orders larger than the number of modes obtained by the test program can be solved. The presentation will start by covering the basis of the approach, then illustrate a trivial application of how it works, and finally illustrate its use on a cantilevered beam using only a fraction of the mode to solve the mass, stiffness, and damping matrices.

Wendler, Bruce H.↗

NASA NTRS1994

Weighted graph based ordering techniques for preconditioned conjugate gradient methods

We describe the basis of a matrix ordering heuristic for improving the incomplete factorization used in preconditioned conjugate gradient techniques applied to anisotropic PDE's. Several new matrix ordering techniques, derived from well-known algorithms in combinatorial graph theory, which attempt to implement this heuristic, are described. These ordering techniques are tested against a number of matrices arising from linear anisotropic PDE's, and compared with other matrix ordering techniques. A variation of RCM is shown to generally improve the quality of incomplete factorization preconditioners.

Clift, Simon S.↗

NASA NTRS1993

Polarimetric radar data decomposition and interpretation

Significant efforts have been made to decompose polarimetric radar data into several simple scattering components. The components which are selected because of their physical significance can be used to classify SAR (Synthetic Aperture Radar) image data. If particular components can be related to forest parameters, inversion procedures may be developed to estimate these parameters from the scattering components. Several methods have been used to decompose an averaged Stoke's matrix or covariance matrix into three components representing odd (surface), even (double-bounce) and diffuse (volume) scatterings. With these decomposition techniques, phenomena, such as canopy-ground interactions, randomness of orientation, and size of scatters can be examined from SAR data. In this study we applied the method recently reported by van Zyl (1992) to decompose averaged backscattering covariance matrices extracted from JPL SAR images over forest stands in Maine, USA. These stands are mostly mixed stands of coniferous and deciduous trees. Biomass data have been derived from field measurements of DBH and tree density using allometric equations. The interpretation of the decompositions and relationships with measured stand biomass are presented in this paper.

Sun, Guoqing↗

NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

NASA NTRS2002

The Delta x B = 0 Constraint Versus Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. It has previously been shown by Powell et al. (1995) that an 'almost' equivalent MHD system in non-conservative form can be derived. This non-conservative system has a better conditioned eigensystem. Aside from Powell et al., the MHD equations can be derived from basic principles in either conservative or non-conservative form. The Delta x B = 0 constraint of the MHD equations is only an initial condition constraint, it is very different from the incompressible Navier-Stokes equations in which the divergence condition is needed to close the system (i.e., to have the same number of equations and the same number of unknown). In the MHD formulations, if Delta x B = 0 initially, all one needs is to construct appropriate numerical schemes that preserve this constraint at later time evolutions. In other words, one does not need the Delta x B condition to close the MHD system. We formulate our new scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the divergence of the magnetic condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗

NASA NTRS1997

Bromine Dioxide, OBrO: The Rotational Spectrum and Molecular Properties

The spectra are well descrie by Hamiltonian which included centrifugal distortion efects for fine and hyperfine terms. The molecular structure and the harmonic force field have been derived, and they as well as fine and hyperfine structure constants, are compared with data of related molecules and ESR data from OBrO isolated in cryogenic salt matrices.

Bromine Dioxide Rotational Spectrum↗

NASA NTRS2001

A Curved, Elastostatic Boundary Element for Plane Anisotropic Structures

The plane-stress equations of linear elasticity are used in conjunction with those of the boundary element method to develop a novel curved, quadratic boundary element applicable to structures composed of anisotropic materials in a state of plane stress or plane strain. The curved boundary element is developed to solve two-dimensional, elastostatic problems of arbitrary shape, connectivity, and material type. As a result of the anisotropy, complex variables are employed in the fundamental solution derivations for a concentrated unit-magnitude force in an infinite elastic anisotropic medium. Once known, the fundamental solutions are evaluated numerically by using the known displacement and traction boundary values in an integral formulation with Gaussian quadrature. All the integral equations of the boundary element method are evaluated using one of two methods: either regular Gaussian quadrature or a combination of regular and logarithmic Gaussian quadrature. The regular Gaussian quadrature is used to evaluate most of the integrals along the boundary, and the combined scheme is employed for integrals that are singular. Individual element contributions are assembled into the global matrices of the standard boundary element method, manipulated to form a system of linear equations, and the resulting system is solved. The interior displacements and stresses are found through a separate set of auxiliary equations that are derived using an Airy-type stress function in terms of complex variables. The capabilities and accuracy of this method are demonstrated for a laminated-composite plate with a central, elliptical cutout that is subjected to uniform tension along one of the straight edges of the plate. Comparison of the boundary element results for this problem with corresponding results from an analytical model show a difference of less than 1%.

Smeltzer, Stanley S.↗

NASA NTRS2004

Microcrystals and Amorphous Material in Comets and Primitive Meteorites: Keys to Understanding Processes in the Early Solar System

Comets, fine-grained matrices of chondrites, and chondritic interplanetary dust particles (IDPs) are each composed of both crystalline and amorphous silicates. The primitive solar nebula, in which comets and asteroids accreted, was formed from the collapsed core of a Giant Molecular Cloud, that, in turn, condensed from materials present in the interstellar medium (ISM). Despite observations that reveal the presence of crystalline magnesium silicate minerals in the shells of very high mass-loss-rate stars [1,2], typical silicate grains in the ISM are most likely to be amorphous, given their relatively long residence time in such a high radiation environment. An upper limit of ~3% crystalline grains can be derived from their non-detection in spectra of ISM solids [3]. If the vast majority of grains that enter the primitive solar nebula are amorphous, then the observation of crystalline dust in comets and primitive chondrite matrices indicates the action of specific processes required to transform the amorphous starting materials into the crystals that are observed.

Nuth, J. A.↗

DOE OSTI2022

Extracting single fiber transverse and shear moduli from off-axis misalignment fiber tensile testing

Small diameter (<100 μm) fibers (e.g. carbon fibers, Kevlar, and fiberglass) and wires (e.g. ultrafine copper and aluminum wires) are frequently used in many different engineering applications, such as for light weighting structures, electromagnetic shielding for aircraft/infrastructure/EVs, vibration damping, biological sensors, aerospace electrical devices, and electric windings just to name a few. Due to the manufacturing process, the fibers and wires are pulled and stretched to produce a preferential alignment. Therefore, thin fibers and wires typically display different properties along the length of the fiber as opposed to their cross section and many fibers/wires are considered transversely isotropic. The axial properties of fibers/wires can be ascertained via tensile testing of single-filaments or fiber tows, but the radial properties require much more effort to measure. Knowing these properties is important for the accurate prediction of micromechanical models and manipulation of fibers during micromanufacturing. In this paper, a new technique was developed to determine the transverse/shear moduli and strength of a material by conducting tensile tests of the material at increasing misalignment angles from the tensile axis. Due to the transversely isotropic nature of the material, the transverse/shear moduli and strength influence the experimental results recorded by the test machine to different degrees based on the amount of misalignment in the test setup. An equation was derived to determine the influence of each of the material properties based on the misalignment angle by manipulating the stiffness matrix for transversely isotropic materials using the transformation matrices. Then, curve fitted coefficients were used to identify the material properties. Here, the strengths were similarly determined by curve fitting an off-axis Tsai-Hill failure criteria to determine the influence of transverse, shear, and tensile strengths based on the complex loading condition provided by the off-axis tensile test. Zoltek Panex 35 carbon fibers were used to demonstrate this new technique and the determined properties were then compared to those obtained from nanoindentation and from literature. Fracture surfaces provide insight into the different failure mechanisms at various misalignment angles.

36 MATERIALS SCIENCE↗

NASA NTRS1975

Problems, questions and results in the use of the BBN-model

A hierarchical structured control model was used to describe pilot performance in a multivariable control problem. The model is based on the assumption that a well trained operator optimized a quadratic optimization criterion. The model considers some special characteristics of the human operator as time delay, neuromuscular time lag, observation and motor noise and the ability to extract the first derivative of a displaced value. Numerical results of a simple human engineering example using the BBN model demonstrate the influence of the variation of different parameters on covariance matrices of the state and observation vectors.

Dey, D.↗

NASA NTRS1987

Ceramic matrix and resin matrix composites: A comparison

The underlying theory of continuous fiber reinforcement of ceramic matrix and resin matrix composites, their fabrication, microstructure, physical and mechanical properties are contrasted. The growing use of organometallic polymers as precursors to ceramic matrices is discussed as a means of providing low temperature processing capability without the fiber degradation encountered with more conventional ceramic processing techniques. Examples of ceramic matrix composites derived from particulate-filled, high char yield polymers and silsesquioxane precursors are provided.

Hurwitz, Frances I.↗

NASA NTRS1986

The use of multidimensional scaling for facilities layout - An application to the design of the Space Station

Before an optimum layout of the facilities for the proposed Space Station can be designed, it is necessary to understand the functions that will be performed by the Space Station crew and the relationships among those functions. Five criteria for assessing functional relationships were identified. For each of these criteria, a matrix representing the degree of association of all pairs of functions was developed. The key to making inferences about the layout of the Space Station from these matrices was the use of multidimensional scaling (MDS). Applying MDS to these matrices resulted in spatial configurations of the crew functions in which smaller distances in the MDS configuration reflected closer associations. An MDS analysis of a composite matrix formed by combining the five individual matrices resulted in two dimensions that describe the configuration: a 'private-public' dimension and a 'group-individual' dimension. Seven specific recommendations for Space Station layout were derived from analyses of the MDS configurations. Although these techniques have been applied to the design of the Space Station, they can be applied to the design of any facility where people live or work.

Tullis, Thomas S.↗

NASA NTRS1987

Ceramic matrix and resin matrix composites - A comparison

Hurwitz, Frances I.↗

NASA NTRS1990

Equations of motion of slung load systems with results for dual lift

General simulation equations are derived for the rigid body motion of slung load systems. These systems are viewed as consisting of several rigid bodies connected by straight-line cables or links. The suspension can be assumed to be elastic or inelastic, both cases being of interest in simulation and control studies. Equations for the general system are obtained via D'Alembert's principle and the introduction of generalized velocity coordinates. Three forms are obtained. Two of these generalize previous case-specific results for single helicopter systems with elastic or inelastic suspensions. The third is a new formulation for inelastic suspensions. It is derived from the elastic suspension equations by choosing the generalized coordinates so as to separate motion due to cable stretching from motion with invariant cable lengths. The result is computationally more efficient than the conventional formulation, and is readily integrated with the elastic suspension formulation and readily applied to the complex dual lift and multilift systems. Equations are derived for dual lift systems. Three proposed suspension arrangements can be integrated in a single equation set. The equations are given in terms of the natural vectors and matrices of three-dimensional rigid body mechanics and are tractable for both analysis and programming.

Cicolani, Luigi S.↗