Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Predictive Control”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Advanced Fuels Campaign Execution Plan

The Advanced Fuels Campaign (AFC) Execution Plan outlines the strategy, mission, scope, near-term and long-term goals, structure, and organization associated with nuclear fuels and materials research, development, and demonstration activities within the Department of Energy’s (DOE) Nuclear Fuel Cycle and Supply Chain (NFCSC) program. NFCSC has been given responsibility to identify and mature advanced fuel technologies for the DOE using a science-based approach, focused on developing a fundamental understanding of nuclear fuels and materials to drive development of integrated nuclear fuel and materials technology. This science-based approach combines theory, experiments, and multiscale modeling and simulation to achieve a predictive understanding of relevant behaviors ranging from fuel fabrication processes (and their resulting fuel microstructures) through fuel/cladding performance under irradiation (in contrast to more empirical, observation-based approaches frequently used in fuel performance modeling and fuel qualification). The traditional scope of AFC includes the evaluation and development of multiple fuel forms to support two fuel cycle options: once-through and full recycle. The word “fuel” is used generically to include conventional fuels, transmutation targets, and any associated cladding or duct materials. The once-through fuel cycle addresses advanced light water reactor fuels with enhanced performance, extended burnup, and reduced waste generation. In fiscal year (FY) 2012, AFC’s scope expanded to include research, development, and demonstration (RD&D) for light water reactor (LWR) fuels with enhanced accident tolerance. Fuel fabrication activities include the development of innovative methods to enhance process efficiencies, reduce waste, and improve control over as-fabricated fuel microstructural properties to achieve desired in-reactor performance. Using modern modeling and simulation approaches, the objective is to predict fresh fuel properties given the feedstock characteristics and fabrication process parameters. The performance-related activities include small-scale, in-reactor, and out-of-reactor phenomenological testing (distinct from, but synergistic with, integral prototypic testing) and extensive, quantitative characterization (focusing on characterization of fuel and cladding materials at the scale of microstructure) both before and after testing. Larger-scale, prototypic experiments are conducted in concert with phenomenological testing to drive a Fuel Development and Qualification program, incorporating a fundamental understanding of fuel behavior performance characteristics. Then, using the tools developed under the productive science-based approach, fuels will be optimized to meet specific performance requirements, thereby minimizing the need to repeatedly perform large-scale, integral experiments over a wide parametric range as a means of experimental exploration. Two significant initiatives are underway within AFC. First, a gap analysis completed in early FY 2019 identified critical irradiation testing needs that are lacking within the national light water reactor (LWR) fuels testbed since the shutdown of the Halden Reactor in 2018. The identified gaps are for instrumented, prototypic testing of LWR fuels, especially under boiling water reactor conditions, ramp conditions, and conditions leading to fuel failure; these needs exist for supporting current LWR fuels and their possible extension to higher burnups, but are especially urgent relative to near-term development and qualification of accident-tolerant fuels. Recommendations that resulted from the Halden Gap Analysis focused on enhancements at Advanced Test Reactor (ATR) and Transient Reactor Test Facility (TREAT) to fill gaps in testing capabilities relative to these needs. Second, a concerted effort to develop and demonstrate a systematic approach to accelerating the development, testing, and qualification of new fuel systems has been initiated. This is highlighted by a test strategy that combines the considerable advances in multiscale, mechanistic fuel modeling of recent years with a MiniFuel separate effects test program in the High Flux Isotope Reactor (HFIR) and a Fission Accelerated Steady-state Testing (FAST) semi-integral accelerated test program in ATR. This approach is being tested/demonstrated using the metallic fuel system, but if successful it is expected to be extensible to multiple fuel types and diverse applications. This document includes an overview of the NFCSC program, a definition of science-based development of nuclear fuels, near-term goals for Advanced LWR fuels (ALFs), and longer-term goals for Advanced Reactor Fuels (ARFs) RD&D. This includes the activities that will be conducted to achieve success toward the grand challenge, as well as the goals and milestones to be achieved over the next few decades of research and development. Long-term goals are based on the DOE Office of Nuclear Energ

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Relation between interfacial shear and friction force in 2D materials

