Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “favorability”

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

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

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Pressure-induced phase transformations in Nb 2 O 5 from ab initio molecular dynamics simulations

Here, the equations of state for the H-, L- and T-Nb 2 O 5 polymorphs were investigated with density functional theory. The shock Hugoniot curves were predicted from ab initio molecular dynamics simulations along 300–2500 K isotherms. Enthalpies derived from simulations at 300 K show that the H → L transition occurs at 6.2 GPa, in agreement with the measured value of ~5.9 GPa, and is more favorable than the H → T transition. Pressure-induced amorphization (PIA) of L-Nb 2 O 5 is in the range ~18-25 GPa, consistent with experiments reporting PIA above 21.4 GPa. The simulated PIA of T-Nb 2 O 5 between ~16 and 26 GPa reproduces PIA observed above 19.2 GPa.

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

Advanced supercritical water-based process concepts for treatment and beneficial reuse of brine in oil/gas production

Produced water generated from oil/gas reservoirs is a cost burden for oil/gas producers, with over 21 billion barrels of the waste generated in 2019. Average disposal costs range from 4.00–8.00·bbl -1 , with up to 20·bbl -1 spent on transportation alone. Dissolved solids content is a contributing factor which limits economical treatment options for this waste stream. This work proposes a novel technique that can handle high salinity waste, employing favorable properties past the critical point of water. In this study, process simulation of two supercritical water desalination (SCWD) scenarios was completed. Zero liquid discharge (ZLD) and brine concentration operating scenarios were compared, weighing the associated economics and benefits for each case. The results were shown to be economically feasible for brines with a high dissolved solids content, ranging from 3.49 to 17.28·bbl -1 in an expanded sensitivity analysis.

42 ENGINEERING↗
DOE OSTI2025

Stirring and peristaltic pumping alter flow electrode particle size and morphology

Flowable carbon-based electrodes are widely studied in electrochemical systems for energy storage and water treatment. Most laboratory- and pilot-scale studies employ active mixing for slurry preparation and handling and use peristaltic pumps for slurry circulation. Magnetic stirring is the most prevalent mixing method at the laboratory scale, and the combination of magnetic mixing and peristaltic pumping is favored for convenience and compatibility with particulate suspensions. However, these methods can strongly affect particle size and morphology (and therefore cell performance). In this study, we evaluate how mixing and pumping influence particle degradation in two commonly used carbon materials: spherical carbon black (CB) beads and powdered activated carbon (AC). We performed 24 h stirring-only experiments using magnetic, orbital, and overhead stirring. Separately, we conducted 32 h circulation experiments using a commercial peristaltic pump. Particle morphology and size were analyzed using optical microscopy and laser diffraction. The initially spherical CB beads (~580 μm, 1.5 wt%) exhibit rapid and extensive breakdown under both peristaltic pumping and magnetic stirring. By comparison, AC particles (initially ~31 μm) tested at 1.5 and 15 wt% showed measurable size reductions of up to ~25 % under either stirring or pumping. In conclusion, these findings highlight the need to account for and quantify morphological changes caused by both stirring and pumping in the design and associated evaluation of flowable electrode systems.

Chemistry↗
DOE OSTI2025

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2022

Capture of novel sp 3 hybridized Z-BN by compressing boron nitride nanotubes with small diameter

Experimental synthesis of new sp 3 hybridized carbon/boron nitride structures remains challenging despite that numerous sp 3 structures have been proposed in theory. Here, we showed that compressed multi-walled boron nitride nanotubes (MWBNNTs) and boron nitride peapods (C 60 @BNNTs) with small diameters could transform into a new sp 3 hybridized boron nitride allotrope (Z-BN). This strategy is considered from the topological transition point of view in boron nitride nanotubes upon compression. Due to the increased curvature in compressed small-diameter MWBNNTs, the uncommon 4- and 8-membered rings in Z-BN could be more favorably formed. And the irreversible tube collapse is proved to be a critical factor for the capture of the formed Z-BN, because of the competition between the resilience of tube before collapse and the stress limitation for the lattice stabilization of Z-BN upon decompression. In this case, Z-BN starts to form above 19.0 GPa, which is fully reversible below 45 GPa and finally becomes quenchable at 93.5 GPa. Consequently, this collapse-induced capture of the high-pressure phase could also be extended to other tubular materials for quenching novel sp 3 structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Record-breaking 2020 summer marine heatwaves in the western North Pacific

