SEARCH · Engineering Papers

Results for “SECONDARY 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

305 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗

DOE OSTI2023

TEAMER: Water Tunnel Data from Testing the Pterofin Skimmer Concept

Pterofin's Skimmer concept relies on a flapping and pitching hydrofoil to extract hydrokinetic energy from water flows. The concept aims to utilize unsteady fluid dynamics phenomena (added mass, shed vorticity, and unsteady boundary layer development) to achieve higher lift coefficients, enabling increased power density of the hydrokinetic device and a fundamental shift in the rpm/torque scaling of the power take off compared with turbines. The Applied Research Laboratory at Penn State, in collaboration with Pterofin, designed and built a proof-of-concept flapping/pitching mechanism which was subsequently tested in ARL's 12-inch water tunnel facility. The mechanical power supplied to or extracted from the mechanism was measured for a range of hydrofoils provided by Pterofin over operating conditions including reduced frequency, Reynolds number, and the ratio between pitching and flapping amplitudes. The power lost to friction in the mechanism was removed from the net power measurement by means of a bare hub tare, with the resultant hydrodynamic power being used to calculate a mechanism-independent and non-dimensional power coefficient. The product of this effort is a dataset describing the power coefficient of a hydrofoil having simultaneous pitching and flapping motions, both of which are approximately sinusoidal. Power coefficients were collected for a range of primary design variables including: - Reduced frequency: 0.01 to 0.95 - Pitching/flapping peak angle ratio: 1.5 to 3.0 - Chord-based Reynolds number: 60,000 to 560,000 Secondary design variables relating to the hydrofoil geometry were explored including: - Aspect ratio - Planform shape - Section thickness distribution - Hydrofoil position relative to the pitching axis - Hydrofoil sweep angle relative to the pitching axis Measured data are provided in mean and time series formats. MATLAB scripts are provided which can be used to generate figures of time-averaged and phase-averaged hydrodynamic power coefficients calculated from the measured data. A complete description of the experiment and data reduction can be found in the Post Access Report for the Pterofin Skimmer test effort which will be available on the TEAMER website. This work was supported by the Pacific Energy Ocean Trust via a TEAMER award.

16 TIDAL AND WAVE POWER↗

DOE OSTI2020

First on-line detection of radioactive fission isotopes produced by laser-accelerated protons

Abstract The on-going developments in laser acceleration of protons and light ions, as well as the production of strong bursts of neutrons and multi- $$\hbox {MeV}$$ MeV photons by secondary processes now provide a basis for novel high-flux nuclear physics experiments. While the maximum energy of protons resulting from Target Normal Sheath Acceleration is presently still limited to around $$100 \, \hbox {MeV}$$ 100 MeV , the generated proton peak flux within the short laser-accelerated bunches can already today exceed the values achievable at the most advanced conventional accelerators by orders of magnitude. This paper consists of two parts covering the scientific motivation and relevance of such experiments and a first proof-of-principle demonstration. In the presented experiment pulses of $$200 \, \hbox {J}$$ 200 J at $$\approx \, 500 \, \hbox {fs}$$ ≈ 500 fs duration from the PHELIX laser produced more than $$10^{12}$$ 10 12 protons with energies above $$15 \, \hbox {MeV}$$ 15 MeV in a bunch of sub-nanosecond duration. They were used to induce fission in foil targets made of natural uranium. To make use of the nonpareil flux, these targets have to be very close to the laser acceleration source, since the particle density within the bunch is strongly affected by Coulomb explosion and the velocity differences between ions of different energy. The main challenge for nuclear detection with high-purity germanium detectors is given by the strong electromagnetic pulse caused by the laser-matter interaction close to the laser acceleration source. This was mitigated by utilizing fast transport of the fission products by a gas flow to a carbon filter, where the $$\upgamma$$ γ -rays were registered. The identified nuclides include those that have half-lives down to $$39 \, \hbox {s}$$ 39 s . These results demonstrate the capability to produce, extract, and detect short-lived reaction products under the demanding experimental condition imposed by the high-power laser interaction. The approach promotes research towards relevant nuclear astrophysical studies at conditions currently only accessible at nuclear high energy density laser facilities.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Pilot-Scale Testing of an Integrated Circuit for the Extraction of Rare Earth Minerals and Elements from Coal and Coal Byproducts Using Advanced Separation Technologies

