SEARCH · Engineering Papers

Results for “Pressure Core”

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.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Synchrotron X-ray metrology of dopant distribution and oxidation state in high pressure CVD grown TM 2+ :ZnSe optical fibers

High pressure chemical vapor deposition (HPCVD) has shown to be a promising method for producing Cr 2+ :ZnSe and Fe 2+ :ZnSe mid-IR optical fiber lasers. The oxidation state of the dopants and their spatial homogeneity are critical for the characterization of effective fiber laser performance; however, this is challenging for small fiber cores and small doping concentrations. We demonstrate using synchrotron micro X-ray fluorescence (XRF) mapping to study the distribution of the dopants throughout the fiber cross-sections. Furthermore, we study the local chemical environment of these unique fiber structures using micro X-ray near edge absorption spectroscopy (XANES). Our study reveals that transition metal doped ZnSe deposited using HPCVD has nanoscale dopant aggregation although the material is chemically identical to the commercially produced diffusion doped laser crystals, as indicated by the presence +2 oxidation state of the dopants. This work demonstrates the power of synchrotron-based X-ray techniques for use as metrology tools for improving the performance of laser materials in fiber geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Water Formation Reaction under Interfacial Confinement: Al0.25Si0.75O2 on O-Ru(0001)

Confined nanosized spaces at the interface between a metal and a seemingly inert material, such as a silicate, have recently been shown to influence the chemistry at the metal surface. In prior work, we observed that a bilayer (BL) silica on Ru(0001) can change the reaction pathway of the water formation reaction (WFR) near room temperature when compared to the bare metal. In this work, we looked at the effect of doping the silicate with Al, resulting in a stoichiometry of Al0.25Si0.75O2. We investigated the kinetics of WFR at elevated H2 pressures and various temperatures under interfacial confinement using ambient pressure X-ray photoelectron spectroscopy. The apparent activation energy was lower than that on bare Ru(0001) but higher than that on the BL-silica/Ru(0001). The apparent reaction order with respect to H2 was also determined. The increased residence time of water at the surface, resulting from the presence of the BL-aluminosilicate (and its subsequent electrostatic stabilization), favors the so-called disproportionation reaction pathway (*H2O + *O ↔ 2 *OH), but with a higher energy barrier than for pure BL-silica.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The Surface Chemistry of Methanol on Cu 3 Pd(111): Effects of Metal Alloying and Reaction with Hydrogen

Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS) was used to study the adsorption and surface chemistry of methanol on a Cu 3 Pd(111) model surface. The composition and morphological properties of the pristine Cu 3 Pd(111) substrate were analyzed using a combination of low-energy electron diffraction (LEED), scanning tunneling microscopy (STM), and low-energy ion scattering (LEIS). The results of ion scattering showed segregation of Pd toward the surface, with a Pd/Cu ratio close to 0.5, a larger value than the ratio of 0.33 expected for a Cu 3 Pd bimetallic structure. The surface of the alloy exhibited a good LEED pattern with three-fold long-range periodicities. In STM, clusters of palladium with hexagonal arrays and Pd- Pd distances of 2.7-2.8 Å were detected. Bonding to copper perturbed the electronic properties of the atoms in the Pd clusters, shifting their 4d states toward higher binding energy with respect to the Fermi level. Further, the valence band spectrum of Cu 3 Pd(111) exhibited a line shape that was very different from those displayed by Cu(111) or Pd(111). At low pressures, the adsorption of methanol on Cu 3 Pd(111) at 300 K mainly produced CH 3 O, CO and CH x species. AP-XPS showed that most Pd atoms in the surface of the Cu 3 Pd(111) alloy interacted with the decomposition products of methanol. No significant changes were observed in the core levels of copper upon the adsorption and dissociation of methanol, suggesting that the molecule mainly interacted with Pd sites of the alloy. Reaction with hydrogen led to fast removal of CH x , C and PdC x species from Cu 3 Pd(111) at moderate (< 450 K) temperatures and prevented a CH x → C transformation. If the stability of adsorbed CH 3 O is used as a descriptor for the hydrogenation of CO 2 to methanol, Cu 3 Pd should be a much better catalyst than monometallic palladium. This may be a consequence of electronic and ensemble effects in the alloy that moderate the reactivity of Pd sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Initial TCV operation with a baffled divertor