Record-breaking marine heatwaves (MHWs) occurred in the western North Pacific during the summer of 2020. These unprecedented MHWs were consistent with favorable large-scale conditions that are linked to an anomalous western North Pacific subtropical high (WNPSH), resulting mainly from sea surface temperature (SST) anomalies across the tropical oceans. In addition, a moderate La Niña-like pattern was also conducive to transporting warm seawater to the western North Pacific. Mixed-layer heat budgets suggest that surface heat flux contributed to the SST anomaly in the western subtropical Pacific. In contrast, oceanic heat advection dominanted in the South China Sea and the western equatorial Pacific. Numerical model experiments indicated that the tropical Indian Ocean SST anomalies were responsible for the enhanced WNPSH. The increased zonal SST gradient across the tropical Pacific also played an important role. Inter-ocean interactions can modulate climate variability through ocean-atmospheric coupling and deserve more attention when predicting MHWs within the context of global warming. Additionally, it is critical to consider MHWs as a powerful tool in detecting acute, intense thermal stress events in the coral bleaching pre-warning system.

Coral bleaching↗
DOE OSTI2021

Molecular packing modulation enabling optimized blend morphology and efficient all small molecule organic solar cells

The molecular planarity of a small molecule donor is one of the key factors affecting its molecular packing behavior, the morphology of the related active blend and thus the photovoltaic properties of all small molecule organic solar cells (ASM-OSCs). Here, a series of new benzodithiophene (BDT)-based A-π-D-π-A type small molecules were developed, namely SD1, SD2 and SD3, by combining the chemical modification of side chain (sulfuration), π-bridge (alkoxy groups introducing) and end group (different type or different alkyl chain). The sulfuration of the BDT side chains and introducing alkoxy group to the π-bridge could induce intermolecular S … S and intramolecular S … O non-covalent interactions for higher planarity of these small molecules. The rhodanine and cyanoacetate end groups showed less effect on the optical and electronic properties but significantly affected the molecular packing behavior and blend morphology. From SD1:Y6-T to SD2:Y6-T and SD3:Y6-T, the face-to-face packing was gradually reduced, while the face-to-edge packing was increased, which resulted in lowered exciton dissociation efficiency, further decreased the device performance. SD1:Y6-T blend exhibited more favorable face-to-face packing and ideal morphology for efficient charge dissociation and transport, contributing to a champion PCE of 10.12%. This work provided a deep insight of the donor and acceptor packing behavior effect on the photovoltaic property by molecular design.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Cross-scale modeling of storm surge, tide, and inundation in Mid-Atlantic Bight and New York City during Hurricane Sandy, 2012

Driven by high-resolution NAM (North America Mesoscale) and ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric fields, a 3-D unstructured grid SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model) was applied to simulate total water level in the Mid-Atlantic Bight during Hurricane Sandy in 2012. The simulated storm surge, tide, significant wave height, and wave peak period results compared favorably with NOAA observations along US East Coast, Long Island Sound, and New York Harbor. The maximum total water level at The Battery in the New York Harbor was accurately simulated with an absolute error of less than 0.08 m (out of 2.9 m peak surge) and a timing difference less than 10 min, which includes the effects on the order of 0.1–0.3 m (5–10%) from the coastal setup induced by the surface waves. The scenario comparisons of (1) “NAM” versus “ECMWF”, (2) “2-D” versus “3-D”, and (3) “with” versus “without” wind wave model were examined. The 3-D barotropic model forced by ECMWF including the effects of wind waves performs the best, attributed to wave-induced radiation stress and reduction of bottom stress when the 3-D version is used. Simultaneously, the ELCIRC-Sub model, using a 5-m topography/bathymetry sub-grid and a regular 200-m resolution finite volume computational model grid, was developed to simulate the street-level inundation in New York City. The momentum and mass fluxes calculated by the coarser base grid model were effectively coupled with the sub-grid so that running a full-blown high-resolution base model is not required. The ELCIRC-Sub model uses MPI (Message Passing Interface) parallel computing to enlarge the coverage of the land surface and an efficient non-linear solver to improve the accuracy of the wetting-and-drying processes. The temporal comparisons of modeled water level with NOAA's tidal stations and USGS' rapid-deployed gauges showed overall performance with an average error of 0.1 m. Particularly, ELCIRC-Sub captured the profile of the highest peak surge (3.9 m), which rose 2 m within 3 hour at Kings Point, NY. The spatial comparisons of the modeled peak water level at 80 surveyed locations showed an average error less than 0.13 m. The modeled maximum inundation extent also matched well with an 80% hit rate against FEMA's Hurricane Sandy maximum flooding extent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A multi-pronged evaluation of aldehyde-based tripeptidyl main protease inhibitors as SARS-CoV-2 antivirals

