Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “injection rate”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A reduced-turbulence regime in the Large Helical Device upon injection of low-Z materials powders

Abstract Recently an improved confinement regime, characterized by reduced turbulent fluctuations has been observed in the Large Helical Device upon the injection of boron powder into the plasma (Nespoli et al 2022 Nat. Phys. 18 350–56). In this article, we report in more detail the experimental observations of increased plasma temperature and the decrease of turbulent fluctuations across the plasma cross section, on an extended database. In particular, we compare powders of different materials (B, C, BN), finding similar temperature improvement and turbulence response for the three cases. Modeling of the powder penetration into the plasma and of neoclassical electric field and fluxes support the interpretation of the experimental results. Additionally, we report evidence of the temperature improvement increasing with powder injection rates and decreasing for both increasing density and heating power. Though, plasma turbulence response varies depending on the initial conditions of the plasma, making it difficult to draw an inclusive description of the phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Power and isotope effects in the ITER baseline scenario with tungsten and tungsten-equivalent radiators in DIII-D

Abstract Experiments in DIII-D document the ITER Baseline Scenario (IBS) at q 95 ∼ 3 and P IN / P LH ∼ 1–2, in both deuterium and hydrogen utilizing Kr and Xe as Tungsten-equivalent radiators. The power threshold for H-mode operation ( P LH ) was determined experimentally without added impurities and found to be about a factor of two higher than the scaling law. In recent IBS experiments in deuterium, intrinsic levels of metals such as Tungsten (W) or molybdenum and inconel are present that reduce the pedestal pressure by 20%–25%. A complete radiative collapse of deuterium IBS plasmas occurs at W core concentrations C W = 10 −5 . Simulations show that for core temperatures expected for ITER, the plasmas would not have a radiative collapse at C W = 1 × 10 −5 , moreover Q = 8–10 would still be achieved for C W up to 3 × 10 −5 . In contrast to deuterium, the IBS in hydrogen is not affected by intrinsic high-Z impurities, indicating that hydrogen H-modes in ITER may not inform the D-T phase with respect to W accumulation and discharge survival. Compared to deuterium, the pedestal pressure in hydrogen is ∼25% lower, with much higher ELM frequency of 150 Hz, decreasing with input power. Krypton was injected in a matrix scan of input power and impurity flow in IBS hydrogen discharges. Krypton impurity density profiles in hydrogen are similar to deuterium plasmas, but at Kr flows that are 2–3 times higher for the same input power. Krypton is transported into the core and affects the whole radius; at the highest injection rates a radiative collapse occurs at core radiation fractions of 0.3–0.35, consistent with the expected maximum W radiation fraction for ITER core plasmas. Comparing the results with previous International Tokamak Physics Activity database studies of the IBS confirms that at higher radiation fraction due to high-Z impurities, a drop in H 98 of >10% is observed. On the other hand, the results using Kr as a W-equivalent radiator indicate that metal (W) devices at lower core temperatures than ITER may provide overly pessimistic performance extrapolations to ITER for deuterium-tritium operation. The new DIII-D results support a more attractive option for the ITER Research Plan with a short hydrogen phase for system commissioning, transitioning to deuterium operations as soon as possible to provide relevant conditions for deuterium-tritium operations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kicking the can down the road: understanding the effects of delaying the deployment of stratospheric aerosol injection

Abstract Climate change is a prevalent threat, and it is unlikely that current mitigation efforts will be enough to avoid unwanted impacts. One potential option to reduce climate change impacts is the use of stratospheric aerosol injection (SAI). Even if SAI is ultimately deployed, it might be initiated only after some temperature target is exceeded. The consequences of such a delay are assessed herein. This study compares two cases, with the same target global mean temperature of ∼1.5° C above preindustrial, but start dates of 2035 or a ‘delayed’ start in 2045. We make use of simulations in the Community Earth System Model version 2 with the Whole Atmosphere Coupled Chemistry Model version 6 (CESM2-WACCM6), using SAI under the SSP2-4.5 emissions pathway. We find that delaying the start of deployment (relative to the target temperature) necessitates lower net radiative forcing (−30%) and thus larger sulfur dioxide injection rates (+20%), even after surface temperatures converge, to compensate for the extra energy absorbed by the Earth system. Southern hemisphere ozone is higher from 2035 to 2050 in the delayed start scenario, but converges to the same value later in the century. However, many of the surface climate differences between the 2035 and 2045 start simulations appear to be small during the 10–25 years following the delayed SAI start, although longer simulations would be needed to assess any longer-term impacts in this model. In addition, irreversibilities and tipping points that might be triggered during the period of increased warming may not be adequately represented in the model but could change this conclusion in the real world.

Brody, Ezra (ORCID:000900030008681X)↗
DOE OSTI2024

Fireball antinucleosynthesis

