Search NASA⌕ Search

SEARCH · Search NASA

Results for “riser flows”

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

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

51 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Development of Methods for In-line Monitoring of Yield Stress During the Transfer of Radioactive Waste - 20435

Radioactive waste slurries at the Savannah River Site (SRS) and Hanford often behave as non-Newtonian fluids with a yield stress. The yield stress increases the size of equipment needed to transport or mix these slurries and can increase the risk of process upsets. This property is often measured by collection of samples and transportation to a laboratory, which can increase the risk of personnel exposure and even change the rheology. In addition, slurry rheology has been observed to change with time. Real-time in-line monitoring (RTIM) of rheological properties of slurries can have substantial benefits to critical processes which are currently in operation or will be part of future plans in Department of Energy (DOE) complex processes. The need for an appropriate RTIM technology that can potentially minimize or eliminate numerous sampling operations was recently identified. This research at Florida International University focuses on addressing the need for a reliable RTIM technology at US Department of Energy sites. The objective of this work is to investigate and develop appropriate and reliable technologies that can perform accurate and easy measurements of yield stress with minimal changes to the flow conditions at US Department of Energy sites. Use of ultrasound in longitudinal and shear modes was considered as a method for monitoring of yield stress variations with no disturbance to the flow. Measurement of shear wave velocity (in shear mode) and monitoring of sound velocity variations with start of flow from rest (longitudinal mode) are the underlying principles. Measurement of pressure loss and liquid rise (due to static differential pressure) during flow of slurries were two other methods that were considered and investigated. In addition, analysis of stability depth of a penetrometer in slurry medium (penetrometer method) revealed a potential suitability. This paper presents efforts associated with investigation of different methods mentioned earlier for kaolin-water mixtures at different concentrations. Theoretical analysis showed that all methods under investigation were capable of monitoring yield stress in a desirable range of 0 to 70 pascal. A bench scale test setup was constructed for investigation of pressure loss and liquid rise methods. Simulants at different concentrations were created for initial testing campaigns. For ultrasonic testing, effects of yield stress on sound propagation speed in longitudinal and shear modes were investigated for simulants placed inside PVC pipes. Parameters such as simulant concentration and pipe diameter were varied to find optimal values. It was found that spacing (or medium thickness) of 1/8 inch or less was necessary to observe echoes with sufficient strength. In addition, a significant effect of yield stress on sound propagation velocity was observed. For the liquid rise method, kaolin-water simulant at 1.17 specific gravity was pumped through a 1/4 inch pipe and liquid height was measured in a vertical riser branching off a three-way connection. Static gauge pressure was varied within [0.3 to 1.1] psi range by increase of pumping pressure and flow. Results showed consistent yield stress data within 10% deviation from values reported by a rheometer. (authors)

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

Loop Thermosyphon Enhanced Solar Collector

A two-phase loop thermosyphon solar collector was developed to transfer thermal energy efficiently and passively from a concentrated solar collector to a thermal desalination process. The two-phase flow in the loop thermosyphon is driven by a difference in gravitational pressure head between the liquid return line and the two-phase riser located between the evaporator and condenser. A numerical model incorporating mass, energy, and pressure drop balances was developed to predict system performance and use as a design tool. Three iterations of loop thermosyphon solar collectors were designed, fabricated, and tested to validate the numerical analysis and demonstrate system integration. Additionally, a novel solar receiver design was evaluated based on volumetric absorption of the concentrated sunlight by an optically absorptive two-phase working fluid. Life tests were performed to evaluate the working fluid thermal stability and resistance to ultraviolet degradation. A glass receiver tube and associated glass-to-metal junctions were developed. Finally, a full-scale loop thermosyphon was integrated with a parabolic trough solar collector and tested on-sun.

14 SOLAR ENERGY↗
DOE OSTI2022

Tank 49H solids disturbance analysis

