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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Microstructural development in DED stainless steels: applying welding models to elucidate the impact of processing and alloy composition

Austenitic stainless steel microstructures produced by directed energy deposition (DED)are analogous to those developed during welding, particularly high energy density welding. To better understand microstructural development during DED, theories of microstructural evolution,which have been established to contextualize weld microstructures, are applied in this study to microstructural development in DED austenitic stainless steels. Phenomenological welding models that describe the development of oxide inclusions, compositional microsegregation, ferrite,matrix austenite grains, and dislocation substructures are utilized to clarify microstructural evolution during deposition of austenitic stainless steels. Two different alloys, 304L and 316L, arecompared to demonstrate the broad applicability of this framework for understanding microstmctural development during the DED process. Despite differences in grain morphology and solidification mode for these two alloys (which can be attributed to compositional differences),similar tensile properties are achieved. It is the fine-scale compositional segregation and dislocation structures that ultimately determine the strength of these materials. The evolution of microsegregation and dislocation structures is shown to be dependent on the rapid solidification and thermomechanical history of the DED processing method and not the composition of the starting material.

36 MATERIALS SCIENCE↗

DOE OSTI2024

On the applicability of the Redlich-Kister framework for viscosity estimation of molten halide salt mixtures

For molten halide salt mixtures already being utilized or under consideration for carbon-free energy production systems, it is crucial that their viscosity is well understood so that system thermal hydraulics can be reliably assessed. Because of the difficulty in accurately measuring molten halide viscosity and the sheer size of the matrix of possible higher order salt mixtures that may be of interest to the energy industry, there are several gaps in the quantified understanding of molten halide viscosity across this matrix. As such, both first-principles and semi-empirical modeling techniques may be crucial for rapidly assessing this broad, complex compositional domain. Herein, the Redlich-Kister framework is applied to assess the feasibility of broadly interpolating and estimating the viscosity of several pseudobinary and pseudoternary molten halide salt systems that may be of key interest to the energy industry. The framework is based on the assumption that an ideal component and a nonideal component collectively describe the viscosity as a function of composition and temperature for a given molten halide system. Three different ideal models were considered for the ideal component, including Grunburg-Nissan, Katti-Chaudhri, and Gambill methods. Regarding the pseudobinary interpolations, the Redlich-Kister models with either the Grunburg-Nissan or Katti-Chaudhri models as the ideal component resulted in either highly (average error less than 5%) or reasonably (average error less than 15%) accurate interpolations of pseudobinary halide viscosity; BeF 2 - or UF 4 -bearing salts tended to result in reasonably accurate interpolations, whereas other pseudobinary mixtures tended to show high accuracy. Regarding the pseudoternary extrapolations, the Redlich-Kister framework shows reasonable success at estimating the extent to which a pseudoternary system may indicate deviations from ideal Grunburg-Nissan mixing, where discrepancies with comparative experimental data generally stay within 30%. Finally, the primary reasons identified for such discrepancies are (1) inaccuracy in the underlying experimental data, (2) different complexation behavior in the higher order systems compared to the pseudobinary subsystems, and (3) extrapolation into temperatures too far out of the domain, which is valid for the underlying experimental data feeding the Redlich-Kister model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Moisture-induced autonomous surface potential oscillations for energy harvesting

A variety of autonomous oscillations in nature such as heartbeats and some biochemical reactions have been widely studied and utilized for applications in the fields of bioscience and engineering. Here, we report a unique phenomenon of moisture-induced electrical potential oscillations on polymers, poly([2-(methacryloyloxy)ethyl] dimethyl-(3-sulfopropyl) ammonium hydroxide-co-acrylic acid), during the diffusion of water molecules. Chemical reactions are modeled by kinetic simulations while system dynamic equations and the stability matrix are analyzed to show the chaotic nature of the system which oscillates with hidden attractors to induce the autonomous surface potential oscillation. Using moisture in the ambient environment as the activation source, this self-excited chemoelectrical reaction could have broad influences and usages in surface-reaction based devices and systems. As a proof-of-concept demonstration, an energy harvester is constructed and achieved the continuous energy production for more than 15,000 seconds with an energy density of 16.8 mJ/cm 2 . A 2-Volts output voltage has been produced to power a liquid crystal display toward practical applications with five energy harvesters connected in series.

