Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Marcus equation”

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.

54 records · Page 3

Cluster bootstrap for cosmological correlators

We show that cosmological wavefunction coefficients associated with n-site chain and loop graphs for a cubic scalar theory in de Sitter spacetime have symbol alphabets given by subsets of A 2n−2 and B 2n−1 cluster variables, respectively, and satisfy the associated cluster adjacency properties. The key step in proving this is identifying a precise connection between graph “tubings” that appear in the kinematic flow equation and polygon “triangulations” that encode the combinatorics of cluster compatibility. Our results imply that cosmological wavefunction coefficients in a general power-law FRW cosmology satisfy cluster adjacency to all orders in the ϵ expansion around the de Sitter limit. We use this information as bootstrap input to show that de Sitter symbols for n ≤ 4 are uniquely determined by simple physical constraints.

differential and algebraic geometry↗

Near-global summer circulation response to the spring surface temperature anomaly in Tibetan Plateau –– the GEWEX/LS4P first phase experiment

Subseasonal to seasonal (S2S) prediction of droughts and floods is one of the major challenges of weather and climate prediction. Recent studies suggest that the springtime land surface temperature/subsurface temperature (LST/SUBT) over the Tibetan Plateau (TP) can be a new source of S2S predictability. The project “Impact of Initialized Land Surface Temperature and Snowpack on Subseasonal to Seasonal Prediction (LS4P)” was initiated to study the impact of springtime LST/SUBT anomalies over high mountain areas on summertime precipitation predictions. The present work explores the simulated global scale response of the atmospheric circulation to the springtime TP land surface cooling by 16 current state-of-the-art Earth System Models (ESMs) participating in the LS4P Phase I (LS4P-I) experiment. The LS4P-I results show, for the first time, that springtime TP surface anomalies can modulate a persistent quasi-barotropic Tibetan Plateau-Rocky Mountain Circumglobal (TRC) wave train from the TP via the northeast Asia and Bering Strait to the western part of the North America, along with the springtime westerly jet from TP across the whole North Pacific basin. The TRC wave train modulated by the TP thermal anomaly play a critical role on the early summer surface air temperature and precipitation anomalies in the regions along the wave train, especially over the northwest North America and the southern Great Plains. The participant models that fail in capturing the TRC wave train greatly under-predict climate anomalies in reference to observations and the successful models. These results suggest that the TP LST/SUBT anomaly via the TRC wave train is the first order source of the S2S variability in the regions mentioned. Furthermore, the TP surface temperature anomaly can influence the Southern Hemispheric circulation by generating cross-equator wave trains. However, the simulated propagation pathways from the TP into the Southern Hemisphere show large inter-model differences. More dynamical understanding of the TRC wave train as well as its cross-equator propagation into the Southern Hemisphere will be explored in the newly launched LS4P phase II experiment.

54 ENVIRONMENTAL SCIENCES↗

Enhanced dilepton emission from a phase transition in dense matter

It is demonstrated that the presence of a phase transition in heavy-ion collisions, at beam energies that probe dense quantum chromodynamics (QCD) matter, leads to a significant enhancement of the dilepton yield at low invariant mass and requires the embedding of low beam energies per produced pion due to the extended emission time. In addition, the temperature of low-mass dileptons shows a modest decrease due to the mixed phase. The emission of dileptons in the SIS18–SIS100 beam energy range is studied by augmenting the ultra-relativistic quantum molecular dynamics (UrQMD) transport model with a realistic density-dependent equation of state, as well as two different phase transitions. This is achieved by extending the molecular dynamics interaction part of the UrQMD model to a density-dependent interaction potential with a high-density minimum, leading to a phase transition and metastable coexisting high-density states. Together with a high-precision measurement, these simulations will be able to constrain the existence of a phase transition in QCD up to densities of several times nuclear saturation density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Systematics of the chemical freeze-out line in the high baryon density regime explored at SIS100

