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 163 records · Page 9

Noninteractive macroscopic reliability model for ceramic matrix composites with orthotropic material symmetry

A macroscopic noninteractive reliability model for ceramic matrix composites is presented. The model is multiaxial and applicable to composites that can be characterized as orthotropic. Tensorial invariant theory is used to create an integrity basis with invariants that correspond to physical mechanisms related to fracture. This integrity basis is then used to construct a failure function per unit volume (or area) of material. It is assumed that the overall strength of the composite is governed by weakest link theory. This leads to a Weibull type model similar in nature to the principle of independent action (PIA) model for isotropic monolithic ceramics. An experimental program to obtain model parameters is briefly discussed. In addition, qualitative features of the model are illustrated by presenting reliability surfaces for various model parameters.

Duffy, Stephen F.↗

Noninteractive macroscopic reliability model for ceramic matrix composites with orthotropic material symmetry

A macroscopic noninteractive reliability model for ceramic matrix composites is presented. The model is multiaxial and applicable to composites that can be characterized as orthotropic. Tensorial invariant theory is used to create an integrity basis with invariants that correspond to physical mechanisms related to fracture. This integrity basis is then used to construct a failure function per unit volume (or area) of material. It is assumed that the overall strength of the composite is governed by weakest link theory. This leads to a Weibull type model similar in nature to the principle of independent action (PIA) model for isotropic monolithic ceramics. An experimental program to obtain model parameters is briefly discussed. In addition, qualitative features of the model are illustrated by presenting reliability surfaces for various model parameters.

Duffy, Stephen F.↗

Bayesian event categorization matrix approach for explosion monitoring

Current efforts to correctly categorize natural events from suspected explosion sources with data that is collected by ground- or space-based sensors presents historical challenges that remain unaddressed by the Event Categorization Matrix (ECM) model. Smaller historical events (lower yield explosions) may have data available from fewer measurement techniques than are available today, and therefore, a historical event record can lack a complete set of discriminants. The covariance structures can also differ between such observations of event (source-type) categories. Both obstacles are problematic for the classic ECM model. Our work addresses this gap and presents a Bayesian update to the previous ECM model, termed the Bayesian Event Categorization Matrix model, which can be trained on partial observations and does not rely on a pooled covariance structure. We further augment the ECM model with Bayesian Decision Theory so that false negative or false positive rates of an event categorization can be reduced in an intuitive manner. To demonstrate improved categorization rates for the Bayesian Event Categorization Matrix model, we compare an array of Bayesian and classic models with multiple performance metrics using Monte Carlo experiments. We use both synthetic and real data. Our Bayesian models show consistent gains in overall accuracy and lower false negative rates relative to the classic ECM model. Here, we propose future avenues to improve Bayesian Event Categorization Matrix models’ decision making and predictive capability.

58 GEOSCIENCES↗

Matrix perturbation for analytical model improvement

A matrix perturbation method is proposed to calculate the Jacobian matrix and to compute the new eigendata for the parameter estimation procedure. The advantages of the method are the applicability to large complex structures without knowing the analytical expressions for the mass and stiffness matrices, and a cost effective approach for the re-computation of the eigendata. This method also allows the use of other measurements such as modal forces, kinetic energy distribution, and strain energy distributions in the estimation procedure. A realistic sample problem is presented to demonstrate the effectiveness of the proposed method.

Chen, J. C.↗

Plasmon coupling in topological insulator multilayers

Topological insulators (TIs) house two-dimensional Dirac plasmons on their surfaces. In TI thin films, plasmons on the top and bottom surfaces couple electrostatically, giving rise to an acoustic and an optical plasmon mode. By extension, a superlattice comprising TI layers and trivially-insulating layers could house multiple complex plasmon modes which could be used to create a new type of Dirac metamaterial. In this paper, we synthesize TI superlattices, fabricate them into stripe arrays, characterize their optical and plasmonic properties, and model the samples using transfer matrices. Furthermore, we excite plasmon modes that couple to the phonons in the superlattice, resulting in hybrid plasmon phonon polaritons. These modes are modeled using an analytical Fano resonance model and the extracted resonant positions are reproduced by transfer matrix modeling. We also excite an epsilon-near-zero mode in the top dielectric material (Bi 0.5 In 0.5 ) 2 Se 3 . Understanding the behavior of the polariton modes in this complex system will lend insight into the many-body interaction in low-dimensional systems.

36 MATERIALS SCIENCE↗

On the Relaxation Dynamics of Lohe Oscillators on Some Riemannian Manifolds