The primary objective of this project was to develop and demonstrate an integrated pilot-scale circuitry for recovering high-value rare earth elements (REEs) from coal and coal byproducts. The target performance was to produce a mixed REE product with content of at least two percent by weight on a dry mass basis in a cost-effective and environmentally benign manner. During the first nine months of the project period (Phase 2 Budget Period 2), pilot plant construction was completed including all field site startup activities such as permitting, engineering design, procurement/bidding, unit fabrication, site construction, equipment installation, module assembly, safety training, and circuit shakedown. During the remaining 21 months of project period (Phase 2 Budget Period 3), detailed field-testing activities were performed including feedstock sample collection and preparation, exploratory testing, circuit modification, detailed parametric study, and performance optimization. A detailed techno-economic analysis was performed based on the pilot plant testing findings which provided various scenarios for REE production. The project successfully accomplished the proposed target performance by producing mixed rare earth oxide (REO) with greater than 90% purity by weight in a continuous pilot scale operation from two distinctly different coarse refuse materials (i.e., West Kentucky No. 13 and Fire Clay coal seams), and at least three secondary sources (i.e., heap leach process and naturally formed acid mine drainage system). Project partners included the University of Kentucky, Virginia Tech, West Virginia University, Alliance Coal, Blackhawk Mining, Mineral Refining Company, and Mineral Separation Technologies. The pilot scale test facility was constructed at a former mining complex owned by Alliance Natural Resource Partners (Alliance Coal). The site was rehabilitated to accommodate the equipment installation, construction and fabrication, electrical power requirement, water line management and containment. The process units constructed and installed included X-ray sorting unit, crushing and grinding unit, physical separation unit, acid leaching unit, solvent extraction unit, and wastewater management unit. A rare earth mineral concentration unit was constructed as a standalone unit for flexible operation. A detailed environmental assessment and control plan was carried out to identify and quantify any potential impacts of the pilot-scale processing circuitry on the human and eco-system health and well-being. Corresponding mitigation strategies and control measures were provided. A conceptual flowsheet was developed to effectively remove thorium and uranium from high purity rare earth oxide mix or any potential radionuclide enriched stream. The two distinct feedstock materials were secured from the Blackhawk Mining Complex in eastern Kentucky where the Fire Clay (Hazard No. 4) seam is processed. The West Kentucky No. 13 (Baker) coarse refuse material was collected from an active process stream at an Alliance coal preparation plant located in western Kentucky. Characterization analysis indicated that both of feed materials generated from the two sources contained >300 ppm of TREEs on a dry whole mass basis which met the requirements for a qualified feed stock. The two feedstocks were further upgraded using a dual x-ray sorter to prepare the feed material for hydrometallurgical circuit. Thermal treatment on feed material prior to leaching was found to: 1) improve the leaching recovery of REEs, 2) increase the leaching kinetics, and 3) allow the leaching reaction to occur at lower acidity. Roasting at 600°C was selected as the pre-treatment condition for both West Kentucky No. 13 and Fire Clay coarse refuse material. Over 40% of leaching recovery was achieved by roasting West Kentucky No. 13 material having a top particle size of 3 mm in the pilot scale operation using 1.2M sulfuric acid leaching at 75OC. Initial pilot scale testing involved continuous operation of the pilot plant for 94 hours. The leaching unit was operated at solid-to-liquid ratio of 1 to 10 (w/v) using 0.5M sulfuric acid solution at a temperature of 75°C. The continuous solvent extraction circuit utilized rougher and cleaner units with DEHPA and TBP as the extractants. An innovative stripping circuit was developed to accumulate the REE concentration in the stripping solution to a level above 600 ppm. A bleed stream from the recycled strip solution was treated using oxalic acid precipitation which produced a high grade rare earth oxalate. The oxalate product was roasted to remove the oxalate which produced a rare earth oxide product having a purity greater 90%. Due to high concentrations of contaminant ions in the pregnant leach solution (PLS), a modified flowsheet was developed that involved pre-concentration of the REEs using multiple stages of precipitation and redissolution. The advantage of this process was improved removal of contamination before the downstream purification process and a significant cost reduction relative to the circuit that utilized the solvent extraction process. The modified circuitry included processes involving leaching, multistage precipitation, redissolution, and oxalate precipitation followed by roasting of the oxalate product. The circuit produced a mixed REO that was 92.96% pure from the initial test. A detailed parametric test plan was carried out which involved varying key parameters including solids feed rate, acid flowrate, acid concentration, multistage precipitation pH, redissolution pH, oxalate precipitation dosage and pH. The response variables included REE recovery, contaminant recovery, REO product grade and overall chemical consumption. Test results indicated that the acid-to-solid ratio is the key parameter to leaching efficiency as performance deteriorated with an increase in solids concentration. The optimal pH determined for REE precipitation and redissolution was 6.5 and 2.5, respectively. Additional tests were conducted to further improve the flowsheet. Recirculating a portion of the PLS to the feed of the leach tanks improved the leaching performance by lowering the pH of the leaching system and reducing the contamination recovery by shortening the residence time. Moreover, the removal of Al prior to REE precipitation significantly reduced the oxalic acid consumption in the oxalate precipitation circuit. The modified circuit produced over 90% grade REO by weight from both West Kentucky No. 13 and Fire Clay coarse refuse material in pilot scale continuous test programs. A case analysis model was developed to project the REE and major contaminants concentration in each PLS stream based on the leaching condition, pH cut point, and oxalic acid dosage. A correlation was established using empirical and semi-empirical models. Using the models, chemical consumption required for each stage was predicted based on the projected performance of the hydrometallurgy circuit. After identifying the optimum conditions, validation tests were carried out for the treatment of both West Kentucky No. 13 and Fire Clay coarse refuse materials in the pilot plant. The actual circuit performance and chemical consumptions were very close to the model predictions. Other than the two coarse refuse sources, several secondary feedstocks were also tested in the pilot plant facility. A “heap leach” system was constructed using the coarse refuse material generated from cleaning the West Kentucky No. 13 seam coal. Using the two stage SX rougher and cleaner circuit, a concentrate with a grade >90% REO was produced while recovering >97% of the REEs from the heap leach PLS. Naturally generated acid mine drainage (AMD) from West Kentucky No.13 mine was processed using the multistage precipitation circuit in the pilot plant in a test conducted for a period of 32 hours. The final grade of the mix RE oxide produced from the AMD was 90.84% with an overall circuit recovery of 64%. The primary source of REE was the selective precipitation steps involving iron and aluminum rejection. The hydrophobic-hydrophilic separation (HHS) process was proven to effectively recover coal from fine waste materials. For REM recovery, the HHS process was able to produce concentrates at grades of approximately 1.8% REE on an ash basis; however, recovery values were typically low, <10%, under the optimal conditions determined in the laboratory-scale testing. Staged testing of the pilot-scale HHS process for coal recovery and semi-continuous laboratory testing for REM testing showed that a total concentration ratio of more than 15x was observed for the REM recovery process. A circuit simulation package was developed for REE extraction and purification using a spreadsheet-based platform (Microsoft Excel). The REESim circuit simulation package is configured to track the mass and volume flows of components passing through a series of unit operations specified and configured by the user. The mass rates can then be utilized by the user to determine important performance indicators such as product mass yields, concentrate purity levels, element-by-element recoveries, and so forth. The techno-economic analysis showed that the roasting and leaching operations were the most expensive capital items, each contributing approximately 30% to the total capital cost. One notable contributor to the high production costs was the low REE recovery observed in the pilot scale trials. The product basket price was shown to have a strong influence on the economic viability of the scenarios, with the scandium price being the most significant influencer. Operating cost was shown to be extremely sensitive to REE recovery, REE feed grade, and leaching acid consumption. An analysis of ten different scenarios for a 500 t/h commercial operation revealed that three were economically favorable, producing internal rates of return varying from 27.7% to 33.1% and payback periods of 4 to 5 years. The project successfully developed and demonstrated a process to recover REEs from coal and coal byproducts in a pilot-plant operation which consistently produced over 90% grade REO mix from varies types of feedstocks. Commercialization analysis showed that the technology readiness level successfully achieved TRL 6 at the end of the project and demonstrated the need and the potential for scaling the process to further advance the technologies toward the goal of providing a domestic supply of REEs at a commercial scale.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

