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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nodal Admittance Matrix Based Area Partition Method for Small-Signal Stability Analysis of Large-Scale Power Electronics Based Power Systems

In power electronics-based power systems (PEPSs), small-signal stability is an important factor for system design and operation, where the impedance-based approach is often used. However, unlike small-scale PEPSs with simple and straightforward impedance models, it would take much more efforts to derive the large-scale PEPSs impedance model, which is very complicated and sometimes may get wrong results due to the elimination of right-half plane (RHP) poles during the impedance aggregation process. To simplify the derivation procedure and analyze the small-signal stability of large-scale PEPSs, this paper proposes a nodal admittance matrix (NAM) based area partition method. In this method, the large-scale PEPS is divided into several sub-areas, and the stability is analyzed within the sub-area first, and then the interconnection stability among these sub-areas is analyzed. The proposed method is scalable and can help to locate the weakest areas/converters that may cause instability in the whole system. In this paper, the concept of the proposed method and its application to an example system are introduced. Experimental results are also given to validate the effectiveness of the proposed method.

Qiao, Liang↗
DOE OSTI2026

Advanced measurement techniques in quantum Monte Carlo: The permutation matrix representation approach

In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also derive more complicated non-local or even dynamic observable estimators. In contrast, we show that arbitrary static observables can be estimated within the permutation matrix representation (PMR) flavor for any Hamiltonian. We then generalize these results to general imaginary-time correlation functions and non-trivial integrated susceptibilities thereof. Finally, we demonstrate the practical versatility of our method by estimating various non-local, random observables for the transverse-field Ising model on a square lattice and a toy random model.

Permutation matrix representation↗
DOE OSTI2023

High Fidelity Computational Model for Fluidized Bed Experiments (Final Report)

MFiX was integrated with preconditioners and linear solver packages in Trilinos via MFIX, Fortran, C, and C++ wrappers. The MFIX wrapper interpreted the structure of the matrix and vector from MFiX while the Fortran wrapper transfers this information to a low-level language, C. C-wrapper transforms the memory semantic between Fortran and C++ language. C++ wrapper passes the matrix and vector to a Trilinos package and gets the solution. This solution was transferred to MFiX via C++, C, Fortran, MFIX wrapper. The framework is used to integrate first, second and third generation linear solvers in Trilinos with MFiX.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Product Consistency Test Results for the High-Chromium Matrix Glasses

This report provides chemical analysis of a series of Product Consistency Test (PCT) leachates from simulated nuclear waste glasses. The resulting data will be used in the development of improved property/composition models. For some of the glass leachates, minor scatter among the triplicate values of some analytes was observed. For other leachates, there were more significant differences among the triplicate values. A review of the PCT data noted that there was little difference between the normalized values on the basis of targeted or measured glass composition. Several of the study glasses, both quenched and canister centerline cooled, have NC i values that are greater than the Hanford Tank Waste Treatment and Immobilization Plant immobilized low-activity waste constraint of 4 g/L for B, Na, and Si. For many of the study glasses, heat treatment had only a marginal impact on the NC i values. The samples of the Environmental Assessment (EA) reference glass included with each PCT set had inconsistent NC i values, although some degree of variation in NC i values for the EA glass is typical. The release rates for boron and sodium were highly correlated for the study glasses, while for the other analytes, less correlation was observed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A constitutive model for glass–ceramic materials

Glass–ceramics have received recent attention for use in glass–ceramic to metal hermetic seals. Due to their heterogeneous microstructure, these materials exhibit a number of advantageous responses over conventional glass based seals. Key amongst them is the possibility of a controllable thermal strain response and apparent coefficient of thermal expansion which may be used to minimize thermally induced residual stresses for aforementioned seals. Furthermore, these behaviors result from an inorganic glass matrix and variety of crystalline ceramic phases including silica polymorph(s) that may undergo reversible solid-to-solid transformations with associated inelastic strain. Correspondingly, these materials exhibit complex thermomechanical responses associated with multiple inelastic mechanisms (viscoelasticity and phase transformation). While modeling these behaviors is essential for developing and analyzing the corresponding applications, no such model exists. Therefore, in this work a three-dimensional continuum constitutive model for glass–ceramic materials combining these various inelastic mechanisms is developed via an internal state variable approach. A corresponding fully implicit three dimensional numerical formulation is also proposed and implemented. The model is used to simulate existing experiments and validate the proposed formalism. As an example, the simple seal problem of a glass–ceramic seal inside a concentric metal shell is explored. Finally, the impact of the cooling rates, viscoelastic shift factors, and inelastic strain on final residual stress state are all investigated and the differing contributions highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Status of Lattice QCD Calculations of One and Two-Nucleon Matrix Elements

