Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “salinity gradient energy”

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.

34 records · Page 2

Tree-inspired lignin microrods-based composite heterogeneous nanochannels for ion transport and osmotic energy harvesting

One of the key processes of tree lignification is that lignin penetrates into the cell wall and fills in the cell wall framework, thereby increasing the hardness and hydrophobic of the tree channels, which is beneficial to consolidate and support the tree cell wall and water transport. Inspired by this natural process, we demonstrated a lignin-based nanofluidic heterogeneous membrane that closely mimics the channels in tree, which can realize ion transport function and effectively capture reverse electrodialysis. The membrane was synthesized by heating dealkaline lignin and PVA at 200 °C and this formed fusiform microrods and a closed-packed membrane. Simultaneously, this membrane composited with anodized aluminum (AAO) channels membrane at 200 °C to form asymmetric heterogeneous nanochannels membrane, which can transport counter-ions and harvest osmotic energy. This membrane implements ion current rectification in 0.1 M KCl electrolyte solution at pH 3 due to the confinement of pores and opposite surface charges in lignin-based heterogeneous nanochannel. An output power density of 0.97 W m- 2 is obtained under a 50-fold salinity gradient, which can be further improved to 1.19 W m -2 by increasing the salinity gradient from 50-fold to 500-fold. Accordingly, this nanofluidic membranes were prepared by using lignin, the key component in tree, which not only mimicked a crucial process of the water and ionic transport process of channels in tree, but also had the prospect in the field of osmotic energy harvesting.

09 BIOMASS FUELS↗

Efficiently resolving the terrestrial-aquatic interface in E3SM with sub-grid methods to improve coastal simulations (Final Technical Report)

The broad goal of this project was to extend the capabilities of MPAS‐Ocean, the ocean component of DOE’s Energy Exascale Earth System Model (E3SM), to better resolve the hydrodynamics of the terrestrial‐aquatic interface without significantly increasing computing resource requirements. Specifically, our goal was to better resolve the small-scale features across the terrestrial‐aquatic interface which includes the low lying coastal floodplain and a dendritic network of hydraulically efficient channels connected to estuarine systems and back bays, barrier islands and high energy inlets. The processes driving flow across the terrestrial‐aquatic interface include upland and floodplain hydrology, tides, wind and atmospheric pressure (including coastally generated storm surge across continental shelves and within estuaries and bays), wind waves and water temperature and salinity gradients.

58 GEOSCIENCES↗

Eco-evolutionary strategies for relieving carbon limitation under salt stress differ across microbial clades

With the continuous expansion of saline soils under climate change, understanding the eco-evolutionary tradeoff between the microbial mitigation of carbon limitation and the maintenance of functional traits in saline soils represents a significant knowledge gap in predicting future soil health and ecological function. Through shotgun metagenomic sequencing of coastal soils along a salinity gradient, we show contrasting eco-evolutionary directions of soil bacteria and archaea that manifest in changes to genome size and the functional potential of the soil microbiome. In salt environments with high carbon requirements, bacteria exhibit reduced genome sizes associated with a depletion of metabolic genes, while archaea display larger genomes and enrichment of salt-resistance, metabolic, and carbon-acquisition genes. This suggests that bacteria conserve energy through genome streamlining when facing salt stress, while archaea invest in carbon-acquisition pathways to broaden their resource usage. These findings suggest divergent directions in eco-evolutionary adaptations to soil saline stress amongst microbial clades and serve as a foundation for understanding the response of soil microbiomes to escalating climate change.

54 ENVIRONMENTAL SCIENCES↗

Performance of a Natural Gas Solid Oxide Fuel Cell System With and Without Carbon Capture

