Integration of the commercial blastFoam software and the ASCR-funded MFEM Finite Element Libraries (blastFEM) (Final Report)
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “multiphase CFD”
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.
A CONDUCTION-BASED HEAT PIPE MODEL FOR ANALYZING THE ENTIRE PROCESS OF LIQUID-METAL HEAT PIPE STARTUP
The Heat Pipe-cooled Microreactor (HPM) is one of the micro nuclear reactor designs under active study at the U.S. Idaho National Laboratory. Among the major concerns of HPM research is to understand the startup behavior of the heat pipe-cooled system associated with the startup of the liquid-metal heat pipes initialing from frozen state. The startup of liquid-metal heat pipes typically involves a number of nonlinear mass and heat transport processes including the phase change from solid to liquid and vapor. Hence, it is still a huge challenge to simulate the liquid-metal heat pipe startup using conventional CFD methods and software. The major difficulties of numerical CFD modeling come from the phase-change process, multiphase interaction, microporous wick flow, and compressible gas dynamics that occur during startup of the liquid-metal heat pipes. This paper proposes a simplified conduction-based method to provide practical insights into the entire startup process of the liquid-metal heat pipes while mitigating the challenges of addressing all the complex physics. We discuss the theoretical basis and modeling assumptions to analyze the liquid-metal heat pipe startup from frozen state based solely on heat-conduction equations. Then, the proposed model is implemented into the commercial CFD software to verify the model performance. The model prediction results are discussed via the comparison with the experimental data obtained from sodium heat-pipe startup experiments.
42 ENGINEERING↗
MFIX-Exa: Performance prediction of multiphase energy conversion devices
MFIX-Exa targets the high-fidelity CFD-DEM model in which particles are unresolved by the fluid grid (typically using a fluid mesh approximately twice the particle diameter) but discrete particle dynamics including collisions are resolved with a simple linear-spring dashpot model. The fluid field is solved with a modern low Mach number formation of a cell-centered, nodal-pressure approximate projection method using a Godunov scheme . Physics capabilities include: complex geometries via embedded boundaries (EBs), open and closed system ideal gas equations of state, species transport and heterogeneous chemi- cal reactions. MFIX-Exa is built on the AMReX software framework, the ECP Block-Structured Adaptive Mesh Refinement (AMR) Co-Design Center, which allows the code to be portable and performant.Further integration into the ECP ecosystem includes HYPRE linear solvers and Ascent for in situ visualization. Current project efforts focus on scaling up realistic simulations to the KPP of the challenge problem: NETs 50kW chemical looping reactor, shakedown and scaling on JSLE and OLCF TDS systems and reaching for any remaining performance improvements.
Musser, Jordan↗
Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics
The U.S. Department of Energy National Energy Technology Laboratory's (NETL) 50 kW th chemical looping reactor (CLR) has an underperforming cyclone, which was designed using empirical correlations. To improve the performance of this cyclone using computational fluid dynamics (CFD)-based modeling simulations, four critical design parameters including the vortex tube radius and length, barrel radius, and the inlet width and height were optimized. Here, NETL's open source multiphase flow with interphase exchange (MFiX) CFD code has been used to model a series of cyclones by systematically varying the geometric design parameters. To perform the optimization process, the surrogate modeling and sensitivity analysis followed by the optimization capability in nodeworks was used. The basic methodology for the process is to employ a statistical design of experiments (DOE) method to generate sampling simulations that fill the design space. Corresponding CFD models are then created, executed, and postprocessed. A response surface is created to characterize the relationship between input parameters and the quantities of interest (QoI). Finally, the CFD-surrogate is used by an optimization method to find the optimal design condition based on the objective and constraints prescribed. The resulting optimal cyclone has a larger diameter and longer vortex tube, a larger diameter barrel, and a taller and narrower solids inlet. The improved design has a predicted pressure drop 11 times lower than the original design while reducing the mass loss by a factor of 2.3.
42 ENGINEERING↗
Droplet-Gas Interactions in Nanoelectrospray Multiphase Flow
Electrospray (ES) is production of a charged droplet plume that transits from an electrically biased supply capillary to a counter electrode under the influence of an electric field. Drag interactions with the surrounding gas decelerate the droplets while transferring the momentum from droplets to accelerate the surrounding gas, resulting in an induced gas flow. This work seeks to characterize the structure of the resulting gas jet and identify the key mechanisms defining the flow structure. The phenomenon of gas jets produced by momentum transfer from nano-electrospray (nES) plumes is explored with schlieren visualization, thermal anemometry, and numerical simulations. Schlieren visualization experiments provide information on the flow structure in support of simulation predictions, and the hot thermistor anemometry measurements of gas velocities outside the spray demonstrate quantitatively validated simulation results. The study reveals the formation of a moderately high velocity coaxial gas jet within the nES plume and provides insight into the overall flow structure of the induced flow. The multiphase electrohydrodynamic simulations enable numerical experimentation to explore the fundamental physics of coupled droplet-gas transport and the resulting flow structure. The simulations, confirmed by the experiments, reveal gas jetting induced by nES with a narrow (sub hundred micrometers in diameter) core originating from the nES liquid-jet breakup region and a surrounding larger-in-extent zone (several hundred micrometers in diameter) of lower velocity gas flow within the electrospray plume. This behavior is due to nES ejecting a stream of droplets from the tip of a narrow liquid cone-jet, where the combined effect of many small droplets transferring momentum to a confined region of gas yields a narrow but high-velocity gas stream. Furthermore, the micro-jets produced by nES have practical utility for mass spectrometry, 3D printing and fabrication, and thermal management.
77 NANOSCIENCE AND NANOTECHNOLOGY↗
Transient simulation of biomass combustion in a circulating fluidized bed riser
Interest in circulating fluidized bed (CFB) boilers as a power generation technology has sky-rocketed 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. Herin this paper summarizes the modeling effort on a 50 kW th CFB combustor designed, built, and operated at CanmetENERGY in Ottawa, Canada. The numerical model employs 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 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 both the physics and chemistry of a CFB combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.
09 BIOMASS FUELS↗
Development and Evaluation of a General Drag Model for Gas-Solid Flows via Deep Learning
This project presents the development and evaluation of a general drag model for gas–solid multiphase flows using deep learning techniques. A comprehensive database of more than 4,000 experimental and numerical data points for spherical and non spherical particles was compiled, incorporating geometric features such as sphericity, aspect ratio, and orientation. Several predictive approaches—including traditional em pirical correlations, machine learning, and deep neural networks—were benchmarked, with the proposed Drag Coefficient Correlation-aided Deep Neural Network (DCC DNN) demonstrating superior accuracy. To account for particle–particle interactions, additional drag data were generated using CFD-based simulations of packed and flu idized beds, leading to the development of a retrained model capable of incorporat ing volume fraction effects. Integration of the trained model with the MFiX CFD solver was achieved using FTorch, enabling drag predictions during discrete element method (DEM) simulations. Validation against experimental data for single particles and fluidized beds confirmed the model’s improved predictive ability, particularly for non-spherical geometries. While the model performed strongly under fluidized con ditions, limitations remained in unfluidized regimes, suggesting a need for expanded datasets. Overall, this study demonstrates the feasibility of combining deep learning with physics-informed CFD to improve drag modeling for gas–solid flows, with promis ing implications for scaling multiphase simulations in industrial applications.
42 ENGINEERING↗
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↗
Experimental measurements of gas-particle flows in large-scale strippers
Measurements of gas–solid particle flow in a large-scale stripper unit are reported. The experiments target industrial (pilot) scale measurements of gas-particle flow that can be used to validate numerical methods for modeling multiphase flows, such as computational fluid dynamics coupled to the discrete element method (CFD-DEM). Specifically, experiments were performed with Geldart Group B (533 µm) glass beads in a stripper unit. In general, previous CFD-DEM validation datasets were performed in benchtop experimental systems and limited to O(10 5 ) particles. Here, in this work, the number of particles in the isolated section of the stripper is O(10 7 ) for all tests. Measurements are reported for a ~2 m tall section of the ~1 m diameter stripper. Radial mass flux profiles measured at five axial positions and three axial pressure drop measurements are reported for six operating conditions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations
Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.
42 ENGINEERING↗
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↗
Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture
First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.
Absorption columns↗
Deterministic Calibration of MFiX-PIC, Part 1: Settling Bed
The Particle-in-cell (PIC) numerical approach for modeling granular solids in fluid flow has gained significant interest in recent years. Valued for its often shorter time-to-solution, the PIC formulation relies on modeling statistical groupings of particles called parcels in cooperation with a solids stress model to affect local solids velocity. This is in contrast to the discrete element model (DEM) where every particle in a system is modelled individually and directly coupled to local solids velocity through Newtonian mechanics. The U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) develops and maintains Multiphase Flow with Interphase eXchanges (MFiX), a collection of open-source computational fluid dynamics (CFD) solvers. Included in the MFiX suite are traditional two-fluid model (TFM) and DEM solvers, and a recently added PIC solver (NETL, 2021). In general, PIC methodologies offer an accuracy trade-off in lieu of computational speed; and therefore, it is important to assess the credibility of MFiX-PIC simulations. For this purpose, a systematic verification, validation and uncertainty quantification (VVUQ) effort was initiated at NETL to assess the new PIC solver
42 ENGINEERING↗
Modeling and simulation of multiphase flows
This presentation provides an overview of the National Energy Technology Laboratory’s (NETL) multiphase computational fluid dynamics codes. The highly successful Multiphase Flows with Interphase eXchanges (MFIX) suite has been used to model a wide range of applications including post-combustion carbon capture, bioreactor optimization, and bio-FCC regeneration. MFIX-Exa, a state-of-the-art CFD code, developed under DOE’s Exascale Computing Project, is built on the AMReX software framework (https://amrex-codes.github.io/) and is designed to leverage modern accelerator-based compute architectures. This presentation further reviews the underlying physical models of both MFIX and MFIX-Exa and contrasts their similarities and differences. Examples of past and present CFD simulations will illustrate how scientific computing at NETL is being used not only for scientific exploration but also for design, optimization and scale-up of multiphase flow devices.
Musser, Jordan [NETL]↗
Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification
The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions
12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
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↗
A filter-dependent granular temperature model from large-scale CFD-DEM data
The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.
AMReX↗
Multiscale Experiments and Multiphysics Simulation of Multiphase Flow for Transportable Small Modular Reactors
A new type of safe, small, transportable nuclear reactor would address the intense and ever-growing global demand for energy produced via a resilient, carbon-free energy source. In this regard, transportable small modular reactors (SMRs) are being designed and developed for electricity generation within small/micro-grid/off-grid isolated systems, as well as for heat generation in industrial/residential applications. These reactors feature the capability to be fully factory fabricated and then directly transported to utilities’ sites as “plug-and-play” systems. Research and development (R&D) programs are underway at Idaho National Laboratory (INL) to successfully design, develop, and demonstrate such safe-by-design mobile reactor technologies, in collaboration with partner organizations. Multiscale experimental facilities and multiphysics simulation tools are required for reactor design verification and validation (V&V), and licensing. These advanced reactors are intended to feature passive safety systems such as passive containment cooling systems (PCCS), which consist of multiphase flows and multispecies distributions. This seminar talk will focus on designing and analyzing transportable SMR PCCS by using multiscale experiments and multiphysics computational fluid dynamics (CFD) simulations to support reactor licensing and safety. The corresponding research challenges are addressed via supportive verification and validation results generated by the models and simulation tools in combination with selective parametric and uncertainty analysis. This solution approach could blaze the trail for commercial adoption of such technologies. The facilities, simulation capabilities, and research opportunities available at INL in regard to such reactors and the integrated energy systems with which they go hand in hand are also discussed briefly, and may spark interest in deeper research as well as new collaborative projects.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