Understanding the interfacial properties between an atomic layer and its substrate is of key interest at both the fundamental and technological levels. From Fermi level pinning to strain engineering and superlubricity, the interaction between a single atomic layer and its substrate governs electronic, mechanical and chemical properties. Here, we measure the hardly accessible interfacial transverse shear modulus of an atomic layer on a substrate. By performing measurements on bulk graphite, and on epitaxial graphene films on SiC with different stacking orders and twisting, as well as in the presence of intercalated hydrogen, we find that the interfacial transverse shear modulus is critically controlled by the stacking order and the atomic layer–substrate interaction. Importantly, we demonstrate that this modulus is a pivotal measurable property to control and predict sliding friction in supported two-dimensional materials. The experiments demonstrate a reciprocal relationship between friction force per unit contact area and interfacial shear modulus. As a result, the same relationship emerges from simulations with simple friction models, where the atomic layer– substrate interaction controls the shear stiffness and therefore the resulting friction dissipation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Prediction of defect properties in concentrated solid solutions using a Langmuir-like model

The alleged existence of sluggish diffusion in high-entropy alloys has drawn controversy. In high-entropy alloys and, in general, in all solids, transport properties are controlled by point defect concentration, which must be known before performing atomistic simulations to compute transport coefficients. In this work, we present a general Langmuir-like model for defect concentration in an arbitrarily complex solid solution and apply this model to generate expressions for concentrations of vacancies and small interstitial atoms. We then calculate the vacancy concentration as a function of temperature in the equiatomic CoNiCrFeMn and FeAl alloys with modified embedded-atom-method potentials for various chemical orderings, showing there is no clear correlation between vacancy thermodynamics and chemical ordering in the CoNiCrFeMn alloy, but clear systematic patterns for FeAl. We believe this is due to the high stability of disordered, random, and ordered intermetallic phases, respectively, in the CoNiCrFeMn and FeAl systems. Finally, this work provides future avenues to the prediction of thermal interstitials and vacancies in solid solutions, which is necessary for models of nonequilibrium behavior of solid solutions.

composition↗
DOE OSTI2023

A hybrid analytical and numerical model for cross-over and performance decay in a unit cell vanadium redox flow battery

Developing an accurate and efficient model for cross-over is critical for improving the long-term performance of redox flow batteries (RFBs). Here, this work presents a hybrid analytical and numerical model that combines a two-dimensional analytical solution to the active species, a one-dimensional analytical model for cross-over mechanisms, and a zero-dimensional numerical model for outlet concentrations of reactants. By comparing with experiment over 41 cycles (ca. 144 h), the model reported a mean voltage difference of 0.0089 V, mean time difference of 14 s per cycle (of 3.5 h), maximum relative difference for capacity and energy of 1.34% and 1.63%, respectively. The predicted mean concentrations for V 2+ , V 3+ , and VO$_{2}^{+}$ in membrane are 65%~77% of measured values from the literature. Upon validation, the model reproduced behaviors in electrolyte imbalance similar to those observed in experiments and numerical models, and revealed the control of cross-over, self-discharge, stoichiometry of side reactions, and Coulombic efficiency on electrolyte imbalance. The model also demonstrates excellent computational efficiency for simulating 41 cycles (around 37,000 points) within 3~4 s. The demonstrated efficiency and accuracy for predicting cross-over and its impacts on voltage, capacity&energy decay, membrane concentrations, and electrolyte imbalance makes it a reliable tool for optimizing RFBs’ long-term performance.

25 ENERGY STORAGE↗
DOE OSTI2023

Chemical Blast Standard (1 kg)

Chemical explosions create blast waves with large overpressure disturbances. It is important to develop a standard blast model based on data to accurately predict acoustic blast-wave amplitudes near detonations and invert for explosion energy from distant observations of blast-wave signals. However, open data from large, controlled chemical explosions with reliable ground truth can be challenging to find. The lack of access to such data could limit the number of contributions to related research and potentially stifle the rate of discoveries or validation of existing models. Here, to address these data scarcity problem, we have curated and compiled a standardized set of 817 blast-wave waveforms from 19 distinct high-explosive events. The blast-wave waveforms are standardized to a 1 kg trinitrotoluene explosion using scaling laws and corrections for location effects. A brief overview of the dataset is presented along with explosion feature models as well as recommendations for extracting explosion features. The resulting dataset is distributed to an open repository in both Seismic Analysis Code and pandas DataFrame formats containing the waveforms, the scaled distances, and the sample rates.

