Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

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 559 records · Page 31

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Effects of Different Types of Entrances on Natural Ventilation in a Subway Station

In this study, the natural ventilation of a horizontal entrance of a typical subway station is investigated based on numerical simulations and experiments. In addition, a renormalization group k-e model (RNG k-e model) is applied to compute both the internal and external airflow patterns. Computational fluid dynamics (CFD) simulations are validated based on the experimental results. Furthermore, a multiple variable regression model is employed to study how the different parameters affect the internal airflow rates of the subway station statistically, according to the experimental results. For this typical model, the parameter importance can be ranked as follows: (1) outdoor wind speed; (2) flow resistance of the subway station; (3) height of the wind catcher; and (4) length of the wind catcher. To understand the detailed pressure distribution of the horizontal entrance of the subway station with and without the wind catcher, 3D numerical simulations are conducted for different scenarios. We attempt to alter the size of the wind catcher (including the length and height) to study the characteristics of the pressure on the horizontal entrance of the subway station under outdoor wind-driven conditions. The pressure distributions for the entrance for different scenarios are compared and analyzed. Finally, the interactions between the internal and external flows are investigated by changing the resistance of the subway station. When the internal flow resistance is changed, the pressure coefficient (Cp) of the entrance is different as well. Hence, the Cp is not only affected by the outdoor environment, but is also influenced by the internal airflows.

CFD↗

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗

Evaluation of Fish-Related Properties of Kaplan Turbines at the Design Phase: Simulation-based outcomes vs. Experimental Data

The development of turbine technology faces growing demands to maximize the survival of migratory fish passing through turbine flows. In this context, two strategies have emerged for quantifying hazardous hydraulic conditions: computer-based evaluations at design stage and recordings with autonomous sensors deployed in prototypes. The former is a desktop evaluation with many modelling assumptions (idealization) and the latter is a field technique that introduces various unknown and uncontrollable factors (uncertainty.) The present work introduces and implements a third method based on test rig measurements and the corresponding computer-based predictions of conditions that negatively affect fish survivability. The experimental work was conducted in a five-bladed Kaplan turbine model in which miniaturized autonomous sensors (SF Mini, developed at the Pacific Northwest National Laboratory, U.S.) measured fish-relevant hydraulic features. The modelling work involved flow simulations according to industry practices and the representation of fish trajectories through the simulated flow conditions. We compared both the experimental measurements and CFD outcomes and discussed the challenges and advantages of the modelling strategies, as well as the benefits for turbine engineers in need of incorporating effective design concepts to mitigate fish mortality through turbines.

Romero-Gomez, Pedro↗

UPDATES FROM THE INVOLUTE WORKING GROUP

The HFIR, RHF, and FRM II reactors represent a particular class of Research and Test Reactors that provide some of the most intense and continuous neutron fluxes for science, industry, and medical applications. These high-performance reactors have achieved compact cores by operating with Highly Enriched Uranium fuel (HEU, 235U/U ≥ 20 wt. %) and utilizing fuel plates curved as an involute. Due to the proliferation risks, the international community aims to reduce or eliminate, when possible, the use of HEU fuel in civilian facilities by converting them to a Low-Enriched Uranium fuel (LEU, 235U/U < 20 wt. %). Conversion of these reactors without significantly compromising their performance or safety is a challenging endeavor that can tremendously benefit from advanced computational tools and thus, eliminate unnecessary conservatism to ensure sufficient thermal margins. Therefore, models are being developed using modern Computational Fluid Dynamics (CFD) and Computational Structural Mechanics (CSM) software to evaluate the steady-state safety margins of various LEU designs instead of being reliant on the more traditional, conservative methods. To gain the confidence and acceptance of high-fidelity modeling by the nuclear regulators, Argonne National Laboratory (ANL) and the involute reactors have formed an informal scientific group, the Involute Working Group (IWG). The IWG facilitates inter-organizational collaboration on experimental benchmarking, code-to-code comparisons, and Verification and Validation (V&V). This paper describes some of the recent IWG efforts in validating software against the existing experimental data, as well as code-to-code comparisons of different software used by the IWG members.

