Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Critical Fluids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Thermal Protection System Evaluation Using Arc-jet Flows: Flight Simulation or Research Tool?

The arc-jet has been used to evaluate thermal protection systems (TPS) and materials for the past forty years. Systems that have been studied in this environmerd include ablators, active, and passive TPS concepts designed for vehicles entering planetary and Earth atmospheres. The question of whether arc-jet flow can simulate a flight environment or is it a research tool that provides an aero-thermodynamic heating environment to obtain critical material properties will be addressed. Stagnation point tests in arc-jets are commonly used to obtain material properties such as mass loss rates, thermal chemical stability data, optical properties, and surface catalytic efficiency. These properties are required in computational fluid dynamic codes to accurately predict the performance of a TPS during flight. Special facilities have been developed at NASA Ames Research Center to approximate the flow environment over the mid-fuselage and body flap regions of proposed space-planes type vehicles. This paper compares flow environments generated in flight over a vehicle with those created over an arc-jet test articles in terms of scale, chemistry, and fluid dynamic properties. Flight experiments are essential in order to validate the material properties obtained from arc-jet tests and used to predict flight performance of any TPS being considered for use on a vehicle entering the Earth atmosphere at hypersonic speed.

Stewart, David A.↗
NASA NTRSDate not supplied

Role of Cryogenic Aerodynamic Testing for Current and Future Technologies

The motivation for cryogenic wind tunnels originated during the 1960s from an internationally recognized need for a high Reynolds number test capability based on experiences with preflight predictions of aerodynamic characteristics and an anticipated need in support of research and development for future aerospace vehicle systems. More specifically, the motivation for a flight (or near flight) Reynolds number ground test facility was captured in the foreword to reference 1 which states:“...AGARD held a Specialists’ Meeting in Paris on ‘Transonic Aerodynamics’ in recognition of the fact that the absence of adequate theoretical methods and wind tunnels of high enough Reynolds number had already led to costly shortcomings in the transonic performance of certain combat and transport aircraft.” Preflight prediction of flight characteristics is a necessary process for the developer of any aerospace vehicle, and introduces significant risk to the success of the vehicle. Whether the vehicle customer is commercial or governmental, the final full-scale vehicle must meet certain requirements to be certified as safe, and certain performance requirements to be economically successful. Aircraft companies strive to know the flight characteristics and performance of their vehicle with high confidence prior to flight, thus enabling optimal design trades and elimination of any costly modifications to the aircraft during and after initial flight testing. The problems of predicting flight characteristics across the full flight envelope prior to flight have been and continue to be challenging. Much has changed since the US National Transonic Facility (NTF) became operational in the 1980s, followed by the European Transonic Windtunnel (ETW) in 1990s, not the least of which has been the advancement and positive impact of computational fluid dynamics (CFD) on aircraft design and development. Today, the need for integrated CFD and ground-based high Reynolds number test and evaluation capability remains as aerodynamics will always be central to defining an aircraft directly as well as providing critical input to other disciplines. Today the economic stakes of being surprised during flight tests are higher than ever. This oral-only presentation describes the past, current, and future role of high Reynolds number aerodynamic testing in aerospace vehicle design and development. Discussion of relevant flow physics, past experiences with preflight prediction, and future needs, challenges, and opportunities is included. 1.“Facilities and Techniques for Aerodynamic Testing at Transonic Speeds and High Reynolds Number,” AGARD CP-83-71, 1971.

Aerodynamics↗
DOE OSTI2022

Boiling and cavitation caused by transient heat transfer in superfluid helium-4

