Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid bubbles”

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

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Investigating biomass composition and size effects on fast pyrolysis using global sensitivity analysis and CFD simulations

It is notoriously difficult to build an accurate universal model for biomass pyrolysis due to its sensitivity to a wide number of critical material attributes such as chemical species and physical sizes. In this work, a biomass pyrolysis kinetics with 32 heterogeneous reactions and 59 species was implemented in an open-source multiphase computational fluid dynamics (CFD) software MFiX and validated against two different experimental pyrolysis data sets that provided detailed data describing chemical component yields. The reaction scheme was then used to build a surrogate model and assess the sensitivity of pyrolysis yields to feedstock compositions. The sensitivity analysis determined that the yield of bio-char showed a strong positive sensitivity to the carbon-rich lignin and tannin pseudo-species in the reaction scheme while the bio-oil and bio-gas were correlated to oxygen-rich lignin pseudo-species. The reaction scheme was then integrated into a coarse-grained discrete element model to simulate fast pyrolysis in a bubbling fluidized bed over a range of feedstock particle sizes. The reactor simulations showed further sensitivity to particle size and hydrodynamics. Notably, particles under 0.5 mm have small heat transfer limitations but left the reactor before completely converting and thus reduced the bio-oil yield. Results from this study can be used to guide future development of highly accurate models for fast pyrolysis reactors with a variety of feedstock properties and operating conditions.

09 BIOMASS FUELS↗
DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Effects of variable deceleration periods on Rayleigh-Taylor instability with acceleration reversals

The dynamics of an interfacial flow that is initially Rayleigh-Taylor unstable but becomes statically stable for some intermediate period due to the reversal of the externally imposed acceleration field is studied. We discuss scenarios that consider both single and double-acceleration reversals. The accel-decel (AD) case consists of a single reversal imposed at an instant after the constant acceleration instability has entered a self-similar regime. The layer of mixed fluid ceases to grow upon acceleration reversal, and the dominant mechanics are due to internal wave oscillations. Variation of mass flux and the Reynolds stress anisotropy is observed due to the action of the internal waves. Here, a second reversal of the AD case that is termed as accel-decel-accel, ADA is then explored; the response of the mixing layer is shown to depend strongly on the duration and the periodicity of the Reynolds stress anisotropy of the mixing layer during the deceleration period. We explore the effect of this variable deceleration period after the second acceleration reversal where the flow once again becomes Rayleigh-Taylor unstable based on metrics that include the integral mixing-layer width, bubble and spike amplitudes, mass flux, Reynolds stress anisotropy tensor, and the molecular mixing parameter.

42 ENGINEERING↗
DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Experimental and Numerical Investigation of Flow Boiling in Additive Manufactured Foam Structures With Vapor Pathways

Abstract The unique properties of metal foams make them potential candidates for a range of applications, including microsystem thermal management. Using additive manufacturing to create foam-type structures can improve upon prior thermal solutions by eliminating thermal interface materials and allowing for customization/local control of parameters. In the present investigation, flow boiling in additive-manufactured metal foams is investigated both experimentally and numerically. Two test samples, one with uniform structure and the other with pathways for vapor removal, are compared both experimentally and numerically. A conjugate computational fluid dynamics and heat transfer (CFD-HT) model utilizing a three-dimensional volume of fluid (VOF) model with accompanying evaporation/condensation model provided in-depth visualization of the boiling flow phenomena. The experiments generated the thermohydraulic performance over a range of heat fluxes, demonstrating that the sample incorporating dedicated vapor pathways performed better in both pressure and heat transfer performance metrics compared to the uniform foam. Additionally, negative system-level effects (i.e., hydraulic oscillations) were shown to be abated using the vapor removal structures. The numerical model yielded further insight into the factors contributing to the improved performance. Results indicated the pathways functioned as vapor removal channels, allowing the generated vapor to vent from the foam structure into the lanes. Further computational investigations demonstrated changes in flow regimes, where the addition of vapor channels caused the flow to change from churn to annular. Bubble behavior unique to the vapor pathway structure was studied, showing stagnant regions that eject vapor into the channel.

