Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transitional flow”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

In-Cu alloy substrates for low-temperature chemical vapor deposition of Mo2C

Two-dimensional transition metal carbides are promising materials because of their potential for combining the favorable properties of transition metal carbides with the high aspect ratio of two-dimensional materials. Though commonly produced by top-down wet-chemical synthesis methods, synthesis by chemical vapor deposition is being considered because of its ability to achieve large areas, controlled layer thickness, and reduced defect density. Typically, liquid Cu is used as a synthesis substrate, though the high melting temperature of Cu (1085 °C) requires high synthesis temperatures. Thus, alternative substrates have been studied in order to reduce the necessary melting temperature. This work systematically studies the impact of synthesis parameters and substrate composition on the growth of ultrathin Mo2C (∼4–120 nm thick) by chemical vapor deposition on an In-Cu alloy. Mo2C flake size increases, and graphene/Mo2C heterostructures form with an increase in the methane flow rates. Increasing the In composition slightly decreases surface coverage and coalescence but does not appreciably impact the Mo2C flake size. Increasing In content also decreases the alloy substrate melting temperature so that a lower temperature synthesis (800 °C) can be performed. However, the necessary high temperatures for pyrolysis of methane lead to a lower limit for the synthesis temperature, similar to graphene.

Young, Katherine T. (ORCID:0000000285845576)↗
DOE OSTI2022

Experimental pub crawl from Rayleigh–Bénard to magnetostrophic convection

The interplay between convective, rotational and magnetic forces defines the dynamics within the electrically conducting regions of planets and stars. Yet their triadic effects are separated from one another in most studies, arguably due to the richness of each subset. In a single laboratory experiment, we apply a fixed heat flux, two different magnetic field strengths and one rotation rate, allowing us to chart a continuous path through Rayleigh–Bénard convection (RBC), two regimes of magnetoconvection, rotating convection and two regimes of rotating magnetoconvection, before finishing back at RBC. Dynamically rapid transitions are determined to exist between jump rope vortex states, thermoelectrically driven magnetoprecessional modes, mixed wall- and oscillatory-mode rotating convection and a novel magnetostrophic wall mode. Thus, our laboratory ‘pub crawl’ provides a coherent intercomparison of the broadly varying responses arising as a function of the magnetorotational forces imposed on a liquid-metal convection system.

42 ENGINEERING↗
DOE OSTI2023

Hydrate-bearing sands in the Terrebonne Basin record the transition from ponded deposition to bypass in the deep-water Gulf of Mexico

Here, we integrated well and seismic data to interpret the stratigraphic architecture and the paleogeographic evolution of two hydrate-bearing sand reservoirs in the Terrebonne Basin in the deep-water Gulf of Mexico. These sands were deposited during the final stage of ponded fill within the Terrebonne Basin. The Green sand filled the final accommodation in the basin as a perched apron; it was formed coincidentally with downstream erosion and sediment bypass which cut through a canyon formed on the southeastern edge of the basin. Subsequently, a low energy turbidite channel continued to bypass sediment through the basin and aggrade. The Orange sand was then deposited as a second perched apron and it was once again incised by the low energy turbidite channel. The Green and Orange sands are regional sheet sands that extend and thicken deep into the Terrebonne Basin. They provide a connected regional aquifer that captured and focused gas migration to the crest of the structure where hydrate was formed. These sands were formed by large, unconfined, turbidity flows. We envision that they will be composed of coarser grained deposits relative to previously explored hydrate reservoirs within levee deposits. As such, they have the potential to be relatively coarse grained and hence have higher permeability: a condition favorable to production of hydrate deposits.

3D seismic interpretation↗
DOE OSTI2025

Solidification cracking of laser melted commercial-purity tungsten

