SEARCH · Engineering Papers

Results for “Multiphase flows”

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

SourceResearch connectionFrom yearThrough year

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

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Online ion mobility spectrometry of nanoparticle formation by non-thermal plasma conversion of metal salts in liquid aerosol droplets

The synthesis of nanoparticles by reaction of liquid aerosol droplets containing precursors in a flow-through, atmospheric-pressure, non-thermal plasma offers a continuous, scalable, substrate- and stabilizer-free approach for direct deposition into liquids or onto soft substrates. However, the combination of multiphase and non-equilibrium chemistry makes the process complicated and poorly understood. Here, we present ion mobility spectrometry measurements of liquid water droplets containing silver nitrate passing through an atmospheric-pressure dielectric barrier discharge reactor that allows us to monitor silver nanoparticle formation online for the first time. Mobility diameter distributions were obtained with the plasma on and off, and exhibited a shift, which was related to the degree of conversion of silver nitrate. The silver nanoparticles were also collected and characterized by UV–visible absorbance spectroscopy and transmission electron microscopy to support the online measurements. Importantly, negligible conversion was found when the water was removed by a diffusion dryer, suggesting that the key reducing species are in the liquid phase, such as solvated electrons. Finally, the study demonstrates how ion mobility spectrometry measurements can be applied to provide insight into this approach to nanoparticle synthesis.

42 ENGINEERING↗

DOE OSTI2024

Machine Learning-Enhanced Multiphase CFD for Carbon Capture Modeling Run Data

Repository for the data generated as part of the 2023-2024 ALCC project "Machine Learning-Enhanced Multiphase CFD for Carbon Capture Modeling." The data was generated with MFIX-Exa's CFD-DEM model. The problem of interest is gravity driven, particle-laden, gas-solid flow in a triply-periodic domain of length 2048 particle diameters with an aspect ratio of 4. The mean particle concentration ranges from 1% to 40% and the Archimedes number ranges from 18 to 90. The particle-to-fluid density ratio, particle-particle restitution and friction coefficients and domain aspect ratio are held constant at values of 1000, 0.9, 0.25 and 4, respectively. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231 using NERSC award ALCC-ERCAP0025948.

AMReX↗

DOE OSTI2019

Investigation of Gravity-Driven Infiltration Instabilities in Smooth and Rough Fractures Using a Pairwise-Force Smoothed Particle Hydrodynamics Model

In this work, we use a three-dimensional multiphase pairwise-force smoothed particle hydrodynamics (PF-SPH) method to simulate infiltration dynamics in smooth and rough fractures. We investigate the influence of roughness and injection rate on the fluid flow profile and flow velocity. Three types of fractures are considered with different degrees of roughness, including a smooth fracture. Both the rough and smooth fractures exhibit flow instability, fingering, and intermittent flow regimes for low infiltration rates. A stable fluid profile only is achieved for high fluxes when fractures are fully saturated. In all considered fractures, the mean flow velocity exponentially increases with increasing infiltration rate. To characterize the flow instability, we calculate standard deviations of the velocity in different sections of the fracture relative to the mean flow velocity. For the considered infiltration rates, we find that roughness decreases the flow velocity and increases the velocity standard deviations. Hence, we observe that the likelihood of discontinuities in flow regimes increases with the degree of roughness.

58 GEOSCIENCES↗

DOE OSTI2016

A spectral cascade-transport turbulence model for wall-bounded flows

A spectral cascade-transport model has been developed and applied to turbulent channel flows (Re{sub τ} = 550, 950, and 2000 based on friction velocity, u{sub τ}; or Re{sub δ} = 8,500; 14,800 and 31,000, based on the mean velocity and channel half-width). This model is an extension of a spectral model previously developed for homogeneous single and two-phase decay of isotropic turbulence and uniform shear flows; and a spectral turbulence model for wall-bounded flows without resolving the boundary layer. A direct numerical simulation (DNS) of turbulent channel flow was used to help develop this model and to assess its performance in the 1D direction across the channel width. The resultant spectral model is capable of predicting the mean velocity, turbulent kinetic energy and energy spectrum distributions for single-phase wall-bounded flows all the way to the wall, where the model source terms have been developed to account for the wall influence. The model is in the process of implementation into the 3D multiphase CFD code NPHASE-CMFD and the latest results are within reasonable error of the 1D predictions. (authors)