42 ENGINEERING↗

DOE OSTI2020

Low-rank Sachdev-Ye-Kitaev models

Motivated by recent proposals of experimental realization of fast scramblers, we study a family of solvable variants of the ($\textit{q}$ = 4) Sachdev-Ye-Kitaev model in which the rank and eigenvalue distribution of the coupling matrix $J_{ ij,kl}$ are tuneable. When the rank is proportional to the number of fermions, the low temperature behavior is sensitive to the eigenvalue distribution. We obtain a complete classification of the possible non-Fermi liquid quantum phases. These include two previously studied phases whose fermion scaling dimension depends continuously on the rank; we show that they are maximally chaotic, but necessitate an extensively degenerate or negative semidefinite coupling matrix. In this work, more generic distributions give rise to “almost Fermi liquids” with a scaling dimension Δ = 1/2, but which differ from a genuine Fermi liquid in quasiparticle decay rate, quantum Lyapunov exponent, and/or specific heat.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Extracting off-diagonal order from diagonal basis measurements

Quantum gas microscopy has developed into a powerful tool to explore strongly correlated quantum systems. However, discerning phases with topological or off-diagonal long range order requires the ability to extract these correlations from site-resolved measurements. Here, we show that a multiscale complexity measure can pinpoint the transition to and from the bond ordered wave phase of the one-dimensional extended Hubbard model with an off-diagonal order parameter, sandwiched between diagonal charge and spin density wave phases, using only diagonal descriptors. We study the model directly in the thermodynamic limit using the recently developed variational uniform matrix product states algorithm, and draw our samples from degenerate ground states related by global spin rotations, emulating the projective measurements that are accessible in experiments. Our results will have important implications for the study of exotic phases using optical lattice experiments. Published by the American Physical Society 2024

1-dimensional systems↗

DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Shedding Light on Primary Donors in Photosynthetic Reaction Centers

Chlorophylls (Chl)s exist in a variety of flavors and are ubiquitous in both the energy and electron transfer processes of photosynthesis. The functions they perform often occur on the ultrafast (fs–ns) time scale and until recently, these have been difficult to measure in real time. Further, the complexity of the binding pockets and the resulting protein-matrix effects that alter the respective electronic properties have rendered theoretical modeling of these states difficult. Recent advances in experimental methodology, computational modeling, and emergence of new reaction center (RC) structures have renewed interest in these processes and allowed researchers to elucidate previously ambiguous functions of Chls and related pheophytins. This is complemented by a wealth of experimental data obtained from decades of prior research. Studying the electronic properties of Chl molecules has advanced our understanding of both the nature of the primary charge separation and subsequent electron transfer processes of RCs. In this review, we examine the structures of primary electron donors in Type I and Type II RCs in relation to the vast body of spectroscopic research that has been performed on them to date. Further, we present density functional theory calculations on each oxidized primary donor to study both their electronic properties and our ability to model experimental spectroscopic data. This allows us to directly compare the electronic properties of hetero- and homodimeric RCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Dark Energy Survey Supernova 5YR Data Release