Lattice QCD calculations of basic hadronic quantities are now often performed with systematic control over the physical pion mass, continuum and infinite volume limits. This progress signals a new era in which lattice QCD will be used to compute basic properties of hadrons, nucleons and light nuclei directly from the Standard Model (SM). In particular, we will be able to use lattice QCD, combined with Effective Field Theory, to make quantitative statements about the interaction of SM matter and potential beyond the SM physics ranging from direct dark matter detection to permanent electric dipole moments. I will briefly review the status and challenges of lattice QCD calculations of select one- and two-nucleon matrix elements. These calculations are significantly more challenging than the basic quantities computed with all systematics controlled, but substantial progress is being made.

Walker-Loud, Andre↗
DOE OSTI2024

Jacobian-based Model Diagnostics and Application to Equation Oriented Modeling of a Carbon Capture System

Equation-oriented (EO) modeling has the potential to enable the effective design and optimization of the operation of advanced energy systems. However, advanced modeling of energy systems results in a large number of variables and non-linear equations, and it can be difficult to search through these to identify the culprit(s) responsible for convergence issues. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms so they can be rescaled. A further singular value decomposition can be per-formed to identify degenerate sets of equations and remaining scaling issues. This work presents an EO model of a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. The IDAES diagnostics tools were successfully applied to this flowsheet to identify problems to improve model robustness and enable the optimization of process design and operating conditions of a carbon capture system.

Allan, Douglas↗
DOE OSTI2025

Analytical and EZmock covariance validation for the DESI 2024 results

The estimation of uncertainties in cosmological parameters is an important challenge in Large-Scale-Structure (LSS) analyses. For standard analyses such as Baryon Acoustic Oscillations (BAO) and Full-Shape two approaches are usually considered. First: analytical estimates of the covariance matrix use Gaussian approximations and (nonlinear) clustering measurements to estimate the matrix, which allows a relatively fast and computationally cheap way to generate matrices that adapt to an arbitrary clustering measurement. On the other hand, sample covariances are an empirical estimate of the matrix based on an ensemble of clustering measurements from fast and approximate simulations. While more computationally expensive due to the large amount of simulations and volume required, these allow us to take into account systematics that are impossible to model analytically. In this work we compare these two approaches in order to enable DESI's key analyses. We find that the configuration space analytical estimate performs satisfactorily in BAO analyses and its flexibility in terms of input clustering makes it the fiducial choice for DESI's 2024 BAO analysis. On the contrary, the analytical computation of the covariance matrix in Fourier space does not reproduce the expected measurements in terms of Full-Shape analyses, which motivates the use of a corrected mock covariance for DESI's 2024 Full Shape analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Power capture and power take-off load of a self-balanced dual-flap oscillating surge wave energy converter

Wave energy converters are an important part of future renewable energy infrastructure. Predicting their power matrix, capture width ratio, and power take-off loads at a targeted site is required for performance assessment before deployment. Because their testing is very expensive, numerical modeling and simulations play a significant role in those assessments. Linear potential flow theory has limited accuracy under large amplitude wave forcing. More accurate predictions can be obtained by using higher-fidelity models, which are computationally expensive. We present a framework for multi-fidelity numerical simulations to determine the hydrodynamic response, wave capture capability, and power take-off load of a full-scale dual-flap oscillating surge wave energy converter. This design exploits out-of-phase motion by setting the distance between the flaps to half the wavelength of the most occurring wave. The simulations are validated using a 1:10 model experiments in a wave tank. Based on these validations, it was determined that Euler simulations provide an acceptable prediction with 90% reduction in computational time with only 11% error. Utilizing Euler simulations at full-scale, the results demonstrate that the annual electrical energy output is 1.79 GWh under regular wave conditions. Here, one significant improvement over single-flap designs is the capture width ratio which exceeds unity.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Pairing properties of the $t$–$t$'–$t$"–$J$ model

