Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase Flow”

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 235 records · Page 13

Stress from long-range interactions in particulate system

In particulate systems, the ensemble average accounting for effects from all particles is expressed as the expected value related to the nearest particle probability. Using this expression, a stress is defined even for long-range particle interactions without the divergence difficulty because of the rapid far field decay of the probability. Similar to the potential part of the virial stress in molecular systems, this stress represents particle-particle interactions through the surrounding field. As examples, the stresses are calculated for systems of charged particles and disperse multiphase flows.

97 MATHEMATICS AND COMPUTING↗

Biomass Combustion in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing resulting in higher combustion efficiency and the ability to use lower quality fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions. SO 2 emissions are conveniently controlled through the addition of Ca-based sulphur sorbents within the combustor. This report summarizes the current modeling effort on a 50 kW th CFB combustor with a diameter of 10 cm and a height of 5 m; designed, built, and operated at CanmetENERGY in Ottawa, Canada employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using 9 kg of olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture the physics and chemistry of a circulating fluidized bed combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

20 FOSSIL-FUELED POWER PLANTS↗

UNESE Data Analysis - Disko Elm Gas Transport Characteristics (NA-22 Quarterly Report)

Two research papers describing gas transport studies in subsurface containment environments appropriate for underground nuclear explosions were revised and accepted for publication in the peer-reviewed journals, Nature – Scientific Reports and Journal of Geophysical Research. The Nature – Scientific Reports paper is a primary deliverable and describes the gas-transport experiment performed in P Tunnel (NNSS) and 1) subsequent analysis of field data as part of the UNESE program. The Journal of Geophysical Research paper describes an experiment, supported mainly by DARPA and led by researchers at Weston Geophysical and New Mexico Tech., with 2) analyses supported by this project to evaluate the ability of SF6 gas tracer to track xenon gas migration under conditions related to a previous UNESE experiment. Two new tasks consider 3) evaluation of detonation heating and multiphase flow on arrival of xenon at the surface and 4) the impact of surface or Langmuir adsorption of gases being transported through pores and along fractures.

58 GEOSCIENCES↗

Computational Fluid Dynamics Analysis of a 12 MW Circulating Fluidized Bed Rise

A three-dimensional (3D) two fluid model (TFM) model of the 12 MW circulating fluidized bed (CFB) riser at Chalmers University was developed using the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The hydrodynamic behavior of the riser was validated by comparing axial pressure profiles from experimental measurements with the results from the simulations. Under some of the operating conditions where there was no external circulation of solids, the pressure profiles from the simulation matched well with the experimental measurements. For the cases where there was external circulation of solids, the pressure profiles with the monodisperse solid phase modeling underpredicted the solids concentration in the freeboard region. With the inclusion of polydispersity in the model, the predicted pressure profiles agreed well with the experimental measurements. A two-dimensional (2D) model of the riser was used to perform simulations of coal combustion. The predicted flue gas outlet compositions of oxygen and carbon dioxide compared reasonably to experiments. The experiments also reported some concentration of carbon monoxide which was overpredicted in the simulations. Additional experimental data such as axial oxygen concentration profile and pressure profiles were also compared with the experiment. It was observed that the oxygen concentration in the bottom bed from the simulations was higher than in the experiments, which suggested that most of the oxidation reactions occurred in the bottom bed. The concentration of volatiles and char in the simulations was lower in the bottom bed than in the experiments, which suggest that there is not sufficient penetration of the solid phase species into the bed.

20 FOSSIL-FUELED POWER PLANTS↗

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗

