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

The interpretation of temperature and salinity variables in numerical ocean model output and the calculation of heat fluxes and heat content

Abstract. The international Thermodynamic Equation of Seawater 2010 (TEOS-10) defined the enthalpy and entropy of seawater, thus enabling the global ocean heat content to be calculated as the volume integral of the product of in situ density, ρ, and potential enthalpy, h0 (with reference sea pressure of 0 dbar). In terms of Conservative Temperature, Θ, ocean heat content is the volume integral of ρcp0Θ, where cp0 is a constant “isobaric heat capacity”. However, many ocean models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) as well as all models that contributed to earlier phases, such as CMIP5, CMIP3, CMIP2, and CMIP1, used EOS-80 (Equation of State – 1980) rather than the updated TEOS-10, so the question arises of how the salinity and temperature variables in these models should be physically interpreted, with a particular focus on comparison to TEOS-10-compliant observations. In this article we address how heat content, surface heat fluxes, and the meridional heat transport are best calculated using output from these models and how these quantities should be compared with those calculated from corresponding observations. We conclude that even though a model uses the EOS-80, which expects potential temperature as its input temperature, the most appropriate interpretation of the model's temperature variable is actually Conservative Temperature. This perhaps unexpected interpretation is needed to ensure that the air–sea heat flux that leaves and arrives in atmosphere and sea ice models is the same as that which arrives in and leaves the ocean model. We also show that the salinity variable carried by present TEOS-10-based models is Preformed Salinity, while the salinity variable of EOS-80-based models is also proportional to Preformed Salinity. These interpretations of the salinity and temperature variables in ocean models are an update on the comprehensive Griffies et al. (2016) paper that discusses the interpretation of many aspects of coupled Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improved coarse-graining methods for two dimensional tensor networks including fermions

We show how to apply renormalization group algorithms incorporating entanglement filtering methods and a loop optimization to a tensor network which includes Grassmann variables which represent fermions in an underlying lattice field theory. As a numerical test a variety of quantities are calculated for two dimensional Wilson-Majorana fermions and for the two flavor Gross-Neveu model. The improved algorithms show much better accuracy for quantities such as the free energy and the determination of Fisher’s zeros.

Phase transitions↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2024

New proxies for second-order cumulants of conserved charges in heavy-ion collisions within the EPOS4 framework

Proxies for cumulants of baryon number 𝐵, electric charge 𝑄, and strangeness 𝑆 are usually measured in heavy-ion collisions via moments of net-number distribution of given hadronic species. Since these cumulants of conserved charges are expected to be sensitive to the existence of a critical point in the phase diagram of nuclear matter, it is crucial to ensure that the proxies used as substitutes are as close to them as possible. Hence, we use the EPOS 4 framework to generate Au + Au collisions at several collision energies of the BNL Relativistic Heavy Ion Collider beam energy scan. We compute second-order net cumulants of 𝜋, 𝐾, and 𝑝, for which experimental data have been published as well as the corresponding conserved charge cumulants. We then compare them with proxies, defined in previous lattice QCD and hadron resonance gas model studies, which are shown to reproduce more accurately their associated conserved charge cumulants. We investigate the impact of hadronic rescatterings occurring in the late evolution of the system on these quantities, as well as the amount of signal actually originating from the bulk medium which endures a phase transition.

Physics↗
DOE OSTI2020

Cooperative Research and Development Agreement between National Energy Technology Laboratory and Korea Institute of Energy Research (KIER) [Abstract]

KIER and NETL will collaborate to study hydrate-based technologies to enhance desalination, gas separation, and hydrocarbon flow through pipelines and wellbores. With optimized selections of hydrate formers and reactors, gas hydrates can be effective and economic mechanisms of separation for solid-water and gas mixtures, and gas hydrate formation or dissociation can be significantly delayed or expedited with suitably selected additives. Comprehensive laboratory studies, including identifying novel hydrate formers or inhibitors to control temperature and pressure requirements and understanding their working mechanism or developing reactor systems for production of large quantities of hydrates, will be performed at both KIER and NETL either parallelly or serially.

03 NATURAL GAS↗
DOE OSTI2024

Final Technical Report: Center for Mechanistic Control of Unconventional Formations

