Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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 451 records · Page 25

Subsurface Uranium Fate and Transport: Integrated Experiments and Modeling of Coupled Biogeochemical Mechanisms of Nanocrystalline Uraninite Oxidation by Fe(III)-(hydr)oxides - Project Final Report

Subsurface bacteria including sulfate reducing bacteria (SRB) reduce soluble U(VI) to insoluble U(IV) with subsequent precipitation of UO 2 . We have shown that SRB reduce U(VI) to nanometer-sized UO 2 particles (1-5 nm) which are both intra- and extracellular, with UO 2 inside the cell likely physically shielded from subsequent oxidation processes. We evaluated the UO 2 nanoparticles produced by Desulfovibrio desulfuricans G20 under growth and non-growth conditions in the presence of lactate or pyruvate and sulfate, thiosulfate, or fumarate, using ultrafiltration and HR-TEM. Results showed that a significant mass fraction of bioreduced U (35-60%) existed as a mobile phase when the initial concentration of U(VI) was 160 µM. Further experiments with different initial U(VI) concentrations (25 - 900 M) in MTM with PIPES or bicarbonate buffers indicated that aggregation of uraninite depended on the initial concentrations of U(VI) and type of buffer. It is known that under some conditions SRB-mediated UO 2 nanocrystals can be reoxidized (and thus remobilized) by Fe(III)-(hydr)oxides, common constituents of soils and sediments. To elucidate the mechanism of UO 2 reoxidation by Fe(III) (hydr)oxides, we studied the impact of Fe and U chelating compounds (citrate, NTA, and EDTA) on reoxidation rates. Experiments were conducted in anaerobic batch systems in PIPES buffer. Results showed EDTA significantly accelerated UO 2 reoxidation with an initial rate of 9.5 M day-1 for ferrihydrite. In all cases, bicarbonate increased the rate and extent of UO 2 reoxidation with ferrihydrite. The highest rate of UO 2 reoxidation occurred when the chelator promoted UO 2 and Fe(III) (hydr)oxide dissolution as demonstrated with EDTA. When UO 2 dissolution did not occur, UO 2 reoxidation likely proceeded through an aqueous Fe(III) intermediate as observed for both NTA and citrate. To complement to these laboratory studies, we collected U-bearing samples from a surface seep at the Rifle field site and have measured elevated U concentrations in oxic iron-rich sediments. To translate experimental results into numerical analysis of U fate and transport, a reaction network was developed based on Sani et al. (2004) to simulate U(VI) bioreduction with concomitant UO 2 reoxidation in the presence of hematite or ferrihydrite. The reduction phase considers SRB reduction (using lactate) with the reductive dissolution of Fe(III) solids, which is set to be microbially mediated as well as abiotically driven by sulfide. Model results show the oxidation of HS– by Fe(III) directly competes with UO 2 reoxidation as Fe(III) oxidizes HS– preferentially over UO 2 . The majority of Fe reduction is predicted to be abiotic, with ferrihydrite becoming fully consumed by reaction with sulfide. Predicted total dissolved carbonate concentrations from the degradation of lactate are elevated (log(pCO 2 ) ~ –1) and, in the hematite system, yield close to two orders-of-magnitude higher U(VI) concentrations than under initial carbonate concentrations of 3 mM. Modeling of U(VI) bioreduction with concomitant reoxidation of UO 2 in the presence of ferrihydrite was also extended to a two-dimensional field-scale groundwater flow and biogeochemically reactive transport model for the South Oyster site in eastern Virginia. This model was developed to simulate the field-scale immobilization and subsequent reoxidation of U by a biologically mediated reaction network.

54 ENVIRONMENTAL SCIENCES↗

2018 FEARCE Development: A Robust and Accurate Engine Modeling Software

