Search NASA⌕ Search

SEARCH · Search NASA

Results for “activity coefficient model”

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

Comment on “Thermodynamic Models for the (HClO 4 + NaClO 4 ){aq} and (HBr + NaBr){aq} Systems at 298.15 K and 0.1 MPa” Authored by Oakes, C. S., Ward, A. L., Chugunov, N. Journal of Chemical & Engineering Data , 68 , 2554–2562

Oakes et al. (2023) published a review article in this journal. In that paper, Oakes et al. (2023) developed thermodynamic models to describe electrolyte solutions for HClO 4 –NaClO 4 –H 2 O and HBr–NaBr–H 2 O systems, based on literature data. In their paper, previously published work from researchers in the field was criticized; some of it is ours. Here, in this brief Comment, we first comment on their models, and then we briefly provide a technical response to that criticism.

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

A Theoretical Model for Computing Freezing Point Depression of Lithium-Ion Battery Electrolytes

Reliable prediction of freezing point depression in liquid electrolytes will accelerate the development of improved Li-ion batteries which can operate in low temperature environments. In this work we establish a computational methodology to calculate activity coefficients and liquidus lines for battery-relevant liquid electrolytes. Electronic structure methods are used in conjuction with classical molecular dynamics simulations and theoretical expressions for Born solvation energy, ion-atmosphere effects from Debye-Hückel theory and solvent entropic effects. The framework uses no a priori knowledge beyond neat solvent properties and the concentration of salt. LiPF 6 in propylene carbonate (PC), LiPF 6 in dimethyl carbonate (DMC) and LiClO 4 in DMC are investigated up to 1 molal with accuracy better than 3 °C when compared to experimental freezing point measurements. Finally, we find that the difference in freezing point depression between the propylene carbonate-based electrolyte and the dimethyl carbonate electrolytes originates from the difference in the solvent dielectric constant.

25 ENERGY STORAGE↗
DOE OSTI2021

Biomass pyrolysis devolatilization kinetics of herbaceous and woody feedstocks

Devolatilization kinetics were determined using a modified micropyrolyzer reactor for several biomass feedstocks: switchgrass, corn stover, red oak, and pine. Here, the micropyrolyzer was directly coupled to a flame ionization detector (FID) to track the release of volatiles from the biomass. Time series data from these experiments was analyzed to determine apparent devolatilization rates. Care was taken to assure the experiments were isothermal and kinetically limited calculating a Biot number less than 0.1 and pyrolysis numbers greater than 10, which simplifies the derivation of devolatilization rates. A single, first order reaction was able to model devolatilization rates at temperatures up to 500 °C. No correlation was found between the inorganic content of the biomass and its rate of devolatilization. Apparent activation energies were in the range of 54.9–88.4 kJ mol –1 . The rate coefficient at 500 °C was calculated as 1.90–5.14 s –1 for the four feedstocks.

09 BIOMASS FUELS↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs): An Efficient Physics-Informed Neural Networks Framework for Interface Problems