The high melting temperature of tungsten (W) makes it an attractive candidate for energy generation applications; however, its use is limited by its poor ductility at low temperatures. This limitation affects even melt-based additive manufacturing (AM) processes such as laser powder bed fusion (PBF-LB), with as-fabricated pure W exhibiting both longitudinal and transverse cracks. Metallurgical and processing factors that affect these cracks are still being explored. This work utilizes powderless single-tracks on pure W plates made over a range of power and velocity combinations in three unique PBF-LB setups: one with flowing argon, one with continuous vacuum, and an in situ high-speed synchrotron X-ray radiography setup. Each had different oxygen activities in their build environments. Longitudinal cracks with oxides appearing to exude from the crack were found for tracks deposited in a build environment with argon shield gas, while no such cracks were observed for builds conducted in vacuum. A calculation showed that direct oxygen ingress into the melt pool from the build environment is negligible when argon shield gas is present. Thus, incorporation of spatter is a likely mechanism for oxygen ingress into the melt pool. Radiography showed extensive keyholing, spatter generation, and crack formation at keyhole porosity. A heat transfer calculation showed the crack formation time was consistent with cooling below its ductile-to-brittle transition (DBT) temperature. This work identifies solidification cracking as a feasible mechanism in pure W beyond the well-known DBT related cracking.

Materials science↗
DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Universal Anomalous Diffusion of Quantized Vortices in Ultraquantum Turbulence

In classical viscous fluids, turbulent eddies are known to be responsible for the rapid spreading of embedded particles. However, in inviscid quantum fluids where the turbulence is induced by a chaotic tangle of quantized vortices, dispersion of the particles can be achieved via a nonclassical mechanism, i.e., their binding to the evolving vortices. However, knowledge on how the vortices diffuse and spread in quantum-fluid turbulence is very limited, especially for the so-called ultraquantum turbulence (UQT) generated by a random tangle of vortices. Here we report a systematic numerical study of the apparent diffusion of vortices in UQT in superfluid helium-4 using the full Biot-Savart simulation. We reveal that the vortices in the superfluid exhibit a universal anomalous diffusion (superdiffusion) at small times, which transits to normal diffusion at large times. This behavior is found to be the result of a generic scaling property of the vortex velocity. Our simulation at finite temperatures also nicely reproduces recent experimental observations. Lastly, the knowledge obtained from this study may form the base for understanding turbulent transport and universal vortex dynamics in various quantum fluids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Train small, model big: Scalable physics simulators via reduced order modeling and domain decomposition

Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the pilot-scale model, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. The ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. Here this method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about 15–40 times faster with only ~1% relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Time-Resolved X-ray Operando Observations of Lithiation Gradients across the Cathode Matrix and Individual Oxide Particles during Fast Cycling of a Li-Ion Cell

Lithiated transition metal oxides serve as active materials in the positive electrode (cathode) of lithium-ion cells. During electrochemical cycling, lithium ions intercalate and deintercalate into these oxide particles. This behavior causes two types of lithiation gradients to emerge: (i) a bulk gradient across the depth of the cathode matrix (averaged over individual oxide particles) and (ii) a microscopic gradient across the particles themselves, which also depends on their location in the electrode. Here we show how both gradients can be studied using operando X-ray diffraction during 4C charge and 4C discharge. The oxide (de)lithiation is estimated from the unit cell parameters by indexing the X-ray diffraction spectra. By fitting the lithiation profiles with orthogonal polynomials, the bulk gradients across the electrode thickness are quantified. These gradients develop as the current flows through the cell and dissipate during open-circuit and potentiostatic-hold periods. Further details of lithiation dynamics can be obtained through shape analysis of the Bragg peaks. In particular, from electrochemical model simulations, we show that the width and skewness of the (003) peak track (de)lithiation fronts moving across the individual oxide particles.

(003) peak↗
DOE OSTI2022

Start-up and Emergency Shutdown Modeling for a Coal-fired 10MWe sCO₂ Power Plant

