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 91 records · Page 5

Improved semiclassical model for real-time evaporation of matrix black holes

In this paper, we study real-time classical matrix mechanics of a simplified 2 × 2 matrix model inspired by the black hole evaporation problem. This is a step towards making a quantitative model of real-time evaporation of a black hole, which is realized as a bound state of D0-branes in string theory. The model we study is the reduction of Yang–Mills in 2 + 1 dimension to 0 + 1 dimensions, which has been corrected with an additional potential that can be interpreted as a zero-point energy for fermions. Our goal is to understand the lifetime of such a classical bound state object in the classical regime. To do so, we pay particular attention to when D-particles separate to check that the “off-diagonal modes” of the matrices become adiabatic and use that information to improve on existing models of evaporation. It turns out that the naive expectation value of the lifetime with the fermionic correction is infinite. This is a logarithmic divergence that arises from very large excursions in the separation between the branes near the threshold for classical evaporation. The adiabatic behavior lets us get some analytic control of the dynamics in this regime to get this estimate. This divergence is cutoff in the quantum theory due to quantization of the adiabatic parameter, resulting in a long lifetime of the bound state, with a parametric dependence of order log(1/ℏ).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fuel Matrix Degradation Model Development Update: Alloy Corrosion Rates and Hydrogen Generation

This work is being performed as part of the DOE NE Spent Fuel and Waste Science and Technology Campaign, Argillite and Crystalline Rock R&D work packages: SF-19AN01030101 and SF-19AN01030201. This document meets the milestone M4SF-19AN010301013 for Argillite R&D and the milestone M4SF-19AN010302013 for Crystalline R&D. The fuel matrix degradation (FMD) model calculates the degradation rate of spent UO2 fuel based on fundamental electrochemical theory and thermodynamics. The model provides the quantitative basis for radionuclide source term predictions in the generic disposal system analysis (GDSA) performance assessment (PA) model. The FMD model has been implemented in a manner that facilitates its integration with the GDSA-PA code and a preliminary, integrated FMD-GDSA model has been successfully tested (Jerden et al., 2017). The specific focus of on-going work is to accurately and quantitatively represent the generation of H2 in a breached waste package and its effect on the spent fuel degradation rate so that this key process can be accurately represented in PA models. It has been shown experimentally that millimolar concentrations of dissolved H2 in contact with spent fuel will inhibit oxidative dissolution and decrease the fuel degradation rates by as much as three orders of magnitude (e.g., Röllin et al., 2001, Ollila, 2008). This will lead to a significant decrease in the radionuclide source term values used in repository PA models. Sensitivity studies using the FMD model have shown that the dissolved H2 concentration is the dominant environmental variable affecting the UO2 spent fuel dissolution rate as long as a source of H2 persists in the breached waste package (Jerden et al., 2015). Therefore, on-going experimental and modeling efforts focus on H2 production, accumulation and interfacial reactions within a breached waste package. The anoxic corrosion of metallic engineering materials within a breached waste package will be the main source of H2 in crystalline and argillite rock repository systems. The most important H2-producing alloys present within a waste package include stainless steel, carbon steel, and aluminum alloys and the Zircaloy fuel cladding. A steel corrosion module was added to the FMD model to account for these H2 sources and couple the H2 generation rate with fuel degradation processes (Jerden et al., 2017). As discussed below, there is a need for experimental data from electrochemical corrosion tests to quantify the H2 generation rates due to the corrosion of these waste package alloys as a function of key environmental variables (Eh, pH, groundwater chemistry). This report presents results from series of electrochemical corrosion tests with materials present in waste packages that provide information needed to parameterize and validate the FMD model. The alloys tested include 316L stainless steel (316SS), AISI 4320 carbon steel (C-steel) and Zircaloy-4 fuel cladding (unirradiated). Potentiodynamic polarization and potentiostatic polarization tests were performed in buffer solutions at pH 4, 7 and 10 with NaCl concentrations of 0, 4 and 100 millimolar (mM). The potentiostatic tests were performed at voltages varying from 0.5 volts down to -0.66 volts, vs the standard hydrogen electrode (SHE), and all tests were performed at laboratory ambient temperatures (~22 °C). The corrosion rates measured during the electrochemical tests were used as direct input to a new FMD/in-package chemistry model that is being developed to provide accurate long-term spent fuel degradation rates as chemical conditions (Eh, pH, speciation) evolve within a breached waste package. This new model was developed by coupling the FMD model with the reactive transport code X1t, which is a module within the Geochemist’s Workbench (GWB) reaction path modeling code (Jerden et al., 2019). The reactive transport model consists of a 1D domain discretized with 21 reaction diffusion cells and includes a single waste packag

