Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

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 937 records · Page 52

Synthesis and Structural Study of Substituted Ternary Nitrides for Ammonia Production

Over the past few decades, inorganic nitride materials have grown in importance in part due to their potential as catalysts for the synthesis of NH 3 , a key ingredient in fertilizer and precursor to industrial chemicals. Of particular interest are the ternary (ABN) or higher-order nitrides with high metal-to-nitrogen ratios that show promise in enhancing NH 3 synthesis reaction rates and yields via heterogeneous catalysis or chemical looping. Although metal nitrides are predicted to be numerous, the stability of nitrogen triple bonds found in N 2 , especially in comparison to the metal–nitrogen bonds, has considerably hindered synthetic efforts to produce complex nitride compounds. In this study, we present an exhaustive down-selection process to identify ternary nitrides for a promising chemical looping NH 3 production mechanism. We also report on a facile and efficient two-step synthesis method that can produce well-characterized η-carbide Co 3 Mo 3 N/Fe 3 Mo 3 N or filled β-manganese Ni 2 Mo 3 N ternaries, as well as their associated quaternary, (Co,Fe) 3 Mo 3 N, (Fe,Ni) 2 Mo 3 N, and (Co,Ni) 2 Mo 3 N, solid solutions. To further explore the quaternary space, syntheses of (Co,Ni) 3 Mo 3 N (Ni ≤ 10 mol %) and Co 3 (Mo,W) 3 N (W ≤ 10 mol %) were also investigated. The structures of the nitrides were characterized via X-ray powder diffraction. The morphology and compositions were characterized with scanning electron microscopy. In conclusion, the multitude of chemically unique, but structurally related, nitrides suggests that properties such as nitrogen activity may be tunable, making the materials of great interest for NH 3 synthesis schemes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical Simulations in Support of a Long-Term Test of Gas Production From Hydrate Accumulations on the Alaska North Slope: Water Production and Associated Design and Management Issues

Here, we investigated numerical simulation strategies for a long-term test of depressurization-induced gas production from the B1 Sand of Unit B at the Hydrate-01 Stratigraphic Test Well. The main objective of this study was to estimate fluid production rates (with emphasis on water production) under a variety of conditions and production scenarios and contribute new insights to the design and management of the field test. In the first part of the study, we investigated the system response to a three-step depressurization process using two limiting sets of flow properties─the expected maximum and minimum intrinsic and effective permeabilities─for the very heterogeneous reservoir. In the second part, we investigated the effect of the production interval length and placement within the formation relative to the boundaries of the hydrate-bearing unit. The best performing well configuration was used in the third part of the study, which used the most representative subsurface flow properties to investigate the effect of the depressurization strategy on the production performance. The best overall performance (largest gas production with modest water production and a strong response at the observation wells) was obtained with a 10 m-long well situated 3 m below the top of the formation and a three-step depressurization scheme at 15-day intervals to a terminal bottomhole pressure of 2.8 MPa. The overall production performance was enhanced by a faster rate of depressurization. Estimated water production rates in all cases were limited and easily manageable. None of the tested well configurations or depressurization strategies significantly reduced water production without also severely reducing gas production. In all the investigated cases, 95% of the long-term fraction of produced water was replenished by inflows from the boundaries and could not be reduced. These substantial water inflows are an unavoidable feature of HU-B and cannot be easily mitigated by a hydraulic control.

02 PETROLEUM↗

Addition of tabulated equation of state and neutrino leakage support to illinoisgrmhd

Here we have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to illinoisgrmhd, an open-source GRMHD code for dynamical spacetimes in the einstein toolkit. These new features are provided by two new, nrpy+-based codes: nrpyeos, which performs highly efficient EOS table lookups and interpolations, and nrpyleakage, which implements a new, adaptive mesh refinement (AMR)-capable neutrino leakage scheme in the einstein toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of illinoisgrmhd, as well as nrpyeos and nrpyleakage, is pedagogically documented in jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the einstein toolkit.

79 ASTRONOMY AND ASTROPHYSICS↗

Multistage mixed precision iterative refinement

