Search NASA⌕ Search

SEARCH · Search NASA

Results for “pressure dependence”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetic structure, excitations and short-range order in honeycomb Na 2 Ni 2 TeO 6

Na 2 Ni 2 TeO 6 has a layered hexagonal structure with a honeycomb lattice constituted by Ni 2+ and a chiral charge distribution of Na + that resides between the Ni layers. In the present work, the antiferromagnetic transition temperature of Na 2 Ni 2 TeO 6 is confirmed at TN˜ 27 K, and further, it is found to be robust up to 8 T magnetic field and 1.2 GPa external pressure; and, without any frequency-dependence. Slight deviations from nominal Na-content (up to 5%) does not seem to influence the magnetic transition temperature, TN. Isothermal magnetization curves remain almost linear up to 13 T. Our analysis of neutron diffraction data shows that the magnetic structure of Na 2 Ni 2 TeO 6 is faithfully described by a model consisting of two phases described by the commensurate wave vectors vec kc, (0.5 0 0) and (0.5 0 0.5), with an additional short-range order component incorporated in to the latter phase. Consequently, a zig-zag long-range ordered magnetic phase of Ni 2+ results in the compound, mixed with a short-range ordered phase, which is supported by our specific heat data. Theoretical computations based on density functional theory (DFT) predict predominantly in-plane magnetic exchange interactions that conform to a J1-J2-J3 model with a strong J3 term. The computationally predicted parameters lead to a reliable estimate for TN and the experimentally observed zig-zag magnetic structure. A spin wave excitation in Na 2 Ni 2 TeO 6 at E ˜ 5 meV at T = 5 K is mapped out through inelastic neutron scattering experiments, which is reproduced by linear spin wave theory calculations using the J values from our computations. Our specific heat data and inelastic neutron scattering data strongly indicate the presence of short-range spin correlations, at T > TN, stemming from incipient antiferromagnetic clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Cloud Atlas: Unraveling the Vertical Cloud Structure with the Time-series Spectrophotometry of an Unusually Red Brown Dwarf

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352—an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR 8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the HST and the Spitzer time-series observations, as well as the time-averaged spectra of WISE 0047. Additionally, in the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/CH 4 chemistry, our models also reproduce the redder $J-{K}_{{\rm{s}}}$ color of WISE 0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Molecular structure and catalytic promotional effect of Mn on supported Na 2 WO 4 /SiO 2 catalysts for oxidative coupling of methane (OCM) reaction

The structure and promotional effect of Mn in supported Mn-Na 2 WO 4 /SiO 2 catalysts for the oxidative coupling of methane (OCM) reaction has been debated for a longtime in the literature. In the current investigation, with the aid of multiple in-situ characterization studies, we show that the freshly calcined supported 1.2Mn-5Na 2 WO 4 /SiO 2 catalyst possesses crystalline Na 2 WO 4 , Mn 2 O 3 and SiO 2 (cristobalite phase) along with surface MnO x and Na-WO x sites at low temperature and oxidizing environments. Under the OCM reaction environment (T>800°C), the crystalline Na 2 WO 4 phase melts and Mn 2 O 3 phase reduces. In contrast, the surface MnO x and Na-WO x sites exhibit excellent thermal and chemical stability. Exposure of the 1.2Mn-5Na 2 WO 4 /SiO 2 catalyst to the OCM reaction environment redisperses the molten Na 2 WO 4 phase on the SiO 2 support to form new surface WO x sites. Interestingly, the stable MnO x species interacts with both molten Na 2 WO 4 phase and surface Na-WO x sites during OCM reaction. Controlled transient kinetic experiments in TAP and detailed steady state OCM fixed-bed reaction studies reveal the role and promotional effect of Mn in the 1.2Mn-5Na 2 WO 4 /SiO 2 catalyst. The W-oxides (both molten Na 2 WO 4 and surface Na-WO x sites) are the active sites for the catalytic OCM reaction and the MnO x species only function as promoters. The promotion of MnO x strongly depends on the gas phase O 2 partial pressure and the MnO x species act as mediators for oxygen exchange between the gas phase molecular O 2 and catalyst lattice oxygen. In conclusion, the temperature dependent MnO x promotion reveals that the MnO x species selectively promote the molten Na 2 WO 4 phase at lower reaction temperature and the surface Na-WO x sites at higher temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Post‐Perovskite Transition in Fe‐ and Al‐Bearing Bridgmanite: Effects on Seismic Observables

