Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physical quantities”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2023

Simultaneous measurements of volume, pressure, optical images, and crystal structure with a dynamic diamond anvil cell: A real-time event monitoring system

The dynamic diamond anvil cell (dDAC) technique has attracted great interest because it possibly provides a bridge between static and dynamic compression studies with fast, repeatable, and controllable compression rates. The dDAC can be a particularly useful tool to study the pathways and kinetics of phase transitions under dynamic pressurization if simultaneous measurements of physical quantities are possible as a function of time. We report the development of a real-time event monitoring (RTEM) system with dDAC, which can simultaneously record the volume, pressure, optical image, and structure of materials during dynamic compression runs. In particular, the volume measurement using both Fabry–Pérot interferogram and optical images facilitates the construction of an equation of state (EoS) using the dDAC in a home-laboratory. We also developed an in-line ruby pressure measurement (IRPM) system to be deployed at a synchrotron x-ray facility. This system provides simultaneous measurements of pressure and x-ray diffraction in low and narrow pressure ranges. The EoSs of ice VI obtained from the RTEM and the x-ray diffraction data with the IRPM are consistent with each other. The complementarity of both RTEM and IRPM systems will provide a great opportunity to scrutinize the detailed kinetic pathways of phase transitions using dDAC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices

In this paper, we report a random sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for under- standing the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbor hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbor hopping and 2D lattices. In addition, constraints on real-space con gurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Theory of optical responses in clean multi-band superconductors

Electromagnetic responses in superconductors provide valuable information on the pairing symmetry as well as physical quantities such as the superfluid density. However, at the superconducting gap energy scale, optical excitations of the Bogoliugov quasiparticles are forbidden in conventional Bardeen-Cooper-Schrieffer superconductors when momentum is conserved. Accordingly, far-infrared optical responses have been understood in the framework of a dirty-limit theory by Mattis and Bardeen for over 60 years. Here we show, by investigating the selection rules imposed by particle-hole symmetry and unitary symmetries, that intrinsic momentum-conserving optical excitations can occur in clean multi-band superconductors when one of the following three conditions is satisfied: (i) inversion symmetry breaking, (ii) symmetry protection of the Bogoliubov Fermi surfaces, or (iii) simply finite spin-orbit coupling with unbroken time reversal and inversion symmetries. This result indicates that clean-limit optical responses are common beyond the straightforward case of broken inversion symmetry. We apply our theory to optical responses in FeSe, a clean multi-band superconductor with inversion symmetry and significant spin-orbit coupling. This result paves the way for studying clean-limit superconductors through optical measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Analysis of streaked images of x-ray self-emission in laser-driven spherical implosions

Imaging of x-ray self-emission provides a powerful in situ measurement of the spatial and temporal evolution of high-energy-density plasmas. However, interpretation of these measurements requires detailed understanding of the data-generating process. This work presents a case study in the interpretation of x-ray self-emission data for the specific application of streaked one-dimensional slit imaging of spherical laser-driven implosions. A comprehensive generative model of the streaked slit-imaging diagnostic is developed including detailed treatments of the radiation transfer, photometrics, and photostatistics associated with the measurement. The model is used to generate realistic synthetic streaked images and to analyze experimental streaked images to extract important physical quantities of interest. An example analysis of streaked images from implosion experiments on the OMEGA laser is presented, where the model developed in this work is used to constrain the trajectory and peak velocity of the implosion using Bayesian inference.

Bayesian inference↗
DOE OSTI2024

$\overline{TKE}$ Parameterization and $\bar{v}$ Uncertainty Analysis for CGMF

