Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Physical mechanism causing rapid changes in ultrarelativistic electron pitch angle distributions right after a shock arrival: Evaluation of an electron dropout event: Drift Shell Splitting on the Dayside

Three mechanisms have been proposed to explain relativistic electron flux depletions (dropouts) in the Earth's outer radiation belt during storm times: adiabatic expansion of electron drift shells due to a decrease in magnetic field strength, magnetopause shadowing and subsequent outward radial diffusion, and precipitation into the atmosphere (driven by EMIC wave scattering). Which mechanism predominates in causing electron dropouts commonly observed in the outer radiation belt is still debatable. In the present study, we evaluate the physical mechanism that may be primarily responsible for causing the sudden change in relativistic electron pitch angle distributions during a dropout event observed by Van Allen Probes during the main phase of the 27 February 2014 storm. During this event, the phase space density of ultrarelativistic (>1MeV) electrons was depleted by more than 1 order of magnitude over the entire radial extent of the outer radiation belt (3 < L* < 5) in less than 6 h after the passage of an interplanetary shock. We model the electron pitch angle distribution under a compressed magnetic field topology based on actual solar wind conditions. Although these ultrarelativistic electrons exhibit highly anisotropic (peaked in 90°), energy-dependent pitch angle distributions, which appear to be associated with the typical EMIC wave scattering, comparison of the modeled electron distribution to electron measurements indicates that drift shell splitting is responsible for this rapid change in electron pitch angle distributions. This further indicates that magnetopause loss is the predominant cause of the electron dropout right after the shock arrival.

58 GEOSCIENCES↗
DOE OSTI2019

Premixed flames subjected to extreme turbulence: Some questions and recent answers

It has been predicted that several changes will occur when premixed flames are subjected to the extreme levels of turbulence that can be found in practical combustors. Our report is a review of recent experimental and DNS results that have been obtained for the range of extreme turbulence, and it includes a discussion of cases that agree or disagree with predictions. “Extreme turbulence” is defined to correspond to a turbulent Reynolds number (Re Τ , based on integral scale) that exceeds 2800 or a turbulent Karlovitz number that exceeds 100, for reasons that are discussed in Section 2.1. Several data bases are described that include measurements made at Lund University, the University of Sydney, the University of Michigan and the U.S. Air Force Research Lab. The data bases also include DNS results from Sandia National Laboratory, the University of New South Wales, Newcastle University, the California Institute of Technology and the University of Cambridge. Several major observations are: (a) DNS now can be achieved for a realistic geometry (of the Lund University jet burner) even for extreme turbulence levels, (b) state relations (conditional mean profiles) from DNS and experiments do tend to agree with laminar profiles, at least for methane-air and hydrogen-air reactants that are not preheated, and (c) regime boundaries have been measured and they do not agree with predicted boundaries. Results indicate that the range of conditions for which flamelet models should be valid is larger than what was previously believed. Additional parameters have been shown to be important; for example, broken reactions occur if the “back-support” is insufficient due to the entrainment of cold gas into the product gas. Turbulent burning velocity measurements have been extended from the previous normalized turbulence levels (u’/S L ) of 24 up to a value of 163. Turbulent burning velocities no longer follow the trend predicted by Shchelkin but they tend to follow the trend predicted by Damköhler. The boundary where flamelet broadening begins was measured to occur at Re Taylor = 13.8, which corresponds to an integral scale Reynolds number (Re T ) of 2800. This measured regime boundary can be explained by the idea that flame structure is altered when the turbulent diffusivity at the Taylor scale exceeds a critical value, rather than the idea that changes occur when Kolmogorov eddies just fit inside a flamelet. A roadmap to extend DNS to complex chemistry and to higher Reynolds numbers is discussed. Exascale computers, machine learning, adaptive mesh refinement and embedded DNS show promise. Some advances are reviewed that have extended the use of line and planar PLIF and CARS laser diagnostics to studies that consider complex hydrocarbon fuels and harsh environments.

42 ENGINEERING↗
DOE OSTI2018

Experimental Investigation of Gas Flow and Hydrate Formation Within the Hydrate Stability Zone

Here, we form methane hydrate by injecting methane gas into a brine-saturated, coarse-grained sample under hydrate-stable thermodynamic conditions. Hydrate forms to a saturation of 11%, which is much lower than that predicted assuming three-phase (gas-hydrate-brine) thermodynamic equilibrium (67%). During hydrate formation, there are temporary flow blockages. We interpret that a hydrate skin forms a physical barrier at the gas-brine interface. The skin fails periodically when the pressure differential exceeds the skin strength. Once the skin is present, further hydrate formation is limited by the rate that methane can diffuse through the solid skin. This process produces distinct thermodynamic states on either side of the skin that allows gas to flow through the sample. This study illuminates how gas can be transported through the hydrate stability zone and thus provides a mechanism for the formation of concentrated hydrate deposits in sand reservoirs. It also illustrates that models that assume local equilibrium at the core-scale and larger may not capture the fundamental behaviors of these gas flow and hydrate formation processes.

