Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid bubbles”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

CFD Analysis in Advance of the NASA Juncture Flow Experiment

NASA through its Transformational Tools and Technologies Project (TTT) under the Advanced Air Vehicle Program, is supporting a substantial effort to further investigate the formation and origin of separation bubbles found on wing-body juncture zones. The flow behavior in these regions is highly complex, difficult to measure experimentally, and challenging to model numerically. Multiple wing configurations were designed and evaluated using Computational Fluid Dynamics (CFD), and a series of wind tunnel risk reduction tests were performed to further down-select the candidates for the final experiment. This paper documents the CFD analysis done in conjunction with the 6 percent scale risk reduction experiment performed in NASA Langley’s 14- by 22-Foot Subsonic Tunnel. The combined CFD and wind tunnel results ultimately help the Juncture Flow committee select the wing configurations for the final experiment.

Lee, H. C.↗
NASA NTRS2012

Microgravity Passive Phase Separator

A new invention disclosure discusses a structure and process for separating gas from liquids in microgravity. The Microgravity Passive Phase Separator consists of two concentric, pleated, woven stainless- steel screens (25-micrometer nominal pore) with an axial inlet, and an annular outlet between both screens (see figure). Water enters at one end of the center screen at high velocity, eventually passing through the inner screen and out through the annular exit. As gas is introduced into the flow stream, the drag force exerted on the bubble pushes it downstream until flow stagnation or until it reaches an equilibrium point between the surface tension holding bubble to the screen and the drag force. Gas bubbles of a given size will form a front that is moved further down the length of the inner screen with increasing velocity. As more bubbles are added, the front location will remain fixed, but additional bubbles will move to the end of the unit, eventually coming to rest in the large cavity between the unit housing and the outer screen (storage area). Owing to the small size of the pores and the hydrophilic nature of the screen material, gas does not pass through the screen and is retained within the unit for emptying during ground processing. If debris is picked up on the screen, the area closest to the inlet will become clogged, so high-velocity flow will persist farther down the length of the center screen, pushing the bubble front further from the inlet of the inner screen. It is desired to keep the velocity high enough so that, for any bubble size, an area of clean screen exists between the bubbles and the debris. The primary benefits of this innovation are the lack of any need for additional power, strip gas, or location for venting the separated gas. As the unit contains no membrane, the transport fluid will not be lost due to evaporation in the process of gas separation. Separation is performed with relatively low pressure drop based on the large surface area of the separating screen. Additionally, there are no moving parts, and there are no failure modes that involve fluid loss. A patent application has been filed.

Paragano, Matthew↗
DOE OSTI2020

Multiphase Reacting Flow Simulations and Optimization of Commercial-Scale Aerobic Bioreactors

Aerobic bioreactors are widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries, wherein microbial bioreactions aid in the conversion of sugars to macromolecules. These reactors are also being actively investigated as a cost-effective pathway for the production of low-value commodities such as biofuels and animal feed. Reactor scale-up is one of the important challenges when designing these systems. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of commercial-scale (~ 500 m3) aerobic bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for oxygen uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM, to perform these simulations, which has been validated against small-scale reactors in our previous work. This work examines the performance of three different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor oxygen mass-transfer coefficient, gas hold-up, and oxygen distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented.

09 BIOMASS FUELS↗
DOE OSTI2022

State of Knowledge on Aerosols and Bubble Transport for Mechanistic Source Term Analysis of Molten Salt Reactors

A comprehensive review of the current state of knowledge surrounding molten salt aerosol behavior and bubble transport in molten salt systems relevant to mechanistic source term (MST) analyses of molten salt reactors (MSR) and fluoride high temperature reactors (FHR) was provided. Although relatively little experimental work has been carried out on investigating these phenomena in molten salts compared to that for other nuclear fluids of interest (e.g., aqueous, sodium), the amount of work that has been reported still provides a baseline understanding to inform and design future experimental campaigns on these topics. The bulk of this data is derived from the Molten Salt Reactor Program (MSRP) which was carried out at Oak Ridge National Laboratory (ORNL) nearly 50 years ago, leading to the successful operation of the Molten Salt Reactor Experiment (MSRE). Many other related research efforts are also highlighted in this review, such as recent efforts under the DOE:NE Advanced Reactor Technology (ART) MSR Campaign. The goal of this report is to provide a foundation for modelers to base assumptions on, as well as to make recommendations to experimentalists attempting to close information gaps. As such, Section 4 contains recommendations for future experimental work that may be needed to close current information gaps related to MST modeling of MSRs. Additionally, given the flowing fuel nature of MSRs, the relevance of aerosol and bubble transport behavior in molten salts goes beyond just MST modeling, evidenced by experiences with the pump bowl and off-gas system of the MSRE.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗
NASA NTRS1987

Characteristics of fluid flow in the combustion synthesis of TiC from the elements

The results of a numerical investigation of finite reservoir effects on capillary spreading at small reservoir dimensions are presently related to wave propagation phenomena in the combustion synthesis of TiC from its two elemental constituents. It is noted that gravitational forces can affect bubble coalescence by nonbuoyant means under the suitable conditions, although these conditions are expected to be rare in combustion synthesis. Finite-curved reservoirs can drive capillary flow due to surface tension and wall contact forces; these cause the wall and the metal to be completely reconfigured during combustion synthesis.