Bergeron, Aurelien↗

Multi-Physics Simulations of Molten Chloride Fast Reactor using Nek5000 and PROTEUS-NODAL

This report documents the FY20 work to demonstrate a specific multi-physics capability within NEAMS tools to model molten chloride fast reactors (MCFR). Specifically, the CFD code Nek5000 and the neutronics code PROTEUS-NODAL were coupled off-line to calculate the delayed neutron precursor distribution within the flowing fuel salts in an MCFR core. An MCFR benchmark problem and a 3D test problem based on a conceptual MSR design were successfully solved through the coupled simulation of Nek5000 and PROTEUS-NODAL. In the coupled simulations for these test problems, Nek5000 provided high-quality flow field information, and PROTEUS-NODAL performed the neutron diffusion calculation along with the precursor drift model. The impact of precursor drift on the core status was also investigated. Additionally, the precursor drift model implemented in PROTEUS-NODAL was verified by making use of the reference precursor distribution generated from Nek5000. Results of these test problems suggest that the multi-physics capability of Nek5000 and PROTEUS-NODAL was properly developed for tracking delayed neutron precursor drift.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation of RANS Models for Predicting Bypass Flow in Prismatic VHTRs

In prismatic Very-High Temperature Reactors (VHTRs), gaps occur between neighboring fuel assemblies referred to as bypass gaps. This work aims to measure the flow distribution between the bypass gaps between neighboring assemblies and the coolant channels within the assemblies. Computational fluid dynamics (CFD) was done using Reynolds-Averaged Navier Stokes (RANS) modeling in STAR-CCM+ to assess the ability of RANS models to accurately capture this phenomena. The RANS results were compared to experimental worked performed by Seoul National University (SNU) and Korea Atomic Energy Research Institute (KAERI) for validation. The CFD results showed good comparison with the experimental data available from SNU and flow split was characterized for multiple configurations, inlet mass flow rates, and bypass gap sizes.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn

The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.

42 ENGINEERING↗

Simulations of Bypass Flow in Prismatic VHTRs

In prismatic Very-High Temperature Reactors (VHTRs), gaps occur between neighboring fuel assemblies referred to as bypass gaps. This work aims to measure the flow distribution between the bypass gaps between neighboring assemblies and the coolant channels within the assemblies. Computational fluid dynamics (CFD) was done using Reynolds-Averaged Navier Stokes (RANS) modeling in STAR-CCM+ to assess the ability of RANS models to accurately capture this phenomena. The RANS results were compared to experimental worked performed by Seoul National University (SNU) and Korea Atomic Energy Research Institute (KAERI) for validation. The CFD results showed good comparison with the experimental data available from SNU and flow split was characterized for multiple configurations, inlet mass flow rates, and bypass gap sizes.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

TRANSIENT THERMAL MODELING OF THE HIGH BURNUP DEMONSTRATION RESEARCH PROJECT CASK USING STAR-CCM+ AND COBRA-SFS

The Department of Energy in collaboration with the Electric Power Research Institute is in the process of conducting the High Burnup Demonstration Project. Where the objective is to characterize the performance of high-burnup fuel in long term storage. As part of this demonstration, a TN-32B dry storage cask was instrumented and loaded with spent nuclear fuel at North Anna Nuclear Generating Station in November 2017. The project cask provides a unique opportunity to gain information on spent fuel and cask performance from an in-service operational system. The cask was instrumented with thermocouples inside of the fuel assemblies, then loaded and dried using normal procedures. After the drying process a thermal soak period where the cask was left indoors was used to obtain steady temperatures for model comparison along with surface temperatures. This paper details thermal modeling validation work that was done to model both the steady state and transient cases. Two modeling tools were used to predict temperatures in the cask. The general purpose CFD and heat transfer code STAR-CCM+ was used with both a detailed pin-by-pin model was used along with a more efficient k-effective simplification. The other code is COBRA-SFS, a purpose-built detailed thermal modeling tool developed and maintained at PNNL. Results from all the modeling tools were compared blind to the test data. Each model utilized general design information and compared reasonably well to the blind steady state data using convection and surface temperature boundary conditions. Post-test adjustments were later made to better reflect the “as built” conditions of the cask. Based on lessons learned from the steady state analysis the vacuum drying process was modeled in both codes. Utilizing the measured loading conditions the transient models were able to compare very well with the measured data. Overall the project showed an ability to model spent fuel storage conditions very well and future work is planned to generalize the methodology used for vacuum drying.