Tank 49H serves as the feed tank for the Salt Waste Processing Facility (SWPF). Transfers into Tank 49H may disturb any solids that have settled to the bottom of the tank, resulting in feed that may exceed the insoluble solids content limit of 1,200 mg/L of the SWPF Waste Acceptance Criteria (WAC). During a transfer into Tank 49H, material that free falls into Tank 49H through a downcomer could potentially disturb any solids on the bottom of the tank and scour or suspend solids from a settled solids layer or turbid region. A previous analysis and report evaluated the potential to disturb solids when transferring into Tank 49H through the B4 riser and recommended a minimum tank level of 120 inches to prevent disturbing any solids in the bottom of the tank. The scope of this task is to perform additional fluid flow analysis to determine whether accounting for disturbed particle settling and particle mixing and dispersion during transfer could allow the minimum liquid level to be reduced below 120 inches. The analysis utilized models from the technical literature to calculate the size and shape of the “plunging jet” as a function of input parameters such as initial velocity, initial jet diameter, elevation of the initial jet, liquid level in the tank, and solids depth. The analysis relied on the M-Star® simulations performed for the previous analysis to provide bounding estimates of the amount of solid particles disturbed and used the MStar® software to calculate the dispersion and mixing of the disturbed solids with other liquid in the tank as the solids are transported to the transfer pump. The analysis showed that with a solid particle size of 5 micron or less, a liquid level of 120 inches should be maintained to prevent significant disturbance of the solid layer at the bottom of Tank 49H.

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

Preliminary Design of Reactor Cavity Cooling System for a Horizontal Compact HTGR

The Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed MIGHTR in 3 years and support its commercialization as a safe and low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. This report documents the preliminary design study of the RCCS for the HC-HTGR. It includes the establishment of the design requirements, a high-level design study by initial scoping calculations, and preliminary performance calculations of the HC-HTGR RCCS design. Design requirements for the HC-HTGR RCCS have been established to guide preliminary design activities and scoping performance calculations. Initial scoping calculations including estimation of the water inventory, estimation of HVAC thermal capability, and a parametric study on loop dimensions by standalone RCCS analysis. Based on scoping calculation results, a set of baseline dimensions of the HC-HTGR RCCS was derived. A water panel modeling approach was investigated to explore various potential design options for the water panel under consideration for the HC-HTGR RCCS using RELAP5-3D. A test case study was performed to assess the prediction capability of two modeling approaches. The results were compared with CFD simulations conducted in constant RPV temperature and heat flux boundary conditions. It confirms the capability of the RELAP5-3D modeling approach to include all important heat transfer mechanisms expected in the HC-HTGR RCCS operation conditions. Then, a reference RELAP5-3D model for the 1/8 th of a compartment of the preliminary design of the HC-HTGR RCCS was developed. A preliminary performance analysis was conducted to evaluate single-phase natural circulation performance with different top tank temperature values and panel conduction performance in various operation conditions. From single-phase natural circulation performance analysis, the system operation mode was investigated in normal operating and limiting design conditions. It showed operation mode in a subcooled state with a proper top tank water cooling system. Parasitic heat loss by both internal air flow and RCCS was estimated, showing it satisfies maintaining below target maximum heat loss of the HC-HTGR RCCS. From the panel conduction performance analysis, two candidate materials for the riser tube such as carbon steel and stainless steel were compared in the thermal performance of HC-HTGR RCCS. From a single water panel test compared with CFD simulation results, it was confirmed that the current capability of the RELAP5-3D modeling approach for the water panel predicts the thermal conduction of two different materials of the water panel. Then, system-level thermal performance analysis was performed for 1/8 th of the compartment of the preliminary HC-HTGR RCCS design. It was first observed that the current preliminary HC-HTGR RCCS design had minimal impact on the overall thermal performance of the water panel by changing pipe material from carbon steel to stainless steel. From Argonne’s effort on the ongoing water-based NSTF testing program, several considerations other than the thermal performance point of view were addressed to be considered in selecting pipe materials.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2005

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗
NASA NTRS2019

