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

Environmental controls on simulated deep moist convection initiation occurring during RELAMPAGO-CACTI

This study synthesizes the results of 13 high-resolution simulations of deep convective updrafts forming over idealized terrain using environments observed during the RELAMPAGO and CACTI field projects. Using composite soundings from multiple observed cases, and variations upon them, we explore the sensitivity of updraft properties (e.g., size, buoyancy, and vertical pressure gradient forces) to influences of environmental relative humidity, wind shear, and mesoscale orographic forcing that support or suppress deep convection initiation (CI). Emphasis is placed on differentiating physical processes affecting the development of updrafts (e.g., entrainment-driven dilution of updrafts) in environments typifying observed successful and null (i.e., no CI despite affirmative operational forecasts) CI events. Thermally-induced mesoscale orographic lift favors the production of deep updrafts originating from ~1–2-km-wide boundary layer thermals. Simulations without terrain forcing required much larger ( ~5-km-wide) thermals to yield precipitating convection. CI outcome was quite sensitive to environmental relative humidity; updrafts with increased buoyancy, depth, and intensity thrived in otherwise inhospitable environments by simply increasing the free tropospheric relative humidity. This implicates the entrainment of free-tropospheric air into updrafts as a prominent governor of CI, consistent with previous studies. Sensitivity of CI to the environmental wind is manifested by: 1) low-level flow affecting the strength and depth of mesoscale convergence along the terrain, and 2) clouds encountering updraft-suppressing pressure gradient forces while interacting with vertical wind shear in the free-troposphere. Among the ensemble of thermals occurring in each simulation, the widest deep updrafts in each simulation were the most sensitive to environmental influences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Hybrid RANS-LES of the Atmospheric Boundary Layer for Wind Farm Simulations: Preprint

Wind farm simulations often do not accurately represent wake-atmospheric boundary layer (ABL) interactions, blade boundary layer (BL) dynamics, and turbine-turbine interactions. In this work, we use Active Model Split (AMS), a new hybrid Reynolds-Averaged Navier Stokes (RANS)-large eddy simulation (LES) model, which is well suited to capture these effects because the model can (i) accurately simulate the ABL with the Coriolis effect, (ii) is accurate in adverse pressure gradients such as those near wind turbine blades, and (iii) has sufficiently low computational cost to simulate multiple turbines while resolving the blade BL. For simplicity and consistency we develop AMS to be used throughout the domain rather than in a zonal method. We implement our work in the massively parallel flow solver, Nalu-Wind, so that our model can access the compute resources needed for blade-resolved simulations of multiple wind turbines. To accomplish these aims, we modify the baseline AMS by changing the RANS contribution to SST k - omega with a length scale limiter, adding the Coriolis effect, and developing an appropriate wall treatment. We show that AMS of the ABL with the Coriolis effect matches LES reference results better than those obtained with RANS. We describe our plans to add buoyancy effects and wind turbines to our AMS simulations.

atmospheric boundary layer↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Investigation of mode activity in NBI-heated experiments of Wendelstein 7-X

