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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Parametric effects in nanobeams and AFM

Vibration dynamics of forced cantilever beams that are used in nanotechnology such as atomic force microscope modeling and carbon nanotubes is considered in terms of a fundamental response within a matrix framework. The modeling equations are written as a matrix differential equation subject to tip-sample general boundary conditions. At the junctions, where there are discontinuities due to different material or beam thickness, compatibility conditions are prescribed. Forced responses are given by convolution of the input load with the time domain Green matrix function. The corresponding matrix transfer function and modes of a multispan cantilever beam are determined in terms of solution basis of the same shape generated by a fundamental solution. Simulations were performed for a three stepped beam with a piezoelectric patch subject to pulse forcing terms and with surface effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Surrogate Model Development of Spent Fuel Degradation for Repository Performance Assessment

In model simulations of deep geologic repositories, UO 2 fuel matrix degradation typically begins as soon as the waste package breaches and groundwater contacts the fuel surface. The initial degradation rate depends on the timing of these events, burnup of the fuel, temperature, and concentrations of dissolved reactants. Estimating the initial rate of degradation is fairly straightforward, but as UO 2 corrosion products precipitate on the fuel surface and the movement of dissolved species between the fuel surface and environment is impeded by the precipitated solids, the rate is more difficult to quantify. At that point, calculating the degradation rate becomes a reactive-transport problem in which a large number of equations must be solved by iteration for a large number of grid cells at each time step. The consequence is that repository simulations, which are already expensive, become much more expensive, especially when hundreds or thousands of waste packages breach. The Fuel Matrix Degradation (FMD) model is the process model of the Spent Fuel and Waste Science and Technology (SFWST) campaign of the US Department of Energy (DOE). It calculates spent fuel degradation rates as a function of radiolysis, redox reactions, electrochemical reactions, alteration layer growth, and diffusion of reactants through the alteration layer. Like other similar fuel degradation process models, it is a complicated model requiring a large number of calculations and iterations at each time step.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Modeling Workloads of a Linear Electromagnetic Code for Load Balancing Matrix Assembly

This report presents our work to model the workloads of a linear electromagnetic application based on the method of moments in the frequency domain to effectively load balance the matrix assembly. This application is particularly challenging to load balance due to its lack of persistent iterative behavior, its operation under tight memory constraint (where the matrix may fill 80% of memory on each node), and the algorithmic complexity of the computational method. This report describes the first step in our work to apply an inspector-executor approach for load balancing workloads where key parameters are exposed during the inspector phase and a pre-trained model is applied to predict relative task weights for the load balancer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
DOE OSTI2016

New type of hill-top inflation

We suggest a new type of hill-top inflation originating from the initial conditions in the form of the microcanonical density matrix for the cosmological model with a large number of quantum fields conformally coupled to gravity. Initial conditions for inflation are set up by cosmological instantons describing underbarrier oscillations in the vicinity of the inflaton potential maximum. These periodic oscillations of the inflaton field and cosmological scale factor are obtained within the approximation of two coupled oscillators subject to the slow roll regime in the Euclidean time. This regime is characterized by rapid oscillations of the scale factor on the background of a slowly varying inflaton, which guarantees smallness of slow roll parameters ϵ and η of the following inflation stage. A hill-like shape of the inflaton potential is shown to be generated by logarithmic loop corrections to the tree-level asymptotically shift-invariant potential in the non-minimal Higgs inflation model and R{sup 2}-gravity. The solution to the problem of hierarchy between the Planckian scale and the inflation scale is discussed within the concept of conformal higher spin fields, which also suggests the mechanism bringing the model below the gravitational cutoff and, thus, protecting it from large graviton loop corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

$B^0_{(s)}$-mixing matrix elements from lattice QCD for the Standard Model and beyond