Engineering↗
DOE OSTI2024

MOSCATO Development and Integration in Fiscal Year 2024

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy evolution for molten salt equipment. In FY24, significant developments were made to the code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: 1. Implementation of tritium transport capabilities and validation with experimental data: To enable modeling of tritium and other fission gases within MSRs, we implemented gas transport within MOSCATO via inclusion of couple mass transport equations within the salt and structural alloys. Comparisons to experimental data from literature showed good agreement with respect to tritium release rates. 2. Preliminary implementation of two-phase flow models in MOSCATO: To model tritium and other gases above their solubility limits, we implemented preliminary two-phase flow models within MOSCATO to account for bubble transport. The first model adopted was the Level-Set approach, which can handle the high void fraction regime, but with a requirement for high mesh resolution thus high computational expense. In this report, we present a verification of the Level-Set method using a simple benchmark case. We also performed a demonstration of the code as applied to an experimental case involving cover gas flow through salt in an experimental vessel. The second model adopted was the Eulerian-Eulerian dispersed flow model, which is computationally cheaper but limited to low void fraction regimes, such as bubbly flow. Validation and verification have not yet been performed for the Eulerian-Eulerian approach, but a preliminary implementation was completed. 3. Validation with static corrosion experiments: Static corrosion experimental data for stainless steel coupons within molten salts was used to further validate the corrosion model in MOSCATO. To do so, we leveraged the existing models in MOSCATO and simulated the sample mass loss and mass gain phenomena. Several ion species, including Cr 2+ , Fe 2+ and H + , were simulated in salt using the PNP solver, while Cr 0 and Fe 0 were simulated with a diffusion solver in stainless steel. The mass loss of the samples was compared with experimental data, and good agreement was achieved. These combined activities served to further expand the capabilities of MOSCATO and make it more generally applicable to the full range of phenomena that can control chemistry and corrosion in molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Comprehensive experimental investigation on biomass-glass beads binary fluidization: A data set for CFD model validation

Fluidization behavior of biomass and glass beads binary mixtures in a bubbling fluidized bed was experimentally investigated. Mixtures containing different mass fraction of Loblolly Pine white wood and glass beads were fluidized at different fluidization velocities. The particle properties were characterized in a QICPIC that uses a dynamic image processing method to measure both particle size and sphericity. The minimum fluidization velocity was determined using the pressure drop method. An image processing method was developed to capture the dynamic expanded bed height at a very high frequency. The effect of biomass mass fraction and inlet gas velocity on mixing and segregation behavior was studied and analyzed through pressure drop measurements. Pressure drop fluctuations and expanded bed height fluctuations via fast Fourier transform were analyzed and compared. Here, the complete and accurate experimental data reported in this study could provide a benchmark data set for various computational fluid dynamics models validation, calibration, and identification.

42 ENGINEERING↗
DOE OSTI2025

Mitigation of Boiling-Induced Thermal Degradation Using Microporous Nickel Inverse Opals Structures