The 2018 operation phase (OP 1.2b) of the stellarator Wendelstein 7-X (W7-X) included, for the first time, neutral beam injection (NBI) to heat the plasma. Since the injection geometry at W7-X is not parallel, this generates both passing and trapped fast particles. During longer phases of NBI injection, with the primary purpose to study the heating efficiency of this system, Alfvén eigenmodes (AEs) were observed by a number of diagnostics, including the phase contrast imaging (PCI) system, the magnetic pick-up coils (Mirnov coils), and the soft x-ray multi-camera tomography system (XMCTS). Alfvén eigenmodes are of great interest for future fusion reactors as it has been shown that the resonant interaction of fast ions with self-excited AEs can lead to enhanced transport of fast ions and potentially to energy losses. This is especially true for so-called gap-modes, Alfvén eigenmodes with frequencies in gaps of the continuous spectrum, since they lack continuum damping. These modes are commonly known to be excited by fast ions, but other destabilizing mechanisms, e.g. the electron-pressure gradient are also possible. In this article we present a first analysis of the experimentally observed frequencies from the theoretical side. The calculation of shear Alfvén wave continua for selected cases and the assignment of observed frequencies to the gaps of the continuous spectra are presented. Using the ideal-MHD code CKA (Könies A. 2007 10th IAEA TM on Energetic Particles in Magnetic Confinement System), we find gap modes that match the experimental measurements in terms of the observed frequencies. We emphasize the crucial roles played by the coupling of sound and Alfvén waves as well as of the Doppler shift arising as a consequence of the radial electric field in W7-X. We employ the perturbative gyrokinetic code CKA-EUTERPE (Feh´er 2013 Simulation of the interaction between Alfv´en waves and fast particles), using a slowing-down distribution function for the fast ions as calculated by the Monte-Carlo particle following code ASCOT (Hirvijoki et al 2014 Comput. Phys. Commun. 185 1310–21) to assess the fast-ion drive. Here we find that the fast-ion drive is insufficient to overcome the background-plasma damping. The fact that unstable modes were observed experimentally may point to problems with the modelling or indicate the existence of other destabilizing mechanisms, e.g. associated with the electron-pressure gradient (Windisch et al 2017 Plasma Phys. Control. Fusion 59 105002) that sensitively depend on the profiles of the background plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Generation of a Strong Parallel Electric Field and Embedded Electron Jet in the Exhaust of Moderate Guide Field Reconnection

Abstract Magnetospheric Multiscale observed an extended current layer concurrent with a strong parallel electric field. The current layer is embedded in the reconnection exhaust and is consistent with an extended magnetized jet expected in roughly symmetric moderate guide field reconnection. The strong parallel electric field observed at the jet balances a strong gradient in the parallel electron pressure caused by a transition between magnetized electrons with T e ‖ > T e ⊥ on one side of the jet and demagnetized, roughly isotropic electrons across the thin current layer. Simulation results show how this transition can occur and how the associated electron pressure gradients in Ohm's law balance the parallel electric field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Structural transitions of 4:1 methanol–ethanol mixture and silicone oil under high pressure

A 4:1 (volume ratio) methanol–ethanol (ME) mixture and silicone oil are two of the most widely used liquid pressure-transmitting media (PTM) in high-pressure studies. Their hydrostatic limits have been extensively studied using various methods; however, the evolution of the atomic structures associated with their emerging nonhydrostaticity remains unclear. Here, we monitor their structures as functions of pressure up to ∼30 GPa at room temperature using in situ high-pressure synchrotron x-ray diffraction (XRD), optical micro-Raman spectroscopy, and ruby fluorescence spectroscopy in a diamond anvil cell. No crystallization is observed for either PTM. The pressure dependence of the principal diffraction peak position and width indicates the existence of a glass transition in the 4:1 ME mixture at ∼12 GPa and in the silicone oil at ∼3 GPa, beyond which a pressure gradient emerges and grows quickly with pressure. There may be another liquid-to-liquid transition in the 4:1 ME mixture at ∼5 GPa and two more glass-to-glass transitions in the silicone oil at ∼10 GPa and ∼16 GPa. By contrast, Raman signals only show peak weakening and broadening for typical structural disordering, and Raman spectroscopy seems to be less sensitive than XRD in catching these structural transitions related to hydrostaticity variations in both PTM. These results uncover rich pressure-induced transitions in the two PTM and clarify their effects on hydrostaticity with direct structural evidence. The high-pressure XRD and Raman data on the two PTM obtained in this work could also be helpful in distinguishing between signals from samples and those from PTM in future high-pressure experiments.

Chen, Xiehang (ORCID:0000000168541506)↗
DOE OSTI2023

Extending the low-recycling, flat temperature profile regime in the lithium tokamak experiment-β (LTX-β) with ohmic and neutral beam heating