Abstract Low precision arithmetic, in particular half precision (16‐bit) floating point arithmetic, is now available in commercial hardware. Using lower precision can offer significant savings in computation and communication costs with proportional savings in energy. Motivated by this, there has been a renewed interest in mixed precision iterative refinement schemes for solving linear systems , and new variants of GMRES‐based iterative refinement have been developed. Each particular variant with a given combination of precisions leads to different condition number‐based constraints for convergence of the backward and forward errors, and each has different performance costs. The constraints for convergence given in the literature are, as an artifact of the analyses, often overly strict in practice, and thus could lead a user to select a more expensive variant when a less expensive one would have sufficed. In this work, we develop a multistage mixed precision iterative refinement solver which aims to combine existing mixed precision approaches to balance performance and accuracy and improve usability. For a user‐specified initial combination of precisions, the algorithm begins with the least expensive approach and convergence is monitored via inexpensive computations with quantities produced during the iteration. If slow convergence or divergence is detected using particular stopping criteria, the algorithm switches to use a more expensive, but more reliable variant. A novel aspect of our approach is that, unlike existing implementations, our algorithm first attempts to use “stronger” GMRES‐based solvers for the solution update before resorting to increasing the precision(s). In some scenarios, this can avoid the need to refactorize the matrix in higher precision. We perform extensive numerical experiments on a variety of random dense problems and problems from real applications which confirm the benefits of the multistage approach.

Oktay, Eda↗

Cooperative Effects Associated with High Electrolyte Concentrations in Driving the Conversion of CO2 to C2H4 on Copper

Compared to a conventional electrolyte concentration of 1 M HCOOK, the use of a highly concentrated 7.1 M HCOOK electrolyte increases the Faradaic efficiency (FE) ratio of C2H4/CO from 2.2 +- 0.3 to 18.3 +- 4.8 at -1.08 V vs. reversible hydrogen electrode (RHE) on a Cu gas-diffusion electrode. Based on electrochemical analysis and ab initio molecular dynamics (AIMD) simulation, the identity and concentration of the cation and anion play more important roles in controlling the CO2R reaction pathway than the bulk CO2 solubility and the bulk pH of electrolytes. In situ attenuated reflectance surface enhanced infrared absorption spectroscopy (ATR-SEIRAS) suggests that, unlike 1 M HCOOK, the *CO-bridge-binding mode on Cu is dominant in 7.1 M HCOOK electrolyte, which potentially results in less CO release and higher yield of C2H4. This study demonstrates that although we can tailor the electrolyte composition to shift product selectivity, the factors that control the product selectivity are numerous and cannot be distilled down into one correlated property-reactivity relationship.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗

Development of simplified model for injection rate prediction of diesel injectors during transient and steady operation

Determining the injection rate of the diesel injector during the transient short injection operations is a critical matter to enhance the combustion process by improving air/fuel mixture formation. In this regard, numerous injection models have been suggested in the literature to predict the injection rate accurately. However, such algorithms have not precisely predicted the rate of injection (ROI) in transient short injection events. In addition, most of those models require many input parameters to characterize hydraulic and mechanical subsystems in the injector, which can make them inconvenient to tune for different injector specifications, in other words not suitable for universal usage. Therefore, our study aims to develop a simple injection model that can predict injection rate precisely in the short and long injection durations considering relatively fewer input parameters. The model traces the pressure variation inside the sac based on the bulk modulus theory to predict nozzle exit velocity, thus the injection rate of the injector. In the algorithm, the bulk modulus of fuel (E) and discharge coefficients of the nozzle hole and needle seat (C d,H and C d,N ) have been considered to vary rather than set constant during the injection process since such flow parameters would have a critical impact on the flow development at the transient stage, thus affecting the injection rate of diesel injectors. Here with this modelling strategy, the accuracy of the new model on ROI prediction was confirmed with acceptable model error for both short and long injection durations by validating against experimental results. To understand the effect of the new model scheme on ROI prediction, ROI results with the new model have been discussed by comparing it with the ROI results predicted with constant flow parameters.

42 ENGINEERING↗

Data-Driven Kinetic Reaction Networks for Separation Chemistry

