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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Nanocomposite structure of two-line ferrihydrite powder from total scattering

Ferrihydrite is one of the most important iron-containing minerals on Earth. Yet determination of its atomic-scale structure has been frustrated by its intrinsically poor crystallinity. The key difficulty is that physically-different models can appear consistent with the same experimental data. Using X-ray total scattering and a nancomposite reverse Monte Carlo approach, we evaluate the two principal contending models—one a multi-phase system without tetrahedral iron(III), and the other a single phase with tetrahedral iron(III). Our methodology is unique in considering explicitly the complex nanocomposite structure the material adopts: namely, crystalline domains embedded in a poorly-ordered matrix. The multi-phase model requires unphysical structural rearrangements to fit the data, whereas the single-phase model accounts for the data straightforwardly. Hence the latter provides the more accurate description of the short- and intermediate-range order of ferrihydrite. We discuss how this approach might allow experiment-driven (in)validation of complex models for important nanostructured phases beyond ferrihydrite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

RE-INTEGRATE EMT Simulation Software: Graph Convolutional Network for Sparse Matrix Pattern Detection

The increasing complexity of power networks, driven by proliferation of inverters, presents analytical challenges that simplified models often fail to capture, necessitating Electromagnetic Transient (EMT) simulations. EMT models are represented as discretized differential-algebraic equations (DAEs), forming a linear system Ax = b that is computationally intensive to solve. Due to inherent sparsity of adjacency matrix A, distinct patterns emerge that, when accurately identified, enable efficient solver selection to minimize computation time. However, identifying ideal pattern is complicated by numerous reordering algorithms and limited structural insights. To address this, we introduce a Graph Convolutional Network (GCN) model for classifying sparse matrix patterns common in power system analysis. The model, achieving 96% test accuracy, is validated using PV plant models of 125 MW capacities connected to New England 39-bus transmission system (TS), and further scaled to a 4,992-bus network with 384 PV plants, yielding 191, 616 × 191, 616 sized A matrix. For all cases, the GCN model accurately identifies the matrix’s intrinsic sparse pattern, demonstrating its potential to enhance solver performance in EMT analysis.

