Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Quality”

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,063 records · Page 59

THE EFFECT OF FUEL REACTIVITY AND EXHAUST GAS RECIRCULATION ON KNOCK PROPENSITY OF NATURAL GAS

The development of high efficiency, spark ignited natural gas engines is currently limited by engine knock at high compression ratio/elevated boost pressures and misfire at lean conditions/high exhaust gas recirculation (EGR) levels. The knock and misfire limits are further confounded by the wide variety in fuel reactivity observed in “pipeline quality” natural gas. In this study, a rapid compression machine was used to characterize the effects of EGR and variation in natural gas fuel reactivity on the homogeneous ignition delay, flame propagation rate, and end- gas autoignition propensity for stoichiometric natural gas/oxidizer/EGR blends. A reduced chemical kinetic mechanism was also developed to accurately model the homogeneous ignition delays measured in the Colorado State University rapid compression machine (CSU RCM). Pipeline quality natural gas with a range of chemical reactivity (68 < Methane Number < 95) was simulated using mixtures of CH4, C2H6, and C3H8. Exhaust gas recirculation gases were simulated with mixtures of Ar, CO2, CO, and NO at substitution rates of 0 to 30 mass percent. Ignition delay period under homogeneous autoignition conditions was measured at compressed pressures of 30.2 to 34.0 bar and compressed temperatures of 667 to 980 K. End-gas autoignition fraction and flame propagation rate were measured by initiating a laser spark in the center of the combustion chamber, after compression, at pressures of 30.7 to 32.7 bar and temperatures of 751 to 795 K. The results indicate that both fuel reactivity and the presence of reactive species (NO and CO) in the exhaust gas recirculation have a strong impact on end-gas autoignition fraction. A chemical kinetic mechanism was developed to predict homogeneous ignition delays for pipeline quality natural gas in a pressure and temperature range of 1-100 bar and 500-1000 K respectively. This mechanism accurately predicted measured homogeneous ignition delay in the RCM with a total average relative error of 11.0%.

Mohr, Jeffrey↗

Recent advances in superconducting-mixer simulations

Over the last few years, considerable progress have been made in the development of techniques for fabricating high-quality superconducting circuits, and this success, together with major advances in the theoretical understanding of quantum detection and mixing at millimeter and submillimeter wavelengths, has made the development of CAD techniques for superconducting nonlinear circuits an important new enterprise. For example, arrays of quasioptical mixers are now being manufactured, where the antennas, matching networks, filters and superconducting tunnel junctions are all fabricated by depositing niobium and a variety of oxides on a single quartz substrate. There are no adjustable tuning elements on these integrated circuits, and therefore, one must be able to predict their electrical behavior precisely. This requirement, together with a general interest in the generic behavior of devices such as direct detectors and harmonic mixers, has lead us to develop a range of CAD tools for simulating the large-signal, small-signal, and noise behavior of superconducting tunnel junction circuits.

Withington, S.↗

Achieving extreme light confinement in low-index dielectric resonators through quasi-bound states in the continuum

Obtaining large field enhancement in low-refractive-index dielectric materials is highly relevant to many photonic and quantum optics applications. However, confining light in these materials is challenging, owing to light leakage through coupling to continuum modes in the surrounding environment. Here, we investigate the possibility of achieving high quality factors in low-index dielectric resonators through the bound states in the continuum (BIC). Our simulations demonstrate that destructive interference between leaky modes can be achieved by tuning the geometrical parameters of the resonator arrays, leading to the emergence of quasi-BIC in resonators that have a small index contrast to the underlying substrates. The resultant large field enhancement gives rise to giant quality factors and Purcell effects. By introducing vertical mirror symmetry, the quasi-BIC can be tuned into an ideal BIC. In addition, the quasi-BIC can modify the emission patterns of the coupled emitters, rendering highly directional and focused far-field emission. These findings may provide a path for the practical implementation of photonic and quantum devices based on low-index dielectric materials.

42 ENGINEERING↗

MOOSE-based Tritium Migration Analysis Program, Version 8 (TMAP8) for advanced open-source tritium transport and fuel cycle modeling