Understanding complex, multistep chemical reactions at the molecular level is a major challenge whose solution would greatly benefit the design and optimization of numerous chemical processes. The separation of rare-earth (4f) and actinide (5f) elements is an example where improving our chemical understanding is important for designing and optimizing new chemistries, even with a limited number of observations. Here, in this work, we leverage data-driven artificial intelligence and machine-learning approaches to develop kinetic reaction networks that describe the liquid–liquid extraction mechanism of uranium using N,N-di-2-ethylhexyl-isobutyramide (DEHiBA). Specifically, we compare and contrast the properties of two classes of models: (1) purely data-driven models that are regularized using chemistry-agnostic, L1 regression and (2) chemistry-informed models that are regularized using relative reaction energies provided by quantum mechanical calculations. We observe that purely data-driven models are unbiased, simple, and accurate in their predictions of experimental measurements when provided with sufficient data but are difficult to fully constrain and interpret. In contrast, chemistry-informed models exhibit significantly improved chemical interpretability and consistency, providing a detailed description of the separation process while achieving high accuracy through ensemble averaging. Overall, the dominant species predicted to be extracted into the organic phase is UO 2 (NO 3 ) 2 (DEHiBA) 2 , agreeing with experimental slope analysis, thermodynamic modeling, EXAFS, and crystal structures. This work demonstrates that leveraging the fundamental structure of the problem can lead to efficient learning schemes that provide both accurate predictions and chemical insights at a low computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dependence of Convective Cloud Microphysical Properties on Environmental Conditions during the TRACER and ESCAPE Field Campaigns: A Synergistic Approach of Observations, Machine Learning and Parcel Models

The sensitivity of convective clouds to aerosols and their interactions with environment, combined with limited observational constraints in parameterizations, introduces significant uncertainties in atmospheric models. Here, this study investigates the dependence of convective cloud microphysical properties on environmental conditions using a synergistic approach that combines unique observations from the TRACER and ESCAPE field campaigns, machine learning techniques, and parcel model simulations with a super-droplet microphysics scheme. A random forest algorithm identifies in-situ vertical velocity (w), temperature (T), and surface fine-mode aerosol mass concentration as the three most important environmental conditions influencing cloud properties including liquid water content (LWC), number concentration for particles with D max < 50 μm (N c ,<50), 50 μm ≤ D max ≤ 3000 μm (N c,50–3000 ), and droplet effective diameter (D e ). Results show that LWC, N c,<50 , and N c,50–3000 significantly increase with w in updrafts. Across w bins, as T decreases, LWC, D e , and N c,50–3000 increase, while N c,<50 decreases, which are closely linked to the distance above cloud bases. Warmer cloud bases yield higher LWC, greater N c,50–3000 , and smaller N c,<50 , while polluted environments produce greater N c,<50 . Parcel model simulations successfully replicate these observed dependencies. The simulation results indicate that warmer cloud bases enhance condensation generating larger droplets, and differences in droplet sizes are then amplified through collision-coalescence, resulting in a greater N c,50–3000 . Polluted conditions result in a greater N c,<50 primarily due to enhanced cloud condensation nuclei activation despite increased collision-coalescence rates compared to pristine conditions. This study provides observed quantitative patterns characterizing cloud microphysical properties as a function of key environmental parameters, offering valuable constraints for improving physics parameterizations and numerical models.

54 ENVIRONMENTAL SCIENCES↗

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids

This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.

graphics processing units↗

Hybrid Basis and Multi-Center Grid Method for Strong-Field Processes

We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗

Responses of summer mesoscale convective systems to irrigation over the North China Plain based on convection-permitting model simulations