Abstract The primary phase of the Earth’s lower mantle, (Al, Fe)‐bearing bridgmanite, transitions to the post‐perovskite (PPv) phase at Earth’s deep mantle conditions. Despite extensive experimental and ab initio investigations, there are still important aspects of this transformation that need clarification. Here, we address this transition in (Al 3+ , Fe 3+ )‐, (Al 3+ )‐, (Fe 2+ )‐, and (Fe 3+ )‐bearing bridgmanite using ab initio calculations and validate our results against experiments on similar compositions. Consistent with experiments, our results show that the onset transition pressure and the width of the two‐phase region depend distinctly on the chemical composition: (a) Fe 3+ ‐, Al 3+ ‐, or (Al 3+ , Fe 3+ )‐alloying increases the transition pressure, while Fe 2+ ‐alloying has the opposite effect; (b) in the absence of coexisting phases, the pressure‐depth range of the Pv‐PPv transition is likely too broad to cause a sharp D” discontinuity (<30 km); (c) the average Clapeyron slope of the two‐phase regions are consistent with previous measurements, calculations in MgSiO 3 , and inferences from seismic data. In addition, (d) we observe a softening of the bulk modulus in the two‐phase region. The consistency between our results and experiments gives us the confidence to proceed and examine this transition in aggregates with different compositions computationally, which will be fundamental for resolving the most likely chemical composition of the D" region by analyses of tomographic images.

Geochemistry & Geophysics↗
DOE OSTI2023

Transient Pressure Interference During CO2 Injection in Saline Aquifers

Abstract CO2 injection in subsurface geological formations (e.g. deep saline aquifers) causes pressure perturbations over a large area surrounding the injection well. Observation wells are widely considered in geologic CO2 storage (GCS) projects where the pressure perturbation induced by CO2 injection is measured. In this work, we use analytical and numerical modeling tools along with field data to examine the pressure behavior in GCS projects before and after CO2 arrival at an observation well. Prior to CO2 arrival, a baseline pressure trend is established which corresponds to single-phase brine flow across the observation well (approximated by Theis solution). Therefore, analysis of early-time pressure data is straightforward, provides the single-phase flow characteristics (mobility and storativity), and helps establishing a baseline pressure change that can be extended beyond the single-phase flow period at the observation well. Upon CO2 arrival, a departure from this baseline trend is expected. For the pressure to detect the CO2 arrival at an observation well, the departure from baseline pressure behavior must be significant and well above the background noise levels. We use existing analytical models to determine the strength of the expected pressure departure signal from the baseline trend upon CO2 arrival. The strength of the expected pressure departure is found to be directly proportional to the mobility ratio. Accordingly, we establish a criterion to determine whether the pressure at an observation well can detect the CO2 arrival. We present an analysis approach through application to synthetic and field data and show the characteristic pressure behavior before and after CO2 arrival. We show that while generally the pressure can be either above or below the expected baseline pressure trend, it would be likely above the baseline upon CO2 arrival. This is because the mobility ratio becomes less than unity after CO2 arrival. We show that depending on the reservoir characteristics, changes in the pressure trend may or may not be sufficient to detect the CO2 arrival.

Engineering↗
DOE OSTI2019

Pressure-driven valence increase and metallization in the Kondo insulator Ce 3 Bi 4 Pt 3

We report the results of high-pressure x-ray diffraction, x-ray absorption, and electrical transport measurements of the Kondo insulator (KI) Ce 3 Bi 4 Pt 3 up to 42 GPa, the highest pressure reached in the study of any Ce-based KI. Here, we observe a smooth decrease in volume and movement toward intermediate Ce valence with pressure, both of which point to increased electron correlations. Despite this, temperature-dependent resistance data show the suppression of the interaction-driven ambient pressure insulating ground state. We also discuss potential ramifications of these results for the predicted topological KI state.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Stable and metastable structures of tin (IV) oxide at high pressure

We have analyzed SnO 2 with a combination of synchrotron X-ray diffraction and X-ray absorption spectroscopy across a pressure range of 0 → 82.9 GPa with thermal annealing by a CO 2 laser allowing access to all of the known high-density polymorphs of SnO 2 , and here report their crystallographic information. The metastability of the post-rutile α -PbO 2 and PdF 2 structures in SnO 2 are investigated by experiment and PW-DFT simulations, revealing a complex energetic landscape and suggesting a significant dependence of the observed phases on the pressure–temperature pathway taken in experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ignition delay times of NH 3 /DME blends at high pressure and low DME fraction: RCM experiments and simulations

