Search NASA⌕ Search

SEARCH · Search NASA

Results for “Injection testing”

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 1,693 records · Page 94

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Isotope ratio characteristics and sensitivity for uranium determinations using a liquid sampling-atmospheric pressure glow discharge ion source coupled to an Orbitrap mass analyzer

Abstract The continued development of the liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma as an ion source for diverse, elemental/isotopic analysis applications continues. To this end, characterization of the capabilities in performing precise and accurate isotope ratio (IR) determinations is essential. Based on past experience with the Thermo Exactive Orbitrap mass analyzer, the LS-APGD was interfaced with this instrument for these tests. While the Orbitrap platform has demonstrated excellent mass resolution and accuracy in “organic” mass spectrometry (MS) applications, work using an Orbitrap for IR analysis is very sparse. These efforts build off previous work in this coupling, where the importance of a few of the LS-APGD discharge parameters and Orbitrap data acquisition methods on IR precision and accuracy were probed. Presented here are the results of a study that evaluated the analytical precision for natural uranium sample (assumed 235U/238U = 0.0072) determinations. The instrumental parameters evaluated include the number of microscans and scans making up a data acquisition set, uranium concentration/signal level, sample make-up, and Fourier transform digitization window. Ultimately, a precision of 0.41% relative standard deviation (RSD) can be achieved for a single determination, with a reproducibility of 1.63 %RSD over 10 separate analytical measurements. A preliminary study of matrix effects on IR measurements of U is presented, highlighting the importance of pre-mass selection before injection into the Orbitrap. The analytical system sensitivity is suggested with the ability to produce a calibration function having an R2 value of >0.99 over a range of 4 orders of magnitude of concentration (~1 – 1000 ng mL-1). These efforts demonstrate the very promising pairing of the LS-APGD ionization source and the Orbitrap, pointing as well to definitive paths forward to better utilize both components in high quality isotope ratio mass spectrometry (IRMS).

Environmental Molecular Sciences Laboratory↗

Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks

Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simplified predictive models for CO 2 sequestration performance assessment

CO2 sequestration in deep saline formations is increasingly being considered as a viable strategy for the mitigation of greenhouse gas emissions from anthropogenic sources. In this context, detailed numerical simulation based models are routinely used to understand key processes and parameters affecting pressure propagation and buoyant plume migration following CO2 injection into the subsurface. As these models are data and computation intensive, the development of computationally-efficient alternatives to conventional numerical simulators has become an active area of research. Such simplified models can be valuable assets during preliminary CO2 injection project screening, serve as a key element of probabilistic system assessment modeling tools, and assist regulators in quickly evaluating geological storage projects. We present three strategies for the development and validation of simplified modeling approaches for CO2 sequestration in deep saline formations: (1) simplified physics-based modeling, (2) statisticallearning based modeling, and (3) reduced-order method based modeling. In the first category, a set of full-physics compositional simulations is used to develop correlations for dimensionless injectivity as a function of the slope of the CO2 fractional-flow curve, variance of layer permeability values, and the nature of vertical permeability arrangement. The same variables, along with a modified gravity number, can be used to develop a correlation for the total storage efficiency within the CO2 plume footprint. Furthermore, the dimensionless average pressure buildup after the onset of boundary effects can be correlated to dimensionless time, CO2 plume footprint, and storativity contrast between the reservoir and caprock. In the second category, statistical “proxy models” are developed using the simulation domain described previously with two approaches: (a) classical Box-Behnken experimental design with a quadratic response surface, and (b) maximin Latin Hypercube sampling (LHS) based design with a multidimensional kriging metamodel fit. For roughly the same number of simulations, the LHS-based metamodel yields a more robust predictive model, as verified by a k-fold cross-validation approach (with data split into training and test sets) as well by validation with an independent dataset. In the third category, a reduced-order modeling procedure is utilized that combines proper orthogonal decomposition (POD) for reducing problem dimensionality with trajectory-piecewise linearization (TPWL) in order to represent system response at new control settings from a limited number of training runs. Significant savings in computational time are observed with reasonable accuracy from the PODTPWL reduced-order model for both vertical and horizontal well problems – which could be important in the context of history matching, uncertainty quantification and optimization problems. The simplified physics and statistical learning based models are also validated using an uncertainty analysis framework. Reference cumulative distribution functions of key model outcomes (i.e., plume radius and reservoir pressure buildup) generated using a 97-run full-physics simulation are successfully validated against the CDF from 10,000 sample probabilistic simulations using the simplified models. The main contribution of this research project is the development and validation of a portfolio of simplified modeling approaches that will enable rapid feasibility and risk assessment for CO2 sequestration in deep saline formations.