Engineered microporous structures have received much attention in high-heat-flux electronics cooling due to their high thermal conductivity and permeability, and large surface area for heat transfer, but are susceptible to boiling-induced thermal degradation. Here, this study investigates the efficacy of nickel inverse opals (NiIOs) in mitigating structural degradation caused by corrosion-assisted erosion during pool boiling with water as the working fluid. First, we compared the reliability of NiIOs to copper inverse opals (CuIOs) for a 3-day pool boiling test at constant heat flux. The NiIOs demonstrated superior resistance to thermal degradation due to their inherent corrosion resistance and mechanical strength. Subsequently, we conducted a more controlled experiment to show the effect of heat flux on the degradation of the NiIOs while excluding the effect of temperature variations. Pool boiling tests of 20-μm-thickness NiIOs covering an area of ∼11 × 11 mm 2 with a 2.5 × 2.5 mm 2 heater at the center were conducted at heat flux levels of 20%, 40%, and 60% of the critical heat flux (CHF) for 3 days. The NiIOs subjected to heat flux levels of 20% and 40% CHF showed minimal degradation, while the sample subjected to 60% CHF showed erosion on the top surface due to higher bubble formation and departure rate. These results show the potential of NiIOs as a promising solution for long-term thermal management in high-power electronic devices, although design considerations for maximum allowable heat flux are necessary for reliable operation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Demonstration of a multi-channel fluidized bed particle–supercritical carbon dioxide heat exchanger for concentrating solar applications

High-temperature thermal energy storage in oxide particles at temperatures above 600°C can couple concentrated solar energy with high-efficiency thermal power cycles to provide dispatchable solar-driven electricity. Challenges remain in developing cost-effective primary heat exchangers, which require expensive alloys, to extract the high-temperature thermal energy from the particles to power cycle fluids, such as supercritical CO 2 (sCO 2 ) in recuperated Brayton cycles. To explore one pathway for cost-effective, high-temperature particle heat exchangers, the current study demonstrates a shell-and-plate, particle–sCO 2 heat exchanger with narrow- channel fluidized beds coupled with micro-channel sCO 2 flows in the heat exchanger walls. This study evaluates the feasibility of multiple parallel, narrow-channel fluidized beds in shell-and-plate particle–sCO 2 HXs, to achieve high bed-wall heat fluxes at elevated temperatures. A reduced-order model simulates the narrow- channel, fluidized-bed particle–sCO 2 heat exchanger to design the fluidized bed geometry, in terms of depth, height, and number of channels,for a nominal 40-kWth heat exchanger at particle and sCO 2 inlet temperatures up to 600 °C and 400 °C respectively. The resulting shell-and-plate heat exchanger design operates with bubbling fluidization of the downward-flowing oxide particles to enhance bed-wall heat transfer. The heat exchanger core is fabricated with etched sCO 2 micro-channels in thin wall plates that are diffusion bonded to spacer frames to form the shell-and-plate structure with 12 parallel, fluidized bed channels, 10.4 mm deep. The heat exchanger is tested at the National Solar Thermal Test Facility at Sandia National Laboratories with CARBOBEAD HSP particles at design particle flow rates of 0.20 kg s –1 and inlet temperatures up to 525 °C. Results show that fluidization across multiple parallel channel beds can maintain uniform particle inventory with a common freeboard zone above the heat exchanger core. Bubbling fluidization improves particle–wall heat transfer coefficients but also increases axial dispersion of particle thermal energy, which lowers the log- mean temperature difference such that total heat transfer remains relatively constant to within ±10% over a broad range of fluidization gas velocities. The axial dispersion required particle and sCO 2 flow rates to be increased by 25% over model-designed conditions to achieve the targeted 40 kWth, which indicates the importance of incorporating axial dispersion into heat exchanger design models and of deploying bed structures to suppress it. Furthermore, this study demonstrates the feasibility and preferred fluidizing gas conditions for particle heat exchangers for releasing high-temperature thermal energy storage systems.

14 SOLAR ENERGY↗
DOE OSTI2021

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2019

Retained Hydrogen Literature Study for DWPF