58 GEOSCIENCES↗
DOE OSTI2020

Power Prediction of Airborne Wind Energy Systems Using Multivariate Machine Learning

Kites can be used to harvest wind energy at higher altitudes while using only a fraction of the material required for conventional wind turbines. In this work, we present the kite system of Kyushu University and demonstrate how experimental data can be used to train machine learning regression models. The system is designed for 7 kW traction power and comprises an inflatable wing with suspended kite control unit that is either tethered to a fixed ground anchor or to a towing vehicle to produce a controlled relative flow environment. A measurement unit was attached to the kite for data acquisition. To predict the generated tether force, we collected input–output samples from a set of well-designed experimental runs to act as our labeled training data in a supervised machine learning setting. We then identified a set of key input parameters which were found to be consistent with our sensitivity analysis using Pearson input–output correlation metrics. Finally, we designed and tested the accuracy of a neural network, among other multivariate regression models. The quality metrics of our models show great promise in accurately predicting the tether force for new input/feature combinations and potentially guide new designs for optimal power generation.

17 WIND ENERGY↗
DOE OSTI2024

Towards accurate prediction of configurational disorder properties in materials using graph neural networks

Abstract The prediction of configurational disorder properties, such as configurational entropy and order-disorder phase transition temperature, of compound materials relies on efficient and accurate evaluations of configurational energies. Previous cluster expansion methods are not applicable to configurationally-complex material systems, including those with atomic distortions and long-range orders. In this work, we propose to leverage the versatile expressive capabilities of graph neural networks (GNNs) for efficient evaluations of configurational energies and present a workflow combining attention-based GNNs and Monte Carlo simulations to calculate the disorder properties. Using the dataset of face-centered tetragonal gold copper without and with local atomic distortions as an example, we demonstrate that the proposed data-driven framework enables the prediction of phase transition temperatures close to experimental values. We also elucidate that the variance of the energy deviations among configurations controls the prediction accuracy of disorder properties and can be used as the target loss function when training and selecting the GNN models. The work serves as a fundamental step toward a data-driven paradigm for the accelerated design of configurationally-complex functional material systems.

Chemistry↗
DOE OSTI2023

Band gap predictions of double perovskite oxides using machine learning

Abstract The compositional and structural variety inherent to oxide perovskites spawn wide-ranging applications. In perovskites, the band gap E g , a key material parameter for these applications, can be optimally controlled by varying the composition. Here, we implement a hierarchical screening process in which two cross-validated and predictive machine learning models for band gap classification and regression, trained using exhaustive datasets that span 68 elements of the periodic table, are applied sequentially. The classification model separates wide band gap materials, with E g ≥ 0.5 eV, from materials which have zero or relatively small band gaps, namely E g < 0.5 eV, and the second regression model quantitatively predicts the gap value of the wide band gap compounds. The study down-selects 13,589 cubic oxide perovskite compositions that are predicted to be experimentally formable, thermodynamically stable, and have a wide band gap. Of these, a subset of 310 compounds, which are predicted to be stable and formable with a confidence greater than 90%, are identified for further investigation. Our models are methodically analyzed via performance metrics and inter-dependence of model features to gain physical insight into the band gap prediction problem. Design maps to identify the variation of band gap with substitution of different elements are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-throughput protein characterization by complementation using DNA barcoded fragment libraries