The systematic uncertainties of chemical freeze-out fits at SIS100 energies (Au+Au reactions at $\sqrt{s_{NN}}$ = 3 – 5 GeV) are studied using UrQMD simulations. Although hadron production in UrQMD does not occur on a sharp chemical freeze-out hyper-surface, the extracted fit quality is shown to be very good. The extracted chemical parameters depend on the selected hadron species as well as the underlying equation of state (EoS) of the matter. Including light nuclei and anti-protons in the fit increases the expected freeze-out temperature, while a stiffer EoS increases the obtained chemical potential. Similarly, the baryon densities extracted by the thermal fits depend on the choice of hadrons as well as the underlying equation of state. These results are important for the upcoming CBM@FAIR physics program and highlight that a degree of caution is advised when one relates the chemical freeze-out curve to features on the QCD phase diagram like the critical endpoint or a possible phase transition.

Hofmann, Emma Lilith [GSI Helmholtzzentrum für Sch↗

Low-lying dipole response of Ni 64

In this work, two complementary real-photon scattering experiments were conducted on the proton-magic 64 Ni nucleus to study the dipole response up to its neutron-separation energy of $S_n$ = 9.7 MeV. By combining both measurements, 87 $E$1 and 23 $M$1 transitions were identified above 4.3 MeV. The results of the observed $M$1 transitions were compared to shell-model calculations using two different model spaces. It was found that the inclusion of excitations across the $Z$ = 28 shell gap in the calculations has a large impact. Furthermore, average cross sections for decays to the ground state (elastic transitions) as well as to lower-lying excited states (inelastic decays) were determined. The corresponding $E$1 channel was compared to calculations within the relativistic equation of motion (REOM) framework. Whereas the calculations of highest possible complexity reproduce the fragmentation and overall behavior of the $E$1 average elastic cross section well, the predicted absolute cross sections are approximately twice as high as the experimental upper limits even though the latter also include an estimate of the inelastic-decay channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dense nuclear matter equation of state from heavy-ion collisions

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. Furthermore, this white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Extreme compression of planetary gases: High-accuracy pressure-density measurements of hydrogen-helium mixtures above fourfold compression

Hydrogen (H 2 ) and helium (He), the most abundant elements in the universe, pose a unique challenge in measuring the equation of state of the mixture, owing to their differing physical properties. There remains a need for data with high enough precision to discriminate between existing equation of state (EOS) mix models in order to understand the internal structure of gas-giant planets. Here, we have measured the EOS of precompressed H 2 - He mixtures at conditions directly relevant to the planetary interiors using hypervelocity gas guns and Sandia’s Z machine with less than 10% uncertainty in density, enabling validation of mixture models. We precompressed 50:50 molar mixtures of H 2 -He to 0.1–0.2 GPa and directly measured particle velocity (in gas-gun experiments) and shock velocities (in Z-machine experiments). To complement the experimental efforts, we also computed the Hugoniots of precompressed H 2 -He mixtures using density-functional-theory-based molecular dynamics. Furthermore, we observe approximately 3- to 4.3-fold density compression at pressures up to 44 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

CTF Validation and Verification: Version 4.3

Coolant-Boiling in Rod Arrays- Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) analysis. It uses a two-fluid, three-field (i.e., fluid film, fluid drops, and vapor) modeling approach. Both subchannel and 3D Cartesian forms of nine conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and has been used and modified by several institutions over the last several decades. COBRA-TF is also used at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently became the CTF. One part of the improvement process includes validating the methods in CTF. This document seeks to provide a certain level of certainty and confidence in the predictive capabilities of the code for the scenarios it was designed to model—rod bundle geometries with operating conditions that are representative of prototypical pressurized water reactor (PWR)s and boiling water reactor (BWR)s in both normal and accident conditions. This is done by modeling a variety of experiments that simulate these scenarios and then presenting a qualitative and quantitative analysis of the results that demonstrates the accuracy to which CTF is capable of capturing specific quantities of interest.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CTF Validation and Verification: Version 4.4