Superfluid helium-4 (He II) has been widely utilized as a coolant in various scientific and engineering applications due to its superior heat transfer capability. An important parameter required in the design of many He II based cooling systems is the peak heat flux $q^*$, which refers to the threshold heat flux above which boiling spontaneously occurs in He II. Past experimental and numerical studies showed that $q^*$ increases when the heating time $t_h$ is reduced, which leads to an intuitive expectation that very high $q^*$ may be achievable at sufficiently small $t_h$. Knowledge on how $q^*$ actually behaves at small $t_h$ is important for applications such as laser ablation in He II. Furthermore we present a numerical study on the evolution of the thermodynamic state of the He II in front of a planar heater by solving the He II two-fluid equations of motion. For an applied heat flux, we determine the heating time beyond which the He II near the heater transits to the vapor phase. As such, a curve correlating $q^*$ and $t_h$ can be obtained, which nicely reproduces some relevant experimental data. Surprisingly, we find that there exists a critical peak heat flux $q^*_c$, above which boiling occurs nearly instantaneously regardless of $t_h$. We reveal that the boiling in this regime is essentially cavitation caused by the combined effects of the first-sound and the second-sound waves in He II. Based on this physical picture, an analytical model for $q^*_c$ is developed, which reproduces the simulated $q^*_c$ values at various He II bath temperatures and hydrostatic head pressures. This work represents a major progress in our understanding of transient heat transfer in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Radial Compressor Design and Off-Design for Trans-critical CO2 Operating Conditions

Presentation for the The 7th International Supercritical CO2 Power Cycles Symposium, February 21 – 24, 2022, San Antonio, Texas. A multi-staged intercooled centrifugal compressor design is performed for an indirect-fired sCO2 cycle and the design point data are used to generate the performance maps using a 2D mean-line analysis code. Fluid and turbomachinery similarity laws are applied to find an appropriate surrogate ideal gas fluid to be used in the design and performance map generation codes. The performance maps are then scaled back for the sCO2 fluid case using the affinity laws.

Uysal, Selcuk Can↗
NASA NTRS2015

Non-Linear Slosh Damping Model Development and Validation

Propellant tank slosh dynamics are typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control (GN&C) analysis. For a partially-filled smooth wall propellant tank, the critical damping based on classical empirical correlation is as low as 0.05%. Due to this low value of damping, propellant slosh is potential sources of disturbance critical to the stability of launch and space vehicles. It is postulated that the commonly quoted slosh damping is valid only under the linear regime where the slosh amplitude is small. With the increase of slosh amplitude, the critical damping value should also increase. If this nonlinearity can be verified and validated, the slosh stability margin can be significantly improved, and the level of conservatism maintained in the GN&C analysis can be lessened. The purpose of this study is to explore and to quantify the dependence of slosh damping with slosh amplitude. Accurately predicting the extremely low damping value of a smooth wall tank is very challenging for any Computational Fluid Dynamics (CFD) tool. One must resolve thin boundary layers near the wall and limit numerical damping to minimum. This computational study demonstrates that with proper grid resolution, CFD can indeed accurately predict the low damping physics from smooth walls under the linear regime. Comparisons of extracted damping values with experimental data for different tank sizes show very good agreements. Numerical simulations confirm that slosh damping is indeed a function of slosh amplitude. When slosh amplitude is low, the damping ratio is essentially constant, which is consistent with the empirical correlation. Once the amplitude reaches a critical value, the damping ratio becomes a linearly increasing function of the slosh amplitude. A follow-on experiment validated the developed nonlinear damping relationship. This discovery can lead to significant savings by reducing the number and size of slosh baffles in liquid propellant tanks.

Yang, H. Q.↗
DOE OSTI2021

On the feasibility of using physics-informed machine learning for underground reservoir pressure management

In this work, we evaluate the feasibility of using physics-informed machine learning (PIML) for underground energy-related pressure management. To this end, we develop a PIML framework to manage underground reservoir pressures by training neural networks to determine fluid extraction rates for dedicated extraction wells during fluid injection operations given a range of reservoir conditions (e.g., transmissivity and storativity). We implement an automatically-differentiable analytical physics model of fluid flow in porous media within the PIML framework as a proxy for more complicated models. This allows us to execute a sufficient number of training scenarios to fully evaluate the feasibility of using PIML to support pressure management activities. We quantify the number of physics-model parameters required for automatic differentiation to become more efficient than finite-difference gradient calculations. We use a simple scenario with a single injector, extractor, and critical location for our feasibility analysis. We evaluate the effect of the size of the training dataset (i.e., the number of reservoir condition samples) on the accuracy and efficiency of the PIML framework. For an equivalent number of model evaluations, the larger training dataset took less time to train and produced a neural network that was able to more accurately manage reservoir pressures. We also evaluate the effect of the training dataset batch size (i.e., number of reservoir condition samples used to update the neural network coefficients during training; i.e., how the training dataset is partitioned). While training ran faster with larger batch sizes, they produced neural networks that managed pressures less accurately. We demonstrate the approach on a more complex scenario involving 10 injectors, 10 extractors, and 4 critical locations (a relatively high well density of 20 wells/km2). We provide the number of forward and adjoint model evaluations required in each case as an indication of the feasibility of using PIML for pressure management when more complicated physics models with longer execution times are used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Probing oil recovery in shale nanopores with small-angle and ultra-small-angle neutron scattering