As an essential enzyme of SARS-CoV-2, the COVID-19 pathogen, main protease (M Pro ) is a viable target to develop antivirals for the treatment of COVID-19. By varying chemical compositions at both P2 and P3 positions and the N-terminal protection group, we synthesized 18 tripeptidyl M Pro inhibitors that contained also an aldehyde warhead and β-(S-2-oxopyrrolidin-3-yl)-alaninal at the P1 position. Systematic characterizations of these inhibitors were conducted, including their in vitro enzymatic inhibition potency, X-ray crystal structures of their complexes with M Pro , their inhibition of M Pro transiently expressed in 293T cells, and cellular toxicity and SARS-CoV-2 antiviral potency of selected inhibitors. These inhibitors have a large variation of determined in vitro enzymatic inhibition IC 50 values that range from 4.8 to 650 nM. Here, the determined in vitro enzymatic inhibition IC 50 values reveal that relatively small side chains at both P2 and P3 positions are favorable for achieving high in vitro M Pro inhibition potency, the P3 position is tolerable toward unnatural amino acids with two alkyl substituents on the α-carbon, and the inhibition potency is sensitive toward the N-terminal protection group. X-ray crystal structures of M Pro bound with 16 inhibitors were determined. In all structures, the M Pro active site cysteine interacts covalently with the aldehyde warhead of the bound inhibitor to form a hemithioacetal that takes an S configuration. For all inhibitors, election density around the N-terminal protection group is weak indicating possible flexible binding of this group to M Pro . In M Pro , large structural variations were observed on residues N142 and Q189. Unlike their high in vitro enzymatic inhibition potency, most inhibitors showed low potency to inhibit M Pro that was transiently expressed in 293T cells. Inhibitors that showed high potency to inhibit M Pro transiently expressed in 293T cells all contain O-tert-butyl-threonine at the P3 position. These inhibitors also exhibited relatively low cytotoxicity and high antiviral potency. Overall, our current and previous studies indicate that O-tert-butyl-threonine at the P3 site is a key component to achieve high cellular and antiviral potency for tripeptidyl aldehyde inhibitors of M Pro .

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Oxygen Vacancies Enriched Bi Based Catalysts for Enhancing Electrocatalytic CO2 Reduction to Formate

The practical application of CO2 conversion through electrochemical reduction is still hindered by the lack of high performance and cost-effective catalysts. It is thus imperative to explore feasible strategies for scalable fabrication of catalysts for CO2 reduction to targeting products with high efficiency. Herein, for the first time, a hybrid of BiOx encapsulated in carbon thin layer is developed as high active electrocatalyst for CO2 conversion, which is synthesized by pyrolysis using Bi MOF as precursor and sacrificial template. The carbon layer overcoating BiOx catalyst exhibits superior performance for CO2 reduction to formate, achieving a maximum Faradaic efficiency of 89.3 % and a maximum current density of 37.8 mA cm-2. It is found that the as-prepared BiOx catalyst is enriched with oxygen vacancies stabilized by the encapsulating carbon layer. Both the carbon layer and oxygen vacancies improve the charge transfer of the catalyst, in favor of the high formate current density. The defects of oxygen vacancies play a vital role in CO2 adsorption and activation and thereby further enhancing the electrocatalytic ability for formate production.

oxygen vacancy, bismuth, carbon coating, CO2 elect↗
DOE OSTI2021

CO 2 electro-reduction on Cu 3 P: Role of Cu(I) oxidation state and surface facet structure in C 1 -formate production and H 2 selectivity

We report the catalytic activity and mechanism of copper(I) phosphide, Cu 3 P, with predominant [00Ι] facet exposure for the electrochemical reduction of CO 2 (CO 2 RR) to formic acid. Crystalline nanosheets of this compound that show a preferential [00Ι] facet orientation exhibit undiminished CO 2 RR activity after 16 hours with full retention of crystal structure and surface chemical speciation and no detectable corrosion. In contrast to the range of products formed on Cu metal, CuO, and Cu 2 O, the CO 2 RR on Cu 3 P [00Ι] produces mainly hydrogen and formate as the sole carbon product in KHCO 3 electrolyte. Analysis of the Cu 3 P [00Ι] facet by HAADF-STEM was used to determine the surface lattice structure, while both XPS and Auger spectroscopies were used to determine the surface chemical speciation from the kinetic energies of ionized electrons. The presented analysis identifies isolated trigonal CuP 3 sites on the Cu 3 P[00Ι]-Cu 3 P 3 terminated surface and the Cu(I) oxidation state as precursor to the active catalyst. The CO 2 RR selectivity to formate and the higher turnover rate for H 2 production on the [00Ι] facet allows a structure-activity analysis and chemical mechanism to be proposed. Formation of a surface hydride at isolated *H-CuP 3 sites is proposed as the catalytic site in forming both H 2 and formate, while the long Cu—Cu separation retards forming C-C coupling products. In conclusion, these results disprove previously published claims of Cu(I) oxidation state as a sufficient criterion to promote CO 2 RR to C 2+ products, show that stronger bonded hydrides, *H-CuP 3 on Cu 3 P, favor production of the C1 product formate over all other carbon products, and predict that stronger formate binding (bidentate) is needed for CO 2 RR currents to compete with H 2 production.

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

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