Coolant-Boiling in Rod Arrays- Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) analysis. It uses a two-fluid, three-field (i.e., fluid film, fluid drops, and vapor) modeling approach. Both subchannel and 3D Cartesian forms of nine conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and has been used and modified by several institutions over the last several decades. COBRA-TF is also used at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently became the CTF. One part of the improvement process includes validating the methods in CTF. This document seeks to provide a certain level of certainty and confidence in the predictive capabilities of the code for the scenarios it was designed to model—rod bundle geometries with operating conditions that are representative of prototypical pressurized water reactor (PWR)s and boiling water reactor (BWR)s in both normal and accident conditions. This is done by modeling a variety of experiments that simulate these scenarios and then presenting a qualitative and quantitative analysis of the results that demonstrates the accuracy to which CTF is capable of capturing specific quantities of interest.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data

Three-dimensional, time-resolved blood flow measurement (4D-flow) is a powerful research and clinical tool, but improved resolution and scan times are needed. Therefore, this study aims to (1) present a postprocessing framework for optimization-driven simulation-based flow imaging, called 4D-flow High-resolution Imaging with a priori Knowledge Incorporating the Navier-Stokes equations and the discontinuous Galerkin method (4D-flow HIKING), (2) investigate the framework in synthetic tests, (3) perform phantom validation using laser particle imaging velocimetry, and (4) demonstrate the use of the framework in vivo. An optimizing computational fluid dynamics solver including adjoint-based optimization was developed to fit computational fluid dynamics solutions to 4D-flow data. Synthetic tests were performed in 2D, and phantom validation was performed with pulsatile flow. Reference velocity data were acquired using particle imaging velocimetry, and 4D-flow data were acquired at 1.5 T. In vivo testing was performed on intracranial arteries in a healthy volunteer at 7 T, with 2D flow as the reference. Results Synthetic tests showed low error (0.4%-0.7%). Phantom validation showed improved agreement with laser particle imaging velocimetry compared with input 4D-flow in the horizontal (mean -0.05 vs -1.11 cm/s, P < .001; SD 1.86 vs 4.26 cm/s, P < .001) and vertical directions (mean 0.05 vs -0.04 cm/s, P = .29; SD 1.36 vs 3.95 cm/s, P < .001). In vivo data show a reduction in flow rate error from 14% to 3.5%. Phantom and in vivo results from 4D-flow HIKING show promise for future applications with higher resolution, shorter scan times, and accurate quantification of physiological parameters.

4D-flow MRI↗

A Platform-Independent X-ray Diffraction Diagnostic for Phase Transition Kinetics in Traditional and Synthetic Microstructure Materials (LDRD Project 213088 Final Report)

Pulsed-power generators using the magnetic loading technique are able to produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength-strain rate dependence, structural phase transitions, and density of crystal defects such as dislocations. Here, we present a cost effective, compact X-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically-driven ramp compression of materials with single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials community. The success in fielding this new XRD diagnostic dramatically improves our predictive capability and understanding of rate-dependent behavior at or near phase transition. As Sandia plans the next-generation pulse-power driver platform, a key element needed to deliver new state-of-the-art experiments will be having the necessary diagnostic tools to probe new regimes and phenomena. These diagnostics need to be as versatile, compact, and portable as they are powerful. The development of a platform-independent XRD diagnostic gives Sandia researchers a new window to study the microstructure and phase dynamics of materials under load. This project has paved the way for phase transition research in a variety of materials with mission interest.

36 MATERIALS SCIENCE↗

Nitrogen fixation rates in the Guinea Dome and the equatorial upwelling regions in the Atlantic Ocean