A minor respiratory process with major global implications: is atmospheric methane oxidation in tree stems driven by stem respiration rather than microbial methanotrophy?

Tree stem surfaces are widely recognized as sites of carbon dioxide (CO₂) efflux and oxygen (O₂) influx, reflecting the dynamics of aerobic respiration of photosynthate substrates, such as sugars, delivered via the phloem. Stems are also largely considered passive conduits for methane (CH₄) produced in anoxic soils via microbial methanogenesis, where CH₄ is thought to be transported upward through the transpiration stream and/or diffusion and emitted through stem surfaces and the canopy. However, recent observations from dynamic stem chambers suggest that stems may also act as active sinks for atmospheric CH₄. Despite these findings, the extent and drivers of stem CH₄ consumption remain poorly characterized across biomes, species, and environmental gradients, and its quantitative relationship to stem respiration has not been established. Moreover, previous studies captured only snapshot fluxes, leaving diurnal patterns of CH₄ exchange uncharacterized. Here, we address these limitations by combining real-time measurements of stem CH₄ and O₂ uptake under ambient conditions in a California cherry tree, using a dynamic stem gas exchange system with three chambers receiving a continuous flow of ambient air and automated chamber and reference air sampling every 10 min. Our results confirm that stems of upland trees can actively consume both atmospheric CH₄ and O₂, but with decreasing temperature sensitivity as daily temperatures increase. Early mornings were marked by rapid influxes of both gases, followed by declining uptake as temperatures rose further. Methane uptake was tightly coupled with O₂ influx and represented a minor (0.012% ± 0.002%) fraction of stem respiratory activity, as determined by concurrent O₂ uptake. These findings suggest that while atmospheric CH₄ oxidation is a minor respiratory process in stems, it is strongly linked with stem physiological activity. This challenges the current assumption that terrestrial CH₄ uptake is driven solely by microbial methanotrophy and raises the possibility that living stem tissues may contribute to CH₄ oxidation through an as-yet-unidentified plant-based mechanism.

