Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal fluid physics”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

GaAs growth rates of 528 μm/h using dynamic-hydride vapor phase epitaxy with a nitrogen carrier gas

Herein we demonstrate record GaAs growth rates approaching 530 μm/h using nitrogen carrier gas and 400 μm/h using hydrogen carrier gas in a dynamic-hydride vapor phase epitaxy reactor. We measured root mean square surface roughness below 1 nm using a 1 μm × 1 μm atomic force microscopy scan for GaAs growth rates up to 483 μm/h using a nitrogen carrier gas and 400 μm/h using a hydrogen carrier gas. We performed computational fluid dynamics modeling to study the effect of the carrier gas choice on the thermal profile within the reactor and how that influences the degree of AsH 3 decomposition. The modeling suggests that the lower thermal conductivity of the nitrogen carrier gas minimizes the amount of AsH 3 that thermally decomposes before reaching the wafer surface and the heavier atomic mass decreases the likelihood that AsH 3 will reach the heated reactor walls, leading to a growth rate enhancement relative to the hydrogen carrier case in a hydride-enhanced growth regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Generation of localized reactor point kinetics parameters using coupled neutronic and thermal fluid models for pebble-bed reactor transient analysis

The systems analysis of anticipated operating occurrences and design basis accidents for pebble-bed reactor systems requires knowledge of neutron point kinetics equations (PKE) parameters. Typically, the generation of PKE parameters is performed in a global manner using standalone neutronics calculations, without the inclusion of thermal fluid distributions. We utilize Griffin and Pronghorn for generating global and local PKE parameters which includes the use of thermal fluid distributions to account for localized effects. This work establishes a methodology for calculating PKE parameters for a pebble bed reactor with a coupled neutronics/thermal fluids analysis. PKE parameters generated on a global and local basis are compared against a diffusion solve for a typical load-following transient. Locally-generated neutron kinetic parameters are able to reduce the maximum error in the transient power level from 5% to below 1.5%; along with this, bulk temperature errors were reduced from 12 K to 4 K.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Observation of shock-front separation in multi-ion-species collisional plasma shocks

We observe shock-front separation and species-dependent shock widths in multi-ion-species collisional plasma shocks, which are produced by obliquely merging plasma jets of a He/Ar mixture (97% He and 3% Ar by initial number density) on the plasma liner experiment. Visible plasma emission near the He-I 587.6-nm and Ar-II 476.5–514.5-nm lines is simultaneously recorded by splitting a single visible image of the shock into two different fast-framing cameras with different narrow bandpass filters (589 ± 5 nm for observing the He-I line and 500 ± 25 nm for the Ar-II lines). For conditions in these experiments (pre-shock ion and electron densities ≈5×10 14 cm –3 , ion and electron temperatures of ≈2.2 eV, and relative plasma-merging speed of 22 km/s), the observationally inferred magnitude of He/Ar shock-front separation and the shock widths themselves are <1 cm, which correspond to ~50 post-shock thermal ion–ion mean free paths. These experimental length scales are in reasonable qualitative and quantitative agreement with results from 1D multi-fluid simulations using the CHICAGO code. However, we note there are differences between the experimentally inferred and simulation-predicted ionization states and line emission intensities, particularly in the post-shock region. Overall, the experimental and simulation results are consistent with theoretical predictions that the lighter He ions diffuse farther ahead within the overall shock front than the heavier Ar ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Preliminary Thermal-Hydraulic Analysis of the UTA–2 Subcritical Linear Accelerator–Driven System