This data set was generated by the Dark Energy Survey (DES) team, and was used for the cosmological analysis in the DES key supernova cosmology paper (arXiv:2401.02929). When using this data please cite both the data release paper (Sanchez et al. 2024, ApJ, arXiv:2406.05046) and the DES key paper (Dark Energy Survey Collaboration 2024, ApJL, arXiv:2401.02929). At DES-SN5YR Github Repo (https://github.com/des-science/DES-SN5YR), we provide a comprehensive list of references for the different data products (SALT3 light-curve model, calibration, DCR corrections, SN classification, SN distances and uncertainties covariance matrix). If you have questions, post an issue on the DES-SN5YR Github repository and/or get in touch with the contact people indicated on the Github repository.

Contractor, Fermilab [Fermilab]↗

DOE OSTI2020

Improvement of SCALE-XSPROC Multigroup Cross Section Processing Based on the CENTRM Pointwise Slowing Down Calculation

The SCALE-XSProc multigroup (MG) cross section processing procedure based on the CENTRM pointwise slowing down calculation is the primary procedure to process problem-dependent self-shielded MG cross sections and scattering matrices for neutron transport calculations.This procedure supports various cell-based geometries including slab, 1-D cylindrical, 1-D spherical and 2-D rectangular configurations and doubly heterogeneous particulate fuels. Recently, this procedure has been significantly improved to be applied to any advanced reactor analysis covering thermal and fast reactor systems, and to be comparable to continuous energy (CE) Monte Carlo calculations. Some reactivity bias and reaction rate differences have been observed compared with CE Monte Carlo calculations, and several areas for improvement have been identified in the SCALE-XSProc MG cross section processing: (1) resonance self-shielding calculations within the unresolved resonance range, (2) 10 eV thermal cut-off energy for the free gas model, (3) on-the-fly adjustments to the thermal scattering matrix, (4) normalization of the pointwise neutron flux, and (5) fine MG energy structure. This procedure ensures very accurate MG cross section processing for high-fidelity deterministic reactor physics analysis for various advanced reactor systems.

Kim, Kang Seog↗

DOE OSTI2020

Analyzing the Effects of Cyberattacks on Distribution System State Estimation: Preprint

Key components of power systems—such as energy management systems, automatic generation control, and state estimation—are under serious vulnerability from cyber attacks. Cyber threats in electric grids have increased significantly because of the increased interconnectivity of supervisory control and data acquisition systems and public network infrastructure. As the penetration level of distributed energy resources increases, it is imperative to employ system-monitoring techniques such as state estimation for the reliable operation of distribution systems. Recently, multiple methods have been developed that exploit the low rank property of distribution system state matrix and are robust to bad data, such as matrix completion. This paper analyzes the impact of various realistic cyber attack scenarios on matrix completion. Realistic cyber attack scenarios are converted into data corruption models that are used in an extensive simulation of a custom IEEE 123-bus system.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2026

Sensitivity of Next-Generation CMB Surveys to Neutrinos and Other Light Relics

Neutrinos and other light relics leave characteristic imprints in the cosmic microwave background anisotropies, making their observation a sensitive probe of the particle content and thermal history of the early universe. The energy density in these relativistic species is parameterized by their effective number $N_\mathrm{eff}$. Measuring this parameter at the percent level, which is a long-standing science goal of CMB-S4 and other experiments, would test a wide range of well-motivated physics within and beyond the Standard Model of particle physics. In this paper, we present Fisher-matrix forecasts of the projected sensitivity to $N_\mathrm{eff}$ of several CMB-S4 survey configurations considered during its extensive design phase. The conceptual design reaches $σ(N_\mathrm{eff}) < 0.03$ over its seven-year observing period, while the revised configuration achieves the same precision over a longer timescale. We complement these results with a cosmic-variance-limited survey over the same multipole range to quantify the room for improvement accessible with additional instrumental, observational, and theoretical efforts. Finally, we discuss the broad implications of precise $N_\mathrm{eff}$ measurements for the radiation sector, big bang nucleosynthesis, light thermal relics, and other early-universe physics. The forecasts presented in this work are performed with the publicly released DRAFT (Dark Radiation Anisotropy Flowdown Team) tool. It provides an end-to-end pipeline from simulated foreground maps and component separation to delensing and projected sensitivities for any cosmological parameter, and it can be directly applied to other cosmic microwave background survey designs.

Trendafilova, Cynthia [Case Western Reserve U., CE↗

DOE OSTI2022

Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl 2 and OsOCl 2

The van der Waals oxide dichlorides MOX 2 (M=V, Ta, Nb, Ru, and Os; X=halogen element), with different electronic densities, are attracting considerable attention. Ferroelectricity, spin-singlet formation, and orbital-selective Peierls phases were reported in this family with d 1 or d 2 electronic configurations, all believed to be caused by the strongly anisotropic electronic orbital degree of freedom. Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl 2 and OsOCl 2 with d 4 electronic configurations. Different from a previous study using VOI 2 with d 1 configuration, these systems with 4d4 or 5d4 do not exhibit a ferroelectric instability along the a axis. Due to the fully occupied dxy orbital in RuOCl 2 and OsOCl 2 , the Peierls instability distortion disappears along the b axis, leading to an undistorted Immm phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the a axis. For this reason, the materials of our focus can be regarded as “effective one-dimensional” systems even when they apparently have a dominant two-dimensional lattice geometry. The large crystal-field splitting energy (between dxz/yz and dxy orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the J = 0 singlet state in MOCl 2 (M=Ru or Os) with electronic density n = 4, resulting in a spin-1 system. Moreover, we find staggered antiferromagnetic order with π wave vector along the M-O chain direction (a axis) while the magnetic coupling along the b axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the t 2g orbitals for the Os atoms to construct a multiorbital Hubbard model for the M-O chains. Staggered AFM with ↑|-|↓-↑|-|↓ spin structure dominates in our density matrix renormalization group calculations, in agreement with density functional theory calculations. Our results for RuOCl 2 and OsOCl 2 provide guidance to experimentalists and theorists working on this interesting family of oxide dichlorides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Triangular lattice Hubbard model physics at intermediate temperatures

Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Enabling fast charging of lithium-ion batteries through secondary- /dual- pore network: Part I - Analytical diffusion model

Battery performance is strongly correlated with electrode microstructural properties. Enabling fast charging of lithium-ion batteries requires improved through-plane ionic diffusion that can be achieved through, among other strategies, structured electrodes with a secondary- or dual-pore network (SPN). In this work, an analytical model investigates the impact of such an SPN on ionic diffusion with a composite electrode, considering various pore-channel geometries and comparing to standard electrodes with identical gravimetric- and volumetric-specific theoretical capacities. Relevant SPN design parameters and tortuosity coefficients are identified according to three optimization objectives that aim to balance the improved overall through-plane diffusion, thanks to the coarse aligned channels, and degraded in-plane diffusion because of the porous matrix densification required to maintain gravimetric- and volumetric-specific theoretical capacities. The model indicates that a relatively low amount of SPN is required and that electrodes with high through-plane tortuosity and low in-plane tortuosity benefit most from such architecture.

25 ENERGY STORAGE↗

DOE OSTI2025

Breakdown of superdiffusion in perturbed quantum integrable spin chains and ladders

Superdiffusive transport with dynamical exponent z=3/2 has been firmly established at finite temperature for a class of integrable systems with a non-Abelian global symmetry G. On the inclusion of integrability-breaking perturbations, diffusive transport with z=2 is generically expected to hold in the limit of late time. Recent studies of the classical Haldane-Ishimori-Skylanin model have found that perturbations that preserve the global symmetry lead to a much slower time scale for the onset of diffusion, albeit with uncertainty over the exact scaling exponent. That is, for perturbations of strength λ, the characteristic timescale for diffusion goes as t*∼λ−α for some α. Using large-scale matrix product state simulations, we investigate this behavior for perturbations to the canonical quantum model showing superdiffusion: the S=1/2 quantum Heisenberg chain. We consider a ladder configuration and look at various perturbations that either break or preserve the SU(2) symmetry, leading to scaling exponents consistent with those observed in one classical study [McCarthy , Phys. Rev. B 110, L180301 (2024)2469-995010.1103/PhysRevB.110.L180301]: α=2 for symmetry-breaking terms and α=6 for symmetry-preserving terms. We also consider perturbations from another integrable point of the ladder model with G=SU(4) and find consistent results. Finally, we consider a generalization to an SU(3) ladder and find that the α=6 scaling appears to be universal across superdiffusive systems when the perturbations preserve the non-Abelian symmetry G.

Wang, Kevin↗

DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗

DOE OSTI2021

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Putting AlphaFold models to work with phenix.process_predicted_model and ISOLDE

AlphaFold has recently become an important tool in providing models for experimental structure determination by X-ray crystallography and cryo-EM. Large parts of the predicted models typically approach the accuracy of experimentally determined structures, although there are frequently local errors and errors in the relative orientations of domains. Importantly, residues in the model of a protein predicted by AlphaFold are tagged with a predicted local distance difference test score, informing users about which regions of the structure are predicted with less confidence. AlphaFold also produces a predicted aligned error matrix indicating its confidence in the relative positions of each pair of residues in the predicted model. The phenix.process_predicted_model tool downweights or removes low-confidence residues and can break a model into confidently predicted domains in preparation for molecular replacement or cryo-EM docking. These confidence metrics are further used in ISOLDE to weight torsion and atom–atom distance restraints, allowing the complete AlphaFold model to be interactively rearranged to match the docked fragments and reducing the need for the rebuilding of connecting regions.

59 BASIC BIOLOGICAL SCIENCES↗