As part of a Phase II Front End Engineering Design (FEED) study for the DOE Fossil Fuel Large Scale Pilots program, investigations of novel start-up and emergency shutdown methods for a stoker-fed, coal-fired, 10MWe sCO₂ power plant proposal were performed. These investigations were performed using 1D transient system models of the plant created within the GT-SUITE system modeling platform. The system model components were created using vendor quotation data (heat exchangers), performance maps (turbomachines), or FMU models (fired heater). Previously, Echogen Power Systems (Echogen) has used a sCO₂ power cycle start-up method that fills the system with liquid CO₂ to facilitate easy use of an electrically-driven start pump to transition the system from initial CO₂ fill to turbo-compressor initialization. The novel start-up method presented attempts this transition with a minimal filling of the system with liquid CO₂, to reduce the total CO₂ inventory required. This method leaves the start pump vulnerable to a two phase inlet condition as system pressure is below the CO₂ saturation pressure for the condenser cooling water temperature (Tcw). System model cases at Tcw of 18°C, 26°C, and 32°C, corresponding to cold, design, and hot ambient days, were analyzed to determine if the SP inlet condition could be kept as a subcooled liquid and estimate the CO₂ inventory reduction amount./p> Stoker-fed, coal-fired heaters continue to emit heat for some time even after emergency shutdown from events such as a power failure. This heat emission would lead to a failure of the fired heater, as the metal overheats, if the CO₂ flow is shut off. The emergency shutdown method presented utilizes CO₂ vented by a controllable vent valve (CRV) to provide CO₂ cooling flow to the fired heater. To determine the effectiveness of this method at keeping the fired heater peak tubing metal temperature below the ASME material temperature limit, for pressures below 5 MPa, of 816 °C, system model cases with CRV diameters ranging from 4” to 10”, and an alternative CO₂ vent routing with CRV size of 6”, were investigated./p>

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Mitigating Data Center Impact on Grid Stability: A Coordinated Control Strategy Using Verrus StabiliGrid Architecture

Large data centers, which now represent a significant and growing share of the total U.S. grid load, can inadvertently destabilize the electrical grid when they disconnect simultaneously during brief voltage disturbances. The July 10, 2024, Eastern Interconnection incident, in which a sub-100-millisecond transmission fault triggered the cascading loss of approximately 1,500 MW of data center load, illustrates this vulnerability. While commercial battery energy storage systems (BESS) deployed in data centers provide device-level fault ride-through per IEEE 1547, they lack coordination with facility protection logic and uninterruptible power supplies (UPS), limiting their effectiveness as grid-stabilizing assets. This report presents the Verrus StabiliGrid architecture, a coordinated control framework that integrates BESS, UPS, and point-of-interconnection (POI) protection settings to enable data centers to ride through both undervoltage and overvoltage grid contingencies without disconnecting. The four-step strategy encompasses: (1) high-resolution power quality monitoring to detect the grid state during events such as undervoltage, overvoltage, underfrequency, and overfrequency; (2) POI protection settings that allow for extended ride-through and grid-connected operation during grid contingencies; (3) grid state-driven autonomous dispatch of assets to improve grid resilience by reducing power draw during undervoltage or absorbing more power during overvoltage events; and (4) coordinated post-recovery dispatch of data center assets to restore firm load to pre-contingency levels. Validated through controller-hardware-in-the-loop (C-HIL) simulations at the National Laboratory of the Rockies, results show grid import restoration to pre-fault levels within 100 milliseconds of voltage recovery. This work advances the ability of data centers to transition from passive, disturbance-sensitive loads to active participants in grid stability, a capability increasingly required by emerging NERC and ERCOT regulatory frameworks.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Influence of Surface Defects on WO 3 Photoelectrodes for Catalyzing Chloride Oxidation in Water