The Tokamak à Configuration Variable (TCV) tokamak is in the midst of an upgrade to further its capability to investigate conventional and alternative divertor configurations. To that end, modular and removable gas baffles have been installed to decrease the coupling between the divertor and the plasma core. The baffles primarily seek to suppress the transit of recycling neutrals to closed flux surfaces. A first experimental campaign with the gas baffles has shown that the baffled divertor remains compatible with a wide range of configurations including snowflake and super-X divertors. Plasma density ramp experiments reveal an increase of the neutral pressure in the divertor by up to a factor ×5 compared to the unbaffled divertor and thereby qualitatively confirm simulations with the SOLPS-ITER code that were used to guide the baffle design. Finally, together with a range of new and upgraded divertor diagnostics, the baffled TCV divertor is now used to validate divertor models for ITER and next step devices with particular emphasis on geometric variations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Overview of results from the 2023 DIII-D negative triangularity campaign

Negative triangularity (NT) is a potentially transformative configuration for tokamak-based fusion energy with its high-performance core, edge localized mode (ELM)-free edge, and low-field-side divertors that could readily scale to an integrated reactor solution. Previous NT work on the TCV and DIII-D tokamaks motivated the installation of graphite-tile armor on the low-field-side lower outer wall of DIII-D. A dedicated multiple-week experimental campaign was conducted to qualify the NT scenario for future reactors. During the DIII-D NT campaign, high confinement ($H_\mathrm{98y,2}\gtrsim$ 1), high current ($q_\mathrm{95} \lt$ 3), and high normalized pressure plasmas ($\beta_\mathrm{N} \gt$ 2.5) were simultaneously attained in strongly NT-shaped discharges with average triangularity $\delta_\mathrm{avg}$ = −0.5 that were stably controlled. Experiments covered a wide range of DIII-D operational space (plasma current, toroidal field, electron density and pressure) and did not trigger an ELM in a single discharge as long as sufficiently strong NT was maintained; in contrast, to other high-performance ELM-suppression scenarios that have narrower operating windows. These strong NT plasmas had a lower outer divertor X-point shape and maintained a non-ELMing edge with an electron temperature pedestal, exceeding that of typical L-mode plasmas. Also, the following was achieved during the campaign: high normalized density ($n_\mathrm{e}$/$n_\mathrm{GW}$ of at least 1.7), particle confinement comparable to energy confinement with $Z_\mathrm{eff}\sim2$ , a detached divertor without impurity seeding, and a mantle radiation scenario using extrinsic impurities. These results are promising for a NT fusion pilot plant but further questions on confinement extrapolation and core-edge integration remain, which motivate future NT studies on DIII-D and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗

DOE OSTI2026

High-Burnup LOCA Burst Susceptibility BISON Analysis in PWRs and BWRs

Accurately assessing high-burnup fuel behavior during loss-of-coolant accidents (LOCAs) is essential for understanding fuel fragmentation, relocation, and dispersal (FFRD) risks across the US light-water reactor fleet. This work updates previous Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program multiphysics LOCA analyses for a pressurized water reactor (PWR) and a boiling water reactor (BWR) by incorporating recent model and material property advancements in the BISON fuel performance code, including a high-burnup structure (HBS) model, revised cladding burst criteria, and updated thermal–mechanical correlations. This update was needed to support ongoing industry initiatives and upcoming regulatory changes. Full-core, rod-resolved operating histories generated using Virtual Environment for Reactor Analysis (VERA) and system-level LOCA conditions obtained from TRACE were applied to statistically representative rod samples in BISON to evaluate burst behavior and FFRD susceptibility. These calculations used two cladding burst correlations and three fuel pulverization models so that the predictions of these models could be compared. The updated PWR simulations show markedly improved numerical stability as the number of crashed simulations decreased by 95% compared to the previous study, and hence higher confidence in results. The updated PWR simulations predicted cladding bursts exclusively among once-burned, high-power rods, with two different cladding burst models identifying the same burst-susceptible population. Resulting FFRD susceptibility estimates are significantly reduced compared with earlier studies, driven by cooler predicted fuel and plenum temperatures, lower hoop strains, and reduced fission gas release in the updated models. In contrast, none of the BWR rods were predicted to burst under either burst criterion, reaffirming minimal BWR FFRD susceptibility even with updated HBS and material models. Comparisons between the PWR and BWR end-of-cycle predictions are made. Comparison with prior work highlights significant shifts in PWR fuel performance metrics and confirmation of earlier BWR conclusions. Overall, the updated results underscore the importance of having high-resolution detailed modeling capability and continuously integrating evolving material models and physics into high-resolution multiphysics simulations. The unified assessment presented here strengthens confidence in predicting high-burnup LOCA behavior by improving agreement between different cladding burst correlations. These results also provide an improved foundation for future BISON model development, FFRD susceptibility calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗

