SEARCH · Search NASA

Results for “Thermal spikes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Thermo-Hydro-Mechanical Modeling of Brine Migration in a Heated Borehole Test in Bedded Salt

Abstract This research paper focuses on the thermo-hydro-mechanical (THM) modeling of brine migration in a heated borehole test conducted as part of the ongoing Brine Availability Test in Salt (BATS) at the Waste Isolation Pilot Plant (WIPP) in New Mexico. It is a component of the international collaboration project DECOVALEX-2023 (DEvelopment of COupled models and their VALidation against EXperiments), which aims to understand the THM processes governing brine flow in heated rock salt repositories through collaborative analysis by multiple research teams. Using the TOUGH–FLAC simulator, THM simulations were performed and compared with data from the BATS phase 1a. This experiment involved two identical horizontal-borehole arrays, one heated and one serving as a control, both equipped with sensor arrays. Analysis of measurements revealed water flow rate surges during heater power transitions, with the highest jump observed during cooling. Acoustic emission activity exhibited distinct patterns in response to heater power changes, suggesting that damage to rock salt is particularly pronounced during the cooling phase. The THM simulations successfully captured these phenomena, highlighting the significance of thermal effects, brine migration, and mechanical behavior in predicting brine availability in heated and damaged rock salt. Our modeling also revealed the critical interplay between heating and cooling-induced damage and its influence on flow properties, particularly affecting brine inflow estimation. Notably, we found that cooling-induced brine inflow spikes result from increased permeability due to tensile dilatancy. These findings have important implications for the development of robust containment strategies and enhance our understanding of the complex processes involved in repository performance.

58 GEOSCIENCES↗

DOE OSTI2023

Proton nuclear magnetic resonance (1H NMR) of flammable organic chemicals in radioactive high–level supernatant waste at the Savannah River Site (SRS)

The Savannah River Site stores approximately 36 million gallons of radioactive and hazardous waste that contains approximately 245 million curies. The waste is sent through various chemical processes to reduce its volume and to separate various components. The facility plans to replace formic acid (a chemical used to reduce soluble mercury) with glycolic acid. Recycle solution with glycolate may flow back to the tank farm, where the glycolate can generate hydrogen gas by thermal and radiolytic mechanisms. The current analytical method for detecting glycolate (ion chromatography) in supernatant requires a large dilution to reduce interference from the nitrate anions. Hydrogen nuclear magnetic resonance is an analytical method that requires less sample dilution. It takes advantage of the CH 2 group in glycolate. Liquid samples were spiked with four different levels of glycolate to build a calibration line, as it is recommended in the standard addition method. The detection and quantitation limits determined were 1 and 5 ppm, respectively, for 32 scans, which is well below the process limit of 10 ppm. In one test, 800 scans of a supernatant spiked with 1 ppm glycolate resulted in a -CH 2 peak with a signal-to-noise ratio of 36.

1H↗

DOE OSTI2022

Switching Dynamics in Vanadium Dioxide-Based Stochastic Thermal Neurons

We report on switching dynamics of individual and coupled vanadium dioxide (VO 2 ) devices subject to voltage pulses as the temperature is systematically varied from room temperature spanning the insulator–metal transition (IMT) temperature. The switching voltage of single devices has a strong relationship with both temperature and voltage pulsewidth. Two-step switching in connected VO 2 devices has been noted in current transient plots and was found to depend on temperature, pulsewidth, and pulse amplitude. Experimental switching behavior measured from VO 2 artificial neurons was implemented into a spiking neural network (SNN). During training, modulating the switching voltage via temperature affords a novel method to implement homeostasis with the coupled devices. Simulation results show the efficacy of the stochastic neuronal characteristics and the proposed homeostasis mechanism on a standard digit recognition task. As a result, these studies contribute to ongoing efforts in neuromorphic computing exploiting collective phase transitions.

42 ENGINEERING↗

DOE OSTI2021

Low-temperature (< 200 °C) degradation of electronic nicotine delivery system liquids generates toxic aldehydes