Jerden, James↗

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗

Topological order in matrix Ising models

We study a family of models for an N_1 \times N_2 N 1 × N 2 matrix worth of Ising spins S_{aB} S a B . In the large N_i N i limit we show that the spins soften, so that the partition function is described by a bosonic matrix integral with a single ‘spherical’ constraint. In this way we generalize the results of to a wide class of Ising Hamiltonians with O(N_1,\mathbb{Z})\times O(N_2,\mathbb{Z}) O ( N 1 , ℤ ) × O ( N 2 , ℤ ) symmetry. The models can undergo topological large N N phase transitions in which the thermal expectation value of the distribution of singular values of the matrix S_{aB} S a B becomes disconnected. This topological transition competes with low temperature glassy and magnetically ordered phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Loop equation and exact soft anomalous dimension in $ \mathcal{N} $ = 4 super Yang-Mills

BPS Wilson loops in supersymmetric gauge theories have been the subjects of active research since they are often amenable to exact computation. So far most of the studies have focused on loops that do not intersect. In this paper, we derive exact results for intersecting 1/8 BPS Wilson loops in $ \mathcal{N} $ = 4 supersymmetric Yang-Mills theory, using a combination of supersymmetric localization and the loop equation in 2d gauge theory. The result is given by a novel matrix-model-like representation which couples multiple contour integrals and a Gaussian matrix model. We evaluate the integral at large N, and make contact with the string worldsheet description at strong coupling. As an application of our results, we compute exactly a small-angle limit (and more generally near-BPS limits) of the cross anomalous dimension which governs the UV divergence of intersecting Wilson lines. The same quantity describes the soft anomalous dimension of scattering amplitudes of W-bosons in the Coulomb branch.

't Hooftand Polyakov loops↗

Real-time monitoring of trace noble gases using laser-induced breakdown spectroscopy—An investigation of the impact of bulk gas on plasma properties and sensitivity

The impact of Ar and He bulk gases on laser-induced breakdown spectroscopy (LIBS) real-time monitoring of trace Xe and Kr was assessed. LIBS is being developed as a monitoring tool for measuring noble gas transport in molten salt systems, in which traditional sensors may face challenges associated with radiation, corrosive materials, and/or mixed phases. The plasma temperature and electron densities of LIBS plasmas were measured in both static and various flowing Ar and He streams (0–5 L min −1 ). The use of an Ar bulk gas resulted in higher plasma temperature, greater electron densities by an order of magnitude, and extended plasma lifetime compared with when He bulk gas was used. Gas flow rate was found to have little impact on plasma temperature; however, its effect on electron density was significant, indicating the need to consider flow rate–specific models. Matrix effects on emission peaks were reported for both bulk gases. Due to these matrix effects, multivariate models were developed for Xe and Kr ranging from 0 to 700 ppm in both bulk gases. Although the predictive behavior was similar (root mean square error of prediction ranging from 11.1 to 20.6 ppm), the limits of detection were superior in He (Xe: 22.9 ppm, Kr: 30.4 ppm). Furthermore, these models were employed in demonstrative real-time tests (>1 h), which showed strong predictive precision (relative standard deviation <5 %) regardless of the bulk gas. Ultimately, this study provides a guide for the considerations required when developing gaseous LIBS models for real-time monitoring.

Gas flow effects↗

Recent progress in molecular engineering to tailor organic–inorganic interfaces in composite membranes

Organic-inorganic composite membranes are of great interest in modern water treatment processes because they offer potentially superior separation efficiency and advanced funtionality by integrating the properties of polymers and inorganics. The biggest challenge in the fabrication and applications of organic-inorganic composite membranes is the incompatibility of organic-inorganic interfaces. In this minireview, we summarize the most recent advances in molecular engineering to tailor the properties of interfaces in composite membranes. Three typical models (i.e. mixed matrix model, interface composite model, and dual-layer composite model) are presented to demonstrate how to regulate these interfaces via molecular engineering and how the interfacial properties ultimately affect the membrane performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LS-DYNA Implementation of Polymer Matrix Composite Model Under High Strain Rate Impact

