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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Parameter identifiability of linear dynamical systems

It is assumed that the system matrices of a stationary linear dynamical system were parametrized by a set of unknown parameters. The question considered here is, when can such a set of unknown parameters be identified from the observed data? Conditions for the local identifiability of a parametrization are derived in three situations: (1) when input/output observations are made, (2) when there exists an unknown feedback matrix in the system and (3) when the system is assumed to be driven by white noise and only output observations are made. Also a sufficient condition for global identifiability is derived.

Glover, K.↗

NASA NTRS1977

A special purpose symbolic manipulation program for the derivation of the equations of motion for large flexible structures

This paper outlines a computer program especially tailored to the task of deriving explicit equations of motion for structures with point-connected substructures. The special purpose program is written in FORTRAN and is designed for performing the specific algebraic operations encountered in the derivation of explicit equations of motion. The derivation is by the Lagrangian approach. Using an orderly kinematical procedure and a discretization and/or truncation scheme, it is possible to write the kinetic and potential energy of each substructure in a compact vector-matrix form. Then, if each element of the matrices and vectors encountered in the kinetic and potential energy is a known algebraic expression, the computer program performs the necessary operations to evaluate the kinetic and potential energy of the system explicitly. Lagrange's equations for small motions about equilibrium can be deduced directly from the explicit form of the system kinetic and potential energy.

Hale, A. L.↗

NASA NTRS1978

Derivation of the equations of motion for complex structures by symbolic manipulation

Hale, A. L.↗

NASA NTRS1990

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗

NASA NTRS2002

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Freed, Alan D.↗

DOE OSTI2024

Evaluation of a practical approach for field scale moisture flow modeling in heterogeneous media at a semiarid site

Abstract A practical approach for modeling field‐scale moisture flow in a highly heterogeneous unsaturated medium is described in this study. The validity of this approach is demonstrated through comparison of the numerical simulations with field observations at a semiarid site located in southcentral Washington State. The methodology is based on upscaling the core scale hydraulic properties and combining power‐law and tensorial connectivity‐tortuosity (PA‐TCT) approaches to derive macroscopic anisotropy parameters for each hydrostratigraphic unit (HSU) identified in the field. Each heterogeneous HSU is approximated by an equivalent homogeneous medium (EHM) model for which PA‐TCT parameters are used in the flow simulations. The available field data on moisture content and matric potential are compared with steady‐state flow simulations based on the mean form of Richards' equation. While the homogenization or averaging of heterogeneities, embedded in the EHM modeling approximation, cannot capture all of the field‐scale variability, the simulated steady‐state moisture and matric potential profiles capture well the central tendency of the field data. This approach is deemed practical for assessing the fate and transport of contaminants in highly heterogeneous unsaturated media at the transport scale of hundreds of meters.

Khaleel, Raziuddin↗

DOE OSTI2022

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗

NASA NTRS1988

Dynamic structural correlation via nonlinear programming techniques

A solution to the correlation between structural dynamic test results and finite element analyses of the same components is presented in this paper. Basically, the method can be categorized as a Levenberg-Marquardt type Gauss-Newton method which requires only the differences between FE modal analyses and test results and their first derivatives with respect to preassigned design variables. With proper variable normalization and equation scaling, the method has been made numerically better-conditioned and the inclusion of the Levenberg-Marquardt technique overcomes any remaining difficulty encountered in inverting singular or near-singular matrices. An important feature is that each iteration requires only one function evaluation along with the associated design sensitivity analysis and so the procedure is computationally efficient.

Ting, T.↗

NASA NTRS2020

Shape Estimation for Elongated Deformable Object using B-spline Chained Multiple Random Matrices Model

In this paper, a B-spline chained multiple random matrix models (RMMs) representation is proposed to model geometric characteristics of an elongated deformable object. The hyper degrees of freedom structure of the elongated deformable object make its shape estimation challenging. Based on the likelihood function of the proposed B-spline chained multiple RMMs, an expectation-maximization (EM) method is derived to estimate the shape of the elongated deformable object. A split and merge method based on the Euclidean minimum spanning tree (EMST) is proposed to provide initialization for the EM algorithm. The proposed algorithm is evaluated for the shape estimation of the elongated deformable objects in scenarios, such as the static rope with various configurations (including configurations with intersection), the continuous manipulation of a rope and a plastic tube, and the assembly of two plastic tubes. The execution time is computed and the accuracy of the shape estimation results is evaluated based on the comparisons between the estimated width values and its ground-truth, and the intersection over union (IoU) metric.