Tritium management is critical for the safety, sustainability, and economics of fusion energy systems, and advanced and reliable modeling tools help accelerate the development of tritium technologies. This paper presents the Tritium Migration Analysis Program, Version 8 (TMAP8), an open-source, MOOSE-based application developed to provide state-of-the-art tritium transport and fuel cycle modeling capabilities. TMAP8 aims to expand the capabilities of previous versions (i.e., TMAP4 and TMAP7) by leveraging modern computational techniques, ensuring high software quality assurance standards (key to building trust), and enabling multispecies, multiscale, and multiphysics simulations for integrated tritium transport modeling in complex geometries. This paper outlines TMAP8’s scope and rigorous development practices, emphasizing its transparency, accessibility, modularity, and reliability. We present the current suite of verification and validation cases based on those from TMAP4, demonstrating TMAP8’s accuracy and reliability against analytical solutions and experimental data. Additionally, the paper showcases TMAP8’s integrated fuel cycle modeling capabilities, highlighting its applicability at various scales and levels. The TMAP8 code and documentation are openly available, promoting collaborative development and widespread adoption within the fusion community. Future work will soon expand TMAP8’s verification and validation suite to include those from TMAP7 and other recent experimental studies for validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Characterizing Detailed Grain Shape and Size Distribution Properties of Lunar Regolith

Introduction: As the nation prepares to return to the Moon, there is an increasing need for testing tools, instruments, and equipment in simulated environments on Earth to ensure successful operations during lunar missions. Regolith will affect all aspects of future lunar missions, from plume interactions during landing to space suit and tool design [1]. Because of this, it is important to understand the grain shape and size properties of lunar regolith and how those influence regolith behavior in order to prepare for these missions. This knowledge is also vital to create more accurate lunar regolith simulants for testing equipment in a lunar environment. While particle size analyses have been performed on most Apollo soils using simple sieving, shape has only been crudely addressed [2]. This work analyzes 4 lunar regolith samples to provide a better understanding of these size and shape parameters and will provide a more accurate baseline of data to create high-fidelity lunar regolith simulants. Methods: New technologies exist today that are capable of measuring size and shape simultaneously for hundreds of thousands of particles in a single measurement. We conducted a rigorous analysis of the particle size distribution (PSD) as well as the size-dependent 2D and 3D shape parameters of lunar regolith samples of different compositions and maturity levels. This analysis was done using a Microtrac SYNC which provides a unique combination of tri-laser diffraction and Dynamic Image Analysis (DIA). Sample Selection. Four regolith samples were selected for analysis based on maturity level and lunar terrain type: • 10084 – Mature high-Ti mare regolith • 15601 – Immature low-Ti mare regolith • 64501 – Mature highland regolith • 67461 – Immature highland regolith Sample Analysis. After receiving the samples, each sample was imaged with an optical microscope (Figure 1). To obtain 2D and 3D particle size and shape, 0.1 g of each sample was then added to the SYNC for analysis, which outputs more than 30 size and shape parameters for each individual grain as well as the complete size distribution from 0.01-2000 μm by blending laser diffraction and DIA together. For the 0.1 g sample masses requested, ~100,000 grains per sample were captured by DIA. Results: From the PSD analysis, it can be seen that the average particle size of 10084 is ~24.5 µm which is much smaller than the average particle size of the Apollo sample collection (~72 µm) [3], while the average particle size of samples 15601, 64501, and 67461 are larger than the Apollo sample average at ~106 µm, ~103 µm, and ~118.5 µm respectively. The size and shape measurements for the samples output ~30 parameters, four of which were focused on for this study: sphericity, aspect ratio, roundness, and concavity (Figure 3). However, after analyzing the plots it was found that only sphericity and aspect ratio showed differences between the samples. These two parameters are measured on a scale from 0 to 1, with 1 being a perfect sphere with equal dimensions. The results show that the sphericity values are slightly high-er for the mature samples 10084 and 64501 (~0.96) than they are for samples 15601 and 67461 (~0.95) (Figure 4). The aspect ratio values are slightly lower for samples 64501 (~0.7) and 67461 (~0.74) than they are for samples 10084 (~0.76) and 15601 (~0.8) (Figure 5). Discussion: Mature regoliths are those that have been exposed to micrometeorite and solar wind bombardment for long periods of time, breaking up particles and causing them to become more rounded [4]. Therefore, the smaller particle sizes and higher sphericity values for samples 10084 and 64501 are expected due to their higher maturity level compared to 15601 and 67461. However, the aspect ratio values are not dependent on maturity and are instead dependent on terrain type. The lower aspect ratio values for the high-land samples is potentially due to the higher plagioclase content, which occurs in elongated particles and does not break down as easily as the pyroxenes or olivines that are present in the mare. Conclusion and Future Work: The results of this work provide a baseline of high-quality data that will contribute to the creation of high-fidelity lunar simulants and will greatly benefit NASA’s efforts of establishing a human presence on the Moon. Future work includes performing unsupervised image classification on the ~105 particle images per sample in order to identify different classes of grains. These grain classes can then be linked to detailed shape properties, and the relative abundance of each class in the samples can be compared. Acknowledgments: We would like to thank the Extraterrestrial Materials Analysis Group (ExMAG) and the Astromaterials Allocation Review Board (AARB) for allocating lunar regolith samples 10084,9010, 15601,365, 64501,249, and 67461,171 for this work. References: [1] Taylor, L. A. et al. (2005) AIAA #2510. [2] Katagiri et al. (2015) ASCE. [3] Carrier III, W. D. (2005) Lunar Geotech Institute Tech Report. [4] McKay, D. S. et al. (1991) The Lunar Sourcebook, Chapter 7.

