Initial Impact of Liquid Jets and Entrainment of Air Bubbles: A Numerical Study

Author(s):  
Xiaoliang Qu ◽  
Lyes Khezzar ◽  
Zhenlin Li

This paper presents a three-dimensional unsteady numerical simulation of a turbulent plunging liquid jet without artificial surface disturbance impinging on a quiecent liquid pool. The focal point of the study is the initial impact and air entrainment process. The multiphase, Volume of Fluid Model is used in combination with the Reynolds Averaged k-ε turbulence model. The process of the initial impact of the jet on the free surface, the subsequent formation of an air cavity and the subsequent break-down of the cavity into small bubbles are captued and analyzed. These simulations show clearly and in detail the process of air carryunder by the liquid-liquid jet. The air cavity caused by the intial jet impact stretches deeply under the pool surface untill break down due to the shear created by a torroidal vortex. The predicted maximum height of the developing air cavity shows very good agreement with existing semi-empirical correlations from the literature and experiments. The velocity of the front of the air cavity is equal to about half the jet valocity at impact as shown by previous works and the predicted penetration depth shows acceptable agreement with previous correlations. The VOF model shows a strong capability of tracking the interface between two phases.

2018 ◽  
Vol 30 (1) ◽  
Author(s):  
Bentalha Chakib

Stepped spillway is a power full hydraulic structure for energy dissipation because ofthe large value of the surface roughness. The performance of the stepped spillway is enhancedwith the presence of air that can prevent or reduce the cavitation damage. This work aims tosimulate air entrainment and determine the characteristics of flow at stepped spillways. Withinthis work flow over stepped chute is simulated by using fluent computational fluid dynamics(CFD). The volume of fluid (VOF) model is used as a tool to simulate air-water interaction onthe free surface thereby the turbulence closure is derived in the k −ε turbulence standard model.The found numerical results agree well with experimental results.


2019 ◽  
Vol 42 (1) ◽  
pp. 42-48
Author(s):  
Chakib Bentalha ◽  
Mohammed Habi

Abstract Stepped spillway is hydraulic structure designed to dissipate the excess in kinetic energy at the downstream of dams and can reduce the size of stilling basin at the toe of the spillway or chute. The flow on a stepped spillway is characterised by the large aeration that can prevent or reduce the cavitation damage. The air entrainment starts where the boundary layer attains the free surface of flow; this point is called “point of inception”. Within this work the inception point is determined by using software Ansys Fluent where the volume of fluid (VOF) model is used as a tool to track the free surface thereby the turbulence closure is derived in the k − ε turbulence standard model. This research aims to find new formulas for describe the variation of water depth at step edge and the positions of the inception point, at the same time the contour map of velocity, turbulent kinetic energy and strain rate are presented. The found numerical results agree well with experimental results like the values of computed and measured water depth at the inception point and the numerical and experimental inception point locations. Also, the dimensionless water depth profile obtained by numerical method agrees well with that of measurement. This study confirmed that the Ansys Fluent is a robust software for simulating air entrainment and exploring more characteristics of flow over stepped spillways.


2020 ◽  
Vol 53 (34) ◽  
pp. 345202
Author(s):  
Feng Gao ◽  
Qiao Wang ◽  
Dong Dai ◽  
Wenjun Ning ◽  
Yuhui Zhang ◽  
...  

2021 ◽  
pp. 174425912098418
Author(s):  
Toivo Säwén ◽  
Martina Stockhaus ◽  
Carl-Eric Hagentoft ◽  
Nora Schjøth Bunkholt ◽  
Paula Wahlgren

Timber roof constructions are commonly ventilated through an air cavity beneath the roof sheathing in order to remove heat and moisture from the construction. The driving forces for this ventilation are wind pressure and thermal buoyancy. The wind driven ventilation has been studied extensively, while models for predicting buoyant flow are less developed. In the present study, a novel analytical model is presented to predict the air flow caused by thermal buoyancy in a ventilated roof construction. The model provides means to calculate the cavity Rayleigh number for the roof construction, which is then correlated with the air flow rate. The model predictions are compared to the results of an experimental and a numerical study examining the effect of different cavity designs and inclinations on the air flow rate in a ventilated roof subjected to varying heat loads. Over 80 different test set-ups, the analytical model was found to replicate both experimental and numerical results within an acceptable margin. The effect of an increased total roof height, air cavity height and solar heat load for a given construction is an increased air flow rate through the air cavity. On average, the analytical model predicts a 3% higher air flow rate than found in the numerical study, and a 20% lower air flow rate than found in the experimental study, for comparable test set-ups. The model provided can be used to predict the air flow rate in cavities of varying design, and to quantify the impact of suggested roof design changes. The result can be used as a basis for estimating the moisture safety of a roof construction.


2021 ◽  
Vol 87 (2) ◽  
Author(s):  
Vikram S. Dharodi ◽  
Amita Das