Increasing oil production from unconventional shale reservoirs is crucial to meet growing energy demands while achieving lower carbon emission than conventional crude oil. Enhanced oil recovery (EOR) has been proposed to improve hydrocarbon recovery rates through the injection of a fluid into the reservoir to facilitate residual oil release from the shale formation. However, economical and sustainable implementation of EOR requires advanced knowledge of fluid behavior in nano-sized pore spaces in shale. In this study, we utilize small-angle neutron scattering (SANS) and ultra-small-angle neutron scattering (USANS) as experimental probes to examine decane removal from a shale nano- to micro-porous matrix, utilizing methane as the injectant. The extent of decane saturation and recovery post-methane pressurization is quantified for clay-rich and carbonate-rich shale samples. A key finding is that extraction of decane by methane on depressurization is related to the methane-decane critical point. Furthermore, we found that although clay-rich shale had a much higher porosity of 5.6%, compared with 1.2% for carbonate-rich shale, decane was more easily removed from the carbonate-rich matrix, leading to similar hydrocarbon yields. These promising results demonstrate the ability of SANS and USANS to provide key insights into oil recovery from nano- to micron-sized pores in shale matrices. Combined with effects of various fractures on fluid behavior in shale, this experimental technique can be used to assess the viability of EOR injectants.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Asc-Seurat: analytical single-cell Seurat-based web application

Abstract Background Single-cell RNA sequencing (scRNA-seq) has revolutionized the study of transcriptomes, arising as a powerful tool for discovering and characterizing cell types and their developmental trajectories. However, scRNA-seq analysis is complex, requiring a continuous, iterative process to refine the data and uncover relevant biological information. A diversity of tools has been developed to address the multiple aspects of scRNA-seq data analysis. However, an easy-to-use web application capable of conducting all critical steps of scRNA-seq data analysis is still lacking. Summary We present Asc-Seurat, a feature-rich workbench, providing an user-friendly and easy-to-install web application encapsulating tools for an all-encompassing and fluid scRNA-seq data analysis. Asc-Seurat implements functions from the Seurat package for quality control, clustering, and genes differential expression. In addition, Asc-Seurat provides a pseudotime module containing dozens of models for the trajectory inference and a functional annotation module that allows recovering gene annotation and detecting gene ontology enriched terms. We showcase Asc-Seurat’s capabilities by analyzing a peripheral blood mononuclear cell dataset. Conclusions Asc-Seurat is a comprehensive workbench providing an accessible graphical interface for scRNA-seq analysis by biologists. Asc-Seurat significantly reduces the time and effort required to analyze and interpret the information in scRNA-seq datasets.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

LES simulations of a vacuum membrane distillation channel with geometric alterations

