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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hadronic vacuum polarization for the muon g - 2 from lattice QCD: Complete short and intermediate windows

We present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment a μ in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to a μ , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain a μ l l , SD ( conn ) = 48.139 ( 11 ) stat ( 91 ) syst [ 92 ] total × 10 - 10 and a μ l l , W ( conn ) = 206.90 ( 14 ) stat ( 61 ) syst [ 63 ] total × 10 - 10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are a μ SD = 69.05 ( 1 ) stat ( 21 ) syst [ 21 ] total × 10 - 10 and a μ W = 236.45 ( 17 ) stat ( 83 ) syst [ 85 ] total × 10 - 10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Multi-channel 1-to-2 transition amplitudes in a finite volume

We derive a model-independent expression for finite-volume matrix elements. Specifically, we present a relativistic, non-perturbative analysis of the matrix element of an external current between a one-scalar in-state and a two-scalar out-state. Our result, which is valid for energies below higher-particle inelastic thresholds, generalizes the Lellouch-Luscher formula in two ways: we allow the external current to inject arbitrary momentum into the system and we allow for the final state to be composed an arbitrary number of strongly coupled two-particle states with arbitrary partial waves (including partial-wave mixing induced by the volume). We also illustrate how our general result can be applied to some key examples, such as heavy meson decays and meson photo production. Finally, we point out complications that arise involving unstable resonance states, such as B to K*+l+l when staggered or mixed-action/partially-quenched calculations are performed.

Briceno, Raul↗
DOE OSTI2017

North Slope of Alaska Snow Intensive Operational Period Field Campaign Report

The campaign was motivated by the need to improve the quantification of measurements of ice-phase precipitation in the Arctic and was by the acquisition and deployment of the new X- and Ka/W-band radars. These radars opened up an opportunity for the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility to obtain spatial estimates of snowfall rates using the polarimetric X-band measurements and dual-frequency measurements (using different combinations of the three wavelengths). However, calculations of X- and Ka-band radar back-scattering of ice crystal aggregates with their complex structure suggest that the commonly used T-matrix approach (Matrosov et al. 2007) for modeling the radar back-scattering underestimates the reflectivity by several decibels, with errors increasing with increasing radar frequency (Botta et al. 2010, 2011). Moreover, the X-band polarimetric measurements and the Ka/W-band measurements are sensitive to the assumed shape of the snow (Botta et al. 2011). One of the five ARM two-dimensional video disdrometers (manufactured by Joanneum Research) were deployed in Barrow at the ARM North Slope of Alaska (NSA) site from 1 October, 2011 to 31 May, 2012 in an attempt to use the instrument in a novel way. The instrument was originally designed to measure the drop size distribution of rain but it seemed worthwhile to explore its capability to quantify ice precipitation particle size and shape distributions in the cold north for scattering calculations and precipitation estimations. Furthermore, this deployment gave us an opportunity to see how reliable it could be in arctic conditions.

Ka-band radar↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Study on generating of thermal neutron scattering cross sections for LiH

LiH is designated as a promising moderator and shielding material because of its low density, high melting point and large fraction of H atoms. However, lack of the thermal neutron cross sections of LiH makes numerical calculation deviate from experimental data to some extent. As a result, it is necessary to study LiH thermal kernel effect. The phonon property of LiH has been investigated by first-principles calculations using the plane-wave pseudo potential method with CASTEP code. The scattering law and the thermal neutron scattering cross sections for Li and H have been generated using this distribution. The results have been compared with zirconium hydride data. The GASKET and NJOY/LEAPR codes have been used in the calculation of scattering law, whose results have been compared with the reference; the discrepancy mainly comes from phonon spectrums and its expansion. LEAPR had the capability to compute scattering through larger energy and momentum transfers than GASKET did. By studying LiH phonon spectrum and constructing the model of LiH thermal kernel and scattering matrix, the ACE format LiH thermal neutron cross sections for MCNP software could be made and used for reactor Neutronics calculation. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗
DOE OSTI2016

High resolution spectroscopy of jet cooled phenyl radical: The ν{sub 1} and ν{sub 2} a{sub 1} symmetry C–H stretching modes