The fuel cell program at the United States Department of Energy (DOE) National Energy Technology Laboratory (NETL) is focused on the development of low-cost, highly efficient, and reliable fossil-fuel-based solid oxide fuel cell (SOFC) power systems that can generate environmentally-friendly electric power with at least 90 percent carbon capture. NETL’s SOFC technology development roadmap is aligned with near-term market opportunities in the distributed generation sector to validate and advance the technology while paving the way for utility-scale natural gas (NG)- and coal-derived synthesis gas-fueled applications via progressively larger system demonstrations. The present study represents a part of a series of system evaluations being carried out at NETL to aid in prioritizing technological advances along research pathways to the realization of utility-scale SOFC systems, a transformational goal of the fuel cell program. In particular, the system performance of utility-scale NG fuel cell (NGFC) systems with and without carbon dioxide (CO2) capture is presented. The NGFC system analyzed features an external auto-thermal reformer (ATR) feeding the fuel to the SOFC system consisting of planar anode-supported SOFC with separated anode and cathode off-gas streams. In systems with CO2 capture, an air separation unit (ASU) is used to provide the oxygen for the ATR and for the combustion of unutilized fuel in the SOFC anode exhaust along with a CO2 purification unit to provide a nearly pure CO2 stream suitable for transport for usage in enhanced oil recovery operations or for storage in underground saline formations. Remaining thermal energy in the exhaust gases is recovered in a bottoming steam Rankine cycle while supplying any process heat requirements. A reduced order model (ROM) developed at the Pacific Northwest National Laboratory (PNNL) is used to predict the SOFC performance. The ROM, while being computationally effective for system studies, provides other detailed information about the state of the stack, such as the internal temperature gradient, generally not available from simple performance models often used to represent the SOFC. Such additional information can be important in system optimization studies to preclude operation under off-design conditions that can adversely impact overall system reliability. The NGFC system performance was analyzed by varying salient system parameters, including the percent of internal (to the SOFC module) NG reformation—ranging from 0 to 100 percent—fuel utilization, and current density. The impact of advances in underlying SOFC technology on electrical performance was also explored.

solid oxide fuel cell (SOFC), natural gas fuel cel↗

Machine learning and deep learning for mineralogy interpretation and CO 2 saturation estimation in geological carbon Storage: A case study in the Illinois Basin

Carbon capture and storage (CCS) is a promising approach to simultaneously maintaining energy security and reducing carbon dioxide (CO 2 ) emissions under the current energy portfolio that is dominated by fossil fuel energy. Pre-injection formation characterization and post-injection CO 2 monitoring are two critical tasks to guarantee storage efficiency in CCS. The CCS projects in the Illinois Basin, the first large-scale CO 2 injection into saline aquifers in the United States, employed conventional and the latest pulsed neutron logging (PNL) tools for mineralogy interpretation and CO 2 saturation estimation, which provide valuable references for future CCS projects. Because of the inherent fuzziness of petrophysical measurements and complex subsurface heterogeneity, interpreting well-logging data is time-consuming, and its accuracy can be user-biased. In recent years, data-driven methods have been widely used to capture the non-linear patterns between input features and interpretation results. This work applied and evaluated four commonly used machine learning (ML) models, including ridge regression (RR), random forest (RF), gradient boosting regression (GBR), support vector regression (SVR), and one deep learning (DL) model, the artificial neural network (ANN). We optimized the hyperparameters of the four ML models and the DL model using the simulated annealing algorithm and the grid search strategy, respectively. The input features of the mineralogy interpretation models were eleven conventional well-logging parameters, and the label data (i.e., ground truth) were the porosity and volumetric fractions of six minerals, including quartz, feldspar, dolomite, calcite, clay, and iron minerals. The results demonstrated that the GBR and RF models were superior in predicting volumetric fractions of minerals and porosity; label data with low coefficient of variation (CV) values tended to yield better performance. For CO 2 saturation estimation, the RF was the best-performing model, followed by SVR, ANN, GBR, and RR. Furthermore, we conducted feature importance ranking using the permutation importance algorithm and found that the formation sigma and well pressure were the most important features in this study. In conclusion, the study of CCS projects in the Illinois Basin bridges the gap between the limited knowledge and understanding of geological carbon storage and the increasing demand for reliable, cost-effective, and sustainable energy solutions.

58 GEOSCIENCES↗

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DISCOVR: Development of Integrated Screening, Cultivar Optimization, and Verification Research

The DISCOVR consortium has developed and is applying an integrated platform and workflow for standardized, deep characterization of highly productive and resilient microalgae strains. New strains are screened using three consecutive tiers in a conceptual down-selection funnel to arrive at the most promising strains. At tier I, strains are screened on gradient incubators to determine their temperature and salinity tolerance range, and to evaluate their resistance to infectious agents (e.g., chytrids). At tier II, the strains’ winter and summer season biomass productivities are quantified in Laboratory Environmental Algae Pond Simulator (LEAPS) photobioreactors, and biomass composition is determined under N-replete and N-deplete conditions. At tier III, the best winter and summer season strains are evaluated in indoor crash ponds to determine their resilience to induced infections, and strains are tested in outdoor raceway ponds to confirm high biomass productivity and culture stability. The top strains are then evaluated in seasonal SOT outdoor pond culture trials at AzCATI, the dedicated Department of Energy SOT (State of Technology) testbed. Additional research is carried out on evaluating growth promoting molecules, reducing the impact of oxygen inhibition, optimizing pond operational conditions, measuring compositional shift dynamics, developing integrated pest management strategies, analysing culture health using 96 well-plate tools, implementing spectroradiometric monitoring and using machine learning to predict culture crashes.