Gang Yao↗

NASA NTRS1973

Structural Aspects of System Identification

The problem of identifying linear dynamical systems is studied by considering structural and deterministic properties of linear systems that have an impact on stochastic identification algorithms. In particular considered is parametrization of linear systems so that there is a unique solution and all systems in appropriate class can be represented. It is assumed that a parametrization of system matrices has been established from a priori knowledge of the system, and the question is considered of when the unknown parameters of this system can be identified from input/output observations. It is assumed that the transfer function can be asymptotically identified, and the conditions are derived for the local, global and partial identifiability of the parametrization. Then it is shown that, with the right formulation, identifiability in the presence of feedback can be treated in the same way. Similarly the identifiability of parametrizations of systems driven by unobserved white noise is considered using the results from the theory of spectral factorization.

Glover, Keith↗

NASA NTRS1988

The eddy transport of nonconserved trace species derived from satellite data

Using the approach of the Garcia and Solomon (1983) model and data obtained by the LIMS instrument on Nimbus 7, the chemical eddy transport matrix for planetary waves was calculated, and the chemical eddy contribution to the components of the matrix obtained from the LIMS satellite observations was computed using specified photochemical damping time scales. The dominant component of the transport matrices for several winter months were obtained for ozone, nitric acid, and quasi-geostrophic potential vorticity (PV), and the parameterized transports of these were compared with the 'exact' transports, computed directly from the eddy LIMS data. The results indicate that the chemical eddy effect can account for most of the observed ozone transport in early winter, decreasing to less than half in late winter. The agreement between the parameterized and observed nitric acid and PV was not as good. Reasons for this are discussed.

Smith, Anne K.↗

NASA NTRS1991

Second-order discrete Kalman filtering equations for control-structure interaction simulations

A general form for the first-order representation of the continuous, second-order linear structural dynamics equations is introduced in order to derive a corresponding form of first-order Kalman filtering equations (KFE). Time integration of the resulting first-order KFE is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete KFE involving only symmetric, N x N solution matrix.

Park, K. C.↗

NASA NTRS1999

A Textbook for a First Course in Computational Fluid Dynamics

This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Zingg, D. W.↗

DOE OSTI2024

A Bayesian Learning Approach to Wireless Outdoor Heatmap Construction using Deep Gaussian Process

We present a novel Bayesian learning approach to outdoor radio heatmap construction utilizing deep Gaussian process (GP). The proposed approach employs a two-layer hierarchy which consists of two cascaded Gaussian processes that are capable of modeling more complex input-output relations than standard single-layer Gaussian processes. Since deriving the exact model likelihood is challenging, a lower bound is optimized instead so that gradient descent-based methods can be performed to find out the optimal model parameters. Typically, inducing points are used in GPs to facilitate low-rank approximation of covariance (kernel) matrices for computation speedup. However, the inaccuracy induced by inducing points can accumulate when stacking multiple layers of GP which may hinder the performance of deep GP. Moreover, since inducing points need to be learned, having them at all layers of deep GP also incurs computational burden. To overcome the above challenges, in contrast to the canonical deep GP model, we use a modified architecture where a full standard GP resides in the first layer and inducing points are only introduced for the second layer. This modified architecture strikes a balance between model accuracy and training complexity. In the proposed model, the noise parameter of the first GP layer is also eliminated to improve the training efficiency as the noise parameter at the output of the second layer suffices to model the uncertainty in the output. The proposed approach is evaluated on real-world datasets, in the form of location-Received Signal Strength (RSS) pairs, collected from the Platform for Open Wireless Data-driven Experimental Research (POWDER) located at the campus of the University of Utah. Experiment results show that the proposed approach can achieve smaller prediction errors on various training and testing data configurations than DNN-based and GP-based methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Enhanced purification coupled with biophysical analyses shows cross-β structure as a core building block for Streptococcus mutans functional amyloids

