Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 397 records · Page 22

PBE-GGA predicts the B8↔B2 phase boundary of FeO at Earth’s core conditions

FeO is a crucial component of the Earth’s core, and its thermodynamic properties are essential to developing more accurate core models. It is also a notorious correlated insulator in the NaCl-type (B1) phase at ambient conditions. It undergoes two polymorphic transitions at 300 K before it becomes metallic in the NiAs-type (B8) structure at ~100 GPa. Although its phase diagram is not fully mapped, it is well established that the B8 phase transforms to the CsCl-type (B2) phase at core pressures and temperatures. Here, we report a successful ab initio calculation of the B8↔B2 phase boundary in FeO at Earth’s core pressures. We show that fully anharmonic free energies computed with the Perdew–Burke–Ernzerhof-generalized gradient approximation coupled with thermal electronic excitations reproduce the experimental phase boundary within uncertainties at P > 255 GPa, including the largely negative Clapeyron slope of –52 MPa/K. This study validates the applicability of a standard density functional theory functional to FeO under Earth’s core conditions and demonstrates the theoretical framework that enables complex predictive studies of this region.

58 GEOSCIENCES↗

Electron Temperature Gradient Driven Transport Model for Tokamak Plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

ETG↗

Efforts to stabilize composite localization elements

Localization finite elements seek to provide a robust framework for modeling ductile failure. They utilize the same constitutive model as the bulk material through the introduction of a length scale in a specialized deformation gradient that regularizes displacement discontinuity. Similar to many other elements, localization elements exhibit locking and associated pressure oscillations under incompressible plastic flow, which is a critical issue when attempting to model pressure-driven damage evolution. These issues can be drastically improved through what are essentially reduced integration techniques for the Jacobian and pressure, but there seem to be pressure-related instabilities that persist and are specific to localization elements. This memo summarizes recent efforts to mitigate and understand this problem, mostly for the 12-node composite wedge localization element in particular. At this point, it remains unclear whether the pressure fields within any localization element can be sufficiently stabilized in order to properly model failures that include softening or fracture.

36 MATERIALS SCIENCE↗

Fiber Optic Sensing Technologies for Battery Management Systems and Energy Storage Applications

Applications of fiber optic sensors to battery monitoring have been increasing due to the growing need of enhanced battery management systems with accurate state estimations. The goal of this review is to discuss the advancements enabling the practical implementation of battery internal parameter measurements including local temperature, strain, pressure, and refractive index for general operation, as well as the external measurements such as temperature gradients and vent gas sensing for thermal runaway imminent detection. A reasonable matching is discussed between fiber optic sensors of different range capabilities with battery systems of three levels of scales, namely electric vehicle and heavy-duty electric truck battery packs, and grid-scale battery systems. The advantages of fiber optic sensors over electrical sensors are discussed, while electrochemical stability issues of fiber-implanted batteries are critically assessed. This review also includes the estimated sensing system costs for typical fiber optic sensors and identifies the high interrogation cost as one of the limitations in their practical deployment into batteries. Finally, future perspectives are considered in the implementation of fiber optics into high-value battery applications such as grid-scale energy storage fault detection and prediction systems.

25 ENERGY STORAGE↗

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗

Expedition UT-GOM2-2 Site H