Huesemann, Michael H↗

Viral niche-partitioning: comparative genomics of giant viruses across environmental gradients in a high Arctic freshwater-saltwater lake

Giant viruses (GVs; Nucleocytoviricota) impact the biology and ecology of a wide range of eukaryotic hosts, with implications for global biogeochemical cycles. Here, we investigated GV niche separation in highly stratified Lake A at the northern coast of Ellesmere Island, Nunavut, Canada. This lake is composed of a layer of ice-covered freshwater that overlies saltwater derived from the ancient Arctic Ocean, and it therefore provides a broad gradient of environmental conditions and ecological habitats, each with a distinct protist community and rich assemblages of associated GVs. The upper layer (mixolimnion) had measurable light and oxygen, and contained diverse GVs linked to photosynthetic protists, indicating adaptation to surface biotic and abiotic conditions. In contrast, the saline lower layer (monimolimnion), lacking oxygen and light, hosted GVs associated with predicted heterotrophic protists, some of which are known for a predatory lifestyle, and with several viral genes suggesting adaptation to deep-water anaerobic conditions. Our observations underscore the coupling between physical and chemical gradients, microeukaryotes and their associated GVs in Lake A, and provide insight into the potential for GVs to directly and indirectly impact host metabolism. There were similarities between the genetic composition of GVs and the metabolic processes of their potential hosts, implying co-evolution and niche-adaptation within the lake habitats. Notably, we found a greater presence of viral rhodopsins in deeper water layers, suggesting an evolutionary relationship with potential hosts capable of supplementing their energetic needs to thrive in low energy, anoxic conditions.

59 BASIC BIOLOGICAL SCIENCES↗

Sensitivities of the West Greenland Current to Greenland Ice Sheet Meltwater in a Mesoscale Ocean/Sea Ice Model

Abstract Meltwater from the Greenland Ice Sheet can alter the continental shelf/slope circulation and cross-shelf freshwater fluxes and limit deep convection in adjacent basins through surface freshening. We explore the impacts on the West Greenland Current and eastern Labrador Sea with different vertical distributions of the meltwater forcing. In this study, we present the results from global coupled ocean/sea ice simulations, forced with atmospheric reanalysis, that are mesoscale eddy-active (∼2–3-km horizontal spacing) and eddy-permitting (∼6–7-km horizontal spacing) in the study region. We compare the West Greenland Current in mesoscale eddy-active and eddy-permitting without meltwater to highlight the role of small-scale features. The mesoscale eddy-active configuration is then used to assess the change in the eastern Labrador Sea when meltwater is added to the surface or vertically distributed to account for mixing within fjords. In both simulations with meltwater, the West Greenland and West Greenland Coastal Currents are faster than in the simulation with no meltwater; their mean surface speeds are the highest in the vertical distribution case. In the latter case, there is enhanced baroclinic conversion at the shelf break compared to the simulation with no meltwater. When meltwater is vertically distributed, there is an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy. In addition, in the eastern Labrador Sea, the salinity is lower and the meltwater volume is greater when meltwater is vertically distributed. Therefore, the West Greenland Current is sensitive to how meltwater is added to the ocean with implications for the freshening of the Labrador Sea. Significance Statement Our goal is to understand how the flux of freshwater across the West Greenland continental slope into the Labrador Sea is modified by meltwater from the Greenland Ice Sheet. We compare the simulations of the ocean that capture key dynamics along the West Greenland continental slope that have no meltwater, meltwater added to the ocean surface, and meltwater distributed vertically to represent the mixing within fjords. When meltwater is added, the currents along the continental slope are faster, with the greatest increase when meltwater is vertically distributed. In that case, there is enhanced freshening of the Labrador Sea because modified density gradients generate more eddies. Proper representation of the vertical structure of meltwater is important for projecting the impact of freshwater on the subpolar North Atlantic.

58 GEOSCIENCES↗

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗

Changes in environmental and engineered conditions alter the plasma membrane lipidome of fractured shale bacteria