Hossain, Md Rifat [Florida International Universit↗
DOE OSTI2019

Evolution of Microstructure and Creep Behavior in an Fe-Ni-Cr-Nb-C Alloy during Service in Hydrocarbon Cracker Tubes

Heat resistant, cast austenitic Fe-Ni-Cr-Nb-C alloy tubes are used by the petrochemical industry to transport hydrocarbons at temperatures between 850 and 1050 °C. Exposure to carbon from the process stream and nitrogen from the air causes severe carburization and nitridation in the tube walls, which influences creep resistance and degrades their ductility. A better understanding of these processes would enable their service life to be extended. To this end, we report the results of an experimental and computational study of the changes in microstructure and creep behavior in a representative alloy during service. The microstructure of the fresh material is compared to that of 8- and 10-year-old tubes. The results confirm the presence of M{sub 23}C{sub 6}, M{sub 7}C{sub 3}, and M{sub 2}(CN) phases reported in prior studies, and also reveal the nucleation and growth of the carbo-nitride phase M{sub 6}(CN). We show that the spatial variations of carbide and nitride particle volume fractions in aged tubes can be predicted by a simple reaction–diffusion equation. In addition, the uniaxial compressive creep response of the fresh and aged tubes was obtained as functions of temperature and strain rate. The materials deform by power-law creep, with strain rate versus stress relation as ε-dot =Aσ{sup n}. Aging reduces the stress exponent from n = 8 in fresh tubes to n = 5 in aged tubes, and also reduces the pre-exponential coefficient A. Tests on model alloys with composition identical to that of the matrix combined with finite element simulations of creep in model microstructures reveal that Cr depletion is responsible for the change in stress exponent, while the change in A is caused by the increase in particle volume fraction. By combining the microstructure-scale simulations of creep with the predictions of the reaction–diffusion model of microstructure evolution, the creep behavior of tubes left in service beyond 10 years is estimated.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Hub flexibility effects on propfan vibration

The significance of hub flexibility in the nonlinear static and dynamic analyses of advanced turboprop blades is assessed. The chosen blade is the 0.175 scale model of the GE-A7-B4 unducted fan blade. A procedure for coupling the effective hub stiffness matrix to an MSC/NASTRAN finite element model is defined and verified. A series of nonlinear static and dynamic analyses are conducted on the blade for both rigid and flexible hug configurations. Results indicate that hub flexibility is significant in the nonlinear static and dynamic analyses of the GE-A7-B4. In order to insure accuracy in analyses of other blades, hub flexibility should always be considered.

Ernst, Michael A.↗
NASA NTRS2016

Evaluating Secondary Inorganic Aerosols in Three Dimensions

The spatial distribution of aerosols and their chemical composition dictates whether aerosols have a cooling or a warming effect on the climate system. Hence, properly modeling the three-dimensional distribution of aerosols is a crucial step for coherent climate simulations. Since surface measurement networks only give 2-D data, and most satellites supply integrated column information, it is thus important to integrate aircraft measurements in climate model evaluations. In this study, the vertical distribution of secondary inorganic aerosol (i.e., sulfate, ammonium, and nitrate) is evaluated against a collection of 14 AMS flight campaigns and surface measurements from 2000 to 2010 in the USA and Europe. GISS ModelE2 is used with multiple aerosol microphysics (MATRIX, OMA) and thermodynamic (ISORROPIA II, EQSAM) configurations. Our results show that the MATRIX microphysical scheme improves the model performance for sulfate, but that there is a systematic underestimation of ammonium and nitrate over the USA and Europe in all model configurations. In terms of gaseous precursors, nitric acid concentrations are largely underestimated at the surface while overestimated in the higher levels of the model. Heterogeneous reactions on dust surfaces are an important sink for nitric acid, even high in the troposphere. At high altitudes, nitrate formation is calculated to be ammonia limited. The underestimation of ammonium and nitrate in polluted regions is most likely caused by a too simplified treatment of the NH3/NH4(+) partitioning which affects the HNO3/NO3(-) partitioning.

aerosols↗
DOE OSTI2014

Exact results for models of multichannel quantum nonadiabatic transitions

We consider nonadiabatic transitions in explicitly time-dependent systems with Hamiltonians of the form Hˆ(t)=Aˆ+Bˆt+Cˆ/t, where t is time and Aˆ,Bˆ,Cˆ are Hermitian N × N matrices. We show that in any model of this type, scattering matrix elements satisfy nontrivial exact constraints that follow from the absence of the Stokes phenomenon for solutions with specific conditions at t→–∞. This allows one to continue such solutions analytically to t→+∞, and connect their asymptotic behavior at t→–∞ and t→+∞. This property becomes particularly useful when a model shows additional discrete symmetries. Specifically, we derive a number of simple exact constraints and explicit expressions for scattering probabilities in such systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗
NASA NTRS1980

Matrix management for aerospace 2000

The martix management approach to program management is an organized effort for attaining program objectives by defining and structuring all elements so as to form a single system whose parts are united by interaction. The objective of the systems approach is uncompromisingly complete coverage of the program management endeavor. Starting with an analysis of the functions necessary to carry out a given program, a model must be defined; a matrix of responsibility assignment must be prepared; and each operational process must be examined to establish how it is to be carried out and how it relates to all other processes.

Mccarthy, J. F., Jr.↗
DOE OSTI2013

On estimating cosmology-dependent covariance matrices

We describe a statistical model to estimate the covariance matrix of matter tracer two-point correlation functions with cosmological simulations. Assuming a fixed number of cosmological simulation runs, we describe how to build a 'statistical emulator' of the two-point function covariance over a specified range of input cosmological parameters. Because the simulation runs with different cosmological models help to constrain the form of the covariance, we predict that the cosmology-dependent covariance may be estimated with a comparable number of simulations as would be needed to estimate the covariance for fixed cosmology. Our framework is a necessary first step in planning a simulations campaign for analyzing the next generation of cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Robust root clustering for linear uncertain systems using generalized Lyapunov theory

Consideration is given to the problem of matrix root clustering in subregions of a complex plane for linear state space models with real parameter uncertainty. The nominal matrix root clustering theory of Gutman & Jury (1981) using the generalized Liapunov equation is extended to the perturbed matrix case, and bounds are derived on the perturbation to maintain root clustering inside a given region. The theory makes it possible to obtain an explicit relationship between the parameters of the root clustering region and the uncertainty range of the parameter space.

Yedavalli, R. K.↗
NASA NTRS1986

Acoustoelastic constants in dilute two-phase alloys

Acoustoelastic constants are calculated for two-phase alloys containing dilute concentrations of precipitates in a solid-solution matrix, on the basis of a model in which the precipitates are represented as a dilute elastic suspension of spherical particle inclusions in an infinite matrix. The longitudinal propagation velocity in the alloy is thereby obtained in terms of the precipitates' concentration and the elastic moduli of the two phases. Results are presented which indicate that the acoustoelastic constant of longitudinal waves in a dilute two-phase alloy varies linearly with the concentration of second-phase precipitates, in agreement with recent measurements in aluminum and steel alloys where the acoustoelastic constants changed linearly with the second phase's volume fraction.

Salama, K.↗
DOE OSTI2024

Significance of Low‐Velocity Zones on Solute Retention in Rough Fractures

Natural fractures are characterized by high internal heterogeneity. This internal variability is the cause of flow channeling, which in turn leads to contaminant transport taking place primarily along the high-velocity channels. Mass exchange between the high-velocity channels and the low-velocity zones has the potential to enhance contaminant retention, due to solute diffusion into the low-velocity zones and subsequent exposure to additional surface area for diffusion into the bordering rock matrix. Here, we derive a random walk particle tracking method for heterogeneous fractures, which includes an additional term to account for the aperture gradient. The method takes into account advection, diffusion in the fracture and matrix diffusion. The developed numerical framework is applied to assess the effect of low-velocity zones in rough self-affine fractures. The results show that diffusion into low-velocity zones has a visible but modest impact on contaminant retention. The magnitude of this impact does not change considerably, regardless of whether diffusion into the rock matrix is considered in the model, and increases for a decreasing average Péclet number of the fracture.

58 GEOSCIENCES↗
DOE OSTI2021

Phase diagram of the anisotropic triangular lattice Hubbard model

In a recent study [Phys. Rev. X 10, 021042 (2020)], we showed using large-scale density matrix renormalization group (DMRG) simulations on infinite cylinders that the triangular lattice Hubbard model has a chiral spin liquid phase. In this work, we introduce hopping anisotropy in the model, making one of the three distinct bonds on the lattice stronger or weaker compared with the other two. We implement the anisotropy in two inequivalent ways, one which respects the mirror symmetry of the cylinder and one which breaks this symmetry. In the full range of anisotropy, from the square lattice to weakly coupled one-dimensional chains, we find a variety of phases. Near the isotropic limit we find the three phases identified in our previous work: metal, chiral spin liquid, and 120$^\circ$ spiral order; we note that a recent paper suggests the apparently metallic phase may actually be a Luther-Emery liquid, which would also be in agreement with our results. Furthermore, when one bond is weakened by a relatively small amount, the ground state quickly becomes the square lattice N\'{e}el order. When one bond is strengthened, the story is much less clear, with the phases that we find depending on the orientation of the anisotropy and on the cylinder circumference. While our work is to our knowledge the first DMRG study of the anisotropic triangular lattice Hubbard model, the overall phase diagram we find is broadly consistent with that found previously using other methods, such as variational Monte Carlo and dynamical mean field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

GPU-optimized Code for Long-term Simulations of Beam-beam Effects in Colliders

We report on the development of the new code for long-term simulation of beam-beam effects in particle colliders. The underlying physical model relies on a matrix-based arbitrary-order symplectic particle tracking for beam transport and the Bassetti-Erskine approximation for beam-beam interaction. The computations are accelerated through a parallel implementation on a hybrid GPU/CPU platform. With the new code, a previously computationally prohibitive long-term simulations become tractable. We use the new code to model the proposed medium-energy electron-ion collider (MEIC) at Jefferson Lab.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Interfacial reaction in SiC{sub f}/Ti-6Al-4V composite by using transmission electron microscopy

The interfacial reactions of continuous SiC fiber reinforced Ti-6Al-4V matrix composite (SiC{sub f}/Ti-6Al-4V composite) and continuous SiC fiber coated by C reinforced Ti-6Al-4V matrix composite (SiC{sub f}/C/Ti-6Al-4V composite) were investigated by using micro-beam electron diffraction (MBED) and energy disperse spectroscopy (EDS) on transmission electron microscopy (TEM). The sequence of the interfacial reactions in the as-processed and exposed at 900°C for 50h SiC{sub f}/Ti-6Al-4V composites can be described as SiC||TiC||Ti{sub 5}Si{sub 3} + TiC||Ti-6Al-4V and SiC||TiC||Ti{sub 5}Si{sub 3}||TiC||Ti{sub 5}Si{sub 3}||TiC||Ti{sub 5}Si{sub 3}||Ti-6Al-4V, respectively. Additionally, both in as-processed and exposed composites, Ti{sub 3}SiC{sub 2} and Ti{sub 3}Si are absent at the interfaces. For the SiC{sub f}/C/Ti-6Al-4V composite exposed at 900 °C for 50 h, the sequence of the interfacial reaction can be described as SiC||C||TiC{sub F}||TiC{sub C}||Ti-6Al-4V before C coating is completely consumed by interfacial reaction. When interfacial reaction consumes C coating completely, the sequence of the interfacial reaction can be described as SiC||TiC||Ti{sub 5}Si{sub 3}||TiC||Ti-6Al-4V. Furthermore, in SiC{sub f}/C/Ti-6Al-4V composite, C coating can absolutely prevent Si diffusion from SiC fiber to matrix. Basing on these results, the model of formation process of the interfacial reaction products in the composites was proposed. - Highlights: • We obtained the sequence of the interfacial reactions in the as-processed and exposed at 900 °C for 50 h SiC{sub f}/Ti-6Al-4 V composites as well as in the SiC{sub f}/C/Ti-6Al-4 V composite exposed at 900 °C for 50 h. • We verified that both in as-processed and exposed SiC{sub f}/Ti-6Al-4 V composites, Ti{sub 3}SiC{sub 2} and Ti{sub 3}Si are absent at the interfaces. • Carbon coating can absolutely prevent silicon diffusion from SiC fiber to matrix. • Basing on these results, the model of formation process of the interfacial reaction products in the composites was proposed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bayesian Optimization Framework for Imperfect Data or Models

Conventional Bayesian optimization methods implicitly assume that the data and model being optimized are “perfect.” This assumption leads to inaccurate posterior probability distribution functions (PDFs) when applied to “imperfect” data or models. The new Bayesian optimization framework presented in this report provides a way to parameterize the effect of imperfections usually encountered in a prior PDF of generalized data or a model on the posterior PDF. The effects of imperfections are parameterized by a set of constraints imposed on the posterior expectation values of deviations between the data and the model and on their covariance matrix elements. A particular set of values for these constraints conveys an evaluator’s best estimate of the effect of imperfections on the corresponding posterior expectation values. When a prior PDF of generalized data is assumed to be normal, an expression for a posterior PDF satisfying an arbitrary set of constraints is derived analytically for linear models. An analogous iterative algorithm is given for nonlinear models. The corresponding posterior PDF should be used to estimate any posterior expectation values in the presence of imperfections parameterized by that set of constraints. A posterior PDF of a conventional Bayesian optimization method is recovered analytically when all evaluator-specified constraints are set to zero (i.e., in the absence of any imperfections). The analytical expressions derived in this report for normal PDFs and linear models were verified numerically by a Metropolis–Hastings Monte Carlo method. The methods presented herein could be applied to any kind of data or models, including differential cross-section data or integral benchmark experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

GPU-Accelerated Long-Term Simulations of Beam-Beam Effects in Colliders

We present an update on the development of the new code for long-term simulation of beam-beameffects in particle colliders. The underlying physical model relies on a matrix-based arbitrary-order particle tracking (including a symplectic option) for beam transport and the generalized Bassetti-Erskine approximation for beam-beam interaction. The computations are accelerated through a parallel implementation on a hybrid GPU/CPU platform. With the new code, previously computationally prohibitive longterm simulations become tractable. The new code will be used to model the proposed Medium-energy Electron-Ion Collider (MEIC) at Jefferson Lab.

Terzic, B.↗
NASA NTRS2011

Improvements to a Response Surface Thermal Model for Orion Mated to the International Space Station

This study is an extension of previous work to evaluate the applicability of Design of Experiments (DOE)/Response Surface Methodology to on-orbit thermal analysis. The goal was to determine if the methodology could produce a Response Surface Equation (RSE) that predicted the thermal model temperature results within +/-10 F. An RSE is a polynomial expression that can then be used to predict temperatures for a defined range of factor combinations. Based on suggestions received from the previous work, this study used a model with simpler geometry, considered polynomials up to fifth order, and evaluated orbital temperature variations to establish a minimum and maximum temperature for each component. A simplified Outer Mold Line (OML) thermal model of the Orion spacecraft was used in this study. The factors chosen were the vehicle's Yaw, Pitch, and Roll (defining the on-orbit attitude), the Beta angle (restricted to positive beta angles from 0 to 75), and the environmental constants (varying from cold to hot). All factors were normalized from their native ranges to a non-dimensional range from -1.0 to 1.0. Twenty-three components from the OML were chosen and the minimum and maximum orbital temperatures were calculated for each to produce forty-six responses for the DOE model. A customized DOE case matrix of 145 analysis cases was developed which used analysis points at the factor corners, mid-points, and center. From this data set, RSE s were developed which consisted of cubic, quartic, and fifth order polynomials. The results presented are for the fifth order RSE. The RSE results were then evaluated for agreement with the analytical model predictions to produce a +/-3(sigma) error band. Forty of the 46 responses had a +/-3(sigma) value of 10 F or less. Encouraged by this initial success, two additional sets of verification cases were selected. One contained 20 cases, the other 50 cases. These cases were evaluated both with the fifth order RSE and with the analytical model. For the maximum temperature predictions, 12 of the 23 components had all predictions within +/-10 F and 17 were within +/-20 F. For the minimum temperature predictions, only 4 of the 23 components (the four radiator temperatures), were within the 10 F goal. The maximum temperature RSEs were then run through 59,049 screening cases. The RSE predictions were then filtered to find 55 cases that produced the hottest temperatures. These 55 cases were then analyzed using the thermal model and the results compared against the RSE predictions. As noted earlier, 12 of the 23 responses were within +/-10 F at 17 within +/-20 F. These results demonstrate that if properly formulated, an RSE can provide a reliable, fast temperature prediction. Despite this progress, additional work is needed to determine why the minimum temperatures responses and 6 of the hot temperature responses did not produce reliable RSEs. Recommend focus areas are the model itself (arithmetic vs. diffusion nodes) and seeking consultations with statistical application experts.

Miller, StephenW.↗