Enhancing EnergyPlus capabilities to model dynamic building envelopes using python plugin

Nearly half of the energy consumption in the United States is related to buildings, resulting in an urgent need to develop innovative technologies to improve building energy efficiency. Dynamic building envelopes, comprising switchable insulation and thermal energy storage materials, have been proposed recently as a promising solution to reduce buildings' heating and cooling loads by thermally coupling the indoor environment with the ambient environment when beneficial while decoupling them when outdoor conditions are not favorable. Although various related technologies are still underway, the whole-building energy modeling tools, like EnergyPlus, do not have the capability to simulate the transient and dynamic nature of dynamic envelope materials and components to accurately predict their impact on building energy use. The objective of this study is to formulate a method in EnergyPlus simulation engine to model multilayer envelopes, comprising dynamic building materials with variable thermophysical properties, and discuss the changes made to the program using a Python plugin. Furthermore, the thermal performance of the dynamic envelopes using the proposed method is compared and verified with the results from a well-established commercial code, COMSOL Multiphysics. A parametric assessment is also conducted to evaluate the energy efficiency benefits of dynamic envelopes in a single-family residential building, demonstrating total annual energy savings up to 11.6 %, when a dynamic envelope operates alone, and up to 18.2 % when it is combined with a thin layer of phase change material as a thermal storage medium. Finally, a United States wide energy efficiency assessment is presented to showcase the geographical spread of the energy savings. The method designed and implemented in this study provides the researchers with the ability to implement their dynamic insulation methods in EnergyPlus and evaluate the whole building energy impact.

25 ENERGY STORAGE↗
DOE OSTI2025

An economic and technical feasibility analysis of a dual-source heat pump using both the air and the ground

The study investigates the economic and technical performance of a novel dual-source heat pump (DSHP) compared with that of air-source heat pumps (ASHPs) and ground-source heat pumps (GSHPs). The DSHP can use both ambient air and the ground as a heat source or heat sink. It uses ambient air when its temperature is favorable for efficient heat pump operation. When the ambient temperature is too hot or cold, the ground source is used to retain high-efficiency heat pump operation. Since the DSHP can alternately use either the ground heat exchanger (GHE) or ambient air to meet the thermal load, the required size of GHE can be smaller than those of GSHPs. This study models the DSHP using a whole building energy simulation tool (EnergyPlus) coupled with a Python plug-in and Heat Pump Design Model (HPDM) to simulate its heating and cooling performance for a typical single-family home in 15 US climate zones. The required GHE size of the DSHP system is determined through simulations and compared with that of GSHPs. DSHP deployment can reduce electricity use compared to ASHPs, especially in cold climates where it shows a reduction of around 50%. When compared to GSHPs, DSHPs use 20%–40% more electricity in warm climates but consume around the same amount in moderate and colder climates. Since the DSHP can use air source when the ambient temperature is mild, the GHE size needed for the DSHP is about 40% less than that needed for GSHPs in hot climates and about 25% less in cold climates. In conclusion, the life cycle cost analysis shows that the DSHP is economically more feasible than ASHPs in colder regions and economically more feasible than GSHPs in hot and cold regions.

Dual-source heat pumps↗
DOE OSTI2020

Enriched hydrogen production over air and air-steam fluidized bed gasification in a bubbling fluidized bed reactor with CaO: Effects of biomass and bed material catalyst