3D LES simulations were carried out to study the performance of a vacuum membrane distillation module. A wiggly wall profile with/without embedded stiffeners was considered to alleviate axial and radial temperature polarization, the cause of performance decrease in membrane distillation. Results of the flow field show the wiggles create unsteady vortex shedding inducing intense mixing in the channel. Vortex shedding intensity increases as the Reynolds number increases or stiffeners are added. 3D results show a 56% improvement in flux when moving from a flat sheet membrane to a wiggly membrane with stiffeners at a constant Reynolds number, corresponding to an improvement from 11.7 to 41.7 in the Nusselt number, showing that the alterations improved the flux performance by enhancing the heat transfer along the membrane surface and therefore alleviating temperature polarization, a critical bottleneck in membrane distillation systems. Further, a merit criterion was defined based on the Nusselt number and friction factor, and a 40% increase in merit was shown switching from a flat to a wiggly channel, while a 97% merit increase was seen going from flat to wiggly with stiffeners. A 42% enhancement in the flux was also seen moving from a straight channel to a wiggly channel at a higher Reynolds number which highlights the importance of strategically choosing the mass flow rate, as well as inducing flow separation and vortex shedding in the channel to promote mixing and dissipate the thermal boundary layer. The variation in results between the wiggly module with stiffeners and the wiggly module at high Reynolds numbers suggests enhancing mixing structures can be more impactful on flux for this geometry than increasing the flow rate to a turbulent/transitional regime. However, both are preferable for peak system flux performance. Furthermore, a 2D approximation was used to perform simulations on more extended channels to examine the length degradation. The modules with wiggly channels performed at the same flux level with a doubling of the length. On the other hand, the flat sheet modules experienced length degradation by temperature polarization and dropped in flux yield by around 13% of the short-channel value. This work illustrates that modeling the system and understanding how the performance decreases as the membrane surface area increases are critical for a larger module (scaling up from a lab to a prototype module) to maintain high flux performance.

42 ENGINEERING↗
DOE OSTI2021

Dynamic Phase Alignment in Inertial Alfven Turbulence

In weakly-collisional plasma environments with sufficiently low electron beta, Alfvenic turbulence transforms into inertial Alfvenic turbulence at scales below the electron skin-depth. We argue that, in inertial Alfvenic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale-dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as proportional to the perpendicular wavenumber to the power of -1. The kinetic and magnetic energy spectra scale as the perpendicular wavenumber to the power of -5/3 and -11/3, respectively. As a result of the dual direct cascade, the generalized-helicity spectrum scales as the perpendicular wavenumber to the power of -5/3, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the sub-electron-skin-depth range. Turbulent eddies exhibit a phase-space anisotropy, with the parallel wavenumber scaling as the perpendicular wavenumber to the power of 5/3, consistent with critically-balanced inertial Alfven fluctuations. Our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth's magnetosheath and ionosphere, solar corona, non-relativistic pair plasmas, as well as to strongly rotating non-ionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Thermal catastrophe in the plasma sheet boundary layer

This letter presents a first step towards a substorm model including particle heating and transport in the plasma sheet boundary layer (PSBL). The heating mechanism discussed is resonant absorption of Alfven waves. For some assumed MHD perturbation incident from the tail lobes onto the plasma sheet, the local heating rate in the PSBL has the form of a resonance function of the one-fluid plasma temperature. Balancing the local heating by convective transport of the heated plasma toward the central plasma sheet, an 'equation of state" is found for the steady-state PSBL whose solution has the form of a mathematical catastrophe: at a critical value of a parameter containing the incident power flux, the local density, and the convection velocity, the equilibrium temperature jumps discontinuously. Associating this temperature increase with the abrupt onset of the substorm expansion phase, the catastrophe model indicates at least three ways in which the onset may be triggered. Several other consequences related to substorm dynamics are suggested by the simple catastrophe model.

Smith, Robert A.↗
NASA NTRS2001

Enhancement of Pool Boiling Heat Transfer and Control of Bubble Motion in Microgravity Using Electric Fields (BCOEL)

The BCOEL project focuses on improving pool boiling heat transfer and bubble control in microgravity by exposing the fluid to electric fields. The electric fields induce a body force that can replace gravity in the low gravity environment, and enhance bubble removal from the heated surface. A better understanding of microgravity effects on boiling with and without electric fields is critical to the proper design of the phase-change-heat-removal equipment for use in spacebased applications. The microgravity experiments will focus on the visualization of bubble formation and shape during boiling. Heat fluxes on the boiling surface will be measured, and, together with the measured driving temperature differences, used to plot boiling curves for different electric field magnitudes. Bubble formation and boiling processes were found to be extremely sensitive to g-jitter. The duration of the experimental run is critical in order to achieve steady state in microgravity experiments. The International Space Station provides conditions suitable for such experiments. The experimental apparatus to be used in the study is described in the paper. The apparatus will be tested in the KC-135 first, and microgravity experiments will be conducted on board of the International Space Station using the Microgravity Science Glovebox as the experimental platform.

