Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERMEDIATE COUPLING”

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

Proton removal from Br 73 , 75 to Se 72 , 74 at intermediate energies

Herein we report new experimental data for excited states of 72,74 Se obtained from proton removal from 73,75 Br secondary beams on a proton target. The experiments were performed with the Ursinus-NSCL Liquid Hydrogen Target and the combined GRETINA+S800 setup at the Coupled Cyclotron Facility of the National Superconducting Cyclotron Laboratory at Michigan State University. Within uncertainties, the inclusive cross sections for proton removal from 73,75 Br on a proton target are identical suggesting that the same single-particle orbitals contribute to the proton-removal reaction. In addition, details of the partial cross section fragmentation are discussed. The data might suggest that l = 1, 2, 3, and 4 angular momentum transfers are important to understand the population of excited states of 72,74 Se in proton removal. Available data for excited states of 74 Ge populated through the 75 As(d, 3 He) 74 Ge proton-removal reaction in normal kinematics suggest indeed that the fp and sd shell as well as the 1g 9/2 orbital contribute. A comparison to data available for odd-A nuclei supports that the bulk of the spectroscopic strengths could be found at lower energies in the even-even Se isotopes than in, for instance, the even-even Ge isotopes. In addition, the population of high-J states seems to indicate that multi-step processes contribute to proton-removal reactions at intermediate energies in these collective nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

End-of-Century Changes in Orographic Precipitation with the Intermediate Complexity Atmospheric Research Model over the Western United States

Abstract Downscaled precipitation projections were created using the Intermediate Complexity Atmospheric Research (ICAR) model over the western United States to increase the physical realism in orographic precipitation changes. End-of-century simulations from eight models in phase 5 of the Coupled Model Intercomparison Project (CMIP5) were downscaled with ICAR and compared to the widely utilized statistically downscaled dataset, localized constructed analogs (LOCAs), to understand where and why projections of cool-season (September–May) precipitation differed. ICAR and LOCA precipitation projections were similar, but their sign differed in hydrologically relevant regions likely due to ICAR’s simulation of microphysics and mesoscale dynamics with high-resolution topography (6 km). In the Pacific Northwest, cool-season precipitation projections from ICAR showed an increase on the windward side of the Cascades and no significant change within the lee. This difference between the windward and leeward side was attributed to reduced zonal wind speeds, allowing more time for microphysical processes within ICAR. This contrast is enhanced by rain’s faster fall speed compared to snow, limiting transport into the lee. Meanwhile, LOCA projected an increase in precipitation across the Cascades. In the Upper Colorado River basin, LOCA projected an increase in precipitation in high elevation regions (>3000 m), but ICAR projected no significant change or a decrease in precipitation. High elevation differences were most evident in the spring and fall and were also attributed to a snow-to-rain transition and dynamical processes that impacted orographic enhancement within ICAR. Idealized, controlled studies are needed to better isolate individual processes, but these results underscore the importance of including microphysics and mesoscale dynamics within regional-scale precipitation projections. Significance Statement A set of global climate model simulations was downscaled using an atmospheric model that contains key physical equations, referred to as Intermediate Complexity Atmospheric Research (ICAR). ICAR was used to examine projected changes in end-of-century cool-season precipitation over mountains in the western United States. Precipitation projections from ICAR were similar to projections that used statistical relationships to downscale climate projections. However, projections differed between ICAR and statistically downscaled datasets in whether they increased, decreased, or stayed the same in specific, hydrologically relevant regions such as the eastern Cascades and high elevation areas of the Upper Colorado River basin. These differences were attributed to the simulation of physical processes in ICAR. The results highlight the importance of kilometer-scale atmospheric processes in regional climate projections.