A recently developed constitutive model is implemented into LS-DYNA as a user defined material model (UMAT) to characterize the nonlinear strain rate dependent behavior of polymers. By utilizing this model within a micromechanics technique based on a laminate analogy, an algorithm to analyze the strain rate dependent, nonlinear deformation of a fiber reinforced polymer matrix composite is then developed as a UMAT to simulate the response of these composites under high strain rate impact. The models are designed for shell elements in order to ensure computational efficiency. Experimental and numerical stress-strain curves are compared for two representative polymers and a representative polymer matrix composite, with the analytical model predicting the experimental response reasonably well.

Zheng, Xia-Hua↗

Solving puzzles in deformed JT gravity: phase transitions and non-perturbative effects

Recent work has shown that certain deformations of the scalar potential in Jackiw-Teitelboim gravity can be written as double-scaled matrix models. However, some of the deformations exhibit an apparent breakdown of unitarity in the form of a negative spectral density at disc order. We show here that the source of the problem is the presence of a multi-valued solution of the leading order matrix model string equation. While for a class of deformations we fix the problem by identifying a first order phase transition, for others we show that the theory is both perturbatively and non-perturbatively inconsistent. Aspects of the phase structure of the deformations are mapped out, using methods known to supply a non-perturbative definition of undeformed JT gravity. Some features are in qualitative agreement with a semi-classical analysis of the phase structure of two-dimensional black holes in these deformed theories.

2D gravity↗

Random matrix approach to three-dimensional QCD with a Chern-Simons term

We propose a random matrix theory for QCD in three dimensions with a Chern-Simons term at level k which spontaneously breaks the flavor symmetry according to U(2Nf) → U(Nf+k)×U(Nf-k). This random matrix model is obtained by adding a complex part to the action for the k = 0 random matrix model. We derive the pattern of spontaneous symmetry breaking from the analytical solution of the model. Additionally, we obtain explicit analytical results for the spectral density and the spectral correlation functions for the Dirac operator at finite matrix dimension, that become complex. In the microscopic domain where the matrix size tends to infinity, they are expected to be universal, and give an exact analytical prediction to the spectral properties of the Dirac operator in the presence of a Chern-Simons term. Here, we calculate the microscopic spectral density. It shows exponentially large (complex) oscillations which cancel the phase of the k = 0 theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of the accuracy of an offline seasonally-varying matrix transport model for simulating ideal age

Newton-Krylov solvers for ocean tracers have the potential to greatly decrease the computational costs of spinning up deep-ocean tracers, which can take several thousand model years to reach equilibrium with surface processes. One version of the algorithm uses offline tracer transport matrices to simulate an annual cycle of tracer concentrations and applies Newton’s method to find concentrations that are periodic in time. Here we present the impact of time-averaging the transport matrices on the equilibrium values of an ideal-age tracer. We compared annually-averaged, monthly-averaged, and 5-day-averaged transport matrices to an online simulation using the ocean component of the Community Earth System Model (CESM) with a nominal horizontal resolution of 1° × 1° and 60 vertical levels. We found that increasing the time resolution of the offline transport model reduced a low age bias from 12% for the annually-averaged transport matrices, to 4% for the monthly-averaged transport matrices, and to less than 2% for the transport matrices constructed from 5-day averages. The largest differences were in areas with strong seasonal changes in the circulation, such as the Northern Indian Ocean. As a result, for many applications the relatively small bias obtained using the offline model makes the offline approach attractive because it uses significantly less computer resources and is simpler to set up and run.

54 ENVIRONMENTAL SCIENCES↗

High Strain Rate Deformation Modeling of a Polymer Matrix Composite: Composite Micromechanical Model - Part 2

Recently applications have exposed polymer matrix composite materials to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under these extreme conditions. In this second paper of a two part report, a three-dimensional composite micromechanical model is described which allows for the analysis of the rate dependent, nonlinear deformation response of a polymer matrix composite. Strain rate dependent inelastic constitutive equations utilized to model the deformation response of a polymer are implemented within the micromechanics method. The deformation response of two representative laminated carbon fiber reinforced composite materials with varying fiber orientation has been predicted using the described technique. The predicted results compare favorably to both experimental values and the response predicted by the Generalized Method of Cells, a well-established micromechanics analysis method.