Electronic cigarette usage has spiked in popularity over recent years. The enhanced prevalence has consequently resulted in new health concerns associated with the use of these devices. Degradation of the liquids used in vaping have been identified as a concern due to the presence of toxic compounds such as aldehydes in the aerosols. Typically, such thermochemical conversions are reported to occur between 300 and 400 °C. Herein, the low-temperature thermal degradation of propylene glycol and glycerol constituents of e-cigarette vapors are explored for the first time by natural abundance 13 C NMR and 1 H NMR, enabling in situ detection of intact molecules from decomposition. The results demonstrate that the degradation of electronic nicotine delivery system (ENDS) liquids is strongly reliant upon the oxygen availability, both in the presence and absence of a material surface. When oxygen is available, propylene glycol and glycerol readily decompose at temperatures between 133 and 175 °C over an extended time period. Among the generated chemical species, formic and acrylic acids are observed which can negatively affect the kidneys and lungs of those who inhale the toxin during ENDS vapor inhalation. Further, the formation of hemi- and formal acetals is noted from both glycerol and propylene glycol, signifying the generation of both formaldehyde and acetaldehyde, highly toxic compounds, which, as a biocide, can lead to numerous health ailments. The results also reveal a retardation in decomposition rate when material surfaces are prevalent with no directly observed unique surface spectator or intermediate species as well as potentially slower conversions in mixtures of the two components. The generation of toxic species in ENDS liquids at low temperatures highlights the dangers of low-temperature ENDS use.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Rheology and engine performance of very low sulfur fuel oil blended with 10% fast pyrolysis and hydrothermal liquefaction oils in a 2-stroke crosshead engine

The performance and emissions for a downscaled single-cylinder 2-stroke crosshead engine were determined for a very low sulfur fuel oil (VLSFO) when blended with 10 wt.% fast pyrolysis (FP) or hydrothermal liquefaction (HTL) bio-intermediates. The FP and HTL oils were derived from biomass and were observed to contain lower molecular weight (MW) hydrocarbons than neat VLSFO (which was evaluated as a baseline comparison). The addition of either biofuel reduced the overall viscosity of the VLSFO. Aging tests at 50, 90, and 120°C showed that the dynamic viscosity of VLSFO increased with exposure time up to two weeks. Similar trends were observed for the FP and HTL blends, but a pronounced spike in viscosity occurred for these fuels during the early period of exposure. None of the viscosity increases exceeded the operational limits of fuel system pumps. Engine performance studies were conducted under low, medium and high load operational settings. The relative performance of the test fuels was highly dependent on operating condition. In general, the engine results for the three test fuels were similar, but modest improvements in brake thermal efficiency and brake specific fuel consumption were observed, which may be attributed to the heightened reactivity of low molecular weight fraction of the FP and HTL oils.

09 BIOMASS FUELS↗

DOE OSTI2023

Intrinsic and Extrinsic Factors Influencing the Dynamics of VO 2 Mott Oscillators

Oscillatory devices have gained significant interest recently as key components of computing systems based on biomimetic neuronal spiking. An understanding of the time scales underlying the spiking is essential for engineering fast, controllable, low energy devices. However, we find that the intrinsic dynamics of these devices are difficult to properly characterize, as they can be heavily influenced by the external circuitry used to measure them. Here we demonstrate these challenges using a VO 2 Mott oscillator with a sub-100 nm effective size, achieved using a nanogap cut in a metallic carbon nanotube electrode. Given the nanoscale thermal volume of this device, it would be expected to exhibit rapid oscillations. However, due to external parasitics present within commonly used current sources, we see orders of magnitude slower dynamics. Here, we outline methods for determining when measurements are dominated by extrinsic factors and discuss the operating conditions under which intrinsic oscillation frequencies may be observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Hydrothermal liquefaction of sewage sludge – effect of four reagents on relevant parameters related to biocrude and $\mathrm{PFAS}$

Catalytic hydrothermal liquefaction (HTL) can convert wet biomass into crude-like oil and has potential to degrade pollutants during the thermal depolymerization process. To avoid negative environmental impacts from per- and polyfluoroalkyl substances (PFAS) in sewage sludge after its final disposal, we comprehensively evaluated catalytic destruction of PFAS in sludge through HTL. Here the results showed that adding red mud to HTL was beneficial for the biocrude yield and significantly increased the carbon conversion efficiency and energy recovery of the hydrothermal process. HTL at 300 °C for 2 h degraded > 96% of spiked perfluorooctanoic acid (PFOA). The mass of spiked perfluorobutanesulfonic acid (PFBS) and perfluorooctanesulfonic acid (PFOS) and pre-existing perfluorohexanoic acid (PFHxA) and perfluoroheptanoic acid (PFHpA), however, increased in the HTL product streams. The remaining PFAS after HTL mainly stayed in the biocrude phase. Our findings indicated that HTL under the studied condition was not able to degrade all PFAS existed in sewage sludge. Although adding red mud could improve the HTL performance with respect to biocrude, this reagent did not improve the removal of perfluoroalkyl acids (PFAAs). As a matter of fact, the hydrothermal process amended with the selected reagents resulted in increased mass of PFBS, PFHxA, and PFHpA in the HTL products due to degradation of PFAA precursors. Thus, other treatment approaches for PFAS removal from sewage sludge need to be identified to prevent possible contamination of environments receiving PFAS-containing sludge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Exponential Escape Rate of Filamentary Incubation in Mott Spiking Neurons