DOE OSTI2023

Thermal Hydraulic Modeling of an Advanced Nuclear Reactor using open-source MOOSE tools

The aim of this study is to develop model of a small High Temperature Gas cooled Reactor (HTGR) including the balance of plant. This microreactor produces electricity using the thermal power of the nuclear reaction. This work utilizes the MOOSE Multiphysics simulation tools, which are mainly developed at Idaho National Laboratory (INL, Idaho, United States of America). Its thermal-hydraulics and Heat conduction Modules are used to study the fluids behavior in the primary loop and power conversion system and their interactions with the heating structures. More specifically, one verifies that the temperatures, pressures and mass flow rates of the fluids in both loops are consistent and that, at the same time, all the power transfers occur as expected. Moreover, the various mechanical components characteristics (turbine, compressor, pump) are adapted to the reactor operating conditions. This first analysis is conducted using simplified and one-dimensional model for the core. Ultimately, the obtained results are used to build a higher fidelity core model. This step aims to verify that the initial simplified simulation is consistent with the three-dimensional core modeling. In addition, the calculated material temperatures are checked. This study provides a fairly complete model of the thermal hydraulic phenomena for a particular design of High Temperature Gas cooled Reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗

DOE OSTI2026

Optimizing pressurized-water reactor equilibrium cycle using a novel loading pattern encoding and rule-based genetic crossover operators

This work presents an extended multi-batch approach applied in shuffling scheme optimization for equilibrium cycle for pressurized water reactors using Genetic Algorithms (GAs). A new ruled based GA crossover operator called Inherited Location and Batch (ILB) was introduced to enhance offsprings reproduction efficiency specialized for equilibrium cycle optimization problem. This approach was implemented within the Plant ReLoad Optimization (PRLO) framework and validated using a generic reactor model based on the AP1000 design, with core parameters calculated via the CASMO/SIMULATE software package. The ILB approach is then applied for both single and multi-objective problems in maximizing cycle length and core average exposure while minimizing the average enrichment of the 57 fresh fuel assemblies (FAs) per cycle. The optimal solutions are selected based on their dominance to the objectives from all feasible solutions. This research identified three optimal solutions satisfied safety constraints: The first solution minimizes feed enrichment costs with a cycle length of 338.8 days and core exposure of 25.39 MWd/MT; the second solution extends cycle length to 361.2 days, with the highest core exposure of 26.84 MWd/MT, using 3.75 wt% average fuel enrichment; the third solution balances both objectives with a cycle length of 349.6 days, core exposure of 25.82 MWd/MT with a slight enrichment increase compared to the first solution. Collectively, these findings underscore the efficiency and effectiveness of the proposed approach in achieving practical multi-objective optimal equilibrium cycle designs using GAs optimizer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

RELAP5-3D modelling of the HTTR-GT/H 2 secondary system and turbomachinery

