Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Volumetric”

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 217 records · Page 12

Gas flow through rotameters

Using data available for small rotameters that use spherical floats in gas flow, a linear relationship is derived. It is noted that the relationship provides a good fit for variable volumetric flow, density, and viscosity at constant flow height. With low Reynolds numbers (Re being less than 1), the product of the variable volumetric flow and the viscosity becomes constant; at high Reynolds numbers (Re being greater than 2000), the product of the variable volumetric flow and the square root of the density becomes constant. It is pointed out that the equation given here can be used to obtain an indirect calibration with any gas of known density and viscosity. The constancy of the product of the variable volumetric flow and viscosity at low variable volumetric flows is seen as suggesting the development of simple, inexpensive gas viscometers using rotameter technology.

Levin, H.↗

Initial Thrust Measurements of Marshall's Ion-ioN Thruster

Electronegative ion thrusters are a variation of traditional gridded ion thruster technology differentiated by the production and acceleration of both positive and negative ions. Benefits of electronegative ion thrusters include the elimination of lifetime-limiting cathodes from the thruster architecture and the ability to generate appreciable thrust from both charge species. While much progress has been made in the development of electronegative ion thruster technology, direct thrust measurements are required to unambiguously demonstrate the efficacy of the concept and support continued development. In the present work, direct thrust measurements of the thrust produced by the MINT (Marshall's Ion-ioN Thruster) are performed using an inverted-pendulum thrust stand in the High-Power Electric Propulsion Laboratory's Vacuum Test Facility-1 at the Georgia Institute of Technology with operating pressures ranging from 4.8 x 10(exp -5) and 5.7 x 10(exp -5) torr. Thrust is recorded while operating with a propellant volumetric mixture ratio of 5:1 argon to nitrogen with total volumetric flow rates of 6, 12, and 24 sccm (0.17, 0.34, and 0.68 mg/s). Plasma is generated using a helical antenna at 13.56 MHz and radio frequency (RF) power levels of 150 and 350 W. The acceleration grid assembly is operated using both sinusoidal and square waveform biases of +/-350 V at frequencies of 4, 10, 25, 125, and 225 kHz. Thrust is recorded for two separate thruster configurations: with and without the magnetic filter. No thrust is discernable during thruster operation without the magnetic filter for any volumetric flow rate, RF forward Power level, or acceleration grid biasing scheme. For the full thruster configuration, with the magnetic filter installed, a brief burst of thrust of approximately 3.75 mN +/- 3 mN of error is observed at the start of grid operation for a volumetric flow rate of 24 sccm at 350 W RF power using a sinusoidal waveform grid bias at 125 kHz and +/- 350 V. Similar bursts in thrust are observed using a square waveform grid bias at 10 kHz and +/- 350 V for volumetric flow rates of 6, 10, and 12 sccm at 150, 350, and 350 W respectively. The only operating condition that exhibits repeated thrust spikes throughout thruster operation is the 24 sccm condition with a 5:1 mixture ratio at 150 W RF power using the 10 kHz square waveform acceleration grid bias. Thrust spikes for this condition measure 3 mN with an error of +/- 2.5 mN. There are no operating conditions tested that show continuous thrust production.

Caruso, Natalie R. S.↗

Geomechanical and Hydrogeological Evaluation of a Shallow Hydraulic Fracture at the Devine Fracture Pilot Site, Medina County, Texas

UT-Austin’s Devine Fracture Pilot Site, 50 miles southwest of San Antonio, Texas, has been targeted for a comprehensive, multidisciplinary development of fracture diagnostic techniques that are cross-validated by ground-truth data acquisition near a recently created, 175-ft-deep, horizontal hydraulic fracture (Ahmadian et al. 2018 Demonstration of proof of concept of electromagnetic geophysical methods for high resolution illumination of induced fracture networks. In Proceedings of the SPE Hydraulic Fracturing Technology Conference and Exhibition, The Woodlands, Texas, USA, 23–25 January 2018. SPE-189858-MS.). To evaluate the fracture diagnostic methods at this site, we conducted injection tests with a predefined volumetric flow-rate profile, resembling a diagnostic fracture injection test on September 2020. Subsequently, we developed hydrogeological and geomechanical models based on flow-rate and bottomhole-pressure measurements. History-matching efforts using a simplified layer-cake hydrogeological model resulted in the field-scale formation permeability of 9.87 × 10 –15 m 2 (10 mD) and Darcy-scale fracture permeability. The analysis of the bottomhole pressure and injection-rate history showed that (1) the newly created horizontal fracture was closed adjacent to the injection well pre-injection and (2) the initial pump-pressure increase at a nominal volumetric injection rate led to near-well fracture reopening, fluid conductivity increase, and abrupt injection-rate increase. To overcome hydrogeological-model limitations of predicting fracture reopening throughout injection, we extended the modeling to a finite-element, poroelastic analysis of horizontal-fracture growth using a cohesive-zone model. Using this fracture-reopening model, we improved the history match of the transient-pressure response during the experiment by adjusting the hydromechanical properties. Furthermore, post-injection pressure transient analyses helped reduce uncertainty in the overburden-stress gradient, and the initial hydraulic-fracturing simulation verified the plausibility of achieving the surveyed propped fracture area.

02 PETROLEUM↗

MFA-DVR: direct volume rendering of MFA models

3D volume rendering is widely used to reveal insightful intrinsic patterns of volumetric datasets across many domains. However, the complex structures and varying scales of volumetric data can make efficiently generating high-quality volume rendering results a challenging task. Multivariate functional approximation (MFA) is a new data model that addresses some of the critical challenges: high-order evaluation of both value and derivative anywhere in the spatial domain, compact representation for large-scale volumetric data, and uniform representation of both structured and unstructured data. Here, in this paper, we present MFA-DVR, the first direct volume rendering pipeline utilizing the MFA model, for both structured and unstructured volumetric datasets. We demonstrate improved rendering quality using MFA-DVR on both synthetic and real datasets through a comparative study. We show that MFA-DVR not only generates more faithful volume rendering than using local filters but also performs faster on high-order interpolations on structured and unstructured datasets. MFA-DVR is implemented in the existing volume rendering pipeline of the Visualization Toolkit (VTK) to be accessible by the scientific visualization community.

97 MATHEMATICS AND COMPUTING↗

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗

Solid 3D Li–S Battery Design via Stacking 2D Conductive Microporous Coordination Polymers and Amorphous Li–S Layers

To make a lithium-sulfur (Li-S) battery practical, not only high gravimetric energy capacity is important, but also high volumetric energy capacity will be required. The currently explored Li-S cathode designs often deploy systems with liquid electrolyte infiltration, hence with relatively low volumetric capacity. In the current study, we theoretically test a compact solid three-dimensional (3D) design (more like a Li-ion battery cathode than a conventional Li-S cathode) consisted of a sandwich structure alternating between the two-dimensional (2D) Mn-hexaaminobenzene-based coordination polymer (2D Mn-HAB-CP) layer and the amorphous Li-S layer. We study the theoretical limits for both its gravimetric and volumetric energy capacity, as well as its structural stability and Li diffusion within the cathode system. To study the Li diffusion within an amorphous system, we also develop a pull-atom molecular dynamics (PA-MD) to calculate the barrier heights of such disordered systems. We reveal the mechanism that determines the Li diffusion in the amorphous layer of the system. Overall, we find such a 3D solid Li-S cathode can be practical, with sufficient large gravimetric and volumetric energy capacity, as well as the Li diffusion constant. It also solves many other common Li-S cathode problems, from Li polysulfide dissolution to electrical insulating, and structure instabilities.

25 ENERGY STORAGE↗

From core collapse to superluminous: the rates of massive stellar explosions from the Palomar Transient Factory

ABSTRACT We present measurements of the local core-collapse supernova (CCSN) rate using SN discoveries from the Palomar Transient Factory (PTF). We use a Monte Carlo simulation of hundreds of millions of SN light-curve realizations coupled with the detailed PTF survey detection efficiencies to forward model the SN rates in PTF. Using a sample of 86 CCSNe, including 26 stripped-envelope SNe (SESNe), we show that the overall CCSN volumetric rate is $r^\mathrm{CC}_v=9.10_{-1.27}^{+1.56}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$ at 〈z〉 = 0.028, and the SESN volumetric rate is $r^\mathrm{SE}_v=2.41_{-0.64}^{+0.81}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$. We further measure a volumetric rate for hydrogen-free superluminous SNe (SLSNe-I) using eight events at z ≤ 0.2 of $r^\mathrm{SLSN-I}_v=35_{-13}^{+25}\, \text{SNe yr}^{-1}\text{Gpc}^{-3}\, h_{70}^{3}$, which represents the most precise SLSN-I rate measurement to date. Using a simple cosmic star formation history to adjust these volumetric rate measurements to the same redshift, we measure a local ratio of SLSN-I to SESN of ${\sim}1/810^{+1500}_{-94}$, and of SLSN-I to all CCSN types of ${\sim}1/3500^{+2800}_{-720}$. However, using host galaxy stellar mass as a proxy for metallicity, we also show that this ratio is strongly metallicity dependent: in low-mass (logM* < 9.5 M⊙) galaxies, which are the only environments that host SLSN-I in our sample, we measure an SLSN-I to SESN fraction of $1/300^{+380}_{-170}$ and $1/1700^{+1800}_{-720}$ for all CCSN. We further investigate the SN rates a function of host galaxy stellar mass, and show that the specific rates of all CCSNe decrease with increasing stellar mass.

79 ASTRONOMY AND ASTROPHYSICS↗

Attribute-Aware RBFs: Interactive Visualization of Time Series Particle Volumes Using RT Core Range Queries

Smoothed-particle hydrodynamics (SPH) is a mesh-free method used to simulate volumetric media in fluids, astrophysics, and solid mechanics. Visualizing these simulations is problematic because these datasets often contain millions, if not billions of particles carrying physical attributes and moving over time. Radial basis functions (RBFs) are used to model particles, and overlapping particles are interpolated to reconstruct a high-quality volumetric field; however, this interpolation process is expensive and makes interactive visualization difficult. Existing RBF interpolation schemes do not account for color-mapped attributes and are instead constrained to visualizing just the density field. To address these challenges, we exploit ray tracing cores in modern GPU architectures to accelerate scalar field reconstruction. We use a novel RBF interpolation scheme to integrate per-particle colors and densities, and leverage GPU-parallel tree construction and refitting to quickly update the tree as the simulation animates over time or when the user manipulates particle radii. We also propose a Hilbert reordering scheme to cluster particles together at the leaves of the tree to reduce tree memory consumption. Finally, we reduce the noise of volumetric shadows by adopting a spatially temporal blue noise sampling scheme. Our method can provide a more detailed and interactive view of these large, volumetric, time-series particle datasets than traditional methods, leading to new insights into these physics simulations.

Particle Volumes↗

F-Hash: Feature-Based Hash Design for Time-Varying Volume Visualization via Multi-Resolution Tesseract Encoding

Interactive time-varying volume visualization is challenging due to its complex spatiotemporal features and sheer size of the dataset. Recent works transform the original discrete time-varying volumetric data into continuous Implicit Neural Representations (INR) to address the issues of compression, rendering, and super-resolution in both spatial and temporal domains. However, training the INR takes a long time to converge, especially when handling large-scale time-varying volumetric datasets. In this work, we proposed F-Hash, a novel feature-based multi-resolution Tesseract encoding architecture to greatly enhance the convergence speed compared with existing input encoding methods for modeling time-varying volumetric data. The proposed design incorporates multi-level collision-free hash functions that map dynamic 4D multi-resolution embedding grids without bucket waste, achieving high encoding capacity with compact encoding parameters. Our encoding method is agnostic to time-varying feature detection methods, making it a unified encoding solution for feature tracking and evolution visualization. Experiments show the F-Hash achieves state-of-the-art convergence speed in training various time-varying volumetric datasets for diverse features. We also proposed an adaptive ray marching algorithm to optimize the sample streaming for faster rendering of the time-varying neural representation.

deep learning↗

Quantifying Trapped Powder in Electron Beam Powder Bed Fusion

Abstract Electron beam powder bed fusion (PBF-EB) shows great potential for manufacturing complex parts including those with internal cavities for heat exchanger, manifold systems, or energy absorption purposes. PBF-EB allows for the manufacture of channel geometries without the need for support structures. Due to the nature of the powder spreading process, powder feedstock is often trapped in intentionally manufactured cavities. This trapped powder can often be difficult to remove and can disturb the intended flow of fluid through the cavity or damage downstream components in its use case. These trapped powder particles present a risk of contamination and component failure if not completely evacuated. Ti6Al4V is a choice material for aerospace applications due to its high strength to weight ratio and its composition as a nonferrous metal; however, in weight sensitive applications excess entrapped powders or powders loosely attached to the surface could cause undesirable weight increases. The inherent spreading process of PBF-EB is different than laser powder bed fusion (PBF-LB) in its operational temperature, sintering. In addition, PBF-EB is less commonly studied in literature compared to its PBF-LB counterpart, and as a result the complexity of the semi-sintered powder and its spreading behavior are not well understood. Prior work has investigated the difficulty in removing trapped powder from PBF-EB, but these studies do not address how to quantify the amount of trapped powder in the cavity. Thus, an accurate method to measure the amount of trapped powder in the cavity must be investigated. In this work, Ti6Al4V coupons were manufactured with horizontal and vertical cavities of three different sizes. Archimedes testing allows for the determination of density differences caused by porosity and trapped powders by measuring mass and volumetric dispersion. Computed tomography (CT) is well suited for segmenting the internal structure and features of a part and has been studied for applications including voids, porosity, and dross. Thus, CT was explored as a method for evaluating trapped powder content in this work. The volumetric representation of the segmentation of the reconstructed CT volume can vary greatly depending on the input filter and thresholding methods. In this study, four different types of segmentation approaches were evaluated to determine the best approach for segmenting the volume as compared to an operator labeled ground truth. The percentage density results from the Archimedes testing were compared to the volumetric percent density from the computed tomography approach. Differences in packing density between two different internal channel features were investigated. Overall, this work sought to validate the use of computed tomography for the detection of trapped powders and present a framework for volumetric segmentation.

Johnstone, Brian↗

Surfactant Effects in Porous Electrodes for Microemulsion Redox Flow Batteries

The effect of surfactant additives on electrochemical behavior in porous electrodes was investigated using vanadium redox flow battery half-cells and the dependence of volumetric kinetics and mass transport on electrolyte, surfactant, and electrode type was explored. Without surfactant added, carbon paper electrodes demonstrated greater kinetics and transport compared to carbon felt, for a given electrolyte. Additionally, posolyte kinetics are greater than negolyte kinetics by one to three orders of magnitude, depending on the electrode type. Addition of surfactant increased electrode wettability and possibly electrochemical surface area. However, this was accompanied by a decrease in volumetric mass transport, due to stronger electrolyte-electrode interactions. The presence of sodium dodecyl sulfate (SDS) influenced posolyte and negolyte kinetics differently. Kinetics showed a dependence on electrode type and surfactant. On carbon felt, volumetric kinetics decreased for both posolyte and negolyte with SDS addition. On carbon paper, SDS decreased volumetric kinetics for the posolyte but increased (>2X) kinetics for the negolyte! This kinetic enhancement depends on surfactant chemistry: cetyltrimethylammonium bromide, a cationic surfactant, failed to increase kinetics. Furthermore, SDS did not increase areal specific resistance. Furthermore, these findings show the superior performance of carbon paper compared to carbon felt and suggest SDS as a possible VRFB negolyte additive.

25 ENERGY STORAGE↗

Data for Resourceful and Economical Designing of Fermentation Medium for Lab and Commercial Strains of Yeast from Alternative Feedstock: Transgenic Oilcane

Sugarcane plant engineered to accumulate lipids in its vegetative tissue is being developed as a new bioenergy crop. The new crop would be a source of juice, oil, and cellulosic sugars. However, limited tolerance of industrially recognized yeasts towards inhibitors generated during the processing of lignocellulosic biomass to produce fermentable sugars is a major challenge in developing scalable processes for second-generation drop-in fuel production. To this end, hydrolysates generated from engineered sugarcane—‘oilcane’ bagasse contain added phenolics and fatty acids that further restrict the growth of fermenting microorganisms and necessitate nutrient supplementation and/or detoxification of hydrolysate which makes the fermentation process expensive. Herein, we propose a resourceful and economical approach for growing lab and commercial strains of S. cerevisiae on unrefined cellulosic sugars aerobically and fermentatively. An equal ratio of hydrolysate and juice was found optimum for growth and fermentation by lab and commercial strains of Saccharomyces cerevisiae engineered for xylose fermentation. The industrial strain grew and fermented efficiently under low aeration conditions having an ethanol titer, yield, specific and volumetric productivities of 46.96 ± 0.19 g/l, 0.51 ± 0.00 g/g, 0.27 ± 0.02 g/g.h and 1.95 ± 0.01 g/l.h, respectively, while the lab strain grew better under higher aeration conditions having the ethanol titer, yield, specific and volumetric productivities of 24.93 ± 0.09, 0.27 ± 0.00 g/g, 0.17 ± 0.00 g/g.h and 1.04 ± 0.00 g/l.h, respectively. Acclimation of cultures in a blended medium significantly improved the performance of the yeast strains. The addition of transgenic oilcane juice, which is inedible and rich in amino acids, to the hydrolysate averted the need for expensive nutrient supplementation and detoxification steps of hydrolysate. The approach provides an economical solution to reduce the cost of fermentation at an industrial scale for second-generation drop-in fuel production.

Biomass Analytics↗

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ultra-High Operation Temperature SiC-matrix Solar Thermal Air Receiver (HOTSSTAR) enabled by additive manufacturing: Test Facility & Performance Evaluations

Solar Heat for Industrial Processes (SHIP) cavity receivers are capable of generating electricity or industrial process heat by absorbing thermal energy from solar radiation, focused on a small area. The concentration of solar radiation on the small area of the receiver enables the achievement of high temperatures (ranging from 400°C to 1,100°C) of a working fluid, thus making the SHIP technology thermodynamically comparable with conventional power plants. A volumetric receiver consists of a porous structure-generally made of silicon carbide or metal, which absorbs solar radiation and converts it into heat energy. Heat energy from the porous materials is then transferred to the fluid following through them. A volumetric receiver acts as a convective heat exchanger, transferring heat to the fluid through convection. Open-loop volumetric receivers work with air at atmospheric pressure and are suitable for single-cycle or multi-cycle energy plants. A Model Based Systems Engineering (MBSE) approach was used to develop a test bed at Sandia national Laboratories (SNL) capable of demonstrating an open-loop volumetric air receiver developed by General Electric Aerospace (GE Aerospace). This paper presents the development of the various MBSE methods, test bed, and testing operations for the GE air receiver, which was experimentally demonstrated to achieve 1,350°C for over 3 hours of operation and an approximate 70% receiver efficiency. By being able to achieve such high temperatures >1,000°C, this work provides the potential to support many SHIP industrial use cases.

14 SOLAR ENERGY↗

TRAILS Output Files

Overview This data repository contains ZIP files that store compressed versions of the output of running the WaterPaths utility planning and management tool in the DU Re-Evaluation mode (to download the tool, please see this GitHub repository). The tool was used to simulate the six-utility North Carolina Research Triangle problem. Details on the contents of each ZIP file can be seen below. Data details Temporal range: Weekly data for 2,344 weeks from 2015 to 2060 (45 years). Spatial range: Six water utilities in the North Carolina Research Triangle region (0: Chapel Hil/OWASA, 1: Durham, 2: Cary, 3: Raleigh, 4: Pittsboro, and 5: Chatham) File types: CSV and OUT Different solutions available The solution numbers correspond to the different pathway strategies (henceforth referred to as "solutions") discussed in paper's main and supporting text (abstract and link to the paper here). They are as follows: Sol92: The Durham-focused pathway strategy Sol132: The Raleigh-focused pathway strategy Sol140: The regionally-robust pathway strategy Objectives files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Objectives_RDMXX_solsXX_to_XX.csv files. Each CSV file will consist of a row representing all the objective values for that specific solution, while every six columns represents the reliability, restriction frequency, infrastructure net present value ($ mil), peak financial cost, worst-case cost, and unit cost ($ per MG; in that order) for each of the six utilities. There will be 1,000 such files, denoting the performance of the six utilities across the 1,000 deeply uncertain states of the world (DU SOWs). Pathway files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Pathways_sXX_RDMXX.out file. Each OUT corresponds to the set of infrastructure being triggered in a specific DU SOW, and each file will have the name file will consist of four tab-delimited columns that are described as follows: Realization: The realization in which an infrastructure options being triggered utility: The utility currently triggering infrastructure week: The week in which a specific infrastructure option is being triggered infra.: The infrastructure option being triggered If the OUT file contains only the header line, no infrastructure was triggered for that specific DU SOW. Policies files These files can be obtained by unzipping Policies.zip. Each of the 1,000 CSV files within the unzipped folder will contain weekly water use restriction policies for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: 0rest_m: restriction multiplier for utility 0 (values between 0 and 1) 1rest_m: restriction multiplier for utility 1 (values between 0 and 1) 2rest_m: restriction multiplier for utility 2 (values between 0 and 1) 3rest_m: restriction multiplier for utility 3 (values between 0 and 1) 4rest_m: restriction multiplier for utility 4 (values between 0 and 1) 5rest_m: restriction multiplier for utility 5 (values between 0 and 1) 0transf: transfer volume for utility 0 (in MGD) 1transf: transfer volume for utility 1 (in MGD) 2transf: transfer volume for utility 2 (in MGD) 3transf: transfer volume for utility 3 (in MGD) 4transf: transfer volume for utility 4 (in MGD) 5transf: transfer volume for utility 5 (in MGD) Water Sources files These files can be obtained by unzipping WaterSources_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each water source for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xvolume: available water volume from source X (in MGD) Xs_area: surface area of source X (in ACF) Xdemand: demand drawn from a water source from source X (in MGD) Xup_spill: upstream spillage from source X (in MGD) Xww_inflow: wastewater inflow from source X (in MGD) Xcatch_inflow: upstream catchment inflow to source X (in MGD) Xevap: evaporation multiplier for source X (values between 0 and 1) Xds_spill: downstream spillage from source X (in MGD) X_Y_alloc_cap: the allocated capacity from source X to utility Y (values between 0 and 1) X_Y_alloc_dem: the allocated demand from source X to utility Y (values between 0 and 1) Xtrmt_alloc_Y: the allocated treatment capacity from source X to utility Y (values between 0 and 1) Utilities files These files can be obtained by unzipping Utilities_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each utility for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xst_vol: total available storage volume of utility X (in MG) Xcapacity: total storage capacity of utility X (in MG) Xnet_inf: : net inflow for all storage infrastructure for utility X (in MGD) Xst_rof: short term ROF for utility X (values between 0 and 1) Xst_stor_rof: short-term storage ROF for utility X (values between 0 and 1) Xst_trmt_rof: short-term treatment ROF for utility X (values between 0 and 1) Xlt_rof: long-term ROF for utility X (values between 0 and 1) Xlt_stor_rof: long-term storage ROF for utility X (values between 0 and 1) Xlt_trmt_rof: long-term treatment ROF for utility X (values between 0 and 1) Xrest_demand: restricted demand for utility X (in MGD) Xunrest_demand: unrestricted demand for utility X (in MGD) Xunfulf_demand: unfulfilled demand for utility X (in MGD) Xwastewater: wastewater return for utility X (in MGD) Xtreat_capacity: total treatment capacity for utility X (in MG) Xcont_fund: reserve (contingency) fund balance for utility X Xins_pout: insurance payout for utility X (% annual volumetric revenue) Xins_price: insurance price for utility X (% annual volumetric revenue) Xinfra_npv: infrastructure net present value for utility ($mil) Xst_vol: total available storage volume of utility X (in MG) Xdebt_serv: debt service for utility X (usually once per year if the infrastructure is triggered; % annual volumetric revenue) Xstor_vol: total stored volume (in MGD) Xobs_ann_dem: observed annual demand for utility X (in MGD) Xproj_dem: projected annual demand for utility X (in MGD) Xpv_debt_serv: present value of debt service payments for utility X (% annual volumetric revenue) Xgross_rev: gross revenue for utility X ($mil) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

Artificial Intelligence↗

Method and device for determining heats of combustion of gaseous hydrocarbons

A method and device is provided for a quick, accurate and on-line determination of heats of combustion of gaseous hydrocarbons. First, the amount of oxygen in the carrier air stream is sensed by an oxygen sensing system. Second, three individual volumetric flow rates of oxygen, carrier stream air, and hydrocrabon test gas are introduced into a burner. The hydrocarbon test gas is fed into the burner at a volumetric flow rate, n, measured by a flowmeter. Third, the amount of oxygen in the resulting combustion products is sensed by an oxygen sensing system. Fourth, the volumetric flow rate of oxygen is adjusted until the amount of oxygen in the combustion product equals the amount of oxygen previously sensed in the carrier air stream. This equalizing volumetric flow rate is m and is measured by a flowmeter. The heat of combustion of the hydrocrabon test gas is then determined from the ratio m/n.

Singh, Jag J.↗

Diaphragm Pump With Resonant Piezoelectric Drive

A diaphragm pump driven by a piezoelectric actuator is undergoing development. This pump is intended to be a prototype of lightweight, highly reliable pumps for circulating cooling liquids in protective garments and high-power electronic circuits, and perhaps for some medical applications. The pump would be highly reliable because it would contain no sliding seals or bearings that could wear, the only parts subject to wear would be two check valves, and the diaphragm and other flexing parts could be designed, by use of proven methods, for extremely long life. Because the pump would be capable of a large volumetric flow rate and would have only a small dead volume, its operation would not be disrupted by ingestion of gas, and it could be started reliably under all conditions. The prior art includes a number piezoelectrically actuated diaphragm pumps. Because of the smallness of the motions of piezoelectric actuators (typical maximum strains only about 0.001), the volumetric flow rates of those pumps are much too small for typical cooling applications. In the pump now undergoing development, mechanical resonance would be utilized to amplify the motion generated by the piezoelectric actuator and thereby multiply the volumetric flow rate. The prime mover in this pump would be a stack of piezoelectric ceramic actuators, one end of which would be connected to a spring that would be part of a spring-and-mass resonator structure. The mass part of the resonator structure would include the pump diaphragm (see Figure 1). Contraction of the spring would draw the diaphragm to the left, causing the volume of the fluid chamber to increase and thereby causing fluid to flow into the chamber. Subsequent expansion of the spring would push the diaphragm to the right, causing the volume of the fluid chamber to decrease, and thereby expelling fluid from the chamber. The fluid would enter and leave the chamber through check valves. The piezoelectric stack would be driven electrically to make it oscillate at the resonance frequency of the spring and- mass structure. This frequency could be made high enough (of the order of 400 Hz) that the masses of all components could be made conveniently small. The resonance would amplify the relatively small motion of the piezoelectric stack (a stroke of the order of 10 m) to a diaphragm stroke of the order of 0.5 mm. The exact amplification factor would depend on the rate of damping of oscillations; this, in turn, would depend on details of design and operation, including (but not limited to) the desired pressure rise and volumetric flow rate. In order to obtain resonance with large displacement, the damping rate must be low enough that the energy imparted to the pumped fluid on each stroke is much less than the kinetic and potential energy exchanged between the mass and spring during each cycle of oscillation. To minimize the power demand of the pump, a highly efficient drive circuit would be used to excite the piezoelectric stack. This circuit (see Figure 2) would amount to a special-purpose regenerative, switching power supply that would operate in a power-source mode during the part of an oscillation cycle when the excitation waveform was positive and in a power-recovery mode during the part of the cycle when the excitation waveform was negative. The circuit would include a voltage-boosting dc-to-dc converter that would convert between a supply potential of 24 Vdc and the high voltage needed to drive the piezoelectric stack. Because of the power-recovery feature, the circuit would consume little power. It should be possible to build the circuit as a compact unit, using readily available components.

Izenson, Michael G.↗

GRACE and GRACE Follow-On Gravity Observations of Intermediate-Depth Earthquakes Contrasted with Those of Shallow Events

Earthquakes involve mass redistribution within the solid Earth and the ocean, and as a result, perturb the Earth's gravitational field. For most of the shallow (<60 km) earthquakes with Mw > 8.0, the GRACE satellite gravity measurements suggest considerable volumetric disturbance of rocks. At a spatial scale of hundreds of km, the effect of volumetric change exceeds gravity change by vertical deformation; for example, negative gravity anomalies associated with volumetric expansion are characteristic patterns after shallow thrust events. In this study, however, we report contrasting observations of gravity change from two intermediate-depth (100–150 km) earthquakes of 2016 & 2017 Mw 8.0 (two combined) Papua New Guinea thrust faulting events and 2019 Mw 8.0 Peru normal faulting and highlight the importance of compressibility in earthquake deformation. The combined 2016/17 thrust events resulted in a positive gravity anomaly of 5–6 microGal around the epicenter, while the 2019 normal faulting produced a negative gravity anomaly of 3–4 microGal. Our modeling found that these gravity changes are manifestation of vertical deformation with limited volumetric change, distinct from gravity changes after the shallow earthquakes. The stronger resistance of rocks to volume change at intermediate-depth results in largely incompressible deformation and thus in a gravity change dominated by vertical deformation. In addition, malleable rocks under high pressure and temperature at depth facilitated substantial afterslip and/or fast viscoelastic relaxation causing additional vertical deformation and gravity change equivalent to the coseismic change. For the Papua New Guinea events, this means that postseismic relaxation enhanced coseismic uplift and relative sea level decrease.

Shin-Chan Han↗