During Defense Waste Processing Facility (DWPF) operation in the Slurry Receipt and Adjustment Tank (SRAT) and the Slurry Mix Evaporator (SME), hydrogen is primarily produced by radiolysis and catalytic generation. The hydrogen is released from the sludge by agitation and removed from the vapor space by purging with air. Upon loss of agitation, the hydrogen begins accumulating in the sludge, and releases rapidly upon the resumption of agitation. If enough hydrogen accumulated in the sludge, and released rapidly upon the start of agitation, the hydrogen concentration in the vapor space could exceed the lower flammability limit (4 vol %). Savannah River Remediation (SRR) currently manages the risk for retained hydrogen in the DWPF vessels by implementing a Retained Hydrogen Program. The program primarily relies on the timely operation of vessel agitators to ensure that the retained hydrogen eruptions that could lead to flammable vapor spaces do not occur. The current program relies on conservative assumptions concerning gas release and retention (e.g., instantaneous release, conservative hydrogen bubble fraction, no release during quiescent periods) to develop allowable vessel quiescent times (Q-times). SRR would like to relax the conservative assumptions to the Retained Hydrogen Program and allow the Retained Hydrogen Program to extend past Sludge Batch 9. Relaxing the assumptions could potentially allow longer Q-times during agitator stoppages. SRR requested Savannah River National Laboratory (SRNL) personnel to perform a literature study to extract industry data pertaining to retained hydrogen and apply it to the DWPF process. The authors conducted a review of the technical literature, including work from the Savannah River Site (SRS), SRNL, the Hanford site, Pacific Northwest National Laboratory (PNNL), and the Sellafield site to identify literature describing the retention and release of gases with non-Newtonian slurries. The study was undertaken with the knowledge that SRS wastes are unique which makes using non-SRS data challenging. The literature search focused on the following topics: Mixing or blend time in Newtonian and non-Newtonian fluids, Quiescent release of gas from an unagitated non-Newtonian fluid, Release of gas from an agitated non-Newtonian fluid, and Changes in fluid rheology as a function of settling time.

08 HYDROGEN↗
DOE OSTI2022

Very small-scale, segregating-fluidized-bed experiments: A dataset for CFD-DEM validation and uncertainty quantification

IWe report fluidization experiments were conducted on a small scale and with a rapid response (short duration) to enable corresponding simulations at low-computational cost. Rise times are reported for four or fewer polyethylene particles (intruders) in an air-fluidized bed of ~5000 group D glass beads. Experimental inputs were completely characterized—particle properties, system dimensions and operating conditions—which is necessary for validating computational fluid mechanics (CFD)-discrete element method (DEM) including a comprehensive uncertainty quantification (UQ) analysis. Input uncertainties are reported as bounds or cumulative distribution functions of measured values. The staggering number of simulations required to complete a UQ analysis (~O[10 4 ] simulations corresponding to ~5 uncertain inputs) motivates this study. These segregating-bed experiments are designed to permit analogous CFD-DEM simulations to complete in less than a day on a single (~2.5 GHz) computational processor unit (CPU). Segregation times are reported for several operating conditions, intruder sizes, and initial configurations, providing a rich dataset for numerical model testing, validation and UQ.

42 ENGINEERING↗
DOE OSTI2022

Speed of Sound Measurements of Binary Mixtures of Hydrofluorocarbons [Pentafluoroethane (R-125), 1,1-Difluoroethane (R-152a), or 1,1,1,2,3,3,3-Heptafluoropropane (R-227ea)] with Hydrofluoroolefins [2,3,3,3-Tetrafluoropropene (R-1234yf) or trans-1,3,3,3-Tetrafluoropropene (R-1234ze(E))]

