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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A reliable algorithm for optimal control synthesis

In recent years, powerful design tools for linear time-invariant multivariable control systems have been developed based on direct parameter optimization. In this report, an algorithm for reliable optimal control synthesis using parameter optimization is presented. Specifically, a robust numerical algorithm is developed for the evaluation of the H(sup 2)-like cost functional and its gradients with respect to the controller design parameters. The method is specifically designed to handle defective degenerate systems and is based on the well-known Pade series approximation of the matrix exponential. Numerical test problems in control synthesis for simple mechanical systems and for a flexible structure with densely packed modes illustrate positively the reliability of this method when compared to a method based on diagonalization. Several types of cost functions have been considered: a cost function for robust control consisting of a linear combination of quadratic objectives for deterministic and random disturbances, and one representing an upper bound on the quadratic objective for worst case initial conditions. Finally, a framework for multivariable control synthesis has been developed combining the concept of closed-loop transfer recovery with numerical parameter optimization. The procedure enables designers to synthesize not only observer-based controllers but also controllers of arbitrary order and structure. Numerical design solutions rely heavily on the robust algorithm due to the high order of the synthesis model and the presence of near-overlapping modes. The design approach is successfully applied to the design of a high-bandwidth control system for a rotorcraft.

Vansteenwyk, Brett↗
NASA NTRS1994

A fiber matrix model for fluid flow and streaming potentials in the canaliculi of an osteon

A theoretical model is developed to predict the fluid shear stress and streaming potential at the surface of osteocytic processes in the lacunar-canalicular porosity of an osteon when the osteon is subject to mechanical loads that are parallel or perpendicular to its axis. The theory developed in Weinbaum et al. (31) for the flow through a proteoglycan matrix in a canaliculus is employed in a poroelastic model for the osteon. Our formulation is a generalization of that of Petrov et al. (17). Our model predicts that, in order to satisfy the measured frequency dependence of the phase and magnitude of the SGP in macroscopic bone samples, the fiber spacing in the fluid annulus must lie in the narrow range 6-7 nm typical of the spacing of GAG sidechains along a protein monomer. The model predictions for the local SGP profiles in the osteon agree with the experimental observations of Starkebaum et al. (24). The theory predicts that the pore pressure relaxation time, tau d, for a 150-300 microns diameter osteon with the foregoing matrix structure is approximately 0.03-0.13 sec, and that the amplitude of the mean fluid shear stress on the membrane of the osteocytic process at the mean areal radius of the osteon has a maximum at 28 Hz if tau d = 0.06 sec. This maximum, which is independent of the magnitude of the loading, could be important in vivo since the recent experiments of Turner et al. (28) and McLeod et al. (15) have a peak in the strain frequency spectrum between 20 and 30 Hz that also appears to be independent of the type (magnitude) of loading. Numerical predictions for the amplitude of the average fluid shear stress on the osteocytic membrane at the mean areal radius of the osteon show that the fluid shear stress associated with the low amplitude 20-30 Hz spectral strain component is at least as large as the average fluid shear stress associated with the high amplitude 1 Hz stride component, although the latter loading is an order of magnitude larger, and has a magnitude that lies within the middle of the range, 6-30 dynes/cm2, where fluid shear stresses in tissue culture studies with osteoblast monolayers have elicited an intracellular Ca++ response (31). The implications of these results for intracellular electrical communication are discussed.

Non-NASA Center↗
NASA NTRS2010

Guidance and Control System for a Satellite Constellation

A distributed guidance and control algorithm was developed for a constellation of satellites. The system repositions satellites as required, regulates satellites to desired orbits, and prevents collisions. 1. Optimal methods are used to compute nominal transfers from orbit to orbit. 2. Satellites are regulated to maintain the desired orbits once the transfers are complete. 3. A simulator is used to predict potential collisions or near-misses. 4. Each satellite computes perturbations to its controls so as to increase any unacceptable distances of nearest approach to other objects. a. The avoidance problem is recast in a distributed and locally-linear form to arrive at a tractable solution. b. Plant matrix values are approximated via simulation at each time step. c. The Linear Quadratic Gaussian (LQG) method is used to compute perturbations to the controls that will result in increased miss distances. 5. Once all danger is passed, the satellites return to their original orbits, all the while avoiding each other as above. 6. The delta-Vs are reasonable. The controller begins maneuvers as soon as practical to minimize delta-V. 7. Despite the inclusion of trajectory simulations within the control loop, the algorithm is sufficiently fast for available satellite computer hardware. 8. The required measurement accuracies are within the capabilities of modern inertial measurement devices and modern positioning devices.