A series of CH stretch modes in phenyl radical (C{sub 6}H{sub 5}) has been investigated via high resolution infrared spectroscopy at sub-Doppler resolution (∼60 MHz) in a supersonic discharge slit jet expansion. Two fundamental vibrations of a{sub 1} symmetry, ν{sub 1} and ν{sub 2}, are observed and rotationally analyzed for the first time, corresponding to in-phase and out-of-phase symmetric CH stretch excitation at the ortho/meta/para and ortho/para C atoms with respect to the radical center. The ν{sub 1} and ν{sub 2} band origins are determined to be 3073.968 50(8) cm{sup −1} and 3062.264 80(7) cm{sup −1}, respectively, which both agree within 5 cm{sup −1} with theoretical anharmonic scaling predictions based on density functional B3LYP/6-311g++(3df,3dp) calculations. Integrated band strengths for each of the CH stretch bands are analyzed, with the relative intensities agreeing remarkably well with theoretical predictions. Frequency comparison with previous low resolution Ar-matrix spectroscopy [A. V. Friderichsen et al., J. Am. Chem. Soc. 123, 1977 (2001)] reveals a nearly uniform Δν ≈ + 10-12 cm{sup −1} blue shift between gas phase and Ar matrix values for ν{sub 1} and ν{sub 2}. This differs substantially from the much smaller red shift (Δν ≈ − 1 cm{sup −1}) reported for the ν{sub 19} mode, and suggests a simple physical model in terms of vibrational mode symmetry and crowding due to the matrix environment. Finally, the infrared phenyl spectra are well described by a simple asymmetric rigid rotor Hamiltonian and show no evidence for spectral congestion due to intramolecular vibrational coupling, which bodes well for high resolution studies of other ring radicals and polycyclic aromatic hydrocarbons. In summary, the combination of slit jet discharge methods with high resolution infrared lasers enables spectroscopic investigation of even highly reactive combustion and interstellar radical intermediates under gas phase, jet-cooled (T{sub rot} ≈ 11 K) conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Next Generation Electromagnetic Pump Analysis Tools (PLM DOC-0005-2188). Final Report

This report provides the broad historical review of EM Pump development and details of MATRIX development under this project. This report summarizes the efforts made to modernize the legacy performance models used in previous EM Pump designs and the improvements made to the analysis tools. This report provides information on Tasks 1, 3, and 4 of the entire project. The research for Task 4 builds upon Task 1: Update EM Pump Databank and Task 3: Modernize the Existing EM Pump Analysis Model, which are summarized within this report. Where research for Task 2: Insulation Materials Development and Evaluation identified parameters applicable to the analysis model with Task 4, the analysis code was updated, and analyses were made for additional materials. The important design variables for the manufacture and operation of an EM Pump that the model improvement can evaluate are: space constraints; voltage capability of insulation system; maximum flux density through iron; flow rate and outlet pressure; efficiency and manufacturability. The development of the next-generation EM Pump analysis tools during this two-year program provides information in three broad areas: Status of analysis model development; Improvements made to older simulations; and Comparison to experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Constraints on texture zero and cofactor zero models for neutrino mass

Imposing a texture or cofactor zero on the neutrino mass matrix reduces the number of independent parameters from nine to seven. Since five parameters have been measured, only two independent parameters would remain in such models. We find the allowed regions for single texture zero and single cofactor zero models. We also find strong similarities between single texture zero models with one mass hierarchy and single cofactor zero models with the opposite mass hierarchy. We show that this correspondence can be generalized to texture-zero and cofactor-zero models with the same homogeneous costraints on the elements and cofactors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Study of Performance Portability of Low-bit Fused Matrix-Vector Multiplication Kernels in SYCL

Understanding the causes of performance gaps between a portable programming model and a vendor-specific programming model is important for improving performance portability. This paper studies performance portability of low-bit fused general matrix-vector multiplication kernels in SYCL on vendors’ graphics processing units (GPUs). This work introduces the use case, explains the kernel implementations in detail, evaluates the performance of the CUDA, HIP, and SYCL kernels on datacenter, desktop, and laptop GPUs, and investigates the causes of performance gaps. The results show that loop unrolling, kernel dispatch overhead, and sum reduction contribute to the gaps.

Jin, Zheming [ORNL] (ORCID:000000027197780X)↗
DOE OSTI2017

Microstructure-Level Modeling of Stage 3 Fission Gas Release in UO 2 Fuel

A new phase-field model of fission gas bubble evolution was developed and applied to gain an improved understanding of the microstructure-level processes leading to fission gas release from nuclear fuel and to inform improvements to the fission gas release model used in BISON. The phase-field model simulates the evolution of multiple fuel grains and intergranular bubbles simultaneously, tracks the local concentration of vacancies and gas atoms, and includes the effect of hydrostatic gas pressure on the surrounding fuel matrix. This represents the most realistic microstructure-level model of fission gas bubbles developed to date. The model was used to simulate the growth of grain face and grain edge bubbles in a hexagonally periodic 3D grain structure, which contains multiple grain boundaries and triple junctions. The coverage of grain boundaries and triple junctions by fission gas bubbles was calculated as a function of time. Preliminary results suggest that grain edge bubbles do not release a significant amount of gas prior to gas release by grain face bubbles, which supports the assumptions of the fission gas release model currently used in BISON. However, the simulations should be improved by including a physically realistic number of grain edge bubbles in the initial conditions before definitive conclusions can be drawn.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Spin-orbit correlations in the nucleon in the large- N c limit