Full-Loop Simulation of the Combustion of Biomass in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency, and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. This report summarizes the development of a comprehensive computational fluid dynamics (CFD) model for biomass combustion in a 50 kW th CFB combustor employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite at the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL). The experimental combustor was designed, built, and operated at CanmetENERGY, Natural Resources Canada in Ottawa, Canada. This work is the culmination of a multi-year collaboration between NETL and CanmetENERGY. The hydrodynamics in the full loop are validated against cold-flow experiments conducted at CanmetENERGY. To mimic the dense packing in the standpipe in the absence of a collisional model in the MFiX-PIC approach, the ram valve in the return leg is partially closed to obstruct the flow. A combustion reaction scheme was developed that considers the volatiles released by the fuel as a lumped species that fully converts to CO 2 and H 2 O. The so-called enhanced simplified combustion scheme produces an excellent match with experimental data in terms of the species concentrations of CO 2 and O 2 at the outlet and the freeboard, demonstrating the utility of the comprehensive model for a range of flow conditions and operating temperatures. Furthermore, the model can be readily adapted to simulate combustion of a variety of biomass fuels as well as co-fired systems under air and oxy-fuel conditions.

09 BIOMASS FUELS↗

High-Accuracy Simulations to Model Pyrometallurgical Processes in a Secondary Lead Reverberatory Furnace

The US manufacturing industry produces about 1.3 million tons of refined lead each year using secondary sources consisting mainly of lead batteries. ORNL is partnering with Gopher resource, the second largest lead recycling company in the United States, and GTI, to develop a high-fidelity CFD model of a directly fired, reverberatory-style, secondary lead furnace. These High Performance Computing (HPC) simulations are aimed to use first principles modeling for combustion and melting processes of the secondary lead feed while accounting for complex interphase interactions between the gas, solid charge (lead) material, slag, and metal phases. Through validation against operating plant data, this effort will enable significant improvements in design, operational parameters, and energy efficiency, thus improving productivity and refractory lifetime of secondary lead melting furnaces. Estimated savings/reduction of, at least, 1 trillion BTU, 1 million ton/year of greenhouse gas emissions, and $\$50$ million/year to the US lead industry can be expected. ORNL resources and expertise in high-performance computing and multicomponent, multiphase flows were utilized to realize this goal while advancing the understanding of the smelting and melting processes occurring within the furnace.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NRAP-Open-IAM: FutureGen2 Component Models

This report describes the development and testing of three component models for NRAP-Open-IAM, the National Risk Assessment Partnership’s open-source integrated assessment model. The FutureGen2 Lookup Table Reservoir component model is based on interpolation of data from a set of lookup tables. The lookup tables contain pressures and saturations predicted by multiphase flow simulations performed with a reservoir simulator. The FutureGen2 Above Zone Monitoring Interval (AZMI) component is a surrogate model that can be used to estimate the impact that carbon dioxide (CO 2 ) and brine leaks from the CO 2 storage reservoir at the FutureGen 2.0 site might have had on overlying aquifers or monitoring units were a leak to occur. The model estimates the size of “impact plumes” according to five metrics: pH, Total Dissolved Solids (TDS), pressure, dissolved CO 2 and temperature. The FutureGen2 Aquifer component is similar, but is limited to four metrics: pH, Total Dissolved Solids (TDS), pressure, and dissolved CO 2 . The input parameters for each model are the same, but the Aquifer component is applicable to depths between 100 m and 700 m and the AZMI component is applicable from depths between 700 m and 1050 m.

42 ENGINEERING↗

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗

NRAP-Open-IAM: Generic Aquifer Component Development and Testing

The Generic Aquifer Model calculates the concentrations of dissolved salt and dissolved CO 2 surrounding a leaking legacy well. The Generic Aquifer model can also estimate the size of an “impact plume” where concentration changes exceed user-specified thresholds. The model is a component of NRAP-Open-IAM, an open-source Integrated Assessment Model (IAM) developed by the National Risk Assessment Partnership (NRAP) to perform risk assessment for geologic CO 2 storage. The input parameters were selected to cover a wide range of groundwater aquifers and leakage rates. The generic aquifer model was developed using a generative adversarial deep learning network, trained using a large synthetic dataset of STOMP multiphase flow simulations. The deep learning model predictions of dissolved salt and dissolved CO 2 in the aquifer compare well to the original STOMP simulation results. The extent of aquifer impacted by leaking CO 2 or brine is calculated using a user-defined mass fraction threshold. The aquifer impact volumes calculated based on STOMP simulation results compare well to those calculated based on the deep learning model. In a provided python script, gridded observation results from the generic aquifer component of NRAP-Open-IAM are converted to HDF5 format files for monitoring design with the DREAM code.