We calculate—for the first time in three-flavor lattice QCD—the hadronic matrix elements of all five local operators that contribute to neutral B 0 - and B s -meson mixing in and beyond the Standard Model. We present a complete error budget for each matrix element and also provide the full set of correlations among the matrix elements. We also present the corresponding bag parameters and their correlations, as well as specific combinations of the mixing matrix elements that enter the expression for the neutral B-meson width difference. We obtain the most precise determination to date of the SU(3)-breaking ratio ξ=1.206(18)(6), where the second error stems from the omission of charm-sea quarks, while the first encompasses all other uncertainties. The threefold reduction in total uncertainty, relative to the 2013 Flavor Lattice Averaging Group results, tightens the constraint from B mixing on the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light-valence quarks and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing effects to extrapolate the five matrix elements to the physical point. We combine our results with experimental measurements of the neutral B-meson oscillation frequencies to determine the CKM matrix elements |V td | = 8.00(34)(8)×10 -3 , |V ts | = 39.0(1.2)(0.4)×10 -3 , and |V td /V ts | = 0.2052(31)(10), which differ from CKM-unitarity expectations by about 2σ. In addition, these results and others from flavor-changing-neutral currents point towards an emerging tension between weak processes that are mediated at the loop and tree levels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Shared Dosimetry Error in Epidemiological Dose-Response Analyses

Radiation dose reconstruction systems for large-scale epidemiological studies are sophisticated both in providing estimates of dose and in representing dosimetry uncertainty. For example, a computer program was used by the Hanford Thyroid Disease Study to provide 100 realizations of possible dose to study participants. The variation in realizations reflected the range of possible dose for each cohort member consistent with the data on dose determinates in the cohort. Another example is the Mayak Worker Dosimetry System 2013 which estimates both external and internal exposures and provides multiple realizations of "possible" dose history to workers given dose determinants. This paper takes up the problem of dealing with complex dosimetry systems that provide multiple realizations of dose in an epidemiologic analysis. In this paper we derive expected scores and the information matrix for a model used widely in radiation epidemiology, namely the linear excess relative risk (ERR) model that allows for a linear dose response (risk in relation to radiation) and distinguishes between modifiers of background rates and of the excess risk due to exposure. We show that treating the mean dose for each individual (calculated by averaging over the realizations) as if it was true dose (ignoring both shared and unshared dosimetry errors) gives asymptotically unbiased estimates (i.e. the score has expectation zero) and valid tests of the null hypothesis that the ERR slope β is zero. Although the score is unbiased the information matrix (and hence the standard errors of the estimate of β) is biased for β≠0 when ignoring errors in dose estimates, and we show how to adjust the information matrix to remove this bias, using the multiple realizations of dose. Use of these methods for several studies, including the Mayak Worker Cohort and the U.S. Atomic Veterans Study, is discussed.

epidemiological analysis↗
DOE OSTI2014

TCDD induces dermal accumulation of keratinocyte-derived matrix metalloproteinase-10 in an organotypic model of human skin

The epidermis of skin is the first line of defense against the environment. A three dimensional model of human skin was used to investigate tissue-specific phenotypes induced by the environmental contaminant, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Continuous treatment of organotypic cultures of human keratinocytes with TCDD resulted in intracellular spaces between keratinocytes of the basal and immediately suprabasal layers as well as thinning of the basement membrane, in addition to the previously reported hyperkeratinization. These tissue remodeling events were preceded temporally by changes in expression of the extracellular matrix degrading enzyme, matrix metalloproteinase-10 (MMP-10). In organotypic cultures MMP-10 mRNA and protein were highly induced following TCDD treatment. Q-PCR and immunoblot results from TCDD-treated monolayer cultures, as well as indirect immunofluorescence and immunoblot analysis of TCDD-treated organotypic cultures, showed that MMP-10 was specifically contributed by the epidermal keratinocytes but not the dermal fibroblasts. Keratinocyte-derived MMP-10 protein accumulated over time in the dermal compartment of organotypic cultures. TCDD-induced epidermal phenotypes in organotypic cultures were attenuated by the keratinocyte-specific expression of tissue inhibitor of metalloproteinase-1, a known inhibitor of MMP-10. These studies suggest that MMP-10 and possibly other MMP-10-activated MMPs are responsible for the phenotypes exhibited in the basement membrane, the basal keratinocyte layer, and the cornified layer of TCDD-treated organotypic cultures. Our studies reveal a novel mechanism by which the epithelial–stromal microenvironment is altered in a tissue-specific manner thereby inducing structural and functional pathology in the interfollicular epidermis of human skin. - Highlights: • TCDD causes hyperkeratosis and basement membrane changes in a model of human skin. • TCDD induces MMP-10 expression in organotypic cultures of human keratinocytes. • Keratinocyte-expressed MMP-10 accumulates in the dermal compartment. • Keratinocyte K14 promoter-driven TIMP-1 expression ablates TCDD-induced phenotypes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Predicting ductility of Mg/SiCp nanocomposite under multiaxial loading conditions based on unit cell modeling