S R Deitrick↗

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part I - Analysis with a Small Sample Size

Computer simulations, especially of complex phenomena, can be expensive, requiring high-performance computing resources. Often, to understand a phenomenon, multiple simulations are run, each with a different set of simulation input parameters. These data are then used to create an interpolant, or surrogate, relating the simulation outputs to the corresponding inputs. When the inputs and outputs are scalars, a simple machine learning model can suffice. However, when the simulation outputs are vector valued, available at locations in two or three spatial dimensions, often with a temporal component, creating a surrogate is more challenging. In this report, we use a two-dimensional problem of a jet interacting with high explosives to understand how we can build high-quality surrogates. The characteristics of our data set are unique - the vector-valued outputs from each simulation are available at over two million spatial locations; each simulation is run for a relatively small number of time steps; the size of the computational domain varies with each simulation; and resource constraints limit the number of simulations we can run. We show how we analyze these extremely large data-sets, set the parameters for the algorithms used in the analysis, and use simple ways to improve the accuracy of the spatio-temporal surrogates without substantially increasing the number of simulations required.

97 MATHEMATICS AND COMPUTING↗

Surrogate-driven design optimization with uncertainty constraints in Monte Carlo simulations

In multi-objective design tasks, the computational cost increases rapidly when high-fidelity simulations are used to evaluate objective functions. Surrogate models help mitigate this cost by approximating the simulation output, simplifying the design process. However, under high uncertainty, surrogate models trained on noisy data can produce inaccurate predictions, as their performance depends heavily on the quality of training data. This study investigates the impact of data uncertainty on two multi-objective design problems modelled using Monte Carlo transport simulations: a neutron moderator and an ion-to-neutron converter. For each, a grid search was performed using five different tally uncertainty levels to generate training data for neural network surrogate models. These models were then optimized using NSGA-III. The recovered Pareto-fronts were analyzed across uncertainty levels: in the moderator problem, normalized hypervolume dropped from 0.886 at 1.0% uncertainty to 0.748 at 10% uncertainty, while in the converter problem it remained near 0.50 for all cases. Average simulation times were also compared to evaluate the trade-off between accuracy and computational cost. Results show that the influence of simulation uncertainty is strongly problem-dependent. In the neutron moderator case, higher uncertainties led to exaggerated objective sensitivities and distorted Pareto-fronts, reducing normalized hypervolume. In contrast, the ion-to-neutron converter task was less affected—low-fidelity simulations produced results similar to those from high-fidelity data. These findings suggest that a fixed-fidelity approach is not optimal. Surrogate models can recover the Pareto-front under noisy conditions, and multi-fidelity studies help identify suitable uncertainty levels for each problem to balance efficiency and accuracy.

07 ISOTOPE AND RADIATION SOURCES↗

High-fidelity wind farm simulation methodology with experimental validation

The complexity and associated uncertainties involved with atmospheric-turbine-wake interactions produce challenges for accurate wind farm predictions of generator power and other important quantities of interest (QoIs), even with state-of-the-art high-fidelity atmospheric and turbine models. A comprehensive computational study was undertaken with consideration of simulation methodology, parameter selection, and mesh refinement on atmospheric, turbine, and wake QoIs to identify capability gaps in the validation process. For neutral atmospheric boundary layer conditions, the massively parallel large eddy simulation (LES) code Nalu-Wind was used to produce high-fidelity computations for experimental validation using high-quality meteorological, turbine, and wake measurement data collected at the Department of Energy/Sandia National Laboratories Scaled Wind Farm Technology (SWiFT) facility located at Texas Tech University’s National Wind Institute. The wake analysis showed the simulated lidar model implemented in Nalu-Wind was successful at capturing wake profile trends observed in the experimental lidar data.