Research and development of Fast, Easy, Accurate and Robust Continuum Engineering (FEARCE, formerly KIVA-hpFE) for turbulent reactive and multiphase flow, particularly as related to engine modeling, is relevant to the DOE Vehicle Technologies Office efforts at addressing national energy security. Less dependence on petroleum products leads to greater energy security. By Environmental Protection Agency standards, some vehicles are now reaching 42–50 mpg mark. These are conventional gasoline engines. With continued investment and research into new technical innovations, the potential exists to save more than 4 million barrels of oil per day or approximately 200 to 400 dollars million per day. This would be a significant decrease in emissions and use of petroleum and a very large stimulus to the U.S. economy.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Innovative Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Numerical simulation of non-isothermal multiphase porous flow combined with reactive transport is used in a wide range of applications, which include nuclear waste repositories, enhanced recovery of petroleum reservoirs, contaminant remediation, geothermal engineering, and carbon sequestration. Understanding and predicting underground phenomena can have enormous impact on how to deal with world issues like climate change, clean water, and renewable energy. The main motivation for this proposal arises from safety assessment of future nuclear waste repositories using the U.S. Department of Energy (DOE) Geologic Disposal Safety Assessment (GDSA) Framework, and performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the nation's only active nuclear waste repository in Carlsbad, New Mexico. I am a technical staff member at Sandia National Laboratories involved in both projects.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Supersonic flameholding by attached oblique shock waves

The object of this investigation was to explore various regimes of shock-induced combustion that may occur as a result of supersonic flow of a reactive gas mixture past a simple compression ramp or wedge. Attention was confined to combinations of wedge angles, fuel-air mixture ratios and freestream velocities that permit the resulting oblique shock or oblique detonation wave to remain attached. Computational fluid dynamic (CFD) simulations, both inviscid and viscous, were performed and analyzed in order to investigate two modes of compression ramp flameholding; namely, complete oblique detonation wave (ODW) and shock-induced combustion. In addition, the possibility of flame persistence or hysteresis was explored. It was preliminarily concluded that inclusion of viscous effects may be essential to accurately predict the occurrence of the various flameholding modes.

Fort, James A.↗

Water Ice and Life's Roots in Space

Nearly three decades ago as Voyager 2 spacecraft raced out of the Solar System. NASA engineers turned its camera arm around (at the request of the American astronomer Carl Sagan) to take a parting snapshot of Earth. Earth's image was a single pale blue pixel, its color caused by the Rayleigh scattering of sunlight in the water of our oceans. Earth is a water planet, and this is the color of life. No matter how far we travel on our planet, no matter how high or deep, if we find liquid water, we find some form of life that manages to survive there. And yet there is a cruel irony. Water in its solid crystalline form is hostile to life. Organisms can roost in geysers, wallow in brine and gulp down acid, but they cowered from ice. The rigid ordering of water molecules in ice crystals expels impurities and tears organic tissue beyond repair. In fact, about the only good thing you can say about ice is that it gets out of the way: Its low density ensures that it floats and leaves the water dwelling creatures in peace. Recent discoveries have caused us to rethink this basic premise. New lines of evidence both observational and experimental - suggest that prebiotic organic compounds are not only comfortable in, but in fact had their origin in a peculiar form of solid water ice that is ubiquitous in interstellar space, but completely absent from Earth. Only recently have we been able to create even submicroscopic quantities of this ice in terrestrial laboratories, yet it constitutes the most abundant form of water in the universe. Interstellar ice is a far cry from the ice we are so familiar with on Earth. This interstellar ice has no crystalline structure, and despite the fact that its temperature is a scant few degrees above absolute zero (where all molecular motion ceases), it is highly reactive and can flow like water when exposed to radiation. It is in fact this ice's similarity to liquid water that allows it to participate in the creation of the very first organic compounds.

Blake, David↗

Using Bayesian Networks for Sensitivity Analysis of Complex Biogeochemical Models