This article presents an investigation on the ductility of Mg/SiCp under nine different multiaxial loading conditions based on a unit cell (or called representative volumetric element, RVE) finite element model. The constitutive relationship of magnesium matrix is characterized by a model combined Lode angle dependent plasticity model with the CPB06 plasticity model, which is able to describe the tension-compression asymmetry in magnesium matrix. The debonding of the matrix-particle interface and matrix fracture are considered by the modified Mohr-Coulomb (MMC) fracture model. Simulation results show that the interface debonding and subsequent matrix fracture are responsible for the ductility and fracture of nanocomposite. The computational results with parametric studies are compared and well correlated with experimental data, and it finally gives a stress-based full fracture locus of Mg/SiCp nanocomposite.

42 ENGINEERING↗
DOE OSTI2016

Model Predictive Control of A Matrix-Converter Based Solid State Transformer for Utility Grid Interaction

The matrix converter solid state transformer (MC-SST), formed from the back-to-back connection of two three-to-single-phase matrix converters, is studied for use in the interconnection of two ac grids. The matrix converter topology provides a light weight and low volume single-stage bidirectional ac-ac power conversion without the need for a dc link. Thus, the lifetime limitations of dc-bus storage capacitors are avoided. However, space vector modulation of this type of MC-SST requires to compute vectors for each of the two MCs, which must be carefully coordinated to avoid commutation failure. An additional controller is also required to control power exchange between the two ac grids. In this paper, model predictive control (MPC) is proposed for an MC-SST connecting two different ac power grids. The proposed MPC predicts the circuit variables based on the discrete model of MC-SST system and the cost function is formulated so that the optimal switch vector for the next sample period is selected, thereby generating the required grid currents for the SST. Simulation and experimental studies are carried out to demonstrate the effectiveness and simplicity of the proposed MPC for such MC-SST-based grid interfacing systems.

Xue, Yaosuo↗
DOE OSTI2016

Role of vertex corrections in the matrix formulation of the random phase approximation for the multiorbital Hubbard model

In the framework of a multiorbital Hubbard model description of superconductivity, a matrix formulation of the superconducting pairing interaction that has been widely used is designed to treat spin, charge, and orbital fluctuations within a random phase approximation (RPA). In terms of Feynman diagrams, this takes into account particle-hole ladder and bubble contributions as expected. It turns out, however, that this matrix formulation also generates additional terms which have the diagrammatic structure of vertex corrections. Furthermore we examine these terms and discuss the relationship between the matrix-RPA superconducting pairing interaction and the Feynman diagrams that it sums.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Numerical Modeling of Fiber-Reinforced Metal Matrix Composite Processing by the Liquid Route: Literature Contribution