Bryson, Jonathan Lamar↗
NASA NTRS2011

A Sensitivity Study on the Effects of Particle Chemistry, Asphericity and Size on the Mass Extinction Efficiency of Mineral Dust in the Earth's Atmosphere: From the Near to Thermal IR

To determine a plausible range of mass extinction efficiencies (MEE) of terrestrial atmospheric dust from the near to thermal IR, sensitivity analyses are performed over an extended range of dust microphysical and chemistry perturbations. The IR values are subsequently compared to those in the near-IR, to evaluate spectral relationships in their optical properties. Synthesized size distributions consistent with measurements, model particle size, while composition is defined by the refractive indices of minerals routinely observed in dust, including the widely used OPAC/Hess parameterization. Single-scattering properties of representative dust particle shapes are calculated using the T-matrix, Discrete Dipole Approximation and Lorenz-Mie light-scattering codes. For the parameterizations examined, MEE ranges from nearly zero to 1.2 square meters per gram, with the higher values associated with non-spheres composed of quartz and gypsum. At near-IR wavelengths, MEE for non-spheres generally exceeds those for spheres, while in the thermal IR, shape-induced changes in MEE strongly depend on volume median diameter (VMD) and wavelength, particularly for MEE evaluated at the mineral resonant frequencies. MEE spectral distributions appear to follow particle geometry and are evidence for shape dependency in the optical properties. It is also shown that non-spheres best reproduce the positions of prominent absorption peaks found in silicates. Generally, angular particles exhibit wider and more symmetric MEE spectral distribution patterns from 8-10 micrometers than those with smooth surfaces, likely due to their edge-effects. Lastly, MEE ratios allow for inferring dust optical properties across the visible-IR spectrum. We conclude the MEE of dust aerosol are significant for the parameter space investigated, and are a key component for remote sensing applications and the study of direct aerosol radiative effects.

Hansell, R. A., Jr.↗
NASA NTRS1999

High Temperature Mechanical Characterization and Analysis of Al2O3 /Al2O3 Composition

Sixteen ply unidirectional zirconia coated single crystal Al2O3 fiber reinforced polycrystalline Al2O3 was tested in uniaxial tension at temperatures to 1400 C in air. Fiber volume fractions ranged from 26 to 31%. The matrix has primarily open porosity of approximately 40%. Theories for predicting the Young's modulus, first matrix cracking stress, and ultimate strength were applied and evaluated for suitability in predicting the mechanical behavior of Al2O3/Al2O3 composites. The composite exhibited pseudo tough behavior (increased area under the stress/strain curve relative to monolithic alumina) from 22 to 1400 C. The rule-of-mixtures provides a good estimate of the Young's modulus of the composite using the constituent properties from room temperature to approximately 1200 C for short term static tensile tests in air. The ACK theory provides the best approximation of the first matrix cracking stress while accounting for residual stresses at room temperature. Difficulties in determining the fiber/matrix interfacial shear stress at high temperatures prevented the accurate prediction of the first matrix cracking stress above room temperature. The theory of Cao and Thouless, based on Weibull statistics, gave the best prediction for the composite ultimate tensile strength.

Gyekenyesi, John Z.↗
NASA NTRS1989

Simultaneous quaternion estimation (QUEST) and bias determination

Tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function are presented. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and Earth-pointing reference attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. They converge to the same estimates at the same rate for inertially-fixed attitude, but the new algorithm converges more slowly than the differential correction for Earth-pointing attitude. The slower convergence of the new method for non-zero attitude rates is believed to be due to the use of an inadequate approximation for a partial derivative matrix. The new method requires about twice the computational effort of the differential correction. Improving the approximation for the partial derivative matrix in the new method is expected to improve its convergence at the cost of increased computational effort.