Goldberg, Robert K.↗

Atomic Layer Deposition of Vanadium Oxide to Reduce Parasitic Absorption and Improve Stability in n-i-p Perovskite Solar Cells for Tandems

Two critical issues associated with semi-transparent, n-i-p perovskite solar cells for 2-terminal tandem devices are parasitic absorption and long-term instability associated with the widely used spiro-OMeTAD and MoOx hole transport and buffer layers, respectively. In this work, we present an alternative hole contact bilayer that consists of a 30 nm undoped layer of spiro-TTB in conjunction with 9 nm of air-stable vanadium oxide (VOx) deposited via atomic layer deposition. The low absorption of UV and visible light in this bilayer results in the fabrication of a semi-transparent perovskite cell with 18.9 mA cm-2 of photocurrent, a 14% increase compared to the 16.6 mA cm-2 generated in a control device with 150 nm of doped spiro-OMeTAD. The ALD VOx buffer layer shows promise as a stable alternative to MoOx; an unencapsulated Cs0.17FA0.83Pb(Br0.17I0.83)3 device with ALD VOx and ITO as the top contact maintains its efficiency following 1000 hours at 85 degrees C in a N2 environment. Lastly, we use transfer matrix modeling of the optimized perovskite stack to predict its optical performance in a monolithic tandem cell with heterojunction silicon.

14 SOLAR ENERGY↗

Nonlinearity and Strain-Rate Dependence in the Deformation Response of Polymer Matrix Composites Modeled

There has been no accurate procedure for modeling the high-speed impact of composite materials, but such an analytical capability will be required in designing reliable lightweight engine-containment systems. The majority of the models in use assume a linear elastic material response that does not vary with strain rate. However, for containment systems, polymer matrix composites incorporating ductile polymers are likely to be used. For such a material, the deformation response is likely to be nonlinear and to vary with strain rate. An analytical model has been developed at the NASA Glenn Research Center at Lewis Field that incorporates both of these features. A set of constitutive equations that was originally developed to analyze the viscoplastic deformation of metals (Ramaswamy-Stouffer equations) was modified to simulate the nonlinear, rate-dependent deformation of polymers. Specifically, the effects of hydrostatic stresses on the inelastic response, which can be significant in polymers, were accounted for by a modification of the definition of the effective stress. The constitutive equations were then incorporated into a composite micromechanics model based on the mechanics of materials theory. This theory predicts the deformation response of a composite material from the properties and behavior of the individual constituents. In this manner, the nonlinear, rate-dependent deformation response of a polymer matrix composite can be predicted.

Goldberg, Robert K.↗

A note on the Bayesian analysis of experiments with correlated multiple responses using a matrix-logarithmic covariance model

In the literature, analysis of multiple responses from experiments with replicates has modeled the covariance matrix directly as linear models of the transformed variances and correlations, i.e., covariance modeling. This article considers models based on the matrix- logarithm of the covariance matrix. This so-called log-covariance modeling is illustrated with data from actual experiments and compared with the traditional covariance modeling.

97 MATHEMATICS AND COMPUTING↗

Accurate prediction of diradical chemistry from a single-reference density-matrix method: Model application to the bicyclobutane to gauche-1,3-butadiene isomerization

Multireference correlation in diradical molecules can be captured by a single-reference 2-electron reduced-density-matrix (2-RDM) calculation with only single and double excitations in the 2-RDM parametrization. The 2-RDM parametrization is determined by N-representability conditions that are non-perturbative in their treatment of the electron correlation. Conventional single-reference wave function methods cannot describe the entanglement within diradical molecules without employing triple- and potentially even higher-order excitations of the mean-field determinant. In the isomerization of bicyclobutane to gauche-1,3-butadiene the parametric 2-RDM (p2-RDM) method predicts that the diradical disrotatory transition state is 58.9 kcal/mol above bicyclobutane. This barrier is in agreement with previous multireference calculations as well as recent Monte Carlo and higher-order coupled cluster calculations. The p2-RDM method predicts the Nth natural-orbital occupation number of the transition state to be 0.635, revealing its diradical character. The optimized geometry from the p2-RDM method differs in important details from the complete-active-space self-consistent-field geometry used in many previous studies including the Monte Carlo calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Krylov complexity in mixed phase space

We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.

chaos & nonlinear dynamics↗

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