Atmospheric greenhouse gases↗

DOE OSTI2021

MBE growth and donor doping of coherent ultrawide bandgap AlGaN alloy layers on single-crystal AlN substrates

Single-crystal Aluminum Nitride (AlN) crystals enable the epitaxial growth of ultrawide bandgap Al(Ga)N alloys with drastically lower extended defect densities. Here, we report the plasma-MBE growth conditions for high Al-composition AlGaN alloys on single-crystal AlN substrates. An AlGaN growth guideline map is developed, leading to pseudomorphic Al x Ga 1–x N epitaxial layers with x ~0.6–1.0 Al contents at a growth rate of ~0.3 μm/h. These epitaxial layers exhibit atomic steps, indicating step flow epitaxial growth, and room-temperature band edge emission from ~4.5 to 5.9 eV. Growth conditions are identified in which the background impurity concentrations of O, C, Si, and H in the MBE layers are found to be very near or below detection limits. An interesting Si segregation and gettering behavior is observed at the epitaxial AlGaN/AlN heterojunction with significant implications for the formation and transport of 2D electron or hole gases. Well-controlled intentional Si doping ranging from ~2 × 10 17 to 3 × 10 19 atoms/cm 3 is obtained, with sharp dopant density transition profiles. In Si-doped Al 0.6 Ga 0.4 N epilayers, a room-temperature free electron concentration of ~3 × 10 19 /cm 3 , an electron mobility of ~27 cm 2 /V s, and an n-type resistivity of ~7.5 m Ω cm are obtained. Here, the implications of these findings on electronic and photonic devices on single-crystal AlN substrates are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Feasibility of Lead Fast Reactor Heat Exchanger Tube Online Monitoring

Within virtually all reactor systems with heat exchangers, the thinnest barriers between primary and secondary fluids are the heat exchanger tubes. Consistent with light water reactor experience, the advanced reactor heat exchangers are the most susceptible regions for corrosion and leakage (NRC, 2019). Moreover, traditional inspections and repairs using robotics will be expensive or impractical. On-line structural health corrosion monitoring (OLSHM) is therefore desirable to detect tube degradation prior to leaks that may allow mixing of heat exchanger fluids. Rhetorically many molten salt reactor (MSR) designers plan for on-line monitoring of critical reactor components including heat exchangers but there are no proven OLSHM sensors or systems for MSR heat exchangers. Guided wave ultrasound from the tube-end or ID using conventional low-temperature piezoelectric sensors has been shown to be sensitive to pits and cracks less than 50% through-wall in long pipes and tubes (Narayanan M.M., 2019) . Mode and frequency selection are important for practical sensor design decisions and ultimate performance. The L-0 wave-mode is easily excited by L-mode transducers but they have surface normal wave motion that will be sensitive to and more attenuative because of fluid interaction with the tube surface (Rose J.L., 1994). Shear Horizontal (SH) wave modes have material motion parallel to the tube surface therefore are virtually unaffected by fluid interaction. For SH-0 to be nondispersive in steel, the frequency-thickness product (MHz-mm) must be less than approximately 1.5. SH mode transducers must be coupled by a material that can support shear stress; a thin layer of honey works well for temporary sensors and an adhesive or a brazed metallurgical bond for permanent sensors. If sensors are incorporated into the fabrication process, the same kinds of SH-0 waves can be generated from the stagnant flow area of the tube OD without compromising the heat-exchange tube function. Piezoelectric materials are also available that can withstand harsh temperatures of molten salt heat exchangers (Stevenson T., 2015) (Machura et al., 2008). Though they may not have as high-performance characteristics, their performance only needs to be good enough to detect flaws of interest. This project attempts to demonstrate feasibility to develop a sensor that can be sensitive to flaws of interest and perform long-term monitoring of a molten salt reactor heat exchanger. Benefits to the advanced reactor system include: • Minimizing the expense of a reactor shut down for periodic inspection—only shutdown for cause. • Reduce risk of unexpected tube rupture by detecting flaws before they grow to through-wall leaks. • Cost reduction by minimizing need to design for manual/robotic inspectability. Transducer design, characterization of the beam profile, wedges vs. direct application of piezoelectric sensors, and ultrasonic pulse-echo responses from a representative heat exchange tube with fabricated flaw types are discussed.