In this article, speed of sound data measured using a dual-path pulse-echo instrument are reported for three binary refrigerant mixtures, R-125/1234yf, R-1234yf/152a, and R-1234ze(E)/227ea, at compositions of (0.33/0.67) and (0.67/0.33) mole fraction. The speed of sound was studied at temperatures ranging from 230 K to 345 K from pressures slightly above the bubble point curve up to 20 MPa for the mixtures containing R-1234yf and 49 MPa for the R-1234ze(E)/227ea mixtures. The relative combined expanded speed of sound uncertainty ranged from 0.039 % to 0.317 % with a mean uncertainty over all state points of less than 0.10 %. The reported data are compared to the most recent mixture models in REFPROP for each blend studied. Comparisons of available mixture models for the R-125/1234yf and R-1234yf/152a blends exhibit average absolute deviation values ranging from 0.10 % to 0.27 %, and the average absolute deviations for R-1234ze(E)/227ea blends range from 0.62 % to 0.94 %. The comparisons show that only minor adjustments are required to the R-125/1234yf and R-1234yf/152a mixture models to represent the speed of sound data within its uncertainty. However, significant adjustments are needed to improve the current mixture models for the R-1234ze(E)/227ea blend. Deficiencies with the R-1234ze(E)/227ea mixture model are not unexpected since it (1) utilizes a pure-fluid EOS R-1234ze(E) that has been shown to inaccurately represent R-1234ze(E) speed of sound values and (2) uses binary interaction parameters for the chemically similar mixture of R-1234yf/227ea.

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

Velocity Profiling of a Gas–Solid Fluidized Bed Using Electrical Capacitance Volume Tomography

In this study, a method of producing velocity profile maps from electrical capacitance volume tomography (ECVT) measurements by reconstructing displacement from measured changes in capacitance is developed and applied to fluidized bed systems. The mapping of the reconstruction leverages the gradient of the sensitivity distribution of the ECVT sensor to circumvent the need for image cross correlation techniques. Experimental data of both bubbling and slugging fluidized beds are collected in a cold flow model. Adaptation of the technique is discussed in detail, and velocity profiles are obtained for a range of gas flow rates. The produced velocity maps are compared against the established methods of cross correlation and against empirical correlations from the literature and are found to agree well in tracking slug and bubble velocity. The exception is when the tracked object is large relative to the ECVT sensor dimensions, a scenario that can be avoided through proper sensor design. The quantities of average velocity, momentum, and solid and gas volume fraction are derived from the image and velocity profiles. The results demonstrate and extend the power of ECVT as a measurement tool for the study and monitoring of gas–solid fluidized beds by providing a computationally cheaper alternative to 3-D cross correlation for deriving velocity profiles.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Development of Integrated Thermal Fluids Modeling Capability for MSRs

The DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad class of light-water and advanced reactor concepts. The research and development approach under thermal fluids technical area synergistically combines three length and time scales in a hierarchal multi-scale approach. To demonstrate the feasibility and capabilities of a multi-scale thermal fluids capability using these codes, a key joint effort has been pursued to develop an integrated system-and engineering-scale thermal fluids analysis capability with the MOOSE-based codes, through integration of SAM and Pronghorn, both based on the MOOSE framework. This report summarizes the progress in developing an integrated system- and engineering-scale thermal fluids analysis capability based on SAM and Pronghorn for molten salt reactors (MSRs), which gained significant interest in recently years. Two coupling approaches were studied, i.e. separate domain or domain-segregated coupling approach and the domain-overlapping approach. A series of coupled multi-physics models have been developed for a common reference molten salt fast reactor (MSFR) concept, ranging from standalone SAM system model to integrated SAM-Pronghorn-Griffin models. Both the steady state and the transient simulations are performed to the state-of-the-art simulation capabilities of NEAMS software suite in MSFR system applications. This report also covers further development and testing of the gas transport model in SAM for MSR modeling support. The presence of noncondensable gas in MSR systems would have strong impacts on fission gas removal and transport of noble metals throughout the system. Fission products removed through the off-gas system can also impact reactivity and can act as an additional point of heat removal. To ensure that the gas transport model is adequately tested, the supported modeling capabilities and features of SAM were identified, and a suite of tests were developed to test the model for each identified feature. Validation and UQ testing were also performed for the model and demonstrated that the buoyancy term, which was originally developed using non-salt/helium gas experimental data, can capture the experimental gas bubble velocity and diameter to within experimental and code uncertainty. An existing MSRE model was also modified in this work to include the gas transport model and to demonstrate the gas model behavior for realistic conditions of interest.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