42 ENGINEERING↗

DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2019

Highly resolved Eulerian simulations of fuel spray transients in single and multi-hole injectors: Nozzle flow and near-exit dynamics

In high pressure fuel injectors, needle opening and closing transients cause complex off-design fluid dynamics behaviors that profoundly impact the spray and mixture formation processes. These dynamics are completely different from what is known to occur in steady state conditions. In this study, diesel spray transients have been investigated in single-hole and 3-hole nozzles, encompassing internal and external nozzle flow and including needle motion, performing highly resolved (2.5 μm) computational fluid dynamics (CFD) simulations. We focused on end-of-injection (EOI) and start-of-injection (SOI) processes, in order to provide insights in to the physics. The liquid fuel, vapor and gas species are modeled with a single-fluid multiphase mixture approach, with diffuse interface, and with large eddy simulations (LES) of the turbulence. Occurrence of phase change due to cavitation is accounted for, and the spray dispersion is described with a turbulent dispersion model. Detailed needle motion data and orifice internal surface are available from xray synchrotron source measurements carried out at Argonne National Laboratory, and shared through the Engine Combustion Network (ECN) community. Simulations are compared against x-ray phase contrast imaging and radiography of the internal and near-exit flow, in addition to optical microscopy data of the near-exit sprays. Simulation results are found to agree well with available experimental data, and are able to realistically capture local and global features. The simulations allow to gain insight into the physics of gas ingestion and dribbles at EOI, for different hole diameters, operating conditions and number of holes. At SOI, timing of liquid appearance out of the injector and spray tip penetration are adequately predicted, by using the EOI flow field as in-nozzle initialization, and by prescribing the measured tip needle displacement with an informed effective valve opening point inferred from the x-ray observations. Lastly, the variation of spreading angle over time is also discussed in detail for the multi-hole case, including hole-to-hole variations. Due to real geometry features and asymmetric needle motion with eccentric components, it is found that the three holes exhibit swirling flows of increasing intensity as the lift decreases, causing the near cone angle to open and spread, in a quasi-hollow cone structure. These features are not observed in axial single-hole injectors because of their relative simplicity and intrinsic symmetry.

42 ENGINEERING↗

DOE OSTI2017

CFD-DEM Onset of Motion Analysis for Application to Bed Scour Risk Assessment

This CFD study with DEM was done as a part of the Federal Highway Administration’s (FHWA’s) effort to improve scour design procedures. The Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) model, available in CD-Adapco’s StarCCM+ software, was used to simulate multiphase systems, mainly those which combine fluids and solids. In this method the motion of discrete solids is accounted for by DEM, which applies Newton's laws of motion to every particle. The flow of the fluid is determined by the local averaged Navier–Stokes equations that can be solved using the traditional CFD approach. The interactions between the fluid phase and solids phase are modeled by use of Newton's third law. The inter-particle contact forces are included in the equations of motion. Soft-particle formulation is used, which allows particles to overlap. In this study DEM was used to model separate sediment grains and spherical particles laying on the bed with the aim to analyze their movement due to flow conditions. Critical shear stress causing the incipient movement of the sediment was established and compared to the available experimental data. An example of scour around a cylindrical pier is considered. Various depths of the scoured bed and flow conditions were taken into account to gain a better understanding of the erosion forces existing around bridge foundations. The decay of these forces with increasing scour depth was quantified with a ‘decay function’, which shows that particles become increasingly less likely to be set in motion by flow forces as a scour hole increases in depth. Computational and experimental examples of the scoured bed around a cylindrical pier are presented.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Wellbore Fracture Imaging Using Inflow Detection Measurements