Abstract Biological nitrogen fixation is a key process balancing the loss of combined nitrogen in the marine nitrogen cycle. Its relevance in upwelling or high nutrient regions is still unclear, with the few available studies in these regions of the ocean reporting rates that vary widely from below detection limit to > 100 nmol N L −1 d −1 . In the eastern tropical Atlantic Ocean, two open ocean upwelling systems are active in boreal summer. One is the seasonal equatorial upwelling, where the residual phosphorus associated with aged upwelled waters is suggested to enhance nitrogen fixation in this season. The other is the Guinea Dome, a thermal upwelling dome. We conducted two surveys along 23° W across the Guinea Dome and the Equator from 15° N to 5° S in September 2015 and August–September 2016 with high latitudinal resolution (20–60 nm between stations). The abundance of Trichodesmium colonies was characterized by an Underwater Vision Profiler 5 and the total biological nitrogen fixation in the euphotic layer was measured using the 15 N 2 technique. The highest abundances of Trichodesmium colonies were found in the area of the Guinea Dome (9°–15° N) with a maximum of 3 colonies L −1 near the surface. By contrast, colonies were almost absent in the Equatorial band between 2° N and 5° S. The highest nitrogen fixation rate was measured at the northern edge of the Guinea Dome in 2016 (ca. 31 nmol N L −1 d −1 ). In this region, where diazotrophs thrived on a sufficient supply of both phosphorus and iron, a patchy distribution was unveiled by our increased spatial resolution scheme. In the Equatorial band, rates were considerably lower, ranging from below detection limit to ca. 4 nmol N L −1 d −1 , with a clear difference in magnitude between 2015 (rates close to zero) and 2016 (average rates around 2 nmol N L −1 d −1 ). This difference seemed triggered by a contrasting supply of phosphorus between years. Our study stresses the importance of surveys with sampling at fine-scale spatial resolution, and shows unexpected high variability in the rates of nitrogen fixation in the Guinea Dome, a region where diazotrophy is a significant process supplying new nitrogen into the euphotic layer.

Fernández-Carrera, Ana (ORCID:0000000183470996)↗

Pressure-temperature equation of state of Al 2 ⁢O 3 up to 14 Mbar and 40 kK

Sapphire (Al 2 ⁢O 3 ), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al 2 ⁢ O 3 as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al 2 ⁢ O 3 dictating a comprehensive overhaul. We develop two advanced EoSs for Al 2 ⁢ O 3 the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al 2 ⁢O 3 up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Kalita, Patricia [Sandia National Laboratories (SN↗

IDAES-PSE 2.6.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.

AS↗

Results of an interlaboratory study on the working curve in vat photopolymerization II: Towards a standardized method

The working curve measurement in photopolymer additive manufacturing is a ubiquitous measure of the cure depth of a printing resin as a function of radiant exposure of light. The fit parameters from this measurement (the depth of light penetration D p and the critical exposure E c ) are used to evaluate and report a resin’s printability, optimize processing parameters, and inform print and resin quality control. Despite its widespread use, the working curve lacks a standard measurement method. Here, following up on our paper “Results of an Interlaboratory Study on the Working Curve in Vat Photopolymerization” from last year, an interlaboratory study on the working curve was performed using calibrated, reproducible, bandpass filtered light sources. With these light sources, the variability between labs in measured working curves was dramatically reduced from the initial interlaboratory study. Aggregate data from this experiment produced reliable D p and E c measurements at 385 nm of 39.2 ± 3.7 µm and 12.3 ± 3.0 mJ cm −2 , respectively. At 405 nm the values of D p and E c are 69.3 ± 3.8 µm and 17.9 ± 2.3 mJ cm −2 , respectively. The results are agnostic to the thickness measurement tool utilized by participants, ensuring broad applicability across laboratories. We also tested the generalizability of the proposed method of using a filtered light source by filtering a commercial 405 nm light source and obtaining a working curve in agreement with the aggregate data from the interlaboratory study. This interlaboratory study provides a basis for a documentary standard for the working curve, so that the entire photopolymer additive manufacturing industry can share reproducible and interoperable working curve data.

36 MATERIALS SCIENCE↗

Emergence of pseudoresonance in high-intensity resonant inelastic x-ray scattering

We report resonant inelastic x-ray scattering (RIXS) spectra of neon atoms interacting with intense x-ray pulses generated using an x-ray free-electron laser (XFEL). We find that an unexpected peak emerges near the 𝐾⁢𝛼 line of Ne, which does not coincide with any physical resonances of neon ions. We perform theoretical calculations based on a quantum-state-resolved rate-equation approach with x-ray-induced processes including possible resonant excitations. Our dynamics simulations demonstrate that a sequence of multiple resonant photoabsorption events are involved and the interplay of those multiple resonances in combination with the relatively large spectral bandwidth of XFEL radiation leads to the emergent resonance-like structure at a position where no resonances exist. Our finding offers critical guidance for future applications of high-intensity RIXS at XFEL facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