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 1,765 records · Page 98

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Almost sure convergence in quantum spin glasses

Recently, Keating, Linden, and Wells [Markov Processes Relat. Fields 21(3), 537-555 (2015)] showed that the density of states measure of a nearest-neighbor quantum spin glass model is approximately Gaussian when the number of particles is large. The density of states measure is the ensemble average of the empirical spectral measure of a random matrix; in this paper, we use concentration of measure and entropy techniques together with the result of Keating, Linden, and Wells to show that in fact the empirical spectral measure of such a random matrix is almost surely approximately Gaussian itself with no ensemble averaging. We also extend this result to a spherical quantum spin glass model and to the more general coupling geometries investigated by Erdős and Schröder [Math. Phys., Anal. Geom. 17(3-4), 441–464 (2014)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Interpretable and flexible non-intrusive reduced-order models using reproducing kernel Hilbert spaces

This paper develops an interpretable, non-intrusive reduced-order modeling technique using regularized kernel interpolation. Existing non-intrusive approaches approximate the dynamics of a reduced-order model (ROM) by solving a data-driven least-squares regression problem for low-dimensional matrix operators. Our approach instead leverages regularized kernel interpolation, which yields an optimal approximation of the ROM dynamics from a user-defined reproducing kernel Hilbert space. We show that our kernel-based approach can produce interpretable ROMs whose structure mirrors full-order model structure by embedding judiciously chosen feature maps into the kernel. The approach is flexible and allows a combination of informed structure through feature maps and closure terms via more general nonlinear terms in the kernel. We also derive a computable a posteriori error bound that combines standard error estimates for intrusive projection-based ROMs and kernel interpolants. In conclusion, the approach is demonstrated in several numerical experiments that include comparisons to operator inference using both proper orthogonal decomposition and quadratic manifold dimension reduction.

Data-driven model reduction↗

DOE OSTI2015

Crustal and uppermost mantle structure in the Middle East: assessing constraints provided by jointly modelling Ps and Sp receiver functions and Rayleigh wave group velocity dispersion curves

Seismic velocity models are found, along with uncertainty estimates, for 11 sites in the Middle East by jointly modelling Ps and Sp receiver functions and surface (Rayleigh) wave group velocity dispersion. The approach performs a search for models that satisfy goodness-of-fit criteria guided by a variant of simulated annealing and uses statistical tools to assess these products of searches. These tools, a parameter correlation matrix and marginal posterior probability density (PPD) function, allow us to evaluate quantitatively the constraints that each data type imposes on model parameters and to identify portions of each model that are well-constrained relative to other portions. This joint modelling technique, which we call ‘multi-objective optimization for seismology’, does not require a good starting solution, although such a model can be incorporated easily, if available, and can reduce the computation time significantly. Applying the process described above to broadband seismic data reveals that crustal thickness varies from 15 km beneath Djibouti (station ATD) to 45 km beneath Saudi Arabia (station RAYN). A pronounced low velocity zone for both Vp and Vs is present at a depth of ~12 km beneath station KIV located in northern part of greater Caucasus, which may be due to the presence of a relatively young volcano. Similarly, we also noticed a 6-km-thick low velocity zone for Vp beginning at 20 km depth beneath seismic station AGIN, on the Anatolian plateau, while positive velocity gradients prevail elsewhere in eastern Turkey. Beneath station CSS, located in Cyprus, an anomalously slow layer is found in the uppermost mantle, which may indicate the presence of altered lithospheric material. Crustal P- and S-wave velocities beneath station D2, located in the northeastern portion of central Zagros, range between 5.2–6.2 and 3.2–3.8 km s -1 , respectively. Finally, in Oman, we find a Moho depth of 34.0 ± 1.0 km and 25.0 ± 1.0 to 30.0 ± 1.0 km beneath stations S02 and S04, respectively.

58 GEOSCIENCES↗

DOE OSTI2014

Verification and Validation of Carbon-Fiber Laminate Low Velocity Impact Simulations.

Presented is a model verification and validation effort using low - velocity impact (LVI) of carbon fiber reinforced polymer laminate experiments. A flat cylindrical indenter impacts the laminate with enough energy to produce delamination, matrix cracks and fiber breaks. Included in the experimental efforts are ultrasonic scans of the damage for qualitative validation of the models. However, the primary quantitative metrics of validation are the force time history measured through the instrumented indenter and initial and final velocities. The simulations, whi ch are run on Sandia's Sierra finite element codes , consist of all physics and material parameters of importance as determined by a sensitivity analysis conducted on the LVI simulation. A novel orthotropic damage and failure constitutive model that is cap able of predicting progressive composite damage and failure is described in detail and material properties are measured, estimated from micromechanics or optimized through calibration. A thorough verification and calibration to the accompanying experiment s are presented. Specia l emphasis is given to the four - point bend experiment. For all simulations of interest, the mesh and material behavior is verified through extensive convergence studies. An ensemble of simulations incorporating model parameter unc ertainties is used to predict a response distribution which is then compared to experimental output. The result is a quantifiable confidence in material characterization and model physics when simulating this phenomenon in structures of interest.

English, Shawn Allen↗

DOE OSTI2021

An initial exploration of Bayesian model calibration for estimating the composition of rocks and soils on Mars

The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Modeling the Long-Term Behavior of Cementitious Barriers - 18188

Cementitious materials are commonly included as backfill/barrier materials in geological repositories for radioactive waste. An understanding of the long-term behavior of cement-based barriers is therefore required to provide confidence in their safety functions within an environmental safety case. Geochemical modelling has been undertaken on several different radioactive waste repository concepts (and their natural/industrial analogues), to explore the nature of cement degradation over timescales relevant to postclosure safety (many tens of thousands of years to c. 1 Ma). In radionuclide transport calculations undertaken as part of performance assessments, it is typically considered that the properties of cementitious engineered barrier materials change as they degrade, potentially resulting in an increase in porosity hydraulic conductivity, and effective diffusion coefficients. However, reactive transport modelling (supported by some analogue data) suggests that, in some repository systems, a degree of clogging of the cement pore space could occur and that some advective flow paths (cracks in the barriers, fractures in the host rock) could become (to some extent) infilled by cement and rock alteration products. This could result in a reduction in water flow and contaminant transport in the vicinity of the cement-rock interface, but some uncertainties remain (especially process couplings and the nature of water flow in complex fracture networks that are becoming clogged). In addition, the armoring of fracture surfaces could result in decreased potential for diffusion of solutes into the rock matrix. It is envisaged that ongoing work will continue to improve reactive-transport model capabilities, such that the implications cement degradation has for contaminant transport can be predicted with greater confidence. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Phenomenology of the Prethermal Many-Body Localized Regime

The dynamical phase diagram of interacting disordered systems has seen substantial revision over the past few years. Theory must now account for a large prethermal many-body localized regime in which thermalization is extremely slow, but not completely arrested. We derive a quantitative description of these dynamics in short-ranged one-dimensional systems using a model of successive many-body resonances. Here, the model explains the decay timescale of mean autocorrelators, the functional form of the decay—a stretched exponential—and relates the value of the stretch exponent to the broad distribution of resonance timescales. The Jacobi method of matrix diagonalization provides numerical access to this distribution, as well as a conceptual framework for our analysis. The resonance model correctly predicts the stretch exponents for several models in the literature. Successive resonances may also underlie slow thermalization in strongly disordered systems in higher dimensions, or with long-range interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Decentralized Low-Rank State Estimation for Power Distribution Systems

This article considers the low-observability state estimation problem in power distribution networks and develops a decentralized state estimation algorithm leveraging the matrix completion methodology. Matrix completion has been shown to be an effective technique in state estimation that exploits the low dimensionality of the power system measurements to recover missing information. This technique can utilize an approximate (linear) load flow model, or it can be used with no physical models in a network where no information about the topology or line admittance is available. The direct application of matrix completion algorithms requires solving a semi-definite programming (SDP) problem, which becomes computationally challenging for large networks. We therefore develop a decentralized algorithm that capitalizes on the popular proximal alternating direction method of multipliers (proximal ADMM). The method allows us to distribute the computation among different areas of the network, leading to a scalable algorithm. By doing all computations at individual control areas and only communicating with neighboring areas, the algorithm eliminates the need for data to be sent to a central processing unit and thus increases efficiency and contributes to the goal of autonomous control of distribution networks. We illustrate the advantages of the proposed algorithm numerically using standard IEEE test cases.

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

DOE OSTI2024

Altered post-fracture systemic bone loss in a mouse model of osteocyte dysfunction

Femur fracture leads to loss of bone at uninjured skeletal sites, which may increase risk of subsequent fracture. Osteocytes, the most abundant bone cells, can directly resorb bone matrix and regulate osteoclast and osteoblast activity, but their role in systemic bone loss after fracture remains poorly understood. In this study we used a transgenic (TG+) mouse model that overexpresses human B-cell lymphoma 2 (BCL-2) in osteoblasts and osteocytes. This causes enhanced osteoblast proliferation, followed by disruption in lacunar-canalicular connectivity and massive osteocyte death by 10 wk of age. We hypothesized that reduced viable osteocyte density would decrease the magnitude of systemic bone loss after femur fracture, reduce perilacunar remodeling, and alter callus formation. Bone remodeling was assessed using serum biomarkers of bone formation and resorption at 5 d post-fracture. We used micro-computed tomography, high resolution x-ray microscopy, mechanical testing, and Raman spectroscopy to quantify the magnitude of systemic bone loss, as well as changes in osteocyte lacunar volume, bone strength, and bone composition 2 wk post-fracture. Fracture was associated with a reduction in circulating markers of bone resorption in non-transgenic (TG-) animals. TG+ mice exhibited high bone mass in the limbs, greater cortical elastic modulus and reduced post-yield displacement. After fracture, TG+ mice lost less trabecular bone than TG- mice, but conversely TG+ mice exhibited trends toward a lower yield point and reduced femoral cortical thickness after fracture, though these were not statistically significant. Lacunar density was greater in TG+ mice, but fracture did not alter lacunar volume in TG+ or TG- mice. These findings suggest that osteocytes potentially play a significant role in the post-traumatic systemic response to fracture, though the effects differ between trabecular and cortical bone.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Coexistence of non-Abelian chiral spin liquid and magnetic order in a spin-1 antiferromagnet

We study the ground-state properties of a spin-1 Heisenberg model on a square lattice with the first- and second-nearest-neighbor antiferromagnetic couplings $J_1$ and $J_2$ and a three-spin scalar chirality term $J_\chi$. Using the density matrix renormalization group calculation, we map out a global phase diagram including various magnetic order phases and an emergent quantum spin liquid phase. The nature of the spin liquid is identified as a bosonic non-Abelian Moore-Read state from the fingerprint of the entanglement spectra and identification of a full set of topological sectors. We further unveil a stripe magnetic order coexisting with this spin liquid. Furthermore, our results not only establish a rare example of non-Abelian spin liquids in simple spin systems but also demonstrate the coexistence of fractionalized excitations and magnetic order beyond mean-field descriptions.

2-dimensional systems↗

DOE OSTI2020

Source identification by non-negative matrix factorization combined with semi-supervised clustering

Machine-learning methods and apparatus are provided to solve blind source separation problems with an unknown number of sources and having a signal propagation model with features such as wave-like propagation, medium-dependent velocity, attenuation, diffusion, and/or advection, between sources and sensors. In exemplary embodiments, multiple trials of non-negative matrix factorization are performed for a fixed number of sources, with selection criteria applied to determine successful trials. A semi-supervised clustering procedure is applied to trial results, and the clustering results are evaluated for robustness using measures for reconstruction quality and cluster separation. The number of sources is determined by comparing these measures for different trial numbers of sources. Source locations and parameters of the signal propagation model can also be determined. Disclosed methods are applicable to a wide range of spatial problems including chemical dispersal, pressure transients, and electromagnetic signals, and also to non-spatial problems such as cancer mutation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Source identification by non-negative matrix factorization combined with semi-supervised clustering

Alexandrov, Boian S.↗

DOE OSTI2026

Simplex‐based model for nanoparticle grain identification in four‐dimensional scanning transmission electron microscopy data

Grain identification in polycrystalline nanoparticles, for example, determining which crystal phases are present at each spatial location, is fundamental to materials characterisation. This is particularly challenging when grains overlap extensively, as commonly occurs in four-dimensional scanning transmission electron microscopy (4D-STEM) datasets. We propose a simplex-based model (SBM) in which each simplex vertex represents the diffraction pattern (DP) of a pure grain, and the simplex edges and interior represent overlapping grains. Our SBM grain identification algorithm operates on the Bragg disk (BD) data matrix distilled from the 4D-STEM data to identify the grain membership at each scan position, together with a BD feature matrix whose columns represent the DPs for each constituent grain, which is important for identifying the crystal structure of each grain. We solve the model using a two-stage algorithm. In Stage 1, we adapt a linear mixing algorithm to estimate an initial BD feature matrix whose columns represent DPs of potentially overlapping grains. Our Stage 2 algorithm incorporates sparsity considerations to transform the initial BD feature matrix so that its columns represent DPs of pure grains. Using simulated datasets with various grain configurations, we demonstrate that SBM recovers both the BD feature matrix and membership maps more accurately than existing methods, even when a grain lacks any pure region and completely overlaps with other grains.

4D-STEM segmentation↗

DOE OSTI2022

An energy-based study of the embedded element method for explicit dynamics

Abstract The embedded finite element technique provides a unique approach for modeling of fiber-reinforced composites. Meshing fibers as distinct bundles represented by truss elements embedded in a matrix material mesh allows for the assignment of more specific material properties for each component rather than homogenization of all of the properties. However, the implementations of the embedded element technique available in commercial software do not replace the material of the matrix elements with the material of the embedded elements. This causes a redundancy in the volume calculation of the overlapping meshes leading to artificially increased stiffness and mass. This paper investigates the consequences in the energy calculations of an explicit dynamic model due to this redundancy. A method for the correction of the edundancy within a finite element code is suggested which removes extra energy and is shown to be effective at correcting the energy calculations for large amounts of redundant volume.

Martin, Valerie A.↗

DOE OSTI2021

Performance of the rigorous renormalization group for first-order phase transitions and topological phases

Expanding and improving the repertoire of numerical methods for studying quantum lattice models is an ongoing focus in many-body physics. While the density matrix renormalization group (DMRG) has been established as a practically useful algorithm for finding the ground state in one-dimensional systems, a provably efficient and accurate algorithm remained elusive until the introduction of the rigorous renormalization group (RRG) by Landau [Nat. Phys. 11, 566 (2015)1745-247310.1038/nphys3345]. In this paper, we study the accuracy and performance of a numerical implementation of RRG at first-order phase transitions and in symmetry-protected topological phases. Our study is motivated by the question of when RRG might provide a useful complement to the more established DMRG technique. In particular, despite its general utility, DMRG can give unreliable results near first-order phase transitions and in topological phases, since its local update procedure can fail to adequately explore (near-)degenerate manifolds. Furthermore, the rigorous theoretical underpinnings of RRG, meanwhile, suggest that it should not suffer from the same difficulties. We show this optimism is justified, and that RRG indeed determines well-ordered, accurate energies even when DMRG does not. Moreover, our performance analysis indicates that in certain circumstances seeding DMRG with states determined by coarse runs of RRG may provide an advantage over simply performing DMRG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Combining the mapping Hamiltonian linearized semiclassical approach with the generalized quantum master equation to simulate electronically nonadiabatic molecular dynamics

The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. Here, the accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Hypernuclear no-core shell model

We extend the No-Core Shell Model (NCSM) methodology to incorporate strangeness degrees of freedom and apply it to single-Λ hypernuclei. After discussing the transformation of the hyperon-nucleon (YN) interaction into Harmonic-Oscillator (HO) basis and the Similarity Renormalization Group transformation applied to it to improve model-space convergence, we introduce two complementary formulations of the NCSM, one that uses relative Jacobi coordinates and symmetry-adapted basis states to fully exploit the symmetries of the hypernuclear Hamiltonian, and one working in a Slater determinant basis of HO states where antisymmetrization and computation of matrix elements is simple and to which an importance-truncation scheme can be applied. For the Jacobi-coordinate formulation, we give an iterative procedure for the construction of the antisymmetric basis for arbitrary particle number and present the formulae used to embed two- and three-baryon interactions into the many-body space. For the Slater-determinant formulation, we report on the conversion of the YN interaction matrix elements from relative to single-particle coordinates, the importance-truncation scheme that tailors the model space to the description of the low-lying spectrum, and the role of the redundant center-of-mass degrees of freedom. We end with a confirmation of both formulations in the four-body system, giving converged ground-state energies for a chiral Hamiltonian, and present a short survey of the A ≤ 7 hyper-helium isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