Previous work was performed on tuning CGMF parameters for 235 U, 238 U, and Plutonium isotopes. Now work is being done to tune minor uranium isotopes. However, uranium isotopes like 232 U and 236 U have almost no experimental data. We are applying cross-isotope models to extrapolate and tune CGMF on isotopes that lack experimental data. There exist several internal CGMF physics quantities that affect the output of CGMF—multi-chance fission probability, excitation energy sharing, spin-cutoff factor, spin scaling, and fragment total kinetic energy to name a few. The mean fragment total kinetic energy, $\overline{TKE}$, is particularly interesting because of its strong anti-correlation with $\bar{v}$. We are most interested in the mean fragment total kinetic energy before neutron emissions. $\overline{TKE}$ is assumed to be pre-neutron emission unless otherwise stated. Currently in CGMF, the $\overline{TKE}$ model for 233,234,235,238 U are tuned independently to reproduce ν for the associated isotopes. In this report, we will tune a cross-isotope $\overline{TKE}$ model to experimental $\overline{TKE}$ data for 232,233,234,235,236,238 U. Because of the unreliable and sparse nature of $\overline{TKE}$ experimental data, future work will use more reliable experimental $\bar{v}$ data to infer the $\overline{TKE}$ model (and likely other internal CGMF parameters) for uranium isotopes. Such work has been performed previously using a sensitivity analysis and Kalman filter methods.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Nonresonant central exclusive production of charged-hadron pairs in proton-proton collisions at s = 13 TeV

The central exclusive production of charged-hadron pairs in p p collisions at a center-of-mass energy of 13 TeV is examined, based on data collected in a special high- β * run of the LHC. The nonresonant continuum processes are studied with the invariant mass of the centrally produced two-pion system in the resonance-free region, m π + π − < 0.7 or m π + π − > 1.8 GeV . Differential cross sections as functions of the azimuthal angle between the surviving protons, squared exchanged four-momenta, and m π + π − are measured in a wide region of scattered proton transverse momenta, between 0.2 and 0.8 GeV, and for pion rapidities | y | < 2 . A rich structure of interactions related to double-pomeron exchange is observed. A parabolic minimum in the distribution of the two-proton azimuthal angle is observed for the first time. It can be interpreted as an effect of additional pomeron exchanges between the protons from the interference between the bare and the rescattered amplitudes. After model tuning, various physical quantities are determined that are related to the pomeron cross section, proton-pomeron and meson-pomeron form factors, pomeron trajectory and intercept, and coefficients of diffractive eigenstates of the proton. © 2024 CERN, for the CMS and TOTEMs Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Capturing soil moisture and salinity changes in flooded coastal forests using electrical resistivity and induced polarization

Geophysical methods provide high-resolution spatial measurements of physical quantities sensitive to changes in soil moisture and salinity and can be used to monitor soil hydrological responses to flooding. However, extracting quantitative hydrological information from geophysical data remains challenging. In this study, we extended existing petrophysical models to estimate soil moisture and salinity from electrical measurements to address this challenge. We manipulated two hydrologically isolated 2000 m 2 experimental plots by simultaneously inundating them with 265 m 3 of either freshwater or estuarine water. Repeated electrical resistivity and induced polarization measurements were used to image the water and solute infiltration along two transects that are 100 and 42 m in length. Petrophysical models derived from laboratory multi-salinity electrical measurements were used to estimate changes in soil moisture and fluid salinity from field measurements of real and imaginary conductivity during the ecosystem-scale flooding experiment. The real conductivity increased by ∼100% in the freshwater plot and ∼570% in the saltwater plot. The change in imaginary conductivity in the freshwater plot was <1 mS/m, whereas that of the estuarine water plot was ∼5 mS/m. Real conductivity shows dependence on soil moisture content with a coefficient of determination (R 2 ) > 0.7, while the imaginary conductivity shows a dependence on soil salinity with R 2 > 0.6. The results validate the use of electrical resistivity for estimating changes in soil moisture content in response to flooding. Combining electrical resistivity imaging with induced polarization measurements provides the possibility to account for changes in pore fluid conductivity.