High Temperature Gas-cooled Reactors (HTGRs) can be used for the cogeneration of electricity and hydrogen. The JAEA-operated High-Temperature engineering Test Reactor (HTTR-GT/H 2 ) is exploring using the reactor for electricity and hydrogen production. Here, a RELAP5-3D model of the HTTR-GT/H 2 secondary system has been developed using design information. The various components and heat exchangers in the secondary system were investigated, modified, and improved where applicable. With these changes made, the model was run, and the results were compared to the design conditions. The results for the sole-power generation mode were shown to fit the design conditions very well. The largest temperature difference was on the order of 7 K, and the largest pressure difference was on the order of 0.05 MPa. The results for the hydrogen cogeneration mode did not fit the design conditions as well as the other case. The largest temperature difference was about 39 K and the largest pressure difference was about 0.27 MPa at the compressor outlet. The larger differences for the hydrogen cogeneration mode are attributed to the various complex components and the flow being split in the secondary loop. A transient reduction in heat removal capability of the secondary system was investigated. Reactor temperatures were predicted to rise as a result. The core reactivity response due to this increase in temperature is investigated and is expected to add negative reactivity to the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Multiscale Assessment of Caprock Integrity for Geologic Carbon Storage in the Pennsylvanian Farnsworth Unit, Texas, USA

Leakage pathways through caprock lithologies for underground storage of CO2 and/or enhanced oil recovery (EOR) include intrusion into nano-pore mudstones, flow within fractures and faults, and larger-scale sedimentary heterogeneity (e.g., stacked channel deposits). To assess multiscale sealing integrity of the caprock system that overlies the Morrow B sandstone reservoir, Farnsworth Unit (FWU), Texas, USA, we combine pore-to-core observations, laboratory testing, well logging results, and noble gas analysis. A cluster analysis combining gamma ray, compressional slowness, and other logs was combined with caliper responses and triaxial rock mechanics testing to define eleven lithologic classes across the upper Morrow shale and Thirteen Finger limestone caprock units, with estimations of dynamic elastic moduli and fracture breakdown pressures (minimum horizontal stress gradients) for each class. Mercury porosimetry determinations of CO2 column heights in sealing formations yield values exceeding reservoir height. Noble gas profiles provide a “geologic time-integrated” assessment of fluid flow across the reservoir-caprock system, with Morrow B reservoir measurements consistent with decades-long EOR water-flooding, and upper Morrow shale and lower Thirteen Finger limestone values being consistent with long-term geohydrologic isolation. Together, these data suggest an excellent sealing capacity for the FWU and provide limits for injection pressure increases accompanying carbon storage activities.

58 GEOSCIENCES↗

DOE OSTI2025

Black Hole Supernovae, Their Equation of State Dependence, and Ejecta Composition

Abstract Recent literature on core-collapse supernovae suggests that a black hole (BH) can form within ∼1 s of shock revival, while still culminating in a successful supernova. We refer to these as BH supernovae, as they are distinct from other BH formation channels in both timescale and impact on the explosion. We simulate these events self-consistently from core collapse until 20–50 days after collapse using three axisymmetric models of a 60 M ⊙ zero-age main-sequence progenitor star and investigate how the composition of the ejecta is impacted by the BH formation. We employ Skyrme-type equations of state (EOSs) and vary the uncertain nucleonic effective mass, which affects the pressure inside the proto–neutron star through the thermal part of the EOS. This results in different BH formation times and explosion energies at BH formation, yielding final explosion energies between 0.06 and 0.72 × 10 51 erg with 21.8–23.3 M ⊙ of ejecta, of which 0–0.018 M ⊙ is 56 Ni. Compared to expectations from 1D simulations, we find more nuanced EOS dependences of the explosion dynamics, the mass of the BH remnant, and the elemental composition of the ejecta. We investigate why the explosions survive despite the massive overburden and link the shape of the diagnostic energy curve and character of the ejecta evolution to the progenitor structure.