Abstract Recent experiments in the lithium tokamak experiment- β (LTX- β ) have extended the duration, performance, operating conditions, and diagnosis of the flat-temperature profile, low-recycling regime first observed in LTX. As expected, Li retains hydrogen and suppresses edge neutral cooling, allowing increased edge electron temperature, roughly equal to the core T e . Flat temperature profiles had been obtained transiently in LTX, as the plasma density decayed following the cessation of edge gas puffing. Careful control over the fueling in LTX- β has now been shown to sustain the flat T e profile and hot edge unique to the low-recycling regime for multiple confinement times in high performance discharges with decaying or steady density. With low density, the flat T e profile is also seen to extend into the scrape-off layer. Neutral beam heating is observed in target discharges with relatively flat electron temperature profiles ( T edge ∼ T core /2), though beam heating is stronger in discharges with higher fueling, higher density, and depressed edge T e . Beam heating produces additional peaking of the T e profile, without degradation of the energy confinement time. Neutral beam heating of target discharges with relatively flat electron temperature profiles similarly results in broad beam heated temperature profiles. Energy confinement in LTX- β generally compares favorably to ohmic and H-mode scalings, frequently exceeding them by factors of 2–4. New and improved diagnostics in LTX- β enable better characterization of this unique regime, including measurements of ion temperature and high field side Thomson scattering profiles. As an initial step toward characterizing turbulence with no T e gradient and roughly equal density and pressure gradient, core fluctuation spectra have been measured in peaked T e discharges using far-forward scattering and fluctuation reflectometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamic Formation of a Transient Jet from Arched Magnetized Laboratory Plasma

Abstract A laboratory plasma experiment was built to explore the eruptive behavior of arched magnetized plasmas with dimensionless parameters relevant to the Sun’s photosphere ( β ≈ 10 −3 , Lundquist number ≈10 4 , plasma radius/ion gyroradius ≈20, ion–neutral collision frequency ≫ion cyclotron frequency). Dynamic formation of a transient plasma jet was observed in the presence of the strapping magnetic field. The eruption leading to the jet is unintuitive because the arched plasma is both kink- and torus-stable. The jet structure erupts within a few Alfvén transit times from the formation of the arched plasma. Extensive measurements of plasma temperature, density, magnetic field, and flows are presented. In its early stages, the jet plasma flows away from the arch with supersonic speeds (Mach 1.5). This high-speed flow persists up to the resistive diffusion time in the arched plasma and is driven by large gradients in the magnetic and thermal pressures near the birthplace of jets. There are two distinct electric current channels within the jet, one consisting of outgoing electrons and another composed of electrons returning to the anode footpoint. Significant current density around the jet is a consequence of the diamagnetic current produced by a large thermal pressure gradient in the jet. Ion–neutral charge-exchange collisions provide an efficient mechanism to produce the cross-field current and control the dynamics of the complex current channels of the jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Blade Designs For Improved Multi-Phase Performance In sCO 2 Compressors; Part II - Optical Diagnostics In sCO 2 And Experimental Evaluation With Particle Image Velocimetry

This paper presents the second part of a study in which the leading-edge and suction surface of a compressor blade was modified to delay onset of phase change for sCO 2 compressors operating near the critical point. Using a first-of-its-kind apparatus for the measurement of sCO 2 flow fields, Particle Image Velocimetry (PIV) is used for local flow field measurements of two compressor blade geometries: the modified “biased-wedge,” and a conventional constant thickness blade. Utilizing the developed hardware, the feasibility of a simple, laser-based diagnostic for qualitatively measuring liquid phase regions, is also presented. The design of the optical diagnostics rig, a discussion of numerous challenges, and necessary considerations involved in performing optical-based measurements like PIV, in sCO 2 , are discussed. Velocity field measurements for the modified compressor profile show a much lower suction peak compared to a conventional blade. Furthermore, these results validate numerical results at the tested conditions, where the suction side profile of the biased wedge works to minimize the local pressure gradient.