Mott materials such as vanadium oxides, when subject to a strong applied voltage, present an inhomogeneous insulator-to-metal transition with formation of metallic filaments within the insulating bulk. This property is enabling the development of compact and power-efficient neuromorphic devices known as Mott neurons. However, the nature of the transition has not been fully understood yet, as it may be attributed to different effects, including Joule self-heating and hot-carrier injection. Moreover, the experimental determination of the threshold voltage needed to induce the transition has proven to be challenging, as the transition becomes increasingly unpredictable when the threshold is approached. The physical understanding of these issues would not only deepen our understanding of Mott insulators, but would also be an important step toward the realization of neuromorphic devices based on such materials. In this work we use numerical simulations based on the Mott resistor network model to study the nature of the filament incubation and formation process. Here, we show that both electronic and thermal effects, in the form of current density focusing and Joule self-heating, respectively, contribute to the filamentary incubation and growth. Remarkably, we find that the percolation of the metallic filaments near the threshold is intrinsically stochastic, qualitatively similar to the familiar Arrhenius activated behavior and to the stochastic firing of biological neurons. More precisely, we characterize the filament percolation as a Poisson point process, which has the same probability distribution as mathematical models of neuronal firing with an exponential escape rate. Finally, we support the numerical simulation results by performing experiments in VO 2 that are in agreement with the exponential escape rate behavior. Thus, we establish a functionality of Mott insulators that opens a path toward implementing neuromorphic hardware with quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

NUMERICAL MODELING OF A SOLID OXIDE FUEL CELL FOR USE IN REAL-TIME SIMULATION AND CYBER-PHYSICAL SYSTEMS

Cyber-physical systems provide a mechanism with which to investigate the physical phenomena and behavior of traditionally cost-prohibitive or otherwise fragile equipment. For the National Energy Technology Laboratory (NETL), this approach resulted in the Hybrid Performance (Hyper) facility which features a gas turbine-SOFC hybrid cycle utilizing real turbomachinery and a simulated SOFC stack. This allows for the investigation of combined cycle performance and control strategies, in an exhaustive manner, both without fear of destroying delicate state-of-the-art fuel cells, and with the full accuracy of real-world turbomachinery. Issues arose between the transient response of the SOFC model being limited to a sample time of 80 milliseconds, due to the calculation time of the SOFC model taking on average 40 milliseconds to calculate for a given timestep with spikes in calculation time reaching the 80 millisecond threshold. In order to be able to match the speed of transients from the turbomachinery and likewise better discern transient behavior, it was determined that the SOFC model must be optimized to operate at a sample time of 5 milliseconds. Therefore, it is necessary to optimize the SOFC model in order to decrease the calculation time from around 40 milliseconds, down to at the most 5 milliseconds. To do this, both the electrochemical algorithm and the thermal algorithm used to simulate the physical behavior of the SOFC are investigated to determine where improvements can be made. To this end the rootfinding numerical recipes of the electrochemical algorithm are investigated as the complex electrochemistry requires a highly iterative nested dual convergence loop to resolve the voltage-current relationship, and likewise the temporal discretization of the thermal algorithm is modified for the sake of higher accuracy and stability. Ultimately the new electrochemical algorithm featuring higher order rootfinding schemes proves to be efficient enough to reach the sub 5 millisecond target, signifying an order of magnitude reduction in calculation time, and when coupled with the new temporal discretization similar calculation time characteristics show that a fully implicit, higher order temporal discretization can also successfully be used if desired. Ultimately this result means that the cyber-physical simulation system can operate at higher sample rates, and resolve transient events at significantly higher resolution and fidelity.

Arias, Jesus↗

DOE OSTI2020

Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks

Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

A domain wall-magnetic tunnel junction artificial synapse with notched geometry for accurate and efficient training of deep neural networks