Here, we present an efficient physics-informed neural networks (PINNs) framework, termed Adaptive Interface-PINNs (AdaI-PINNs), to improve the modeling of interface problems with discontinuous coefficients and/or interfacial jumps. This framework is an enhanced version of its predecessor, Interface PINNs or I-PINNs (Sarma et al.; https://doi.org/10.1016/j.cma.2024.117135), which involves domain decomposition and assignment of different predefined activation functions to the neural networks in each subdomain across a sharp interface, while keeping all other parameters of the neural networks identical. In AdaI-PINNs, the activation functions vary solely in their slopes, which are trained along with the other parameters of the neural networks. This makes the AdaI-PINNs framework fully automated without requiring preset activation functions. Comparative studies on one-dimensional, two-dimensional, and three-dimensional benchmark elliptic interface problems reveal that AdaI-PINNs outperform I-PINNs, reducing computational costs by 2-6 times while producing similar or better accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Interpretation of Fracture Initiation Points by In-Well Low-Frequency Distributed Acoustic Sensing in Horizontal Wells

Summary Low-frequency distributed acoustic sensing (LF-DAS) exploits the optical phase shift of Rayleigh backscatter in fiber-optic cables to obtain distributed measurements of changes in strain and temperature. Fiber-optic cables are often installed for multistage hydraulic fracture diagnostics in horizontal wells. LF-DAS in an untreated well provides far-field strain measurements, while offset wells are hydraulically fractured. Such a configuration is called crosswell LF-DAS sensing. Crosswell LF-DAS measurements have proved useful in diagnosing fracture hits, fracture azimuth, planarity, cluster efficiency, fracture propagation rates, and the dynamic distance to the fracture front. In contrast, in-well LF-DAS is conducted on the actively fractured well. Due to cool fracture fluid being injected at high injection rates, the strain component of the LF-DAS response is largely obscured by temperature changes. In permanent fiber-optic cable installations, distributed temperature sensing (DTS) is often conducted simultaneously with LF-DAS. An opportunity exists to decouple the temperature and strain components of LF-DAS sensors to observe strain changes on in-well LF-DAS. The LF-DAS response is modeled as linearly dependent on strain and temperature changes. Theoretical LF-DAS temperature and strain sensitivity coefficients are derived based on the changes to the index of refraction and length of the fiber. Using the DTS measurements, temperature changes are computed, smoothed, filtered, and compared to the LF-DAS response. Crossplots of the in-well LF-DAS measurements and temperature changes from DTS measurements far from the actively fractured region are used to validate the theoretical sensitivity coefficients. Uncertainty in the temperature component of the LF-DAS response is quantified. The difficulty in corresponding the different spatial and temporal resolutions of the DTS and LF-DAS measurements is overcome by comparing the responses over a moving temporal and spatial window. If the LF-DAS response at the center of the window agrees with the DTS response within uncertainty, the measurement is filtered out. After filtering, the remaining nonzero in-well LF-DAS measurements are due to changes in strain. The data are then visualized in waterfall plots. The results indicate that the theoretical and observed strain and temperature coefficients agree within 10%. After the temperature component of the in-well LF-DAS response is extracted, the remaining nonzero measurements are located primarily within the actively treated region. Locations with peaks in the strain response are interpreted to indicate fracture initiation points. These fracture initiation points are compared with in-well DTS and high-frequency DAS noise measurements across multiple stages to better understand fracture initiation along the horizontal well.

Engineering↗
DOE OSTI2023

Toward Improved Regional Hydrological Model Performance Using State-Of-The-Science Data-Informed Soil Parameters

Accurate soil moisture and streamflow data are an aspirational need of many hydrologically relevant fields. Model simulated soil moisture and streamflow hold promise but models require validation prior to application. Calibration methods are commonly used to improve model fidelity but misrepresentation of the true dynamics remains a challenge. In this study, we leverage soil parameter estimates from the Soil Survey Geographic (SSURGO) database and the probability mapping of SSURGO (POLARIS) to improve the representation of hydrologic processes in the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) over a central California domain. Our results show WRF-Hydro soil moisture exhibits increased correlation coefficients ( r ), reduced biases, and increased Kling-Gupta Efficiencies (KGEs) across seven in situ soil moisture observing stations after updating the model's soil parameters according to POLARIS. Compared to four well-established soil moisture data sets including Soil Moisture Active Passive data and three Phase 2 North American Land Data Assimilation System land surface models, our POLARIS-adjusted WRF-Hydro simulations produce the highest mean KGE (0.69) across the seven stations. More importantly, WRF-Hydro streamflow fidelity also increases, especially in the case where the model domain is set up with SSURGO-informed total soil thickness. The magnitude and timing of peak flow events are better captured, r increases across nine United States Geological Survey stream gages, and the mean KGE across seven of the nine gages increases from 0.12 to 0.66. Our pre-calibration parameter estimate approach, which is transferable to other spatially distributed hydrological models, can substantially improve a model's performance, helping reduce calibration efforts and computational costs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

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

A Multi-Instrument Cloud Condensation Nuclei Spectrum Product (Final Technical Report)