One of the most striking measurements taken during DOE’s EGS Collab project at the 4850-foot depth location was the so-called ‘sewer cam’, which enabled direct visualization of the flow of water into the production well through fractures during the stimulation. The ability to see directly which fractures were flowing and (roughly) how much was a breakthrough in understanding the topology of the created fracture network. Achieving this kind of fracture flow imaging at FORGE would be more challenging because of the 225°C temperature, but equally or even more valuable if it could be achieved. In 2017, a joint project between Sandia and Stanford developed a downhole tool concept to measure the enthalpy of multiphase fluid entering a geothermal well from individual fractures (Gao et al., 2017). For the FORGE project, measuring enthalpy is of less interest because the fluid is expected to be single-phase liquid water. However, the foundation of the device was the measurement of chloride ion concentration, which could form the basis for a direct measurement of inflow from fractures. During the 2017 project, this novel chloride sensing system was implemented into a laboratory test instrument, and we confirmed the capability of the system to measure the ion concentration of fluid entering a model wellbore through a small entry port. The wellbore was a 6-inch diameter model well, and the port was approximately 0.08 inch (2mm) in diameter. The device could measure the chloride concentration accurately even when the well was flowing in a bubbly flow. Given its accuracy, the tool should be able to identify locations of water entering the wellbore even if the ion concentration differs only slightly from that of the water in the well. It is likely that different fractures may flow slightly different chloride concentrations, which would make it feasible to detect individual fractures as well as to estimate the volume of their flow. Ultimately, we could also recognize different fractures flowing back significantly different ion concentrations after fracturing in the FORGE wells. This could be realized by adding different ions in the fracturing fluids in different fractures created at different stages of stimulation (and modifying the tool to include different ion specificity). Sandia’s tool was shown during the study to have the capability to withstand the 225°C temperature, and the electrochemical sensing elements were tested in the laboratory to 225°C at 1500 psia for 24 hours. An early implementation of the fully integrated downhole electrochemical tool, including high-temperature electronics, robust housing, and wireline truck interface, had previously been constructed and tested successfully at Sandia; thus, hardware development tasks focused on advancing the technology readiness level (TRL) of this promising technology for FORGE deployment, rather than on developing a new scientific basis for its operation. The data collection electronics in this tool allowed for several other sensors (pressure, temperature, flow spinner) to be implemented in parallel as well. The research was a new collaboration between Stanford and Sandia to modify and refine the tool for FORGE deployment, to make the downhole measurements, and to characterize the evolving fractures.

15 GEOTHERMAL ENERGY↗

DOE OSTI2015

Lattice Boltzmann simulation of immiscible fluid displacement in porous media: Homogeneous versus heterogeneous pore network

Injection of anthropogenic carbon dioxide (CO{sub 2}) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO{sub 2} that can be captured and its long-term storage stability in subsurface requires a fundamental understanding of multiphase displacement phenomena at the pore scale. In this paper, the lattice Boltzmann method is employed to simulate the immiscible displacement of a wetting fluid by a non-wetting one in two microfluidic flow cells, one with a homogeneous pore network and the other with a randomly heterogeneous pore network. We have identified three different displacement patterns, namely, stable displacement, capillary fingering, and viscous fingering, all of which are strongly dependent upon the capillary number (Ca), viscosity ratio (M), and the media heterogeneity. The non-wetting fluid saturation (S{sub nw}) is found to increase nearly linearly with logCa for each constant M. Increasing M (viscosity ratio of non-wetting fluid to wetting fluid) or decreasing the media heterogeneity can enhance the stability of the displacement process, resulting in an increase in S{sub nw}. In either pore networks, the specific interfacial length is linearly proportional to S{sub nw} during drainage with equal proportionality constant for all cases excluding those revealing considerable viscous fingering. Our numerical results confirm the previous experimental finding that the steady state specific interfacial length exhibits a linear dependence on S{sub nw} for either favorable (M ≥ 1) or unfavorable (M < 1) displacement, and the slope is slightly higher for the unfavorable displacement.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1972

Research on turbojet combustor pollutant formation