Eggenberger Andersen, Oliver (ORCID:00000002966079↗

DOE OSTI2021

System Modeling of the HTTR and Economic Dispatch Model of the Secondary System

High Temperature Gas-cooled Reactors (HTGRs) can be used for the generation of electricity and their process heat can be used to improve the efficiency of chemical processes such as hydrogen production. The JAEA-operated High-Temperature engineering Test Reactor (HTTR-GT/H2) is exploring using the reactor for electricity and hydrogen production. A RELAP5-3D model of the HTTR-GT/H2 secondary system has been developed using design information. The various components and heat exchangers in the secondary system were modeled and results were compared to the design conditions. The results for the sole-power generation mode were shown to fit the design conditions very well. The largest temperature difference was on the order of 7 K, and the largest pressure difference was on the order of 0.05 MPa. The results for the hydrogen cogeneration mode did not match the design conditions nearly as well. The largest temperature difference was about 39 K and the largest pressure difference was about 0.27 MPa at the compressor outlet. The larger differences for the hydrogen cogeneration mode are attributed to the various complex components and the flow being split in the secondary loop. A transient reduction in heat removal capability of the secondary system was investigated. Reactor temperatures are anticipated to rise as a result. The core reactivity response due to this increase in temperature is investigated and is expected to add negative reactivity to the reactor. An economic dispatch model was developed for a nuclear-driven iodine-sulfur cycle system to determine hydrogen sale prices that would make such a system profitable. The study focuses on the development of the economic model and the role that input data plays on final calculated values. It was found that the input electricity prices, whether using historical data or a host of synthetic time histories, produce significantly different breakeven hydrogen sale prices. As such, great care should be used in these economic dispatch analyses to select reasonable input assumptions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Beyond the Surface: Non-Invasive Low-Field NMR Analysis of Microbially-Induced Calcium Carbonate Precipitation in Shale Fractures

Microbially-induced calcium carbonate precipitation (MICP) is a biological process in which microbially-produced urease enzymes convert urea and calcium into solid calcium carbonate (CaCO 3 ) deposits. MICP has been demonstrated to reduce permeability in shale fractures under elevated pressures, raising the possibility of applying this technology to enhance shale reservoir storage safety. For this and other applications to become a reality, non-invasive tools are needed to determine how effectively MICP seals shale fractures at subsurface temperatures. In this study, two different MICP strategies were tested on 2.54 cm diameter and 5.08 cm long shale cores with a single fracture at 60 °C. Flow-through, pulsed-flow MICP-treatment was repeatedly applied to Marcellus shale fractures with and without sand (“proppant”) until reaching approximately four orders of magnitude reduction in apparent permeability, while a single application of polymer-based “immersion” MICP-treatment was applied to an Eagle Ford shale fracture with proppant. Low-field nuclear magnetic resonance (LF-NMR) and X-Ray computed microtomography (micro-CT) techniques were used to assess the degree of biomineralization. With the flow-through approach, these tools revealed that while CaCO 3 precipitation occurred throughout the fracture, there was preferential precipitation around proppant. Without proppant, the same approach led to premature sealing at the inlet side of the core. In contrast, immersion MICP-treatment sealed off the fracture edges and showed less mineral precipitation overall. This study highlights the use of LF-NMR relaxometry in characterizing fracture sealing and can help guide NMR logging tools in subsurface remediation efforts.

MICP↗

DOE OSTI2024

On discriminating tokamak sawtooth crash models via localized density and temperature measurements

The core electron temperature drops rapidly during the sawtooth crash in tokamak plasmas, which causes heat loss and may lead to fast particle losses or even a disruption. Several models have been proposed for the periodic crash, including the Kadomtsev model with magnetic reconnection and the quasi-interchange model with the growth of higher-mode-number pressure-driven instabilities. 3D MHD simulations were performed for these two models with a goal to develop intuition and to predict qualitatively how different types of sawtooth will appear in various diagnostics. The structures of electron density ne and electron temperature Te show a dominant (1, 1) mode for the Kadomtsev case and a dominant (4, 4) mode for the quasi-interchange case. The oscillations of ne and Te have a positive correlation near the inversion layer for both cases, while their frequencies and amplitudes are different depending on the dominant modes. Particularly, for the Kadomtsev case, we find a relation between the amount of flux reconnected during a sawtooth event and ne or Te oscillations. Therefore, we connect recently developed measurement capabilities for ne and Te to the internal sawtooth behavior. We propose that this method of analysis can help in identifying the type of sawtooth in future experiments augmented by simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