On-Line Tube Monitor, Guided Wave Ultrasound, SH-0↗

DOE OSTI2020

Progress Toward Molten Salt Reactor Heat Exchanger On-Line Monitoring

Within virtually all reactor systems with heat exchangers, the thinnest barriers between primary and secondary fluids are the heat exchanger tubes. Consistent with light water reactor experience, the advanced reactor heat exchangers are the most susceptible regions for corrosion and leakage (NRC, 2019). Moreover, traditional inspections and repairs using robotics will be expensive or impractical. On-line structural health corrosion monitoring (OLSHM) is therefore desirable to detect tube degradation prior to leaks that may allow mixing of heat exchanger fluids. Rhetorically many molten salt reactor (MSR) designers plan for on-line monitoring of critical reactor components including heat exchangers but there are no proven OLSHM sensors or systems for MSR heat exchangers. Guided wave ultrasound from the tube-end or ID using conventional low-temperature piezoelectric sensors has been shown to be sensitive to pits and cracks less than 50% through-wall in long pipes and tubes (Narayanan M.M., 2019) . Mode and frequency selection are important for practical sensor design decisions and ultimate performance. The L-0 wave-mode is easily excited by L-mode transducers but they have surface normal wave motion that will be sensitive to and more attenuative because of fluid interaction with the tube surface (Rose J.L., 1994). Shear Horizontal (SH) wave modes have material motion parallel to the tube surface therefore are virtually unaffected by fluid interaction. For SH-0 to be nondispersive in steel, the frequency-thickness product (MHz-mm) must be less than approximately 1.5. SH mode transducers must be coupled by a material that can support shear stress; a thin layer of honey works well for temporary sensors and an adhesive or a brazed metallurgical bond for permanent sensors. If sensors are incorporated into the fabrication process, the same kinds of SH-0 waves can be generated from the stagnant flow area of the tube OD without compromising the heat-exchange tube function. Piezoelectric materials are also available that can withstand harsh temperatures of molten salt heat exchangers (Stevenson T., 2015) (Machura et al., 2008). Though they may not have as high-performance characteristics, their performance only needs to be good enough to detect flaws of interest. This project attempts to demonstrate feasibility to develop a sensor that can be sensitive to flaws of interest and perform long-term monitoring of a molten salt reactor heat exchanger. Benefits to the advanced reactor system include: • Minimizing the expense of a reactor shut down for periodic inspection—only shutdown for cause. • Reduce risk of unexpected tube rupture by detecting flaws before they grow to through-wall leaks. • Cost reduction by minimizing need to design for manual/robotic inspectability. Transducer design, characterization of the beam profile, and ultrasonic pulse-echo responses from a representative heat exchange tube with fabricated flaw types are discussed.

On-Line Tube Monitor, Guided Wave Ultrasound, SH-0↗

DOE OSTI2023

Optimal gas cooling with minimal aero-opticaldistortion for next generation high-average-power lasers (Progress report, Year 2)

Next-generation lasers that combine high-peak-power ultrashort pulses with high repetition rates have the potential to efficiently deliver the high-average-power outputs required for inertial fusion energy generation, laser-based particle acceleration, and to generate compact secondary x- and 𝛾-ray radiation sources, among others. The combination of long gain length with high surface-to-volume ratio to remove residual heat makes gas-cooled multi-slab amplifier heads an inherently aperture- and average power scalable technology that enables this new class of solid-state lasers, and their thermal management becomes the most pressing technical challenge. As the index-of refraction depends on the gas density, heat removal from the slab surfaces unavoidably entails scattering of propagated light from regions of high temperature fluctuations. The requirement to maximize convective heat extraction from the gain medium under the constraint of minimum optical aberration of the beam renders amplifier slab cooling a technological challenge at the intersection of aero-optics, fluid-mechanics and heat transfer. The goal of this research is to leverage computational fluid dynamics and aero optics simulations of different levels of fidelity to devise a numerical modeling and optimization framework for gas-cooled multi-slab amplifiers. The project will be conducted in close collaboration with the Advanced Photon Technologies Group at Lawrence Livermore National Laboratory, where a gas-cooled test stand for laser materials and thermal management structures will be developed and installed. The final goal of the collaborative effort is to manufacture and benchmark a physical prototype of an numerically optimized amplifier head design. Measurement data on solid-state and thermo-optical distortions will be used for the validation of the numerical models. The key technical challenge in this multi-objective design problem is that of balancing optimal heat removal with minimal optical distortion. To facilitate computationally efficient design optimization and, at the same time, deepen our basic understanding of the interaction of light with turbulence, a complementary study that combines high and low-fidelity flow simulations and aero-optical models is proposed. Along the lines of these two objectives and solution strategies, the project is organized into two thrusts.