The National Nuclear Security Administration’s mission includes establishing a reliable supply of 99 Mo without highly enriched uranium. Oak Ridge National Laboratory (ORNL) supports this objective through collaborative research and development with industrial partners. Niowave Inc., a current partner, is currently designing a subcritical linear accelerator-driven system (ADS) and preparing for the US Nuclear Regulatory Commission’s licensing process. Niowave’s system has the potential to efficiently supply medical radioisotopes. The technology includes a superconducting electron accelerator and a pile of both natural and low-enriched uranium (LEU) targets. The process fissions uranium and many valuable isotopes can then be extracted from the targets. Niowave is currently iterating through conceptual and detailed design processes for several system sizes. This report discusses UTA–2, which is at the demonstration stage. UTA–2 will validate numerical modeling results with experimental measurements before progressing to the detailed design of UTA–3, the commercial-sized ADS. The thermal-hydraulic behavior of the UTA–2 core design was numerically investigated using STAR-CCM+, as described in this report. STAR-CCM+, a state-of-the-art computational fluid dynamics (CFD) software that was commercially developed by Siemens, has an extensive user base and a set of validation studies. It is also compliant with the American Society of Mechanical Engineers’ Nuclear Quality Assurance 1 standard. Three cases are investigated in this report. Case 1 quantified the temperature field within the UTA–2 assembly and water tank using only conduction as the method of thermal energy transport. This simplified approach was overly conservative and yielded wetted cladding temperatures above the coolant saturation temperature. In this case, the maximum temperature of the wetted cladding surface of the highest power rod exceeded the saturation temperature of water by 63.8°C. Because of the overly conservative approach taken in Case 1 and its negative subcooled margin, buoyancy-driven natural circulation flow physics were implemented in Case 2. Adding coolant motion significantly distributed the thermal energy of the system through convective heat transfer. This relatively small amount of convective heat transfer significantly reduced system temperatures and increased the subcooled margin from –63.8 to 61.3°C. This margin confirmed that no boiling was expected during normal operation of UTA–2 at 230 W. Case 3 had no additional physics models but considered an overpower event during which the power of each LEU and natural uranium rod was at its respective peak values. This resulted in a study with the total assembly power equal to 176% of the nominal power of 230 W considered in Cases 1 and 2. The resulting natural circulation flows were slightly enhanced. The subcooled margin decreased slightly to 46.3°C, which is still a significant margin to local boiling of the water within the tank. This margin confirmed that no boiling was expected during an abnormal operation of UTA–2 at 406 W.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Did we hear the sound of the Universe boiling? Analysis using the full fluid velocity profiles and NANOGrav 15-year data

Abstract In this paper, we analyse sound waves arising from a cosmic phase transition where the full velocity profile is taken into account as an explanation for the gravitational wave spectrum observed by multiple pulsar timing array groups. Unlike the broken power law used in the literature, in this scenario the power law after the peak depends on the macroscopic properties of the phase transition, allowing for a better fit with pulsar timing array (PTA) data. We compare the best fit with that obtained using the usual broken power law and, unsurprisingly, find a better fit with the gravitational wave (GW) spectrum that utilizes the full velocity profile. Even more importantly, the thermal parameters that produce the best fit are quite different. We then discuss models that can produce the best-fit point and complementary probes using CMB experiments and searches for light particles in DUNE, IceCUBE-Gen2, neutrinoless doubleβ-decay, and forward physics facilities (FPF) at the LHC like FASERν, etc.

Astronomy & Astrophysics↗
DOE OSTI2025

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Gas-Phase Temperature Mapping of Evaporating Microdroplets

Evaporation is a ubiquitous and complex phenomenon of importance to many natural and industrial systems. Evaporation occurs when molecules near the free interface overcome intermolecular attractions with the bulk liquid. As molecules escape the liquid phase, heat is removed, causing evaporative cooling. The influence of evaporative cooling on inducing a temperature difference with the surrounding atmosphere as well as within the liquid is poorly understood. Here, we develop a technique to overcome past difficulties encountered during the study of heterogeneous droplet evaporation by coupling a piezo-driven droplet generation mechanism to a controlled micro-thermocouple to probe microdroplet evaporation. The technique allowed us to probe the gas-phase temperature distribution using a micro-thermocouple (50 μm) in the vicinity of the liquid–vapor interface with high spatial (±10 μm) and temporal (±100 ms) resolution. We experimentally map the temperature gradient formed surrounding sessile water droplets having varying curvature dictated by the apparent advancing contact angle (100° ≲ θ ≲ 165°). The experiments were carried out at temperatures below and above ambient for a range of fixed droplet radii (130 μm ≲ R ≲ 330 μm). Our results provide a primary validation of the centuries-old theoretical framework underpinning heterogeneous droplet evaporation mediated by the working fluid, substrate, and gas thermophysical properties, droplet apparent contact angle, and droplet size. We show that microscale droplets residing on low-thermal-conductivity substrates such as glass absorb up to 8× more heat from the surrounding gas compared to droplets residing on high-thermal-conductivity substrates such as copper. Finally, our work not only develops an experimental understanding of the heat transfer mechanisms governing droplet evaporation but also presents a powerful platform for the study and characterization of liquid–vapor transport at curved interfaces wetting and nonwetting advanced functional surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Temperature-Dependent Speciation of Ni(II) in Molten Chloride, Bromide, and Iodide Salts