54 ENVIRONMENTAL SCIENCES↗

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model: Improvement and Quality Assurance

The multisegmented wellbore model (MSW) semi-analytically estimates the amount of CO 2 and brine leakage from a leaking legacy well by segmenting it into intervals to simulate site-specific stratigraphic and hydrogeologic properties. The model is a component of the National Risk Assessment Partnership Open-Source Integrated Assessment Model (NRAP-Open-IAM), which was developed to perform risk assessment for geologic CO 2 storage. The new wellbore leakage model, which uses deep learning networks for a caprock segment, was developed to enhance the analytical MSW. The model was trained and validated using a synthetic data set of Subsurface Transport Over Multiple Phases (STOMP) multiphase flow simulations from various geological, well attribute, and operational conditions to ensure its quality. The results demonstrate that the model is more accurate than the existing model in predicting the transport of two-phase fluids (brine and injected CO 2 ) through the well. This report provides a detailed explanation of the model development and quality assurance.

58 GEOSCIENCES↗

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗

MFiX Development Updates

This presentation discusses recent developments of the Multiphase Flow with Interphase eXchanges (MFiX) software. A brief review of all modeling approaches along with their cost versus accuracy is provided to guide users when selecting a model for a given application. The main new features of the past five official releases of MFiX are described. Improvement in chemistry management allow for easier simulation setup, and faster simulation speed. Progress in the implementation of a thin-wall boundary condition is presented. The major new model development over the past year is the release of the Glued Sphere Particle model (GSP), where component spheres are combined together to represent non-spherical particles. The integration of Machine Learning (ML) workflow in the CFD process is discussed with two applications: a surrogate model for stiff chemistry and the development of a PIC stress model using ML.

Dietiker, Jeff↗

EFRC-MUSE: Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials

Phase interactions and fluid properties in geological and other environments are critical in applications ranging from hydrogen production and geologic storage and recovery, carbon dioxide storage and sequestration, and the sustainable use of water resources. The four goals of EFRC-MUSE: Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials were based on the priority directions articulated in the Basic Research Needs documents, and the scientific needs in nanoscience. 1. Develop a fundamental understanding of confinement and surface interactions in mesoscale media with nanometer-sized pores on the phase behavior, thermodynamic and multiphase flow properties of multicomponent fluid mixtures. 2. Examine the impact of mineralogy and material heterogeneity on mechanical properties to better understand chemo-mechanical interactions in material failure. 3. Determine in-operando cross-scale structural and nanostructural material properties with fluids in confinement and flow under realistic condi

58 GEOSCIENCES↗

Using AMReX's Embedded Boundaries to Support MFIX-Exa’s Geometry Capabilities

MFIX-Exa is an AMReX-based application used for the simulation of chemically reacting multiphase flows inside reactors. The reactor geometry is mostly a combination of primitive shapes such as cylinders, cubes and spheres of varying dimensions connected together with potentially a lot of dead space between them. Additionally, the reactors may have particle-laden jets introduced at different parts of their geometry to enhance reactant mixing. This talk presents how MFIX-Exa uses OpenSCAD’s csg format to allow users to define the embedded boundary (EB) geometry and how it supports specifying inflow through certain sections of the EB. The talk will describe how MFIX-Exa leverages the EB features of grid-pruning, restart files and redistribution. Finally, the talk goes over some capabilities which are currently not available in MFIX-Exa but planned for the future, such as moving parts of the reactor geometry.

embedded boundaries↗

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↗

Nuclear Safety [Vol. 30, No. 3, July-September 1989]