We study the pairing properties of the two-dimensional t–t'–t"–J model, where t' and t" are second and third neighbor hoppings, at a doping level x ≈ 0.1. Recent studies of the t–t'–J model find strong pairing for t' > 0, associated with electron doping, but an absence of pairing for t'<0, associated with hole doping. This is in contrast to the cuprates, where the highest transition temperatures appear for hole doping. Model parametrizations for the cuprates estimate a t" comparable to t', which, in principle, might fix this discrepancy. However, we find that it does not; we observe a suppression of pairing for the hole-doped system (t' < 0,t" > 0), while for the electron-doped system (t' > 0,t" < 0) d-wave pairing is robust. As a result, extended hoppings appear to be insufficient to make the one-band t–t'–J model capable of describing the pairing in the hole-doped system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Engineering Graphene-Ceramic 3D Composite Foams by Freeze Drying

A 3D graphene foam produced by the chemical vapor deposition (CVD) technique is recognized as an effective nanofiller material, as it does not agglomerate in the matrix. Though CVD facilitates pristine graphene foam production, the method poses limitations in developing large-scale 3D graphene composite foams as advanced engineering materials. In this paper, a freeze-drying (FD) process is used to produce a composite 3D foam, which is a mixture of graphene nanoplatelet (GNP) and a low-temperature co-fired ceramic (LTCC). The freeze-dried GNP-LTCC reticulated 3D composite foam has an average pore size ranging from 70 to 100 μm. Pore size is varied by regulating the heat transfer rate during the freezing process using a thermally conductive aluminum (Al) mold and a thermally insulating acrylonitrile butadiene styrene (ABS) mold. Computational thermal modeling is used to visualize the heat transfer and its effect on foam pore size. Subsequently, a freeze-dried GNP-LTCC composite foam is embedded into the LTCC matrix to form a hierarchical assembly by a spark plasma technique without compromising the 3D structure of the FD foam. This study established that a simple, eco-friendly, and scalable processing methodology can produce advanced graphene-based 3D composite foams as the future tailorable nanofillers for designing multimatrix materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluation of data driven low-rank matrix factorization for accelerated solutions of the Vlasov equation

Low-rank methods have shown success in accelerating simulations of a collisionless plasma described by the Vlasov equation, but still rely on computationally costly linear algebra every time step. We propose a data-driven factorization method using artificial neural networks, specifically with convolutional layer architecture, that trains on existing simulation data. At inference time, the model outputs a low-rank decomposition of the distribution field of the charged particles, and we demonstrate that this step is faster than the standard linear algebra technique. Numerical experiments show that the method achieves comparable reconstruction accuracy for interpolation tasks, generalizing to unseen test data in a manner beyond just memorizing training data; patterns in factorization also inherently followed the same numerical trend as those within algebraic methods (e.g., truncated singular-value decomposition). However, when training on the first 70% of a time-series data and testing on the remaining 30%, the method fails to meaningfully extrapolate. Despite this limiting result, the technique may have benefits for simulations in a statistical steady-state or otherwise showing temporal stability. These results suggest that while the model offers a computationally efficient alternative for datasets with temporal stability, its current formulation is best suited for interpolation rather than for predicting future states in time-evolving systems. This study thus lays the groundwork for further refinement of neural network-based approaches to low-rank matrix factorization in high-dimensional plasma simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Prompt-delayed γ-ray spectroscopy of neutron-rich 119.121 In isotopes

