Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “attainable yield”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

104 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Steady-state mixing state of black carbon aerosols from a particle-resolved model

Abstract. Black carbon (BC) exerts a notable warming effect due to its strong light absorption, largely influenced by its “mixing state”. However, due to computational constraints, the mixing state is challenging to accurately represent in large-scale models. In this study, we employ a particle-resolved model to simulate the evolution of BC mixing state based on field observation. Our result shows that aerosol compositions, coating thickness (CT) distribution, and optical properties of BC aerosols all exhibit a tendency toward a steady state with a characteristic timescale of less than 1 d, considerably shorter than the BC atmospheric lifetime. The rapid attainment of a steady state suggests that it is reasonable to disregard this pre-steady-state period and instead concentrate on the average properties of BC across extensive spatial and temporal scales. The distribution of CT follows an exponential linear distribution and can be characterized by a single slope parameter k. This distribution is independent of the BC core's distribution. In the model simulation, the mean CT, equivalent to the 1/k, is 62 nm, which is consistent with the statistical results indicating a mean CT of 63 nm. Utilizing the slope parameter k, which effectively characterizes the CT distribution under the steady-state simplifying assumption, the BC absorption enhancement closely corresponds to the results obtained via the particle-resolved method. This study simplifies the BC mixing state description and yields a precise evaluation of the BC optical properties, which has the potential utility for modeling efforts in the refinement of the assessment of BC's radiative effects.

Zhang, Zhouyang↗
DOE OSTI2021

Hydrogen Bonding and Its Effect on the Orientational Dynamics of Water Molecules inside Polyelectrolyte Brush-Induced Soft and Active Nanoconfinement

Despite being the most ubiquitous compound on Earth, the fundamental properties of water are not fully understood, especially in nanoconfinement. Densely grafted polyelectrolyte (PE) molecules attain the configuration of a “brush”: these PE brushes, due to their ability to form hydrogen bonds (HBs) with water via the PE functional groups, act as a source of soft and active nanoconfinement for the brush-trapped water molecules. In this paper, we study the effects of PE brush-induced confinement on the structure, dynamics and energetics of the water-water and water-PE HBs. Our results indicate a significant weakening of the HBs from bulk to sparsely grafted to densely grafted brushes. i.e., by increasing the degree of brush-induced nanoconfinement. We explain that this weakening of water-water HBs is caused by the disruption of the extended network of water molecules within the brush-induced nanoconfinement. This is confirmed by performing a ring structure analysis of the water molecules, which yields a reduction in the average ring size at higher degrees of brush-induced nanoconfinement (i.e., at higher brush grafting densities). Furthermore, we investigate the role of HB on the orientational dynamics of the water molecules. Here, we observe that the rotational motion of the water molecules becomes sluggish inside the PE brushes. Recent findings have indicated that the water and counterions trapped in brush-induced nanoconfinement demonstrate structures (in combination with the PE functional group) analogous to that in “water-in-salt” electrolytes that have seen extensive recent uses for Li-ion battery applications. However, the rotational dynamics of water molecules inside the brush-induced nanoconfinement is found to be distinctly different from that of conventional “water-in-salt” electrolytes in absence of any confinement; therefore, the present study will provide the necessary platform towards conceptualizing polymer-based nanoconfinement for battery applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Effect of Molecular Weight and Layer Thickness on Dielectric Breakdown Strength of Neat and Homopolymer Swollen Lamellar Block Copolymer Films