Orion Capsule Parachute Assembly System (CPAS) Overload Testing Approach and Results

The qualification of the Orion Capsule Parachute Assembly System (CPAS) system includes exposure to loads and dynamic pressures above the required values as validation of the parachutes' structural integrity. As outlined in the certification plan, each of the four parachutes of the system are to be subjected to 110% of their respective maximum dynamic pressure requirements. The Main and Drogue parachutes have satisfied this overload condition in drop testing and due to cost and schedule constraints, the Forward Bay Cover Parachute (FBCP) and Pilot parachute were subjected to the overload condition in the ground testing described in this document. The test objectives and pass/fail criteria were established and require the parachutes to achieve and maintain a target riser load (associated with a minimum of 110% dynamic pressure overload) for a minimum of three seconds while sustaining no failures of any structural members (vent hoop, radials, suspension line) or any damage which propagates into catastrophic failure of the canopy. Considering the assumptions and limitations associated with the ground testing (primarily non-uniform flow field of the ground test system and variations in parachute manufacturing), a method of establishing the desired overload condition was determined by the technical community and covers the 2-sigma bounds of the drag area distribution derived from drop testing. On June 27, 2017 the testing was executed at the High Velocity Airflow System (HiVAS) facility located at the Weapon Survivability Laboratory (WSL) at the Naval Air Warfare Center, Weapons Division (NAWCWD) China Lake. Engineering Development Unit (EDU) parachutes were used as pathfinders to gain experience with achieving the test conditions. Additional runs were executed to measure the airflow at the same location as the canopy skirt, although this data is not required to satisfy the test objectives and parachute pass/fail criteria. The qualification parachutes were successfully exposed to the target conditions and sustained only minor damage.

Daum, Jared S.↗
DOE OSTI2020

Validation and application of a multiphase CFD model for hydrodynamics, temperature field and RTD simulation in a pilot-scale biomass pyrolysis vapor phase upgrading reactor

Accurate prediction of transport phenomena is critical for VPU reactor design, optimization, and scale-up. The current study focused on the validation and application of a multiphase CFD model within an open-source code MFiX for hydrodynamics, temperature field, and residence time distribution (RTD) simulation in a non-reacting circulating fluidized bed riser for biomass pyrolysis vapor phase upgrading (VPU). First, an Eulerian-Eulerian approach three-dimensional CFD model was employed to simulate the pilot-scale VPU riser on the supercomputer Joule. Excellent quantitative agreement between experimental and simulated results was achieved for pressure drops and temperature field in a range of operating conditions. Then the validated multiphase CFD model was applied to predict gas and solid residence time distributions (RTDs) since prediction and analysis of RTD is an important tool to study the complex multiphase flow behavior and mixing inside chemical reactors. The predictions show that solid mean residence time is 3.5 times the gas residence time; the solid RTD is more sensitive to the process gas flow rate than the solids circulation rate.

09 BIOMASS FUELS↗
DOE OSTI2025

Steam generator model design parameter sensitivity study for small modular reactor system

Here, this study focuses on design parameter sensitivity studies pertaining to several Once-Through Steam Generator (OTSG) model cases both with and without a riser using python and advanced risk assessment and optimization tool, i.e. Risk Analysis Virtual Environment (RAVEN) developed at Idaho National Laboratory (INL), to support a Small Modular Reactor (SMR) system. The presented Steam Generator (SG) python-based model is a mathematical representation of a steam-generating unit for a Pressurized Water Reactor (PWR)-type SMR system, including fluid flow and heat transfer equations, models, and correlations. Design studies involve changing the model’s input design parameters (e.g., temperature, pressure, mass flow rate) to observe the resulting effects on the output of the system, such as the Heat Transfer Coefficient (HTC), Reynolds number, Nusselt number, and heat transfer performance. Sensitivity studies analyze the degree to which system output and/or desired parameters (e.g., HTC or heat transfer performance) are sensitive to changes in the input parameters. By using RAVEN, detailed design parametric sensitivity studies. Six input parameters—namely, the pressure, temperature, and mass flow rate for the inlet of the primary-side (hot fluid) and secondary-side (cold fluid) of the SG—were randomly perturbed via RAVEN’s Monte Carlo Sampler module, using uniform distributions (i.e., ±1%, ±5% and ±10 % relative changes) for 600 samples. The analysis results give valuable insights into SG system performance, and provide justification for further research and development such as optimized sensor placement, design verification, validation, and optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Blend Pump Mixing in Tank 27 Evaluation