The Z = 50 shell closure, near N = 82, is unique in the sense that it is the only shell closure with the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , enclosing the magic gap. The interaction of the proton hole/particle in the above-mentioned orbitals with neutrons in the νh 11/2 orbital is an important prerequisite to the understanding of the nuclear structure near N = 82 and the νπ interaction. To explore the structural similarity between the high-spin isomeric states in In ( Z = 49), Sn ( Z = 50), and Sb ( Z = 51) isotopes from a microscopic point of view. In addition, to understand the role of a proton hole or particle in the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , respectively, with neutron holes in the ν h 11 / 2 orbital on these aforementioned isomers. The fusion and transfer induced fission reaction Be 9 ( U 238 , f ) with 6.2 MeV/u beam energy, using a unique setup consisting of AGATA, VAMOS++, and EXOGAM detectors, was used to populate through the fission process and study the neutron-rich In 119 , 121 isotopes. This setup enabled the prompt-delayed γ -ray spectroscopy of isotopes in the time range of 100 ns – 200 μ s . In the odd- A In 119 , 121 isotopes, indications of a short half-life 19/2 - isomeric state, in addition to the previously known 25/2 + isomeric state, were observed from the present data. Further, new prompt transitions above the 25/2 + isomer in In 121 were identified along with reevaluation of its half-life. The experimental data were compared with the theoretical results obtained in the framework of large-scale shell-model calculations in a restricted model space. Furthermore, the ( π g 9 / 2 ν h 11 / 2 ; I | H ^ | π g 9 / 2 ν h 11 / 2 ; I ) two-body matrix elements of residual interaction were modified to explain the excitation energies and the B ( E 2) transition probabilities in the neutron-rich In isotopes. The (i) decreasing trend of E (29/2 + ) - E (25/2 + ) in odd-In (with dominant configuration π g 9 / 2 - 1 ν h 11 / 2 - 2 and maximum aligned spin of 29/2 + ) and (ii) increasing trend of E (27/2 + ) - E (23/2 + ) in odd-Sb (with dominant configuration π g 7 / 2 + 1 ν h 11 / 2 - 2 and maximum aligned spin of 27/2 + ) with increasing neutron number could be understood as a consequence of hole-hole and particle-hole interactions, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗
DOE OSTI2024

Hybrid analysis of radiative corrections to neutron decay with current algebra and effective field theory

We introduce a useful framework for high-precision studies of the neutron beta decay by merging the current algebra description and the fixed-order effective field theory calculation of the electroweak radiative corrections to the neutron axial form factor. We discuss the advantages of this hybrid method and show that it only requires a minimal amount of lattice QCD inputs to achieve a 10 –4 theory accuracy for the Standard Model prediction of the neutron lifetime and the axial-to-vector coupling ratio λ, both important to the search for physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Measurement of the t-channel single top quark cross section in proton-proton collisions at $\sqrt{s}=5.02$ TeV

The single top quark t-channel production cross section is measured in proton-proton collisions at the CERN LHC at $\sqrt{s}=5.02$ TeV, using data recorded with the CMS detector in 2017, corresponding to an integrated luminosity of 302 pb −1 , and resulting in the first CMS measurement of the process at that energy. Events with one electron or muon and two or more jets, among which at least one is identified as originating from a b quark fragmentation, are analyzed. The combined cross section of single top quark (tq) and single top antiquark ($\overline{\mathrm{t}}\mathrm{q}$) production is $\sigma \left(\mathrm{tq}+\overline{\mathrm{t}}\mathrm{q}\right)={25.5}_{-3.5}^{+3.6}{\left(\mathrm{stat}\right)}_{-3.9}^{+4.2}\left(\mathrm{syst}\right)\pm 0.5\left(\mathrm{lumi}\right)$ pb. The individual cross sections are measured to be $\sigma \left(\mathrm{tq}\right)={17.7}_{-2.7}^{+2.9}{\left(\mathrm{stat}\right)}_{-2.4}^{+2.6}\left(\mathrm{syst}\right)\pm 0.3\left(\mathrm{lumi}\right)$ pb and $\sigma \left(\overline{\mathrm{t}}\mathrm{q}\right)={6.7}_{-1.6}^{+2.4}{\left(\mathrm{stat}\right)}_{-2.5}^{+2.1}\left(\mathrm{syst}\right)\pm 0.1\left(\mathrm{lumi}\right)$ pb. Their ratio is measured to be ${\mathcal{R}}_{\mathrm{t}-\mathrm{ch}}={2.7}_{-0.8}^{+1.5}{\left(\mathrm{stat}\right)}_{-0.3}^{+1.2}\left(\mathrm{syst}\right)$. The absolute value of the Cabibbo-Kobayashi-Maskawa matrix element is found to be |f LV V tb | = 0.92 ± 0.09(exp) ± 0.01(theo). The measurements are in good agreement with the standard model predictions at next-to-next-to-leading order accuracy in quantum chromodynamics.

Hadron-Hadron Scattering↗