Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational analysis”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An Experimental Comparison of Cyclic Variations in Diesel–Natural Gas and POMDME–Natural Gas Dual Fuel Combustion

Abstract Cyclic variations in internal combustion engines are caused by various factors, including combustion mixture stratification, in-cylinder flows, local fluctuations in air-fuel ratio, etc. Cyclic variations have a profound impact on engine performance and emissions. In this study, cyclic variations in dual fuel combustion are analyzed, comparing diesel–natural gas (NG) and polyoxymethylene dimethyl ether (POMDME)-NG dual fuel combustion. Cyclic variability was initially quantified using the coefficient of variation of gross indicated mean effective pressure (IMEPg) computed from experimental cylinder pressure data. The cases analyzed in this study had a coefficient of variation (COV) of IMEPg greater than or around 5%, which was the lower limit of onset of instability for this engine. Experiments were performed at two fixed start of injection (SOI) of high-cetane fuel: 310 CAD and 350 CAD. For all experiments, a constant load of 5 bar IMEPg was maintained, and the intake boost pressure and rail pressure were fixed at 1.5 bar and 500 bar, respectively. For each case, 1000 cycles of cylinder pressure data were recorded, filtered, and processed using an in-house heat release analysis code for each cycle. A comparison between individual cycles and the “ensemble averaged cycle” was made for both diesel–NG and POMDME–NG combustion. For the early SOI of 310 CAD, the peak cylinder pressure fluctuations of individual cycle were found to be ± 15 bar for both fuel combinations, compared to the ensemble averaged cycle, and < 1/10th of the cycles had an IMEPg lower than 0.05 bar of the ensemble averaged cycle. However, the peak pressure fluctuations were found to be lower for POMDME–NG (±3 bar) than diesel–NG dual fuel combustion at 350 CAD SOI, indicating lower cyclic variations. The higher reactivity of POMDME helped reduce fluctuations in combustion phasing at the retarded SOI. The presence of cycles of deterioration and cycles of recovery were also observed with diesel–NG combustion for 310 CAD SOI, and the scatter in the IMEPg return map was similar for both fuel combinations. The IMEPg return map for POMDME–NG combustion was less scattered at the 350 CAD SOI.

Engineering↗
DOE OSTI2022

Heteroepitaxy of FCC-on-FCC systems of large misfit

Here, to understand the effects of lattice mismatch on heteroepitaxial growth, we have studied the equilibrium structure and orientation relationships (ORs) of FCC films grown epitaxially on FCC substrates, using molecular dynamics simulations in conjunction with embedded atom method potentials. Three film/substrate systems have been investigated, namely: Ag on Cu, Ag on Ni and Pb on Al. These systems cover a significant range of lattice mismatch, from 12.6% for Ag/Cu to 21.8% for Pb/Al. For each system, the ORs of films on six different substrate orientations, namely: (100), (511), (311), (211), (322) and (111), have been investigated. Films on these susbstrates cover a gradual transition from the oct-cube orientation relationship, which occurs only on (100) substrates, to the heterotwin orientation relationship, which often occurs on (111) substrates. It is found that the resulting ORs vary systematically with substrate orientation, but that the pattern of variation is almost identical for all three systems, and therefore largely independent of mismatch. However, the manner in which mismatch is accommodated does depend on the magnitude of mismatch. Simulations point to an important role for edge-to-edge matching and defects such as stacking faults. An analysis of these results in terms of transformation strains highlights the distinction between the ORs, which are largely independent of mismatch, and the local interfacial structure, which changes directly with mismatch.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Numerical Simulations of High Prandtl Number Liquid Jets Impinging on a Flat Plate

In this work, 3D simulations of oil jets impinging on a flat, heated wall are presented. The numerical setup uses the Volume of Fluid (VoF) method to model the two-phase flow. A careful grid definition across the liquid film, along with the use of the Conjugate Heat Transfer (CHT) approach allowed local heat transfer to be solved with fine resolution at the wall. Variations of liquid flow rate, liquid temperature and surface temperature allow to cover a wide range of local Reynolds and Prandtl numbers (226 < Re < 2850, 77 < Pr < 161). Resulting surface-averaged heat transfer compares very well with experimental measurements conducted in a previous study. In-depth analysis of the flow has identified expected features from the literature. In particular, the impact of jet axial velocity profiles on the heat transfer distribution in the stagnation zone was clearly stated. The increase in heat transfer when warming the liquid film was also reproduced and explained by a decrease in oil viscosity and an increase in film velocity. All those effects were taken into account in correlations for stagnation and local values of Nusselt number. A grid sensitivity study was also conducted, showing that if the grid solving the thermal boundary layer in the stagnation zone can be coarsened without impacting local and surface-averaged predictions of heat transfer, a minimum resolution (2 to 3 cells) within the thermal boundary layer is however required for an accurate prediction of heat transfer.

