Search NASA⌕ Search

SEARCH · Search NASA

Results for “bubbly flow”

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

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Processing Methodology for Double Sensor Probe in a Bubbly Flow

Bubbly flow phenomenon is the most common flow pattern in nuclear plant system. In this paper, according to Hibiki and Bartel, the experiment was designed to observe the fluid conditions surrounding nuclear rod. In recent years, double sensor probe is one of the measuring techniques for the bubbly flow. It can measure local void fraction and bubble velocity. After mathematical statistical calculation, interfacial area concentration is obtained by Kataoka. The double sensor probe method, with the advantage for measuring local interfacial area concentration, has become widely used in many measurement methods. The paper of Liu and Zhao were referred to establish the post-processing for the raw data obtained from measuring instrument. After drawing the signal of the probe, it is important to analyze the raw data for the elimination of unnecessary noise signal and highlight the characteristic signal. Based on the same experimental condition, the different filtering methods, Simple Moving Average (SMA) and Gaussian average, were compared and discussed by the double sensor probe.

42 ENGINEERING↗
DOE OSTI2017

Interface tracking simulations of bubbly flows in PWR relevant geometries

[en] Highlights: • Simulations were performed for turbulent bubbly flows in PWR subchannel geometry. • Liquid turbulence is fully resolved by direct numerical simulation approach. • Bubble behavior is captured using level-set interface tracking method. • Time-averaged single- and two-phase turbulent flow statistical quantities are obtained. - Abstract: The advances in high performance computing (HPC) have allowed direct numerical simulation (DNS) approach coupled with interface tracking methods (ITM) to perform high fidelity simulations of turbulent bubbly flows in various complex geometries. In this work, we have chosen the geometry of the pressurized water reactor (PWR) core subchannel to perform a set of interface tracking simulations (ITS) with fully resolved liquid turbulence. The presented research utilizes a massively parallel finite-element based code, PHASTA, for the subchannel geometry simulations of bubbly flow turbulence. The main objective for this research is to demonstrate the ITS capabilities in gaining new insight into bubble/turbulence interactions and assisting the development of improved closure laws for multiphase computational fluid dynamics (M-CFD). Both single- and two-phase turbulent flows were studied within a single PWR subchannel. The analysis of numerical results includes the mean gas and liquid velocity profiles, void fraction distribution and turbulent kinetic energy profiles. Two sets of flow rates and bubble sizes were used in the simulations. The chosen flow rates corresponded to the Reynolds numbers of 29,079 and 80,775 based on channel hydraulic diameter (D_h) and mean velocity. The finite element unstructured grids utilized for these simulations include 53.8 million and 1.11 billion elements, respectively. This has allowed to fully resolve all the turbulence scales and the deformable interfaces of individual bubbles. For the two-phase flow simulations, a 1% bubble volume fraction was used which resulted in 17 bubbles in the smaller case and 262 bubbles in the larger case. In the larger simulation case the size of the resolved bubbles is 0.65 mm in diameter, and the bulk mesh cell size is about 30 microns. Those large-scale simulations provide new level of details previously unavailable and were enabled by the excellent scaling performance of our two-phase flow solver and access to the state-of-the-art supercomputing resources. The presented simulations used up to 256 thousand processing threads on the IBM BG/Q supercomputer “Mira” (Argonne National Laboratory).

Nuclear Science & Technology↗
DOE OSTI2020

Uncertainty analysis of PIV measurements in bubbly flows considering sampling and bubble effects with ray optics modeling

In this report, a modeling-based approach is developed to investigate the external sources of uncertainty in particle image velocimetry (PIV) measurements for two-phase bubbly flows. Primarily relied on ray optics modeling, this study focuses on the uncertainty in PIV measurements due to bubble-induced light distortion. In addition, the uncertainty due to sampling is also considered through synthetic PIV images generated from Monte Carlo sampling. Based on a PIV-planar laser induced florescence (PLIF) measurement system, three case studies are performed and discussed, considering respectively the sampling effect, bubble effect, and the combination of these two effects. These case studies confirm the validity of the ray optics modeling in obtaining accurate particle images through a benchmark test. Furthermore, it is also found that the bubble effect could significantly influence the particle images. On the other hand, the result of the case study with synthetic PIV images indicates that the bubble effect on particle images can be minimized by removing “outliers” in the obtained velocity vector field through a multi-pass post-processing procedure. However, the uncertainty of the velocity measurements will still increase.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Eulerian-Eulerian large eddy simulation of two-phase dilute bubbly flows