Streptococcus mutans is an etiologic agent of human dental caries that forms dental plaque biofilms containing functional amyloids. Three amyloidogenic proteins, P1, WapA, and Smu_63c were previously identified. C123 and AgA are naturally occurring amyloid-forming fragments of P1 and WapA, respectively. We determined that four amyloidophilic dyes, ThT, CDy11, BD-oligo, and MK-H4, differentiate C123, AgA, and Smu_63c amyloid from monomers, but non-specific binding to bacterial cells in the absence of amyloid precludes their utility for identifying amyloid in biofilms. Congo redinduced birefringence is a more specific indicator of amyloid formation and differentiates biofilms formed by wild-type S. mutans from a triple ΔP1/WapA/Smu_63c mutant with reduced bioiflm forming capabilities. Amyloid accumulation is a late event, appearing in older S. mutans biofilms after 60hours of growth. Amyloid derived from pure preparations of all three proteins is visualized by electron microscopy as mat-like structures. Typical amyloid fibers become evident following protease digestion to eliminate non-specific aggregates and monomers. Amyloid mats, similar in appearance to those reported in S. mutans biofilm extracellular matrices, are reconstituted by co-incubation of monomers and amyloid fibers. X-ray fiber diffraction of amyloid mats and fibers from all three proteins demonstrate patterns reflective of a cross-β amyloid structure.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1981

An influence coefficient method for the application of the modal technique to wing flutter suppression of the DAST ARW-1 wing

The methods used to compute the mass, structural stiffness, and aerodynamic forces in the form of influence coefficient matrices as applied to a flutter analysis of the Drones for Aerodynamic and Structural Testing (DAST) Aeroelastic Research Wing. The DAST wing was chosen because wind tunnel flutter test data and zero speed vibration data of the modes and frequencies exist and are available for comparison. A derivation of the equations of motion that can be used to apply the modal method for flutter suppression is included. A comparison of the open loop flutter predictions with both wind tunnel data and other analytical methods is presented.

Pines, S.↗

NASA NTRS1983

Identification of general linear mechanical systems

Previous work in identification theory has been concerned with the general first order time derivative form. Linear mechanical systems, a large and important class, naturally have a second order form. This paper utilizes this additional structural information for the purpose of identification. A realization is obtained from input-output data, and then knowledge of the system input, output, and inertia matrices is used to determine a set of linear equations whereby we identify the remaining unknown system matrices. Necessary and sufficient conditions on the number, type and placement of sensors and actuators are given which guarantee identificability, and less stringent conditions are given which guarantee generic identifiability. Both a priori identifiability and a posteriori identifiability are considered, i.e., identifiability being insured prior to obtaining data, and identifiability being assured with a given data set.

Sirlin, S. W.↗

DOE OSTI2020

Effects of He Ion Irradiation in Simple Borosilicate Glasses for HLW Immobilisation - 20324

High-level radioactive waste (HLW) consists of fission products and minor actinides. Borosilicate glass is the host material of choice in most countries to safely immobilize HLW, followed by eventual geological disposal. HLW glasses will experience both self-heating and self-irradiation effects due to the decay of fission products and minor actinides present therein. During the first 500 years of geological disposal, beta-gamma decay will be the main source of radiation due to short half-lives of fission products, as compared to very long-lived actinides which will be responsible for irradiation over much longer timescales. High-energy alpha particles and alpha recoil nuclei can potentially have significant impact (stored energy, volume change, accumulated dose, chemical durability) on the performance of the wasteform over longer timescales. Thus, to ensure safe and economic disposal over geological periods (10{sup 4} - 10{sup 6} years), the assessment of performance and properties in terms of radiation effects becomes indispensable. This work aims to simulate the effects of alpha decay in radioactive waste glasses. Two international glasses, namely UK 'MW' lithium-sodium borosilicate and Indian sodium-barium borosilicate, which are being used as HLW immobilization matrices by the UK and India respectively, were prepared by a melt-quench method and annealed below their glass transition temperatures. Both glasses were irradiated with 650 KeV He{sup 2+} ions with doses ranging from 10{sup 15} to 10{sup 17} ions/cm{sup 2}. Structural analysis was carried out using electron paramagnetic resonance (EPR) and Raman spectroscopies. First derivative continuous wave EPR spectra were recorded for all pristine and irradiated glass monoliths at room temperature using an X-band (frequency 9.6 GHz) EPR Bruker spectrometer. Results from EPR analyses showed the formation of irradiation induced paramagnetic centers at g ∼ 2.00 in all irradiated glasses. Raman spectra of pristine and irradiated samples were acquired, and showed notable changes in glass structure post-irradiation, with the emergence of a broad Raman band at 1560 cm{sup -1} upon irradiation with high fluences. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