Designing next-generation lightweight pulsed power devices hinges on understanding the factors influencing the energy storage performance of dielectric materials. Polymer dielectric films have a quadratic dependence of energy storage on the voltage breakdown strength and strategies to enhance the breakdown strength are expected to yield a path toward high energy storage densities. Highly stratified lamellar block copolymer (L-BCP) films of model polystyrene-b-polymethylmethacrylate (PS-b-PMMA) exhibited as much as ~50% enhancement in breakdown voltage (E BD ), (225 % increase in stored energy density, U~ E BD 2 ), compared to unordered as-cast L-BCP films. Such an energy density using amorphous polymer is on par with industry-standard semi-crystalline biaxially oriented polypropylene (BOPP), and as such a notable development in the field. This work develops a deeper understanding of the molecular mechanisms of E BD enhancement in L-BCP films, due to combination of interface formation and spatial isolation of the chain ends into segregated zones. We further confirm the role of chain ends in the breakdown process blending a low M n L-BCP with matched M n homopolymers to attain same layer spacing as neat L-BCP of higher M n . E BD shows a significant decrease at low homopolymer fractions due to increased net chain end density within swollen ordered L-BCP domains in wet-brush regime, followed by increased E BD due to layer thickness increase via segregated “interphase layer” formation by excess homopolymers. Notably, E BD of homopolymer swollen L-BCPs is always lower than neat L-BCP of the same domain spacing, due to overall adverse chain ends contribution from homopolymers. These findings provide important selection rules for L-BCPs for designing next-generation flexible electronics with high energy density solid-state BCP film capacitors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Commodity Thermoplastic Elastomer-Enabled Templated Synthesis of Large-Pore Ordered Mesoporous Materials

Fabrication of ordered mesoporous materials (OMMs) has predominantly relied on templating-based methods. However, these methods are constrained by several limitations, especially the limited pore sizes attainable with commercially available surfactants used as structure-directing agents. To unlock the full potential of the OMMs, it is essential to develop synthetic strategies that facilitate the production of large-pore OMMs using scalable processes and cost-effective precursors. This work demonstrates the use of thermoplastic elastomer (TPE)-derived carbon replicas for synthesizing ordered mesoporous silica (OMS) and metal oxides (OMMOs) via precursor infiltration and template removal. The nanostructural evolution of the resulting inorganic materials was systematically investigated. Specifically, using tetraethyl orthosilicate (TEOS) as a silica precursor, this method can produce an OMS with relatively large pores. To establish the generalizability of this process, the fabrication approach was extended to other commercially available TPEs with varied chemical compositions and molecular weights while consistently resulting in ordered structures. Additionally, this synthetic strategy can be successfully applied to the production of OMMOs, including tin and titanium oxide matrix chemistries, yielding pore sizes of 16.0 and 19.2 nm, respectively. By developing a general method and using low-cost precursors, this work presents a scalable approach for fabricating large-pore OMMs with tunable pore textures and matrix chemistries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗
DOE OSTI2021

Experimental and theoretical determinations of hydrogen isotopic equilibrium in the system CH 4 —H 2 —H 2 O from 3 to 200°C

