Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pressure Gradient”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Relevance of E × B drifts for particle and heat transport in divertors

Radial electric fields up to ~4 kV m –1 are observed in the boundary between the private flux region (PFR) and the scrape-off layer (SOL) driving E × B drifts between the inner and outer targets at speeds up to 2.8 km s –1 in the Tokamak à configuration variable divertor. The resulting E × B fluxes, located in a narrow region ($\Delta {\rho _\Psi } < 0.012$ in normalized radius or $\Delta $ R – R sep <4 mm mapped to the outer midplane) are equivalent to around 20% of the total heat and particle flux to the divertor targets (inner + outer). At the peak E r , the E × B poloidal transport is equivalent to parallel flows with M ∥ ~ 3. In the snowflake divertor with a second X-point in the outer SOL, the drifts in the PFR-SOL boundary were equivalent to around 30% of the total heat and particle flux to the divertor targets and cover a region ~50% wider than in the single null ($\Delta {\rho _\Psi }$ ~ 0.018, $\Delta $R – R sep ~ 6 mm). Furthermore, the location of the PFR-SOL boundary drift shifts radially in the E ∥ × B direction when reversing the toroidal field direction. Peaks in density and electron pressure have been identified near the primary X-point along with large gradients in density, temperature, and potential, the latter resulting in a local electric field ~2.7 kV m –1 which drives a drift (1.9 km s –1 ) upwards towards the closed flux surfaces. Floating potential (V f ) magnitudes up to 75 V (~2 kT e ) were measured, indicating that V f and parallel currents should not be neglected when estimating plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Isotope impact on Alfvén eigenmodes and fast ion transport in DIII-D

Abstract Measurements of beam driven Alfvén Eigenmode (AE) activity in matched deuterium (D) and hydrogen (H) DIII-D plasmas show a dramatic difference in unstable mode activity and fast ion transport for a given injected beam power. The dependence of the unstable AE spectrum in reversed magnetic shear plasmas on beam and thermal species is investigated in the current ramp by varying beam power in a sequence of discharges for fixed thermal and beam species at fixed density. In general, a spectrum of Reversed Shear Alfvén Eigenmodes (RSAEs) and Toroidal Alfvén Eigenmodes (TAEs) are driven unstable with sub-Alfvénic D beam injection while primarily only RSAEs are driven unstable for the H beam cases investigated. Further, for a given beam power, the driven AE amplitude is always reduced with H beams relative to D and for H thermal plasma relative to pure D or mixed D/H plasmas. Estimates of the fast ion stored energy combined with modeling using the hybrid kinetic-MHD code MEGA indicate that the dominant mechanism contributing to the difference between H and D beam drive is the faster classical slowing down of H beam ions relative to D and the resultant lower beam ion pressure. Calculations of the AE induced stored energy deficits using the reduced critical gradient model TGLFEP show quantitative agreement with the observed dependencies on injected power, isotope and minimum safety factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Prospects for the expansion of standing wave ambient pressure photoemission spectroscopy to reactions at elevated temperatures

Standing wave ambient pressure photoemission spectroscopy (SWAPPS) is a promising method to investigate chemical and potential gradients across solid-vapor and solid-liquid interfaces under close-to-realistic environmental conditions, far away from high vacuum. Until now, these investigations have been performed only near room temperature, but for a wide range of interfacial processes, chief among them being heterogeneous catalysis, measurements at elevated temperatures are required. One concern in these investigations is the temperature stability of the multilayer mirrors, which generate the standing wave field. At elevated temperatures, degradation of the multilayer mirror due to, for example, interdiffusion between the adjacent layers, decreases the modulation of the standing wave field, thus rendering SWAPPS experiments much harder to perform. Here, we show that multilayer mirrors consisting of alternate B4C and W layers are stable at temperatures exceeding 600 °C and are, thus, promising candidates for future studies of surface and subsurface species in heterogeneous catalytic reactions using SWAPPS.

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

In-situ thermomechanical measurements during burst testing of Zr claddings using digital image correlation and infrared thermography techniques