Herman, Cila↗
NASA NTRS2001

Enhancement of Pool Boiling Heat Transfer and Control of Bubble Motion in Microgravity Using Electric Fields - BCOEL

The BCOEL project focuses on improving pool boiling heat transfer and bubble control in microgravity by exposing the fluid to electric fields. The electric fields induce a body force that can replace gravity in the low gravity environment, and enhance bubble removal from thc heated surface. A better understanding of microgravity effects on boiling with and without electric fields is critical to the proper design of the phase-change-heat-removal equipment for use in space-based applications. The microgravity experiments will focus on the visualization of bubble formation and shape during boiling. Heat fluxes on the boiling surface will be measured, and, together with the measured driving temperature differences, used to plot boiling curvcs for different electric field magnitudes. Bubble formation and boiling processes were found to be extremely sensitive to g-jitter. The duration of the experimental run is critical in order to achieve steady state in microgravity experiments. The International Space Station provides conditions suitable for such experiments. The experimental appararus to be used in the study is described in the paper. The apparatus will be tested in the KC-135 first, and microgravity experiments will be conducted on board of the International Space Station using the Microgravity Science Glovebox as the experimental platform.

Herman, Cila↗
NASA NTRSDate not supplied

Predicting Quadcopter Noise With The Lattice-Boltzmann Method

From small delivery Unmanned Aircraft Systems (UAS) to short-haul building-to-building aircraft, urban air mobility is gaining traction as a new approach to tackle the problem of transportation in densely populated areas. NASA’s vision for vertical lift vehicles is to capitalize on and improve unique capabilities to greatly benefit the United States’ growing civil flight requirements. This vision is embodied in the Revolutionary Vertical Lift Technology (RVLT)project [1]. Beyond safety, one of the chief concerns of communities where drones are becoming more popular is the noise they generate. Noise will undoubtedly be one of the major obstacles to public acceptance of any new urban air mobility technology. The ability to predict the acoustics of new conceptual aircraft with multiple rotors and complex fuselages is critical to enable the creation of quieter designs. The objective of this research is to build up a better physical understanding of the noise generated by a typical quadcopter drone and what it takes to predict it from first principles using computational fluid dynamics (CFD) with the Lattice-Boltzmann method (LBM). The specific goals are to establish best practices to predict multi-rotor and vehicle interaction noise with LBM, validate these predictions by comparing to wind tunnel measurements, and assess the computational cost necessary to obtain accurate predictions.

Francois Cadieux↗
NASA NTRS2023

Neural Radiance Fields for Shadow Removal and Relighting of Satellite Imagery

Multi-view stereo photogrammetric techniques have been conventionally utilized to generate Global Digital Elevation Models (GDEM) of planetary and lunar surfaces. However, these methods, relying on conventional feature detectors, are often subject to inaccuracies caused by changes in lighting conditions, including diffuse reflection and harsh shading. This has limited the ability of these methods to accurately reconstruct shadowed regions in orbital imagery, such as highly shaded urban areas and the permanently shadowed regions (PSRs) located in the lunar surface, which are critical targets for NASA's Artemis program. Neural Radiance Fields (NeRFs) offer a novel solution to these limitations by breaking away from traditional photogrammetric assumptions of ridged, opaque surfaces. NeRFs are capable of reconstructing 3D objects with variably transmissive properties, such as fluids, PSRs, and reflective surfaces. In this work, we demonstrate the robustness of NeRFs in generating high-fidelity 3D models of terrain from highly shaded orbital imagery acquired from small satellites in low earth orbit (LEO). Our results showcase the ability of NeRF to overcome the limitations of traditional photogrammetric methods and provide a promising solution for reconstructing complex scenes in challenging lighting conditions.

Caleb Ashmore Adams↗
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↗