A wealth of observational data exists on the characteristics of atmospheric particulate matter, over multiple years, at the DOE ARM Southern Great Plains (SGP) Central Facility site. This site is located in a region of the country that frequently experiences weather extremes, and that is removed from many local sources of pollution but is affected by transported smoke, dust, and urban emissions. The relationships between particulate matter, cloud formation and evolution, and precipitation are therefore of strong interest, and are being explored via modeling on a variety of scales. These models require as input detailed information on the characteristics of particles capable of serving as the nuclei for cloud formation. Sufficient data exist to be able to put together a picture of the nature of the total aerosol and the cloud condensation nuclei (CCN) subset, and their variability, through merged data products. This study was aimed at exploiting the multiple measurement types at SGP to develop the first such multi-year estimates. Further, the resulting data were analyzed to understand temporal patterns ranging from hourly to seasonal, thereby gaining insights into the particle sources affecting the atmosphere in this region. DOE-funded datasets that were analyzed in this study include total particle number concentrations, submicron aerosol scattering coefficients, dry aerosol size distributions, and more recently, time-resolved submicron aerosol chemical composition. Data are also available for the number concentrations of particles that are activated in a cloud condensation nucleus instrument at a series of setpoint supersaturations, providing direct observations of the number concentrations of “CCN”. This variable is the quantity that is generally desired for inclusion in numerical models that seek to represent and predict the impacts of varying aerosol characteristics on the formation and microphysical properties of clouds. One limitation of the use of direct CCN observations is that they are not available for supersaturations larger than about 1%, which is insufficient for deep convection and may be insufficient even for shallow convection, depending on the nature of the available CCN and the dynamics of the cloud. We developed a data-based approach to representing the full aerosol size spectrum with size-dependent hygroscopicity, and used this to extrapolate CCN spectra beyond the limited measurements. Five years of SGP aerosol data (2009 -2013) were analyzed. As a side product of our work, we identified and communicated several previously-unflagged data quality issues. The resulting merged aerosol distributions, along with fits for seasonal averages, were published and submitted to the ARM archive as a special value-added product (VAP; submitted as a PI product). CCN spectra were computed for the same data period and will similarly be published and submitted to the archive for use by the community. We also note that our methodologies and findings have been discussed at several Joint ARM User Facility/Atmospheric System Research (ASR) Principal Investigators Meetings and that recent ARM/ASR aerosol data reporting strategies have included similar ideas for data merging, indicating that this work has had a lasting impact on ARM aerosol data acquisition and reporting. The proposed work advances the science of the interactions of aerosols, clouds and precipitation, with direct application to improve representation of such interactions for clouds in regional and global climate models. The archived data will continue to serve research studies in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Radionuclide-specific Parameters Dataset

The radionuclide-specific parameters dataset is searchable for radiological information for multiple isotopes simultaneously. After selecting radionuclides of interest and the desired parameters, the RAIS will generate a table containing the values, chosen according to an established hierarchy. Results can be downloaded in Excel format. 50 parameters are available, including atomic number, soil to animal transfer coefficients, plant uptake coefficients, half-life, specific activity, and water solubility. Seven primary sources are used to populate the dataset of radiological-specific parameters. These values should be used in cancer risk assessments for the calculation of preliminary remediation goals (PRGs), hazard characterization, and transport modeling. Users can select up to 1000 radionuclides per query. The dataset supports environmental risk assessments, regulatory decision-making, and environmental planning with tools for benchmarking against risk-based standards. This structured approach ensures a robust evaluation of environmental risks tailored to regulatory needs.