The stable isotopic composition of methane (CH 4 ) is commonly used to fingerprint natural gas origins. Over the past 50 years, there have been numerous proposals that both microbial and thermogenic CH 4 can form in or later attain hydrogen isotopic equilibrium with water (H 2 O) and carbon isotopic equilibrium with carbon dioxide (CO 2 ). Evaluation of such proposals requires knowledge of the equilibrium fractionation factors between CH 4 and H 2 O or CO 2 at the temperatures where microbial and thermogenic CH 4 form in or are found in the environment, which is generally less than 200°C. Experimental determinations of these fractionation factors are only available above 200°C, requiring extrapolation of these results beyond the calibrated range or the use of theoretical calculations at lower temperatures. Here, we provide a calibration of the equilibrium hydrogen isotopic fractionation factor for CH 4 and hydrogen gas (H 2 ) ( D α CH4(g)–H2(g) ) based on experiments using γ-Al 2 O 3 and Ni catalysts from 3 to 200°C. Results were regressed as a 2 nd order polynomial of 1000 × ln D α CH4(g)–H2(g) vs. 1/T (K -1 ) yielding: 1000 × l n D α C H 4 ( g ) - H 2 ( g ) = 3.5317 × 10 7 T 2 + 2.7749 × 10 5 T - 179.48 We combine this calibration with previous experimental determinations of hydrogen isotope equilibrium between H 2 , H 2 O(g), and H 2 O(l) and we provide an interpolatable experimental calibration of 1000 × ln D α CH4(g)–H2O(l) from 3 to 200°C. Our resulting 4th order polynomial is the following equation: 1000 × l n D α C H 4 ( g ) - H 2 O l = - 7.9443 × 10 12 T 4 + 8.7772 × 10 10 T 3 - 3.4973 × 10 8 T 2 + 5.4398 × 10 5 T - 382.05 At 3°C, the value from our calibration differs by 93‰ relative to what would be calculated based on the extrapolation of the only experimental calibration currently available to temperatures below its calibrated range (lowest temperature of 200°C; Horibe and Craig, 1995). We additionally provide new theoretical estimates of hydrogen isotopic equilibrium between CH 4 (g), H 2 (g), and H 2 O(g) and carbon isotopic equilibrium between CH 4 (g) and CO 2 (g) using Path Integral Monte Carlo (PIMC) calculations. Our PIMC calculations for hydrogen isotopic equilibrium between CH 4 and H 2 agree 1:1 with our experiments. Finally, we compile carbon and hydrogen isotopic measurements of CH 4 , CO 2 , and H 2 O from various environmental systems and compare observed differences between carbon and hydrogen isotopes to those expected based on isotopic equilibrium. We find that isotopic compositions of some microbial gases from marine sedimentary, coalbed, and shale environments are consistent with those expected for CH 4 H 2 O(l) hydrogen and CH 4 CO 2 carbon isotopic equilibrium. In contrast, microbial terrestrial and pure culture gases are not consistent with both CH 4 H 2 O(l) hydrogen and CH 4 CO 2 carbon isotopic equilibrium. Overall, these results are explained qualitatively using previously developed conceptual models that link free energy gradients available to microorganisms to the degree that their enzymes can promote isotope-exchange reactions between CH 4 , CO 2 , and H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Geostrophic to magnetically-damped turbulence in liquid metal rotating magnetoconvection

Understanding planetary core convection dynamics requires the study of convective flows in which the Coriolis and Lorentz forces attain a leading-order, so-called magnetostrophic balance. Experimental investigations of rotating magnetoconvection (RMC) in the magnetostrophic regime are therefore essential to broadly characterize the properties of local-scale planetary core flow. Toward this end, we present here the first thermovelocimetric measurements of magnetostrophic, liquid metal convection, which are made using liquid gallium as the working fluid, at moderate rotation rates (Ekman numbers 10 -4 ≤ Ek ≤ 10 -5 ) and in the presence of dynamically strong magnetic fields (Elsasser number Λ = 1). Complementary rotating convection (RC) experiments are performed at the same rotation rates to serve as reference cases. Our RMC velocity measurements adequately follow a geostrophic turbulent scaling for cases in which local-scale convective inertial forces exceed the Lorentz forces in the fluid bulk. In cases where Lorentz forces exceed local-scale inertia (local or convective interaction parameter N ℓ or N C ≳ 3), the root-mean-square RMC velocities are magnetically damped, yielding values below the geostrophic turbulent RC scaling prediction. An enhancement in heat transfer is observed, which we attribute to the increased coherence of vertically aligned magnetostrophic convective flow. Finally, our results and comparisons to dynamo simulations suggest that the system is in the strong-field regime, where convection-scale flows are magnetically damped with N C ≳ 3).

Convective turbulence↗
DOE OSTI2020

Wave Prediction using X-Band Radar (Final Summary Report)