17 WIND ENERGY↗

Hybrid Grid Techniques for Propulsion Applications

During the past decade, computational simulation of fluid flow for propulsion activities has progressed significantly, and many notable successes have been reported in the literature. However, the generation of a high quality mesh for such problems has often been reported as a pacing item. Hence, much effort has been expended to speed this portion of the simulation process. Several approaches have evolved for grid generation. Two of the most common are structured multi-block, and unstructured based procedures. Structured grids tend to be computationally efficient, and have high aspect ratio cells necessary for efficently resolving viscous layers. Structured multi-block grids may or may not exhibit grid line continuity across the block interface. This relaxation of the continuity constraint at the interface is intended to ease the grid generation process, which is still time consuming. Flow solvers supporting non-contiguous interfaces require specialized interpolation procedures which may not ensure conservation at the interface. Unstructured or generalized indexing data structures offer greater flexibility, but require explicit connectivity information and are not easy to generate for three dimensional configurations. In addition, unstructured mesh based schemes tend to be less efficient and it is difficult to resolve viscous layers. Recently hybrid or generalized element solution and grid generation techniques have been developed with the objective of combining the attractive features of both structured and unstructured techniques. In the present work, recently developed procedures for hybrid grid generation and flow simulation are critically evaluated, and compared to existing structured and unstructured procedures in terms of accuracy and computational requirements.

Koomullil, Roy P.↗

Using an Advanced Distribution Management System Test Bed to Evaluate the Impact of Model Quality on Volt/VAR Optimization: Preprint

In this paper, we present a test bed for evaluating existing and future advanced distribution management system (ADMS) applications in a realistic laboratory setting, including other utility management systems and field equipment. We present an example of using it to evaluate the impact of the ADMS network model quality on a Volt/VAR optimization (VVO) application. The test bed integrates a commercial ADMS with a real-time simulation model of a utility distribution feeder. Representative power and controller hardware are integrated through hardware-in-the-loop (HIL) techniques. The performance of the ADMS VVO application is also evaluated for different levels of measurement density. Initial results indicate that a higher model quality achieves the highest possible energy savings while avoiding voltage violations, whereas a lower model quality results in increased energy savings but at the expense of more voltage violations.

ADMS↗

Simulation of HyMETS Flowfield Around Baby-SPRITE Entry Probe

Arc-jets are unique facilities used in research, development, and evaluation of high-temperature thermal protection systems for hypersonic vehicles and planetary entry systems. Thermochemical non-equilibrium computational fluid dynamics simulations have been carried out for the Hypersonic Materials Environmental Test System arc-jet facility to determine the size of a capsule model before arc-jet testing by better understanding of the physical phenomena. The results show the effect of the test article geometry and the importance of high-quality grids for accurate solutions. Accurate computational modeling of hypersonic flow fields inside arc-jets under simulated planetary entry conditions would help improve the design of thermal protection systems that may enable human exploration of the Moon, Mars, and beyond.

Ventura Diaz, Patricia↗

Between-Simulator Comparison of Motion Filter Order and Break Frequency Interaction Effects on Pilot Control Behavior

This paper investigates the interaction effects of motion filter order and break frequency on pilots' manual control behavior and control performance using two simulators. Eighteen pilots performed the experiment in the Vertical Motion Simulator (VMS) at NASA Ames Research Center and twenty pilots in the SIMONA Research Simulator at Delft University of Technology. The experiment used a full-factorial design with three motion filter orders (first-, second-, and third-order) and two filter break frequencies (0.5 and 2.0 rad/s), in addition to reference no-motion and full-motion conditions. Key task variables, such as the quality of the motion and visual cues and the characteristics of the sidestick, were matched across both simulators. Overall, the expected effects of filter order and break frequency variations were found, with both increasing order and increasing break frequency causing pilots to use less motion feedback in their control strategy, resulting in lower tracking performance. Furthermore, across the wide range of filter orders tested in the experiment, the existing Sinacori-Schroeder motion fidelity criterion was found to be a good predictor of the interaction effects of both filter settings on pilot control behavior. For the same motion condition, there was a consistent offset in the results between simulators, due to the more high-gain control strategy adopted by a number of the VMS pilots. Still, the observed relative trends in pilot control behavior and performance between motion conditions were equivalent in both simulators and thus accurately replicated.