Burst testing has long been a relied-on method for assessing fuel cladding materials in light-water reactoraccident scenarios. However, burst testing historically provided minimal in-situ observations beyond a few thermocouple readings, pressure from transducers, and post-test deformation measurements. Acquisition of full cladding thermomechanical data throughout a burst test would be invaluable to fuel system qualification, material model development and fuel performance code validation, such as the BISON fuel performance code. Recently, methodologies for simultaneously applying 2-dimensional digital image correlation and infrared thermography techniques to quantify the relationship between temperature, internal cladding pressure, and deformation during burst testing were developed. This work details results from the application of these techniques to burst testing of Zircaloy-4 cladding segments subjected to loss-of-coolant accident transients in an air environment. Two types of tests were performed with: (a) cladding segments connected to a fixed pressure reservoir as well as (b) closed systems where internal pressure increased during transients. In-situ thermomechanical measurements including strain and thermal gradient data are presented, demonstrating the effectiveness of applying combined surface measurement techniques to burst testing

Bell, Sam↗
DOE OSTI2020

Technology Enabling Zero-EPZ Micro Modular Reactors: Fabrication of Zirconium Hydride with Controlled Hydrogen Loading

This report completes Milestone M2.2.2 - ZrH pellet fabricated with controlled hydride loading: ZrH1.5±0.1 single phase pellets fabricated with controlled and thermally stable hydride loading. Hydrogen content to be quantified and single-phase confirmed by X-ray diffraction. In this milestone, we report the successful fabrication of delta-phase zirconium hydride pellets using the ORNL bulk metal hydriding system. The successful deployment of zirconium hydride moderator in advanced reactors requires development of a consistent and affordable production pathway along with implementation of a hydrogen retention solution throughout the reactor life. Fabrication of delta-phase zirconium hydride is challenging since the absorption of a large amount of hydrogen into alpha-zirconium induces significant volume expansion and the present hydrogen concentration gradient results in cracking. A fully programmable hydriding system with continuous hydrogen partial pressure and flow control to facilitate processing of massive metal hydride has been developed at ORNL. In this report, the working principle of the hydriding system will be introduced. Characterization of the produced zirconium hydride includes X-ray powder diffraction to identify the present phases and X-ray Computed Tomography to visualize the internal microstructure (including cracks). The results indicate that single delta-phase zirconium hydride pellets with various sizes have been successfully produced.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mid-infrared Pulsed Upconversion Imaging in a Rotating Detonation Combustor

To meet the challenges associated with performing mid-InfraRed (IR) imaging in Rotating Detonation Combustors (RDCs), a novel pulsed mid-IR UpConversion Imaging (UCI) diagnostic has been implemented. UCI is an alternative to direct mid-IR detection that uses nonlinear optical frequency mixing to shift mid-IR wavelengths carrying a target image to shorter wavelengths that can be imaged with high-performance silicon-based CCD/CMOS cameras. This approach offers several favorable properties including high spectral selectivity, high temporal resolution, and superior low-light detectivity. A hydrogen-air research RDC was operated with carbon dioxide addition to allow pulsed UCI imaging of mid-IR luminosity within the combustion channel from spontaneous thermal emissions. The resulting measurements demonstrate high spatiotemporal resolution capable of imaging small structures near the supersonically propagating detonation wave front. The results show how this technique can be used to observe sharp gradients and millimeter-scale structures in the high-temperature, high-pressure zones RDC flow fields.

White, Logan W.↗
DOE OSTI2020

Genomics of sorghum local adaptation to a parasitic plant