An efficient mathematical model is presented for conducting large eddy simulation (LES) of turbulent bubbly flows at low gas hold-ups. The mathematical model is based on the filtered multi-fluid equations derived by application of the component-weighted volume-averaging process of each phase. A two-fluid model is applied to the filtered set of equations, yielding a more computationally tractable set of filtered equations. The filtered equations are formulated in generalized curvilinear coordinates to simulate turbulent gas-liquid bubbly flows in complex geometries. The Germano model is used to account for sub-grid turbulent stresses. The gas and liquid momentum equations are discretized with the hybrid staggered/non-staggered method. The performance of the developed mathematical model is demonstrated by simulating the benchmark test case of a turbulent gas-liquid bubbly flow in a rectangular bubble column. The accuracy of the proposed method is assessed via comparison with the experimental data and other numerical results in the literature.

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

Evaluation of interfacial area transport models for horizontal bubbly flow

This paper presents evaluation results of interfacial area transport equation (IATE) for horizontal air-water bubbly flow. The steady-state one-dimensional, one-group IATE for adiabatic air-water bubbly flow developed by Talley [1] is employed in the current evaluation study. The models evaluated include: (a) drift-flux-based closure relations for void-weighted bubble velocity and void fraction; (b) frictional pressure drop model based on the Lockhart-Martinelli approach; (c) bubble interaction mechanisms; and (d) covariance associated with asymmetric bubble distribution in horizontal two-phase flow. In order to evaluate the models, a local two-phase flow database obtained by the four-sensor conductivity probe in fourteen different test conditions is employed, which covers a wide range of flow conditions in horizontal bubbly flow regime, with superficial liquid and gas velocities ranging from 3.50-6.00 m/s and 0.09- 1.00 m/s, respectively. It is demonstrated that both the void weighted bubble velocity and void fraction are predicted well with an average percent difference of ±2.7%. The fictional pressure loss can be predicted with an average percent difference of ±1.12%. The covariance parameters calculated from the modeling approach generally compare well with those calculated from experimental data within ±20% difference. It is found that the IATE generally predicts the measured interfacial area concentration with an average percent difference of ±5.7%. Furthermore, the effect of liquid velocity on individual bubble interaction mechanisms is investigated. It is found that both the contributions of random collision and turbulent impact decrease as liquid velocity decreases. (authors)

42 ENGINEERING↗
DOE OSTI2018

Interface Tracking Investigation of Geometric Effects on the Bubbly Flow in PWR Subchannels

We report that absorbing heat from the fuel rod surface, water as coolant can undergo subcooled boiling within a pressurized water reactor (PWR) fuel rod bundle. Because of the buoyancy effect, the vapor bubbles generated will then rise along and interact with the subchannel geometries. Reliable prediction of bubble behavior is of immense importance to ensure safe and stable reactor operation. However, given a complex engineering system like a nuclear reactor, it is very challenging (if not impossible) to conduct high-resolution measurements to study bubbly flows under reactor operation conditions. The lack of a fundamental two-phase-flow database is hindering the development of accurate two-phase-flow models required in more advanced reactor designs. In response to this challenge, first-principles–based numerical simulations are emerging as an attractive alternative to produce a complementary data source along with experiments. Leveraged by the unprecedented computing power offered by state-of-the-art supercomputers, direct numerical simulation (DNS), coupled with interface tracking methods, is becoming a practical tool to investigate some of the most challenging engineering flow problems. In the presented research, turbulent bubbly flow is simulated via DNS in single PWR subchannel geometries with auxiliary structures (e.g., supporting spacer grid and mixing vanes). The geometric effects these structures exert on the bubbly flow are studied with both a conventional time-averaging approach and a novel dynamic bubble tracking method. Lastly, the new insights obtained will help inform better two-phase models that can contribute to safer and more efficient nuclear reactor systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Two-Phase Turbulent Kinetic Energy Budget Computation in Co-Current Taylor Bubble Flow

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