The overarching mission of the Center for Mechanistic Control of Unconventional Formations (CMC-UF) was to garner cross-cutting, fundamental, geoscience knowledge to achieve mechanistic control over the strongly coupled nonequilibrium physical and geochemical processes in extreme geological environments including shale, mudstone, marls, and other tight rocks with nanoscale pores. Collectively, these are referred to as unconventional formations and they often play the role of seals for other subsurface storage formations. The fundamental knowledge garnered by CMC-UF enabled science-based management of US unconventional formations for subsurface storage of carbon dioxide and TWh quantities of renewable energy as hydrogen and/or compressed air over longer timeframes as well as for natural gas production, with reduced environmental impacts, in the short term.

58 GEOSCIENCES↗
DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗
DOE OSTI2020

Experimental validation of scintillator detector response to correlated neutrons with MCNP and associated packages

Time, energy, and angle correlated neutrons are emitted in nuclear fission. New, physics-based fission event generators that capture these correlations have been incorporated into MCNP6.2 because correlations can impact nuclear non-proliferation and safeguards systems response. DRiFT is a detector response code that takes MCNP outputs and reproduces laboratory analogous detector quantities. Previous work has compared DRiFT to single detector experiment data and has discussed recent updates to DRiFT specific to correlated fission experiments. In this work, we demonstrate a validation of recent updates to DRiFT for organic scintillator response in correlated fission neutron experiments using a 252Cf source. Furthermore, outputs from MCNP with the FREYA fission event generator were post-processed with DRiFT and compared to experiment data. Simulation and experiment results are in reasonably good agreement for multiple correlated neutron quantities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Dilute Combustion Control Using Spiking Neural Networks

Dilute combustion with exhaust gas recirculation (EGR) in spark-ignition engines presents a cost-effective method for achieving higher levels of engine efficiency. At high levels of EGR, however, cycle-to-cycle variability (CCV) of the combustion process is exacerbated by sporadic occurrences of misfires and partial burns. Previous studies have shown that temporal deterministic patterns emerge at such conditions and certain combustion cycles have a significant influence over future events. Due to the complexity of the combustion process and the nature of CCV, harnessing all the deterministic information for control purposes has remained challenging even with physics based 0-D, 1-D, and high-fidelity computational fluid dynamics (CFD) models. In this study, we present a data-driven approach to optimize the combustion process by controlling CCV adjusting the cycle-to-cycle fuel injection quantity. Readily available data from in-cylinder pressure was used to train a spiking neural network (SNN) which learns the optimal way to manage fuel injection in order to reduce CCV while maintaining acceptable levels of fuel consumption. SNNs are particularly well suited for powertrain control applications due to their ability to be deployed on FPGA-based neuromorphic hardware which are small, inexpensive, and have a low power demand. The high-performance computing (HPC) resources of Oak Ridge National Laboratory were used to run an evolutionary-based training approach for choosing the best SNN configuration that minimizes the size of the network while achieving the desired goal. The neuromorphic hardware with the optimized SNN deployed was connected to the rapid prototyping engine control system for real-time control implementation and tested on a single cylinder version of a GM LNF 4-cylinder engine. The results show a significant reduction of CCV with a small percentage of additional fuel used to stabilize the charge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

FLAG review 2024

We review lattice results related to pion, kaon, 𝐷-meson, 𝐵-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor 𝑓+⁡(0) arising in the semileptonic 𝐾 → 𝜋 transition at zero momentum transfer, as well as the decay-constant ratio 𝑓 𝐾 ⁡/𝑓 𝜋 and its consequences for the Cabibbo–Kobayashi–Maskawa (CKM) matrix elements 𝑉 𝑢⁢𝑠 and 𝑉 𝑢⁢𝑑 . We review the determination of the 𝐵 𝐾 parameter of neutral kaon mixing as well as the additional four 𝐵 parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for 𝑚 𝑐 and 𝑚 𝑏 as well as those for the decay constants, form factors, and mixing parameters of charmed and bottom mesons and baryons. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant 𝛼 𝑠 . We review the determinations of nucleon charges from the matrix elements of both isovector and flavor-diagonal axial, scalar and tensor local quark bilinears, and momentum fraction, helicity moment and the transversity moment from one-link quark bilinears. We also review determinations of scale-setting quantities. Finally, in this review we have added a new section on the general definition of the low-energy limit of the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Radiation Transport Simulation of the Sulfur Counting System in the Radiation Metrology Laboratory (RML)