Markley, F. Landis↗
NASA NTRS1967

Renormalized Born approximation.

Renormalizing approximate wave functions so that amplitude is correct by means of matrix, with applications to Born approximation

BORN APPROXIMATION↗
DOE OSTI2023

Aggregation of Leslie matrix models with application to ten diverse animal species

Some grouping is necessary when constructing a Leslie matrix model because it involves discretizing a continuous process of births and deaths. Here, the level of grouping is determined by the number of age classes and frequency of sampling. It is largely unknown what is lost or gained by using fewer age classes, and I address this question using aggregation theory. I derive an aggregator for a Leslie matrix model using weighted least squares, determine what properties an aggregated matrix inherits from the original matrix, evaluate aggregation error, and measure the influence of aggregation on asymptotic and transient behaviors. To gauge transient dynamics, I employ reactivity of the standardized Leslie matrix. I apply the aggregator to 10 Leslie models developed for animal populations drawn from a diverse set of species. Several properties are inherited by the aggregated matrix: (a) it is a Leslie matrix; (b) it is irreducible whenever the original matrix is irreducible; (c) it is primitive whenever the original matrix is primitive; and (d) its stable population growth rate and stable age distribution are consistent with those of the original matrix if the least squares weights are equal to the original stable age distribution. In the application, depending on the population modeled, when the least squares weights do not follow the stable age distribution, the stable population growth rate of the aggregated matrix may or may not be approximately consistent with that of the original matrix. Transient behavior is lost with high aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

140,142 Ce Neutron Cross Section Resolved Resonance Region Evaluation

A resolved resonance region evaluation of 140,142 Ce was conducted by Oak Ridge National Laboratory. Requested by the US Nuclear Criticality Safety Program, this evaluation is based on recent high-resolution transmission and capture high-resolution measurements of nat Ce and 142 Ce conducted at JRC-Geel at the Geel Linear Accelerator facility. It is also based on recently measured thermal constants available from the EX FOR database. Starting from the resonance parameters from the ENDF/B-VIII.0 library and following a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance in formation was derived by the fit of these experimental datasets using the Reich–Moore approximation of the R-matrix theory, as implemented in the SAMMY code system. The resolved resonance region upper energy limit for 140 Ce was kept at 200 keV, whereas the 142 Ce resonance region was extended from 13 to 26 keV. This new evaluation was found to be in good agreement not only with several integral quantities of interest to the reactor physics community, but also with the stellar Maxwellian-averaged cross section.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1995

Metamorphism and aqueous alteration in low petrographic type ordinary chondrites

In order to investigate the relative importance of dry metamorphism and aqueous alteration in the history of chondruies, chondruies were hand-picked from the Semarkona (petrographic type 3.0), Bishunpur (3. 1), Chainpur (3.4), Dhajala (3.8) and Allegan (5) chondrites, and matrix samples were extracted from the first three ordinary chondrites. The thermoluminescence (TL) properties of all the samples were measured, and appropriate subsets of the samples were analyzed by electron-microprobe and radiochemical neutron activation and the water and H-isotopic composition determined. The TL data for chondrules from Semarkona and Bishunpur scatter widely showing no unambiguous trends, although group B1 chondrules tend to have lower sensitivities and lower peak temperatures compared with group A5 chondrules. It is argued that these data reflect the variety of processes accompanying chondrule formation. The chondrules show remarkably uniform contents of the highly labile elements, indicating mineralogical control on abundance and volatile loss from silicates and loss and recondensation of mobile chalcophiles and siderophiles in some cases. Very high D/H values (up to approx. 8000% SMOW) are observed in certain Semarkona chondrules, a confirmation of earlier work. With increasing petrographic type, mean TL sensitivities of the chondrules increase, the spread of values within an individual meteorite decreases, and peak temperatures and peak widths show trends indicating that the TL is mainly produced by feldspar and that dry, thermal metamorphism is the dominant secondary process experienced by the chondrules. The TL sensitivities of matrix samples also increase with petrographic type. Chainpur matrix samples show the same spread of peak temperatures and peak widths as Chainpur chondruies, indicating metamorphism-related changes in the feldspar are responsible for the TL of the matrix. The TL data for the Semarkona and Bishunpur matrix samples provide, at best, only weak evidence for aqueous alteration, but the matrix contains H with approximately terrestrial D/H values, even though it contains much water. Secondary processes (probably aqueous alteration) presumably lowered the D/H of the matrix and certain chondrules. While chondrule properties appear to be governed primarily by formation processes and subsequent metamorphism, the matrix of Semarkona has a more complex history involving aqueous alteration as a meteorite-wide process.