DEMONSTRATION OF A DATA-DRIVEN PHYSICS-BASED APPROACH FOR COMPUTATIONALLY EFFICIENT CFD PREDICTION OF TWO-PHASE BUBBLY FLOW

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale physics-guided data-driven approach, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. By defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

42 ENGINEERING↗
DOE OSTI2020

Computationally efficient CFD prediction of bubbly flow using physics-guided deep learning

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale framework was proposed in this paper by applying a physics-guided data-driven approach. Instrumental to this framework, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. In this work, by defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Surfactant stabilized bubbles flowing in a Newtonian fluid

Bubbles suspended in a fluid cause the suspension to have different rheological properties than the base fluid. In general, the viscosity of the suspension increases as the volume fraction of the bubbles is increased. A current application, and motivation for this study, is in wellbore cements used for hydrocarbon extraction and carbon sequestration. In these settings, the gas bubbles are dispersed into the cement to reduce the density as well as improve the properties for specific conditions or wellbore issues. In this paper, we use Stokesian dynamics to numerically simulate the behavior of a large number of bubbles suspended in a Newtonian fluid. Going beyond prior work on simulating particles in suspension, we account for the nature of bubbles by allowing for slip on the bubble surface, the deflection on the bubble surface, and a bubble–bubble pairwise interaction that represents the surfactant physics; we do not account for bubble compressibility. We incorporate these interactions and simulate bubble suspensions of monodisperse size at several volume fractions. We find that the bubbles remain better dispersed compared with hard spherical particles that show a greater tendency to structure or cluster.

Rosenbaum, Eilis↗
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2022

Quantitatively Monitoring Bubble-Flow at a Seep Site Offshore Oregon: Field Trials and Methodological Advances for Parallel Optical and Hydroacoustical Measurements

Two lander-based devices, the Bubble-Box and GasQuant-II, were used to investigate the spatial and temporal variability and total gas flow rates of a seep area offshore Oregon, United States. The Bubble-Box is a stereo camera–equipped lander that records bubbles inside a rising corridor with 80 Hz, allowing for automated image analyses of bubble size distributions and rising speeds. GasQuant is a hydroacoustic lander using a horizontally oriented multibeam swath to record the backscatter intensity of bubble streams passing the swath plain. The experimental set up at the Astoria Canyon site at a water depth of about 500 m aimed at calibrating the hydroacoustic GasQuant data with the visual Bubble-Box data for a spatial and temporal flow rate quantification of the site. For about 90 h in total, both systems were deployed simultaneously and pressure and temperature data were recorded using a CTD as well. Detailed image analyses show a Gaussian-like bubble size distribution of bubbles with a radius of 0.6–6 mm (mean 2.5 mm, std. dev. 0.25 mm); this is very similar to other measurements reported in the literature. Rising speeds ranged from 15 to 37 cm/s between 1- and 5-mm bubble sizes and are thus, in parts, slightly faster than reported elsewhere. Bubble sizes and calculated flow rates are rather constant over time at the two monitored bubble streams. Flow rates of these individual bubble streams are in the range of 544–1,278 mm 3 /s. One Bubble-Box data set was used to calibrate the acoustic backscatter response of the GasQuant data, enabling us to calculate a flow rate of the ensonified seep area (~1,700 m 2 ) that ranged from 4.98 to 8.33 L/min (5.38 × 10 6 to 9.01 × 10 6 CH 4 mol/year). Such flow rates are common for seep areas of similar size, and as such, this location is classified as a normally active seep area. For deriving these acoustically based flow rates, the detailed data pre-processing considered echogram gridding methods of the swath data and bubble responses at the respective water depth. The described method uses the inverse gas flow quantification approach and gives an in-depth example of the benefits of using acoustic and optical methods in tandem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation

Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