Host–parasite coevolution can maintain high levels of genetic diversity in traits involved in species interactions. In many systems, host traits exploited by parasites are constrained by use in other functions, leading to complex selective pressures across space and time. Here, we study genome-wide variation in the staple crop Sorghum bicolor (L.) Moench and its association with the parasitic weed Striga hermonthica (Delile) Benth., a major constraint to food security in Africa. We hypothesize that geographic selection mosaics across gradients of parasite occurrence maintain genetic diversity in sorghum landrace resistance. Suggesting a role in local adaptation to parasite pressure, multiple independent loss-of-function alleles at sorghum LOW GERMINATION STIMULANT 1 (LGS1) are broadly distributed among African landraces and geographically associated with S. hermonthica occurrence. However, low frequency of these alleles within S. hermonthica-prone regions and their absence elsewhere implicate potential trade-offs restricting their fixation. LGS1 is thought to cause resistance by changing stereochemistry of strigolactones, hormones that control plant architecture and below-ground signaling to mycorrhizae and are required to stimulate parasite germination. Consistent with trade-offs, we find signatures of balancing selection surrounding LGS1 and other candidates from analysis of genome-wide associations with parasite distribution. Experiments with CRISPR–Cas9-edited sorghum further indicate that the benefit of LGS1-mediated resistance strongly depends on parasite genotype and abiotic environment and comes at the cost of reduced photosystem gene expression. Our study demonstrates long-term maintenance of diversity in host resistance genes across smallholder agroecosystems, providing a valuable comparison to both industrial farming systems and natural communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Low-pressure inductively coupled plasmas in hydrogen: impact of gas heating on the spatial distribution of atomic hydrogen and vibrationally excited states

Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10–20 Pa or 0.075–0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ten-moment fluid model for low-temperature magnetized plasmas

In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.

Kuldinow, Derek Amur (ORCID:0000000319730196)↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2021

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. Finally, the approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Technique for constructing high gradient insulators

A process for constructing a high-tensile strength, high-gradient insulator (HGI) may include stacking alternating layers of conductors and insulators, and vacuum pressure potting the stacked layers onto an insulating rod. The process may also include post machining the stacked layers to form a complete assembly of the HGI.

Watson, Scott Avery↗
DOE OSTI2022

Closed-loop pressure retarded osmosis draw solutions and their regeneration processes: A review

Pressure-Retarded Osmosis (PRO) is an osmotic process that has been used to harvest energy from salinity gradients using a semi permeable membrane. A comparison between open-loop PRO (OLPRO) and closed-loop PRO (CLPRO) was made regarding their performance and costs. In CLPRO, where the diluted draw solution is re-concentrated in the regeneration system to be reutilized in the process, has recently received an intensive focus as the most viable configuration for a standalone power plant. The choice of the PRO draw solution in CLPRO is crucial to garner a high osmotic pressure as the key for the feasibility of the process. Here, in this review, the draw solutions are critically evaluated in the literature in terms of energy output as well as the method of regeneration used to recirculate them. A set of practical criteria has been suggested to appraise the adequacy of the solution for CLPRO application. It was concluded that NH 3 – CO 2 theoretically can produce 170 W/m 2 of power density. Inorganic draw solutes such as NaCl can generate high power density up to 87 W/m 2 . Organic draw solutes with their remarkably low reverse salt flux (RSF) have promising potential for future application in PRO. Similarly, the regeneration systems of the diluted draw solutions have also been reviewed and discussed. How the energy consumption of the regeneration process affects the feasibility of CLPRO is explained. For the specific case of osmotic heat engines (OHEs), when the energy of the regeneration process is supplied by heat waste, the range of applicability of the heat waste in CLPRO in terms of efficiency is defined and compared to Organic Rankine Cycle (ORC). The results showed that CLPRO has better efficiency than ORC for temperatures T < 80 °C, which makes it a promising process or low-grade heat energy recovery. In addition, a PRO-RO hybrid system coupled with solar power can reduce the net specific energy consumption (SEC) to 0.39 kWh/m 3 . The conditions that regeneration processes should operate under to make PRO viable are discussed in the last section. Overall, the study indicates the key factors for optimizing the performance of CLPRO process.

42 ENGINEERING↗
DOE OSTI2022

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗
DOE OSTI2020

Setting the H-mode pedestal structure: variations of particle source location using gas puff and pellet fueling