54 ENVIRONMENTAL SCIENCES↗

Femtosecond timing distribution and control for next generation accelerators and light sources

Femtosecond Timing Distribution At LCLS Free-electron-lasers (FEL) have the capability of producing high photon flux from the IR to the hard x-ray wavelength range and to emit femtosecond and eventually even attosecond pulses. This makes them an ideal tool for fundamental as well as applied re-search. Timing precision at the Stanford Linear Coherent Light Source (LCLS) between the x-ray FEL (XFEL) and ultrafast optical lasers is currently no better than 100 fs RMS. Ideally this precision should be much better and could be limited only by the x-ray pulse duration, which can be as short as a few femtoseconds. An increasing variety of science problems involving electron and nuclear dynamics in chemical and material systems will become accessible as the timing improves to a few femtoseconds. Advanced methods of electron beam conditioning or pulse injection could allow the FEL to achieve pulse durations less than one femtosecond. The objective of the work described in this proposal is to set up an optical timing distribution system based on mode locked Erbium doped fiber lasers at LCLS facility to improve the timing precision in the facility and allow time stamping with a 10 fs precision. The primary commercial applications for optical timing distributions systems are seen in the worldwide accelerator facilities and next generation light sources community. It is reasonable to expect that at least three major XFELs will be built in the next decade. In addition there will be up to 10 smaller machines, such as FERMI in Italy and Maxlab in Sweden, plus the market for upgrading already existing facilities like Jefferson Lab. The total market is estimated to be on the order of a 100 Million US Dollars. The company owns the exclusive rights to the IP covering the technology enabling sub-10 fs synchronization systems. Testing this technology, which has set records in a lab environment, at LCLS, hence in a real world scenario, is an important corner stone of bringing the technology to market.

43 PARTICLE ACCELERATORS↗

Advanced Accelerators: Particle, Photon and Plasma Wave Interactions

The overall objective of this project was to study the acceleration of electrons to very high energies over very short distances based on trapping slowly moving electrons in the fast moving potential wells of large amplitude plasma waves, which have relativistic phase velocities. These relativistic plasma waves, or wakefields, are the basis of table-top accelerators that have been shown to accelerate electrons to the same high energies as kilometer-length linear particle colliders operating using traditional decades-old acceleration techniques. The accelerating electrostatic fields of the relativistic plasma wave accelerators can be as large as GigaVolts/meter, and our goal was to study techniques for remotely measuring these large fields by injecting low energy probe electron beams across the plasma wave and measuring the beam’s deflection. Our method of study was via computer simulations, and these results suggested that the deflection of the probe electron beam was directly proportional to the amplitude of the plasma wave. This is the basis of a proposed diagnostic technique, and numerous studies were performed to determine the effects of changing the electron beam, plasma wave and laser beam parameters. Further simulation studies included copropagating laser beams with the relativistic plasma waves. New interesting results came out of these studies including the prediction that very small scale electron beam bunching occurs, and an anomalous line focusing of the electron beam occurs under certain conditions. These studies were summarized in the dissertation of a graduate student who obtained the Ph.D. in physics. This past research program has motivated ideas for further research to corroborate these results using particle-in-cell simulation tools which will help design a test-of-concept experiment in our laboratory and a scaled up version for testing at a major wakefield accelerator facility.