We study the twist-3 spin-orbit correlations of quarks described by the nucleon matrix elements of the parity-odd rank-2 tensor QCD operator (the parity-odd partner of the QCD energy-momentum tensor). Our treatment is based on the effective dynamics emerging from the spontaneous breaking of chiral symmetry and the mean-field picture of the nucleon in the large- N c limit. The twist-3 QCD operators are converted to effective operators, in which the QCD interactions are replaced by spin-flavor-dependent chiral interactions of the quarks with the pion field. We compute the nucleon matrix elements of the twist-3 effective operators and discuss the role of the chiral interactions in the spin-orbit correlations. We derive the first-quantized representation in the mean-field picture and develop a quantum-mechanical interpretation. The chiral interactions give rise to new spin-orbit couplings and qualitatively change the correlations compared to the quark model picture. We also derive the twist-3 matrix elements in the topological soliton picture where the quarks are integrated out (skyrmion). The methods used here can be extended to other QCD operators describing higher-twist nucleon structure and generalized parton distributions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Electron transport in nanoscale junctions with local anharmonic modes

We study electron transport in nanojunctions in which an electron on a quantum dot or a molecule is interacting with an N-state local impurity, a harmonic (“Holstein”) mode, or a two-state system (“spin”). These two models, the Anderson-Holstein model and the spin-fermion model, can be conveniently transformed by a shift transformation into a form suitable for a perturbative expansion in the tunneling matrix element. We explore the current-voltage characteristics of the two models in the limit of high temperature and weak electron-metal coupling using a kinetic rate equation formalism, considering both the case of an equilibrated impurity, and the unequilibrated case. Specifically, we show that the analog of the Franck-Condon blockade physics is missing in the spin-fermion model. We complement this study by considering the low-temperature quantum adiabatic limit of the dissipative spin-fermion model, with fast tunneling electrons and a slow impurity. While a mean-field analysis of the Anderson-Holstein model suggests that nonlinear functionalities, bistability and hysteresis may develop, such effects are missing in the spin-fermion model at the mean-field level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Influence of helium bubbles on the bulk equation of state of gold under static compression

Nuclear materials often evolve into two-phase systems comprising a bulk matrix with dispersed inert-gas bubbles. The presence of these bubbles can have consequences to the thermomechanical response of materials and is a key life-limiting factor in some nuclear fuel forms. Understanding the behavior of these two-phase, bubble-matrix systems is, thus, important to improved predictive models and frameworks for many nuclear materials applications. While temperature excursions of these two-phase systems have been characterized, fewer studies have focused on the evolution of inert-gas bubbles under pressure. Here, in this paper, we use x-ray tools to interrogate a He-implanted gold foil to determine the pressure-dependent evolution of the individual components (Au matrix + bubbles), and we compare that total pressure dependence to theoretical equation-of-state descriptions based on mixing rules.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Approximate Hessian for accelerating ab initio structure relaxation by force fitting

We present a method to approximate the Hessian matrix of the Born-Oppenheimer energy landscape by using a simple force field model whose parameters are fitted to on-the-flight ab initio results. The inversed Hessian matrix is used as the preconditioner of conjugate gradient algorithms to speed up the atomic structure relaxation, resulting in a speedup factor of 2 to 5 on systems of bulk, slab, sheets, and atomic clusters. Because the force field model employed is simple and general, the parameter fitting is straightforward; the method is applicable to a variety of complicated systems for minimum structure relaxation. In the metal cluster new structure search, the new method yields better structures than the one obtained before with conventional algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels

In this study, a predictive modeling tool is developed for damage analysis and design of hydrogen (H 2 ) storage composite pressure vessels. It integrates micromechanics of matrix cracking into a continuum damage mechanics (CDM) description for damage evolution, and three-dimensional (3D) finite element (FE) modeling of the vessel structural response. At the scale of the composite layer (mesoscale), the temperature-dependent stiffness reduction law in terms of the damage variable for transverse matrix cracking is computed using an Eshelby-Mori-Tanaka approach for the initial composite thermoelastic properties and a self-consistent model for the stiffness reduction as a function of the damage variable. While transverse matrix cracking obeying a damage evolution relation can progressively evolve from an initiation to a saturation state, fiber failure is predicted by a micromechanical fiber rupture criterion that accounts for the fiber strength and matrix stress. The implementation of this integrated multiscale modeling model into a 3D FE formulation enables damage analysis and design of H 2 storage composite pressure vessels. The developed tool is illustrated through 3D damage analyses of a cryogenically compressed H 2 storage vessel model subjected to thermomechanical loadings to investigate effects of the helical layer fiber orientation and loading scenario on damage development, vessel integrity and burst pressure.

08 HYDROGEN↗
DOE OSTI2024

Posterior Covariance Matrix Approximations

Here, the Davis equation of state (EOS) is commonly used to model thermodynamic relationships for high explosive (HE) reactants. Typically, the parameters in the EOS are calibrated, with uncertainty, using a Bayesian framework and Markov Chain Monte Carlo (MCMC) methods. However, MCMC methods are computationally expensive, especially for complex models with many parameters. This paper provides a comparison between MCMC and less computationally expensive Variational methods (Variational Bayesian and Hessian Variational Bayesian) for computing the posterior distribution and approximating the posterior covariance matrix based on heterogeneous experimental data. All three methods recover similar posterior distributions and posterior covariance matrices. This study demonstrates that for this EOS parameter calibration application, the assumptions made in the two Variational methods significantly reduce the computational cost but do not substantially change the results compared to MCMC.

97 MATHEMATICS AND COMPUTING↗