Gasification is one of the methods of generating biopower or biofuels from biomass waste. In this study, a benchscale fluidized bed reactor was used for biomass air and air-steam gasification. Gasification was performed under constant operating conditions (~780 °C, equivalence ratio = ~0.32) to investigate the effect of biomass (switchgrass, pine residues) and bed materials (sand, CaO+ sand, Al 2 O 3 , and CaO + Al 2 O 3 ). All gasification products, such as synthesis gas (syngas), contaminant gases, tar, and biochar (solid) were comprehensively analyzed. The composition of biomass significantly impacted CO and H 2 yield from volatile combustible matter and fixed carbon. Further, the presence of CaO made the condition favorable for the water-gas shift (WGS) reaction combined with the CO 2 carbonation reaction, which increased H 2 concentration. Additional steam with CaO increased H 2 concentration closer to 50% (N 2 free condition) through the combination reactions of steam hydrocarbon reforming and WGS by producing 44 g H 2 /kg dry biomass and 143 g CO /kg dry biomass . The usage of steam reduced the overall yield of contaminant gases, whereas the usage of CaO or Al 2 O 3 decreased the amount of gasification tar by approximately 5.8–6.5 g tar /kg dry biomass . In conclusion, this study can provide valuable experimental data for biomass waste to produce better quality syngas.

08 HYDROGEN↗
DOE OSTI2025

Optimal sizing of battery energy storage systems for peak shaving and demand response using a degradation-aware Bayesian Optimization-Mixed-Integer Linear Programming framework

The increasing integration of renewable energy and rising electricity demand highlight the importance of battery energy storage systems for peak shaving and demand response. Unlike prior approaches that overlook operational impacts on degradation, this study proposes a Bayesian Optimization–Mixed Integer Linear Programming framework for optimal battery energy storage system sizing. In this framework, Mixed Integer Linear Programming determines short-term scheduling while a calibrated electrochemical model iteratively evaluates degradation. The central hypothesis is that the framework can efficiently identify optimal sizes that yield realistic and economically robust outcomes. The method is tested across three scenarios: peak shaving, peak shaving with energy-reduction demand response, and peak shaving with power-reduction demand response. Results show that the framework converge to the optimum within 20 iterations out of 150 possible sizes. Under baseline conditions, the framework consistently selects the smallest feasible system, minimizing unnecessary degradation costs from oversized storage. Sensitivity analyses reveal that larger systems are favored as demand rates or incentives increase. Comparisons of demand response programs indicate that power-reduction demand response offers greater economic benefits than energy-reduction demand response, although demand savings from peak shaving remain the dominant contributor to overall performance. This study demonstrates that the proposed framework balances computational tractability with degradation fidelity, identifies critical economic thresholds for investment, and offers a practical, flexible tool to guide industrial stakeholders in cost-effective battery energy storage system deployment.

Batteries↗
DOE OSTI2021

Maximizing oil production from water alternating gas (CO 2 ) injection into residual oil zones: The impact of oil saturation and heterogeneity

Residual oil zones (ROZs) are widespread reservoirs, characterized by oil at residual saturation, either underlying oil fields (brownfield) or lateral (greenfield) to such fields. These reservoirs have the potential to produce volumes of oil sufficiently significant to make appreciable impacts on the US’s oil reserves and associated incidental CO 2 sequestration. The objective of this study is to improve our understanding the impact of heterogeneous and low oil saturations, in brownfield ROZs, on the effectiveness of water alternating gas (WAG) injection strategies. ROZs occur in the Permian Basin and elsewhere, and operators are using CO 2 injection for enhanced oil recovery (EOR) in these zones. The consensus model for the formation of ROZs is that they were formed by the effect of faster regional aquifer flow, acting over millions of years. Both the magnitude of oil saturation and the spatial distribution of oil differ from water-flooded main pay zones (MPZs). To explore the most effective injection strategies, we conducted simulations of CO 2 injection into synthetic geologic reservoirs. These simulations focused on injection into reservoirs subject to either man-made waterflooding or long-term natural waterflooding. By exploring the impact of varying: oil saturation; well patterns; reservoir heterogeneity; and permeability anisotropy, we attempt to quantify the factors that most influence the effectiveness of WAG injection. WAG ratios (the ratio of injected water and CO 2 , in reservoir volumes) of interest are those that either minimize the net CO 2 utilization ratios or maximize oil production rates. In general, the most effective WAG ratios for ROZs, are consistently less than those observed undergoing CO 2 injection in the same geologic reservoir models after traditional (man-made) waterflooding. This work demonstrates that most favorable WAG ratios for oil production in ROZs are different from those in traditional MPZs because of oil saturation differences. Therefore, CO 2 injection into both zones or directly copying WAG injection designs from MPZs to ROZs might not maximize oil production.

04 OIL SHALES AND TAR SANDS↗