42 ENGINEERING↗

DOE OSTI2021

Integrated parameter and process learning for hydrologic and biogeochemical modules in Earth System Models

Focus area: Primary focal area #2; secondary focal area #3: Learning about parameters and processes of land surface hydrologic and biogeochemical models in Earth System models by integrating machine learning, physics, and big data. Science challenges: How do we maximally leverage big-data observations to improve hydrobiogeochemical process description and parameterization so that such modules more realistically capture hydrologic and vegetation responses and feedbacks under the future climate? For example, how can we leverage physics, limited observations of vegetation and streamflow to better estimate evapotranspiration, and, relatedly, net primary productivity, especially for drought areas? Vegetation plays a critical role in regional and global water cycles; however, existing vegetation models have failed to predict vegetation response to droughts (McDowell & Xu, 2017) , arctic greening (Keenan & Riley, 2018) , and critical transitions between forest and savanna (Hirota et al., 2011) . These studies suggest that when we build process-based models (PBM) parameterized from regional and global plant traits, we tend to poorly describe plant adaptation and local-scale competition processes. The models and their associated parameters assigned for different regions in the world are not capturing essential heterogeneity in vegetation responses at finer spatial scales. Many parameters of the land surface models control hydrology and vegetation dynamics at the same time. The heterogeneity in vegetation response is a function of (i) plant type, (ii) plant size, (iii) competition and succession, (iv) environmental controls, and (v) local variations due to the unique ecological community that are very difficult to describe (e.g., the size of gaps resulting from fire that facilitated the coexistence of pioneering species). In the demographic models, only factors (i) and (iv) were captured, and plant types were generally described only by leaf phenology and climate zones. With current demographic models, we generally consider more traits to define plant types (i) and calibrate these traits to consider factors (ii), (iii) and (iv); however, it is substantially challenging to scale to regional and global simulations due to trait variations across space (Ali et al., 2016). Moreover, it has been noted that hillslope processes, including ridge-to-valley flow and sunny vs. shady slopes are primary organizers of water, energy, and vegetation (Clark et al., 2015; Fan et al., 2019) . Although gradual improvements in the hydrologic model component in earth system models may reduce this error (at a remarkably slow pace), the long-term, gradual impact of hydrology on plant traits are not well captured. Recent work showed that the hydrologic controls exerted by groundwater and lateral flow are primary regulators of rooting depth (Fan et al., 2017) . Such hydrologic controls have seldom been reflected in vegetation model parameterizations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Transforming Aeration Energy in Water Resource Recovery Facilities (WRRFs) through Suboxic Nitrogen Removal (Final Report)

The objective of this project was to advance two key technological components—aeration control strategies and process design methodologies—to support the development and broader adoption of suboxic biological nitrogen removal (SBNR). The project focused on achieving the following three goals: • Enhance Model Predictive Control (MPC) Technology: Advance the DO/Nmaster MPC platform from its initial 2018 pilot deployment at the Chico Water Resource Recovery Facility in California to full-scale integration. This included partnering with a blower technology commercialization partner and incorporating machine learning (ML) capabilities to enable nationwide deployment. • Bridge Knowledge Gaps in SBNR Process Design: Address fundamental gaps in SBNR process understanding through controlled pilot-scale testing at a dedicated pilot facility. These efforts supported the development of robust kinetic models to inform reliable SBNR control, operational strategies, and design frameworks. • Demonstrate Full-Scale Implementation of Low DO/SBNR with ML: Transition low dissolved oxygen (DO)/SBNR coupled with ML from pilot-scale trials to full-scale demonstration in flow-through biological nutrient removal (BNR) systems, with the goal of enabling scalable, nationwide adoption in activated sludge treatment processes. The project included demonstration of SBNR at the pilot scale as performed by Hampton Roads Sanitation District (HRSD) and at the full-scale as performed by the Los Angeles County Sanitation Districts' (LACSD) Pomona Water Reclamation Plant (POWRP).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗

DOE OSTI2023

Lime slurry treatment of soils developing on abandoned coal mine spoil: Linking contaminant transport from the micrometer to pedon-scale