Abstract Our ability to predict, control, or design biological function is fundamentally limited by poorly annotated gene function. This can be particularly challenging in non-model systems. Accordingly, there is motivation for new high-throughput methods for accurate functional annotation. Here, we used co mplementation of aux otrophs and DNA barcode seq uencing (Coaux-Seq) to enable high-throughput characterization of protein function. Fragment libraries from eleven genetically diverse bacteria were tested in twenty different auxotrophic strains of Escherichia coli to identify genes that complement missing biochemical activity. We recovered 41% of expected hits, with effectiveness ranging per source genome, and observed success even with distant E. coli relatives like Bacillus subtilis and Bacteroides thetaiotaomicron . Coaux-Seq provided the first experimental validation for 53 proteins, of which 11 are less than 40% identical to an experimentally characterized protein. Among the unexpected function identified was a sulfate uptake transporter, an O-succinylhomoserine sulfhydrylase for methionine synthesis, and an aminotransferase. We also identified instances of cross-feeding wherein protein overexpression and nearby non-auxotrophic strains enabled growth. Altogether, Coaux-Seq’s utility is demonstrated, with future applications in ecology, health, and engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Computational Fluid Dynamics-Based Modeling of Methane Flows Around Oil and Gas Equipment

Recent studies estimate that emissions from oil and gas production facilities contribute between 20 and 50% of the total methane ( CH 4 ) emitted in the US; therefore, quantifying and reducing these emissions are crucial for achieving climate goals. Methane quantification depends on both measuring methane concentrations and converting them to emissions through a modeling framework. Currently, simple atmospheric dispersion models are primarily used to quantify emissions and concentrations, but these estimates are highly uncertain when quantifying emissions from complex aerodynamic sources, such as oil and gas facilities. This investigation used a CFD modeling approach, which can account for aerodynamic complexity but has hitherto not been used to model methane concentrations downwind of a methane release of a known rate, and compared it against in situ measurements. High-time-resolution (1 Hz) methane concentration and meteorological data were measured during experiments conducted at the METEC on 21 March and 11 July 2024. The METEC site configuration, measured wind data, and controlled emission rates were used as input for the CONVERGE CFD model to model downwind CH 4 concentration. The modeling was carried out between 20 and 70 m, from two different points of release in two separate controlled-release experiments, one from a separator and another from a wellhead. In these experiments, we found that the CFD model could predict the CH 4 concentrations downwind of the release to a good degree. The model was evaluated on multiple metrics to assess its performance in estimating methane concentrations at typical fence line distances (∼30 m). These results help us to understand external flows and the ability of CFD models to predict downwind concentrations in aerodynamically complex environments.

03 NATURAL GAS↗
DOE OSTI2025

Operation at Reduced Atmospheric Pressure and Concept of Reliability Redundancy for Optimized Design of Insulation Systems

Electrified transportation is calling for insulation design criteria that is adequate to provide elevated levels of power density, power dynamics and reliability. Increasing voltage levels are expected to cause accelerated intrinsic and extrinsic aging effects which will not be easily predictable at the design stage due to a lack of suitable modeling. Designing reliable insulation systems would require finding solutions able to control accelerated aging due to an unpredictable increase of intrinsic stresses and the onset of extrinsic stresses as partial discharges. This paper proposes the concept of reliability redundancy for the insulation design of aerospace electrical asset components, which is also validated at lower-than-standard atmospheric pressure. The principle is that extrinsic-aging-free design might be achieved upon determining the aging stress or abnormal service stresses distribution and being sure that aging will not generate conditions that can incept extrinsic aging (partial discharges) during operation life. However, such information is never, in practice, fully available to insulation system designers. Hence, especially in critical applications such as electrified aircraft, aerospace, and combat ships a further level of reliability should be added to a partial-discharge-free design, which can consist of the use of corona-resistant materials and/or of life models able to consider the accelerated aging effect of partial discharges (or any other type of extrinsic-accelerated aging factor). Innovative life modeling considering both extrinsic and intrinsic aging stresses, insulating material testing to estimate model parameters, and a metric for quantifying the extent of corona (or partial discharge) resistance can lead to establishing feasibility and limit conditions for optimized or fully reliability-redundant design. It is shown in the paper that if an extrinsic-aging-free design is not feasible, and it is therefore replaced by a redundant design, a further level of reliability redundancy can be provided by effective condition monitoring plans.

Montanari, Gian Carlo↗
DOE OSTI2024

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗
DOE OSTI2024