Valone, S. M.↗
NASA NTRS1986

Calculations of the stability of some axisymmetric flows proposed as a model of vortex breakdown

The term vortex breakdown refers to the abrupt and drastic changes of structure that can sometimes occur in swirling flows. It was conjectured that the bubble type of breakdown can be viewed as an axisymmetric wave traveling upstream in a primarily columnar vortex flow. In this scenario the wave's upstream progress is impeded only when it reaches a critical amplitude and it loses stability to some nonaxisymmetric disturbance. The stability of some axisymmetric wavy flows to three dimensional disturbances, viewing the amplitude of the wave as a bifurcation parameter is examined. The stability of a set of related columnar vortex flows, constructed by taking the two dimensional flow at a single axial location and extending it throughout the domain without variation, is investigated. The method used will be to expand the perturbation velocity in a series of divergence free vectors which ensures that the continuity equation for the incompressible fluid is satisfied exactly by the computed velocity field. Projections of the stability equation onto the space of inviscid vector fields eliminated the pressure term from the equation and reduces the differential eigen problem to a generalized matrix eigen problem. Results are presented both for the one dimensional, columnar vortex flows and also for the wavy bubble flow.

Mhuiris, N. M. G.↗
DOE OSTI2022

High-fidelity pool boiling simulations on multiple nucleation sites using interface capturing method

Boiling has proved to be one of the most efficient means for heat transfer and is a very important phenomenon during severe accident scenarios in light water reactors. High-fidelity pool boiling simulations can provide a numerical database for improving mechanistic boiling models by allowing for specific evaluation of interactions among bubbles. Previously published pool boiling simulations investigated two nucleation sites in which bubble growth at one site suppressed nucleation at the other site. Based on previous study results, more complicated interface-capturing simulations on pool boiling were conducted using PHASTA code with locally refined unstructured mesh. First, different boundary conditions (BCs) were assessed to support robustness and reproducibility of the boiling model. Then, a scale study was conducted at a larger domain with nine nucleation sites where either nine or four nucleation sites are activated. Involving more nucleation sites increased the complexity of bubble interactions from surrounding sites. Finally, bubble departure behavior influenced by wall heat flux was investigated. In this work, when heat flux was increased, the order of bubble departure changed, but diagonal bubbles always departed after one another. The departure time interval between the first and second bubble reduced as heat flux increased. The corresponding frequency was almost linearly proportional to the heat flux. In addition, bubble departure behavior was found to be greatly influenced by the nucleation site pattern. Multiple nucleation sites resulted in superimposed inhibitive effects from surrounding sites to each bubble, which extensively delayed the departure. This new observation was not discussed in previously published works. The work presented here provides new insight on the fundamental understanding of boiling phenomena, contributes to the development of a 3D multiphase computational fluid dynamics (M-CFD) model, and provides a more comprehensive database for data-driven pool boiling studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Bubble Point Measurements of cis-1,1,1,4,4,4-Hexafluorobutene [R-1336mzz(Z)] + trans-1,2-Dichloroethene [R-1130(E)] mixtures

Saturation pressures of pure R-1336mzz(Z) and R-1130(E) and bubble point pressures of three R-1336mzz(Z)/1130(E) blends were measured from 265 K to 360 K. For each pure refrigerant or refrigerant blend, a total of twenty unique saturation pressures or bubble points were measured. In total 100 unique state points were obtained. Presently, no Helmholtz-energy-explicit type equation of state (EoS) is available for R-1130(E). While an extended corresponding states EoS for R-1130(E) is available to estimate the properties of the R-1336mzz(Z)/1130(E) blend, this model does not resolve the azeotropic behavior of the mixture. Therefore, the perturbed-chain statistical-associating fluid theory (PC-SAFT) EoS is used to model the vapor–liquid equilibria of the R-1336mzz(Z)/1130(E) blend. PC-SAFT model parameters are reported, and the overall performance of the model is characterized by deviations from the experimental data.

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

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Calibration of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation with Machine Learning Surrogate Models

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is one of the most powerful accelerator-driven neutron sources in the world. The intense protons strike on SNS's mercury target to provide bright neutron beams, which also leads to severe fluid-structure interactions inside the target. Prediction of resultant loading on the target is difficult particularly when helium gas is injected into mercury to reduce the loading and mitigate the pitting damage on vessel walls. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this multi-physics problem. However, several uncertain parameters in the R-P model require intensive simulations to determine their optimal values. With the help of machine learning and the measured target strain, we have studied the major uncertain parameters in this R-P model and developed a framework to identify optimal parameters that significantly reduce the discrepancy between simulations and experimental strains. The preliminary results show the possibility of using this mercury/helium mixture and surrogate models to predict a better match of target strain response when the helium gas is injected.

Lin, Lianshan↗
NASA NTRS1993

Characterization of voids formed during liquid impregnation of nonwoven multifilament glass networks as related to composite processing