Pressure and conventional cores were collected at Site H of the Walker Ridge Protracted Area Block 313 in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico) during the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2). Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were recovered. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of about 200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. These ooze intervals also correspond to lighter sediment color, increased Ca content based on X-ray florescence (XRF) core scanning, and increased calcareous nannofossil abundance. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years), with a pronounced increase in sedimentation rate with depth. Below 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C1 /C2 ) and the methane to ethane plus propane (C1 /(C2 +C3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ13C isotopic signature of methane ranges between -69.9 and -78.5 ‰ relative to the Vienna Pee Dee Belemnite (VPDB) standard. Pressure core recovery in sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One pressure core was degassed and the average hydrate saturation in the core was determined to be 24%. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

Enhanced pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal high-βp scenario for the National Spherical Torus Experiment Upgrade (NSTX-U) and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (βN > 5), and significant bootstrap fraction (fBS > 0.6) at Ip/BT = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. The EP H-mode was most often observed as a transition following a large edge-localized mode in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with an elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the kinetic ballooning mode and kink-peeling mode, can be responsible for the increased particle transport in EP H-mode.

Battaglia, D. J. (ORCID:0000000188979740)↗

Operation above the Greenwald density limit in high performance DIII-D negative triangularity discharges

The density limit in strongly-shaped negative triangularity (NT) discharges is studied experimentally in the DIII-D tokamak. Record-high Greenwald fractions f G are obtained, using gas puff injection only, with values up to near 2, where f G is defined as the ratio of the line-averaged density over n G = I p /(π a 2 ), with I p [MA] the plasma current and a[m] the plasma minor radius. A clear higher operational limit with higher auxiliary power is also demonstrated, with the ohmic density limit about two times lower than with additional neutral beam injection heating. The evolution of the electron density, temperature and pressure profiles are analyzed as well. The core density can be up to twice the Greenwald density and keeps increasing, while the value at the separatrix remains essentially constant and slightly below n G . The edge temperature gradient collapses to near zero and NT plasmas are shown to be resilient to such profiles in terms of disruptivity. We also present the time evolution of the inverse electron pressure scale length with the value at the last closed flux surface (LCFS) decreasing below the value at the normalized radius 0.9 near the density limit, demonstrating the clear drop of confinement starting from the edge. This inverse scale length ‘collapse’ at the LCFS also defines well the characteristic behavior of the kinetic profiles approaching a density limit.

density limit↗

Developing a platform for Fresnel diffractive radiography with 1 μ m spatial resolution at the National Ignition Facility

An x-ray Fresnel diffractive radiography platform was designed for use at the National Ignition Facility. It will enable measurements of micron-scale changes in the density gradients across an interface between isochorically heated warm dense matter materials, the evolution of which is driven primarily through thermal conductivity and mutual diffusion. Here we use 4.75 keV Ti K-shell x-ray emission to heat a 1000 μm diameter plastic cylinder, with a central 30 μm diameter channel filled with liquid D 2 , up to 8 eV. This leads to a cylindrical implosion of the liquid D 2 column, compressing it to ~2.3 g/cm 3 . After pressure equilibration, the location of the D 2 /plastic interface remains steady for several nanoseconds, which enables us to track density gradient changes across the material interface with high precision. For radiography, we use Cu He- α x rays at 8.3 keV. Using a slit aperture of only 1 μm width increases the spatial coherence of the source, giving rise to significant diffraction features in the radiography signal, in addition to the refraction enhancement, which further increases its sensitivity to density scale length changes at the D 2 /plastic interface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗

Electron temperature gradient driven transport model for tokamak plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermoacoustic Soret separation

In a fluid mixture in a channel with an axial time-averaged temperature gradient, high-amplitude oscillating flow can greatly increase the axial flux of thermal diffusion (Soret) separation of the components of the mixture. The enhancement occurs when the oscillating lateral temperature gradient greatly exceeds the axial gradient, causing a large oscillating concentration that can be favorably time-phased with the oscillating flow. This process can occur even with a negligible pressure oscillation or with a negligible temperature response to pressure, as is the case in most liquid solutions. The thermal boundary condition imposed by realistic solids on thermoacoustic liquids is imperfect, adding mathematical complications that are absent for typical gases, for which the solid surface is temporally isothermal. Compared with gas mixtures, the high Lewis number in typical liquid solutions reduces the separation flux associated with the time-averaged temperature gradient, but it also reduces the remixing associated with the time-averaged mole-fraction gradient. For large enough channels, the second-law separation efficiency is only slightly reduced from that of steady liquid Soret separation.

36 MATERIALS SCIENCE↗

On the feasibility of using physics-informed machine learning for underground reservoir pressure management

In this work, we evaluate the feasibility of using physics-informed machine learning (PIML) for underground energy-related pressure management. To this end, we develop a PIML framework to manage underground reservoir pressures by training neural networks to determine fluid extraction rates for dedicated extraction wells during fluid injection operations given a range of reservoir conditions (e.g., transmissivity and storativity). We implement an automatically-differentiable analytical physics model of fluid flow in porous media within the PIML framework as a proxy for more complicated models. This allows us to execute a sufficient number of training scenarios to fully evaluate the feasibility of using PIML to support pressure management activities. We quantify the number of physics-model parameters required for automatic differentiation to become more efficient than finite-difference gradient calculations. We use a simple scenario with a single injector, extractor, and critical location for our feasibility analysis. We evaluate the effect of the size of the training dataset (i.e., the number of reservoir condition samples) on the accuracy and efficiency of the PIML framework. For an equivalent number of model evaluations, the larger training dataset took less time to train and produced a neural network that was able to more accurately manage reservoir pressures. We also evaluate the effect of the training dataset batch size (i.e., number of reservoir condition samples used to update the neural network coefficients during training; i.e., how the training dataset is partitioned). While training ran faster with larger batch sizes, they produced neural networks that managed pressures less accurately. We demonstrate the approach on a more complex scenario involving 10 injectors, 10 extractors, and 4 critical locations (a relatively high well density of 20 wells/km2). We provide the number of forward and adjoint model evaluations required in each case as an indication of the feasibility of using PIML for pressure management when more complicated physics models with longer execution times are used.

54 ENVIRONMENTAL SCIENCES↗

Coherent fluctuations in a partially magnetized small magnetron discharge operating on molecular nitrogen

We report experimental data on rotating plasma disturbances in a quasi-steady direct-current magnetron discharge operating with molecular nitrogen at pressures of 27, 40, 67, and 133 Pa. Experiments with a segmented anode measuring the local current flow reveal up to 12 coherent structures (“spokes”) propagating in the +E×B and −E×B azimuthal directions. The current–voltage characteristics are studied and found to follow a negative resistance profile at the lower discharge currents. The existence of coherent structures and the number of spokes observed depend on the discharge current and operating pressure. Lower pressures are more favorable to the retrograde −E×B spoke rotation, which is observed in the 27, 40, and 67 Pa discharges. The results are interpreted with a gradient-driven drift wave model, and using the dispersion relation, we show that the propagation direction depends on the orientation of the local electric field within the plasma. Using an electron fluid model, we show how experimental conditions can lead to field reversals that generate −E×B spoke rotation in the lower pressure discharges.

Physics↗

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Preliminary Cryogenic Cold Test Results of the First 9-Cell LSF Shape Cavity

Following successful prototyping and testing of single- & 5-cell LSF shape cavities *, **, the first 9-cell LSF shape cavity LSF9-1 was successfully constructed using an innovative process at JLab with the in-house facilities. The cavity was then shipped to KEK for post-fabrication mechanical adjustment and ILC TDR style treatment and surface processing. Cold testing was carried out at the JLab VTA facility, instrumented with a suite of Kyoto instruments. Favorable values for the bath pressure detuning sensitivity and Lorentz force detuning coefficient were experimentally measured, validating the design improvement in cell stiffeners. Pass-band measurements indicate 4 out of 9 cells reaching gradient capability of > 45 MV/m, including 2 cells reaching 51 MV/m. Cornell OST detectors identified the cell and location responsible for the current hard quench limit. Multipacting-like barriers observed in end cells are investigated both analytically and numerically. The cavity was shipped to FNAL and received a light EP at the joint ANL/FNAL facility for further cold testing at Jlab. Two new 9-cell LSF cavities are being constructed including one made of large-grain niobium material.

43 PARTICLE ACCELERATORS↗

Focusing Ions at Atmospheric Pressure using Nonlinear DC Voltage Sequences Applied to a Stacked Ring Ion Guide

Many modern ion mobility (IM) and mass spectrometers (MS) operate under low pressure (= 10 Torr) and employ high voltage radiofrequencies (RF) to provide ion confinement. Unfortunately, RF effectiveness drastically decreases as pressure increases, and few techniques for focusing ions at elevated pressures exist. Here we demonstrate a new approach for focusing ions at atmospheric pressure (AP) by applying nonlinear DC voltage sequences following quadratic and power (exponential) functions to a stacked ring ion guide. We used ion trajectory simulations to rigorously explore how ions react to nonlinear electric fields and validate the simulations with a set of ion current measurements performed at AP. Ion trajectory simulations show that ions initially defocus near the entrance of the device but then become intensely focused as they travel through the device. Contour plots for both nonlinear voltage sequences show electric field lines that increasingly curve inwards as a function of distance, resulting in spatial ion focusing. Experimental ion current and spot size measurements were performed at AP using a 10-cm stacked ring ion guide and a segmented Faraday cup detector. Quadratic sequences produced ~5% smaller spot sizes (~22.8 mm) and ~25% higher ion current compared to a linear voltage sequence (~24.0 mm). Alternatively, power sequences produced ~64% smaller spot sizes (~8.7 mm), albeit with ~10x lower ion current. However, both nonlinear voltage sequences produced similar ion currents at the center of the Faraday cup detector, indicating that higher ion densities are achieved when using nonlinear voltage gradients. These results demonstrate a new way to focus ions at AP, and the capabilities demonstrated here provide fundamental insights on how to keep ions inside analytical devices at elevated pressures without RF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measuring a Low Horizontal Hydraulic Gradient in a High Transmissivity Aquifer

Abstract High transmissivity aquifers typically have low hydraulic gradients (i.e., a flat water table). Measuring low gradients using water levels can be problematic because measurement error may be greater than the true difference in water levels (i.e., a low signal‐to‐noise ratio). In this study, the feasibility of measuring a hydraulic gradient in the range of 10 −6 to 10 −5 m/m was demonstrated. The study was performed at a site where the depth to water from land surface ranged from 40.1 to 94.2 m and the aquifer transmissivity was estimated at 41,300 m 2 /d (hydraulic conductivity of 18,800 m/d). The goals of the study were to reduce measurement error as much as practicable and assess the importance of factors affecting water level measurement accuracy. Well verticality was the largest source of error (0.000 to 0.168 m; median of 0.014 m), and geodetic survey of casing elevations was the next most important source of error (0.002 to 0.013 m; median of 0.005 m). Variability due to barometric pressure fluctuations was not an important factor at the site. Hydraulic heads were measured to an accuracy of ±0.0065 m, and the average hydraulic gradient was estimated to be 8.0 × 10 −6 (±0.9 × 10 −6 ) m/m. The improvement in accuracy allowed for two reversals in the groundwater flow direction to be identified, after which the gradient averaged 2.5 × 10 −5 (±0.4 × 10 −5 ) m/m. This study showed it is possible to sufficiently control sources of error to measure hydraulic gradients in the 10 −6 to 10 −5 m/m range.

McDonald, John P.↗