conjugate heat transfer↗
DOE OSTI2026

Composition, Activity, and Stability of IrO x Oxygen Evolution Reaction Electrocatalysts

The oxygen evolution reaction (OER) is integral to several electrochemical energy conversion and storage technologies, including carbon dioxide reduction to value added fuels, nitrogen reduction to ammonia, reversible fuel cells, rechargeable metal−air batteries, and water electrolysis to produce hydrogen. Iridium oxide (IrO x ) is widely recognized as the benchmark OER catalyst for acidic environments. Despite widespread use of IrO x catalysts, most notably in proton-exchange membrane water electrolyzers (PEMWEs), a comprehensive understanding of the physicochemical properties of commercial catalysts and the impact of these properties on both the activity and stability of these catalysts is lacking. Here, we study commercial IrO x catalysts with different physicochemical properties, three nominally considered amorphous and three rutile, to elucidate how structural and compositional variations affect OER activity and stability. Utilizing standardized aqueous electrochemical protocols, time-resolved dissolution quantification using inductively-coupled plasma mass spectrometry, and physicochemical characterization, including multiple synchrotron X-ray techniques, we systematically correlate catalyst properties with OER performance and degradation behavior aided by principal component analysis (PCA). Our results demonstrate the general trend of amorphous IrO x having higher intrinsic activity but limited stability and crystalline rutile IrO 2 having lower activity but enhanced stability against dissolution. The trends within the amorphous and rutile catalyst groups correlate with inherent material properties, including phase composition and structure, crystallinity, particle size, surface area, and surface structure/chemistry. Notably, we identify a rutile catalyst with the largest crystallite/ domain sizes, moderate surface area, a small fraction of hydrous phase, and a favorable pore structure (trimodal distributions of pore sizes ranging from 2−5 nm) that exhibits the best balance between activity and stability among the six catalysts studied here. These findings illustrate a fundamental structure-governed trade-off between activity and stability and highlight the critical role of surface chemistry modification and structure engineering in IrO x catalyst optimization.

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

A Bayesian Deep Learning Approach to Near-Term Climate Prediction

Since model bias and associated initialization shock are serious shortcomings that reduce prediction skills in state-of-the-art decadal climate prediction efforts, we pursue a complementary machine-learning-based approach to climate prediction. The example problem setting we consider consists of predicting natural variability of the North Atlantic sea surface temperature on the interannual timescale in the pre-industrial control simulation of the Community Earth System Model. While previous works have considered the use of recurrent networks such as convolutional LSTMs and reservoir computing networks in this and other similar problem settings, we currently focus on the use of feedforward convolutional networks. In particular, we find that a feedforward convolutional network with a Densenet architecture is able to outperform a convolutional LSTM in terms of predictive skill. Next, we go on to consider a probabilistic formulation of the same network based on Stein variational gradient descent and find that in addition to providing useful measures of predictive uncertainty, the probabilistic (Bayesian) version improves on its deterministic counterpart in terms of predictive skill. Finally, we characterize the reliability of the ensemble of machine learning models obtained in the probabilistic setting by using analysis tools developed in the context of ensemble numerical weather prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Arctic extreme precipitation changes from 1980 to 2022 in response to sea ice decline and enhanced atmospheric rivers