Thermal Analysis, Spent Fuel, Dry Storage↗

Modeling cell venting and gas-phase reactions in 18650 lithium ion batteries during thermal runaway

A numerical model is developed to study cell venting, internal pressure, and gas-phase dynamics behavior of 18650 Li-ion cells undergoing thermal runaway. A k-e Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer's equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. Then, a series of computational fluid dynamics (CFD) simulations are conducted on a single 18650 cell at various states-of-charge (100%, 50%, 25%) to study detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is found that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

25 ENERGY STORAGE↗

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗

Impact of Different Thermal Gradients on the Dynamics of Cylindrical Lithium-ion Cells Subject to Accelerated Aging and on Module Performance

This study investigates the impacts of applying different thermal gradient patterns to cylindrical lithium-ion cells in a module on cell dynamics (temperatures, current flows, state of charge), module performance (evolution of resistance, capacity, and energy versus cycle number), and module lifetime. The thermal gradients were generated using cooling plates (CPs) with three different flow-field designs, namely, straight, perpendicular, and U-turn. The study uses computational fluid dynamics (CFD), the pseudo-two-dimensional (P2D) battery model, capacity loss and increased impedance due to the growth of a solid-electrolyte-interphase, and the electric current distribution from module terminals to cells that depends on the series-parallel electrical connections among the cells. The impact of the thermal gradient (resulting from the CP designs) on the variability in resistance, current, state of charge, and voltage among the cells was analyzed and linked to differences in the module's performance. Applying a thermal gradient to parallel-connected strings of series-connected cells led to variation in the current through each parallel string and an imbalance in the voltage of series-connected cells. Module performance is poorer when the thermal gradient causes a voltage imbalance than when it causes a current imbalance. Module performance becomes the worst when both current variation and voltage imbalance happen together. For instance, the module's lifetime (estimated as reaching 80% of its initial capacity) varied by 5% to 17.5%, depending on the magnitude and pattern of the imposed thermal gradient. As the relative orientation between thermal gradients and cells' electrical connectivity influences the module's performance, appropriate consideration should be given to the choice of the CP, especially if large thermal gradients are allowed.

Battery thermal management↗

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↗

Characterizing Impacts of Dry Coal Feeding in High Pressure Oxy-Coal Combustion Systems

Reaction Engineering International (REI) has managed a team of experts from the University of Utah, Southeast University (SEU) in Nanjing, China, Electric Power Research Institute (EPRI), Corrosion Management Ltd. (C-M), Praxair, and Brigham Young University (BYU) to investigate dry pulverized coal feeding for pressurized oxy-coal combustion. Dry feed firing systems for entrained flow, pressurized, oxy-coal combustors have not been well developed, although related technologies have been used in the Shell Gasification Process and for pressurized fluid bed combustion. DOE-funded research recently completed at REI and the University of Utah focused on characterizing impacts of high temperatures and pressures in oxy-coal combustion systems. For high pressure combustion, that research used a coal slurry feed into a 17 bar pressurized combustor. As a consequence of that research, it was identified that fuel feeding and firing system flexibility are challenges that require attention. Based on that experience, the approach of using a coal slurry feed system leads to challenges in producing consistent atomization of the slurry, which causes burnout problems, especially at high pressures. In addition, slurry atomization processes may be difficult to scale to sizes appropriate for practical commercial use. Dry pressurized coal burner systems, on the other hand, have the potential to yield efficiency gains, improve flexibility and facilitate applications at larger scales. Experimental work was conducted at the University of Utah Industrial Combustion and Gasification Research Facility as well as the 100 kW pressurized oxy-coal combustor (POC) facility at Brigham Young University. Mechanism development and CFD-based combustion and dense-phase flow modeling were performed at REI. Successful completion of the project objectives has resulted in the following key deliverables: 1) Design and prototype of a pressurized pulverized coal feeding and oxy-firing system 2) Data from a 100kW, pressurized (15 bar) entrained flow reactor with a dry feeding delivery and burner system that describes flame characteristics, radiative heat flux profiles, carbon burnout, along with characteristics of ash aerosols, fouling, and slagging. 3) Validated and transportable models that describe the relevant conditions in pressurized oxy-combustion systems and that can be used for scale-up and optimization. 4) Principles to guide design of high pressure, pilot-scale and full-scale coal oxy-firing systems. 5) Assessment of pressurized oxy-combustion impacts on key parameters relevant to oxy-coal fired utility boilers such as coal devolatilization, char oxidation, mineral matter transformation, deposition, and corrosion. The experimental data, pressurized oxy-firing system principles, and process mechanisms provided by this work can be used by electric utilities, boiler OEMs, equipment suppliers, design firms, software vendors, consultants and government agencies to assess the use of high temperature and high pressure oxy-combustion in current research and to guide development of new oxy-coal boiler designs.

