SEARCH · Engineering Papers

Results for “flux balance 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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2020

A machine learning approach to emulation and biophysical parameter estimation with the Community Land Model, version 5

Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibration targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Energetics of surface melt in West Antarctica

We use reanalysis data and satellite remote sensing of cloud properties to examine how meteorological conditions alter the surface energy balance to cause surface melt that is detectable in satellite passive microwave imagery over West Antarctica. This analysis can detect each of the three primary mechanisms for inducing surface melt at a specific location: thermal blanketing involving sensible heat flux and/or longwave heating by optically thick cloud cover, all-wave radiative enhancement by optically thin cloud cover, and föhn winds. We examine case studies over Pine Island and Thwaites glaciers, which are of interest for ice shelf and ice sheet stability, and over Siple Dome, which is more readily accessible for field work. During January 2015 over Siple Dome we identified a melt event whose origin is an all-wave radiative enhancement by optically thin clouds. During December 2011 over Pine Island and Thwaites glaciers, we identified a melt event caused mainly by thermal blanketing from optically thick clouds. Over Siple Dome, those same 2011 synoptic conditions yielded a thermal-blanketing-driven melt event that was initiated by an impulse of sensible heat flux and then prolonged by cloud longwave heating. The December 2011 synoptic conditions also generated föhn winds at a location on the Ross Ice Shelf adjacent to the Transantarctic Mountains, and we analyze this case with additional support from automatic weather station data. In contrast, a late-summer thermal blanketing period over Pine Island and Thwaites glaciers during February 2013 showed surface melt initiated by cloud longwave heating and then prolonged by enhanced sensible heat flux. One limitation thus far with this type of analysis involves uncertainties in the cloud optical properties. Nevertheless, with improvements this type of analysis can enable quantitative prediction of atmospheric stress on the vulnerable Antarctic ice shelves in a steadily warming climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Multidisciplinary Constraints on the Thermal‐Chemical Boundary Between Earth's Core and Mantle

Abstract Heat flux from the core to the mantle provides driving energy for mantle convection thus powering plate tectonics, and contributes a significant fraction of the geothermal heat budget. Indirect estimates of core‐mantle boundary heat flow are typically based on petrological evidence of mantle temperature, interpretations of temperatures indicated by seismic travel times, experimental measurements of mineral melting points, physical mantle convection models, or physical core convection models. However, previous estimates have not consistently integrated these lines of evidence. In this work, an interdisciplinary analysis is applied to co‐constrain core‐mantle boundary heat flow and test the thermal boundary layer (TBL) theory. The concurrence of TBL models, energy balance to support geomagnetism, seismology, and review of petrologic evidence for historic mantle temperatures supports Q CMB ∼15 TW, with all except geomagnetism supporting as high as ∼20 TW. These values provide a tighter constraint on core heat flux relative to previous work. Our work describes the seismic properties consistent with a TBL, and supports a long‐lived basal mantle molten layer through much of Earth's history.

58 GEOSCIENCES↗

DOE OSTI2026

Simulation of Divertor Performance in ST40 Under Dynamic Double-Null Plasmas

A power fraction model was implemented for the simultaneous prediction of 3-D surface temperature evolution at all four divertor targets in near-double-null (DN) tokamak configurations, which is especially important for compact high-field devices that may not have the ability to dissipate large amounts of power on the high-field side. Evaluating the power-sharing between the four divertor strike points in a disconnected DN configuration is important for understanding the overall power balance, as well as for optimizing the power exhaust performance and prolonging the survivability of the plasma facing components (PFCs). This power-sharing is typically evaluated in terms of the separation between the primary and secondary separatrices at the outboard midplane, $\textit {dR}_{\text {sep}}$. The Heat flux Engineering Analysis Toolkit (HEAT) is coupled with Brunner’s power fraction model to simulate the deposited heat flux and resultant temperature change on 3-D divertor targets in a dynamic DN (DDN) pulse operation in ST40, a high-field spherical tokamak. The simulation results showed that with DDN operation, the operation time has significantly increased compared with single-null geometry configurations.