43 PARTICLE ACCELERATORS↗

Exploring the Relationship Between Octane Sensitivity and Heat-of-Vaporization

The latent heat-of-vaporization (HoV) of blends of biofuel and hydrocarbon components into gasolines has recently experienced expanded interest because of the potential for increased HoV to increase fuel knock resistance in direct-injection (DI) engines. Several studies have been conducted, with some studies identifying an additional anti-knock benefit from HoV and others failing to arrive at the same conclusion. Consideration of these studies holistically shows that they can be grouped according to the level of fuel octane sensitivity variation within their fuel matrices. When comparing fuels of different octane sensitivity significant additional anti-knock benefits associated with HoV are sometimes observed. Studies that fix the octane sensitivity find that HoV does not produce additional anti-knock benefit. New studies were performed at ORNL and NREL to further investigate the relationship between HoV and octane sensitivity. Three fuels were formulated for the ORNL study with matched RON and octane sensitivity, but with differing HoV. Experiments with these fuels in a 1.6-liter GTDI engine showed that the fuels exhibited very similar combustion phasing under knock-limited spark advance (KLSA) conditions. Fuels having a range of RON, octane sensitivity, and HoV were tested at NREL in a single-cylinder GDI engine under conditions where octane sensitivity has little effect on knock resistance. KLSA was found to be well correlated with RON. These results reinforce the concept that HoV anti-knock effects can be viewed as a contributor to octane sensitivity. From this viewpoint, HoV effects manifest themselves as increases in octane sensitivity.

30 DIRECT ENERGY CONVERSION↗

Effect of Combustion Products on Detonation Velocities

Pressure gain combustion has received increased interest in recent years because of the potential for increased thermodynamic efficiencies in propulsion and power generation applications. In pressure gain devices, mixing of fuel and oxidizer with combustion products can occur in the combustion chamber as the reactants are injected and the exhaust products are purged. Most research, however, considers only fresh reactants. This is significant because mixing can result in decreased detonation velocities. Laminar flame speed calculations show that dilution with combustion products can reduce the flame speed as a consequence of both dilution and chemistry effects. The objective of this study is to identify the effects of dilution with combustion products on detonation velocities, expanding on the results for laminar flame speed calculations. Experiments were performed using a pulse detonation engine instrumented with photodiodes to measure detonation velocities. Sonic orifices were used to meter mass flow rates of fuel (C 3 H 8 ), oxidizer (N 2 O), and combustion products (CO 2 ) or inert diluents (N 2 ). The experiments were performed such that the sensible enthalpy of the diluents was matched, thus observed velocity sensitivities indicate a chemical kinetic effect. The data indicates that non-inert diluents demonstrate a suppression in detonation velocities resulting from chemical kinetic effects. Specifically, testing with CO 2 showed a significant velocity deficit when compared to inert N 2 .

Bean, Derek↗

Ability of Particulate Matter Index to describe sooting tendency of various gasoline formulations in a stratified-charge spark-ignition engine

This study investigates the ability of Particulate Matter Index (PMI) to describe the sooting behavior of various gasoline formulations in a stratified-charge (SC) spark-ignition engine. Specifically, the engine was operated at 2000 rpm with an intake pressure of 130 kPa where soot formation is known to primarily occur in the bulk gases. Exhaust soot emissions were measured for nine test fuels at various exhaust gas recirculation levels. A comparison between measured soot levels and PMI shows that PMI is a relatively poor predictor of the sooting tendency of the tested fuels under lean SC combustion. Among the fuels, three fuels, namely the diisobutylene blend, High Olefin, and E30 fuels exhibit measured soot behavior opposite of that predicted by PMI. Optical diagnostics were utilized to further investigate the in-cylinder phenomena for these three fuels. Analysis of natural luminosity and diffused back-illumination extinction imaging suggests that fuel-induced differences in the amount of soot formed are responsible for a majority of the discrepancy in measured versus predicted sooting tendency. Fuel-induced differences in soot oxidation and spray development seem to play minor roles. Because the combustion and air-fuel mixing processes for lean SC combustion are different from conventional stoichiometric operation it is hypothesized that the PMI correlation needs to be modified to account for differences in stoichiometric air-fuel ratio and level of oxygenation between fuels. Furthermore, the role of fuel volatility in PMI possibly needs to be de-emphasized for SC operation with fuel injection into compression-heated gases.