One of the technologies used to produce metal matrix composites (MMCs) is liquid route processing. One solution is to inject a liquid metal under pressure or at constant rate through a fibrous preform. This foundry technique overcomes the problem of the wettability of ceramic fibers by liquid metal. The liquid route can also be used to produce semiproducts by coating a filament with a molten metal. These processes involve physical phenomena combined with mass and heat transfer and phase change. The phase change phenomena related to solidification and also to the melting of the metal during the process notably result in modifications to the permeability of porous media, in gaps in impregnation, in the appearance of defects (porosities), and in segregation in the final product. In this article, we provide a state-of-the-art review of numerical models and simulation developed to study these physical phenomena involved in MMC processing by the liquid route.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Nonequilibrium Steady-State Transport in Quantum Impurity Models: A Thermofield and Quantum Quench Approach Using Matrix Product States

Here, the numerical renormalization group (NRG) is tailored to describe interacting impurity models in equilibrium, but it faces limitations for steady-state nonequilibrium, arising, e.g., due to an applied bias voltage. We show that these limitations can be overcome by describing the thermal leads using a thermofield approach, integrating out high energy modes using NRG, and then treating the nonequilibrium dynamics at low energies using a quench protocol, implemented using the time-dependent density matrix renormalization group. This yields quantitatively reliable results for the current (with errors ≲3%) down to the exponentially small energy scales characteristic of impurity models. We present results of benchmark quality for the temperature and magnetic field dependence of the zero-bias conductance peak for the single-impurity Anderson model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Modular Growth and Nucleation Formulations in Modular Plasticity Models

Recent investigations like the second and third Sandia Fracture Challenges have characterized and demonstrated the performance of a variety of failure techniques and models. These surveys have considered a wide breadth of models encapsulating both general failure criteria as well as those focusing on pore nucleation and growth. Extensive reviews exist on both topics. The former category generally consists of classic models like the Johnson-Cook or Wilkins criteria. These models were recently added to modular plasticity models in the Library of Advanced Materials for Engineering (LAME) as criteria for use with element death capabilities. The latter category was not treated in that effort. There exists a large class of failure models based on predicting the evolution of pores and failure associated with such microstructures. While the exact mechanisms and corresponding impact on the macroscale behavior remain an active area of research, a large suite of formulations have been proposed combining different features of both pore nucleation and subsequent growth. The most famous of these are based on the popular Gurson model of pore growth derived via micromechanical analysis assuming a plastically incompressible matrix. Numerous other models exist for both growth and nucleation and the Cocks-Ashby growth and Horstemeyer-Gokhale nucleation models have been used successfully in recent Sandia Fracture Challenges. This specific combination is colloquially referred to as the "BCJ-failure model as it has been frequently used with the Bammann-Chisea-Johnson plasticity model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bayesian Inference of State-Level COVID-19 Basic Reproduction Numbers across the United States

Although many persons in the United States have acquired immunity to COVID-19, either through vaccination or infection with SARS-CoV-2, COVID-19 will pose an ongoing threat to non-immune persons so long as disease transmission continues. We can estimate when sustained disease transmission will end in a population by calculating the population-specific basic reproduction number ℛ0, the expected number of secondary cases generated by an infected person in the absence of any interventions. The value of ℛ0 relates to a herd immunity threshold (HIT), which is given by 1−1/ℛ0. When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely (barring mutations allowing SARS-CoV-2 to escape immunity). Here, we report state-level ℛ0 estimates obtained using Bayesian inference. Maximum a posteriori estimates range from 7.1 for New Jersey to 2.3 for Wyoming, indicating that disease transmission varies considerably across states and that reaching herd immunity will be more difficult in some states than others. ℛ0 estimates were obtained from compartmental models via the next-generation matrix approach after each model was parameterized using regional daily confirmed case reports of COVID-19 from 21 January 2020 to 21 June 2020. Our ℛ0 estimates characterize the infectiousness of ancestral strains, but they can be used to determine HITs for a distinct, currently dominant circulating strain, such as SARS-CoV-2 variant Delta (lineage B.1.617.2), if the relative infectiousness of the strain can be ascertained. On the basis of Delta-adjusted HITs, vaccination data, and seroprevalence survey data, we found that no state had achieved herd immunity as of 20 September 2021.

60 APPLIED LIFE SCIENCES↗