manual control↗

Application of the Rossi-alpha method to simulations of HEU and organic scintillators

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook contains thousands of high-quality benchmark measurements that include detailed descriptions of experiment designs and exhaustive accounting of uncertainties. These benchmarks are used extensively by the criticality safety and nuclear data communities. A subcritical benchmark involving bare, highly enriched uranium (HEU) has yet to be accepted to the handbook. The Measurement of Uranium Subcritical and Critical (MUSiC) was proposed to address this gap. MUSiC will be performed at the National Criticality Experiments Research Center (NCERC) and will leverage bare HEU hemisphere shells that can be assembled in a variety of configurations. The measurement has the following key benefits: it will be the first set of fast, subcritical measurements of HEU to be submitted to ICSBEP; ten configurations will be measured spanning from nearly no multiplication to critical (e.g., it is a measurement of similar systems with different re activities); measurements of an all-HEU system will be of value for nuclear data validation, providing results that assist in determining biases in both measurements and simulations. A variety of detector systems will be used including a Neutron Multiplicity Array Detector (NoMAD) system (similar to the MC-15), four small volume 0.635 cm (#8960;) × 7.59 cm 3He detectors (ideal for measuring prompt neutron decay constants due to their fast recovery speed), and an array of eight 5.08 cm (#8960;) × 5.08 cm EJ-309 organic scintillator detectors. The focus of this work is the simulation of the subcritical MUSiC configurations including the EJ-309 detectors and the application of the Rossi-α method. The EJ-309 organic scintillator detector array is a new detection system that will be deployed at NCERC. Organic scintillators are ideal for measurements of fast, subcritical systems due to their short time response on the order of nanoseconds.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Analysis of Handling Qualities and Power Consumption for Urban Air Mobility (UAM) eVTOL Quadrotors with Degraded Heave Disturbance Rejection and Control Response

A piloted handling qualities study of urban air mobility (UAM) electric vertical take-off and landing (eVTOL) quadrotors was performed utilizing the Vertical Motion Simulator (VMS) facility at NASA Ames Research Center. Rotor speed and variable pitch-controlled variants of a six-passenger conceptual design vehicle were assessed with different levels of degradation to control response and disturbance rejection bandwidth (DRB) in the heave axis. In previous work, preliminary trends across several handling quality rating categories reflected the effects of these degradations. Additionally, the impact of using different test standards and turbulence on the ratings were discussed. This paper elaborates on those results, but also provides insight into unexpected trends observed during the study including: a disharmony in attitude response, subpar ratings for the baseline Level 1 performance vehicle, and excessive drift and yaw couplings observed in a lateral reposition maneuver. Moreover, shortcomings of the handling quality scales and comparisons of power consumption among the vehicles in the various test conditions are presented.

Handling Qualities↗

Optimization for Bioenergy Systems

The Sustainable Aviation Fuel (SAF) Grand Challenge (Langholtz, 2024 ) seeks to generate 35 billion gallons of SAF each year by 2050, with corn stover, an agricultural byproduct, playing a key role as a feedstock. This study develops an optimization framework to enhance the quality and quantity of corn stover while ensuring economic and environmental viability. Using the Decision Support System for Agrotechnology Transfer (DSSAT) crop model, we simulate the effects of cover crops on rotation yield, soil moisture balance, and nitrogen cycling across diverse climates and soils. The model outputs, including yield data and soil quality changes, inform a Mixed-Integer Linear Programming (MILP) optimization model. This model aims to maximize economic and environmental returns by incorporating production costs, direct and indirect income, and environmental incentives. The optimization model evaluates 280 agriculture management plans composed of various crop management strategies, including corn stover removal rates, cover crop adoption, and fertilization practices. It seeks to identify the optimal combination of crop and tillage decisions for each subfield, maximizing profits while enhancing soil carbon sequestration and reducing greenhouse gas emissions. Outputs include detailed subfield locations, optimal management plans, and profits per hectare and per acre, allowing for comparison with literature values on farm profits. This study provides a robust optimization framework supporting the SAF Grand Challenge by proposing economically viable and environmentally sustainable strategies for corn stover utilization. The findings highlight corn stover's potential as a sustainable feedstock for SAF production, offering practical solutions to enhance its quality and quantity while maintaining soil health. Idaho is used as a case study to demonstrate the framework's applicability and effectiveness in real-world scenarios. Langholtz, M. H., Davis, M., Hellwinckel, C., De La Torre Ugarte, D., Efroymson, R., Jacobson, R., Milbrandt, A., Coleman, A., Davis, R., Kline, K. L., Badgett, A., Curran, S., Schmidt, E., Theiss, T., Fried, J., English, B., Lambert, L., Cook, H., Field, J., ... Walker, L. (2024). 2023 Billion-Ton Report: An Assessment of U.S. Renewable Carbon Resources. https://doi.org/10.2172/2441098 DSSAT Foundation. (2025). Decision Support System for Agrotechnology Transfer (DSSAT). Retrieved from https://dssat.net/

09 - BIOMASS FUELS↗

Magnetic Designs of New First Target Beamline Magnets for the ORNL SNS Upgrade

The Spallation Neutron Source (SNS) at ORNL currently is being upgraded from 1.0 GeV to 1.3 GeV. Several water-cooled magnets should be upgraded to transport 30% of higher beam energy. New chicane, injection/extraction septum, and Lambertson magnets were designed. Designing the magnets was a challenging task because the new magnets required good combined integrated field quality and needed to occupy the old magnets space but with about 20% greater integrated magnetic field. Additional strong requirements applied to the magnets fringe field do not disturb the circulating beam. The special field profiles had to be provided in foil areas between magnets. The analysis described here was based on OPERA3D simulations. A special technique was used for analyzing the integrated field harmonics. Initially, the particle track was simulated, and integrated field components were calculated along this track for the reference radius, which were used for the harmonics analysis. In addition, 3D field maps were provided for beam optics simulations. The final beamline analysis confirmed good beam transmission and low losses.

3d field maps↗

Reference H Cycle 3 Stability, Control, and Flying Qualities Batch Assessments

This work is an update of the assessment completed in February of 1996, when a preliminary assessment report was issued for the Cycle 2B simulation model. The primary purpose of the final assessment was to re-evaluate each assessment against the flight control system (FCS) requirements document using the updated model. Only a limited number of final assessments were completed due to the close proximity of the release of the Langley model and the assessment deliverable date. The assessment used the nonlinear Cycle 3 simulation model because it combines nonlinear aeroelastic (quasi-static) aerodynamic with hinge moment and rate limited control surface deflections. Both Configuration Aerodynamics (Task 32) and Flight Controls (Task 36) were funded in 1996 to conduct the final stability and control assessments of the unaugmented Reference H configuration in FY96. Because the two tasks had similar output requirements, the work was divided such that Flight Controls would be responsible for the implementation and checkout of the simulation model and Configuration Aerodynamics for writing Madab "script' files, conducting the batch assessments and writing the assessment report. Additionally, Flight Controls was to investigate control surface allocations schemes different from the baseline Reference H in an effort to fulfill flying qualities criteria.

Henderson, Dennis K.↗

Mobilization of trace metals from caprock and formation rocks at the Illinois Basin – Decatur Project demonstration site under geological carbon dioxide sequestration conditions

One concern for geologic CO 2 sequestration is the potential leakage of CO 2 or CO 2 -saturated brine containing trace metals into overlying aquifers, which poses the risk of adversely affecting underground sources of drinking water. In this work, rock and brine samples collected in the Illinois Basin – Decatur Project (IBDP), a large-scale geological CO 2 sequestration demonstration project in Decatur, Illinois, USA, as well as synthetic brine samples were used to understand the mobilization of trace metals as a result of CO 2 –rock–brine interactions under IBDP-specific conditions (50 °C and 20.6 MPa). Rock sample characterization indicated that, at the IBDP site, trace metal concentrations were greatest in cap rock samples from the Eau Claire Formation, compared with those from formations above and below this Formation. The natural brine samples collected from the Mt. Simon Sandstone reservoir were highly reducing and saline, with trace metal concentrations up to 680 times greater than the U.S. Environmental Protection Agency-prescribed drinking water standards. Batch leaching experiments indicated that both trace metal mobilization from rock and immobilization from the brine occurred when high-pressure CO 2 was introduced into the rock-brine system. The amounts of metals mobilized from the rock generally accounted for <5% of the total metals in the rock, but for some metals, including Ni, Pb, and Tl, up to 63% of the metals in the rock were mobilized to the brine. Leaching of trace metals into synthetic brines was different from that into the natural brine. The results of this study provided current information on the trace element sources and will help in risk simulations and experimental design to further evaluate potential impact of leakage on groundwater quality, which is important for future GCS projects to consider for monitoring and containment assurance purposes.

58 GEOSCIENCES↗