ABSTRACT Microorganisms that persist in fractured shale reservoirs cause several problems including secreting foul gases and forming biofilms. Current biocontrol measures often fail due to limited knowledge of their in situ activities. The plasma membrane protects the cell, mediates many of its critical functions, and responds to intracellular cues and ecological perturbations through physicochemical modifications. As such, it provides valuable insight into the physiological adaptation of microorganisms in disturbed environmental systems. Here, we (i) demonstrate how changes in salinity and hydraulic retention time (HRT) influence the plasma membrane intact polar lipid (IPL) chemistry of model bacterium, Halanaerobium congolense WG10, and mixed microbial consortia enriched from shale-produced fluids and (ii) elucidate adjustments in membrane IPL chemistry during biofilm growth relative to planktonic cells. We incubated H. congolense WG10 in chemostats under three salinities (7%, 13%, and 20% NaCl), operated under three HRTs (19.2, 24, and 48 h), and in drip flow biofilm reactors under the same salinity gradients. Also, mixed microbial consortia in produced fluids were enriched in triplicate chemostat vessels under three HRTs (19.2, 24, and 72 h) and biofilm reactors. Lipids were analyzed by ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Our results show that phosphatidylglycerols, cardiolipins, and phosphatidylethanolamines were predominantly enriched in planktonic H. congolense WG10 cells grown at hypersalinity (20%) compared to optimum (13%). In addition, several zwitterionic phosphatidylcholines and phosphatidylethanolamines were higher in abundance during biofilm growth. These observations suggest that microbial adaptation and biofilm formation in fractured shale are enabled by strategic plasma membrane IPL chemistry adjustments. IMPORTANCE Microorganisms inadvertently introduced into the shale reservoir during fracturing face multiple stressors including brine-level salinities and starvation. However, some anaerobic halotolerant bacteria adapt and persist for long periods of time. They produce hydrogen sulfide, which sours the reservoir and corrodes engineering infrastructure. In addition, they form biofilms on rock matrices, which decrease shale permeability and clog fracture networks. These reduce well productivity and increase extraction costs. Under stress, microbes remodel their plasma membrane to optimize its roles in protection and mediating cellular processes such as signaling, transport, and energy metabolism. Hence, by observing changes in the membrane lipidome of model shale bacteria, Halanaerobium congolense WG10, and mixed consortia enriched from produced fluids under varying subsurface conditions and growth modes, we provide insight that advances our knowledge of the fractured shale biosystem. We also offer data-driven recommendations for improving biocontrol efficacy and the efficiency of energy recovery from unconventional formations.

03 NATURAL GAS↗

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗

A hierarchy of global ocean models coupled to CESM1

We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed-layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three-dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed-layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid-latitude westerly jet and of the Inter-tropical Convergence Zone (ITCZ) in response to sea-ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed-layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.

54 ENVIRONMENTAL SCIENCES↗

DISCOVR strain pipeline screening – Part I: Maximum specific growth rate as a function of temperature and salinity for 38 candidate microalgae for biofuels production

Here, to identify high productivity strains for microalgal biofuels generation, the maximum specific growth rate of 38 strains was measured as a function of salinity (i.e., 5, 15, and 35 PSU) and temperature (i.e., at 8 temperatures along a linear gradient from ca. 5 to 45°C) to determine the most suitable growth medium salinity and best growing season, respectively, for outdoor raceway pond cultivation. The following strains were evaluated: Agmenellum quadruplicatum UTEX 2268, Anabaena sp. ATCC 33081, Arthrospira fusiformis UTEX 2721, Arthrospira platensis UTEX 3086, Chlorella vulgaris NREL 4-C12, Chlorella autotrophica CCMP 243, Chlorella sorokiniana DOE1044, Chlorella sorokiniana DOE 1116, Chlorella sorokiniana DOE 1412 (UTEXB3016), Chlorella vulgaris LRB AZ-1201, Chlorococcum littorale UTEX 117, Chlorococcum sp. UTEX-B P7, Chloromonas reticulata CCALA 870, Coelastrella sp. DOE 0202, Cyanobacterium sp. AB1, Micractinium reisseri NREL 14-F2, Microchloropsis gaditana CCMP1894, Microchloropsis salina CCMP 1776, Monoraphidium sp. MONOR1, Monoraphidium minutum 26B-AM, Nannochloropsis oceanica CCAP 849/10, Oscillatoria cf. priestleyi CCMEE 5020.1-1, Picochlorum celeri TG2-WT-CSM/EMRE, Picochlorum oklahomensis CCMP 2329, Picochlorum renovo NREL 39-A8, Picochlorum soloecismus DOE 101, Porphyridium cruentum CCMP 675, Scenedesmus acutus LRB-AP-0401, Scenedesmus obliquus DOE 0152.z, Scenedesmus obliquus UTEX393,Scenedesmus rubescens NREL 46B-D3, Scenedesmus sp. IITRIND2, Stichococcus minor CCMP 819, Stichococcus minutus CCALA 727, Synechococcus elongatus UTEX2973.1, Tetraselmis striata LANL 1001, Tisochrysis lutea CCMP 1324, and Tribonema minus UTEXB3156. For each strain, the identity and the presence of bacterial cohorts was determined using 18S and 16S rDNA sequencing, respectively. The maximum specific growth rate versus temperature data were also used to determine the activation energies (Arrhenius equation) for most strains. For all strains, the measured salinity and temperature tolerance data were compared to those reported in the literature. The fastest growing strains were down-selected for subsequent biomass productivity measurements in climate-simulation photobioreactors, as reported in the next paper in the issue.