ST40↗

DOE OSTI2022

Divertor plasma detachment: roles of plasma momentum, energy, and particle balances

Abstract Contrary to statements from some recent papers, the bifurcation of the scrape-off layer and divertor plasma parameters, as described within a 1D model, does exist. The analysis of both numerical simulations and experimental data shows that the real physical reasons that can cause a strong reduction of the plasma particle flux to the target j d are impurity radiation loss and volumetric plasma recombination. Fitting the parameter describing the ‘plasma momentum removal’, based on the results of 2D numerical simulations, then using this fit in subsequent analytic expressions for j d , and finally interpreting the evolution of j d as the result of elastic plasma-neutral interactions break the causality principle and lead to a wrong physical picture of divertor plasma detachment.

Physics↗

DOE OSTI2025

Air-sea interactions yielding rapid Beaufort Sea ice losses during the 2021 ONR THINICE Pilot Field Campaign

Seasonal sea ice coverage has declined in the Pacific Arctic region. Losses of the ice cover have modified the region’s surface heat and momentum fluxes providing weather systems more direct interaction with the upper ocean. For example, recent cyclones that have passed over these waters have initiated multiple air-sea interactions and have impacted mixed layer depth, wave properties, and sea ice mass budgets due to diabatic and/or advective processes. As a United States Office of Naval Research (ONR) Departmental Research initiative, one key component of the ONR THINICE program was to advance knowledge on Arctic cyclone coupling with the ocean surface, its ice cover, and understand how such coupling can lead to rapid ice loss. Leading up to the primary field campaign in August 2022, a pilot field campaign from 19 August to 13 September 2021 focused in Utqiagvik, Alaska and Longyearbyen, Svalbard was carried out. This period was marked by frequent Arctic cyclones, periods of strong low-level winds, and stints of extreme short-term Beaufort Sea ice area losses. Cyclones traversed these waters much more frequently than usual during this 26-day period, but did not exhibit extreme central pressures indicative of high intensity and impactful storms. Using ERA5 atmospheric reanalysis fields, CERES and AIRS surface flux retrievals, and in situ sensor data collected from ice-tethered profilers and ice mass balance buoys, we highlight air-sea-ice interactions in the Beaufort Sea during and between these cyclones leading up to extreme ice loss events. Our case study analysis is aimed at understanding how Arctic cyclones, including air-sea interactions prompted by their passage over the marginal ice zone, influence the end-of-summer ice conditions and ice loss events. This study provides an example where frequent, yet weak storms drive rather large marginal sea-scale losses albeit with heterogenous sea ice effects at local scales.