The tentative identification of approximately ten relativistic antihelium ( He ¯ ) cosmic-ray events at AMS-02 would, if confirmed, challenge our understanding of the astrophysical synthesis of heavy antinuclei. We propose a novel scenario for the enhanced production of such antinuclei that is triggered by isolated, catastrophic injections of large quantities of energetic Standard Model (SM) antiquarks in our galaxy by physics beyond the Standard Model (BSM). We demonstrate that SM antinucleosynthetic processes that occur in the resulting rapidly expanding, thermalized fireballs of SM plasma can, for a reasonable range of parameters, produce the reported tentative ∼ 2 : 1 ratio of He ¯ 3 to He ¯ 4 events at AMS-02, as well as their relativistic boosts. Moreover, we show that this can be achieved without violating antideuterium or antiproton flux constraints for the appropriate antihelium fluxes. A plausible BSM paradigm for the catastrophic injections is the collision of macroscopic composite dark-matter objects carrying large net antibaryon number. Such a scenario would require these objects to be cosmologically stable, but to destabilize upon collision, promptly releasing a fraction of their mass energy into SM antiparticles within a tiny volume. We show that, in principle, the injection rate needed to attain the necessary antihelium fluxes and the energetic conditions required to seed the fireballs appear possible to obtain in such a paradigm. We leave open the question of constructing a BSM particle physics model to realize this, but we suggest two concrete scenarios as promising targets for further investigation. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Unstable Invasion of Sedimenting Granular Suspensions

Here we investigate the development of mobility inversion and fingering when a granular suspension is injected radially between horizontal parallel plates of a cell filled with a miscible fluid. While the suspension spreads uniformly when the suspension and the displaced fluid densities are exactly matched, even a small density difference is found to result in a dense granular front which develops fingers with angular spacing that increase with granular volume fraction and decrease with injection rate. We show that the timescale over which the instability develops is given by the volume fraction dependent settling timescale of the grains in the cell. We then show that the mobility inversion and the nonequilibrium Korteweg surface tension due to granular volume fraction gradients determine the number of fingers at the onset of the instability in these miscible suspensions.

02 PETROLEUM↗
DOE OSTI2021

Stabilization of Alfvén Eigenmodes in DIII-D via Controlled Energetic Ion Density Ramp and Validation of Theory and Simulations

In this work, fast-ion driven Alfvén waves with frequency close to the ion cyclotron frequency (ƒ = 0:58ƒ ci ) excited by energetic ions from a neutral beam are stabilized via a controlled energetic ion density ramp for the first time in a fusion research plasma. The scaling of wave amplitude with injection rate is consistent with theory for single mode collisional saturation near marginal stability. The wave is identified as a shear-polarized global Alfvén eigenmode excited by Doppler-shifted cyclotron resonance with fast-ions with sub-Alfvénic energetic ions, a first in fusion research plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

Herein we formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results [1]-[4], indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Dual Nature of Entrainment-Mixing Signatures Revealed through Large-Eddy Simulations of a Convection-Cloud Chamber

Abstract Entrainment of subsaturated air into a cloud can influence its optical and microphysical properties in various ways, depending on the droplet evaporation and turbulent mixing time scales. Previous experiments in the Pi convection-cloud chamber have revealed that, given a fixed entrained air property, the mixing of entrained subsaturated air results in complete evaporation of some cloud droplets, with the rest remaining unchanged. This is a signature of inhomogeneous mixing. While comparing the results of entrainment with varying air properties, the mixing signature appears as if the subsaturated air is well mixed with the cloud to evenly reduce the droplets’ size. In other words, taken together, the experiments appear to have the signature of homogeneous mixing. To explore these results in a greater depth, we conduct large-eddy simulations combined with a bin microphysics scheme. Our results reproduce the similar signatures of inhomogeneous and homogeneous mixing, implying that LES can resolve the inhomogeneous mixing when the grid spacing is smaller than the entrained air parcel. Additionally, we observe that increasing the aerosol injection rate enhances the signature of inhomogeneous mixing, while coarser grid spacing diminishes it. Finally, the change in wall fluxes in response to various entrained air properties confirms that the homogeneous signature seen in the analysis of an ensemble of simulations is the result of various equilibrium states. This further strengthens the suggestion that the homogeneous mixing signature found in aircraft observations near the cloud top may result from combining entrainment events of different intensities, possibly caused by various-sized eddies. Significance Statement Large-eddy simulation and size-resolved microphysics can resolve time scales for turbulent mixing and evaporation and, therefore, are well suited for reproducing, extending, and interpreting the entrainment experiment in the Pi convection-cloud chamber. Our simulation results confirm (i) the inhomogeneous mixing signature for an individual entrainment event and (ii) the appearance of homogeneous mixing in an ensemble of entrainment episodes. Furthermore, we demonstrate that the inhomogeneous mixing signature is more pronounced in a polluted cloud, but coarser grid spacing in simulations may compromise the accuracy of this signature. Last, the homogeneous mixing signature results from various equilibrium states established for different entrainment intensities and adjusted wall fluxes, which are challenging to measure experimentally but can be easily analyzed in the simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