Inspired by the parallelism and efficiency of the brain, several candidates for artificial synapse devices have been developed for neuromorphic computing, yet a nonlinear and asymmetric synaptic response curve precludes their use for backpropagation, the foundation of modern supervised learning. Spintronic devices—which benefit from high endurance, low power consumption, low latency, and CMOS compatibility—are a promising technology for memory, and domain-wall magnetic tunnel junction (DW-MTJ) devices have been shown to implement synaptic functions such as long-term potentiation and spike-timing dependent plasticity. In this work, we propose a notched DW-MTJ synapse as a candidate for supervised learning. Using micromagnetic simulations at room temperature, we show that notched synapses ensure the non-volatility of the synaptic weight and allow for highly linear, symmetric, and reproducible weight updates using either spin transfer torque (STT) or spin–orbit torque (SOT) mechanisms of DW propagation. We use lookup tables constructed from micromagnetics simulations to model the training of neural networks built with DW-MTJ synapses on both the MNIST and Fashion-MNIST image classification tasks. Accounting for thermal noise and realistic process variations, the DW-MTJ devices achieve classification accuracy close to ideal floating-point updates using both STT and SOT devices at room temperature and at 400 K. Our work establishes the basis for a magnetic artificial synapse that can eventually lead to hardware neural networks with fully spintronic matrix operations implementing machine learning.

42 ENGINEERING↗

DOE OSTI2022

Transmission Electron Microscopy Study on the Effect of Thermal and Electrical Stimuli on Ge2Te3 Based Memristor Devices

Memristor devices fabricated using the chalcogenide Ge 2 Te 3 phase change thin films in a metal-insulator-metal structure are characterized using thermal and electrical stimuli in this study. Once the thermal and electrical stimuli are applied, cross-sectional transmission electron microscopy (TEM) and X-ray energy-dispersive spectroscopy (XEDS) analyses are performed to determine structural and compositional changes in the devices. Electrical measurements on these devices showed a need for increasing compliance current between cycles to initiate switching from low resistance state (LRS) to high resistance state (HRS). The measured resistance in HRS also exhibited a steady decrease with increase in the compliance current. High resolution TEM studies on devices in HRS showed the presence of residual crystalline phase at the top-electrode/dielectric interface, which may explain the observed dependence on compliance current. XEDS study revealed diffusion related processes at dielectric-electrode interface characterized, by the separation of Ge 2 Te 3 into Ge- and Te- enriched interfacial layers. This was also accompanied by spikes in O level at these regions. Furthermore, in-situ heating experiments on as-grown thin films revealed a deleterious effect of Ti adhesive layer, wherein the in-diffusion of Ti leads to further degradation of the dielectric layer. This experimental physics-based study shows that the large HRS/LRS ratio below the current compliance limit of 1 mA and the ability to control the HRS and LRS by varying the compliance current are attractive for memristor and neuromorphic computing applications.

Shallcross, Austin↗

DOE OSTI2025

Heating of Cs2Te photocathode via field emission and radiofrequency pulsed heating: Implication toward breakdown

The occurrence of radiofrequency (RF) breakdown limits operational electromagnetic gradients in accelerator structures. Experimental evidence often suggests that breakdown events are associated with temperature and dark current spikes on the surface of RF devices. In the past decade, there has been increased interest in unveiling the mechanism behind breakdown initiation in metal copper and copper alloys; however, efforts regarding a breakdown phenomenon in photocathode-relevant semiconductors have been more limited. In this work, we explore field-emission-assisted heating via Nottingham and Joule processes, as a possible candidate for breakdown initiation. For this, field emission from an intrinsic Cs2Te ultrathin film coated on a copper substrate was modeled within the Stratton–Baskin–Lvov–Fursey formalism, describing the processes and effects in the bulk and on the surface of a photocathode exposed to high RF electromagnetic fields. It is shown that a field-emission characteristic deviates significantly from the classical Fowler–Nordheim (FN) theory, whereby predicting that dark current is orders of magnitude lower than the one expected by FN law. Conventional pulsed heating was also found to impose negligible heating to the photocathode. Both conclusions suggest that a Cs2Te photocathode coated on a metal substrate would be insensitive to catastrophic thermal-material runaway breakdown, unlike what is observed for metal surfaces. Finally, a few unconventional breakdown candidate scenarios are identified and discussed, including thermoelastic deformation and avalanche breakdown.

Shinohara, Ryo (ORCID:0000000276995732)↗

DOE OSTI2021

Chemiresistor sensor based on ion-imprinted polymer (IIP)-functionalized rGO for Cd(II) ions in water