The research for developing a modular computer program for describing gas-turbine combustor flow and chemical fields is reported. The status of investigations is presented for three key elements of the unified model: turbulent mixing with swirl, turbulent mixing with swirl including finite-rate kinetics, and multiphase recirculation zone analysis

Source record↗

NASA NTRS1993

In-flight pressure distributions and skin-friction measurements on a subsonic transport high-lift wing section

Flight experiments are being conducted as part of a multiphased subsonic transport high-lift research program for correlation with wind-tunnel and computational results. The NASA Langley Transport Systems Research Vehicle (B737-100 aircraft) is used to obtain in-flight flow characteristics at full-scale Reynolds numbers to contribute to the understanding of 3-D high-lift, multi-element flows including attachment-line transition and relaminarization, confluent boundary-layer development, and flow separation characteristics. Flight test results of pressure distributions and skin friction measurements were obtained for a full-chord wing section including the slat, main-wing, and triple-slotted, Fowler flap elements. Test conditions included a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.40). Pressure distributions were obtained at 144 chordwise locations of a wing section (53-percent wing span) using thin pressure belts over the slat, main-wing, and flap elements. Flow characteristics observed in the chordwise pressure distributions included leading-edge regions of high subsonic flows, leading-edge attachment-line locations, slat and main-wing cove-flow separation and reattachment, and trailing-edge flap separation. In addition to the pressure distributions, limited skin-friction measurements were made using Preston-tube probes. Preston-tube measurements on the slat upper surface suggested relaminarization of the turbulent flow introduced by the pressure belt on the slat leading-edge surface when the slat attachment line was laminar. Computational analysis of the in-flight pressure measurements using two-dimensional, viscous multielement methods modified with simple-sweep theory showed reasonable agreement. However, overprediction of the pressures on the flap elements suggests a need for better detailed measurements and improved modeling of confluent boundary layers as well as inclusion of three-dimensional viscous effects in the analysis.

Yip, Long P.↗

DOE OSTI2020

Relative Permeability for Offshore Enhanced Oil Recovery

Relative permeability (kr ) is the description of multiphase transport through porous media that is most widely accepted and utilized to scale relationships up to the field scale through simulations. Previous research at NETL has shown a dependence of kr on the flow rates and porous media structure that is poorly captured in most descriptions of this process. Our project was originally scoped to be a 3-year project to (2019-2021): (1) determine if this poor literature description is true for offshore EOR; (2) collect data on the generation of relevant kr curves for offshore EOR; and (3) distribute this collected data, methodology and resultant curves.

Crandall, Dustin↗

DOE OSTI2015

TOUGH+ v1.5 Core Code

TOUGH+ v1.5 is a numerical code for the simulation of multi-phase, multi-component flow and transport of mass and heat through porous and fractured media, and represents the third update of the code since its first release [Moridis et al., 2008]. TOUGH+ is a successor to the TOUGH2 [Pruess et al., 1991; 2012] family of codes for multi-component, multiphase ﬂuid and heat ﬂow developed at the Lawrence Berkeley National Laboratory. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstations, PC, Macintosh). TOUGH+ v1.5 employs dynamic memory allocation, thus minimizing storage requirements. It has a completely modular structure, follows the tenets of Object-Oriented Programming (OOP), and involves the advanced features of FORTRAN 95/2003, i.e., modules, derived data types, the use of pointers, lists and trees, data encapsulation, defined operators and assignments, operator extension and overloading, use of generic procedures, and maximum use of the powerful intrinsic vector and matrix processing operations. TOUGH+ v1.5 is the core code for its family of applications, i.e., the part of the code that is common to all its applications. It provides a description of the underlying physics and thermodynamics of non-isothermal flow, of the mathematical and numerical approaches, as well as a detailed explanation of the general (common to all applications) input requirements, options, capabilities and output specifications. The core code cannot run by itself: it needs to be coupled with the code for the specific TOUGH+ application option that describes a particular type of problem. The additional input requirements specific to a particular TOUGH+ application options and related illustrative examples can be found in the corresponding User's Manual.

Moridis, GeorgeJ↗

DOE OSTI2018