01 COAL, LIGNITE, AND PEAT↗

Characterizing Impacts of Dry Coal Feeding in High Pressure Oxy-Coal Combustion Systems

Reaction Engineering International (REI) has managed a team of experts from the University of Utah, Southeast University (SEU) in Nanjing, China, Electric Power Research Institute (EPRI), Corrosion Management Ltd. (C-M), Praxair, and Brigham Young University (BYU) to investigate dry pulverized coal feeding for pressurized oxy-coal combustion. Dry feed firing systems for entrained flow, pressurized, oxy-coal combustors have not been well developed, although related technologies have been used in the Shell Gasification Process and for pressurized fluid bed combustion. DOE-funded research recently completed at REI and the University of Utah focused on characterizing impacts of high temperatures and pressures in oxy-coal combustion systems. For high pressure combustion, that research used a coal slurry feed into a 17 bar pressurized combustor. As a consequence of that research, it was identified that fuel feeding and firing system flexibility are challenges that require attention. Based on that experience, the approach of using a coal slurry feed system leads to challenges in producing consistent atomization of the slurry, which causes burnout problems, especially at high pressures. In addition, slurry atomization processes may be difficult to scale to sizes appropriate for practical commercial use. Dry pressurized coal burner systems, on the other hand, have the potential to yield efficiency gains, improve flexibility and facilitate applications at larger scales. Experimental work was conducted at the University of Utah Industrial Combustion and Gasification Research Facility as well as the 100 kW pressurized oxy-coal combustor (POC) facility at Brigham Young University. Mechanism development and CFD-based combustion and dense-phase flow modeling were performed at REI. Successful completion of the project objectives has resulted in the following key deliverables: 1) Design and prototype of a pressurized pulverized coal feeding and oxy-firing system 2) Data from a 100kW, pressurized (15 bar) entrained flow reactor with a dry feeding delivery and burner system that describes flame characteristics, radiative heat flux profiles, carbon burnout, along with characteristics of ash aerosols, fouling, and slagging. 3) Validated and transportable models that describe the relevant conditions in pressurized oxy-combustion systems and that can be used for scale-up and optimization. 4) Principles to guide design of high pressure, pilot-scale and full-scale coal oxy-firing systems. 5) Assessment of pressurized oxy-combustion impacts on key parameters relevant to oxy-coal fired utility boilers such as coal devolatilization, char oxidation, mineral matter transformation, deposition, and corrosion. The experimental data, pressurized oxy-firing system principles, and process mechanisms provided by this work can be used by electric utilities, boiler OEMs, equipment suppliers, design firms, software vendors, consultants and government agencies to assess the use of high temperature and high pressure oxy-combustion in current research and to guide development of new oxy-coal boiler designs.

pressurized oxy-coal combustion, sub-micron ash ae↗