Sensitivity analysis is a vital tool in numerical modeling to identify important parameters and processes that contribute to the overall uncertainty in model outputs. We developed a new sensitivity analysis method to quantify the relative importance of uncertain model processes that contain multiple uncertain parameters. The method is based on the concepts of Bayesian networks (BNs) to account for complex hierarchical uncertainty structure of a model system. We derived a new set of sensitivity indices using the methodology of variance-based global sensitivity analysis with the Bayesian inference. The framework is capable of representing the detailed uncertainty information of a complex model system using BNs and affords flexible grouping of different uncertain inputs given their characteristics and dependency structures. We have implemented the method on a real-world biogeochemical model at the groundwater-surface water interface within the Hanford Site’s 300 Area. The uncertainty sources of the model were first grouped into forcing scenario and three different processes based on our understanding of the complex system. The sensitivity analysis results indicate that both the reactive transport and groundwater flow processes are important sources of uncertainty for carbon-consumption predictions. Within the groundwater flow process, the structure of geological formations is more important than the permeability heterogeneity within a given geological formation. Our new sensitivity analysis framework based on BNs offers substantial flexibility for investigating the importance of combinations of interacting uncertainty sources in a hierarchical order, and it is expected to be applicable to a wide range of multi-physics models for complex systems.

Dai, Heng↗

Analysis of the Reactive CO 2 Surface Flux in Electrocatalytic Aqueous Flow Reactors

We study how the mass transfer of a chemical species in shear flow is suppressed by the production of a second species that partially reacts away the first species. The second species is produced at the surface with the first species as a reactant. Our work is directly motivated by electrochemical CO 2 reduction in aqueous flow reactors, where OH– molecules generated by the CO evolution reaction react away CO 2 molecules, ultimately inhibiting the mass transfer of CO 2 to the cathode surface. We derive a simple approximation of the Sherwood number, a dimensionless measure of the mass flux into the surface, as a function of the Péclet number, surface Damköhler number, and bulk Damköhler number, and we compare our approximation to numerical solutions of the governing equations. We find that in the diffusion-limited regime, the Sherwood number is well described by the classical Graetz-Lévêque result reduced by a reaction factor due to the competing bulk reaction; with the stoichiometry relevant to electrochemical CO 2 reduction, this reaction factor is 1/2. While the model problem we solve provides insight into how OH– production affects CO 2 mass transfer, experimental systems often involve more complex chemistry. We thus also show how a common buffered electrolyte, KHCO 3 , affects our results.

30 DIRECT ENERGY CONVERSION↗

Predictive performance modeling framework for a novel enclosed particle receiver configuration and application for thermochemical energy storage

Concentrating solar power (CSP) plants can provide dispatchable power with a thermal energy storage capability for increased renewable-energy grid penetration. Particle-based CSP systems permit higher temperatures, and thus, potentially higher solar-to-electric efficiency than state-of-the-art molten-salt heat-transfer systems. This paper describes a detailed numerical analysis framework for estimating the performance of a novel, geometrically complex, enclosed particle receiver design. The receiver configuration uses arrays of small tubular absorbers to collect and subsequently transfer solar energy to a flowing particulate medium. The enclosed nature of the receiver design renders it amenable to either an inert heat-transfer medium, or a reactive heat-transfer medium that requires a controllable ambient environment. The numerical analysis framework described in this study is demonstrated for the case of thermal reduction of CaCr 0.1 Mn 0.9 O 3-$\delta$ for thermochemical energy storage. The modeling strategy consists of Monte Carlo ray tracing for absorbed solar-energy distributions from a surround heliostat field, computational fluid dynamics modeling of small-scale local tubular arrays, surrogate response surfaces that approximately capture simulated tubular array performance, a quasi-two-dimensional reduced-order description of counter-flow reactive solids and purge gas, and a radiative exchange model applied to embedded-cavity structures at the size scale of the full receiver. In this work we apply the numerical analysis strategy to a single receiver configuration, but the framework can be generically applicable to alternative enclosed designs. In conclusion, we assess sensitivity of receiver performance to surface optical properties, heat-transfer coefficients, solids outlet temperature, and purge-gas feed rates, and discuss the significance of model assumptions and results for future receiver development.

14 SOLAR ENERGY↗

Geologic hydrogen as an emerging fuel: experimental insights, thermodynamics, kinetics, and reactive transport modeling