Here, autoignition delay times of ammonia/dimethyl ether (NH 3 /DME) mixtures were measured in a rapid compression machine with DME fractions of 0, 2 and 5 and 100% in the fuel. The measurements were performed at equivalence ratios $\varphi$= 0.5, 1.0 and 2.0 and pressures in the range 10–70 bar; depending on the fuel composition, the temperatures after compression varied from 610 K to 1180 K. Admixture of DME is seen to have a dramatic effect on the ignition delay time, effectively shifting the curves of ignition delay vs. temperature to lower temperatures, up to ~250 K compared to pure ammonia. Two-stage ignition is observed at $\varphi$= 1.0 and 2.0 with 2% and 5% DME in the fuel, despite the pressure being higher than that at which pure DME shows two-stage ignition. At $\varphi$= 0.5, a reproducible pre-ignition pressure rise is observed for both DME fractions, which is not observed in the pure fuel components. A novel NH 3 /DME mechanism was developed, including modifications in the NH 3 subset and addition of the NH 2 + CH 3 OCH 3 reaction, with rate coefficients calculated from ab initio theory. Simulations faithfully reproduce the observed pre- ignition pressure rise. While the mechanism also exhibits two-stage ignition for NH 3 /DME mixtures, both qualitative and quantitative improvement is recommended. The overall ignition delay times for ammonia/DME mixtures are predicted well, generally being within 50% of the experimental values, although reduced performance is observed for pure ammonia at $\varphi$= 2.0. Simulating the ignition process, we observe that the DME is oxidized much more rapidly than ammonia. Analysis of the mechanism indicates that this ‘early DME oxidation’ generates reactive species that initiate the oxidation of ammonia, which in turn begins heat release that raises the temperature and accelerates the oxidation process towards ignition. The reaction path analysis shows that the low-temperature chain-branching reactions of DME are important in the early oxidation of the fuel, while the sensitivity analysis indicates that several reactions in the oxidation of DME, including cross reactions between DME and NH 3 species presented here, are critical to ignition, even at fractions of 2% DME in the fuel.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Mode-specific pressure effects on the relaxation of an excited nitromethane molecule in an argon bath

The vibrational and rotational mode-specific relaxations of CH 3 NO 2 with 50 kcal/mol of initial internal energy in an argon bath is computed at 300 K at pressures of 10-400 atm. This work uses archived information from our previously published [J. Chem. Phys. 142, 014303 (2015)] molecular dynamics simulations and employs our previous published [J. Chem. Phys. 151, 034303 (2019)] method for projecting time-dependent Cartesian velocities onto normal mode eigenvectors. The computed relaxations cover three types of energies: vibrational, rotational, and Coriolis. In general, rotational and Coriolis relaxations in all modes are initially fast followed by an orders of magnitude slower relaxation. For all modes, that slower relaxation rate is approximately comparable to the vibrational relaxation rate. For all three types of energies, there are small-scale mode-to-mode variations. Of particular prominence is the exceptionally fast relaxation shared in common by the external rotation about the C-N axis, the internal hindered rotation of the CH 3 group relative to the NO 2 group, and the symmetric stretch of the CH 3 group.

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

Candidate structure for the H 2 -PRE phase of solid hydrogen

Experimental progress finally reached the metallic solid hydrogen phase, which was predicted by Wigner and Huntington over 80 years ago. However, the different structures in the phase diagram are still being debated due to the difficulty of diffraction experiments for high-pressured hydrogen. The determination of crystal structures under extreme condition is both of the basic condensed matter physics, and in planetary science: the behavior of giant gaseous planets (e.g., Jupiter and Saturn) strongly depends on the properties of inner high-pressured hydrogen. This work describes possible structures appearing under high pressures of 400~600 GPa. We applied a structural search using particle swarm optimization with density functional theory (DFT) to propose several candidate structures. For these structures, we performed fixed-node diffusion Monte Carlo simulations combined with DFT zero-point energy corrections to confirm their relative stability. We found P2 1 /c-8 as a promising candidate structure for the H 2 -PRE phase. P2 1 /c-8 is predicted the most stable at 400 and 500 GPa. P2 1 /c-8 reproduces qualitatively the IR spectrum peaks observed in the H 2 -PRE phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Polymorphism and Negative Linear Compressibility in Pyrazine-d4