Arctic extreme precipitation (EP) broadly impacts permafrost degradation, glacier and snow cover changes, and ice sheet mass balance as well as ecosystems. However, investigation of EP spatiotemporal variations over the Arctic remains challenging, and their primary drivers are still poorly understood. Performance estimation of three state-of-the-art reanalysis products (Climate Forecast System Reanalysis (CFSR), European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA-5), and Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2)) against gauge-based precipitation observations reveals that MERRA-2 outperforms other reanalysis for Arctic EP changes. Based on MERRA-2 data, both annual EP amount and occurrence days averaged over the Arctic show statistically significant positive trends during 1980–2022 (3.37 ± 1.03 mm dec −1 and 0.42 ± 0.17 d dec −1 , respectively), with the most pronounced increase in the autumn. Spatial heterogeneity in annual and seasonal EP trends is found across the Arctic, with the largest positive annual trends of 30 mm dec −1 over the Bering Sea and the Denmark Strait. The significant EP increase is closely associated with intensified atmospheric rivers (ARs) and widespread decline in sea ice concentration (SIC). Specially, SIC and ARs are responsible for 12% and 50% of Arctic EP inter-annual variance, respectively, while ARs directly contribute 28.3% of the total annual EP amounts. These findings explain the mechanistic controls on Arctic EP, providing critical insights for projecting polar weather and climate extremes and their impacts on the Arctic environment.

atmospheric rivers↗
DOE OSTI2020

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Verification of a Modeling Toolkit for the Design of Building Electrical Distribution Systems

DC electrical distribution systems offer many potential advantages over their AC counterparts. They can facilitate easier integration with distributed energy resources, improve system energy efficiency by eliminating AC/DC converters at end-use devices (e.g., laptop chargers), and reduce installation material, time, and cost. However, DC electrical distribution systems present additional design considerations, largely resulting from potentially greater magnitude and variation in cable losses. Modeling and simulation are rarely used to design such systems. However, the greater dependency of DC system energy efficiency on design choices such as distribution voltages, architecture, and integration of PV and BESS suggests that modeling and simulation may be required. Such system performance analysis is currently not a standard practice, in part due to limited availability and validation of capable software tools. This paper characterizes the accuracy of a Modelica-based Building Electrical Efficiency Analysis Model (BEEAM) toolkit, as a precursor for validating its use to perform system performance analysis and inform design decisions. The study builds upon previous verification research by characterizing complete systems comprised of commercially available equipment, and providing a more detailed analysis of simulation results. Five lighting systems with varying electrical distribution architectures were designed using market-available equipment, installed in a laboratory environment, modeled using BEEAM, and simulated using three Modelica integrated development environments (IDEs). Simulated and measured results were compared to characterize toolkit accuracy. Initial results revealed that simulated performance was mostly within ±5% of measured system-level and device-level performance. While simulation results were not found to be dependent on the IDE, some Modelica compiler interoperability issues were identified. Although the BEEAM toolkit showed promise for the targeted use case, further work is needed to determine whether the demonstrated 5% accuracy is sufficient for making real-world design decisions, and for BEEAM to advance from an interesting research tool to one that can impact real-world building projects.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Grid impact analysis using controller-hardware-in-the-loop for high-power vehicle charging stations

A controller-hardware-in-the-loop (CHIL) architecture for the evaluation of grid impacts arising from high-power vehicle charging stations is presented in this paper. Unlike simulation-based studies, the proposed method can be used to capture the interactions of the grid and the charging load along with charger controllers in real time. The proposed method can be used to evaluate the impact of charging load on the grid in terms of voltage variations and line congestion. The proposed CHIL platform allows for de-risking the vehicle charging station deployment by simulating the complex interactions among all the components of a vehicle charging station - i.e., the grid, vehicle, and charger controller - in a realistic manner before using the charging station in a grid. Further, the proposed CHIL approach can be used to evaluate the voltage regulation causalities of the vehicle charging station. Experimental results are presented in the paper to illustrate the applicability of the proposed method in a laboratory environment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Reliability Assessment of Cooling Fans for PV Inverters: Testing, Modeling, and Case Studies

The reliability of photovoltaic (PV) inverters is critical for long-term solar system performance, with cooling fan failures frequently leading to costly downtime. While much research exists on general cooling fan reliability, little attention has been given to fans operating within PV inverters and their unique environmental challenges. Here, this article proposes a comprehensive methodology to address this gap. First, a failure mode and effects analysis is performed on fans to identify the key failure mechanisms in PV applications, their corresponding stressors, and the models necessary for lifetime prediction. Second, an accelerated life test is designed and conducted to collect valuable experimental data for PV inverter fans in a reasonable amount of time. Third, a mathematical conversion of dynamic mission profiles into effective constant stress levels is derived. Fourth, case studies are given, showcasing lifetime estimates that account for geographic variations in mission profile data. The results demonstrate that this integrated approach leads to an accurate reliability assessment for PV inverter cooling fans.