58 GEOSCIENCES↗
DOE OSTI2015

Where is the cosmic-ray modulation boundary of the heliosphere?

The intensity of Galactic cosmic rays in the heliosphere is modulated by solar activities. The outer boundary where the solar modulation begins has always been a subject matter of debate in the cosmic-ray and heliophysics community. Various experimental methods and theoretical model calculations have been used to determine the boundary. Although the heliopause was always suspected to be the boundary, it is only until very recently after Voyager 1 had crossed the heliopause did we confirm that the boundary is indeed the heliopause. In this paper, we use a model simulation and detailed Voyager observation of cosmic rays at the heliopause crossing to show that the modulation boundary, in fact, is a fraction of an AU beyond the heliopause. Such a conclusion requires a very low turbulence level of the interstellar magnetic field in the outer heliosheath. According to the quasi-linear theory, a low level of turbulence should result in a very large diffusion coefficient parallel to the magnetic field and a very small perpendicular diffusion coefficient. For the first time, we are confident that Voyager 1 has obtained the truly pristine local interstellar cosmic-ray spectra down to the energies below 1 MeV. The cosmic-ray intensity is rapidly filtered by a thin layer of the interstellar magnetic field immediately outside of the heliopause. Its filtration amount depends on the conditions of magnetic field turbulence on the both sides of the heliopause, thus making it solar-cycle dependent as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2018

SELF-LIMITING TRANSIENT PULSE SIMULATION METHOD EXHIBITING TIME LAG PHENOMENON USING MAMMOTH

Coupling of multi-physics and multi-scale allows the exploration of phenomenon not previously studied with high fidelity. The heat diffusion time lag in TREAT-like fuel requires both multi-physics and multi-scale coupling to model the full transient behavior. The reactor physics code MAMMOTH, which has the ability to couple the heat conduction equation and the neutron diffusion equation, simulated a self-limiting transient pulse meeting these requirements. The MOOSE MultiApp system allowed the coupling of many micro-scale simulations with a single macro-scale simulation. The micro-scale geometry was a 0.0044 cm diameter UO2 fuel grain surrounded by a graphite moderator. The macro-scale geometry was a homogeneous fuel with the same isotope ratio of the micro-scale material surrounded by a graphite reflector. Two temperature-dependent two energy group cross-section libraries were computed using Serpent 2 for use in MAMMOTH. Both the homogeneous unit cell and heterogeneous unit cell cross-section library had infinite reactor boundary conditions. The macro-scale simulation had finite boundary conditions. By varying the feedback temperature, macro or micro, given to the cross-sections, the effect of the time lag was explored using a Picard coupling method. The time lag increases the peak power density by 9 % to 10 % and the energy deposited by 8 % to 9 %. The effect of differing grain and moderator temperature at the micro-scale was also explored showing differences in time lag behavior. Future research will improve the fidelity of the method through representative macro-scale geometries, temperature-dependent cross-sections, and greater energy group resolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
DOE OSTI2014

Generic Argillite/Shale Disposal Reference Case