Experimental and Computational Studies of Crystal Nucleation in Composition Gradients (Final Report)

The major goals of this project were to develop and validate a predictive nucleation model that incorporates composition gradients and is applicable to a large range of metallic systems which form both stochiometric and non-stochiometric compounds. Computationally we planned to expand the extent of thermodynamics-based theories and lay the groundwork to improve the general understanding and predictive capabilities of the role that gradients play in phase formation, glass formability and stability, and nucleation and growth events. To address these goals, we used Molecular Dynamics (MD) simulations and both isothermal and isochronal nanocalorimetric experiments on amorphous phases with a controlled composition gradient with the hopes of validating and improving the classical nucleation model and its use in solid solutions. While the computations were successful and identified ways to improve the classical nucleation model, the in situ nanocalorimetry studies proved very challenging due to unexpected difficulties in fabricating effective calorimeters and samples. Thus, we could not experimentally validate our predicted influence of composition gradients on nucleation. Nonetheless, important insights were gained and modifications to the classical nucleation theory were suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simultaneously improving accuracy and computational cost under parametric constraints in materials property prediction tasks

Abstract Modern data mining techniques using machine learning (ML) and deep learning (DL) algorithms have been shown to excel in the regression-based task of materials property prediction using various materials representations. In an attempt to improve the predictive performance of the deep neural network model, researchers have tried to add more layers as well as develop new architectural components to create sophisticated and deep neural network models that can aid in the training process and improve the predictive ability of the final model. However, usually, these modifications require a lot of computational resources, thereby further increasing the already large model training time, which is often not feasible, thereby limiting usage for most researchers. In this paper, we study and propose a deep neural network framework for regression-based problems comprising of fully connected layers that can work with any numerical vector-based materials representations as model input. We present a novel deep regression neural network, iBRNet, with branched skip connections and multiple schedulers, which can reduce the number of parameters used to construct the model, improve the accuracy, and decrease the training time of the predictive model. We perform the model training using composition-based numerical vectors representing the elemental fractions of the respective materials and compare their performance against other traditional ML and several known DL architectures. Using multiple datasets with varying data sizes for training and testing, We show that the proposed iBRNet models outperform the state-of-the-art ML and DL models for all data sizes. We also show that the branched structure and usage of multiple schedulers lead to fewer parameters and faster model training time with better convergence than other neural networks. Scientific contribution: The combination of multiple callback functions in deep neural networks minimizes training time and maximizes accuracy in a controlled computational environment with parametric constraints for the task of materials property prediction.

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

Strain‐Driven Stabilization of a Room‐Temperature Chiral Multiferroic with Coupled Ferroaxial and Ferroelectric Order

Noncollinear ferroic materials are sought after as testbeds to explore the intimate connections between topology and symmetry, which result in electronic, optical, and magnetic functionalities not observed in collinear ferroic materials. For example, ferroaxial materials have rotational structural distortions that break mirror symmetry and induce chirality. When ferroaxial order is coupled with ferroelectricity arising from a broken inversion symmetry, it offers the prospect of electric-field-control of the ferroaxial distortions and opens up new tunable functionalities. However, chiral multiferroics, especially ones stable at room temperature, are rare. A strain-stabilized, room-temperature chiral multiferroic phase in single crystals of BaTiS 3 is reported here. Using first-principles calculations, the stabilization of this multiferroic phase having P6 3 space group for biaxial tensile strains exceeding 1.5% applied on the basal ab-plane of the room temperature P6 3 cm phase of BaTiS 3 is predicted. The chiral multiferroic phase is characterized by rotational distortions of TiS 6 octahedra around the long c-axis and polar displacement of Ti atoms along the c-axis. Furthermore, the ferroaxial and ferroelectric distortions and their domains in P6 3 -BaTiS 3 are directly resolved using atomic resolution scanning transmission electron microscopy. Landau-based phenomenological modeling predicts a strong coupling between the ferroelectric and the ferroaxial order making P6 3 -BaTiS 3 an attractive test bed for achieving electric-field-control of chirality.

Chalcogenide↗
DOE OSTI2026