18S and 16S rDNA sequencing↗

Managing Subsurface Pressure Buildup and Interference in Commercial-Scale CO 2 Storage Project with Proximal Injection Wells

Large-scale decarbonization using carbon capture and storage (CCS) is likely to involve many commercial-scale CO 2 storage projects located in close proximity to each other. This close proximity raises concerns over pressure interference among the storage projects. Pressure interference between injection and storage efforts can reduce the practicable CO 2 storage resource and force wells to inject CO 2 at a lower rate to avoid the fracture pressure thresholds per United States Environmental Protection Agency (EPA) Class VI well regulations to preserve injection and confining zone integrity and potentially mitigate against inducing seismic activity. These analyses employ numerical full-physics reservoir modeling to evaluate how pressure buildup fronts and CO 2 plumes evolve under commercial-scale injection volumes of CO 2 in which multiple storage sites located in close proximity occur in tandem. The simulation models mimic injection at pseudo basin-scale and assume homogeneous saline formation(s) as storage targets with a pair of upper and lower homogeneous seal layer/s. These analyses specifically investigate the efficacy of two basin-wide reservoir pressure management strategies in addressing the technical challenges associated with pressure buildup and CO 2 plume commingling. The strategies explored include: 1) enlarging the area of injection well spacing (WS) and 2) storing CO 2 in a stacked sequence (SSS) of saline formations compared to a single formation. The storage and confining zones properties assumed were held common across the scenarios, unless specified otherwise. Analyses results show that after injecting 4 million tons per year for 30 years using 4 separate wells (each injecting 1 million metric tons per year), the radius of CO 2 plume extends to a mere 3 km or less from injection wells. Meanwhile, the radius of pressure buildup ranges on the order of tens to a few hundreds of kilometers, depending on the magnitude of pressure buildup threshold that one would use to define the front. CO 2 plume commingling from different injection wells appears to occur 50 years post-injection, especially under scenarios with narrowly spaced (i.e., < 5 km apart) injection well locations. Findings from sensitivity cases on the well spacing suggest that storage formations modeled would require different well spacing to avoid fracture pressure thresholds. For instance, modeled storage formations with high fracture gradients (i.e., 0.8 psi/ft) would need less than 5–km well spacing, whereas those with lower fracture gradients (i.e., 0.7 psi/ft) would need approximately 20–km well spacing, based on assumed modeling parameters. Under stacked injection, the pressure challenges (described above) still exist but are more alleviated due to distributing the same injection volume across more available reservoir volume. These analyses demonstrate that stacked-sequence storage can effectively address the challenges, while still providing the same target CO 2 storage volumes and allowing a large number of storage projects to be deployed in the same basin by better utilizing the available storage resource across different reservoir depths. Among cases modeled, the resulting pressure buildup front is most suppressed when each storage project distributes injection volumes over several wells, each of which injects a portion of the total CO 2 across the stacked sequence. This strategy results in the smallest CO 2 aerial footprint amongst scenarios evaluated but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 42 percent compared to the commercial-scale single-formation storage), the result of which is crucial to maintain caprock integrity. The findings presented by this research draw attention to the importance of greater coordination among storage operators and regulatory stakeholders to foster the upscaling and deployment of CCS. These analyses provide insights into required decision-making when considering multi-project deployment in a shared basin. Because these analyses evaluate a very specific geologic situation, they bear further investigations across other geologic situations.

42 ENGINEERING↗