42 ENGINEERING↗

Enterprise: a reduced-scale, flexible fuel, single-cylinder crosshead marine diesel research engine

A custom-built 1/10 scale, flexible-fuel, single-cylinder crosshead marine diesel research engine, dubbed the Enterprise, was designed and built as a platform for conducting marine lubricant research. The engine design was focused on maintaining lubricant boundary conditions similar those of full-scale engines. The Enterprise has a rated speed of 625 rpm to match the linear mean piston speed of full-scale marine engines with the reduced stroke and is capable of 18 bar brake mean effective pressure (BMEP). A high-swirl, uniflow-scavenged, two-stroke combustion system was designed to operate on either distillate diesel or heavy fuel oil (HFO). A hydraulic electronic unit injector (HEUI) fuel injection system with two injectors having customized three-hole nozzles provides fueling and is driven by the MAHLE Flexible ECU (MFE). Intake air boost is provided using a variable speed industrial air compressor along with closed-loop electronically controlled pressure regulation. Exhaust backpressure is also controlled electronically using an electro-pneumatic butterfly valve. A customized, miniaturized Hans Jensen Lubetronic cylinder lubrication system, which is capable of delivering 0.88 mg of cylinder oil per lubricator per injection to the ring pack at up to the maximum engine overspeed rating of 656 rpm, was commissioned for the Enterprise. Four independently temperature-controlled cooling zones are used to provide control of the liner temperature profile. The four coolant zones as well as the system oil and cylinder oil supplies are all controlled within 1°C for consistent thermal boundary conditions. The engine is instrumented with an in-cylinder pressure transducer as well as an inline torque transducer and numerous process temperature and pressure measurements to allow for detailed thermodynamic efficiency analysis, and experimental facilities allow for precise measurement of liner wear and characterization of deposits on the piston and rings. A radio telemetry system from IR Telemetrics is used to measure piston crown temperatures during engine operation. Engine design considerations and experimental facilities are described in detail, along with the results of a study investigating the impact of both system oil and cylinder oil viscosity on friction and efficiency. A test matrix was designed to compare base stock chemistries and viscosities for both system and cylinder oils. Experiments were completed using a fuel economy test protocol that was defined, consisting of repeating six speed/load operating conditions. An automated dynamometer control system was implemented to ensure repeatability of the test procedure. Preliminary testing indicated that the coolant, oil, and exhaust temperatures reach a steady state within 5 minutes of changing operating conditions. The system oil viscosity was found to have a significant impact on friction, with a reduction in viscosity from SAE 30 to SAE 20 resulting in over 6% reduction in friction mean effective pressure (FMEP) and over a 2% reduction in brake-specific fuel consumption (BSFC). Chemistry effects were also evident in the system and cylinder oils. Repeatability comparisons between the individual repetitions of each operating condition were used to ensure that results are statistically valid. For the initial cylinder oil viscosity experiments, observed changes were within the established/normal repeatability; ongoing experiments will determine whether improvements in boundary condition control will enable cylinder friction changes of this magnitude to be measured on this engine platform.

Kaul, Brian↗

Time-lapse seismic data inversion for estimating reservoir parameters using deep learning