Geologic hydrogen (GeoH 2 ) is emerging as a viable clean energy source. It is largely produced through serpentinization, a geological process in which ultramafic rocks react with water under suitable temperature and pressure. Here, this review synthesizes the current understanding of H 2 generation by serpentinization, with an emphasis on reaction mechanisms, kinetics, and thermodynamics, as well as on modeling flow and transport of reacting fluids in geological formations. We describe the role of mineral assemblages, such as olivine and pyroxene, fluid-rock interactions, and catalytic surfaces, in influencing GeoH 2 yield and reaction rates. By integrating models of reaction kinetics, subsurface reactive flow and transport, and the serpentinization process, and by accounting for the thermodynamic state of the system, this review aims to guide future GeoH 2 research and to evaluate the potential of natural hydrogen as a sustainable clean energy source.

Moradi, Rasoul [Univ. of Southern California, Los ↗

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗

Computer study of emergency shutdowns of a 60-kilowatt reactor Brayton space power system

A digital computer study of emergency shutdowns of a 60-kWe reactor Brayton power system was conducted. Malfunctions considered were (1) loss of reactor coolant flow, (2) loss of Brayton system gas flow, (3)turbine overspeed, and (4) a reactivity insertion error. Loss of reactor coolant flow was the most serious malfunction for the reactor. Methods for moderating the reactor transients due to this malfunction are considered.

Tew, R. C.↗

Modelling Nuclear Thermal Propulsion Reactor Startup Transients

The National Aeronautics and Space Administration (NASA) has set the goal of a manned mission to Mars by the year 2030 [1] and charged the national academy of sci- ences "to identify primary technical and programmatic chal- lenges, merits, and risks for maturing space nuclear propulsion technologies of interest to a future human Mars exploration mission" [2]. One relevant technology, nuclear thermal propul- sion (NTP), has notable advantages over traditional chemical rockets; most important among them is the ability to produce larger specific impulse on the order of 900s. The reduction of mission time is crucial for a manned mission to Mars to reduce the risk for the crew. Due to its higher specific impulse, NTP satisfies this need and is pursued as one technology to get humans to Mars [3, 4]. The construction of an NTP sys- tem has to negotiate several challenges laid out in Ref. [2]; one of these challenges is the need to startup the NTP sys- tem from essentially cold conditions to full power within one minute. This paper focuses on studying the neutronics and thermal-hydraulics behavior of a simplified NTP model dur- ing prescribed reactivity insertions and mass flow rate (MFR) ramps. It is the goal of this paper to investigate startup, peak material temperatures, and average specific impulse for a low enriched Uranium (LEU), ceramic and metal material (CER- MET) NTP system when varying reactivity insertion and MFR ramps.

33 ADVANCED PROPULSION SYSTEMS↗

Guided plasma jets directed onto wet surfaces: angular dependence and control

The optimal use of atmospheric pressure plasma jets (APPJs) for treatment of surfaces—inorganic, organic and liquid—depends on being able to control the flow of plasma-generated reactive species onto the surface. The typical APPJ is a rare gas mixture (RGM) flowed through a tube to which voltage is applied, producing an RGM plasma plume that extends into the ambient air. The RGM plasma plume is guided by a surrounding shroud of air due to the higher electric field required for an ionization wave (IW) to propagate into the air. Additionally, the mixing of the ambient air with the RGM plasma plume then determines the production of reactive oxygen and nitrogen species (RONS). The APPJ is usually oriented perpendicular to the surface being treated. However, the angle of the APPJ with respect to the surface may be a method to control the production of reactive species to the surface due to the change in APPJ propagation properties and the resulting gas dynamics. In this paper, we discuss results from computational and experimental investigations addressing two points—propagation of IWs in APPJs with and without a guiding gas shroud as a function of angle of the APPJ with respect to the surface; and the use of this angle to control plasma activation of thin water layers. We found that APPJs propagating out of the plasma tube into a same-gas environment lack any of the directional properties of shroud-guided jets, and largely follow electric field lines as the angle of the plasma tube is changed. Guided APPJs propagate coaxially with the tube as the angle is changed, and turn perpendicularly towards the surface only a few mm above the surface. The angle of the APPJ produces different gas dynamic distributions, which enable some degree of control over the content of RONS transferred to thin water layers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CFD Simulations of a Regenerative Process for Carbon Dioxide Capture in Advanced Gasification Based Power Systems

This project describes the work carried out to prepare a highly reactive and mechanically strong MgO based sorbents and to develop a Population Balance Equations (PBE) approach to describe the evolution of the particle porosity distribution that is linked with Computational Fluid Dynamics (CFD) to perform simulations of the CO2 capture and sorbent regeneration. A large number of MgO-based regenerable sorbents were prepared using low cost and abundant dolomite as the base material. Among various preparation parameters investigated the potassium/magnesium (K/Mg) ratio was identified as the key variable affecting the reactivity and CO2 capacity of the sorbent. The optimum K/Mg ratio is about 0.15. The sorbent formulation HD52-P2 was identified as the “best” sorbent formulation and a large batch (one kg) of the sorbent was prepared for the detailed study. The results of parametric study indicate the optimum carbonation and regeneration temperatures are 360° and 500°C, respectively. The results also indicate that steam has a beneficial effect on the rate of carbonation and regeneration of the sorbent and that the reactivity and capacity of the sorbent decreases in the cycling process (sorbent deactivation). The results indicate that to achieve a high CO2 removal efficiency, the bed of sorbent should be operated at a temperature range of 370-410°C which also favors production of hydrogen through the WGS reaction. To describe the carbonation reaction kinetics of the MgO, the Variable Diffusivity shrinking core Model (VDM) was developed in this project, which was shown to accurately fit the experimental data. An important advantage of this model is that the changes in the sorbent conversion with time can be expressed in an explicit manner, which will significantly reduce the CFD computation time. A Computational Fluid Dynamic/Population Balance Equations (CFD/PBE) model was developed that accounts for the particle (sorbent) porosity distribution and a new version of the method of moments, called Finite size domain Complete set of trial functions Method Of Moments (FCMOM) was used to solve the population balance equations. The PBE model was implemented in a commercial CFD code, Ansys Fluent 13.0. The code was used to test the model in some simple cases and the results were verified against available analytical solution in the literature. Furthermore, the code was used to simulate CO2 capture in a packed-bed and the results were in excellent agreement with the experimental data obtained in the packed bed. The National Energy Laboratory (NETL) Carbon Capture Unit (C2U) design was used in simulate of the hydrodynamics of the cold flow gas/solid system (Clark et al.58). The results indicate that the pressure drop predicted by the model is in good agreement with the experimental data. Furthermore, the model was shown to be able to predict chugging behavior, which was observed during the experiment. The model was used as a base-case for simulations of reactive flow at elevated pressure and temperatures. The results indicate that by controlling the solid circulation rate, up to 70% CO2 removal can be achieved and that the solid hold up in the riser is one of the main factors controlling the extent of CO2 removal. The CFD/PBE simulation model indicates that by using a simulated syngas with a composition of 20% CO2, 20% H2O, 30% CO, and 30% H2, the composition (wet basis) in the reactor outlet corresponded to about 60% CO2 capture with and exit gas containing 65% H2. A preliminary base-case-design was developed for a regenerative MgO-based pre-combustion carbon capture process for a 500 MW IGCC power plant. To minimize the external energy requirement, an extensive heat integration network was developed in Aspen/HYSYS® to produce the steam required in the regenerator and heat integration. In this process, liquid CO2 produced at 50 atm can easily be pumped and sequestered or stored. The preliminary economic analyses indicate that the estimated cost of carbon v capture is in the range of $31-$44/ton, suggesting that a regenerative MgO-Based process can be a viable option for pre-combustion carbon dioxide capture in advanced gasification based power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

First principles simulation of reacting hypersonic flow over a blunt wedge

This article presents molecular-level analysis of a reactive, near-continuum, Mach 21 nitrogen flow over a blunt wedge using the direct molecular simulation (DMS) method. The flow conditions lead to internal energy excitation and dissociation in the flow field, resulting in thermal and chemical nonequilibrium in the flow. Thermal nonequilibrium in the vibrational mode is observed to extend to the molecular level, where the vibrational energy distributions at various points in the flow field are observed to be non-Boltzmann. Furthermore, this is the first reactive DMS calculation where the wall is assumed to be isothermal and full momentum accommodation of the particles is enforced, hence incorporating viscous wall effects. Since the DMS method uses a quantum mechanically generated interaction potential as its only modeling input, all thermochemical and transport properties of the flow field can directly be attributed to the ab initio potential energy surface. Using the DMS solution as a benchmark, this article assesses the performance of Navier–Stokes computational fluid dynamics solutions using lower fidelity two-temperature models. Two models are chosen as points of comparison: the well-known Park two-temperature model and the recently developed modified Marrone and Treanor model.

Mechanics↗

Non-Stationary Power System Forced Oscillation Analysis using Synchrosqueezing Transform

Non-stationary forced oscillations (FOs) have been observed in power system operations. However, most detection methods assume that the frequency of FOs is stationary. In this paper, we present a methodology for the analysis of nonstationary FOs. Firstly, Fourier synchrosqueezing transform (FSST) is used to provide a concentrated time-frequency representation of the signals that allows identification and retrieval of non-stationary signal components. To continue, the Dissipating Energy Flow (DEF) method is applied to the extracted components to locate the source of forced oscillations. The methodology is tested using simulated as well as real PMU data. In conclusion, the results show that the proposed FSST-based signal decomposition provides a systematic framework for the application of DEF Method to non-stationary FOs.

42 ENGINEERING↗

Integrated Experimental and Modeling Studies of Mineral Carbonation as a Mechanism for Permanent Carbon Sequestration in Mafic/Ultramafic Rocks

A program of laboratory experiments, modeling and fieldwork was carried out at Yale University, University of Maryland, and University of Hawai‘i, under a DOE Award (DE-FE0004375) to study mineral carbonation as a practical method of geologic carbon sequestration. Mineral carbonation, also called carbon mineralization, is the conversion of (fluid) carbon dioxide into (solid) carbonate minerals in rocks, by way of naturally occurring chemical reactions. Mafic and ultramafic rocks, such as volcanic basalt, are natural candidates for carbonation, because the magnesium and iron silicate minerals in these rocks react with brines of dissolved carbon dioxide to form carbonate minerals. By trapping carbon dioxide (CO 2 ) underground as a constituent of solid rock, carbonation of natural basalt formations would be a secure method of sequestering CO 2 captured at power plants in efforts to mitigate climate change. Geochemical laboratory experiments at Yale, carried out in a batch reactor at 200°C and 150 bar (15 MPa), studied carbonation of the olivine mineral forsterite (Mg 2 SiO 4 ) reacting with CO 2 brines in the form of sodium bicarbonate (NaHCO 3 ) solutions. The main carbonation product in these reactions is the carbonate mineral magnesite (MgCO 3 ). A series of 32 runs varied the reaction time, the reactive surface area of olivine grains and powders, the concentration of the reacting fluid, and the starting ratio of fluid to olivine mass. These experiments were the first to study the rate of olivine carbonation under passive conditions approaching equilibrium. The results show that, in a simple batch reaction, olivine carbonation is fastest during the first 24 hours and then slows significantly and even reverses. A natural measure of the extent of carbonation is a quantity called the carbonation fraction, which compares the amount of carbon removed from solution, during a run, to the maximum amount that could have been removed if the olivine initially present had fully dissolved and the cations released had subsequently precipitated in carbonate minerals. The carbonation fractions observed in batch experiments with olivine grains and powders varied significantly, from less than 0.01 (1%) to more than 0.5 (50%). Over time, the carbonation fractions reached an upper limit after about 24 to 72 hours of reaction, then stayed constant or decreased. The peak Final Scientific/Technical Report DE-FE0004275 | Mineral Carbonation | 4 coincided with the appearance of secondary magnesium-bearing silicate minerals, whose formation competes for magnesium ions in solution and can even promote conditions that dissolve magnesite. The highest carbonation fractions resulted from experiments with low ratios of concentrated solution to olivine, during which amorphous silica spheres or meshes formed, instead of secondary silicate minerals. The highest carbonation fractions appear to result from competing effects. Precipitation of silica layers on olivine reduces the reactive surface area and, thus, the rate of olivine dissolution (which ultimately limits the carbonation rate), but these same silica layers can also inhibit the formation of secondary silicate minerals that consume magnesite formed in earlier stages of carbonation. Simulation of these experiments with simple geochemical models using the software program EQ3/6 reproduces the general trends observed—especially the results for the carbonation fraction in short-run experiments. Although further experimentation and better models are needed, this study nevertheless provides a framework for understanding the optimal conditions for sequestering carbon dioxide by reacting CO 2 -bearing fluids with rocks containing olivine minerals. A series of experiments at the Rock Physics Laboratory at the University of Maryland studied the carbonation process during deformation of thermally cracked olivine-rich rock samples (dunite) saturated with CO 2 brines of varying compositions. A goal of these geomechanical experiments was to see if flow and deformation processes, which accompany natural carbonation reactions in underground settings, work to enhance or inhibit the reactions. The experiments involved hydrostatic compaction, followed by deformation at a constant rate of strain. Sample permeability was monitored during the reactions. Comparison of the samples’ volume changes to their axial strains (shortening) during deformation indicates that samples reacted with CO 2 -saturated brines accommodate more axial compaction, before the onset of dilation (a swelling that precedes rock failure), than samples reacted with distilled water. Analyses of the reacted samples with scanning electron microscope (SEM) images indicate, first, that dissolution of olivine occurring in the initial stages of carbonation can provide pathways to fluid flow that sustain the reaction, and, second, that carbonate minerals precipitated along existing fractures in the rocks may serve as asperities, or roughness on a crack’s surface that restricts its closure. Final Scientific/Technical Report DE-FE0004275 | Mineral Carbonation | 5 In a related study undertaken by one of the principal investigators as a spin-off of the main project, a simple model of (magnesite) crystal growth in the pore space of basalts undergoing carbonation was developed. The model suggests that, under a carefully controlled program of CO 2 injection, carbonate mineral growth can harden the rock formation against earthquakes that might otherwise be induced by the injection of large fluid volumes (Yarushina and Bercovici, 2013). The overall conclusion of the research project is that mineral carbonation of underground mafic and ultramafic rock formations is a viable candidate for long-term sequestration of man-made carbon dioxide. No results obtained during the project indicate that the method is inherently intractable in its implementation; moreover, enormous volumes of basalt near Earth’s surface are candidate locations for large-scale injection programs. The geochemical experiments do indicate, however, that there will be significant engineering challenges in maintaining high rates of carbonation, by delaying the onset of chemical conditions that promote formation of secondary silicate minerals and, therefore, slow down, or even reverse, the carbonation process. It remains an open question as to whether carbonation processes can be sustained for many years in an engineered system operating on a large scale—a scale capable of accommodating millions of tons of CO 2 annually. The development of realistic theoretical models that can systematically describe the combined effects of reactive flow, precipitation and geomechanical deformation is a major barrier to further understanding of the practical viability of mineral carbonation as large-scale method of carbon sequestration.

58 GEOSCIENCES↗

Control system design for a small pressurized water reactor - 141

Small Pressurized Water Reactors (SPWR) have unique designs from the commercial large Pressurized Water Reactors (PWRs). Integral design is adopted. There are no hot legs and cold legs between the reactor core and the steam generators like in the PWR. Strong interactions exist among different loops. The pressure vessel is large in volume and the coolant inventory is also in a large amount. The inertia of the coolant is large and it takes more time for the primary system to respond to disturbances. The steam generator adopts the once-through design and the water inventory is small. The steam generator is very sensitive to disturbances. These unique characteristics challenge the control system design of SPWR. Relap5 is used to model the SPWR. Steady-state simulation is performed to verify the correctness of the model. In the reactor power control system, both the reactor fission power and the coolant average temperature are regulated by the control rod reactivity. In the feedwater flow control system, the coordination between the reactor and the turbine is considered and coolant average temperature is adopted as one measurable disturbance for the control system. The coolant pressure is adjusted based on the heaters and spray in the pressurizer. The water level in the pressurizer is controlled by the charging flow and the letdown flow is compensated to calculate the charging flow needed. Transient simulations are carried out to evaluate the control system performance. When the reactor is perturbed, the reactor can be stabilized under the control system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