At the request of staff members from the Radiation Metrology Laboratory (RML), a series of Monte Carlo radiation transport calculations were performed using two different models of the detector geometry of the RMLs sulfur counting system. The fraction of electrons from each β-decay of 32 P in the sulfur pellet that enter the window of the sulfur counting system was calculated with both MCNP and ITS. In addition, the differential energy distributions of the electrons entering the counting system window were computed. There was significant agreement between the integral and differential quantities calculated by the two transport codes. Summary tables are for the integral efficiency values are found in the body of the report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Isotope Production at the Jefferson Lab LERF Facility (Final Report)

The project initially proceeded well, meeting all of the Milestones. However, the need to create a new irradiation area in the LERF (the initial proposal was to use the CEBAF Injector) used up a larger fraction of the budget than foreseen. A further problem arose when the hot-cell became unavailable at VCU, followed by VCU dropping out of the collaboration. This meant that we were unable to separate the 67 Cu from the gallium. We replaced this by carrying out detailed gamma spectroscopy to determine the quantity of 67 Cu produced.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Error-Bounded Learned Scientific Data Compression with Preservation of Derived Quantities

Scientific applications continue to grow and produce extremely large amounts of data, which require efficient compression algorithms for long-term storage. Compression errors in scientific applications can have a deleterious impact on downstream processing. Thus, it is crucial to preserve all the “known” Quantities of Interest (QoI) during compression. To address this issue, most existing approaches guarantee the reconstruction error of the original data or primary data (PD), but cannot directly control the problem of preserving the QoI. In this work, we propose a physics-informed compression technique that is composed of two parts: (i) reduction of the PD with bounded errors and (ii) preservation of the QoI. In the first step, we combine tensor decompositions, autoencoders, product quantizers, and error-bounded lossy compressors to bound the reconstruction error at high levels of compression. In the second step, we use constraint satisfaction post-processing followed by quantization to preserve the QoI. To illustrate the challenges of reducing the reconstruction errors of the PD and QoI, we focus on simulation data generated by a large-scale fusion code, XGC, which can produce tens of petabytes in a single day. The results show that our approach can achieve a high compression amount while accurately preserving the QoI within scientifically acceptable bounds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Towards advanced polarized electron sources

Polarized electrons play an important role in high-energy and nuclear physics, and their properties have also been exploited in ultrafast electron microscopy. Currently, gallium arsenide crystals illuminated by circular polarized infrared laser light are commonly used for generating polarized electrons. However, the achievable accelerating voltage and the gradient of these electrostatic sources limit the beam quality and quantity. A solution could be to combine gallium arsenide photocathodes with radio-frequency electron guns, which are capable of accelerating beams with significantly higher gradients and voltage. Here we report the successful operation of a gallium arsenide photocathode in a superconducting radio-frequency gun. Our findings are relevant for future sources of polarized electrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Trailhead for quantum simulation of SU(3) Yang-Mills lattice gauge theory in the local multiplet basis

Maintaining local interactions in the quantum simulation of gauge field theories relegates most states in the Hilbert space to be unphysical—theoretically benign, but experimentally difficult to avoid. Reformulations of the gauge fields can modify the ratio of physical to gauge-variant states often through classically preprocessing the Hilbert space and modifying the representation of the field on qubit degrees of freedom. This paper considers the implications of representing SU(3) Yang-Mills gauge theory on a lattice of irreducible representations in both a global basis of projected global quantum numbers and a local basis in which controlled-plaquette operators support efficient time evolution. Classically integrating over the internal gauge space at each vertex (e.g., color isospin and color hypercharge) significantly reduces both the qubit requirements and the dimensionality of the unphysical Hilbert space. Initiating tuning procedures that may inform future calculations at scale, the time evolution of one and two plaquettes are implemented on one of IBM’s superconducting quantum devices, and early benchmark quantities are identified. The potential advantages of qudit environments, with either constrained two-dimensional hexagonal or one dimensional nearest-neighbor internal state connectivity, are discussed for future large-scale calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