accelerated life testing (ALT)↗
DOE OSTI2021

Deep Analysis Net with Causal Embedding for Coal-fired Power Plant Fault Detection and Diagnosis (DANCE4CFDD)

Fault detection and diagnosis is critical to power plant operation to ensure attaining high reliability while reducing operation cost. As more renewable power is introduced to the power grid, traditional fossil power plants take on the extra burden of excessive load cycling to compensate the generation variability from renewable power. Such load cycling will pose more reliability challenges to power plant operation. There are a number of challenges faced by today’s asset health management system in coal- fired (or gas) power plants: 1) high-dimensional nonlinear interaction among multiple time series measurements; 2) high measurement variance induced by operational conditions/modes; 3) variation among asset types and plant configurations; and 4) a small number of faulty events to learn from. To cope with these challenges, today’s fielded asset health management systems rely heavily on manual efforts from domain experts and hand-crafted features or rules based on domain knowledge. Despite its role in plant reliability, such a practice is costly and hinders its scalability and sustainability, particularly when a plant undergoes modifications. The objective of this project is to develop a novel end-to-end AI learning system that is trainable (i.e., the AI representation of a complex system behavior can be directly learned from properly labeled data) for accurate fault detection and root cause analysis. The ability to create a fault detection model directly from time series could alleviate the efforts associated with today’s asset management solution development. In the course of this project, we have achieved the following: Created an AI model development environment incorporating state-of-the-art neural network architectures for rapid model development and evaluation; Developed novel learning strategies for training of fault detection model; Developed special-purpose neural network architecture embedded with variable association graph aiming for better interpretability; Developed a learning strategy to leverage a small number of faulty events for enhanced fault detection capability; Conducted detailed experimental study based on public benchmark datasets and demonstrated the effectiveness of the proposed solution; and Validated the developed system with data from both a coal-fired plant boiler dynamic simulation model and real-world coal-fired power plant covering multiple asset and fault types. Overall, the project attained a technology readiness level of TRL 5 from TRL 2 at the beginning of the project.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Optimizing the Synthetic Potential of O 2 : Implications of Overpotential in Homogeneous Aerobic Oxidation Catalysis

Molecular oxygen is the quintessential oxidant for organic chemical synthesis, but many challenges continue to limit its utility and breadth of applications. Extensive historical research has focused on overcoming kinetic challenges presented by the ground-state-triplet electronic structure of O 2 and the various reactivity and selectivity challenges associated with reactive oxygen species derived from O 2 reduction. This Perspective will analyze thermodynamic principles underlying catalytic aerobic oxidation reactions, borrowing concepts from the study of the oxygen reduction reaction (ORR) in fuel cells. This analysis is especially important for "oxidase"-type liquid-phase catalytic aerobic oxidation reactions, which proceed by a mechanism that couples two sequential redox half-reactions: (1) substrate oxidation, and (2) oxygen reduction, typically affording H 2 O 2 or H 2 O. The catalysts for these reactions feature redox potentials that lie between the potentials associated with the substrate oxidation and oxygen reduction reactions, and changes in the catalyst potential lead to variations in effective overpotentials for the two half reactions. Catalysts that operate at low ORR overpotential retain more thermodynamic driving force for the substrate oxidation step, enabling O 2 to be used in more challenging oxidations. While catalysts that operate at high ORR overpotential have less driving force available for substrate oxidation, they often exhibit different or improved chemoselectivity relative to the high-potential catalysts. The concepts are elaborated in a series of case studies to highlight their implications for chemical synthesis. Examples include comparisons of (a) NO x /oxoammonium and Cu/nitroxyl catalysts, (b) high-potential quinones and amine oxidase biomimetic quinones, and (c) Pd aerobic oxidation catalysts, with or without NO x cocatalysts. In addition, we show how reductive activation of O 2 provides a means to access potentials not accessible with conventional oxidase-type mechanisms. Altogether, this analysis highlights the central role of catalyst overpotential in guiding the development of aerobic oxidation reactions.

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

Economic analysis of the benefits to petroleum refiners for low carbon boosted spark ignition biofuels