Rayleigh–Taylor (RT) and buoyancy-driven (BD) instabilities are driven by gravity in a fluid system with inhomogeneous density. The paper investigates these instabilities for a strongly coupled dusty plasma medium. This medium has been represented here in the framework of the generalized hydrodynamics (GHD) fluid model which treats it as a viscoelastic medium. The incompressible limit of the GHD model is considered here. The RT instability is explored both for gradual and sharp density gradients stratified against gravity. The BD instability is discussed by studying the evolution of a rising bubble (a localized low-density region) and a falling droplet (a localized high-density region) in the presence of gravity. Since both the rising bubble and falling droplet have symmetry in spatial distribution, we observe that a falling droplet process is equivalent to a rising bubble. We also find that both the gravity-driven instabilities get suppressed with increasing coupling strength of the medium. These observations have been illustrated analytically as well as by carrying out two-dimensional nonlinear simulations. Part 2 of this paper is planned to extend the present study of the individual evolution of a bubble and a droplet to their combined evolution in order to understand the interaction between them.


Author(s):  
Jean Franc¸ois Sigrist ◽  
Christian Laine ◽  
Dominique Lemoine ◽  
Bernard Peseux

This paper is related to the study of a nuclear propulsion reactor prototype for the French Navy. This prototype is built on ground and is to be dimensioned toward seismic loading. The dynamic analysis takes the coupled fluid structure analysis into account. The basic fluid models used by design engineers are inviscid incompressible or compressible. The fluid can be described in a bidimensional by slice or a three-dimensional approach. A numerical study is carried out on a generic problem for the linear FSI dynamic problem. The results of this study are presented and discussed. As a conclusion, the three-dimensional inviscid incompressible fluid appears to be the best compromise between the description of physical phenomena and the cost of modeling. The geometry of the reactor is such that large displacements of the structure in the fluid can occur. Therefore, the linearity hypothesis might not be longer valid. The case of large amplitude imposed oscillating motion of a cylinder in a confined fluid is numerically studied. A CFD code is used to investigate the fluid behavior solving the NAVIER-STOKES equations. The forces induced on the cylinder by the fluid are computed and compared to the linear solution. The limit of the linear model can then be exhibited.


Processes ◽  
2020 ◽  
Vol 8 (7) ◽  
pp. 810
Author(s):  
Jade Gesare Abuga ◽  
Tiri Chinyoka

The flow of viscoelastic fluids may, under certain conditions, exhibit shear-banding characteristics that result from their susceptibility to unusual flow instabilities. In this work, we explore both the existing shear banding mechanisms in the literature, namely; constitutive instabilities and flow-induced inhomogeneities. Shear banding due to constitutive instabilities is modelled via either the Johnson–Segalman or the Giesekus constitutive models. Shear banding due to flow-induced inhomogeneities is modelled via the Rolie–Poly constitutive model. The Rolie–Poly constitutive equation is especially chosen because it expresses, precisely, the shear rheometry of polymer solutions for a large number of strain rates. For the Rolie–Poly approach, we use the two-fluid model wherein the stress dynamics are coupled with concentration equations. We follow a computational analysis approach via an efficient and versatile numerical algorithm. The numerical algorithm is based on the Finite Volume Method (FVM) and it is implemented in the open-source software package, OpenFOAM. The efficiency of our numerical algorithms is enhanced via two possible stabilization techniques, namely; the Log-Conformation Reformulation (LCR) and the Discrete Elastic Viscous Stress Splitting (DEVSS) methodologies. We demonstrate that our stabilized numerical algorithms accurately simulate these complex (shear banded) flows of complex (viscoelastic) fluids. Verification of the shear-banding results via both the Giesekus and Johnson-Segalman models show good agreement with existing literature using the DEVSS technique. A comparison of the Rolie–Poly two-fluid model results with existing literature for the concentration and velocity profiles is also in good agreement.


Author(s):  
M. Mustafa ◽  
A. Mushtaq ◽  
T. Hayat ◽  
A. Alsaedi

Abstract Here we address the influence of heat/mass transfer on MHD axisymmetric viscoelastic fluid flow developed by an elastic sheet stretching linearly in the radial direction. Constitutive relations of Maxwell fluid model are utilized in mathematical formulation of the problem. Non-linear radiation heat flux is factored in the model which accounts for both small and large temperature differences. Chemical reaction effects with modified Arrhenius energy function are analyzed which are not yet explored for viscoelastic fluid flows. Highly accurate numerical computations are performed. Our computations show S-shaped profiles of temperature function in case of sufficiently large temperature differences. Species concentration increases when activation energy for chemical reaction is increased. However, both chemical reaction rate and temperature gradient tend to reduce the solute concentration.


2019 ◽  
Vol 40 (6) ◽  
pp. 763-768 ◽  
Author(s):  
N. A. Khismatullina
Keyword(s):  

Sign in / Sign up

Export Citation Format

Share Document