Xie, T.↗
DOE OSTI2023

R-matrix Resolved Resonance Region Evaluation of 140,142 Ce

Oak Ridge National Laboratory completed the resolved resonance region (RRR) evaluation of the two most abundant cerium isotopes, 140 Ce (88.45%) and 142 Ce (11.11%), as requested by the US Nuclear Criticality Safety Program. These evaluations are based on recent high-resolution transmission and capture measurements performed on nat Ce and highly enriched 142 Ce samples at the JRC-Geel Linear Accelerator facility, as well as measured thermal constants available from the EXFOR database. Starting from the resonance parameters of the ENDF/B-VIII.0 library followed by a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance information were derived by fitting these measured data using the Reich–Moore approximation of the R -matrix theory, as implemented in the SAMMY code system. The 140 Ce RRR upper energy limit was kept at 200 keV, whereas the 142 Ce resonance region was extended from 13 to 26 keV. Updated statistical properties were obtained for the new evaluations and compared to those derived from the ENDF/B-VIII.0 nuclear data library. The new evaluation work improved some of the discrepancies found in previous work, such as the capture resonance integral and stellar Maxwellian-averaged cross sections. These integral quantities were mainly derived from the fit of the latest measured data, especially the neutron capture yield data for 142 Ce isotope.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures

An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.

Gibson, J. S.↗
NASA NTRS1988

An invariant asymptotic formula for solutions of second-order linear ODE's

An invariant-matrix technique for the approximate solution of second-order ordinary differential equations (ODEs) of form y-double-prime = phi(x)y is developed analytically and demonstrated. A set of linear transformations for the companion matrix differential system is proposed; the diagonalization procedure employed in the final stage of the asymptotic decomposition is explained; and a scalar formulation of solutions for the ODEs is obtained. Several typical ODEs are analyzed, and it is shown that the Liouville-Green or WKB approximation is a special case of the present formula, which provides an approximation which is valid for the entire interval (0, infinity).

Gingold, H.↗
DOE OSTI2022

Multirate Exponential Rosenbrock Methods

In this paper we propose a novel class of methods for high-order accurate integration of multirate systems of ordinary differential equation initial-value problems. The proposed methods construct multirate schemes by approximating the action of matrix φ functions within explicit exponential Rosenbrock (ExpRB) methods, thereby called multirate ExpRB (MERB) methods. They consist of the solution to a sequence of modified “fast” initial-value problems, which may themselves be approximated through subcycling any desired initial-value problem solver. In addition to proving how to construct MERB methods from certain classes of ExpRB methods, we provide rigorous convergence analysis of these methods and derive efficient MERB schemes of orders 2 through 6 (the highest-order infinitesimal multirate methods to date). Lastly, we then present numerical simulations to confirm these theoretical convergence rates and to compare the efficiency of MERB methods against other recently introduced high-order multirate methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A new 181 Ta neutron resolved resonance region evaluation

A new 181 Ta neutron resolved resonance region evaluation has been performed from the thermal energy range up to approximately 2.5 keV. The R-matrix SAMMY code was used with the Reich–Moore approximation to evaluate resonance parameters from several experimental data sets. A Monte Carlo approach was used for resonance spin assignments and generating 59 small fictitious resonance levels which were shown to improve the cumulative level, Porter-Thomas, and Wigner distributions as compared to theoretical predictions. Covariance information was also generated for the entire resolved resonance region. Finally, the positive impact of the new evaluation was validated through benchmark calculations which were sensitive to the 181 Ta cross section and showed improvement in the reactivity bias for several benchmark cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2001

Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

NASA Discipline Plant Biology↗