A refinery modeling framework is developed to estimate the benefits of blending high-quality biofuels directly with refinery gasoline components for attaining a premium grade fuel (also termed as Co-Optima Boosted SI gasoline here). Our results change the paradigm of bio-blendstocks (BBs) being competitors to fossil components, by identifying opportunities for refineries to add value to their product slate, from some favorable BB properties. This potential value can be characterized by calculating the breakeven value (BEV), as defined down below. The proposed modeling framework incorporates extensive data from (1) projected product over the next few decades, (2) crude oil and refinery products pricing, and (3) fuel specifications. The complete refinery models serve as a basis for assessing the value of biofuels, assuming profitability remains the same for representative petroleum refinery configurations. Our assessment showed wide range of variation of biofuels BEV from $\$$20-$\$$120/bbl, within the considered blending level and crude prices. Further, the BEV was correlated with the fuel octane ratings such as octane numbers (research, RON and motor octane numbers, MON) and both antiknock index (AKI, average of RON and MON) and sensitivity (S, difference between RON and MON), with a slightly higher correlation with the sensitivity. However, the expected decrease in gasoline demand for the upcoming years could negatively impact biofuels demand and value, in a business-as-usual scenario. Our analysis also showed a more valuable bio-blendstocks incorporation in smaller refineries since they can enhance the capabilities for producing specialty, high-value fuels/products, and introduce high octane-barrels into otherwise constrained blending operations. Additional implications towards refiners include opportunities to rebalance operations, access to high-value fuel markets, and synchronization with broader transportation industry trends. Furthermore, results indicate the value of Co-Optima boosted spark ignition (BSI) efficiency gains can extend to refiners to incentivize decarbonization and diversified feedstock production.

09 BIOMASS FUELS↗
DOE OSTI2026

An International Laboratory Comparison Study on Approximating the Enthalpy of Adsorption via the Clausius‐Clapeyron Approach

Materials-based gas capture and storage is an increasingly important area of research. Robust and accurate determination of material properties is required for judicial selection of materials for specific applications and for engineering materials–based systems at scale. One key property is the strength of the adsorbate–adsorbent interaction often quantified via the isosteric enthalpy of adsorption. The heat of adsorption can be measured directly through calorimetry; however, a more widely used approach is to apply the Clausius-Clapeyron (CC) equation to adsorption isotherms collected at different temperatures. While this approach appears to be straightforward, there exist multiple variants in the application of the methodologies employed. This raises the question on how these variations may or may not affect the determined results. Presented here is a discussion of the most common methodologies and a comparison of indirect determinations (via CC) of the isosteric enthalpy of adsorption by different laboratories on identical material. Included in that comparison are discussions on the measurement and analysis reproducibility. Importantly, details of the methodologies are shown to be critical when comparing enthalpies among laboratories, and different methodologies contribute to significant discrepancies and artifacts in the results. Recommendations are provided to promote robust determination and the reporting thereof.

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

Monitoring Seasonal Shear Wave Velocity Changes in the Top 6 m at Garner Valley in Southern California With Borehole Data

Abstract Subsurface structures play important roles in seismic ground motion, crustal hydrology, stability of the built environment, and more. Constraining temporal changes of subsurface shear wave velocity ( V S ) can provide useful information to all these topics and the growing field of hydrological monitoring with seismic velocity. Using borehole records at Garner Valley, CA, we estimate seasonal subsurface V S variations from impulse response functions (IRFs) of earthquake data (2005–2018) along with IRFs and cross‐correlation of cross‐hole experiment data (2015–2018). The inferred V S variations are up to ∼25% in the top 6 m and ∼10% at 2–5 m in depth. The V S variations correlate strongly with the water table depth changes, suggesting that the changes are mostly due to fluctuations of pore pressure in the shallow material. The shallow velocity changes alter the near‐surface conditions, can affect seismic hazard estimation, and may be improperly attributed to deeper processes without careful analysis.

58 GEOSCIENCES↗
DOE OSTI2024

Random projection using random quantum circuits