For future waste disposal operations, Savannah River Mission Completion (SRMC) plans to operate Tank 27 as a blend tank. It is desirable to avoid additional worker exposure by using the planned commercial submersible mixing pump (CSMP) installations in Risers B2 and B4 to function as the blend pumps for prolonged blend tank operations. The purpose of this task is to evaluate the impacts to mixing across a range of pump installation heights from 1 to 201 inches. Specifications include: A maximum pump diameter of 22.5 inches, dual nozzles placed tangentially opposed with nominal inside diameter of 2 inches, a total flow range (both nozzles combined) of 500 to 1950 gallons per minute (gpm) and a maximum revolutions per minute (RPM) of 1800. Tank 27 is expected to have maximum fill height of 360 inches.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1959

Flow Studies of Decelerators at Supersonic Speeds

Wind tunnel tests recorded the effect of decelerators on flow at various supersonic speeds. Rigid parachute models were tested for the effects of porosity, shroud length, and number of shrouds. Flexible model parachutes were tested for effects of porosity and conical-shaped canopy. Ribbon dive brakes on a missile-shaped body were tested for effect of tension cable type and ribbon flare type. The final test involved a plastic sphere on riser lines.

Source record↗
DOE OSTI2024

Multi-physics Modeling of Radiative Heat Transfer and Fluid Flow for the Reactor Cavity Cooling System

High-temperature gas-cooled reactors (HTGRs) are notable for their high thermal efficiency and potential for combined heat and power applications. These reactors are particularly appealing due to their advanced passive safety features. HTGRs utilize passive safety systems that function without requiring active components like pumps or compressors during emergencies. These reactor designs depend on a Reactor Cavity Cooling System (RCCS) to manage decay heat removal from the reactor pressure vessel (RPV) during accident conditions. The RCCS consists of vertical rectangular channels known as "risers" or riser ducts positioned around the RPV. These risers receive heat from the RPV through both convective and radiative heat transfer mechanisms. Understanding the interplay of multiple physical phenomena, such as fluid dynamics, heat transfer, and neutron interactions, is essential for the effective design and operation of nuclear reactors, particularly for systems like the RCCS. Multi-physics simulations provide a comprehensive approach to studying these interactions, offering detailed insights and enhancing accuracy. They are especially important in RCCS designs, where the interaction between radiative and convective heat transfer can significantly impact system performance. By leveraging multi-physics simulations, complex reactor behaviors can be modeled without compromising the fidelity of the underlying physical processes. This work aims to establish a robust methodology for coupling multiple physical processes in an air-cooled RCCS. By focusing on validating this multi-physics approach, the study involves designing test cases that simulate various conditions to verify the numerical models employed. The outcomes of this research will provide critical insights for accurately modeling and optimizing complex nuclear systems like the RCCS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

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