Experiments in DIII-D show that a particle source location inside the top of the H-mode pedestal (pellets) maintains a higher pedestal pressure than an edge source (gas fueling) through a widening of the electron temperature pedestal with reduction of the temperature gradient. The effect of these two fueling schemes on the H-mode pedestal structure was examined in DIII-D by comparing controlled pellet-fueled and gas-fueled discharges across a fueling scan up to 40 torr 1 s –1 . High resolution electron profiles reveal that gas fueling lowers the pedestal pressure as the density profile shifts radially outwards and the separatrix density increases, while pellet fueling maintains a constant pedestal pressure. The neutral source locations from pellets and gas are determined with the PELLET and UEDGE codes, respectively, and quantify the particle source localization. Pellets provide significant ionization inside the pedestal top while gas puffing localizes ionization in the scrape-off-layer and pedestal foot, broadly consistent with the density profile structure influenced by the source. ELMs are observed to increase in frequency and reduce impurity content as fueling is increased. Stability analysis with ELITE shows that both conditions are near the type-I ELM corner of the peeling–ballooning stability diagram, which is altered significantly by the introduction of pellets. Since transport mechanisms are not observed to change substantially with particle source location, wider pedestals allow the pellet-fueled discharges to retain higher pedestal temperatures at similar pedestal densities. EPED1 is tested to capture the pedestal pressure, under-predicting the height with pellets and over-predicting the height with gas by ~15%. Overall, these results have important implications for future reactors where pellet fueling will be the primary particle source due to an opaque scrape-off-layer by showing that the neutral source location plays a role in setting the structure of the H-mode pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetics and Reactor Modeling of VFA Ketonization for Sustainable Aviation Fuel Production

Ketonization of volatile fatty acids (VFAs) produced via arrested methanogenesis of wet waste represents the only unit operation of wet waste upgrading to sustainable (net-zero or negative life cycle CO2 emissions) aviation fuel not currently in industrial practice. Ketone product yields of close to 100% have been obtained during gas-phase reactions of VFAs over oxide catalysts at the laboratory scale, but design of an industrial ketonization reactor requires understanding of the impacts of reactant and product partial pressures, deactivation, and heat and mass transport phenomena on observed ketonization rates. Our work leverages rigorous kinetic analysis of ketonization processes to inform reactor scale-up efforts through packed-bed reactor modeling. We first present results of a kinetic study of ketonization of a model VFA, hexanoic acid, over an industrial ZrO2 catalyst performed in a packed-bed microreactor in conditions free of significant heat or mass transfer gradients. Major findings of the analysis include: (i) hexanoic acid saturates all catalyst active sites at relatively low partial pressure (~10 kPa) and (ii) ketonization products 6-undecanone, water, and CO2 inhibit reaction rates. Kinetic data are used to fit a rate expression quantifying the functional dependence of ketonization rate on partial pressures of VFA reactants and ketone, water, and CO2 products. A packed-bed reactor model describing vapor-phase hexanoic acid ketonization over ZrO2 extrudate catalysts is developed based on the kinetic model. The effects of (i) bed- and pellet-scale mass- and heat-transfer limitations and (ii) axial pressure drop guide our development of recommendations for optimal sizing, temperature, and influent flow composition of an industrial-scale ketonization reactor. The quantitative understanding of VFA ketonization developed in this study represents an advance toward derisking the VFA ketonization step of wet waste upgrading to sustainable aviation fuel.

biojet↗
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 OSTI2020

Direct Measurement of Mass-Temperature Distributions in the Core of Implosions on the National Ignition Facility (Full Technical Final Report)

Inertial confinement fusion plasmas are characterized by several kilo-electron-volt (keV) electron temperatures and ion and electron densities approaching 100 to 1000 times solid. In these implosions which simulate the dynamics of our sun, spectroscopic dopants can lead to detailed exploration of the plasma gradients that are not possible in other places. This is both in the central core pressure and temperature as well as any “mix” resultant from ablation and hydrodynamic interactions. This work will show new Bayesian atomic physics models applied to ICF plasmas to give both the temperature and density of the internal hot spot to better describe the physical nature of these implosions along with trying to spatially resolve the mix coming from the ablator region. The mission impact is to demonstrate this effect and to lead to a more quantitative description of mix in inertial confinement plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