EGS Collab Experiment 1 Stimulation Data

Stimulation data from Experiment 1 of EGS Collab, which occurred on the 4850 ft level of the Sanford Underground Research Facility (SURF). A detailed description of the stimulation data is provided in the StimulationDataNotes.docx and is also available on the EGS Collab Wiki. A Meta Data Cheat Sheet, which describes all of the channels in the Raw CSV files, is available as well. Note that this cheat sheet is a comprehensive meta data descriptor and channels were added as the experiment evolved. This means that some columns may not be populated in early data. Additionally, we have included the chat logs from these experiments. The experiments were broadcast over teleconferencing software and real-time data displays were available to remote observers. The logs contain important observations from those personnel performing the experiment and the remote contributors. Finally, we have included summary and individual plots of all of the data for the user to compare to.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

EGS Collab Experiment 1: Circulation Testing Processed data

This submission includes processed and reduced data for circulation testing that was conducted at the 164' fracture on the 4850 ft level of the Sanford Underground Research Facility. The circulation tests were done to test the flow through the 164' fracture in the EGS Collab Experiment 1 testbed on the 4850 ft level of the Sanford Underground Research Facility in order to validate computer models concerned with flow processes and heat exchange processes in EGS. This data was processed from raw data that is available in the GDR submission linked below: "EGS Collab Circulation Testing Raw Data and Documentation." Python scripts used to process and reduce the data are included. The scripts also serve as "documentation" as well as examples for data processing. More information about the test itself as well as the rationale and process of data processing is available at wiki page linked below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Microwave-Assisted Hydrogen Generation from Hydrocarbon-Bearing Reservoir Rocks: Stage-Dependent Thermal Runaway and In-Situ Carbonate Engineering

Microwave-assisted hydrogen generation from hydrocarbon-bearing reservoir rocks is strongly influenced by mineralogy, methane activation, carbonate reactions, and thermal runaway behavior. This study investigates a new approach in which carbonate phases are generated in-situ through the reaction of internally produced CO2 with Ca(OH)2 under microwave heating conditions. The objective is to evaluate how rock mineralogy, methane injection, and Ca(OH)2 addition influence hydrogen generation, carbon redistribution, and stage-dependent reaction pathways during microwave exposure. Microwave heating experiments were conducted using Permian Basin reservoir rocks under three experimental conditions: rock-only experiments under Ar atmosphere, CH4–Ar experiments without additive, and CH4–Ar experiments containing 5 wt% Ca(OH)2. Methane-assisted experiments were performed under continuous injection of 30 standard cubic centimeters per minute (sccm) CH4 and 30 sccm Ar. Based on thermal runaway behavior, each experiment was divided into three operational stages: Before Thermal Runaway (BR), After Thermal Runaway–Decrease in Microwave Power (ARD), and After Thermal Runaway– Increase in Microwave Power (ARI). Temperature and gas composition were continuously monitored throughout the experiments. The rock-only experiments demonstrated that hydrogen generation can occur intrinsically from hydrocarbon-bearing rocks under microwave heating, even without externally injected methane. However, hydrogen production did not correlate solely with kerogen content, indicating that mineralogy strongly influences hydrogen-generation pathways. Correlation analyses suggested that kerogen decomposition initially generated CH4, CO, and CO2, followed by secondary hydrocarbon reactions associated with H2 and C2 hydrocarbon formation. Methane-assisted experiments substantially increased hydrogen production; however, identical methane injection rates produced significantly different hydrogen yields among the rock samples, confirming that mineralogical composition controls methaneconversion behavior under microwave heating conditions. The addition of Ca(OH)2 significantly altered carbon evolution behavior in a stage-dependent manner. During the BR stage, Ca(OH)2 reduced gas-phase CO2 production, particularly in carbonate-rich rocks, indicating favorable conditions for in-situ carbonation and carbonate deposition prior to extensive thermal decomposition. The suppression of CO2 during BR became more pronounced with increasing carbonate content of the rock system. After thermal runaway, carbonate-containing systems exhibited enhanced hydrogen generation behavior, suggesting that carbonatederived mineral transformations and carbonate-mediated reactions contribute to high-temperature hydrogen-generation pathways. In carbonate-poor rocks, Ca(OH)2 addition enabled simultaneous URTeC 4493775 2 enhancement of hydrogen production and partial suppression of CO2 release during the post-runaway stages. Overall, the results demonstrate that microwave-assisted hydrogen generation is governed by dynamically evolving interactions among kerogen decomposition, methane activation, mineral transformations, carbonate formation/decomposition, and thermal runaway behavior. This work introduces in-situ carbonate engineering with Ca(OH)2 as a strategy for coupling hydrogen generation with partial insitu carbon management under microwave heating conditions.