A technique based on matching the refractive index of an invading liquid to that of a fiber mat was used to study entrapment of air ('voids') that occurs during forced in-plane radial flow into nonwoven multifilament glass networks. The usefulness of this technique is demonstrated in quantifying and mapping the air pockets. Experiments with a series of fluids with surface tensions varying from 28 x 10(exp -3) to 36 x 10(exp -3) N/m, viscosities from 45 x 10(exp -3) to 290 x 10(exp -3) Pa.s, and inlet flow rates from 0.15 x 10(exp -6) to 0.75 x 10(exp -6) m(exp 3)/s, showed that void content is a function of the capillary number characterizing the flow process. A critical value of capillary number, Ca = 2.5 x 10(exp -3), identifies a zone below which void content increases exponentially with decreasing capillary number. Above this critical value, negligible entrapment of voids is observed. Similar experiments carried out on surface treated nonwoven mats spanning a range of equilibrium contact angles from 20 deg to 78 deg showed that there is a critical contact angle above which negligible entrapment is observed. Below this value, there is no apparent effect of contact angle on the void fraction - capillary number relationship described earlier. Studies on the effect of filament wettability, and fluid velocity and viscosity on the size of the entrapment (voids) were also carried out. These indicate that larger sized entrapments which envelop more than one pore are favored by a low capillary number in comparison to smaller, pore level bubbles. Experiments were carried out on deformed mats - imposing high permeability spots at regular intervals on a background of low permeability. The effect of these spatial fluctuations in heterogeneity of the mat on entrapment is currently being studied.

Mahale, Anant D.↗
NASA NTRS2001

SHIVA-(Spaceflight Holography in a Virtual Apparatus)

SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus) will expand our understanding of the fundamental physics of particle movement in fluids by exploiting the power of holography in a spaceflight experiment'. In addition, the study will exploit the movement of particles in fluids to observe and quantify microgravity phenomena that are extremely important in materials sciences with applications both in space and on earth. The regime under scrutiny is the low Reynolds number, Stokes regime or creeping flow, which covers particles and bubbles moving at very low velocity. The equations describing this important regime have been under development and investigation for over 100 years and yet a complete analytical solution of the general equation had remained elusive yielding only approximations and numerical solutions. In the course of the ongoing NASA NRA, the first analytical solution of the general equation was produced by members of the investigator team using the mathematics of fractional derivatives. This opened the way to an even more insightful and important investigation of the phenomena in microgravity.

Trolinger, James D.↗
DOE OSTI2025

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗
DOE OSTI2024

Testing of a 40-kWth Counterflow Particle-Supercritical Carbon Dioxide Narrow-Channel, Fluidized Bed Heat Exchanger

Particle-based primary heat exchangers (HXs) must deliver sCO2 fluid temperatures above 700°C to couple particle-based concentrating solar receivers and thermal energy storage (TES) sub-systems with efficient sCO2 power cycles. Particle-sCO2 HX designs have struggled to meet DOE cost targets (≤ $150/kWth) due to the amount of expensive nickel alloys necessary for manufacturing full-scale, particle-sCO2 HXs. Our team has demonstrated that mild bubbling fluidization of falling particles in a counterflow narrow-channel fluidized bed can reduce required HX surface area and thus, costs by increasing particle-wall heat transfer coefficients hT,w > 800 W m-2 K-1. This paper reports on the fabrication and testing of a stainless steel, particle-sCO2 HX with 12 fluidized-bed channels approximately 10.5 mm deep spaced between diffusion-bonded, micro-channel sCO2 plates. The HX with a core length of ≈0.56 m is fed with CARBOBEAD HSP particles through a short, fluidized freeboard zone just above the core. Testing to date in the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories has shown that parallel bed fluidization maintains uniform particle inventory across the instrumented channels. Heat transfer thermal duty between the particle and sCO2 flows exceeds 30 kWth with sCO2 inlet temperatures of 200ºC and particle inlet temperatures up to 440ºC and mass flow rates of 0.2 kg s-1 fluidized by counterflowing gas flow rates of 0.005 kg s-1. Tests at higher particle and sCO2 inlet temperatures (600ºC and 400ºC respectively) are targeted to achieve > 40 kWth with model-predicted overall heat transfer coefficients U > 400 W m-2 K-1.

14 SOLAR ENERGY↗
DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
NASA NTRS1991

Spiral vane bioreactor

A spiral vane bioreactor of a perfusion type is described in which a vertical chamber, intended for use in a microgravity condition, has a central rotating filter assembly and has flexible membranes disposed to rotate annularly about the filter assembly. The flexible members have end portions disposed angularly with respect to one another. A fluid replenishment medium is input from a closed loop liquid system to a completely liquid filled chamber containing microcarrier beads, cells and a fluid medium. Output of spent medium is to the closed loop. In the closed loop, the output and input parameters are sensed by sensors. A manifold permits recharging of the nutrients and pH adjustment. Oxygen is supplied and carbon dioxide and bubbles are removed and the system is monitored and controlled by a microprocessor.

Morrison, Dennis R.↗
DOE OSTI2021

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