14 SOLAR ENERGY↗
DOE OSTI2023

Flow friction in a wide rectangular duct with a beveled inlet: Informal Notes of Fundamental Fluid Physics Studies

Fundamental streamwise pressure distributions were measured at Idaho State University to address a primary objective of obtaining insight into the apparent evolution of air flow regimes at moderate Reynolds numbers in thin, wide rectangular ducts with converging inlets. Ducts of 3.030 and 4.077 mm spacings s (called three and four mm), having non-dimensional lengths of (L/s) ≈ 500.8 and 372.2 and widths (W/s) ≈ 114.6 and 83.2, respectively, were compared. Flow rates ranged from Reynolds numbers based on hydraulic diameter (Re Dh ) of about 980 to 10,000 and corresponding non-dimensional downstream pressure gradient parameters (-K p ) covered the range 0.038 to 0.0061, respectively.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Dusty Gas Model for Solid Oxide Fuel Cell Fuel Electrode

This model applies the Dusty Gas Model simulate multi-component species transport in SOFC (solid oxide fuel cell) anodes which considers the pressure gradient across the fuel electrode. This studyhas been verified with the analytical solution for different fuel electrode thicknesses and with literature values. The model was developed using the VoronoiFVM platform in Julia which is a built in implicit and semi implicit solver that integrates electrochemical behavior, microstructural effects, and transient analysis for accurate prediction of species transport under varying conditions.

dusty gas model (DGM)↗
DOE OSTI2020

Generic Post-Closure Safety Assessment for Deep Borehole Disposal of Heat-Generating Waste - 20451

The deep borehole disposal (DBD) concept for high-activity waste consists of drilling a large-diameter borehole into crystalline basement rock to a depth of about 5,000 m, placing waste packages in the lower portion of the borehole - the waste emplacement zone - and sealing and plugging the upper portion of the borehole with a combination of bentonite, cement plugs, and sand or crushed rock ballast. The waste emplacement zone is several times deeper than typical mined repositories and is well below the typical maximum depth of fresh groundwater resources. Recent work describes a generic reference case for deep borehole disposal of Cs and Sr capsules stored at the Hanford Site. Two scenarios are described: an undisturbed (nominal) scenario, and disturbed scenario, in which a waste package is stuck adjacent to a transmissive fracture zone in crystalline basement. Simulations of the processes affecting radionuclide release and transport in the two scenarios provide the basis for quantitative post-closure performance assessment (PA) of the reference case. This paper presents a probabilistic PA of the stuck waste package scenario, i.e. uncertainties affecting radionuclide mobilization and transport are propagated through 200 realizations of a predictive model to quantify uncertainties in predicted radionuclide concentrations within the borehole and the fracture zone. Simulations predict no radionuclide transport in the borehole to distances greater than 100 m above the stuck waste package. An assumed lateral pressure gradient results in predictions of limited radionuclide transport in the fracture zone at a lateral distance of 200 m from the stuck waste package. These results complement previous results for the nominal scenario and contribute to the development of a generic safety case for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Collinear dual-pulse laser optical breakdown and energy deposition

A non-equilibrium model for laser-generated plasmas is used to represent collinear dual nano-second-pulse interactions. In the specific case considered, the breakdown is initiated with an ultraviolet (UV) laser pulse at 266 nm, which pre-ionizes the gas, and it is followed by a near-infrared (NIR) pulse at 1064 nm, which deposits significant energy into the ionized mixture. The model is validated against corresponding experiments, and simulation results are interrogated to understand key features of the plasma-kernel dynamics and the post-discharge hydrodynamics. The hydrodynamics of the non-equilibrium plasma is governed by the two-temperature Navier–Stokes equations accounting for both multiphoton ionization and inverse bremsstrahlung. The interaction between the laser beam and the plasma is modeled based on the radiative transfer equation. The temporal and spatial offsets of the two pulses can generate ionization kernels with different topology and dynamics. It is shown that the UV pre-ionization pulse can tailor the plasma region, leading to a larger ionized volume than would occur for a single-pulse breakdown, and it increases the efficiency of the energy deposition for the following NIR discharge. Vorticity in the early post-discharge phase is generated via baroclinic torque from the misalignment of the radial gradient of density (sudden gas expansion) and the strong pressure gradient (initiated by energy deposition from the laser beam). It is predicted that the collinear dual pulse can be designed to specify the plasma kernel size, shape, and maximum temperature in the breakdown phase, and the initial post-breakdown vorticity and plasma-core decay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Emergence of a nocturnal low-level jet from a broad baroclinic zone