This study reports the design and development of a novel chemiresistor (CR) sensor using ion imprinted polymer (IIP)-functionalized reduced graphene oxide (rGO) [IIP/rGO-CR] for cadmium ions (Cd(II)) determination in water. The sensor consisted of a CR transducer made of rGO channel bridging source and drain electrodes prepared by self-assembly and thermal reduction of graphene oxide (GO) on Au interdigitated electrodes chip fabricated on Si/SiO 2 substrate. The IIP was then grafted on rGO using surface-initiated reversible addition-fragmentation chain transfer (RAFT) polymerization with polyethylenimine (PEI) and methylacrylic acid (MAA) as dual functional monomers and Cd(II) ions as template through UV light-initiated copolymerization. The IIP functionalized on rGO acted as an effective recognition element that modulated the resistance of rGO-CR upon binding of Cd(II), enabling Cd(II) detection at ppb level in aqueous solutions. The prepared IIP/rGO-CR sensor worked effectively in the linear range of 2~200ppb and achieved a limit of detection (LOD) of 0.83 ppb, which is lower than the World Health Organization guidelines of 3ppb for drinking water quality. The developed sensor of IIP/rGO-CR showed a high selectivity against a variety of trace and heavy metal ions found in water and good stability for up to 60 days when stored at room temperature for Cd(II) determination in water. Further, the sensor was successfully applied to analyzing Cd(II) spiked in tap, lake and river waters with a 94.5%–113.5% recovery, demonstrating a high degree of accuracy even in complex water samples. Finally, our results illustrated that the CR sensor of IIP functionalized rGO provides a potential platform for sensitive, robust and low-cost environmental analysis of Cd(II) in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗

DOE OSTI2024

Electron Heating in 2D Particle-in-cell Simulations of Quasi-perpendicular Low-beta Shocks

Abstract We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta ( β p = 0.25) collisionless ion–electron shocks with mass ratio m i / m e = 200, fast Mach number M ms = 1 –4, and upstream magnetic field angle θ Bn = 55°–85° from the shock normal n ˆ . It is known that shock electron heating is described by an ambipolar, B -parallel electric potential jump, Δ ϕ ∥ , that scales roughly linearly with the electron temperature jump. Our simulations have Δ ϕ ∥ / ( 0.5 m i u sh 2 ) ∼ 0.1 –0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measure ϕ ∥ , including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate of ϕ ∥ in our low- β p shocks. We further focus on two θ Bn = 65° shocks: a M s = 4 ( M A = 1.8 ) case with a long, 30 d i precursor of whistler waves along n ˆ , and a M s = 7 ( M A = 3.2 ) case with a shorter, 5 d i precursor of whistlers oblique to both n ˆ and B ; d i is the ion skin depth. Within the precursors, ϕ ∥ has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the M s = 4 , θ Bn = 65° case, ϕ ∥ shows a weak dependence on the electron plasma-to-cyclotron frequency ratio ω pe /Ω ce , and ϕ ∥ decreases by a factor of 2 as m i / m e is raised to the true proton–electron value of 1836.

Tran, Aaron (ORCID:0000000334834890)↗

DOE OSTI2020

Effects of Turbulence and Temperature Fluctuations on Knock Development in an Ethanol/Air Mixture

The effects of turbulence on knock development and intensity for a thermally inhomogeneous stoichiometric ethanol/air mixture at a representative end-gas autoignition condition in internal combustion engines are investigated using direct numerical simulations with a skeletal reaction mechanism. Two- and three-dimensional simulations are performed by varying the most energetic length scale of temperature, $l_T$ , and its relative ratio with the most energetic length scale of turbulence, $l_T/l_e$ , together with two different levels of the turbulent velocity fluctuation, $u'$ . It is found that $l_T$ / $l_e$ and the ratio of ignition delay time to eddy-turnover time, $\tau _{ig}/\tau _t$ , are the key parameters that control the detonation development. An increase in either $l_T$ or $l_e$ enhances the detonation propensity by allowing a longer run-up distance for the detonation development. The characteristic length scale of the temperature field, $l_T$ , is significantly modified by high turbulence intensity achieved by a large $l_e$ and $u'$ . The intense turbulence mixing effectively distributes the initial temperature field to broader scales to support the developing detonation waves, thereby increasing the likelihood of the detonation formation. On the contrary, high turbulence intensity with a short mixing time scale, achieved by a small $l_e$ and a large $u'$ , reduces the super-knock intensity attributed to the finer broken-up structures of detonation waves. Either $\tau _{ig}/\tau _t$ less than unity or $l_e = l_T$ even with a large $u'$ is found to have no significant effect on super-knock mitigation. Finally, high turbulent intensity may induce high-pressure spikes comparable to the von Neumann spike. Overall, increased temperature and pressure by combustion heating, noticeably after the peak of heat release rate, significantly enhance the collision and interaction of multiple emerging autoignition fronts near the ending combustion process, resulting in localized high-pressure spikes.

42 ENGINEERING↗