Multiphase Model of Semisolid Slurry Generation and Isothermal Holding During Cooling Slope Rheoprocessing of A356 Al Alloy

In the present paper, we present an experimentally validated 3D multiphase and multiscale solidification model to understand the transport processes involved during slurry generation with a cooling slope. In this process, superheated liquid alloy is poured at the top of the cooling slope and allowed to flow along the slope under the influence of gravity. As the melt flows down the slope, it progressively loses its superheat, starts solidifying at the melt/slope interface with formation of solid crystals, and eventually exits the slope as semisolid slurry. In the present simulation, the three phases considered are the parent melt as the primary phase, and the solid grains and air as secondary phases. The air phase forms a definable air/liquid melt interface as the free surface. After exiting the slope, the slurry fills an isothermal holding bath maintained at the slope exit temperature, which promotes further globularization of microstructure. The outcomes of the present model include prediction of volume fractions of the three different phases considered, grain evolution, grain growth, size, sphericity and distribution of solid grains, temperature field, velocity field, macrosegregation and microsegregation. In addition, the model is found to be capable of making predictions of morphological evolution of primary grains at the onset of isothermal coarsening. The results obtained from the present simulations are validated by performing quantitative image analysis of micrographs of the rapidly oil-quenched semisolid slurry samples, collected from strategic locations along the slope and from the isothermal slurry holding bath.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Modeling Multiphase Effects in CO 2 Compressors at Subcritical Inlet Conditions

An advanced numerical framework to model CO 2 compressors over a wide range of subcritical conditions is presented in this paper. Thermodynamic and transport properties are obtained through a table look-up procedure with specialized features for subcritical conditions. Phase change is triggered by the difference between the local values of pressure and saturation pressure, and both vaporization and condensation can be modeled. Rigorous validation of the framework is presented for condensation in high pressure CO 2 using test data in a De Laval nozzle. The comparisons between computations and test data include: condensation onset locations, Wilson line, and nozzle pressure profiles as a function of inlet pressures. The framework is applied to the Sandia compressor that has been modeled over broad range of conditions spanning the saturation dome including: near critical inlet conditions (305.4 K, and 7.843 MPa), pure liquid inlet conditions (at 295 K), pure vapor inlet conditions (at 302 K), and two-phase inlet conditions (at 290 K). Multiphase effects ranging from cavitation at the liquid line to condensation at the vapor line have been simulated. Here, the role of real fluid effects in enhancing multiphase effects at elevated temperatures closer to the critical point has been identified. The performance of the compressor has been compared with test data; the computational fluid dynamics (CFD) results also show that the head-flow coefficient curve collapses with relatively minor scatter, similar to the test data, when the flow coefficient is defined on the impeller exit meridional velocity.

14 SOLAR ENERGY↗

DOE OSTI2022

An Open-Source Framework for the Computational Analysis and Design of Autothermal Chemical Processes

The main project object was to develop software tools for simulating non-equilibrium autothermal processes, improving the prospects for identifying and designing such systems. The project demonstrates the use of these tools to simulate autothermal pyrolysis, a process recently developed at the pilot-plant scale at Iowa State University. In such process, instead of externally heating a reactor to pro-vide the enthalpy of pyrolysis, sufficient oxygen in the form of air is introduced into the reactor to support partial oxidation of reactants and products with the exothermic energy released supporting endothermic pyrolysis reactions. A fluidized bed is used to assure good mixing of biomass and oxidant and provide an isothermal reaction environment. The amount of oxygen required depends upon the kind of biomass being pyrolyzed and parasitic heat losses from the reactor. For example, for woody biomass pyrolyzed under conditions that simulate adiabatic operation, equivalence ratios can be as low as 0.06, compared to 0.20 or higher for autothermal gasifiers. By removing the heat transfer bottleneck of conventional pyrolysis, operation in autothermal mode allowed a significant increase in reactor throughput process, approaching five times the throughput of the conventionally operated pyrolyzer. Different simulation strategies were considered and developed: a zero-dimensional chemistry model was used to verify the applicability of kinetic schemes to predict biomass fast pyrolysis in autothermal conditions. Conventional chemical reactor models such as the plug flow reactor and the partially stirred reactor were used to investigate the role of mixing in the fluidized bed pyrolyzer and to establish the impact of mixing time on the gas-phase reactions. A comprehensive multiphase computational fluid dynamics (mCFD) framework, including polydisperse granular phase modeling and detailed chemical kinetics was formulated and used to model the experimental setup for autothermal biomass fast pyrolysis at ISU. Multiphase CFD was also used to investigate the role of biomass feed positioning on the mixing of biomass in the pyrolizer. Finally, a reduced order model (ROM), suitable to be implemented in process simulators was obtained. Both the mCFD and the ROM were validated against experiments.

