Fluids

Research

18 pages, 22383 KiB

Open AccessArticle

Numerical Study on an Interface Compression Method for the Volume of Fluid Approach

by Yuria Okagaki, Taisuke Yonomoto, Masahiro Ishigaki and Yoshiyasu Hirose

Fluids 2021, 6(2), 80; https://doi.org/10.3390/fluids6020080 - 10 Feb 2021

Cited by 17 | Viewed by 3574

Many thermohydraulic issues about the safety of light water reactors are related to complicated two-phase flow phenomena. In these phenomena, computational fluid dynamics (CFD) analysis using the volume of fluid (VOF) method causes numerical diffusion generated by the first-order upwind scheme used in [...] Read more.

Many thermohydraulic issues about the safety of light water reactors are related to complicated two-phase flow phenomena. In these phenomena, computational fluid dynamics (CFD) analysis using the volume of fluid (VOF) method causes numerical diffusion generated by the first-order upwind scheme used in the convection term of the volume fraction equation. Thus, in this study, we focused on an interface compression (IC) method for such a VOF approach; this technique prevents numerical diffusion issues and maintains boundedness and conservation with negative diffusion. First, on a sufficiently high mesh resolution and without the IC method, the validation process was considered by comparing the amplitude growth of the interfacial wave between a two-dimensional gas sheet and a quiescent liquid using the linear theory. The disturbance growth rates were consistent with the linear theory, and the validation process was considered appropriate. Then, this validation process confirmed the effects of the IC method on numerical diffusion, and we derived the optimum value of the IC coefficient, which is the parameter that controls the numerical diffusion. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Figure 1

17 pages, 2273 KiB

Open AccessArticle

Monolithic Solvers for Incompressible Two-Phase Flows at Large Density and Viscosity Ratios

by Mohamed El Ouafa, Stephane Vincent and Vincent Le Chenadec

Fluids 2021, 6(1), 23; https://doi.org/10.3390/fluids6010023 - 5 Jan 2021

Cited by 8 | Viewed by 3093

Abstract