Currier, William Ryan [NOAA/Physical Sciences Labo↗
DOE OSTI2021

Test of the universality of τ and μ lepton couplings in W-boson decays with the ATLAS detector

The standard model of particle physics encapsulates our best current understanding of physics at the smallest scales. A fundamental axiom of this theory is the universality of the couplings of the different generations of leptons to the electroweak gauge bosons. The measurement of the ratio of the decay rate of W bosons to τ leptons and muons, R ( τ/μ ), constitutes an important test of this axiom. Using 139 fb –1 of proton–proton collisions recorded with the ATLAS detector at a centre-of-mass energy of 13 TeV, we report a measurement of this quantity from di-leptonic $t\overline{t}$ events where the top quarks decay into a W boson and a bottom quark. We can distinguish muons originating from W bosons and those originating from an intermediate τ lepton through the muon transverse impact parameter and differences in the muon transverse momentum spectra. The measured value of R ( τ/μ ) is 0.992 ± 0.013 [± 0.007(stat) ± 0.011(syst)] and is in agreement with the hypothesis of universal lepton couplings as postulated in the standard model. This is the only such measurement from the Large Hadron Collider, so far, and obtains twice the precision of previous measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2004

Sensitivity of Precipitation in Coupled Land-Atmosphere Models

The project objective was to understand mechanisms by which atmosphere-land-ocean processes impact precipitation in the mean climate and interannual variations, focusing on tropical and subtropical regions. A combination of modeling tools was used: an intermediate complexity land-atmosphere model developed at UCLA known as the QTCM and the NASA Seasonal-to-Interannual Prediction Program general circulation model (NSIPP GCM). The intermediate complexity model was used to develop hypotheses regarding the physical mechanisms and theory for the interplay of large-scale dynamics, convective heating, cloud radiative effects and land surface feedbacks. The theoretical developments were to be confronted with diagnostics from the more complex GCM to validate or modify the theory.

Neelin, David↗
DOE OSTI2016

High resolution spectroscopy of jet cooled phenyl radical: The ν{sub 1} and ν{sub 2} a{sub 1} symmetry C–H stretching modes

A series of CH stretch modes in phenyl radical (C{sub 6}H{sub 5}) has been investigated via high resolution infrared spectroscopy at sub-Doppler resolution (∼60 MHz) in a supersonic discharge slit jet expansion. Two fundamental vibrations of a{sub 1} symmetry, ν{sub 1} and ν{sub 2}, are observed and rotationally analyzed for the first time, corresponding to in-phase and out-of-phase symmetric CH stretch excitation at the ortho/meta/para and ortho/para C atoms with respect to the radical center. The ν{sub 1} and ν{sub 2} band origins are determined to be 3073.968 50(8) cm{sup −1} and 3062.264 80(7) cm{sup −1}, respectively, which both agree within 5 cm{sup −1} with theoretical anharmonic scaling predictions based on density functional B3LYP/6-311g++(3df,3dp) calculations. Integrated band strengths for each of the CH stretch bands are analyzed, with the relative intensities agreeing remarkably well with theoretical predictions. Frequency comparison with previous low resolution Ar-matrix spectroscopy [A. V. Friderichsen et al., J. Am. Chem. Soc. 123, 1977 (2001)] reveals a nearly uniform Δν ≈ + 10-12 cm{sup −1} blue shift between gas phase and Ar matrix values for ν{sub 1} and ν{sub 2}. This differs substantially from the much smaller red shift (Δν ≈ − 1 cm{sup −1}) reported for the ν{sub 19} mode, and suggests a simple physical model in terms of vibrational mode symmetry and crowding due to the matrix environment. Finally, the infrared phenyl spectra are well described by a simple asymmetric rigid rotor Hamiltonian and show no evidence for spectral congestion due to intramolecular vibrational coupling, which bodes well for high resolution studies of other ring radicals and polycyclic aromatic hydrocarbons. In summary, the combination of slit jet discharge methods with high resolution infrared lasers enables spectroscopic investigation of even highly reactive combustion and interstellar radical intermediates under gas phase, jet-cooled (T{sub rot} ≈ 11 K) conditions.

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

Unconventional Pathway in the Gas‐Phase Synthesis of 9 H ‐Fluorene (C 13 H 10 ) via the Radical–Radical Reaction of Benzyl (C 7 H 7 ) with Phenyl (C 6 H 5 )

Abstract The simplest polycyclic aromatic hydrocarbon (PAH) carrying a five‐membered ring—9 H ‐fluorene (C 13 H 10 )—is produced isomer‐specifically in the gas phase by reacting benzyl (C 7 H 7 ⋅) with phenyl (C 6 H 5 ⋅) radicals in a pyrolytic reactor coupled with single photon ionization mass spectrometry. The unconventional mechanism of reaction is supported by theoretical calculations, which first produces diphenylmethane and unexpected 1‐(6‐methylenecyclohexa‐2,4‐dienyl)benzene intermediates (C 13 H 12 ) accessed via addition of the phenyl radical to the ortho position of the benzyl radical. These findings offer convincing evidence for molecular mass growth processes defying conventional wisdom that radical‐radical reactions are initiated through recombination at their radical centers. The structure of 9 H ‐fluorene acts as a molecular building block for complex curved nanostructures like fullerenes and nanobowls providing fundamental insights into the hydrocarbon evolution in high temperature settings.

He, Chao↗
DOE OSTI2022

Unconventional Pathway in the Gas-Phase Synthesis of 9 H -Fluorene (C 13 H 10 ) via the Radical–Radical Reaction of Benzyl (C 7 H 7 ) with Phenyl (C 6 H 5 )

The simplest polycyclic aromatic hydrocarbon (PAH) carrying a five-membered ring - 9H-fluorene (C 13 H 10 ) - is produced isomer-specifically in the gas phase by reacting benzyl (C 7 H 7 •) with phenyl (C 6 H 5 •) radicals in a pyrolytic reactor coupled with single photon ionization mass spectrometry. The unconventional mechanism of reaction is supported by theoretical calculations, which first produces diphenylmethane and unexpectedly 1-(6-methylenecyclohexa-2,4-dienyl)benzene intermediates (C 13 H 12 ) accessed via addition of the phenyl radical to the ortho position of the benzyl radical. These findings offer convincing evidence for molecular mass growth processes defying conventional wisdom that radical–radical reactions are initiated through recombination at their radical centers. Furthermore, the structure of 9H-fluorene acts as a molecular building block for complex curved nanostructures like fullerenes and nanobowls providing fundamental insights into the hydrocarbon evolution in high temperature settings.

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

{sup 45}Sc Solid State NMR studies of the silicides ScTSi (T=Co, Ni, Cu, Ru, Rh, Pd, Ir, Pt)

The silicides ScTSi (T=Fe, Co, Ni, Cu, Ru, Rh, Pd, Ir, Pt) were synthesized by arc-melting and characterized by X-ray powder diffraction. The structures of ScCoSi, ScRuSi, ScPdSi, and ScIrSi were refined from single crystal diffractometer data. These silicides crystallize with the TiNiSi type, space group Pnma. No systematic influences of the {sup 45}Sc isotropic magnetic shift and nuclear electric quadrupolar coupling parameters on various structural distortion parameters calculated from the crystal structure data can be detected. {sup 45}Sc MAS-NMR data suggest systematic trends in the local electronic structure probed by the scandium atoms: both the electric field gradients and the isotropic magnetic shifts relative to a 0.2 M aqueous Sc(NO{sub 3}){sub 3} solution decrease with increasing valence electron concentration and within each T group the isotropic magnetic shift decreases monotonically with increasing atomic number. The {sup 45}Sc nuclear electric quadrupolar coupling constants are generally well reproduced by quantum mechanical electric field gradient calculations using the WIEN2k code. Highlights: Black-Right-Pointing-Pointer Arc-melting synthesis of silicides ScTSi. Black-Right-Pointing-Pointer Single crystal X-ray data of ScCoSi, ScRuSi, ScPdSi, and ScIrSi. Black-Right-Pointing-Pointer {sup 45}Sc solid state NMR of silicides ScTSi.

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

Universal molecular-confined synthesis of interconnected porous metal oxides-N-C frameworks for electrocatalytic water splitting

The rational synthesis of high performance electrocatalysts at low cost for water splitting (hydrogen and oxygen evolution reaction, HER and OER) is highly desirable but remains a key challenge. Herein, we report a versatile molecular-confining route to construct strongly coupled metal oxides-N-C frameworks with interconnected configuration. By simply chelating various transition metal ions with ethylenediaminetetraacetic acid disodium salt (EDTA) in agarose hydrogel and subsequent carbonization, the high specific surface area porous N-C frameworks with strongly coupled metal oxides (e.g., manganese, iron, cobalt, nickel and their mixed oxides) can be prepared as advanced nonprecious water splitting electrocatalysts. For instance, the resulting Co 3 O 4 -N-C frameworks with the high surface area of 153 m 2 g -1 , and the moderate nitrogen content of 1.23%, require an overpotential of 324 mV to afford a current density of 10 mA cm -2 in 0.1 M KOH for OER, superior to commercial RuO 2 catalysts. Density functional theory (DFT) calculations reveal that the strong coupling between N-C and Co 3 O 4 tunes the local electronic structure of Co for high-valence active sites and assures optimal adsorption energies of OER intermediates. Moreover, Co 2 P-N-C electrocatalysts derived from Co 3 O 4 -N-C also exhibit excellent HER performance under 1.0 M KOH with a low overpotential of 139 mV to reach 10 mA cm -2 , a small Tafel slope of 45 mV dec -1 and impressive stability, underscoring the versatility of our synthetic strategy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Extraordinary Magnetic Response of an Anisotropic 2D Antiferromagnet via Site Dilution

A prominent characteristic of 2D magnetic systems is the enhanced spin fluctuations, which reduce the ordering temperature. We report that a magnetic field of only 1000th of the Heisenberg superexchange interaction can induce a crossover, which for practical purposes is the effective ordering transition, at temperatures about 6 times the Néel transition in a site-diluted two-dimensional anisotropic quantum antiferromagnet. Such a strong magnetic response is enabled because the system directly enters the antiferromagnetically ordered state from the isotropic disordered state, skipping the intermediate anisotropic stage. The underlying mechanism is achieved on a pseudospin-half square lattice realized in the [(SrIrO3)1/(SrTiO3)2] superlattice thin film that is designed to linearly couple the staggered magnetization to external magnetic fields by virtue of the rotational symmetry-preserving Dzyaloshinskii–Moriya interaction. Furthermore, our model analysis shows that the skipping of the anisotropic regime despite finite anisotropy is due to the enhanced isotropic fluctuations under moderate dilution.

2D antiferromagnet↗
DOE OSTI2014

Near-resonant rotational energy transfer in HCl–H{sub 2} inelastic collisions

We present a new four-dimensional (4D) potential energy surface for the HCl–H{sub 2} van der Waals system. Both molecules were treated as rigid rotors. Potential energy surface was obtained from electronic structure calculations using a coupled cluster with single, double, and perturbative triple excitations method. The four atoms were described using the augmented correlation-consistent quadruple zeta basis set and bond functions were placed at mid-distance between the HCl and H{sub 2} centers of mass for a better description of the van der Waals interaction. The global minimum is characterized by the well depth of 213.38 cm{sup −1} corresponding to the T-shape structure with H{sub 2} molecule on the H side of the HCl molecule. The dissociation energies D{sub 0} are 34.7 cm{sup −1} and 42.3 cm{sup −1} for the complex with para- and ortho-H{sub 2}, respectively. These theoretical results obtained using our new PES are in good agreement with experimental values [D. T. Anderson, M. Schuder, and D. J. Nesbitt, Chem. Phys. 239, 253 (1998)]. Close coupling calculations of the inelastic integral rotational cross sections of HCl in collisions with para-H{sub 2} and ortho-H{sub 2} were performed at low and intermediate collisional energies. Significant differences exist between para- and ortho-H{sub 2} results. The strongest collision-induced rotational HCl transitions are the transitions with Δj = 1 for collisions with both para-H{sub 2} and ortho-H{sub 2}. Rotational relaxation of HCl in collision with para-H{sub 2} in the rotationally excited states j = 2 is dominated by near-resonant energy transfer.

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

Anomalous magnetic structure and spin dynamics in magnetoelectric LiFePO 4

We report significant details of the magnetic structure and spin dynamics of LiFePO 4 obtained by single-crystal neutron scattering. Our results confirm a previously reported collinear rotation of the spins away from the principal b axis, and they determine that the rotation is toward the a axis. In addition, we find a significant spin-canting component along c. Furthermore, the possible causes of these components are discussed, and their significance for the magnetoelectric effect is analyzed. Inelastic neutron scattering along the three principal directions reveals a highly anisotropic hard plane consistent with earlier susceptibility measurements. While using a spin Hamiltonian, we show that the spin dimensionality is intermediate between XY- and Ising-like, with an easy b axis and a hard c axis. As a result, it is shown that both next-nearest neighbor exchange couplings in the bc plane are in competition with the strongest nearest neighbor coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Deciphering competing elementary steps to correlate electrocatalyst chemical state with activity

The overpotential in multielectron transfer heterogeneous electrocatalysis fundamentally arises from thermodynamic and kinetic disparities among elementary steps; however, deciphering coupled and competing steps has long remained a challenge. Here, we establish an electrochemical deconvolution paradigm based on key processes in electrocatalytic reactions, such as charge accumulation, electron/proton transfer, and intermediate evolution, to resolve competing elementary steps. Taking the oxygen evolution reaction as a prototypical reaction, we design a model catalyst featuring a precise isolated cation-anion vacancy pair and track the previously elusive electrochemical behavior of lattice oxygen by disentangling interference from adsorbed oxygen intermediates. Mechanistically, the lattice oxygen oxidation pathway originates from the spontaneous, nonelectrochemical deprotonation of replenished water molecules coordinated to unsaturated cation sites. Alternating current techniques further reveal that although lattice oxygen oxidation requires a higher potential than metal oxidation, it exhibits faster kinetics, providing insight into its superior catalytic activity. These findings establish a direct experimental correlation between the initial chemical state and the catalytic activity and prove surface-confined lattice oxygen cycling. Furthermore, expanding conventional potential-current analysis into a multidimensional framework enables disentanglement of thermodynamic and kinetic contributions of key elementary steps, thereby guiding the rational optimization of various complex multielectron transfer reactions.

OER↗
DOE OSTI2023

The genetics of aerotolerant growth in an alphaproteobacterium with a naturally reduced genome

Reduced genome bacteria are genetically simplified systems that facilitate biological study and industrial use. The free-living alphaproteobacterium Zymomonas mobilis has a naturally reduced genome containing fewer than 2,000 protein-coding genes. Despite its small genome, Z. mobilis thrives in diverse conditions including the presence or absence of atmospheric oxygen. However, insufficient characterization of essential and conditionally essential genes has limited broader adoption of Z. mobilis as a model alphaproteobacterium. Here, we use genome-scale CRISPRi-seq (clustered regularly interspaced short palindromic repeats interference sequencing) to systematically identify and characterize Z. mobilis genes that are conditionally essential for aerotolerant or anaerobic growth or are generally essential across both conditions. Comparative genomics revealed that the essentiality of most “generally essential” genes was shared between Z. mobilis and other Alphaproteobacteria, validating Z. mobilis as a reduced genome model. Among conditionally essential genes, we found that the DNA repair gene, recJ, was critical only for aerobic growth but reduced the mutation rate under both conditions. Further, we show that genes encoding the F1FO ATP synthase and Rhodobacter nitrogen fixation (Rnf) respiratory complex are required for the anaerobic growth of Z. mobilis. Combining CRISPRi partial knockdowns with metabolomics and membrane potential measurements, we determined that the ATP synthase generates membrane potential that is consumed by Rnf to power downstream processes. Rnf knockdown strains accumulated isoprenoid biosynthesis intermediates, suggesting a key role for Rnf in powering essential biosynthetic reactions. Our work establishes Z. mobilis as a streamlined model for alphaproteobacterial genetics, has broad implications in bacterial energy coupling, and informs Z. mobilis genome manipulation for optimized production of valuable isoprenoid-based bioproducts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The Genetics of Aerotolerant Growth in Zymomonas mobilis

Reduced genome bacteria are genetically simplified systems that facilitate biological study and industrial use. The free-living Alphaproteobacterium, Zymomonas mobilis, has a naturally reduced genome containing fewer than 2000 protein coding genes. Despite its small genome, Z. mobilis thrives in diverse conditions including the presence or absence of atmospheric oxygen. However, insufficient characterization of essential and conditionally essential genes has limited broader adoption of Z. mobilis as a model Alphaproteobacterium. Here, we use genome-scale CRISPRi-seq to systematically identify and characterize Z. mobilis genes that are conditionally essential for aerotolerant or anaerobic growth, or are generally essential across both conditions. Comparative genomics revealed that the essentiality of most "generally essential" genes was shared between Z. mobilis and other Alphaproteobacteria, validating Z. mobilis as reduced genome model. Among conditionally essential genes, we found that the DNA repair gene, recJ, was critical only for aerobic growth but reduced the mutation rate under both conditions. Further, we show that genes encoding the F1FO ATP synthase and Rnf respiratory complex are required for anaerobic growth of Z. mobilis. Combining CRISPRi partial knockdowns with metabolomics and membrane potential measurements, we determined that the ATP synthase generates membrane potential that is consumed by Rnf to power downstream processes. Rnf knockdown strains accumulated isoprenoid biosynthesis intermediates, suggesting a key role for Rnf in powering essential biosynthetic reactions. Our work establishes Z. mobilis as a streamlined model for alphaproteobacterial genetics, has broad implications in bacterial energy coupling, and informs Z. mobilis genome manipulation for optimized production of valuable isoprenoid-based bioproducts

Alphaproteobacteria↗
NASA NTRS1985

Cyanobacterial mats: Microanalysis of community metabolism

The microbial communities in two sites were studied using several approaches: (1) light microscopy; (2) the measurement of microprofiles of oxygen and sulfide at the surface of the microbial mat; (3) the study of diurnal variation of oxygen and sulfides; (4) in situ measurement of photosynthesis and sulfate reduction and study of the coupling of these two processes; (5) measurement of glutathione in the upper layers of the microbial mat as a possible oxygen quencher; (6) measurement of reduced iron as a possible intermediate electron donor along the established redoxcline in the mats; (7) measurement of dissolved phosphate as an indicator of processes of break down of organic matter in these systems; and (8) measurement of carbon dioxide in the interstitial water and its delta C-13 in an attempt to understand the flow of CO2 through the systems. Microbial processes of primary production and initial degradation at the most active zone of the microbial mat were analyzed.

Cohen, Y.↗
NASA NTRS1994

Gravitational instabilities in a proto-planetary disk including the effects of magnetic fields

We investigate the gravitational instability of a thin, Keplerian protoplanetary disk including the effects of a largely azimuthal magnetic field. The model follows that of our previous work (Noh, Vishniac, & Cochran 1991) except for the inclusion of a magnetic field. The disk is assumed to consist of neutral and ionized gas and neutral dust which are coupled by gravity and friction. The growth rates and eigenfunctions are calculated numerically using nonaxisymmetric linear perturbation methods. The results show that the growth rate has a maximum at some intermediate azimuthal number m, but for each value of m it is reduced relative to the unmagnetized case. The effects of the magnetic field appear more strongly on small scales. As the strength of the equilibrium magnetic field increases the growth rates decrease, and the maximum instability occurs at a lower value on m due to the increasing magnetic pressure. The response of each component to the magnetic field is discussed using the behavior of the eigenfunctions in the radial direction. With the inclusion of the magnetic field, the effects of the ionization fraction and friction on the growth rates also appear to be important for high m modes. Increasing the ionization fraction or the friction suppresses instability, but only slightly changes the maximally unstable azimuthal scales. The enhanced growth rates due to a dust component for which thermal pressure is negligible are somewhat reduced by the inclusion of a magnetic field. The effects of different boundary conditions (reflecting and transmitting) on the growth rates are also shown.

Noh, Hyerim↗
DOE OSTI2021

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