Pyrazine (1,4-diazine) is a widely used linker in coordination polymers, molecular magnets, and metal–organic frameworks, yet its solid-state polymorphism is controlled by weak C–H···N interactions. Using variable-temperature and high-pressure neutron diffraction, a previously unrecognized low-temperature transition to phase IV is identified, showing that the form stabilized under pressure can also be accessed below ∼90 K at ambient pressure in the perdeuterated compound. The transition from phase III → IV involves a symmetry-breaking lattice distortion driven by the reorganization of C–H···N hydrogen bonds, accompanied by negative linear compressibility, which is primarily accommodated through the reduction in void-space within the structure. Structural analysis reveals that while the overall monoclinic framework of phase IV is preserved under both temperature and pressure routes, the evolution of hydrogen-bonded network is distinct depending on the pathway taken. These results establish pyrazine as a model molecular solid for studying pressure–temperature polymorphism and demonstrate how weak directional interactions can produce anomalous mechanical responses in organic crystals.

Ridley, Christopher J. [Oak Ridge National Laborat↗
DOE OSTI2024

Viscosity Measurements at High Pressures: A Critical Appraisal of Corrections to Stokes' Law

Abstract Fluids and melts in planetary interiors significantly influence geodynamic processes from volcanism to global‐scale differentiation. The roles of these geofluids depend on their viscosities ( η ). Constraining geofluid η at relevant pressures and temperatures relies on laboratory‐based measurements and is most widely done using Stokes' Law viscometry with falling spheres. Yet small sample chambers required by high‐pressure experiments introduce significant drag on the spheres. Several correction schemes are available for Stokes' Law but there is no consensus on the best scheme(s) for high‐pressure experiments. We completed high‐pressure experiments to test the effects of (a) the relative size of the sphere diameter to the chamber diameter and (b) the top and bottom of the chamber, that is, the ends, on the sphere velocities. We examined the influence of current correction schemes on the estimated viscosity using Monte Carlo simulations. We also compared previous viscometry work on various geofluids in different experimental setups/geometries. We find the common schemes for Stokes' Law produce statistically distinct values of η . When inertia of the sphere is negligible, the most appropriate scheme may be the Faxén correction for the chamber walls. Correction for drag due to the chamber ends depends on the precision in the sinking distance and may be ineffective with decreasing sphere size. Combining the wall and end corrections may overcorrect η . We also suggest the uncertainty in η is best captured by the correction rather than propagated errors from experimental parameters. We develop an overlying view of Stokes' Law viscometry at high pressures.

Geochemistry & Geophysics↗
DOE OSTI2025

Effects of equilibrium pressure on plasma response to RMPs in a spherical tokamak

This study presents a comprehensive analysis of the equilibrium pressure on the plasma response to resonant magnetic perturbations (RMPs) in the spherical tokamak (ST) MAST-U, employing both single-fluid and MHD-kinetic hybrid models (implemented via the MARS-F/K codes). As a key finding, the study identifies two different pressure-driven eigenmodes, exhibiting Sturmian property, that affect the Troyon no-wall limits for the onset of the n = 1 and n = 2 ( n is the toroidal mode number) ideal external kink instabilities as well as the corresponding plasma response to the applied RMP. With increasing equilibrium pressure, the plasma response to RMPs is significantly enhanced in the ST plasma, particularly in the high-pressure regime where kinetic effects strongly stabilize the external kink instability. The Troyon no-wall limit divides the plasma response into two regions: well below the limit, the response amplitudes and trends (versus pressure) are similar between the fluid and kinetic models; as the equilibrium pressure approaches the Troyon limit, the kinetic model predicts significant amplification of the RMP field, up to 30 times for cases considered. A relatively weak dependence of the optimal coil phasing on the equilibrium pressure is computed in this ST plasma, similar to the trend obtained for the conventional aspect ratio devices. These findings underscore the importance of incorporating kinetic effects in accurate prediction of the plasma response to RMPs in high-pressure ST tokamak plasmas and provide a theoretical basis for optimizing RMP-based control of the edge-localized modes in future ST devices.

RMP↗