In this paper, we investigate the accuracy and robustness of three classes of methods for solving two-phase incompressible flows on a staggered grid. Here, the unsteady two-phase flow equations are simulated by finite volumes and penalty methods using implicit and monolithic approaches (such [...] Read more.

In this paper, we investigate the accuracy and robustness of three classes of methods for solving two-phase incompressible flows on a staggered grid. Here, the unsteady two-phase flow equations are simulated by finite volumes and penalty methods using implicit and monolithic approaches (such as the augmented Lagrangian and the fully coupled methods), where all velocity components and pressure variables are solved simultaneously (as opposed to segregated methods). The interface tracking is performed with a Volume-of-Fluid (VOF) method, using the Piecewise Linear Interface Construction (PLIC) technique. The home code Fugu is used for implementing the various methods. Our target application is the simulation of two-phase flows at high density and viscosity ratios, which are known to be challenging to simulate. The resulting strategies of monolithic approaches will be proven to be considerably better suited for these two-phase cases, they also allow to use larger time step than segregated methods. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Figure 1

18 pages, 7518 KiB

Open AccessArticle

Assessment and Prediction of Air Entrainment and Geyser Formation in a Bottom Outlet: Field Observations and CFD Simulation

by James Yang, Penghua Teng, Junhu Nan, Shicheng Li and Anders Ansell

Fluids 2020, 5(4), 203; https://doi.org/10.3390/fluids5040203 - 9 Nov 2020

Cited by 1 | Viewed by 2761

Abstract

Air entrainment at the intake of a bottom outlet often gives rise to air pockets in its conduit and formation of geysers. The outlet in question comprises a bulkhead gate, gate shaft, horizontal conduit, and exit. Operations show that it suffers from appreciable [...] Read more.

Air entrainment at the intake of a bottom outlet often gives rise to air pockets in its conduit and formation of geysers. The outlet in question comprises a bulkhead gate, gate shaft, horizontal conduit, and exit. Operations show that it suffers from appreciable flow fluctuations and blowouts in the tailwater, which leads to gate operation restrictions. For the purpose of understanding the hydraulic phenomenon, both prototype discharge tests and three-dimensional computational fluid dynamics (CFD) modeling of two-phase flows are performed. The operational focus of the facility are small and large gate openings. The CFD results reveal that, with air entrained in the gate shaft, continual breakup and coalescence of air bubbles in the conduit typify the flow. At small openings below 1 meter, the air–water flow is characterized by either distinct blowouts of regular frequency or continuous air release. In terms of geyser behaviors inclusive of frequency, the agreement is good between field and numerical studies. At large openings, the gate becomes fully submerged, and the flow is discharged without air entrainment and blowouts. The paper showcases the air–water flow features in a typical bottom outlet layout in Sweden, which is intended to serve as an illustration of the study procedure for other similar outlets. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Graphical abstract

18 pages, 833 KiB

Open AccessArticle

Control and Optimization of Interfacial Flows Using Adjoint-Based Techniques

by Alexandru Fikl, Vincent Le Chenadec and Taraneh Sayadi

Fluids 2020, 5(3), 156; https://doi.org/10.3390/fluids5030156 - 10 Sep 2020

Cited by 7 | Viewed by 2594

Abstract

The applicability of adjoint-based gradient computation is investigated in the context of interfacial flows. Emphasis is set on the approximation of the transport of a characteristic function in a potential flow by means of an algebraic volume-of-fluid method. A class of optimisation problems [...] Read more.

The applicability of adjoint-based gradient computation is investigated in the context of interfacial flows. Emphasis is set on the approximation of the transport of a characteristic function in a potential flow by means of an algebraic volume-of-fluid method. A class of optimisation problems with tracking-type functionals is proposed. Continuous (differentiate-then-discretize) and discrete (discretize-then-differentiate) adjoint-based gradient computations are formulated and compared in a one-dimensional configuration, the latter being ultimately used to perform optimisation in two dimensions. The gradient is used in truncated Newton and steepest descent optimisers, and the algorithms are shown to recover optimal solutions. These validations raise a number of open questions, which are finally discussed with directions for future work. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Figure 1

24 pages, 6023 KiB

Open AccessArticle

Numerical Study on the Flow and Heat Transfer Coupled in a Rectangular Mini-Channel by Finite Element Method for Industrial Micro-Cooling Technologies

by Claude Aurélien Kamdem Kamdem and Xiaolu Zhu

Fluids 2020, 5(3), 151; https://doi.org/10.3390/fluids5030151 - 3 Sep 2020

Cited by 3 | Viewed by 3723

Abstract

Nowadays, cooling high thermal flows in compact volumes continues to be one of the crucial problems in the industry. With the advent of advanced technologies, much more attention has been paid to how to improve the performance of cooling systems in the area [...] Read more.

Nowadays, cooling high thermal flows in compact volumes continues to be one of the crucial problems in the industry. With the advent of advanced technologies, much more attention has been paid to how to improve the performance of cooling systems in the area of micro-technologies. Rectangular mini-channels are typical representatives which commonly used for cooling applications. However, micro-technologies still face the problem of low performance due to the low productivity of cooling related to unbefitting physical parameter values. Here, this work studies the applicability of the heat transfer scheme of convective flow and flow boiling in a rectangular mini-channel for satisfying the cooling requirement of industrial micro-technologies, through a simulation model governed by the coupled mechanism from Navier-Stokes (N-S) equation and heat transfer equations with phase change effect. In this work, various hydraulic diameters and different inlet fluid speed are used to calculate the different velocity profiles, pressure drops, coefficients of friction and finally, the distribution of the temperatures and dissipated heat flux. The simulation results show the applicability of the rectangular mini-channel in diverse applications such as engine cooling, electronic components, automotive on-board electronics and aerospace engineering. Flow boiling simulation results reveal that the obtained patterns were smooth mixture flow and discrete flow. The dissipated heat flux can reach 1.02–5.34 MW/m² for a hydraulic diameter of 0.5 mm. We show that the system with the gradient temperature that evolves increasingly along the top and bottom walls of the channels presents the highest heat flux dissipated in flow boiling. Additionally, the fin efficiency of the system is 0.88 and the coefficient value of convective heat transfer is in the range between 5000 < h < 100,000, which indicates the flow boiling heat transfer is effective in the mini-channel when the Reynolds number is less than 400. It provides a significant heat exchange for cooling in these application areas. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Graphical abstract

20 pages, 2342 KiB

Open AccessArticle

Computational Modeling of Bubbles Growth Using the Coupled Level Set—Volume of Fluid Method

by Amir Taqieddin, Yuxuan Liu, Akram N. Alshawabkeh and Michael R. Allshouse

Fluids 2020, 5(3), 120; https://doi.org/10.3390/fluids5030120 - 23 Jul 2020

Cited by 7 | Viewed by 4880

Abstract

Understanding the generation, growth, and dynamics of bubbles as they absorb or release dissolved gas in reactive flows is crucial for optimizing the efficiency of electrochemically gas-evolving systems like alkaline water electrolysis or hydrogen production. To better model these bubbly flow systems, we [...] Read more.

Understanding the generation, growth, and dynamics of bubbles as they absorb or release dissolved gas in reactive flows is crucial for optimizing the efficiency of electrochemically gas-evolving systems like alkaline water electrolysis or hydrogen production. To better model these bubbly flow systems, we use a coupled level set and volume of fluid approach integrated with a one-fluid transport of species model to study the dynamics of stationary and rising bubbles in reactive two-phase flows. To accomplish this, source terms are incorporated into the continuity and phase conservation equations to allow the bubble to grow or shrink as the species moves through the interface. Verification of the hydrodynamics of the solver for non-reactive systems demonstrates the requisite high fidelity interface capturing and mass conservation necessary to incorporate transport of species. In reactive systems where the species impacts the bubble volume, the model reproduces the theoretically predicted and experimentally measured diffusion-controlled growth rate (i.e.,

R (t) \propto t^{0.5}

). The model is then applied to rising bubbles to demonstrate the impact of transport of species on both the bubble velocity and shape as well as the concentration field in its wake. This improved model enables the incorporation of electric fields and chemical reactions that are essential for studying the physicochemical hydrodynamics in multiphysics systems. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Graphical abstract

19 pages, 4594 KiB

Open AccessEditor’s ChoiceArticle

CFD Modeling of Hydrocyclones—A Study of Efficiency of Hydrodynamic Reservoirs

by Marvin Durango-Cogollo, Jose Garcia-Bravo, Brittany Newell and Andres Gonzalez-Mancera

Fluids 2020, 5(3), 118; https://doi.org/10.3390/fluids5030118 - 21 Jul 2020

Cited by 22 | Viewed by 8752

Abstract

The dynamics of hydrocyclones is complex, because it is a multiphase flow problem that involves interaction between a discrete phase and multiple continuum phases. The performance of hydrocyclones is evaluated by using Computational Fluid Dynamics (CFD), and it is characterized by the pressure [...] Read more.

The dynamics of hydrocyclones is complex, because it is a multiphase flow problem that involves interaction between a discrete phase and multiple continuum phases. The performance of hydrocyclones is evaluated by using Computational Fluid Dynamics (CFD), and it is characterized by the pressure drop, split water ratio, and particle collection efficiency. In this paper, a computational model to improve and evaluate hydrocyclone performance is proposed. Four known computational turbulence models (renormalization group (RNG) k-

ε

, Reynolds stress model (RSM), and large-eddy simulation (LES)) are implemented, and the accuracy of each for predicting the hydrocyclone behavior is assessed. Four hydrocyclone configurations were analyzed using the RSM model. By analyzing the streamlines resulting from those simulations, it was found that the formation of some vortices and saddle points affect the separation efficiency. Furthermore, the effects of inlet width, cone length, and vortex finder diameter were found to be significant. The cut-size diameter was decreased by 33% compared to the Hsieh experimental hydrocyclone. An increase in the pressure drop leads to high values of cut-size and classification sharpness. If the pressure drop increases to twice its original value, the cut-size and the sharpness of classification are reduced to less than 63% and 55% of their initial values, respectively. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Graphical abstract

17 pages, 5534 KiB

Open AccessArticle

Assessment of CFD Solvers and Turbulent Models for Water Free Jets in Spillways

by António Muralha, José F. Melo and Helena M. Ramos

Fluids 2020, 5(3), 104; https://doi.org/10.3390/fluids5030104 - 30 Jun 2020

Cited by 9 | Viewed by 3684

Abstract

The capability of two different OpenFOAM^® solvers, namely interFoam and twoPhaseEulerFoam, in reproducing the behavior of a free water jet was investigated. Numerical simulations were performed in order to obtain the velocity and air concentration profiles along the jet. The turbulence intensity [...] Read more.

The capability of two different OpenFOAM^® solvers, namely interFoam and twoPhaseEulerFoam, in reproducing the behavior of a free water jet was investigated. Numerical simulations were performed in order to obtain the velocity and air concentration profiles along the jet. The turbulence intensity was also analyzed. The obtained results were compared with published experimental data and, in general, similar velocity and air concentration profiles were found. InterFoam solver is able to reproduce the velocity field of the free jet but has limitations in the simulation of the air concentration. TwoPhaseEulerFoam performs better in reproducing the air concentration along the jet, the results being in agreement with the experimental data, although the computational runs are less stable and more time consuming. The sensitivity analysis of the inlet turbulent intensity showed that it has no influence in the characteristics of the jet core. With this research it is possible to conclude that: interFoam with k-Epsilon (k-ε) turbulence model is the best choice if the goal of the numerical simulations is the simulation of the velocity field of the jet. Meanwhile, twoPhaseEulerFoam with mixturek-Epsilon (mk-ε) shall be considered if the objective is the simulation of the velocity field and the air concentration. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Figure 1

17 pages, 8458 KiB

Open AccessArticle

Numerical Simulation of Breathing Mode Oscillation on Bubble Detachment

by Takao Oku, Hiroyuki Hirahara, Tomohiro Akimoto and Daiki Tsuchida

Fluids 2020, 5(2), 96; https://doi.org/10.3390/fluids5020096 - 16 Jun 2020

Cited by 4 | Viewed by 3523

Abstract

When a bubble detaches from a nozzle immersed in water, a sound is emitted owing to the detachment. The bubble deformation and sound emission generated after detachment has been investigated in many studies, in which the breathing mode with a natural frequency was [...] Read more.

When a bubble detaches from a nozzle immersed in water, a sound is emitted owing to the detachment. The bubble deformation and sound emission generated after detachment has been investigated in many studies, in which the breathing mode with a natural frequency was discussed based on the dynamics of the interface between the air and water. In this study, the deformation of a bubble was observed, and the sound emitted upon detachment was measured experimentally. To analyze the bubble deformation process, a computational fluid dynamics (CFD) simulation was conducted using the volume of fluid (VOF) method to predict the sound emission. In the analysis, the deformation behavior, the oscillation frequencies, sound pressure, and radius variation were discussed by comparing the numerical and experimental data. Furthermore, the natural frequency and low frequency vibrations were discussed based on the interference between the detached bubbles and the air column vibrations. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Figure 1

23 pages, 2099 KiB

Open AccessFeature PaperArticle

Consistent Velocity–Pressure Coupling for Second-Order L²-Penalty and Direct-Forcing Methods

by Arthur Sarthou, Stéphane Vincent and Jean-Paul Caltagirone

Fluids 2020, 5(2), 92; https://doi.org/10.3390/fluids5020092 - 10 Jun 2020

Cited by 1 | Viewed by 2161

Abstract

The present work studies the interactions between fictitious-domain methods on structured grids and velocity–pressure coupling for the resolution of the Navier–Stokes equations. The pressure-correction approaches are widely used in this context but the corrector step is generally not modified consistently to take into [...] Read more.

The present work studies the interactions between fictitious-domain methods on structured grids and velocity–pressure coupling for the resolution of the Navier–Stokes equations. The pressure-correction approaches are widely used in this context but the corrector step is generally not modified consistently to take into account the fictitious domain. A consistent modification of the pressure-projection for a high-order penalty (or penalization) method close to the Ikeno–Kajishima modification for the Immersed Boundary Method is presented here. Compared to the first-order correction required for the

L^{2}

-penalty methods, the small values of the penalty parameters do not lead to numerical instabilities in solving the Poisson equation. A comparison of the corrected rotational pressure-correction method with the augmented Lagrangian approach which does not require a correction is carried out. Full article

(This article belongs to the Special Issue Advances in Numerical Methods for Multiphase Flows)

► Show Figures

Figure 1

Journal Menu

Journal Browser

Advances in Numerical Methods for Multiphase Flows

Share This Special Issue

Special Issue Editor

Special Issue Information

Keywords

Benefits of Publishing in a Special Issue

Published Papers (10 papers)

Research

Further Information

Guidelines

MDPI Initiatives

Follow MDPI