The random sampling task performed by Google's Sycamore processor gave us a glimpse of the “quantum supremacy era.” This has definitely shed some light on the power of random quantum circuits in this abstract task of sampling outputs from the (pseudo)random circuits. In this paper, we explore a practical near-term use of local random quantum circuits in dimensional reduction of large low-rank data sets. We make use of the well-studied dimensionality reduction technique called the random projection method. This method has been extensively used in various applications such as image processing, logistic regression, entropy computation of low-rank matrices, etc. We prove that the matrix representations of local random quantum circuits with sufficiently shorter depths [ ∼ O ( n ) ] serve as good candidates for random projection. We demonstrate numerically that their projection abilities are not far off from the computationally expensive classical principal components analysis on MNIST and CIFAR-100 image datasets. We also benchmark the performance of quantum random projection against the commonly used classical random projection in the tasks of dimensionality reduction of image data sets and computing von Neumann entropies of large low-rank density matrices. And finally, using variational quantum singular value decomposition, we demonstrate a near-term implementation of extracting the singular vectors with dominant singular values after quantum random projecting a large low-rank matrix to lower dimensions. All such numerical experiments unequivocally demonstrate the ability of local random circuits to randomize a large Hilbert space at sufficiently shorter depths with robust retention of properties of large data sets in reduced dimensions. Published by the American Physical Society 2024

Kumaran, Keerthi (ORCID:0009000949125721)↗
DOE OSTI2025

Cryptic cycling by electroactive bacterioplankton in Trout Bog Lake

The potential for extracellular electron transfer (EET) is a prevailing genomic feature of humic lake bacterioplankton. However, there has been little evidence for the substantial ecological contribution predicted by genetics. We hypothesized that anoxygenic phototrophic electrotrophs and accompanying heterotrophic electrogens cycle dissolved organic matter (DOM) between oxidized and reduced states. We predicted that such bacterioplankton would exhibit diel-scale oscillations due to the light dependency of photosynthesis. Using Trout Bog Lake in Wisconsin, USA, as our model ecosystem, we profiled the water column with depth-discrete metagenomic, physiochemical, and electrochemical analyses. We observed variation in oxidation reduction potential (ORP) in response to sunlight, initiating at depths populated by anoxygenic phototrophs with EET genes. We developed an automated buoy to measure electric current flow between many pairs of electrodes simultaneously, observing correlation in electron consumption to sunlight. Our results, combined with published metatranscriptomic analysis, indicate the occurrence of electron cycling between phototrophic oxidation (electrotrophic metabolism) by Chlorobium and anaerobic respiration (electrogenic metabolism) by Geothrix, involving DOM. We also repeatedly observed gradual seasonal increases in hypolimnion ORP throughout summer. These diel and seasonal patterns imply that electroactive DOM mediates the ecology of electroactive bacteria in lakes, controlling humic lake methane emissions.IMPORTANCEWe investigated the physical, chemical, and redox characteristics of a bog lake and electrodes hung therein to test the hypothesis that dissolved organic matter is being cycled between oxidized and reduced states by electroactive bacterioplankton powered by phototrophy. To do so, we performed field-based analyses on multiple timescales using both established and novel instrumentation. We paired these analyses with recently developed bioinformatics pipelines for metagenomics data to investigate genes that enable electroactive metabolism and accompanying metabolisms. Our results are consistent with our hypothesis and yet upend some of our other expectations. Our findings have implications for understanding greenhouse gas emissions from lakes, including electroactivity as an integral part of lake metabolism throughout more of the anoxic parts of lakes and for a longer portion of the summer than expected. Our results also give a sense of what electroactivity occurs at given depths and provide a strong basis for future studies.

carbon emissions↗
DOE OSTI2022

Polymer Retention “Tailing” Phenomenon Associated with the Milne Point Polymer Flood

Summary Associated with the Milne Point polymer flood (on Alaska’s North Slope), this paper explores the unusual shape of hydrolyzed polyacrylamides (HPAM) breakout/propagation during dynamic polymer retention measurements in Milne Point core material. In contrast to conventional expectations, polymer retention does not delay the initial polymer arrival at the end of a Milne Point core. However, after effluent polymer concentrations rapidly rise to at least 50% of the injected value, the concentration gradually “tails” up over many pore volumes (PVs) before it finally achieves the injected value. To understand the origin and significance of this behavior, a wide range of core experiments were performed, including substantial variations in polymer concentration and molecular weight, core length, preservation state, sand grain size, and mineral composition. Illite was identified as primarily responsible for the tailing phenomenon. This phenomenon has important consequences that must be considered when projecting the performance of the field project. This work suggests that mineralogy analysis (especially for illite and kaolinite) may reveal whether tailing should be accounted for during simulations of polymer propagation/retention in a given field application.

Engineering↗