Historical and abandoned coal mine spoil continues to generate acid- and metal(loid)-rich porewaters and represents a geographically large and diffuse non-point source of contamination to local watersheds. A potentially inexpensive approach to treat these materials and soils developing on them is through the application of lime slurries, to neutralize acidity and encourage the (co)precipitation of metal(loids) with Fe(III)-(oxy)hydroxides, and potentially other metal-oxides and/or Ca-bearing phases. Here, the efficacy of this approach was evaluated through parallel field application and laboratory-based flow-through column experiments. The field site is Huff Run sub-watershed 25 located in Tuscarawas County, Ohio and was chosen in part because it was previously classified as one of the most highly AMD-impacted sub-watersheds in the region. Two locations with historical spoil were chosen and suction lysimeters were installed at 25 and 75 cm depth to monitor porewater composition on two sides at the base of each pile. Half of each slope received seven lime slurry treatments from June through October of 2017. A suite of aqueous (ICP-OES, IC, and TOC-L) and solid phase geochemical and mineralogical approaches (quantitative SEM-EDS and synchrotron μ-XRF) were used to determine how composition, texture, morphology, and spatial distribution of mineral coatings differ in pre- and post-lime treated soils, and how that impacts the distribution and transport of trace metal(loid)s. Mine spoil porewater at site 1 was slightly less alkaline (pH ranging from 7.04 to 7.37) than at site 2 (ranging from 7.55 to 7.71), and average electrical conductivity values at site 1 (316–405 μS cm —1 ) were slightly lower than at site 2 (358–464 μS cm —1 ), although differences between the sites were not significant. Porewater pH and electrical conductivity in all lysimeters decreased over the course of the field season but there was no obvious response to lime treatment at either site or any depth. At site 1, both treatment and depth were significant factors affecting Ca, K, Ni, SO 4 2— , and DOC concentrations while only treatment effects were significant for dissolved Al and Cu (p < 0.05). For all soils, there were no trends in metal concentration observed over time although DOC and SO 4 2— decreased over the field season. Pedon-scale changes in metal porewater concentrations in response to treatment were linked to micrometer-scale changes in mineral surface coatings; specifically, higher concentrations of Ca, Fe, Mn, and Zn were observed in the coatings and no changes were observed in Fe redox speciation, whereas total S decreased likely due to oxidation of S in coal fragments. In contrast to the field experiment, the column experiments exhibited a much greater response in effluent composition with respect to lime treatment. The untreated columns had approximately an order of magnitude more H + leached over the course of the experiment (p < 0.001) and resulted in greater Ca, Al, Cu, and DIC leached and less Mn, Zn, and SO 4 2— . Soils treated with the lime slurry in the column experiments exhibited larger and thicker secondary Fe-coatings, including the addition of Fe-sulfates. Despite clear trends in the laboratory-based column experiment where the lime-to-soil ratio was higher, the effects were either muted or undetected in the field pilot project, suggesting that a higher application rate of lime in the field is needed to achieve a similar effect. This work provides evidence that a less alkaline lime slurry could be a practical and inexpensive method of treating coal mine spoil-impacted soils and represents an important step in linking laboratory-based remediation studies to implemented field-based studies.

58 GEOSCIENCES↗

DOE OSTI2022

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Quantifying Properties for a Mechanistic, Predictive Understanding of Aqueous Impact on Aging of Medium and Low Voltage AC and DC Cabling in Nuclear Power Plants