Radioactive waste disposal in a deep subsurface repository hosted in clay/shale/argillite is a subject of widespread interest given the desirable isolation properties, geochemically reduced conditions, and widespread geologic occurrence of this rock type (Hansen 2010; Bianchi et al. 2013). Bianchi et al. (2013) provides a description of diffusion in a clay-hosted repository based on single-phase flow and full saturation using parametric data from documented studies in Europe (e.g., ANDRA 2005). The predominance of diffusive transport and sorption phenomena in this clay media are key attributes to impede radionuclide mobility making clay rock formations target sites for disposal of high-level radioactive waste. The reports by Hansen et al. (2010) and those from numerous studies in clay-hosted underground research laboratories (URLs) in Belgium, France and Switzerland outline the extensive scientific knowledge obtained to assess long-term clay/shale/argillite repository isolation performance of nuclear waste. In the past several years under the UFDC, various kinds of models have been developed for argillite repository to demonstrate the model capability, understand the spatial and temporal alteration of the repository, and evaluate different scenarios. These models include the coupled Thermal-Hydrological-Mechanical (THM) and Thermal-Hydrological-Mechanical-Chemical (THMC) models (e.g. Liu et al. 2013; Rutqvist et al. 2014a, Zheng et al. 2014a) that focus on THMC processes in the Engineered Barrier System (EBS) bentonite and argillite host hock, the large scale hydrogeologic model (Bianchi et al. 2014) that investigates the hydraulic connection between an emplacement drift and surrounding hydrogeological units, and Disposal Systems Evaluation Framework (DSEF) models (Greenberg et al. 2013) that evaluate thermal evolution in the host rock approximated as a thermal conduction process to facilitate the analysis of design options. However, the assumptions and the properties (parameters) used in these models are different, which not only make inter-model comparisons difficult, but also compromise the applicability of the lessons learned from one model to another model. The establishment of a reference case would therefore be helpful to set up a baseline for model development. A generic salt repository reference case was developed in Freeze et al. (2013) and the generic argillite repository reference case is presented in this report. The definition of a reference case requires the characterization of the waste inventory, waste form, waste package, repository layout, EBS backfill, host rock, and biosphere. This report mainly documents the processes in EBS bentonite and host rock that are potentially important for performance assessment and properties that are needed to describe these processes, with brief description other components such as waste inventory, waste form, waste package, repository layout, aquifer, and biosphere. A thorough description of the generic argillite repository reference case will be given in Jové Colon et al. (2014).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Oxidation of a model alkane aerosol by OH radical: the emergent nature of reactive uptake

An accurate description of the evolution of organic aerosol in the Earth’s atmosphere is essential for climate models. However, the complexity of multiphase chemical and physical transformations has been challenging to describe at the level required to predict aerosol lifetimes and changes in chemical composition. In this work a model is presented that reproduces experimental data for the early stages of oxidative aging of squalane aerosol by hydroxyl radical (OH), a process governed by reactive uptake of gas phase species onto the particle surface. Simulations coupling free radical reactions and Fickian diffusion are used to elucidate how the measured uptake coefficient reflects the elementary steps of sticking of OH to the aerosol as a result of a gas-surface collision, followed by very rapid abstraction of hydrogen and subsequent free radical reactions. It is found that the uptake coefficient is not equivalent to a sticking coefficient or an accommodation coefficient: it is an intrinsically emergent process that depends upon particle size, viscosity, and OH concentration. An expression is derived to examine how these factors control reactive uptake over a broad range of atmospheric and laboratory conditions, and is shown to be consistent with simulation results. Well-mixed, liquid behavior is found to depend on the reaction conditions in addition to the nature of the organic species in the aerosol particle.

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

Influence of skin surface roughness degree on energy characteristics of light scattered by a biological tissue

We present the results of modelling of photometric characteristics of light in soft tissues illuminated by a parallel beam along the normal to the surface, obtained with allowance for the skin roughness parameters and the angular structure of radiation approaching the surface from within the tissue. The depth structure of the fluence rate and the spectra of the diffuse reflection of light by the tissue in the interval of wavelengths 300 – 1000 nm are considered. We discuss the influence of the tilt angle variance of rough surface microelements and light refraction on the studied characteristics. It is shown that these factors lead to the reduction of the radiation flux only in the near-surface tissue layer and practically do not affect the depth of light penetration into the tissue. On the other hand, the degree of the surface roughness and the conditions of its illumination from within the tissue essentially affect the coefficient of diffuse reflection of light and lead to its considerable growth compared to the cases of a smooth interface and completely diffuse illumination, often considered to simplify the theoretical problem solution. The role of the roughness of skin surface is assessed in application to the solution of different direct and inverse problems of biomedical optics. (paper)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Partnership for Edge Physics Simulation