Nuclear Safety is a review journal that covers significant developments in the field of nuclear safety. Its scope includes the analysis and control of hazards associated with nuclear energy, operations involving fissionable materials, and the products of nuclear fission and their effects on the environment. Primary emphasis is on safety in reactor design, construction, and operation; however, the safety aspects of the entire fuel cycle, including fuel fabrication, spent-fuel processing, nuclear waste disposal, handling of radioisotopes, and environmental effects of these operations, are also treated. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 325 Safety of Framatome Advanced Nuclear Steam Supply Systems Designs by J. A. Charles and D. Lange, 333 Book Review of Nuclear Accidents: Intervention Levels for the Protection of the Public by H. B. Piper; ACCIDENT ANALYSIS: 335 Living PRA Computer Systems by S. C. Dinsmore and H.-P. Balfanz, 343 Summary of ICAP Assessments of RELAP5/MOD2 by W. E. Driskell and R. G. Hanson; CONTROL AND INSTRUMENTATION: 352 Thermal Performance Monitoring System at Maanshan Nuclear Power Plant by H.-J. Chao, Y.-P. Lin, G.-H. Jou, L.-Y. Liao, and Y.-B. Chen; DESIGN FEATURES: 358 Warning Systems for Nuclear Power Plant Emergencies by J. H. Sorensen and D. S. Mileti; WASTE AND SPENT FUEL MANAGEMENT: 371 Activities Related to Waste Management Compiled by E G. Silver; OPERATING EXPERIENCES: 382 Steam Generator Tube Performance: Experience with Water-Cooled Nuclear Power Reactors During 1985 by O. S. Tatone and R. L. Tapping, 400 Systems Interaction Analyses: Concepts and Techniques (Part II) by M. D. Muhlheim and G. A. Murphy, 413 Reactor Shutdown Experience Compiled by J. W. Cletcher, 416 Operating U.S. Power Reactors Compiled by E G. Silver; RECENT DEVELOPMENTS: 440 General Administrative Activities Compiled by E G. Silver, 460 Reports, Standards, and Safety Guides by D. S. Queener, 466 Status of Power-Reactor Licensing Activities Compiled by E G. Silver, 470 Proposed Rule Changes as of Mar. 31, 1989; ANNOUNCEMENTS: 334 Proceedings Published, 351 CEC Seminar on Methods and Codes for Assessing the Off-Site Consequences of Nuclear Accidents, 357 Short Course on Multiphase Flow and Heat Transfer: Bases and Applications in A: The Nuclear Power Industry B: The Process Industries, 357 International Conference on Probabilistic Safety Assessment and Management, 478 International Topical Meeting on the Safety, Status, and Future of Non-Commercial Reactors and Irradiation Facilities, 475 The Authors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unveiling the Essential Parameters Driving Mineral Reactions during CO2 Storage in Carbonate Aquifers through Proxy Models

Numerical simulation is a commonly employed technique for studying carbon dioxide (CO2) storage processes in porous media, particularly saline aquifers. It enables the representation of diverse trapping mechanisms and the assessment of CO2 retention capacity within the subsurface. The intricate physicochemical phenomena involved necessitate the incorporation of multiphase flow, accurate depiction of fluid and rock properties, and their interactions. Among these factors, geochemical reaction rates and mechanisms are pivotal for successful CO2 trapping in carbonate reactive rocks. However, research on kinetic parameters and the influence of lithology on CO2 storage remains limited. This limitation is partly due to the challenges faced in laboratory experiments, where the time scale of the reactions and the lack of in situ conditions hinder accurate measurement of mineral reaction rates. This study employs proxy models constructed using response surfaces calibrated with simulation results to address uncertainties associated with geochemical reactions. Monte Carlo simulation is utilized to explore a broader range of parameters and identify influential factors affecting CO2 mineralization. The findings indicate that an open database containing kinetic parameters can support uncertainty assessment. Additionally, the proxy models effectively represent objective functions related to CO2 injectivity and mineralization, with calcite dissolution playing a predominant role. pH, calcite concentration, and CO2 injection rate significantly impact dolomite precipitation, while quartz content remains unaffected.

58 GEOSCIENCES↗