To address the gap in knowledge in understanding degradation in relevant service conditions, this project aimed to develop a mechanistic, predictive model of medium and low voltage cable failure based on the primary environmental degradation parameters of aqueous immersion time, temperature, and the oxidation extent. To do so, we took a two-fold approach toward evaluating degradation as related to the aqueous condition. First, we evaluated the chemical, mechanical, and electrical properties of polymers that comprise the cabling insulation under varied aqueous conditions at different temperatures. Secondly, we evaluated the performance of insulation materials under similar accelerated aging conditions in addition to the dielectric breakdown of the cabling under submersion conditions. Thermal oxidation and immersion of LDPE thin films were done using a Parr vessel and the extent of oxidation was monitored using the carbonyl index from ATR-FTIR spectrum, a ratio quantifying the carbonyl functionality produced on the LDPE film surface. Increasing the temperature, oxidation time, and oxygen pressure increased the water-vapor permeability and the carbonyl content. Thermogravimetric analysis of the films confirmed that there is indeed an initial weight loss of the smaller molecular weight molecules and volatiles. At higher temperatures, there is a secondary mass loss. Up to 70°C, the mass loss curves were similar until 80°C, when the aged films decreased at a much greater rate with respect to increasing temperature. This was seen in the 10, 50, 90% weight loss temperatures. The mechanical properties of sheet PP and HDPE as well as injection molded LDPE, HDPE, HDPP dog bones were plotted with respect to the time spent in the Parr vessel during thermal oxidation. For the sheet PP and HDPE, the UTS and the modulus of elasticity decreased while the elongation increased. The injection-molded LDPE, HDPE, and HDPP dog bones experienced a similar trend where the changes in the properties were within experimental uncertainty. This analysis of LDPE in dry and immersive oxygen-rich environments is an important baseline for future experiments looking at other degradation mechanisms. A predictive aging model for low-density polyethylene insulative cable housings was developed. This model revealed a combination of physical and chemical processes which lead to an increase in permeability and a decrease in strength of the insulator. In creating this predictive model, a gap of knowledge in the literature was revealed in two areas: understanding how the crystallinity of a polymer changes with time and visualizing the predicted pores formed through the insulator which leads to failure. Developing and using an adapted ATR-FTIR method, crystallinity was monitored and compared to bulk crystallinity found via DSC. Inhomogeneity in crystallinity changes were observed but the FTIR method but was found to be less reliable. Pore visualization was found utilizing electrical impedance spectroscopy saturated aged polyethylene films with synthesized citrate-capped gold nanoparticles, and chloroauric acid precursor. results indicate a decreasing impedance of the polyethylene films caused by the transport of ions through the film. SEM imaging was then performed for elemental analysis of the films and counter electrodes for the presence of ions. Chloride and gold ions were detected; however, no nanoparticles were found. This gives an estimated pore size of at least 0.3 nm through the aged film. Cyclic submergence of polyethylene and polypropylene in aqueous solutions of copper sulfate and Harrison’s solution were examined. It was discovered that aging in these mixed conditions and at 90°C did cause a small, yet significant increase to tensile strength for both unaged polyethylene and polypropylene. The PE and PP in cycled and submerged conditions did not have tensile strengths significantly different from dry-aged after 16 weeks. For both solutions, it was determined that capacitance increases with both water tree depth and cable temperature and is relatively independent of changing water tree AR. As for resistance, there is no apparent change between depth percentages of 10 and 80 for both solutions. This is because the water tree has not yet entered the conductor and thus there is no shorting current flowing through the tree yet. Between 80 and 100% water tree depth, it is observed that for both solutions there is a rapid decrease in resistance because the tree is now impacting the conductor shield and allowing shorting current to flow through it. The relationship between resistance and temperature was found to be inversed. As temperature increases, resistance decreases for both water tree solutions with distilled water having the most apparent difference across the range of temperatures simulated. With regards to water tree AR, resistance was found to be relatively independent at depth percentages of 10, 50, and 70, but there was some relationship at 100% depth for both solutions. At this depth, as water tree AR increased, resistance was found to also increase. Through the examination of voltage and electric field distribution plots, it was confirmed for EPR cables that increasing water tree depth leads to increasing distortion. The increase in shorting current at 100% depth for aqueous copper sulfate compared to distilled water was also visualized. Lastly, the relationship between shorting current and temperature for narrow water trees was visualized and validated using the distribution plots. This study has led to a better understanding of the effect of cable temperature, water tree solution, and geometry of the tree region on cable degradation. Furthermore, this study has documented simulation results of water tree degradation in EPR MV cables, which had been previously lacking in information. Future work in this area will add frequency into the mix and examine what effect it has on the rate of cable degradation as a result of water treeing. To compliment the work being carried by UMD on different polymer chemistries, harvested medium voltage (MV) cables were selected for accelerated aging study at ORNL. These medium voltage cables are representative of the vintages and materials that are currently in use at existing nuclear power plants (NPPs). These cables were immersed in water at 90°C and energized for a period of two years at elevated voltage. Partial discharge was measured during this period of time to track potential degradation. Unfortunately, due to the absence of failure, the extent of integration between the UMD techniques and the harvested MV insulation was limited. However, based on the findings and from the UMD research in the previous sections on polyethylene and polypropylene, follow-on characterization for harvested MV cable insulation should focus on sample preparation to take advantage of the UMD techniques to better understand mechanistic degradation in MV cable insulations. This would include: 1.) Utilization of solutions and impedance measurement insulation to track permeability in harvested and accelerated aging insulation samples, 2.) Adaptation of gold nanoparticle synthesis for electrical impedance spectroscopy with supporting SEM characterization to study pore formation in harvested and accelerated aged insulation samples, & 3.) ATR-FTIR crystallinity with multiple temperature dependent harvested MV cable insulation under different air and submerged accelerated aging. In addition, low voltage dielectric spectroscopy and high voltage dissipation factor, tan δ, could be implemented as a condition monitoring tool for harvested MV cable insulation in submerged environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