Extensive irrigation activities in the North China Plain (NCP) significantly influence regional weather and climate. However, previous studies focusing on the NCP were primarily based on coarse-resolution models, which are unable to explicitly resolve convection systems, causing large uncertainty in precipitation simulations. In this study, a convection-permitting model coupled with a dynamic irrigation scheme is utilized to investigate the impacts of irrigation on summertime mesoscale convective systems (MCSs) over the NCP. Sensitivity experiments with irrigation off and on are conducted for 5 summers and an MCS identification and tracking algorithm is applied to both satellite observations and model simulations. We find that incorporating irrigation in the model increases MCS precipitation, which agrees more with observations. The probability distributions of MCS lifetime, area, propagation speed, and intensity are all better simulated with irrigation. Irrigation increases the occurrence frequency of MCSs throughout the entire day. The nighttime increase is partly because of more frequent local initiation of MCS developed from isolated deep convection, while the daytime increase is mainly attributed to the changes in MCSs initiating elsewhere and then propagating to the NCP. On average, irrigation induces additional moisture that is more thermodynamically favorable for precipitation, but this effect is partially offset by the weakened ascending air motion primarily caused by irrigation surface cooling. Compared to weak MCS precipitation events, strong MCS precipitation events experience greater enhancement in precipitation intensity when including irrigation because the offset effect from the change in large-scale ascending air motion is insignificant. In addition, irrigation makes the variation of MCS precipitation intensity more correlated with the variation in ascending motion but less correlated with that in atmospheric moisture. Our results suggest the pronounced impacts of irrigation on MCSs over the NCP which should be included in numerical models to improve regional precipitation simulation and prediction.

54 ENVIRONMENTAL SCIENCES↗

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

Data-driven surrogates for high dimensional models using Gaussian process regression on the Grassmann manifold

This paper introduces a surrogate modeling scheme based on Grassmannian manifold learning to be used for cost-efficient predictions of high-dimensional stochastic systems. The method exploits subspace-structured features of each solution by projecting it onto a Grassmann manifold. This point-wise linear dimensionality reduction harnesses the structural information to assess the similarity between solutions at different points in the input parameter space. The method utilizes a solution clustering approach in order to identify regions of the parameter space over which solutions are sufficiently similarly such that they can be interpolated on the Grassmannian. In this clustering, the reduced-order solutions are partitioned into disjoint clusters on the Grassmann manifold using the eigen-structure of properly defined Grassmannian kernels and, the Karcher mean of each cluster is estimated. Then, the points in each cluster are projected onto the tangent space with origin at the corresponding Karcher mean using the exponential mapping. For each cluster, a Gaussian process regression model is trained that maps the input parameters of the system to the reduced solution points of the corresponding cluster projected onto the tangent space. Using this Gaussian process model, the full-field solution can be efficiently predicted at any new point in the parameter space. In certain cases, the solution clusters will span disjoint regions of the parameter space. In such cases, for each of the solution clusters we utilize a second, density-based spatial clustering to group their corresponding input parameter points in the Euclidean space. The proposed method is applied to two numerical examples. Here, the first is a nonlinear stochastic ordinary differential equation with uncertain initial conditions where the surrogate is used to predict the time history solution. The second involves modeling of plastic deformation in a model amorphous solid using the Shear Transformation Zone theory of plasticity, where the proposed surrogate is used to predict the full strain field of a material specimen under large shear strains.

42 ENGINEERING↗

Numerical Simulations in Support of a Long-Term Test of Gas Production from Hydrate Accumulations on the Alaska North Slope: Reservoir Response to Interruptions of Production (Shut-Ins)