We study the collective relaxation dynamics appearing in weakly coupled Lohe oscillators in a large coupling regime. The Lohe models on the unit sphere and unitary group were proposed as a nonabelian generalization of the Kuramoto model on the unit circle and their emergent dynamics has been extensively studied in previous literature for some restricted class of initial data based on the Lyapunov functional approach and order parameter approach. In this paper, we extend the previous partial results to cover a generic initial configuration via the detailed analysis on the order parameter measuring the modulus of the centroid. In particular, we present a detailed relaxation dynamics and structure of the resulting asymptotic states for the Lohe sphere model. We also present new gradient flow formulations for the Lohe matrix models with the same one-body Hamiltonians on some group manifolds. As a direct application of this new formulation, we show that every bounded Lohe flow which originated from any initial configuration converges asymptotically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantitative comparison of opacities calculated using the R -matrix and distorted-wave methods: Fe xvii

ABSTRACT We present here a detailed calculation of opacities for Fe xvii at the physical conditions corresponding to the base of the Solar convection zone. Many ingredients are involved in the calculation of opacities. We review the impact of each ingredient on the final monochromatic and mean opacities (Rosseland and Planck). The necessary atomic data were calculated with the R-matrix and the distorted-wave (DW) methods. We study the effect of broadening, of resolution, of the extent of configuration sets and of configuration interaction to understand the differences between several theoretical predictions as well as the existing large disagreement with measurements. New Dirac R-matrix calculations including all configurations up to the n = 4, 5, and 6 complexes have been performed as well as corresponding Breit–Pauli DW calculations. The DW calculations have been extended to include autoionizing initial levels. A quantitative contrast is made between comparable DW and R-matrix models. We have reached self-convergence with n = 6 R-matrix and DW calculations. Populations in autoionizing initial levels contribute significantly to the opacities and should not be neglected. The R-matrix and DW results are consistent under the similar treatment of resonance broadening. The comparison with the experiment shows a persistent difference in the continuum while the filling of the windows shows some improvement. This study defines our path to the next generation of opacities and opacity tables for stellar modelling.

Delahaye, F.↗

Oxygen Diffusion and Reaction Kinetics in Continuous Fiber Ceramic Matrix Composites

Previous stressed oxidation tests of C/SiC composites at elevated temperatures (350 C to 1500 C) and sustained stresses (69 MPa and 172 MPa) have led to the development of a finite difference cracked matrix model. The times to failure in the samples suggest oxidation occurred in two kinetic regimes defined by the rate controlling mechanisms (i.e. diffusion controlled and reaction controlled kinetics). Microstructural analysis revealed preferential oxidation along as-fabricated, matrix microcracks and also suggested two regimes of oxidation kinetics dependent on the oxidation temperature. Based on experimental results, observation, and theory, a finite difference model was developed. The model simulates the diffusion of oxygen into a matrix crack bridged by carbon fibers. The model facilitates the study of the relative importance of temperature, the reaction rate constant, and the diffusion coefficient on the overall oxidation kinetics.

Halbig, Michael C.↗

Dark Energy Survey year 3 results: covariance modelling and its impact on parameter estimation and quality of fit

ABSTRACT We describe and test the fiducial covariance matrix model for the combined two-point function analysis of the Dark Energy Survey Year 3 (DES-Y3) data set. Using a variety of new ansatzes for covariance modelling and testing, we validate the assumptions and approximations of this model. These include the assumption of Gaussian likelihood, the trispectrum contribution to the covariance, the impact of evaluating the model at a wrong set of parameters, the impact of masking and survey geometry, deviations from Poissonian shot noise, galaxy weighting schemes, and other sub-dominant effects. We find that our covariance model is robust and that its approximations have little impact on goodness of fit and parameter estimation. The largest impact on best-fitting figure-of-merit arises from the so-called fsky approximation for dealing with finite survey area, which on average increases the χ2 between maximum posterior model and measurement by $3.7{{\ \rm per\ cent}}$ (Δχ2 ≈ 18.9). Standard methods to go beyond this approximation fail for DES-Y3, but we derive an approximate scheme to deal with these features. For parameter estimation, our ignorance of the exact parameters at which to evaluate our covariance model causes the dominant effect. We find that it increases the scatter of maximum posterior values for Ωm and σ8 by about $3{{\ \rm per\ cent}}$ and for the dark energy equation-of-state parameter by about $5{{\ \rm per\ cent}}$.

79 ASTRONOMY AND ASTROPHYSICS↗

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗

FEOTS v0.0.0: a new offline code for the fast equilibration of tracers in the ocean