Geologic carbon sequestration involves the injection of captured carbon dioxide ([Formula: see text]) into subsurface formations for long-term storage. The movement and fate of the injected [Formula: see text] plume is of great concern to regulators because monitoring helps to identify potential leakage zones and determines the possibility of safe long-term storage. To address this concern, we design a deep-learning framework for [Formula: see text] saturation monitoring to determine the geologic controls on the storage of the injected [Formula: see text]. We use different combinations of porosities and permeabilities for a given reservoir to generate saturation and velocity models. We train the deep-learning model with a few time-lapse seismic images and their corresponding changes in saturation values for a particular [Formula: see text] injection site. The deep-learning model learns the mapping from the change in the time-lapse seismic response to the change in [Formula: see text] saturation during the training phase. We then apply the trained model to data sets comprising different time-lapse seismic image slices (corresponding to different time instances) generated using different porosity and permeability distributions that are not part of the training to estimate the [Formula: see text] saturation values along with the plume extent. Our algorithm provides a deep-learning assisted framework for the direct estimation of [Formula: see text] saturation values and plume migration in heterogeneous formations using the time-lapse seismic data. Our method improves the efficiency of time-lapse inversion by streamlining the large number of intermediate steps in the conventional time-lapse inversion workflow. This method also helps to incorporate the geologic uncertainty for a given reservoir by accounting for the statistical distribution of porosity and permeability during the training phase. Tests on different examples verify the effectiveness of our approach.

Geochemistry & Geophysics↗

Electron Cyclotron Heating system status and upgrades on DIII-D

The Electron Cyclotron Heating (ECH) system on the DIII-D tokamak consists of six 110 GHz gyrotrons with corrugated coaxial 31.75 mm waveguide transmission lines and steerable launching mirrors. The system has been gradually updated, leading to increased experimental flexibility and a high system reliability of 91% in the past year. Operationally, the gyrotrons can generate up to a total of 4.8 MW of rf power for pulses up to 5 seconds in length. The maximum ECH energy injected into the DIII-D has been 16.6 MJ. The HE11 mode content is over 85% for all the lines, and the transmission coefficient is better than -1.1 dB for all the transmission lines, close to the theoretical value. A new depressed collector gyrotron was recently installed and was injecting up to 640 kW of power into the plasma during 2014-2015 tokamak operations. Three dual waveguide launchers, which can steer the RF beams ±20 degrees poloidally and toroidally, were used for real-time neoclassical tearing mode control and suppression. The launchers now have increased poloidal scanning speed and beam positioning accuracy of ~±2 mm at the plasma center. A new method of in-situ calibration of the mirror angle was used in conjunction with the upgrading of the encoders and motors for the launchers. Two more gyrotrons are expected to be installed and operational in 2015-2016. The first is a repaired 110 GHz, 1 MW gyrotron that had a gun failure after more than 11 years of operation at DIII-D. The second is a newly designed depressed collector tube in the 1.5 MW class, operating at 117.5 GHz, manufactured by Communications and Power Industries (CPI). It operates in the TE20,9 mode and has achieved 1.8 MW for short pulses during factory testing. Furthermore, this gyrotron is undergoing rework to address a high voltage standoff problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Immediate effects of acute Mars mission equivalent doses of SEP and GCR radiation on the murine gastrointestinal system-protective effects of curcumin-loaded nanolipoprotein particles (cNLPs)