Understanding the fundamental speciation and local structure of metal cations─such as nickel, which is present as a corrosion product─in molten salt media is essential for assessing their impact on the thermal and physical properties of salts used for advanced nuclear energy applications. In this study, we employed an integrated approach of combining X-ray absorption spectroscopy (XAS), UV–vis spectroscopy, and ab initio molecular dynamics (AIMD) to investigate the halide-dependent speciation of divalent nickel ions, Ni(II), in LiCl–KCl, LiBr–KBr, and LiI–KI eutectic salt mixtures from room temperature to 600 °C. Our findings show that both halide type and temperature significantly influence the local coordination of nickel in these molten salt systems. A notable aspect of Ni(II) speciation in LiI–KI salts is the formation of a polyiodide species. In conclusion, our findings provide critical insights into metal ion speciation in molten salt systems, enhancing our understanding of the factors that govern the physicochemical properties of these complex fluids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Recent progress in the design of the K-DEMO divertor

The preliminary conceptual design of the Korean fusion demonstration reactor (K-DEMO) with a major radius of 6.8 m and the fusion power of 2200 MW has been studied since 2012. The overall configuration of the K-DEMO divetor system based on the ITER-like water-cooled tungsten technology is a double-null type symmetric divertor subdivided into 32 toroidal modules for the vertical maintenance. A detached divertor scenario with impurity seeding was considered as the primary approach for the power exhaust to reduce the peak heat flux lower than the engineering limit of 10 MW/m 2 . The power exhaust performance at the scrap off layer was estimated by using UEDGE-2D code, a two-dimensional fluid transport code for collisional edge plasma and neutral species like N, Ne, and Ar. Particle and heat flux on inboard and outboard divertor targets were calculated for the detached cased depending on parameters such as the impurity seeding rate, pumping rate, and the pedestal density. On the other hand, a magnetic solution like X-divertor, snowflake divertor, and super X-divertor to expand the plasma wet area was considered for K-DEMO since the detached divertor increasing a radiation fraction by impurity seeding might be able to be unstable. However, the extremely high current of poloidal coils was required more than the engineering limit, 20 MA, to form magnetic field lines for the alternative divertors. Based on the physical calculation of the edge plasma, engineering analyses were carried out to find out the thermal and structural reliability. The thermo-hydraulic analysis confirmed thermal stability, whether all comprising materials are operating within their allowable temperature windows when the case of the peak heat flux is set to 10 MW/m 2 on the outboard divertor target. The response surface optimization method derived two optimal design candidates employing two kinds of heat sink materials, respectively: the reduced activation ferritic martensitic (RAFM) steel and CuCrZr alloy. Additionally, the drawbacks and merits of the two materials were definite. The optimal design with applying RAFM steel was vulnerable to withstand thermal and mechanical loads since low thermal conductivity caused too thin thickness of the heat sink. On the other hand, the CuCrZr alloy has critical drawbacks in terms of activation and radioactive waste despite its high thermal conductivity. Meanwhile, preliminary electromagnetic (EM) analysis was carried out to estimate the EM loads caused by the abnormal behaviors of plasma since EM loads are one of the most critical external loads for designing a DEMO divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
DOE OSTI2021

Direct numerical simulation of single-species and binary-species boundary layers at high pressure

Direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate flow physics in a binary-species boundary layer at high pressure, with special attention to mass diffusion. The working fluids are nitrogen in the single-species cases, and a mixture of nitrogen and methane in the binary-species cases. An investigation of mean profiles shows that velocity and temperature profiles have steep gradients near the wall, whereas the mass fraction profiles do not have a steep gradient. This result indicates that the similarity law between velocity, temperature and mass fraction does not hold true in the binary-species boundary layer. The comparison of turbulent fluctuations shows that the qualitative characteristics of velocity and temperature are similar to each other. In contrast, profiles of mass fraction fluctuations are largely different from those of velocity and temperature. This result indicates that mass diffusion is not similar to momentum and thermal diffusions in binary-species boundary layers.

Toki, Takahiko↗
DOE OSTI2022

The timbre of Hawking gravitons: an effective description of energy transport from holography

Planar black holes in AdS, which are holographically dual to compressible relativistic fluids, have a long-lived phonon mode that captures the physics of attenuated sound propagation and transports energy in the plasma. We describe the open effective field theory of this fluctuating phonon degree of freedom. The dynamics of the phonon is encoded in a single scalar field whose gravitational coupling has non-trivial spatial momentum dependence. This description fits neatly into the paradigm of classifying gravitational modes by their Markovianity index, depending on whether they are long-lived. The sound scalar is a non-Markovian field with index 3 - d for a d-dimensional fluid. We reproduce (and extend) the dispersion relation of the holographic sound mode to quartic order in derivatives, constructing in the process the effective field theory governing its attenuated dynamics and associated stochastic fluctuations. We also remark on the presence of additional spatially homogeneous zero modes in the gravitational problem, which remain disconnected from the phonon Goldstone mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