Model Predictive Control (MPC) is the best controls framework available to wave energy today. It allows for the maximization of energy yield from a WEC system while respecting various constraints and losses in the system. This type of constrained optimization is a critical capability for the techno-economic optimization of WEC devices and is difficult to do effectively with causal controls frameworks. MPC requires a future prediction of wave excitation forces, requiring wave prediction on a future time horizon of up to 30 seconds. The future horizon required is highly dependent on the WEC and PTO type available. To our knowledge no-one has been able to implement MPC on a WEC device at sea, due to the fact that phase-resolved wave-prediction is not a capability that has been sufficiently developed to date. The key objectives of this project are focused on developing “industry-ready” wave prediction technology building blocks that can be applied to a wide range of different wave energy conversion topologies. Our collaborative work with device developers has demonstrated that the level of improvement attainable for any particular device is heavily dependent on the device and PTO topology chosen. The range of annual average power capture improvements is on the order of 10% to over 100%. However, the important aspect is that we can get within about 10% of theoretical upper limits even when considering errors in the wave prediction and more importantly, we can optimize performance while respecting various constraints such as motion amplitude or peak structural loads. The combination of improved performance while limiting structural loads has a net effect of reducing the levelized cost of energy from these emerging technologies and allows us to design control laws that meet an economic optimum. Efforts under a previous project on controls has focused on developing a set of wave-prediction algorithms that work by leveraging a network of wave sensing measurement buoys that are equipped with real-time telemetry. The present work is focusing primarily on using X-band radar to measure the wave-field around the WEC to predict the future excitation forces on the structure and combining the measurements of these different sources to create improvements in the wave-prediction accuracy.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Application of Electromagnetic Methods for Reservoir Monitoring with Emphasis on Carbon Capture, Utilization, and Storage

The Controlled-Source ElectroMagnetic (CSEM) method provides crucial information about reservoir fluids and their spatial distribution. Carbon dioxide (CO2) storage, enhanced oil recovery (EOR), geothermal exploration, and lithium exploration are ideal applications for the CSEM method. The versatility of CSEM permits its customization to specific reservoir objectives by selecting the appropriate components of a multi-component system. To effectively tailor the CSEM approach, it is essential to determine whether the primary target reservoir is resistive or conductive. This task is relatively straightforward in CO2 monitoring, where the injected fluid is resistive. However, for scenarios involving brine-saturated (water-wet) or oil-wet (carbon capture, utilization, and storage—CCUS) reservoirs, consideration must also be given to conductive reservoir components. The optimization of data acquisition before the survey involves analyzing target parameters and the sensitivity of multi-component CSEM. This optimization process typically includes on-site noise measurements and 3D anisotropic modeling. Based on our experience, subsequent surveys tend to proceed smoothly, yielding robust measurements that align with scientific objectives. Other critical aspects to be considered are using magnetotelluric (MT) measurements to define the overall background resistivities and integrating real-time quality assurance during data acquisition with 3D modeling. This integration allows the fine tuning of acquisition parameters such as acquisition time and necessary repeats. As a result, data can be examined in real-time to assess subsurface information content while the acquisition is ongoing. Consequently, high-quality data sets are usually obtained for subsequent processing and initial interpretation with minimal user intervention. The implementation of sensitivity analysis during the inversion process plays a pivotal role in ensuring that the acquired data accurately respond to the target reservoirs’ expected depth range. To elucidate these concepts, we present an illustrative example from a CO2 storage site in North Dakota, USA, wherein the long-offset transient electromagnetic method (LOTEM), a variation of the CSEM method, and the MT method were utilized. This example showcases how surface measurements attain appropriately upscaled log-scale sensitivity. Furthermore, the sensitivity of the CSEM and MT methods was examined in other case histories, where the target reservoirs exhibited conductive properties, such as those encountered in enhanced oil recovery (EOR), geothermal, and lithium exploration applications. The same equipment specifications were utilized for CSEM and MT surveys across all case studies.

Barajas-Olalde, César↗
DOE OSTI2021

Super-Resolution Photothermal Patterning in Conductive Polymers Enabled by Thermally Activated Solubility

