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.
Results of the 2013 UT modeling benchmark obtained with models implemented in CIVA
The 2013 Ultrasonic Testing (UT) modeling benchmark concerns direct echoes from side drilled holes (SDH), flat bottom holes (FBH) and corner echoes from backwall breaking artificial notches inspected with a matrix phased array probe. This communication presents the results obtained with the models implemented in the CIVA software: the pencilmodel is used to compute the field radiated by the probe, the Kirchhoff approximation is applied to predict the response of FBH and notches and the SOV (Separation Of Variables) model is used for the SDH responses. The comparison between simulated and experimental results are presented and discussed.
Feeder Power Disaggregation: A Data-Efficient Matrix Completion Approach
This paper presents a data-driven algorithm for the feeder power disaggregation problem in distribution systems. Leveraging spatio-temporal power patterns in residential homes, residential power is discomposed into three components: sparse-switching loads, periodic loads, and photovoltaic (PV) generation, which are characterized through the design of two sparse matrices and a low-rank matrix. The matrix completion process is data-efficient because of the matrix sparsity and low rankness, along with the use of power system models. The proposed approach is tested using real-world residential data set on a 33-bus distribution system, demonstrating accurate power disaggregation with efficient matrix completion.
Multi-scale modeling of wastage layer formation in metallic fuel cladding
Fuel-cladding chemical interaction (FCCI) is a major concern for U-Zr metallic fuels' performance, primarily due to the formation of a brittle layer (wastage) in the cladding. This brittle layer, resulting from intermetallic compounds between cladding constituents Fe, Cr, and lanthanide fission products, significantly impacts the cladding's mechanical integrity. Recent efforts focus on developing a mechanistic modeling framework to understand lanthanide production, transport to the fuel-cladding interface, and phase transformation to intermetallic phases. A multi-scale computational approach has been used to calculate lanthanide transport rates, with atomistic calculations determining Nd diffusivities through the solid fuel matrix and along pore surfaces. These diffusivities inform a mesoscale model to determine an effective diffusion coefficient, accounting for porosity and infiltration with bond sodium. This effective diffusivity is used in engineering-scale simulations via the BISON fuel performance code, which has been validated against EBR-II and FFTF reactor experiments.
Interlayer fusion bonding of semi-crystalline polymer composites in extrusion deposition additive manufacturing
This work focuses on the evolution of interlayer fracture toughness properties of fiber-reinforced, semi-crystalline polymers in the extrusion deposition additive manufacturing (EDAM) process. Further, this work bridges the gap between the additive process conditions (time-temperature history) and the effective layer-to-layer fracture properties developed within a printed component. This is the first step to predict delamination that can occur during printing, during cooling to room temperature after printing, and during service performance of an additively manufactured geometry. A phenomenological model is developed for fusion bonding of semi-crystalline polymer matrix composites by coupling the interdiffusion of polymer chains with the evolution of polymer crystallinity. While the interdiffusion is captured by reptation theory of polymer dynamics, the evolution of crystallinity is modeled by phenomenological crystallization kinetics and crystal melting dynamics. Further, a methodology is developed to determine the critical strain energy release rate, G IC of the interlayer interface and experiments are conducted utilizing the double cantilever beam fracture test geometry. In conclusion, predictions of GIC as a function of thermal history are compared with experiments.
Computer code for double beta decay QRPA based calculations
The computer code developed by our group some years ago for the evaluation of nuclear matrix elements, within the QRPA and PQRPA nuclear structure models, involved in neutrino-nucleus reactions, muon capture and β{sup ±} processes, is extended to include also the nuclear double beta decay.
Energy levels and zero field splitting parameter for Fe{sup 2+} doped in ZnS
The aim of present paper is to report the results on the modeling of the crystal field parameters of Fe{sup 2+} doped in host matrix ZnS, simulate the energy levels scheme and calculate the zero field splitting parameter D of such system. The crystal field parameters were modeled in the frame of the superposition model of crystal field and the simulation of the energy levels scheme and calculation of the zero field splitting parameters done by diagonalization the Hamiltonian of Fe{sup 2+}:ZnS system. The obtained results were disscused and compared with experimental data. Satisfactory agreement have been obtained.
Improved Statistical Analysis for the Neutrinoless Double-Beta Decay Matrix Element of 136Xe
Neutrinoless double beta decay nuclear matrix element (M0ν) for 136Xe was recently analyzed using a statistical approach (Phys. Rev. C 107, 045501 (2023)). In the analysis, three initial shell model effective Hamiltonians were randomly altered, and their results for 23 measured observables were used to infer credibility for the M0ν nuclear matrix element (NME) based on a Bayesian Model Averaging approach. In that analysis, a reasonable Gamow-Teller quenching factor of 0.7 was assumed for each starting effective Hamiltonian. Given that the result of the statistical analysis was sensible to this choice, we are here improving that analysis by assuming that the Gamow-Teller quenching factor is also randomly chosen within reasonabe limits for all three starting Hamiltonians. The outcomes are slightly higher expectation values and uncertainties for the M0ν NME.
Hybrid programming-model strategies for GPU offloading of electronic structure calculation kernels
To address the challenge of performance portability and facilitate the implementation of electronic structure solvers, we developed the basic matrix library (BML) and Parallel, Rapid O(N), and Graph-based Recursive Electronic Structure Solver (PROGRESS) library. The BML implements linear algebra operations necessary for electronic structure kernels using a unified user interface for various matrix formats (dense and sparse) and architectures (CPUs and GPUs). Focusing on density functional theory and tight-binding models, PROGRESS implements several solvers for computing the single-particle density matrix and relies on BML. In this paper, we describe the general strategies used for these implementations on various computer architectures, using OpenMP target functionalities on GPUs, in conjunction with third-party libraries to handle performance critical numerical kernels. In this study, we demonstrate the portability of this approach and its performance in benchmark problems.
Short-distance matrix elements for $D$-meson mixing for 2+1 lattice QCD
We study the short-distance hadronic matrix elements for D-meson mixing with partially quenched N f = 2+1 lattice QCD. We use a large set of the MIMD Lattice Computation Collaboration's gauge configurations with a 2 tadpole-improved staggered sea quarks and tadpole-improved Lüscher-Weisz gluons. We use the a 2 tadpole-improved action for valence light quarks and the Sheikoleslami-Wohlert action with the Fermilab interpretation for the valence charm quark. Our calculation covers the complete set of five operators needed to constrain new physics models for D-meson mixing. We match our matrix elements to the MS-NDR scheme evaluated at 3 GeV. We report values for the Beneke-Buchalla-Greub-Lenz-Nierste choice of evanescent operators.
Investigating the Difference in Mechanical Stability of Retained Austenite in Bainitic and Martensitic High-Carbon Bearing Steels using in situ Neutron Diffraction and Crystal Plasticity Modeling
In situ neutron diffraction of the uniaxial tension test was used to study the effect of the surrounding matrix microstructure on the mechanical stability of retained austenite in high-carbon bearing steels. Comparing the samples with bainitic microstructures to those with martensitic ones, it was found that the retained austenite in a bainitic matrix starts transforming into martensite at a lower strain compared to that within a martensitic matrix. On the other hand, the rate of transformation of the austenite was found to be higher within a martensitic microstructure. Crystal plasticity modeling was used to analyze the transformation phenomenon in these two microstructures and determine the effect of the surrounding microstructure on elastic, plastic, and transformation components of the strain. The results showed that the predominant difference in the deformation accumulated was from the transformation strain and the critical transformation driving force within the two microstructures. The retained austenite was more stable for identical loading conditions in case of martensitic matrix compared to the bainitic one. It was also observed that the initial volume fraction of retained austenite within the bainitic matrix would alter the onset of transformation to martensite, but not the rate of transformation.
Interaction-Driven Spontaneous Quantum Hall Effect on a Kagome Lattice
Topological states of matter have been widely studied as being driven by an external magnetic field, intrinsic spin-orbital coupling, or magnetic doping. With this work, we unveil an interaction-driven spontaneous quantum Hall effect (a Chern insulator) emerging in an extended fermion-Hubbard model on a kagome lattice, based on a state-of-the-art density-matrix renormalization group on cylinder geometry and an exact diagonalization in torus geometry.We first demonstrate that the proposed model exhibits an incompressible liquid phase with doublet degenerate ground states as time-reversal partners. The explicit spontaneous time-reversal symmetry breaking is determined by emergent uniform circulating loop currents between nearest neighbors. Moreover, the fingerprint topological nature of the ground state is characterized by quantized Hall conductance. Hence, we identify the liquid phase as a quantum Hall phase, which provides a “proof-of-principle” demonstration of the interaction-driven topological phase in a topologically trivial noninteracting band.
Unifying soil organic matter formation and persistence frameworks: the MEMS model
Soil organic matter (SOM) dynamics in ecosystem-scale biogeochemical models have traditionally been simulated as immeasurable fluxes between conceptually defined pools. This greatly limits how empirical data can be used to improve model performance and reduce the uncertainty associated with their predictions of carbon (C) cycling. Recent advances in our understanding of the biogeochemical processes that govern SOM formation and persistence demand a new mathematical model with a structure built around key mechanisms and biogeochemically relevant pools. Here, we present one approach that aims to address this need. Our new model (MEMS v1.0) is developed from the Microbial Efficiency-Matrix Stabilization framework, which emphasizes the importance of linking the chemistry of organic matter inputs with efficiency of microbial processing and ultimately with the soil mineral matrix, when studying SOM formation and stabilization. Building on this framework, MEMS v1.0 is also capable of simulating the concept of C saturation and represents decomposition processes and mechanisms of physico-chemical stabilization to define SOM formation into four primary fractions. After describing the model in detail, we optimize four key parameters identified through a variance-based sensitivity analysis. Optimization employed soil fractionation data from 154 sites with diverse environmental conditions, directly equating mineral-associated organic matter and particulate organic matter fractions with corresponding model pools. Finally, model performance was evaluated using total topsoil (0–20 cm) C data from 8192 forest and grassland sites across Europe. Despite the relative simplicity of the model, it was able to accurately capture general trends in soil C stocks across extensive gradients of temperature, precipitation, annual C inputs and soil texture. The novel approach that MEMS v1.0 takes to simulate SOM dynamics has the potential to improve our forecasts of how soils respond to management and environmental perturbation. Ensuring these forecasts are accurate is key to effectively informing policy that can address the sustainability of ecosystem services and help mitigate climate change.
Pairing tendencies in a two-orbital Hubbard model in one dimension
The recent discovery of superconductivity under high pressure in the ladder compound BaFe 2 S 3 has opened a new field of research in iron-based superconductors with focus on quasi-one-dimensional geometries. In this publication, using the density matrix renormalization group technique, we study a two-orbital Hubbard model defined in one-dimensional chains. Our main result is the presence of hole binding tendencies at intermediate Hubbard U repulsion and robust Hund coupling J H / U = 0.25. Binding does not occur either in weak coupling or at very strong coupling. The pair-pair correlations that are dominant near half-filling, or of similar strength as the charge and spin correlation channels, involve hole-pair operators that are spin singlets, use nearest-neighbor sites, and employ different orbitals for each hole. As a result, the Hund coupling strength, presence of robust magnetic moments, and antiferromagnetic correlations among them are important for the binding tendencies found here.
Evaluation of the cytotoxic and cellular proteome impacts of food-grade TiO2 (E171) using simulated gastrointestinal digestions and a tri-culture small intestinal epithelial model
Engineered nanomaterials (ENMs) are widely used in the food industry; however, regulations for ENMs in food are still in the early stages of development due to insufficient health data. This study investigated the cytotoxicity and changes to the proteomic profile in an in vitro small intestinal epithelium model after exposure to digested food models containing the ubiquitous engineered particulate food additive, TiO2 (E171). TiO2 at 0.75% or 1.5% (w/w) concentrations in either a fasting food model (FFM) or a standardized food model (SFM) based on American diet were digested using an in vitro oral-gastric-small intestinal simulator, and the resulting digestas were applied to a small intestinal epithelium tri-culture cellular model. Effects on cell layer integrity, cytotoxicity, and oxidative stress were assessed. In order to explore impact on cellular processes beyond basic cytotoxicity, mass spectrometrybased quantitative proteomic analyses of control and exposed tri-culture cells was performed. TiO2 in FFM, but not in SFM, produced significant, dose-dependent cytotoxicity, and at the higher dose caused significant oxidative stress, indicative of a food matrix effect. No significant perturbations of the cellular proteome were observed with TiO2 in either FFM or SFM food models. However, proteins involved in the tricarboxylic acid (TCA) cycle, gluconeogenesis, and translation were up-regulated by digestas from SFM compared to those from FFM, also indicative of food matrix effects in the case of food models alone. Interestingly, the differences in profiles between the two food models was more pronounced in the presence of TiO2 indicative of food matrix effects. Together, these results indicate that TiO2 in a fasting diet may be slightly cytotoxic, and that ingested TiO2 does not significantly alter the epithelial proteome, whereas the food matrix alone can have a dramatic effect on the proteome.
Plasma production and ion heating systems for the Material Plasma Exposure eXperiment
Plasma production and ion heating in the Material Plasma Exposure eXperiment (MPEX), whose design is nearing completion, is accomplished using continuous wave RF power with average power density up to 1.6 MW/m2 at the plasma interface. Plasma is produced using helicon waves coupled through a single helical antenna at 13.56 MHz with power supplied by three 100 kW fixed-frequency RF generators feeding a power combiner network followed by the matching network and launcher. The helicon source can utilize various gasses including hydrogen, deuterium, and helium, with a magnetic field strength in the source region up to 0.2 T, and maximum device |B| of 2.5 T. Power is coupled to ions via ion cyclotron heating at the fundamental resonance using a pair of phased helical antennas operating in the frequency range 4-9 MHz, that launch waves towards the resonance from the antenna region where ω > ωci. ICH power is supplied by a single 500 kW tunable RF transmitter through a 90° power splitter and matching/decoupling network. In the case of both the helicon and ICH systems the antennas are located external to the vacuum, with power transferred through novel water-cooled coaxial vacuum windows consisting of fused quartz outer cylinders and silicon nitride inner cylinders with forced convection water cooling between them. The antenna enclosures are pressurized with dry air to 3 bar absolute for voltage standoff.Several 3-D COMSOL models have been created to simulate the two systems. A model of the helicon region utilizing a cold plasma dielectric tensor with accurate magnetic field and realistic plasma density profiles has been used to calculate the plasma loading/complex antenna input impedance at the launcher feed for various ne and |B| values, for the purpose of estimating power handling. It also incorporates the geometry of all launcher structures relevant to this determination. A still more detailed model of the launcher together with a lossy dielectric plasma surrogate has been used to determine RF electric field values and power losses in the device components.Similar models have been produced to predict the performance and power handling of the ICH launcher. For this device the impedance matrix of the two-element antenna array is calculated using a warm plasma model, necessary to properly determine the wave propagation and absorption.
Superposition of M{sub 5}X{sub 5} superstructures and X-ray diffraction in TiO{sub 1.0} titanium monoxide
The interpretation of diffraction spectra of ordered high-temperature phases of solid solutions and strongly nonstoichiometric compounds is discussed. It has been shown that variations of the intensities of superstructure reflections, which cannot be explained within simple ordering models, can be due to the superposition of superstructures with different symmetries in the matrix of the basis crystal structure. Using an example of atom–vacancy ordering in TiO{sub 1.0} titanium monoxide, a model of the order–order transition state formed by the superposition of low-temperature monoclinic (space group A2/m (C2/m)) and high-temperature cubic (space group Pm3̅m) M{sub 5}X{sub 5} superstructures has been proposed. It has been shown that the transition state is thermodynamically equilibrium and should be implemented instead of the M{sub 5}X{sub 5} cubic superstructure. The transition state model can be considered as an M{sub (5–i)}X{sub (5–i)} superstructure (i = 1, 14/18, 11/18) with the monoclinic symmetry (space group P1m1).
On optimizing Jacobi–Davidson method for calculating eigenvalues in low dimensional structures using eight band k · p model
The paper presents two ways of improving the Jacobi–Davidson method for calculating the eigenvalues and eigenvectors described by eight-band k · p model for quantum dots and other low dimensional structures. First, the method is extended by the application of time reversal symmetry operator. This extension allows efficient calculations of the twofold degeneracy present in the multiband k · p model and other interior eigenvalues. Second, the preconditioner for the indefinite matrix which comes from the discretization of the eight band k · p Hamiltonian is presented. The construction of this preconditioner is based on physical considerations about energy band structure in the k · p model. On the basis of two real examples, it is shown that the preconditioner can significantly shorten the time needed to calculate the interior eigenvalues, despite the fact that the memory usage of the preconditioner and Hamiltionian is comparable. Finally, some technical details for implementing the eight band k · p Hamiltonian and the eigensolver are provided.