Ballinger, Thomas [International Arctic Research C↗

DOE OSTI2025

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗

DOE OSTI2022

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗

DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Rubisco regulation in response to altered carbon status in the cyanobacterium Synechococcus elongatus PCC 7942

Abstract Photosynthetic organisms possess a variety of mechanisms to achieve balance between absorbed light (source) and the capacity to metabolically utilize or dissipate this energy (sink). While regulatory processes that detect changes in metabolic status/balance are relatively well studied in plants, analogous pathways remain poorly characterized in photosynthetic microbes. Here, we explored systemic changes that result from alterations in carbon availability in the model cyanobacterium Synechococcus elongatus PCC 7942 by taking advantage of an engineered strain where influx/efflux of a central carbon metabolite, sucrose, can be regulated experimentally. We observed that induction of a high-flux sucrose export pathway leads to depletion of internal carbon storage pools (glycogen) and concurrent increases in estimates of photosynthetic activity. Further, a proteome-wide analysis and fluorescence reporter-based analysis revealed that upregulated factors following the activation of the metabolic sink are concentrated on ribulose-1,5-bisphosphate carboxylase-oxygenase (Rubisco) and auxiliary modules involved in Rubisco maturation. Carboxysome number and Rubisco activity also increased following engagement of sucrose secretion. Conversely, reversing the flux of sucrose by feeding exogenous sucrose through the heterologous transporter resulted in increased glycogen pools, decreased Rubisco abundance, and carboxysome reorganization. Our data suggest that Rubisco activity and organization are key variables connected to regulatory pathways involved in metabolic balancing in cyanobacteria.

Plant Sciences↗

DOE OSTI2024

Decadal increases in carbon uptake offset by respiratory losses across northern permafrost ecosystems

Tundra and boreal ecosystems encompass the northern circumpolar permafrost region and are experiencing rapid environmental change with important implications for the global carbon (C) budget. We analysed multi-decadal time series containing 302 annual estimates of carbon dioxide (CO 2 ) flux across 70 permafrost and non-permafrost ecosystems, and 672 estimates of summer CO 2 flux across 181 ecosystems. We find an increase in the annual CO 2 sink across non-permafrost ecosystems but not permafrost ecosystems, despite similar increases in summer uptake. Thus, recent non-growing-season CO 2 losses have substantially impacted the CO 2 balance of permafrost ecosystems. Furthermore, analysis of interannual variability reveals warmer summers amplify the C cycle (increase productivity and respiration) at putatively nitrogen-limited sites and at sites less reliant on summer precipitation for water use. Our findings suggest that water and nutrient availability will be important predictors of the C-cycle response of these ecosystems to future warming.

58 GEOSCIENCES↗

DOE OSTI2025

Cyanobacterial circadian regulation enhances bioproduction under subjective nighttime through rewiring of carbon partitioning dynamics, redox balance orchestration, and cell cycle modulation

Abstract Background The industrial feasibility of photosynthetic bioproduction using cyanobacterial platforms remains challenging due to insufficient yields, particularly due to competition between product formation and cellular carbon demands across different temporal phases of growth. This study investigates how circadian clock regulation impacts carbon partitioning between storage, growth, and product synthesis in Synechococcus elongatus PCC 7942, and provides insights that suggest potential strategies for enhanced bioproduction. Results After entrainment to light-dark cycles, PCC 7942 cultures transitioned to constant light revealed distinct temporal patterns in sucrose production, exhibiting three-fold higher productivity during subjective night compared to subjective day despite moderate down-regulation of genes from the photosynthetic apparatus. This enhanced productivity coincided with reduced glycogen accumulation and halted cell division at subjective night time, suggesting temporal separation of competing processes. Transcriptome analysis revealed coordinated circadian clock-driven adjustment of the cell cycle and rewiring of energy and carbon metabolism, with over 300 genes showing differential expression across four time points. The subjective night was characterized by altered expression of cell division-related genes and reduced expression of genes involved in glycogen synthesis, while showing upregulation of glycogen degradation pathways, alternative electron flow components, the pentose phosphate pathway, and oxidative decarboxylation of pyruvate. These molecular changes created favorable conditions for product formation through enhanced availability of major sucrose precursors (glucose-1-phosphate and fructose-6-phosphate) and maintained redox balance through multiple mechanisms. Conclusions Our analysis of circadian regulatory rewiring of carbon metabolism and redox balancing suggests two potential approaches that could be developed for improving cyanobacterial bioproduction: leveraging natural circadian rhythms for optimizing cultivation conditions and timing of pathway induction, and engineering strains that mimic circadian-driven metabolic shifts through controlled carbon flux redistribution and redox rebalancing. While these strategies remain to be tested, they could theoretically improve the efficiency of photosynthetic bioproduction by enabling better temporal separation between cell growth, carbon storage accumulation, and product synthesis phases.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

The role of plasma–atom and molecule interactions on power & particle balance during detachment on the MAST Upgrade Super-X divertor

This paper shows first quantitative analysis of the detachment processes in the MAST Upgrade Super-X divertor (SXD). We identify an unprecedented impact of plasma-molecular interactions involving molecular ions (likely $D^+_2$), resulting in strong ion sinks (Molecular Activated Recombination—MAR), leading to a reduction of ion target flux. The MAR ion sinks exceed the divertor ion sources before electron-ion recombination (EIR) starts to occur, suggesting that significant ionisation occurs outside of the divertor chamber. In the EIR region, $T_e \ll 0.2$ eV is observed and MAR remains significant in these deep detached phases. The total ion sink strength demonstrates the capability for particle (ion) exhaust in the Super-X Configuration. Molecular Activated Dissociation is the dominant volumetric neutral atom creation process can lead to an electron cooling of 20% of $P_\mathrm{SOL}$. The measured total radiative power losses in the divertor chamber are consistent with inferred hydrogenic radiative power losses. This suggests that intrinsic divertor impurity radiation, despite the carbon walls, is minor in the divertor chamber. This contrasts previous TCV results, which may be associated with enhanced plasma-neutral interactions and reduced chemical erosion in the detached, tightly baffled SXD. The above observations have also been observed in higher heat flux (narrower SOL width) type I ELMy H-mode discharges. This provides evidence that the characterisation in this paper may be general.

plasma-molecular interactions↗

DOE OSTI2025

Dynamic soil columns simulate Arctic redox biogeochemistry and carbon release during changes in water saturation

Thawing Arctic permafrost can induce hydrologic change and alter redox conditions, shifting the balance of soil organic matter (SOM) decomposition. There remains uncertainty about how soil saturation and redox transitions impact dissolved and gas phase carbon fluxes, and efforts to link hydrobiogeochemical processes to ecosystem-scale models are limited. This study evaluates SOM decomposition of Arctic tundra soils using column experiments, water chemistry measurements, microbial community analysis, and a PFLOTRAN reactive transport model. Soil columns from a thermokarst channel (TC) and an upland tundra (UC) were exposed to cycles of saturation and drainage, which controlled carbon emissions. During saturation, an outflow of dissolved organic carbon from the UC soil correlated with elevated reduced iron and decreased pH; during drainage, UC carbon dioxide fluxes were 70% higher than TC fluxes. Intermittent methane release was observed for TC, consistent with higher methanogen abundance. Slower drainage in the TC soil correlated with more subtle biogeochemical changes. PFLOTRAN simulations captured experimental trends in soil carbon fluxes, oxygen concentrations, and water contents. The model was then used to evaluate additional soil water drainage rates. This study emphasizes the importance of considering hydrologic change when evaluating and simulating SOM decomposition in dynamic Arctic tundra environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Particle balance of deuterium during deuterium shattered pellet injection shutdown in DIII-D

A particle balance analysis was conducted during a deuterium (D 2 ) shattered pellet injection-induced plasma shutdown on the DIII-D tokamak to determine why less than 20% of the pellet material is assimilated into the core plasma by the mid-current quench (CQ). Initially, most of the D 2 is injected as frozen shards and ionized upon entering the vessel. During the thermal quench, ionized particles move to the divertors and subsequently to the center post (CP) walls, where they rapidly recycle and partially accumulate as neutrals without assimilating into the core plasma. In contrast, the particle flux to the outer midplane walls is negligible, despite being accompanied by hot plasma with electron temperatures exceeding 100 eV. During mid-CQ, volume recombination effects, although not large enough to impact overall particle balance, were significant enough to require accounting for accurate interpretation of fast-framing camera D-alpha signals and the estimation of the CP wall particle flux. In addition, toroidal asymmetries, observed in measurements of toroidal electron density perturbations and the phase of magnetohydrodynamic modes, are present throughout the shutdown and can account for a discrepancy in the assimilation rate for up to 50% of the observed D 2 particle inventory. These sources and sinks of particles and fluxes were identified using absolutely calibrated D-alpha brightness and Langmuir probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

The effects of photosynthetic rate on respiration in light, starch/sucrose partitioning, and other metabolic fluxes within photosynthesis

In the future, plants may encounter increased light and elevated CO 2 levels. How consequent alterations in photosynthetic rates will impact fluxes in photosynthetic carbon metabolism remains uncertain. Respiration in light ( R L ) is pivotal in plant carbon balance and a key parameter in photosynthesis models. Understanding the dynamics of photosynthetic metabolism and R L under varying environmental conditions is essential for optimizing plant growth and agricultural productivity. However, measuring R L under high light and high CO 2 (HLHC) conditions poses challenges using traditional gas exchange methods. In this study, we employed isotopically nonstationary metabolic flux analysis (INST-MFA) to estimate RL and investigate photosynthetic carbon flux, unveiling nuanced adjustments in Camelina sativa under HLHC. Despite numerous flux alterations in HLHC, RL remained stable. HLHC affects several factors influencing RL, such as starch and sucrose partitioning, v o /v c ratio, triose phosphate partitioning, and hexose kinase activity. Analysis of A/C i curve operational points reveals that HLHC’s major changes primarily stem from CO 2 suppressing photorespiration. Integration of these fluxes into a simplified model predicts changes in CBC labeling under HLHC. This study extends our prior discovery that incomplete CBC labeling is due to unlabeled carbon reimported during R L , offering insights into manipulating labeling through adjustments in photosynthetic rates.

Elevated CO2↗

DOE OSTI2025

A model for selecting the best sustainable airport technology alternatives

Air transport is a continually expanding industry, a fact that has become even more evident with the recovery of the aviation industry post-COVID-19. This expansion amplified energy consumption at airports, which was already significantly high. Airport decision-makers are increasingly focusing on improving sustainability and addressing social, economic, and environmental criteria across airports worldwide. In this article, we propose a decision-making model for identifying a sustainable airport technology solution that minimizes the energy consumption of airport lighting while considering economic, emission, and life-cycle criteria. The proposed model combines data envelopment analysis and multi-criteria decision making techniques. We applied the model to the Dallas/Fort Worth International Airport (DFW) as a case study, considering eight lighting technology solutions, each with five luminous flux alternatives. Our model identified the best lighting technology solution for DFW outdoor and indoor environments based on the following criteria: luminous flux, capital costs, life-cycle costs, energy consumption, and emissions (CO2e, NOx, SO2, and PM2.5). The designed model is customizable to any airport and is applicable to a wide range of airport lighting technologies. In our analysis, Light-Emitting Diode lighting emerged as the most sustainable technology option. It ranked first in most cases due to its balance of high efficacy, long lifespan, low life-cycle cost, low capital cost, and lower emissions across all pollutants.

Tchivwila, Moise B↗

DOE OSTI2021

Computational Framework for Machine-Learning-Enabled 13 C Fluxomics

13 C metabolic flux analysis (MFA) has emerged as a powerful tool for synthetic biology. This optimization-based approach suffers long computation time and unstable solutions depending on the initial guess. Here, we develop a machine-learning-based framework for 13 C fluxomics. Specifically, training and test data sets are generated by metabolic network decomposition and flux sampling, in which flux ratios at metabolic nodes and simulated labeling patterns of metabolites are used as training targets and features, respectively. To improve prediction accuracy and simplify the model, automated processes are developed for flux ratio selection based on solvability and feature screening based on importance. We found that predictive performance can be significantly improved using both amino acids and central carbon metabolites in comparison with amino acids alone. Together with measured external fluxes, the predicted flux ratios determine the mass balance system, yielding global flux distributions. This approach is validated by flux estimation using both simulated and experimental data in comparison with canonical 13 C MFA. The approach represents a reliable fluxomics method readily applicable to high-throughput metabolic phenotyping, which highlights the advances of intelligent learning algorithms in synthetic biology, specifically in the Test and Learn stage of the Design-Build-Test-Learn cycle.

13C metabolic flux analysis↗