Doping induced solubility control (DISC) patterning is a recently developed technique which uses the change in polymer solubility upon doping, along with an optical dedoping process, to achieve high resolution optical patterning. DISC patterning is remarkable because the process resolution exceeds the linear Abbe’s diffraction limit, however, no mechanism has been proposed to explain such high resolution. Here, we use diffraction to spatially modulate the light intensity and determine the dissolution rate, revealing a superlinear dependence on light intensity, which yields features that are sharper than the profile of the writing photomask. This rate law is independent of wavelength, indicating that patterning resolution is not dominated by an optical dedoping reaction, as was previously proposed. Instead we show here that the optical patterning mechanism is primarily controlled by the thermal profile generated by the laser. To quantify this effect, the thermal profile and dissolution rate are modeled using a finite-element model and compared against patterned line cross sections as a function of wavelength, laser intensity, and dwell time. Our model reveals that although the laser-generated thermal profile is broadened considerably beyond the lasers resolution, the highly temperature dependent dissolution rate results in selective dissolution near the peak of the thermal profile. Therefore, the key factor in achieving super-resolution patterning is a strongly temperature dependent dissolution rate|a common feature of many polymers. In addition to suggesting several routes to improved resolution, our model also demonstrates that doping is not required for optical patterning of conjugated polymers, as was previously believed. Instead, we demonstrate that superlinear resolution optical patterning is attainable in any conjugated polymer simply by tuning the solvent quality during patterning, thus extending the applicability of our method to a wide class of materials. Here, we demonstrate the generality of photothermal patterning by writing sub-400 nm features into undoped Pff-BT4T-2OD.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2022

An all optical approach for comprehensive in-operando analysis of radiative and nonradiative recombination processes in GaAs double heterostructures

We demonstrate an all optical approach that can surprisingly offer the possibility of yielding much more information than one would expect, pertinent to the carrier recombination dynamics via both radiative and nonradiative processes when only one dominant deep defect level is present in a semiconductor material. By applying a band-defect state coupling model that explicitly treats the inter-band radiative recombination and Shockley–Read–Hall (SRH) recombination via the deep defect states on an equal footing for any defect center occupation fraction, and analyzing photoluminescence (PL) as a function of excitation density over a wide range of the excitation density (e.g., 5–6 orders in magnitude), in conjunction with Raman measurements of the LO-phonon plasmon (LOPP) coupled mode, nearly all of the key parameters relevant to the recombination processes can be obtained. They include internal quantum efficiency (IQE), minority and majority carrier density, inter-band radiative recombination rate (W r ), minority carrier nonradiative recombination rate (W nr ), defect center occupation fraction (f), defect center density (N t ), and minority and majority carrier capture cross-sections (σ t and σ tM ). While some of this information is thought to be obtainable optically, such as IQE and the W r /W nr ratio, most of the other parameters are generally considered to be attainable only through electrical techniques, such as current-voltage (I-V) characteristics and deep level transient spectroscopy (DLTS). Following a procedure developed herein, this approach has been successfully applied to three GaAs double-heterostructures that exhibit two distinctly different nonradiative recombination characteristics. The method greatly enhances the usefulness of the simple PL technique to an unprecedented level, facilitating comprehensive material and device characterization without the need for any device processing.

14 SOLAR ENERGY↗
DOE OSTI2021

Edge-Propagation Discharge Mechanism in CFx Batteries—A First-Principles and Experimental Study

Graphite fluoride (CF x ) cathodes coupled with lithium anodes yield one of the highest theoretical specific capacities (>860 mAh/g) among primary batteries. In practice, the observed discharge voltage (~2.5 V) is significantly lower than thermodynamic limits (>4.5 V), the discharge rate is low, and so far Li/CF x has only been used in primary batteries. Understanding the discharge mechanism at atomic length scales will improve practical CFx energy density, rate capability, and rechargeability. So far, purely experimental techniques have not identified the correct discharge mechanism or explained the discharge voltage. In this work, we apply density functional theory calculations to demonstrate that a CF x -edge propagation discharge mechanism based on lithium insertion at the CF/C boundary in partially discharged CF x exhibits a voltage range of 2.5 to 2.9 V—depending on whether solvent molecules are involved. The voltages and solvent dependence agree with our discharge and galvanostatic intermittent titration technique measurements. The predicted discharge kinetics are consistent with CF x operations. Finally, we predict some Li/CF x rechargeability under the application of high potentials, along a charging pathway different from that of discharge. Our work represents a general, quasi-kinetic framework to understand the discharge of conversion cathodes, circumventing the widely used phase diagram approach which most likely does not apply to Li/CF x because equilibrium conditions are not attained in this system.

25 ENERGY STORAGE↗
DOE OSTI2020

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