An analytical model is presented for the generation of a Blackadar-like nocturnal low-level jet in a broad baroclinic zone. The flow is forced from below (flat ground) by a surface buoyancy gradient and from above (free atmosphere) by a constant pressure gradient force. Diurnally-varying mixing coefficients are specified to increase abruptly at sunrise and decrease abruptly at sunset. With attention restricted to a surface buoyancy that varies linearly with a horizontal coordinate, the Boussinesq-approximated equations of motion, thermal energy, and mass conservation reduce to a system of one-dimensional equations that can be solved analytically. Sensitivity tests with southerly jets suggest that (i) stronger jets are associated with larger decreases of the eddy viscosity at sunset (as in Blackadar theory), (ii) the nighttime surface buoyancy gradient has little impact on jet strength, and (iii) for pure baroclinic forcing (no free-atmosphere geostrophic wind), the nighttime eddy diffusivity has little impact on jet strength, but the daytime eddy diffusivity is very important and has a larger impact than the daytime eddy viscosity. The model was applied to a jet that developed in fair weather conditions over the Great Plains from southern Texas to northern South Dakota on 1 May 2020. The ECMWF Reanalysis v5 (ERA5) for the afternoon prior to jet formation showed that a broad north-south-oriented baroclinic zone covered much of the region. The peak model-predicted winds were in good agreement with ERA5 winds and lidar data from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) central facility in north-central Oklahoma.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hard, transparent, sp3-containing 2D phase formed from few-layer graphene under compression

Despite several theoretically proposed two-dimensional (2D) diamond structures, experimental efforts to obtain such structures are in initial stage. Recent high-pressure experiments provided significant advancements in the field, however, expected properties of a 2D-like diamond such as sp 3 content, transparency and hardness, have not been observed together in a compressed graphene system. Here, we compress few-layer graphene samples on SiO 2 /Si substrate in water and provide experimental evidence for the formation of a quenchable hard, transparent, sp 3 -containing 2D phase. Our Raman spectroscopy data indicates phase transition and a surprisingly similar critical pressure for two-, five-layer graphene and graphite in the 4e6 GPa range, as evidenced by changes in several Raman features, combined with a lack of evidence of significant pressure gradients or local non-hydrostatic stress components of the pressure medium up to ≈ 8 GPa. The new phase is transparent and hard, as evidenced from indentation marks on the SiO 2 substrate, a material considerably harder than graphene systems. Furthermore, we report the lowest critical pressure (≈ 4 GPa) in graphite, which we attribute to the role of water in facilitating the phase transition. Overall, theoretical calculations and experimental data indicate a novel, surface to- bulk phase transition mechanism that gives hint of diamondene formation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Microbarograph - ESRL Hi-Res Microbarograph, Troutdale - Raw Data

High-precision barometers (Paroscientific 6000-16B-IS) are combined with Nishiyama-Bedard Quad Disk pressure probes, measuring pressure (mb) at the surface, nominally 2 m above ground level. Data are sampled at 20 Hz for potential studies of turbulence. The sensors' high accuracy makes them useful for determining horizontal pressure gradients and their relation to wind ramp events, as well as the temporal variability of pressure associated with mountain wakes and waves. **Note different ASCII file formats for Goldendale (z04) and Walla Walla (z09) sites.**

17 WIND ENERGY↗