A major goal of our participation in the Edge Physics Simulation project has been to contribute to the understanding of the self-organization of tokamak turbulence fluctuations resulting in the formation of a staircase structure in the ion temperature. A second important goal is to demonstrate how small scale turbulence in plasmas self-organizes with dynamically driven quasi-stationary flow shear. These goals have been accomplished through the analyses of the statistical properties of XGC1 flux driven Gyrokinetic electrostatic ion temperature gradient (ITG) turbulence simulation data in which neutrals are included. The ITG turbulence data, and in particular fluctuation data, were obtained from a massively parallel flux-driven gyrokinetic full-f particle-in-cell simulation of a DIII-D like equilibrium. Below some the findings are summarized. It was observed that the emergence of staircase structure is related to the variations in the normalized temperature gradient length (R/LT) and the poloidal flow shear. Average turbulence intensity is found to be large in the vicinity of minima in R/LTi, where ITG growth is expected to be lower. The distributions of the occurrences of potential fluctuation are found to be Gaussian away from the staircase-step locations, but they are found to be non-Gaussian in the vicinity of staircase-step locations. The results of analytically derived expressions for the distribution of the occurrences of turbulence intensity and intensity flux were compared with the corresponding quantities computed using XGC1 simulation data and good agreement is found. The derived expressions predicts inward and outward propagation of turbulence intensity flux in an intermittent fashion. The outward propagation of turbulence intensity flux occurs at staircase-step locations and is related to the change in poloidal flow velocity shear and to the change in the ion temperature gradient. The standard deviation, skewness and kurtosis for turbulence quantities were computed and found to be large in the vicinity of the staircase-step structures. Large values of skewness and kurtosis can be explained by a temporary opening and closing of the structure which allows turbulence intensity events to propagate. The staircase patterns may reduce the ion heat transport and a manipulation of these patterns may be used to optimize heat transport in tokamaks. An additional objective of the research in support of the Edge Physics Simulation initiative has been to improve the understanding of scrape-off layer thermal transport. In planning experiments and designing future tokamaks, it is important to understand the physical effects that contribute to divertor heat-load fluxes. The research accomplished will contribute to developing new models for the scrape-off layer region. The XGC0 code was used to compute the heat fluxes and the heat-load width in the outer divertor plates of C-Mod and DIII-D tokamaks. It was observed that the width of the XGC0 neoclassical heat-load was approximately inversely proportional to the total plasma current. Anomalous transport in the H-mode pedestal region of five Alcator C-Mod discharges, representing a collisionality scan, was analyzed. The understanding of anomalous transport in the pedestal region is important for the development of a comprehensive model for the H-mode pedestal slope. It was found that the electron thermal anomalous diffusivities at the pedestal top increase with the electron collisionality. This dependence can point to the DRIBM as the modes that drive the anomalous transport in the plasma edge of highly collisional discharges. The effects of plasma shaping on the H-mode pedestal structure was also investigated. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities were discussed. For the discharges with higher elongation, it was found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continued to evolve to higher pedestal pressures and bootstrap currents until the peeling ballooning stability conditions were satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes were triggered at earlier times. The plasma elongation was found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency was large, and the H-mode pedestal evolves rapidly. It was found that the temperature of neutrals in the scrape-off-layer region can affect the dynamics of the H-mode pedestal buildup. However, the final pedestal profiles were nearly independent of the neutral temperature. The elongation and triangularity affected the pedestal widths of plasma density and electron temperature profiles differently. This study illustrated a new mechanism for controlling the pedestal bootstrap current and the pedestal stability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Effects of temperature and carbon source on the isotopic fractionations associated with O 2 respiration for 17 O/ 16 O and 18 O/ 16 O ratios in E. coli

18 O/ 16 O and 17 O/ 16 O ratios of atmospheric and dissolved oceanic O 2 are used as biogeochemical tracers of photosynthesis and respiration. Critical to this approach is a quantitative understanding of the isotopic fractionations associated with production, consumption, and transport of O 2 in the ocean both at the surface and at depth. We made measurements of isotopic fractionations associated with O 2 respiration by E. coli . Our study included wild-type strains and mutants with only a single respiratory O 2 reductase in their electron transport chains (either a heme-copper oxygen reductase or a bd oxygen reductase). We tested two common assumptions made in interpretations of O 2 isotope variations and in isotope-enabled models of the O 2 cycle: (i) laboratory-measured respiratory 18 O/ 16 O isotopic fractionation factors ( 18 α) of microorganisms are independent of environmental and experimental conditions including temperature, carbon source, and growth rate; And ( ii ) the respiratory ‘mass law’ exponent, θ, between 18 O/ 16 O and 17 O/ 16 O, 17α = ( 18 α) θ , is universal for aerobic respiration. Results demonstrated that experimental temperatures have an effect on both 18α and θ for aerobic respiration. Specifically, lowering temperatures from 37 to 15 °C decreased the absolute magnitude of 18 α by 0.0025 (2.5‰), and caused the mass law slope to decrease by 0.005. We propose a possible biochemical basis for these variations using a model of O 2 reduction that incorporates two isotopically discriminating steps: the reversible binding and unbinding of O 2 to a terminal reductase, and the irreversible reduction of that O 2 to water. In conclusion, we cast our results in a one-dimensional isopycnal reaction-advection-diffusion model, which demonstrates that enigmatic δ 18 O and Δ 17 O variations of dissolved O 2 in the dark ocean can be understood by invoking the observed temperature dependence of these isotope effects.