Adebayo, Moses B. [Univ. of Toledo, OH (United Sta↗
DOE OSTI2026

A sweeping positivity-preserving high-order finite difference WENO scheme for Euler equations

We develop a simple, high-order, conservative and robust positivity-preserving sweeping procedure for the density and the nonlinear pressure function in the compressible Euler equations. Using the scaling limiter in Zhang and Shu (J Comput Phys 229:3091–3120, 2010), we obtain a non-trivial extension of the scalar sweeping technique in Liu et al. (J Sci Comput 73:1028–1071, 2017) for the positivity of pressure. The sweeping procedure developed in this paper is a post-processing technique, which can be applied to any concave functions of the conserved variables in hyperbolic conservation law systems. Thus, it has applications beyond the Euler equations. This procedure preserves positivity and conservation of physical quantities without destroying the accuracy of the underlying scheme. The algorithm works for general schemes including finite difference, finite volume and discontinuous Galerkin methods; however, in this paper we focus on finite difference weighted essentially non-oscillatory (WENO) methods. As a result, we provide numerical tests of the fifth-order finite difference WENO scheme to demonstrate the accuracy and robustness of the technique.

Compressible Euler equations↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

SO(3)-invariance of informed-graph-based deep neural network for anisotropic elastoplastic materials

This work examines the frame-invariance (and the lack thereof) exhibited in simulated anisotropic elasto-plastic responses generated from supervised machine learning of classical multi-layer and informed-graph-based neural networks, and proposes different remedies to fix this drawback. The inherent hierarchical relations among physical quantities and state variables in an elasto-plasticity model are first represented as informed, directed graphs, where three variations of the graph are tested. While feed-forward neural networks are used to train path-independent constitutive relations (e.g., elasticity), recurrent neural networks are used to replicate responses that depends on the deformation history, i.e. or path dependent. In dealing with the objectivity deficiency, we use the spectral form to represent tensors and, subsequently, three metrics, the Euclidean distance between the Euler Angles, the distance from the identity matrix, and geodesic on the unit sphere in Lie algebra, can be employed to constitute objective functions for the supervised machine learning. In this, the aim is to minimize the measured distance between the true and the predicted 3D rotation entities. Following this, we conduct numerical experiments on how these metrics, which are theoretically equivalent, may lead to differences in the efficiency of the supervised machine learning as well as the accuracy and robustness of the resultant models. Neural network models trained with tensors represented in component form for a given Cartesian coordinate system are used as a benchmark. Our numerical tests show that, even given the same amount of information and data, the quality of the anisotropic elasto-plasticity model is highly sensitive to the way tensors are represented and measured. The results reveal that using a loss function based on geodesic on the unit sphere in Lie algebra together with an informed, directed graph yield significantly more accurate rotation prediction than the other tested approaches.

42 ENGINEERING↗
DOE OSTI2021

Errors introduced in fission neutron spectrum measurements using a single reference

Measurements of physical quantities are often made relative to a reference quantity, simplifying the interpretation of the measurement and its uncertainties. However, if applied without careful considerations, this technique can carry with it a significant and complex systematic error that is usually ignored. For example, measurements of continuous neutron spectra relative to a reference spectrum yield a result that does not properly include the important effects of neutron scattering in the experimental environment. These effects act to distort the measured spectrum shape and create notable errors in spectrum-integrated quantities. While this effect is one of a series of known potential sources of systematic error including, but not limited to target impurities, angular distributions, incident beam flux, and others, the error introduced from differences in environmental scattering effects between the reference and desired spectra are frequently overlooked. In this work, we demonstrate the origin of this error using simulated data corresponding to measurements of the neutron-induced prompt fission neutron spectra of 233 U, 235 U, 238 U, and 239 Pu as well as the spontaneous fission neutron spectrum of 252 Cf. Errors in the spectrum shape can be as high as 10%–15% depending on the combination of incident and outgoing neutron energy. The spectrum integral and average spectrum energy also contain errors that change with incident neutron energy and can be as high as ~2%. Although this study is based on simulated data, we also show a method for proving that this error exists in an experimental data set and we show a correction that can be applied to results of this type to obtain a more accurate result, thereby largely removing systematic errors from environmental scattering effects. The potential for these kinds of errors should be considered in measurements producing continuous neutron spectra.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

The upgraded summing NaI(Tl) (SuN++) absorption spectrometer

Simulations of astrophysical processes require a plethora of nuclear physics input. In particular, models of neutron-capture nucleosynthesis like the s, i, and r processes require β-decay information and experimentally constrained neutron-capture reaction rates. Past experiments with the 4π Summing NaI(Tl) (SuN) total absorption spectrometer have provided these physics quantities. Here, we outline an upgrade of SuN to SuN++, where 20 new segments (12 NaI(Tl) and 8 CeBr 3 ) have been integrated into the pre-existing SuN total absorption spectrometer to provide increased energy and time resolution in β-decay experiments. The details of the newly upgraded SuN++ total absorption spectrometer are discussed with results from the commissioning experiment at the Facility for Rare Isotope Beams (FRIB) utilizing a 70 Cu beam.

CeBr3↗
DOE OSTI2025

A Machine Learning Framework for Modeling Ensemble Properties of Atomically Disordered Materials

Atomic disorder can strongly influence material properties such as charge transport, optical response, and catalytic activity. However, efficiently modeling these disorder effects remains challenging for first-principles methods due to the cost of sampling large configurational spaces and computing complex physical quantities. Recent advances of machine learning techniques, particularly graph neural networks (GNNs), has enabled the efficient and accurate predictions of complex material properties, offering promising tools for studying disordered systems. In this work, we present a general machine-learning-assisted computational framework that integrates equivariant GNNs with Monte Carlo simulations to compute the thermodynamic and ensemble-averaged functional properties of disordered materials. Using the surface-termination-disordered MXene monolayer Ti 3 C 2 T 2–x as a representative system, we find that electrical conductivity exhibits an emergent peak near the order–disorder phase transition temperature due to the interplay between electron scattering and doping. In contrast, optical conductivity remains largely insensitive to local atomic disorder and reflects the global surface chemical composition. These results highlight the role of atomic disorder in affecting material properties and demonstrate the potential of our approach for statistically modeling disorder effects in a wide range of materials such as high-entropy alloys and spin liquids.

MXene↗
DOE OSTI2026

Models and Algorithms for Equilibrium Analysis of Mixed-Material Nucleic Acid Systems

Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.

2′OMe-RNA↗
DOE OSTI2023

Enhanced active-site electric field accelerates enzyme catalysis

In this study, the design and improvement of enzymes based on physical principles remain challenging. Here we demonstrate that the principle of electrostatic catalysis can be leveraged to substantially improve a natural enzyme’s activity. We enhanced the active-site electric field in horse liver alcohol dehydrogenase by replacing the serine hydrogen-bond donor with threonine and replacing the catalytic Zn 2+ with Co 2+ . Based on the electric field enhancement, we make a quantitative prediction of rate acceleration—50-fold faster than the wild-type enzyme—which was in close agreement with experimental measurements. The effects of the hydrogen bonding and metal coordination, two distinct chemical forces, are described by a unified physical quantity—electric field, which is quantitative, and shown here to be additive and predictive. These results suggest a new design paradigm for both biological and non-biological catalysts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Cost function for low-dimensional manifold topology assessment

Abstract In reduced-order modeling, complex systems that exhibit high state-space dimensionality are described and evolved using a small number of parameters. These parameters can be obtained in a data-driven way, where a high-dimensional dataset is projected onto a lower-dimensional basis. A complex system is then restricted to states on a low-dimensional manifold where it can be efficiently modeled. While this approach brings computational benefits, obtaining a good quality of the manifold topology becomes a crucial aspect when models, such as nonlinear regression, are built on top of the manifold. Here, we present a quantitative metric for characterizing manifold topologies. Our metric pays attention to non-uniqueness and spatial gradients in physical quantities of interest, and can be applied to manifolds of arbitrary dimensionality. Using the metric as a cost function in optimization algorithms, we show that optimized low-dimensional projections can be found. We delineate a few applications of the cost function to datasets representing argon plasma, reacting flows and atmospheric pollutant dispersion. We demonstrate how the cost function can assess various dimensionality reduction and manifold learning techniques as well as data preprocessing strategies in their capacity to yield quality low-dimensional projections. We show that improved manifold topologies can facilitate building nonlinear regression models.

42 ENGINEERING↗
DOE OSTI2024

Microscopic Theory of Nonlinear Hall Effect in Three-Dimensional Magnetic Systems

Abstract Nonlinear Hall effect (NLHE) has been detected in various of condensed matter systems. Unlike linear Hall effect, NLHE may exist in physical systems with broken inversion symmetry in crystals. On the other hand, real space spin texture may also break inversion symmetry and result in NLHE. We employ the Feynman diagrammatic technique to calculate non-linear Hall conductivity (NLHC) in three-dimensional magnetic systems. The results connect NLHE with the physical quantity of emergent electrodynamics which originates from the magnetic texture. The leading order contribution of NLHC, χabb , is proportional to the emergent toroidal moment T a e , which reflects how the spin textures wind in three dimensions.

Hou 侯, Wen-Tao 文涛↗