N 2 Onet: a global collaborative network facilitating advances in measurement, modeling, and mitigation of agricultural soil nitrous oxide emissions

Nitrogen (N) fertilizer supports global food production, but its use and overuse drive emissions of nitrous oxide (N 2 O), a potent and long-lived greenhouse gas. Understanding the drivers of N 2 O fluxes remains elusive, making it difficult to predict emissions in time and space and to develop and evaluate ways to lower emissions through management. Major scientific uncertainties underlying the understanding of the drivers of N 2 O fluxes identified in a workshop of N 2 O emissions experts include poor process-based understanding of controls on soil N 2 O emissions in the field; insufficient data to reduce uncertainty in N 2 O budgets from the field to regional scales, including N 2 O emission measurements and importantly, field-scale N balances; and high uncertainty in model predictions of soil N 2 O emissions across environmental and management conditions. To reduce these uncertainties, we present the concept of N 2 Onet, a global collaborative initiative to accelerate advances in N 2 O measurement, analyses, and mitigation. N 2 Onet will serve as an observational network of supersites with multi-scale measurements; a database hub for N 2 O flux and ancillary data; and a catalyst for community building, information sharing, and training. By coalescing and coordinating the global community of researchers, N 2 Onet will provide a roadmap for reducing N 2 O emissions from agriculture worldwide.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Enhancing nanoscale charged colloid crystallization near a metastable liquid binodal

Achieving predictive control over crystallization using non-classical nucleation while avoiding kinetic traps would be a step towards designing materials with new functionalities. We address these challenges by inducing the bottom-up assembly of nanocrystals into ordered arrays, or superlattices. Using electrostatics—rather than density—to tune the interactions between particles, we watch self-assembly proceed through a metastable liquid phase. Here, we systematically investigate the phase behaviour as a function of quench conditions in situ and in real time using small-angle X-ray scattering. By fitting to colloid, liquid and superlattice models, we extract the time evolution of each phase and the system phase diagram, which we find to be consistent with short-range attractive interactions. Using the predictive power of the phase diagram, we establish control of the self-assembly rate over three orders of magnitude, and we identify one- and two-step self-assembly regimes, with only the latter implicating the metastable liquid as an intermediate. The presence of the metastable liquid increases the superlattice formation rate relative to the equivalent one-step pathway, and the superlattice order increases with the rate, revealing a generalizable kinetic strategy for promoting and enhancing ordered assembly.

Tanner, Christian P. N. [University of California,↗
DOE OSTI2022

Spatiotemporal Metabolic Network Models Reveal Complex Autotroph-Heterotroph Biofilm Interactions Governed by Photon Incidences

Autotroph-heterotroph interactions are ubiquitous in natural environment and play a key role in controlling various essential ecosystem functions, such as production and utilization of organic matter, cycling of nitrogen, sulfur, and other chemical elements. Understanding how these biofilm metabolic interactions are constrained in space and time remains challenging because fully predictive models designed for this purpose are currently limited. Toward filling this gap, here we developed community metabolic network models for two autotroph-heterotroph biofilm consortia (termed UCC-A and UCC-O), which share a suite of common heterotrophic members but have a single distinct photoautotrophic cyanobacterium (Phormidesmis priestleyi str. ANA and Phormidium sp. OSCR) that provides organic carbon and nitrogen sources to support the growth of heterotrophic partners. After determining model parameters by data fitting using the spatiotemporal distributions of microbial abundances, we comparatively analyzed the resulting biofilm models to examine any fundamental differences in microbial interactions between the two consortia under the variation of key environmental variables: CO2 and photon levels. The UCC-A model predicted generally expected responses, i.e., the autotroph population increased in response to elevated levels of CO2 and photon, followed by increase in the heterotroph population. In contrast, the UCC-O model showed somewhat complicated dynamics, e.g., higher photon incidence rates resulted in the increase in autotroph population but decrease in heterotroph population due to the lowered provision of glucose from the autotroph. A further analysis showed that species coexistence was governed by the photon incidences rather than the carbon availability for UCC-O, which was the opposite for UCC-A.

Phalak, Poonam↗