58 GEOSCIENCES↗
DOE OSTI2019

Hillock formation in co-deposited thin films of immiscible metal alloy systems

Hillocks are protruding sections of a surface that can influence the electronic and mechanical properties of a thin film. In monolithic films, these hillocks are formed due to thermally induced stress gradients during deposition. In this study, the surface morphology of co-sputtered immiscible Cu-X (X is a body centered cubic (BCC) group V or VI metal) thin films is characterized to elucidate the conditions that lead to hillock formation. Cu-Ta and Cu-Mo-Ag films were co-deposited with physical vapor deposition (PVD) magnetron sputtering at 25, 400, 600, and 800°C. A significant number of homogeneously distributed surface hillocks were observed only in the 600 and 800°C films. High-angle annular dark field (HAADF) cross-sectional imaging revealed significant Cu agglomerations underneath the protruding features. For the Cu-Ta films, the Cu was enveloped in a Cu-Ta nanocrystalline matrix. For Cu-Mo-Ag, the Cu was surrounded by Mo-Ag concentration modulations. While thermal stress gradients arise during deposition of immiscible metal films, biaxial stress calculations and literature cases reveal that they are not solely responsible for hillock formation. The observed morphologies align with a surface diffusion kinetic model that evaluates diffusion length of adatoms during deposition as a function of deposition temperature. Finally, this indicates a phase separation driving force paired with the constituent elements’ dissimilar mobilities at elevated deposition temperatures contribute to the presence of hillocks.

36 MATERIALS SCIENCE↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2015

Excitation of gravity waves by fingering convection, and the formation of compositional staircases in stellar interiors

Fingering convection (or thermohaline convection) is a weak yet important kind of mixing that occurs in stably stratified stellar radiation zones in the presence of an inverse mean molecular weight gradient. Brown et al. recently proposed a new model for mixing by fingering convection, which contains no free parameter and was found to fit the results of direct numerical simulations in almost all cases. Notably, however, they found that mixing was substantially enhanced above their predicted values in the few cases where large-scale gravity waves, followed by thermo-compositional layering, grew spontaneously from the fingering convection. This effect is well known in the oceanographic context and is attributed to the excitation of the so-called collective instability. In this work, we build on the results of Brown et al. and of Traxler et al. to determine the conditions under which the collective instability may be expected. We find that it is only relevant in stellar regions that have a relatively large Prandtl number (the ratio of the kinematic viscosity to the thermal diffusivity), O(10{sup −3}) or larger. This implies that the collective instability cannot occur in main-sequence stars, where the Prandtl number is always much smaller than this (except in the outer layers of surface convection zones, where fingering is irrelevant anyway). It could in principle be excited in regions of high electron degeneracy, during He core flash, or in the interiors of white dwarfs. We discuss the implications of our findings for these objects, from both a theoretical and an observational point of view.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Effects of surface processes on hydrogen outgassing from metal in desorption experiments

In this study, effects of surface processes on hydrogen outgassing during desorption experiments are investigated using a standard reaction-diffusion model describing hydrogen transport, retention and desorption in material. Three mechanisms for hydrogen migration and desorption on material surface are considered: hydrogen migration from material bulk to material surface (readsorption), hydrogen migration from material surface to material bulk (reabsorption) and hydrogen desorption from material surface by molecular recombination. Three hydrogen outgassing regimes are identified: i) recombination-limited outgassing when hydrogen recombination and desorption is fast compared to hydrogen reabsorption and slow compared to hydrogen transport from material bulk onto material surface ii) reabsorption-limited outgassing when hydrogen recombination and desorption is slow compared to hydrogen reabsorption and when the effective hydrogen recombination and desorption is slow compared to hydrogen transport from material bulk onto material surface iii) bulk-limited outgassing otherwise. Regimes of hydrogen outgassing from tungsten and beryllium are estimated for various experimental conditions. Analytical expressions of the outgassing flux are then derived for each outgassing regime, and used to characterize TDS spectra obtained in TDS experiments. It is shown that TDS spectra in recombination-limited regime are skewed toward high temperature, while TDS spectra in reabsorptionlimited and bulk-limited regimes are skewed toward low temperature. Furthermore, the temperature of the desorption peak in TDS spectra is shown to decrease in recombination-limited regime and to increase in reabsorption-limited and bulk-limited regimes as the total amount of hydrogen stored in material increases. Lastly, it is observed that the effective hydrogen recombination rate measured in reabsorption-limited permeations experiments may be not be reliably used to model and predict hydrogen retention and recycling from tungsten during plasma operations in tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