Manning, Karessa [Oak Ridge National Laboratory (O↗
DOE OSTI2020

Time reversal symmetry breaking and odd viscosity in active fluids: Green–Kubo and NEMD results

Active fluids, which are driven at the microscale by non-conservative forces, are known to exhibit novel transport phenomena due to the breaking of time reversal symmetry. Recently, Epstein and Mandadapu [arXiv:1907.10041 (2019)] obtained Green-Kubo relations for the full set of viscous coefficients governing isotropic chiral active fluids, including the so-called odd viscosity, invoking Onsager's regression hypothesis for the decay of fluctuations in active non-equilibrium steady states. In this Communication, we test these Green-Kubo relations using molecular dynamics simulations of a canonical model system consisting of actively torqued dumbbells. We find the resulting odd and shear viscosity values from the Green-Kubo relations to be in good agreement with values measured independently through non-equilibrium molecular dynamics flow simulations. This provides a test of the Green-Kubo relations and lends support to the application of the Onsager regression hypothesis in relation to viscous behaviors of active matter systems.

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

RAIS Slope Factors for Radionuclides

For a given radionuclide, Slope Factors (SFs) represent the ELCR equivalent per unit intake (i.e., ingestion or inhalation) or external exposure of that radionuclide. These SFs are used to convert a radionuclide concentration in soil, air, water, or foodstuffs to a radiation ELCR. The primary use of SFs, also called risk coefficients, is to compute the ELCR resulting from site-related exposures. This is accomplished by multiplying the route-specific SF by the CDI of each radionuclide of potential concern for each route of exposure. Ingestion and inhalation slope factors are central estimates in a linear model of the age-averaged, lifetime attributable radiation cancer incidence (fatal and nonfatal cancer) risk per unit of activity inhaled or ingested, expressed as risk/picocurie (pCi). External exposure SFs are central estimates of lifetime attributable radiation cancer incidence risk for each year of exposure to external radiation from photon-emitting radionuclides distributed uniformly in a thick layer of soil, expressed as risk/year per pCi/gram soil. External exposure SF units can also be risk/year per pCi/square centimeters (cm2) of soil. When combined with site-specific media concentration data and appropriate exposure assumptions, SFs can be used to estimate lifetime cancer risks to members of the general population due to radionuclide exposures. The SFs presented on are updated values from Federal Guidance Report (FGR) 13 supplement using the International Commission on Radiological Protection (ICRP) 107 decay data. The derivation and values are found in “Calculations of Slope Factors and Dose Coefficients” (ORNL, 2014).

Noto, Katie [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2024

Near wake evolution of a tidal stream turbine due to asymmetric sheared turbulent inflow with different integral length scales

Tidal stream turbines deployed at highly energetic open water sites are subjected to sheared inflow in the rotor plane. The inflow shear is expected to cause asymmetric loading on the rotor blades and affect the downstream wake. In the current study, two different turbulent inflow conditions, static-high shear and dynamic shear, were generated via an active-grid turbulence generator. A 1:20 scaled three-bladed horizontal axis tidal turbine model was tested in those conditions. The results were compared to a quasi-laminar case with no imposed turbulence or shear. The results show that the high shear reduces the average performance, with a drop of up to 16% in the optimal power coefficient. Besides, the shear profiles increase torque fluctuations and induce significant differences in wake hydrodynamics between the high-speed (upper) and low-speed (lower) regions. The large integral length scales further enhance the load fluctuations perceived by the rotor but have a negligible effect on the mean wake field quantities and the wake recovery. The lower half region featured a faster breakdown of tip vortex structure and a rapid drop of swirl number, a phenomenon conjectured to be a consequence of the strong turbulence intensities and Reynolds stresses in the lower half region. Furthermore, the sheared turbulent inflow also results in a very intensive energy redistribution process towards large-scale, low-frequency motions, which is important to the downstream turbines.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Viscous flow properties and hydrodynamic diameter of phenothiazine-based redox-active molecules in different supporting salt environments

We report viscous flow properties of a redox-active organic molecule, N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT), a candidate for non-aqueous redox flow batteries, and two of its radical cation salts. A microfluidic viscometer enabled the use of small sample volumes in determining viscosity as a function of shear rate and concentration in the non-aqueous solvent, acetonitrile, both with and without supporting salts. All solutions tested show Newtonian behavior over shear rates of up to 30 000 s−1, which was rationalized by scaling arguments for the diffusion-based relaxation time of a single MEEPT molecule without aggregation. Neat MEEPT is flowable but with a large viscosity (412 mPa⋅s at room temperature), which is ∼1000 times larger than that of acetonitrile. MEEPT solutions in acetonitrile have low viscosities; at concentrations up to 0.5 M, the viscosity increases by less than a factor of two. From concentration-dependent viscosity measurements, molecular information was inferred from intrinsic viscosity (hydrodynamic diameter) and the Huggins coefficient (interactions). Model fit credibility was assessed using the Bayesian Information Criterion. It is found that the MEEPT and its charged cations are “flowable” and do not flocculate at concentrations up to 0.5 M. MEEPT has a hydrodynamic diameter of around 8.5 Å, which is almost insensitive to supporting salt and state of charge. This size is comparable to molecular dimensions of single molecules obtained from optimized structures using density functional theory calculations. The results suggest that MEEPT is a promising candidate for redox flow batteries in terms of its viscous flow properties.

Wang, Yilin (ORCID:0000000316102995)↗
DOE OSTI2021

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

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

Determination of Infinite Dilution Activity Coefficients of Molecular Solutes in Ionic Liquids and Deep Eutectic Solvents by Factorization-Machine-Based Neural Networks

Widely known as “green solvents,” ionic liquids (ILs) and deep eutectic solvents (DESs) have been used as substitutes for traditional organic solvents in separation science. To achieve better separations using ILs and DESs, this work aimed to predict the infinite dilution activity coefficients (IDACs) of molecular solutes in these solvents with the state-of-the-art factorization-machine-based neural network (DeepFM). DeepFM combines the benefits of factorization machines and deep neural networks to learn both low-order and high-order interactions among features. The IDAC prediction model was established with 52,372 experimental IDAC datapoints including 260 solvents (252 ILs and 8 DESs) and 112 molecular solutes collected at various temperatures from 288.15 to 428.15 K. Chemical information describing the ILs, DESs, and molecular solutes was included in the IDAC prediction model, including chemical functional groups, molecular weights, and Abraham solvation parameters. The IDAC prediction model showed an improved accuracy compared with alternative models; additionally, DESs were included in the IDAC prediction model for the first time. Here, the model will reduce the energy and resources needed to optimize the selection of ILs and DESs for specific separations, which will promote the development of these green solvents for sustainable chemical processes.

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

One-step lignocellulose depolymerization and saccharification to high sugar yield and less condensed isolated lignin

The cost of sugar production remains a key challenge in future lignocellulosic biorefineries. We demonstrate that ZnBr 2 , an inexpensive inorganic salt, provides nearly theoretical yields of glucose and xylose in one-step from poplar wood at 85 °C and short reaction times at molten salt hydrate (MSH) conditions without an acid. Catalytic depolymerization of the isolated MSH lignin, using a CoS 2 catalyst, yields 17% phenol-like monomers compared to only 1% produced from the acidified MSH lignin. Reductive catalytic fractionation of MSH lignin over Ru/C resulted in two times higher total monomer yield compared to the AMSH lignin. Both the lignin samples were characterized using 2D HSQC NMR and the thioacidolysis method. Thioacidolysis studies reveal 8.4% and 1.8% of β-O-4 linkages in MSH and acidified MSH lignin, respectively. Thermodynamic modeling and 13 C NMR spectroscopy indicate that the effectiveness of this catalyst arises from the strong interaction of the Lewis acidic zinc cation (Zn 2+ ) with the coordinated water molecules resulting in hydrolysis of the metal aquo complex and to the salt-driven increase in the H + activity coefficient. Techno-economic analysis demonstrates that despite being slower, the ZnBr 2 MSH media has cost advantages, compared to conventional hydrolysis and even to the LiBr and ZnBr 2 AMSH, due to the higher quality of lignin.

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

Evaluating the limitations of Bayesian metabolic control analysis

Bayesian Metabolic Control Analysis (BMCA) is a promising framework for inferring metabolic control coefficients in data-limited scenarios, combining Bayesian inference with linear-logarithmic (lin-log) rate laws. These metabolic control coefficients quantify how changes in enzyme activities affect steady-state fluxes and metabolite concentrations across a metabolic network. However, its predictive accuracy and limitations remain underexplored. This study systematically evaluates BMCA’s ability to infer elasticity values, flux control coefficients (FCC), and concentration control coefficients (CCC) under varying data availability conditions using three synthetic metabolic network models. We demonstrate that BMCA predictions are highly dependent on the inclusion of flux and enzyme concentration data, with the omission of these datasets leading to severe inaccuracies. In our synthetic, enzyme-perturbation datasets, external metabolite concentrations had minimal impact and, in some cases, their exclusion improved predictions; when external-nutrient perturbations were introduced and those concentrations were observed, gains were at most modest. Additionally, we find that posterior estimation with both ADVI and HMC can underestimate large-magnitude elasticities in our synthetic settings, with ADVI showing somewhat higher variance under strong up-regulation; thus, recovering |elasticity| ≳ 1.5 remains challenging regardless of the inference engine. ADVI also fails to accurately infer allosteric interactions, even when regulatory effects are strong. While BMCA maintains reasonable accuracy in partially recovering the rankings of the highest FCC values, its estimates of absolute values remain constrained by prior assumptions and data limitations. Our findings reveal the BMCA algorithm’s strengths and weaknesses, providing guidance on its application in metabolic engineering, and highlighting the need for methodological refinements to enhance its predictive capabilities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Modeling and Analysis of Polarization Losses in Solid Oxide Fuel Cells with Siloxane Contamination

In this study, the degradation of the solid-oxide fuel cell (SOFC) nickel-yttria stabilized zirconia anode under decamethyltetrasiloxane (L4) contamination is examined with experiments and modeling. A model is developed for the polarization losses based on the charge transfer coefficient,α, and diffusion layer thickness,δ, and fitted to the experimental data to understand how the siloxane degrades the SOFC performance with time. The results of the model indicate that the total polarization losses increase approximately 44% over the course of the 180 min experiment at 350 mA cm −2 . Activation losses dominate the polarization losses initially but decrease in their total contribution while concentration losses increase. Scanning electron microscopy (SEM) with wavelength dispersive X-ray spectroscopy (WDS) elemental mapping indicates that silicon deposition is highest at the outer edge of the anode and forms a barrier layer to fuel diffusion, increasing concentration losses. When the model is applied to other previous D4 and L4 siloxane experiments conducted over a period of 40 h, similar trends in polarization losses are observed. Polarization losses increase more rapidly with D4 compared to L4 siloxane contamination, with concentration losses increasing the fastest with both types of siloxane.

Electrochemistry↗