03 NATURAL GAS↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2022

Optimization of Water-Alternating-CO2 Injection Field Operations Using a Machine-Learning-Assisted Workflow

Summary This paper will present a robust workflow to address multiobjective optimization (MOO) of carbon dioxide (CO2)-enhanced oil recovery (EOR)-sequestration projects with a large number of operational control parameters. Farnsworth unit (FWU) field, a mature oil reservoir undergoing CO2 alternating water injection (CO2-WAG) EOR, will be used as a field case to validate the proposed optimization protocol. The expected outcome of this work would be a repository of Pareto-optimal solutions of multiple objective functions, including oil recovery, carbon storage volume, and project economics. FWU’s numerical model is used to demonstrate the proposed optimization workflow. Because using MOO requires computationally intensive procedures, machine-learning-based proxies are introduced to substitute for the high-fidelity model, thus reducing the total computation overhead. The vector machine regression combined with the Gaussian kernel (Gaussian-SVR) is used to construct proxies. An iterative self-adjusting process prepares the training knowledge base to develop robust proxies and minimizes computational time. The proxies’ hyperparameters will be optimally designed using Bayesian optimization to achieve better generalization performance. Trained proxies will be coupled with multiobjective particle swarm Optimization (MOPSO) protocol to construct the Pareto-front solution repository. The outcomes of this workflow will be a repository containing Pareto-optimal solutions of multiple objectives considered in the CO2-WAG project. The proposed optimization workflow will be compared with another established methodology using a multilayer neural network (MLNN) to validate its feasibility in handling MOO with a large number of parameters to control. Optimization parameters used include operational variables that might be used to control the CO2-WAG process, such as the duration of the water/gas injection period, producer bottomhole pressure (BHP) control, and water injection rate of each well included in the numerical model. It is proved that the workflow coupling Gaussian-SVR proxies and the iterative self-adjusting protocol is more computationally efficient. The MOO process is made more rapid by squeezing the size of the required training knowledge base while maintaining the high accuracy of the optimized results. The outcomes of the optimization study show promising results in successfully establishing the solution repository considering multiple objective functions. Results are also verified by validating the Pareto fronts with simulation results using obtained optimized control parameters. The outcome from this work could provide field operators an opportunity to design a CO2-WAG project using as many inputs as possible from the reservoir models. The proposed work introduces a novel concept that couples Gaussian-SVR proxies with a self-adjusting protocol to increase the computational efficiency of the proposed workflow and to guarantee the high accuracy of the obtained optimized results. More importantly, the workflow can optimize a large number of control parameters used in a complex CO2-WAG process, which greatly extends its utility in solving large-scale MOO problems in various projects with similar desired outcomes.

Energy & Fuels↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2020

Distributed Temperature Sensing (DTS) to Monitor CO 2 Migration in an Enhanced Oil Recovery Field in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. MRCSP conducted analysis of CO 2 flow patterns using distributed temperature sensing (DTS) data acquired from the Chester 16 pinnacle reef located in Otsego County, Michigan. The reef has one injection well Chester 6-16 and one monitoring well Chester 8-16 where DTS is installed. DTS provides a means to track temperature within the wellbores in real time and characterize the migration patterns of CO 2 in the subsurface when DTS data is combined with other operational data such as the injection rates, bottomhole pressures and temperature. The primary objective of this study was to determine the distribution of CO 2 inflow zones in the A1 Carbonate and Brown Niagaran formations.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

SECARB-USA: Preliminary Regional Assessment of Potential Storage Complexes, Including Capacities and Costs, Based on Application of NRAP and SCO2T

This subtask focused on producing estimates of CO 2 injection rates, storage capacity, and sequestrations costs across the entire Southeast Regional CO 2 Utilization and Storage Acceleration Partnership (SECARB-USA) region. These estimates will then be used in subsequent tasks that rely on understanding CO 2 storage and sequestration cost estimates. The data are generated and reported using 10x10 km (100 km 2 ) grid cells across based on the NATCARB dataset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Plasma-assisted surface modification of biopolymers

The interaction of non-thermal plasmas with surfaces is an important problem of interest in various applications. This work focuses on two main aspects of the problem - simulation and experiments. The development of a low-temperature plasma modeling is described along with experiments focusing on the interaction of a low-temperature plasma with liquid water. The simulation framework is capable of handling arbitrary geometries, chemistry and a fully-coupled boundary condition for plasma-dielectric interface. Representative examples for the verification and validation of the framework are provided. The experimental work deals with the study of low-temperature plasmas with water with the nitrate injection rate compared for three different configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