09 BIOMASS FUELS↗

DOE OSTI2023

Sizing and Selection of Pressure Relief Valves for High-Pressure Thermal–Hydraulic Systems

This study covers the critical concerns related to the sizing, selection, installation, maintenance, and testing of pressure safety valves (PSVs). The aim is to ensure the safety of pressurized systems, hydrostatic transmission systems, and hydraulic plants, including process plants, thermal power plants, and nuclear reactor systems. PSVs are devices that ensure the safety and reliability of pressurized vessels, lines, and systems during overpressure events. The task of selecting which PSV features are of greatest value for a specific purpose is complex—especially in the design of a high-pressure experimental thermal–hydraulic facility for hydrostatic and transient testing of the reactor system—when the systems are in the design and development phases and require qualification and demonstration to prove that they have reached a given level of technological readiness. The present study highlights the required steps for users to follow the associated rules, guidelines, and recommendations. As a part of this research, case studies are presented to help readers better understand the applicable strategy and standards. A discussion and a review of PSV performance degradation and failure are summarized to provide a better understanding of varied process applications and conditions, including fluid flow dynamics, boundary-layer formation and pressure drops, gas bubble formation and collapse, geometric configurations, inlet/outlet piping, abrupt pressure fluctuations, and acoustic resonance. Moreover, this study discusses the servicing and testing of PSVs in a multiphase pressurized system. Overall, it provides a basic overview of how PSVs ensure the safety of pressurized systems, supported by case studies and industrial practices.

02 PETROLEUM↗

DOE OSTI2016

CFD Study of Full-Scale Aerobic Bioreactors: Evaluation of Dynamic O2 Distribution, Gas-Liquid Mass Transfer and Reaction

If advanced biofuels are to measurably displace fossil fuels in the near term, they will have to operate at levels of scale, efficiency, and margin unprecedented in the current biotech industry. For aerobically-grown products in particular, scale-up is complex and the practical size, cost, and operability of extremely large reactors is not well understood. Put simply, the problem of how to attain fuel-class production scales comes down to cost-effective delivery of oxygen at high mass transfer rates and low capital and operating costs. To that end, very large reactor vessels (>500 m3) are proposed in order to achieve favorable economies of scale. Additionally, techno-economic evaluation indicates that bubble-column reactors are more cost-effective than stirred-tank reactors in many low-viscosity cultures. In order to advance the design of extremely large aerobic bioreactors, we have performed computational fluid dynamics (CFD) simulations of bubble-column reactors. A multiphase Euler-Euler model is used to explicitly account for the spatial distribution of air (i.e., gas bubbles) in the reactor. Expanding on the existing bioreactor CFD literature (typically focused on the hydrodynamics of bubbly flows), our simulations include interphase mass transfer of oxygen and a simple phenomenological reaction representing the uptake and consumption of dissolved oxygen by submerged cells. The simulations reproduce the expected flow profiles, with net upward flow in the center of column and downward flow near the wall. At high simulated oxygen uptake rates (OUR), oxygen-depleted regions can be observed in the reactor. By increasing the gas flow to enhance mixing and eliminate depleted areas, a maximum oxygen transfer (OTR) rate is obtained as a function of superficial velocity. These insights regarding minimum superficial velocity and maximum reactor size are incorporated into NREL's larger techno-economic models to supplement standard reactor design equations.

09 BIOMASS FUELS↗