Tungsten oxide (WO 3 ) is an n-type semiconductor due to oxygen vacancies (□ O •• in Kroger-Vink notation) or surface protonation as H x WO 3 . It is one of the few acid-stable oxides under large positive bias, which makes WO 3 ideal for interrogating the mechanism of the chloride oxidation reaction (COR). The large, positive valence band edge of ∼3 eV provides the overpotential necessary to carry out the COR, but the reaction competes with the oxygen-evolution reaction in water. The □ O •• defect density can be controlled by the atmosphere under which the material is annealed, so WO 3 films were prepared by a spin-coating method from an ammonium metatungstate precursor annealed at 500 °C under air, flowing O 2 , and flowing argon. Annealing the films in a flowing O 2 atmosphere hinders the formation of □ O ••, and annealing in Ar leads to greater surface W 6+ , likely due to expelling intercalated H + . The saturated photocurrent density (j ph ) is highest in films with the greatest concentration of W 5+ and greatest concentration of oxide defects: (0.66 mA/cm 2 annealed in air, 0.58 mA/cm 2 annealed in Ar, and 0.49 mA/cm 2 annealed in O 2 , reported at 1.5 V vs Ag/AgCl, pH 3 (before the onset of a dark reaction). The defect concentrations are determined by X-ray photoelectron spectroscopy. In all cases the Faradaic efficiency for the COR is near unity. Finally, we demonstrate that W 5d states can be probed by ligand K-edge X-ray absorption near-edge spectroscopy via pre-edge (Cl 1s → W 5d) transitions, lower in energy than the ligand-centered (Cl 1s → 4p) transition. We use this analysis to show the presence of W─Cl covalent bonds on the WO 3 films post-COR, corroborated by DFT calculations. Furthermore, this result stands in contrast to the commonly assumed mechanistic proposal invoking outer-sphere electron transfer to a physisorbed chloride ion.

Cl K-edge↗
DOE OSTI2021

Decoupling Conductivity and Solubility in Electrolytes Using Microemulsions

Redox flow batteries have recently received considerable attention as possible large-scale energy storage devices, but their low energy density has inhibited widespread application. In this work, a novel strategy of decoupling conductivity and solubility of electrolytes using microemulsion is put forward to enhance ionic conduction of non-aqueous electrolytes, increase the selectivity of active species, improve the battery voltage, and eventually achieve the possibility of high energy density. We report a study of the electrochemistry of ferrocene in single phase Tween® 20/1-butanol/H 2 O/toluene microemulsion system at 20 °C. At low and intermediate surfactant to water weight ratios (<0.5/0.5), the voltammogram exhibits reversible electrochemical behavior, while at high surfactant levels the curves show lower levels of reversibility. The latter voltammograms have a form typically associated with high resistance in solution, consistent with a gradual transition in microstructure as surfactant levels increase. This change in structure is supported by correlations with conductivity results based on the literature. The voltammograms show little evidence of anomalies in double layer capacitance or electrode "blocking" by droplets, suggesting that the electron transfer is indeed occurring in a facile manner between the electrode and the ferrocene inside the oil phase.

25 ENERGY STORAGE↗
DOE OSTI2024

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗
DOE OSTI2024

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modified equilibrium distributions for Cooper-Frye particlization

We introduce a positive-definite single-particle distribution that is suitable for describing the transition from a macroscopic hydrodynamic to a microscopic kinetic description during the late stages of heavy-ion collisions in the presence of moderately large viscous corrections. The modified equilibrium distribution function can be constructed with hydrodynamic input from either relativistic viscous fluid dynamics or anisotropic fluid dynamics. We test the modified equilibrium distribution's hydrodynamic output for a stationary hadron resonance gas subject to either shear stress, bulk pressure, or baryon diffusion current at a given freeze-out temperature and baryon chemical potential. While it does not reproduce all components of the net baryon current and energy-momentum tensor exactly, it significantly improves upon the customary linearized approximations for the nonequilibrium correction δf n which typically lead to unphysical negative distribution functions at large particle momenta. A comparison of particle spectra and p T –differential elliptic flow coefficients from the Cooper–Frye formula computed with the modified equilibrium distribution and with linearized δf n corrections is presented, for two different (2 + 1)–dimensional hypersurfaces corresponding to central and noncentral Pb + Pb collisions at the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