In this work, we investigate by means of numerical simulation a planned year-long field test of depressurization-induced production from a permafrost-associated hydrate reservoir on the Alaska North Slope at the site of the recently drilled Hydrate-01 Stratigraphic Test Well. The main objective of this study is to assess quantitatively the impact of temporary interruptions (well shut-ins) on the expected fluid production performance from the B1 Sand of the stratigraphic Unit B during controlled depressurization over different time scales, as well as on other relevant aspects of the system response that have the potential to significantly affect the design of the field test. We consider eight different cases of depressurization, including (a) rapid depressurization over a 60-day period to a terminal bottomhole pressure P W of 2.8 MPa and (b) a slower depressurization rate to a final P W of 0.6 MPa at the end of the year-long production test, in addition to (c) a multi-step depressurization regime and (d) a quasi-linear continuous depressurization strategy. The results of the study indicate that shut-ins obviously reduce gas release and production during and immediately after their occurrence, but their longer-term effects are strongly dependent on the depressurization regime and on the time of observation, covering the entire range of potential outcomes. Shut-ins (a) have a universally strong negative effect when quasi-linear depressurization is involved regardless of the length of the production period, and (b) have a strong positive effect in multi-step depressurization schemes that becomes apparent earlier for large initial pressure drops, but (c) can also appear to have practically no effect for slow stepwise depressurization at the end of the year-long production test. Shut-ins lead to a rapid reformation of hydrates, even to the point of disappearance of a free gas phase in the reservoir. Rapid depressurization regimes lead to early maximum rates of hydrate dissociation and gas production, while the maximum rates occur at the end of the production test for the cases of slower depressurization. Shut-ins do not appear to have a significant impact on water production, as the cessation of production is followed by higher rates production when depressurization resumes. Similarly, (a) the fraction of produced CH 4 originating from exsolution from the water, (b) the water-to-gas ratio, and (c) the rate of replenishment of produced water by boundary inflows do not appear significantly affected by shut-ins, the effects of which seem to be temporary in the majority of the cases. The study confirmed the superiority of multi-step depressurization methods as the most effective strategies for hydrate dissociation and gas production and showed that two observation wells (located at distances of 30 and 50 m from the production well) are appropriately positioned and both able to capture the P, T, and S G behavior during the fluid production and shut-ins in any of the eight cases we investigated.

03 NATURAL GAS↗

HarDWR - Harmonized Water Rights Records

For a detailed description of the database of which this record is only one part, please see the HarDWR meta-record. Here we present a new dataset of western U.S. water rights records. This dataset provides consistent unique identifiers for each spatial unit of water management across the domain, unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of 7 broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of intersectoral water allocation changes through water markets (Grogan et al., in review). Specifically, the data were formatted for use as input to a process-based hydrologic model, WBM, with a water rights module (Grogan et al., in review). While this specific study motivated the development of the database presented here, U.S. west water management is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al, 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. west. The raw downloaded data for each state is described in Lisk et al. (in review), as well as here. The dataset is a series of various files organized by state sub-directories. The first two characters of each file name is the abbreviation for the state the in which the file contains data for. After the abbreviation is the text which describes the contents of the file. Here is each file type described in detail: XXFullHarmonizedRights.csv: A file of the combined groundwater and surface water records for each state. Essentially, this file is the merging of XXGroundwaterHarmonizedRights.csv and XXSurfaceWaterHarmonizedRights.csv by state. The column headers for each of this type of file are: state - The name of the state the data comes from. FIPS - The two digit numeric state ID code. waterRightID - The unique identifying ID of the water right, the same identifier as its state uses. priorityDate - The priority date associated with the right. origWaterUse - The original stated water use from the state. waterUse - The water use category under the unified use categories established here. source - Whether the right is for surface water or groundwater. basinNum - The alpha-numeric identifier of the WMA the record belongs to. CFS - The maximum flow of the allocation in cubic feet per second (ft3s-1). Arizona is unique among the states, as its surface and groundwater resources are managed with two different sets of boundaries. So, for Arizona, the basinNum column is missing and instead there are two columns: surBasinNum - The alpha-numeric identifier of the surface water WMA the record belongs to. grdBasinNum - The alpha-numeric identifier of the groundwater WMA the record belongs to. XXStatePOD.shp: A shapefile which identifies the location of the Points of Diversion for the state's water rights. It should be noted that not all water right records in XXFullHarmonizedRights.csv have coordinates, and therefore may be missing from this file. XXStatePOU.shp: A shapefile which contains the area(s) in which each water right is claimed to be used. Currently, only Idaho and Washington provided valid data to included within this file. XXGroundwaterHarmonizedRights.csv: A file which contains only harmonized groundwater rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. XXSurfaceWaterHarmonizedRights.csv: A file which contains only harmonized surface water rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. Additionally, one file, stateWMALabels.csv, is not stored within a sub-directory. While we have referred to the spatial boundaries that each state uses to manage its water resources as WMAs, this term is not shared across all states. This file lists the proper name for each boundary set, by state.

Economics↗