Introduction: Missions beyond low Earth orbit (LEO) will expose astronauts to ionizing radiation (IR) in the form of solar energetic particles (SEP) and galactic cosmic rays (GCR) including high atomic number and energy (HZE) nuclei. The gastrointestinal (GI) system is documented to be highly radiosensitive with even relatively low dose IR exposures capable of inducing mucosal lesions and disrupting epithelial barrier function. IR is also an established risk factor for colorectal cancer (CRC) with several studies examining long-term GI effects of SEP/GCR exposure using tumor-prone APC mouse models. Studies of acute short-term effects of modeled space radiation exposures in wildtype mouse models are more limited and necessary to better define charged particle- induced GI pathologies and test novel medical countermeasures (MCMs) to promote astronaut safety. Methods: In this study, we performed ground-based studies where male and female C57BL/6J mice were exposed to γ-rays, 50 MeV protons, or 1 GeV/n Fe-56 ions at the NASA Space Radiation Laboratory (NSRL) with histology and immunohistochemistry endpoints measured in the first 24 h post-irradiation to define immediate SEP/GCR-induced GI alterations. Results: Our data show that unlike matched γ-ray controls, acute exposures to protons and iron ions disrupts intestinal function and induces mucosal lesions, vascular congestion, epithelial barrier breakdown, and marked enlargement of mucosa-associated lymphoid tissue. We also measured kinetics of DNA double-strand break (DSB) repair using gamma-H2AX- specific antibodies and apoptosis via TUNEL labeling, noting the induction and disappearance of extranuclear cytoplasmic DNA marked by gamma-H2AX only in the charged particle-irradiated samples. We show that 18 h pre-treatment with curcumin-loaded nanolipoprotein particles (cNLPs) delivered via IV injection reduces DSB-associated foci levels and apoptosis and restore crypt villi lengths. Discussion: These data improve our understanding of physiological alterations in the GI tract immediately following exposures to modeled space radiations and demonstrates effectiveness of a promising space radiation MCM.

59 BASIC BIOLOGICAL SCIENCES↗

Hard X-Ray Emission from the Eastern Jet of SS 433 Powering the W50 “Manatee” Nebula: Evidence for Particle Reacceleration

We present a broadband X-ray study of W50 (the “Manatee” nebula), the complex region powered by the microquasar SS 433, that provides a test bed for several important astrophysical processes. The W50 nebula, a Galactic PeVatron candidate, is classified as a supernova remnant but has an unusual double-lobed morphology likely associated with the jets from SS 433. Using NuSTAR, XMM-Newton, and Chandra observations of the inner eastern lobe of W50, we have detected hard nonthermal X-ray emission up to ~30 keV, originating from a few-arcminute-sized knotty region (“Head”) located ≲18' (29 pc for a distance of 5.5 kpc) east of SS 433, and constrained its photon index to 1.58 ± 0.05 (0.5–30 keV band). The index gradually steepens eastward out to the radio “ear” where thermal soft X-ray emission with a temperature kT ~ 0.2 keV dominates. The hard X-ray knots mark the location of acceleration sites within the jet and require an equipartition magnetic field of the order of ≳12 μG. The unusually hard spectral index from the “Head” region challenges classical particle acceleration processes and points to particle injection and reacceleration in the subrelativistic SS 433 jet, as seen in blazars and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗

Pore architecture controls on mineral reactivity

Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.

58 GEOSCIENCES↗

The Cost and Benefit of Enhancing Cybersecurity for Hybrid AC/DC Grids

As critical interfaces of AC grids and DC grids inside a hybrid AC/DC grid, the voltage-sourced-converter (VSC) has been demonstrated to be vulnerable to false data injection (FDI) cyber-attacks. As a result, the cyber-attack-induced AC grid frequency deviations and DC grid voltage deviations threaten the secure operation. Here, to enhance cybersecurity in a not only feasible but also cost-effective manner, this paper proposes a cost-benefit-based cyber-defense strategy for a hybrid AC/DC grid. First, this paper establishes a spatial-temporal dual cyber-attack evaluation model, in which the cyber-attack-induced frequency and voltage deviations are modelled in both a spatially and temporally dual manner. Then, the proposed cost-benefit-based cyber-defense strategy is modelled as a VSC commitment problem to achieve the trade-off between maximizing the cyber-defense benefits and minimizing the cyber-defense costs. The VSC commitment problem is then mathematically convexified into a mixed-integer second-order cone programming (MISOCP) problem, which could be efficiently solved in an event-triggered manner against unfolding cyber-attack events. Simulation results on a test hybrid AC/DC grid verified the feasibility and the cost-effectiveness of the proposed cyber-defense strategy.

cost-benefit analysis↗

Final Report: Optimal Model Complexity in Geological Carbon Sequestration: A Response Surface Uncertainty Analysis