Abstract. In this paper we introduce a new software framework for the offline calculation of tracer transport in the ocean. The Fast Equilibration of Ocean Tracers Software (FEOTS) is an end-to-end set of tools to efficiently calculate tracer distributions on a global or regional sub-domain using transport operators diagnosed from a comprehensive ocean model. To the best of our knowledge, this is the first application of a transport matrix model to an eddying ocean state. While a Newton–Krylov-based equilibration capability is still under development and not presented here, we demonstrate in this paper the transient modeling capabilities of FEOTS in an application focused on the Argentine Basin, where intense eddy activity and the Zapiola Anticyclone lead to strong mixing of water masses. The demonstration illustrates progress in developing offline passive tracer simulation capabilities, while highlighting the challenges of the impulse response functions approach in capturing tracer transports by a non-linear advection scheme. Our future work will focus on improving the computational efficiency of the code to reduce time-to-solution, using different basis functions to better represent non-linear advection operators, applying FEOTS to a parent model with unstructured grids (Ocean Model for Prediction Across Scales, MPAS-Ocean), and fully implementing a Newton–Krylov steady-state solver.

54 ENVIRONMENTAL SCIENCES↗

A Unified Level of Service Model for NASA Earth Science Data Stewardship

During the past year, the Interagency Implementation and Concepts Team (IMPACT) reviewed existing service models in use at various NASA data centers in an effort to produce a unified, cohesive, and comprehensive Level of Service Model for all of NASA Distributed Active Archive Centers (DAACs). NASA DAACs are responsible for ensuring NASA Earth Science data are accurately and securely ingested, distributed, supported, and preserved. The term “Service" as used here refers to the spectrum of data management activities and outputs provided by DAACs in support of the data cared fo by each data cente. The unified Level-of- Service (LoS) model described in this presentation utilizes both the NASA-defined data product category and the data processing level to easily identify an appropriate level-of-service to be applied to a data product throughout the full data life cycle. This LoS model is to be used by DAACs when appraising incoming data in order to determine the appropriate and required services to provide. The LoS model utilizes a 3-level system in which services build upon previous levels and thereby require greater commitment and effort both on the part of the DAAC personnel and the data producer at the highest level. The LoS model description also contains examples of ways to communicate with data producers and data users what services can be expected, thereby bringing more consistent user experiences across the enterprise. In this presentation, we will outline the features of the LoS model and describe how it relates to the FAIR data practices and the NOAA Maturity Matrix model.

Smith, Deborah↗

An optimum organizational structure for a large earth-orbiting multidisciplinary space base

An optimum hypothetical organizational structure was studied for a large earth-orbiting, multidisciplinary research and applications space base manned by a crew of technologists. Because such a facility does not presently exist, in situ empirical testing was not possible. Study activity was, therefore, concerned with the identification of a desired organizational structural model rather than with the empirical testing of the model. The essential finding of this research was that a four-level project type total matrix model will optimize the efficiency and effectiveness of space base technologists.

Ragusa, J. M.↗

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗

Dual Phase Synergy Enabled Large Elastic Strains of Nanoinclusions in a Dislocation Slip Matrix Composite

Freestanding nanomaterials (such as nanowires, nanoribbons and nanotubes) are known to exhibit ultra-large elastic strains and ultra-high strengths. However, harnessing their superior intrinsic mechanical properties in bulk composites has proven to be difficult. A recent breakthrough has overcome this difficulty by using a martensitic phase transforming matrix, in which ultra-large elastic strains approaching the theoretical limit is achieved in Nb nanowires embedded in the matrix. This discovery, whereas breaking a long-standing challenge, still limits our ability of harnessing the exceptional properties of nanomaterials and developing ultrahigh strength bulk materials to a narrow selection of phase transforming alloy matrices. In this study, we investigated the possibility to harness the intrinsic mechanical properties of nanoinclusions in conventional dislocation slip matrix based on a principle of synergy between the inclusion and the matrix. The small spacing between the densely populated hard and dislocation-impenetrable nanoinclusions departmentalize the plastic matrix into small domains to effectively impede dislocation motion within the matrix, inducing significant strengthening and large local elastic strains of the matrix, which in turn induced large elastic strains in the nanoinclusions. This dual phase synergy is verified in a Ti 3 Sn inclusions/B2-NiTi(Fe) plastic matrix model materials system. The maximum elastic strain of Ti 3 Sn inclusion obtained in the dislocation slip matrix is comparable to that achieved in a phase transforming matrix. As a result, this finding opens new opportunities for the development of high-strength nanocomposites.

42 ENGINEERING↗

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗

Improvement of a large dynamic analytical model using ground vibration test data

A method has been developed which uses measured normal modes and natural frequencies to improve an analytical mass and stiffness matrix model of a structure. The method directly identifies, without iteration, a set of minimum changes in the analytical matrices which force the eigensolutions to agree with the test measurements. An application is presented in which the analytical model had 508 degrees of freedom and 19 modes were measured at 101 locations on the structure. The resulting changes in the model are judged to be small compared to expectations of errors in the analysis. Thus the improved model is accepted as a reasonable model of the structure with improved dynamic response characteristics. In addition, it is shown that the procedure may be a useful tool in identifying apparent measured modes which are not true normal modes of the structures.

Berman, A.↗