FY23 Report on Water NSTF Testing at Two-Phase Conditions: Off-normal Scenarios

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report details the experimental activities and upgrades performed during the program’s fifth year of water-based operation. The theme of this year of testing was examining off-normal scenarios. In practice, that meant examining scenarios with blockages or other design basis conditions. From a maintenance and capability standpoint that meant assessing valving options for implementing blockages, assessing drain and refill capabilities, and ensuring in all cases that operators could safely, efficiently, and repeatably bring the facility to an off-normal condition and return it to the normal operating state. In one notable case, no solution existed, nor was a solution commercially available, that met the facility needs for implementing blockages within the two-phase region of the chimney. A valve and feedthrough were custom-designed for NSTF to operate submerged in saturated water while being safely actuated from outside the tank. The fifth year of testing included eight matrix tests consisting of 162 hours of active heating, 9,735 kWh of electrical heating, with seven tests classified as Accepted per NQA-1 and one as Trending. All matrix tests were performed at two-phase flow conditions. The power parametric series from previous years was extended to include another scenario with a decay load equivalent to 1.75 MWt, full scale, further resolving the system flow oscillation response as a function of input power. A depletion test was performed to extend previous work and examine the system’s approach to stagnation and geysering that were not previously observed. This test began at an initial inventory level of 50%, prototypic power of 2.1 MWt, and employed an accelerated drain of 0.75 gpm plus boiloff until stagnation conditions were achieved and multiple geysering events were recorded. Further exploring the geysering phenomenon was a separate effects test examining static boiling in the risers and chimney. This static boiling test was performed at 48 kWe, a higher power than testing in the previous year, and included enhancements to the drain system to remove inventory in a better-controlled and more reliable manner. Two tests were performed, at prototypic conditions of 1.4 MWt and 2.4 MWt, respectively, examining the system response to throttled conditions in the singlephase region of the piping. A ball valve immediately before he riser inlet header was increasingly throttled, progressing the system through a repeatable pattern of two-phase oscillations and stable regimes before ultimately inducing stagnation and geysering. Concluding the year, two tests were performed at baseline conditions of 70% initial inventory level and 2.1 MWt of prototypic decay heat to examine both repeatability and the response to throttling within the two-phase region, at the inventory tank inlet using the custom valve mentioned above. The first such test repeated transient operating conditions and revealed significant sensitivity t

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Tank 27 Solids Disturbance Testing Using M-Star for the Impact of CSMP Usage on Blend Tank Mixing Effectiveness

Savannah River Mission Completion (SRMC) plans to operate Tank 27as a blend tank for future waste disposal operations. It is desirable to avoid additional worker exposure by using the planned commercial submersible mixing pump (CSMP) installations in Risers B2 and B4 to function as the blend pumps for prolonged blend tank operations. The purpose of this task is to evaluate the impacts to mixing across a range of pump installation heights from 1 to 201 inches above the tank floor. The CSMPs shall maintain adequate blending within a 24 hour maximum timeframe and an adequate nozzle discharge velocity as outlined in U-ESR-G-00030, using the specifications of the planned CSMPs listed in PO SRRA149971 and M-DS-G-00081. Specifications include: maximum pump diameter of D=22.5 in., dual nozzles placed tangentially opposed with nominal inside diameter of 2 in., total flow range (both nozzles combined) of 500 to 1950 GPM and maximum RPM of 1800. Tank 27 is expected to have a fill height maximum of 360". This analysis began by using results from the 2010 testing and applying them to the Tank 27 design to calculate a mixing time as a function of pump flow rate. These results were modified to account for the higher liquid level in Tank 27. Once the liquid level was increased to 360 inches, the results were modified to account for a higher viscosity in Tank 27 than in the testing. Once the viscosity was increased, the effect of changing pump elevation was applied to calculations of the mixing time in Tank 27. Finally, the uncertainty from the testing was included in the analysis to recommend a blend time as a function of pump flow rate and elevation. The conclusion from this study is that the particles in tank 27 will be suspended at both the 600 and 900 gpm flow rates. Only particles above 90 micron will not be suspended in significant amounts when the flow rate is 600 gpm. At 900 gpm even 100 micron particles were thoroughly suspended in the tank. From the series of simulations it has been concluded that it will not be possible to mix tank 27 at the proposed fluid flow rates and not suspend the solid particles at the bottom of the tank.

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

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