The critical component of a risk assessment study in evaluating GCS is an analysis of uncertainty in CO2 modeling. In such analyses, direct numerical simulation of CO2 flow and leakage requires many time-consuming model runs. Alternatively, analytical methods have been developed which allow fast and efficient estimation of CO2 storage and leakage, although restrictive assumptions on formation rock and fluid properties are employed. In this study, an intermediate approach is proposed based on the Design of Experiment and Response Surface methodology, which consists of using a limited number of numerical simulations to estimate a prediction outcome as a combination of the most influential uncertain site properties. The methodology can be implemented within a Monte Carlo framework to efficiently assess parameter and prediction uncertainty while honoring the accuracy of numerical simulations. The choice of the uncertain properties is flexible and can include geologic parameters that influence reservoir heterogeneity, engineering parameters that influence gas trapping and migration, and reactive parameters that influence the extent of fluid/rock reactions. The method was tested and verified on modeling long-term CO2 flow, non-isothermal heat transport, and CO2 dissolution storage by coupling two-phase flow with explicit miscibility calculation using an accurate equation of state that gives rise to convective mixing of formation brine variably saturated with CO2. All simulations were performed using three-dimensional high-resolution models including a target deep saline aquifer, overlying caprock, and a shallow aquifer. To evaluate the uncertainty in representing reservoir permeability, sediment hierarchy of a heterogeneous digital stratigraphy was mapped to create multiple irregularly shape stratigraphic models of decreasing geologic resolutions: heterogeneous (reference), lithofacies, depositional environment, and a (homogeneous) geologic formation. To ensure model equivalency, all the stratigraphic models were successfully upscaled from the reference heterogeneous model for bulk flow and transport predictions (Zhang & Zhang, 2015). GCS simulation was then simulated with all models, yielding insights into the level of parameterization complexity that is needed for the accurate simulation of reservoir pore pressure, CO2 storage, leakage, footprint, and dissolution over both short (i.e., injection) and longer (monitoring) time scales. Important uncertainty parameters that impact these key performance metrics were identified for the stratigraphic models as well as for the heterogeneous model, leading to the development of reduced/simplified models at lower characterization cost that can be used for the reservoir uncertainty analysis. All the CO2 modeling was conducted using PFLOTRAN – a massively parallel, multiphase, multi-component, and reactive transport simulator developed by a multi-laboratory DOE/SciDAC (Scientific Discovery through Advanced Computing) project (Zhang et al., 2017, in review). Within the uncertainty analysis framework, increasing reservoir depth were investigated to explore its effect on the uncertainty outcomes and the potential for developing gravity-stable injection with increased storage security (Dai et al., 20126; Dai et al., 2017, in review). Finally, to accurately model CO2 fluid-rock reactions and resulting long-term storage as secondary carbonate minerals, a modified kinetic rate law for general mineral dissolution and precipitation was proposed and verified that is invariant to a scale transformation of the mineral formula weight. This new formulation will lead to more accurate assessment of mineral storage over geologic time scales (Lichtner, 2016).

58 GEOSCIENCES↗

Modeling and optimization of steady flow of natural gas and hydrogen mixtures in pipeline networks

Here, we extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of energy systems on fossil fuels. Flow in a pipe is related to endpoint pressures by a basic Weymouth equation model, with an ideal gas equation of state, where the wave speed depends on the hydrogen concentration. At vertices, in addition to mass balance, we also consider mixing of incoming flows of varying hydrogen concentrations. The problems of interest are the heterogeneous gas flow simulation (HGFS), which determines system pressures and flows given fixed boundary conditions and compressor settings, as well as the heterogeneous gas flow optimization (HGFO), which extremizes an objective by determining optimal boundary conditions and compressor settings. We examine conditions for uniqueness of solutions to the HGFS, as well as compare and contrast mixed-integer and continuous nonlinear programming formulations for the HGFO. We develop computational methods to solve both problems, and examine their performance using four test networks of increasing complexity.

08 HYDROGEN↗