Elementary Numerical Analysis of Wet Foam Formation and Study of Its Flow Structures and Physical Behavior

2021 ◽  
Author(s):  
Ebrahim Shirani ◽  
Sima Nasirzade ◽  
Fethi Aloui

Abstract The purpose of this study is to analyze the flow of wet foam and to study the effect of volume fraction, velocity and surface tension and other physical parameters on the foam flow. The most numerical researches done in this area are for single-phase flows. The numerical simulation in this study is the first simulation in the foam flow, in which both the bubbles and the water are simulated as two-phase flow. In this study, fluid containing a surfactant and bubbles are flowing in a duct. The dimensions of the duct cross section is 15 × 60 in millimeters. The numerical solution is performed for three Reynolds numbers of 50, 100 and 1000, three volume fractions of 48, 41 and 28, and three Weber numbers of 0.405, 0.27 and 0.203 (27 different modes), and the effect of the above parameters on the flow behavior and its physical properties have been investigated. It was found that in foam flow, the velocity fluctuations, due to the movement of bubbles in the flow, is in the order of magnitude of the mean velocity. The same is true for wall shear stress. By increasing the Reynolds number, the pressure loss increases, the magnitude of the velocity fluctuations decreases and the frequency of the velocity fluctuations increases. By increasing the Weber number, the pressure loss and the magnitude of the velocity fluctuations decrease and the mean shear stress increases. By increasing the foam quality, pressure loss increases, the mean shear stress and the magnitude of the velocity fluctuations decrease and its frequency increases. And the phenomenon of coalescence causes a sudden increase in momentum speed.

Author(s):  
Shinji Honami ◽  
Wataru Tsuboi ◽  
Takaaki Shizawa

This paper presents the effect of flame dome depth on the total pressure performance and flow behavior in a sudden expansion region of the combustor diffuser without flow entering the dome head. The mean velocity and turbulent Reynolds stress profiles in the sudden expansion region were measured by a Laser Doppler Velocitmetry (LDV) system. The experiments show that total pressure loss is increased, when flame dome depth is increased. Installation of an inclined combuster wall in the sudden expansion region is suggested from the viewpoint of a control of the reattaching flow. The inclined combustor wall is found to be effective in improvement of the diffuser performance. Better characteristics of the flow rate distribution into the branched channels are obtained in the inclined wall configuration, even if the distorted velocity profile is provided at the diffuser inlet.


2015 ◽  
Vol 770 ◽  
pp. 210-246 ◽  
Author(s):  
M. Mehrabadi ◽  
S. Tenneti ◽  
R. Garg ◽  
S. Subramaniam

Gas-phase velocity fluctuations due to mean slip velocity between the gas and solid phases are quantified using particle-resolved direct numerical simulation. These fluctuations are termed pseudo-turbulent because they arise from the interaction of particles with the mean slip even in ‘laminar’ gas–solid flows. The contribution of turbulent and pseudo-turbulent fluctuations to the level of gas-phase velocity fluctuations is quantified in initially ‘laminar’ and turbulent flow past fixed random particle assemblies of monodisperse spheres. The pseudo-turbulent kinetic energy $k^{(f)}$ in steady flow is then characterized as a function of solid volume fraction ${\it\phi}$ and the Reynolds number based on the mean slip velocity $\mathit{Re}_{m}$. Anisotropy in the Reynolds stress is quantified by decomposing it into isotropic and deviatoric parts, and its dependence on ${\it\phi}$ and $Re_{m}$ is explained. An algebraic stress model is proposed that captures the dependence of the Reynolds stress on ${\it\phi}$ and $Re_{m}$. Gas-phase velocity fluctuations in freely evolving suspensions undergoing elastic and inelastic particle collisions are also quantified. The flow corresponds to homogeneous gas–solid systems, with high solid-to-gas density ratio and particle diameter greater than dissipative length scales. It is found that for the parameter values considered here, the level of pseudo-turbulence differs by only 15 % from the values for equivalent fixed beds. The principle of conservation of interphase turbulent kinetic energy transfer is validated by quantifying the interphase transfer terms in the evolution equations of kinetic energy for the gas-phase and solid-phase fluctuating velocity. It is found that the collisional dissipation is negligible compared with the viscous dissipation for the cases considered in this study where the freely evolving suspensions attain a steady state starting from an initial condition where the particles are at rest.


2013 ◽  
Vol 136 (2) ◽  
Author(s):  
A. M. Rashad ◽  
S. Abbasbandy ◽  
Ali J. Chamkha

In recent years, nanofluids have attracted attention as a new generation of heat transfer fluids in building heating, heat exchangers, plants, and automotive cooling applications because of their excellent thermal performance. Various benefits of the application of nanofluids include improved heat transfer, heat transfer system size reduction, minimal clogging, microchannel cooling, and miniaturization of systems. In this paper, a study of steady, laminar, natural convection boundary-layer flow adjacent to a vertical cylinder embedded in a thermally stratified nanofluid-saturated non-Darcy porous medium is investigated. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis, and a generalized porous media model, which includes inertia and boundary effects, is employed. The cylinder surface is maintained at a constant nanoparticles volume fraction, and the wall temperature is assumed to vary with the vertical distance according to the power law form. The resulting governing equations are nondimensionalized and transformed into a nonsimilar form and then solved by Keller box method. A comparison is made with the available results in the literature, and our results are in very good agreement with the known results. A parametric study of the physical parameters is made, and a representative set of numerical results for the velocity, temperature, and volume fraction, as well as local shear stress and local Nusselt and Sherwood numbers, are presented graphically. The salient features of the results are analyzed and discussed. The results indicate that, when the buoyancy ratio or modified Grashof number increases, all of the local shear stress, local Nusselt number, and the local Sherwood number enhance while the opposite behaviors are predicted when the thermophoresis parameter increases. Moreover, increasing the value of the surface curvature parameter leads to increases in all of the local shear stress and the local Nusselt and Sherwood numbers while the opposite behaviors are obtained when either of the thermal stratification parameter or the boundary effect parameter increases.


Author(s):  
Rogelio Chovet ◽  
Fethi Aloui ◽  
Laurent Keirsbulck

Some industrial processes are associated with the flow of aqueous foams inside horizontal channels. Examples are found in the oil, food and cosmetic industries. This type of flow presents an important pressure loss, originated from the shear stress exerted by the channel walls. Foam flow is one of the most complex fluids. In a macroscopic point of view, the physical-chemical interaction between the bubbles can be related to some non-Newtonian models (Bingham law, power law, etc.) or an apparent viscosity. These last can represent the internal deformations of fluid elements when shear stress is applied. An experimental facility able to create this type of flow is not so easy to design. Many parameters must be taken into consideration. So, Computational Fluid Dynamics (CFD) constitutes an ideal technique for analyzing this kind of problem. The aim of this study is to validate the use of Computational Fluid Dynamics in order to correctly predict the pressure losses and the velocity fields of a foam flowing through a straight channel and singularities (fence and half-sudden expansion). Simulations for a realistic scenario: two-phase flow, change in the surface tension, bubble size, were undertaken. Obtained results showed that simulations are not able to accurately reproduce for such a complex fluid, the important aspects of this study, such as the pressure losses and the velocity fields. Therefore, an approximation to a Bingham fluid was made. For a foam flow quality of 70% and a velocity of 2 cm/s, the numerical results are justified by experimental evidence. Experiments have been done and predictions for the flow behavior are extrapolated. Results show that the software is able to recreate the behavior of foam flow through a straight channel and singularities. However, this approach is extremely sensitive to the choice of several parameters, like the apparent viscosity, the yield stress, the viscosity consistence, etc.


Author(s):  
Francesco Soranna ◽  
Yi-Chih Chow ◽  
Oguz Uzol ◽  
Joseph Katz

Stereoscopic PIV measurements examine the flow structure and turbulence within a rotor near wake located in a non-uniform field generated by a row of Inlet Guide Vanes (IGVs). The experiments are performed in a refractive index matched facility that provides unobstructed view of the entire flow field. The data are acquired at 10 closely spaced radial planes located near mid-span, which enable measurements of all the components of the mean strain rate and Reynolds stress tensors. Chopping and variations of advection speed of the upstream IGV wakes, as they pass along the rotor blade, create a non-uniform flow that shears the rotor wake. However, the phase averaged flow at mid span remains almost two-dimensional. Due to the overwhelming effects of the non-uniform strain field, the presently observed trends of the Reynolds stresses within the sheared wake differ from those measured in previous studies of curved wakes. The axial velocity fluctuations increase along the suction/outer side of the wake, while the other components decay. On the pressure/inner part of the wake the circumferential velocity fluctuations are higher. The shear stress has a complex distribution, but is also higher on the suction side. To explain these trends the stresses and production rates are examined in coordinate systems aligned with the principal strain directions. As expected, the production is high along the compressive directions and low, even negative, in the extensive directions. Accordingly, the compressed normal stress component increases along the wake, while the extended component decays. The in-plane shear stress component and its associated production remain very high in the principal coordinate system of the strain. Projecting the stresses back to the laboratory coordinate system explains the observed inhomogeneous anisotropic distribution of Reynolds stresses within the kinked wake.


2011 ◽  
Vol 687 ◽  
pp. 41-71 ◽  
Author(s):  
Partha S. Goswami ◽  
V. Kumaran

AbstractThe particle and fluid velocity fluctuations in a turbulent gas–particle suspension are studied experimentally using two-dimensional particle image velocimetry with the objective of comparing the experiments with the predictions of fluctuating force simulations. Since the fluctuating force simulations employ force distributions which do not incorporate the modification of fluid turbulence due to the particles, it is of importance to quantify the turbulence modification in the experiments. For experiments carried out at a low volume fraction of $9. 15\ensuremath{\times} 1{0}^{\ensuremath{-} 5} $ (mass loading is 0.19), where the viscous relaxation time is small compared with the time between collisions, it is found that the gas-phase turbulence is not significantly modified by the presence of particles. Owing to this, quantitative agreement is obtained between the results of experiments and fluctuating force simulations for the mean velocity and the root mean square of the fluctuating velocity, provided that the polydispersity in the particle size is incorporated in the simulations. This is because the polydispersity results in a variation in the terminal velocity of the particles which could induce collisions and generate fluctuations; this mechanism is absent if all of the particles are of equal size. It is found that there is some variation in the particle mean velocity very close to the wall depending on the wall-collision model used in the simulations, and agreement with experiments is obtained only when the tangential wall–particle coefficient of restitution is 0.7. The mean particle velocity is in quantitative agreement for locations more than 10 wall units from the wall of the channel. However, there are systematic differences between the simulations and theory for the particle concentrations, possibly due to inadequate control over the particle feeding at the entrance. The particle velocity distributions are compared both at the centre of the channel and near the wall, and the shape of the distribution function near the wall obtained in experiments is accurately predicted by the simulations. At the centre, there is some discrepancy between simulations and experiment for the distribution of the fluctuating velocity in the flow direction, where the simulations predict a bi-modal distribution whereas only a single maximum is observed in the experiments, although both distributions are skewed towards negative fluctuating velocities. At a much higher particle mass loading of 1.7, where the time between collisions is smaller than the viscous relaxation time, there is a significant increase in the turbulent velocity fluctuations by ${\ensuremath{\sim} }1$–2 orders of magnitude. Therefore, it becomes necessary to incorporate the modified fluid-phase intensity in the fluctuating force simulation; with this modification, the mean and mean-square fluctuating velocities are within 20–30 % of the experimental values.


2021 ◽  
Vol 249 ◽  
pp. 03039
Author(s):  
James Jenkins ◽  
Diego Berzi

We first phrase a boundary-value problem for a dense, steady, fully-developed, gravitational flow of identical inelastic spheres over in inclined bumpy base in the absence of sidewalls. We then obtain approximate analytical solutions for the profiles of the solid volume fraction, the strength of the velocity fluctuations, and the mean velocity of the flow. We compare these with those obtained in numerical solutions of the exact equations.


2001 ◽  
Vol 429 ◽  
pp. 307-342 ◽  
Author(s):  
ROBERTO ZENIT ◽  
DONALD L. KOCH ◽  
ASHOK S. SANGANI

Experiments were performed in a vertical channel to study the behaviour of a monodisperse bubble suspension for which the dual limit of large Reynolds number and small Weber number was satisfied. Measurements of the liquid-phase velocity fluctuations were obtained with a hot-wire anemometer. The gas volume fraction, bubble velocity, bubble velocity fluctuations and bubble collision rate were measured using a dual impedance probe. Digital image analysis was performed to quantify the small polydispersity of the bubbles as well as the bubble shape.A rapid decrease in bubble velocity with bubble concentration in very dilute suspensions is attributed to the effects of bubble–wall collisions. The more gradual subsequent hindering of bubble motion is in qualitative agreement with the predictions of Spelt & Sangani (1998) for the effects of potential-flow bubble–bubble interactions on the mean velocity. The ratio of the bubble velocity variance to the square of the mean is O(0.1). For these conditions Spelt & Sangani predict that the homogeneous suspension will be unstable and clustering into horizontal rafts will take place. Evidence for bubble clustering is obtained by analysis of video images. The fluid velocity variance is larger than would be expected for a homogeneous suspension and the fluid velocity frequency spectrum indicates the presence of velocity fluctuations that are slow compared with the time for the passage of an individual bubble. These observations provide further evidence for bubble clustering.


2018 ◽  
Vol 6 (4) ◽  
Author(s):  
Varit Kunopagarnwong ◽  
Thongchai Rohitatisha Srinophakun ◽  
Wisarut Manasthammakul

The flow behavior of liquid-solid particles in mixing tanks using a modified Rushton impeller, called a Rushton V-cut impeller, was studied. Both the Rushton and Rushton V-cut impellers were compared at a 300 rpm stirring speed and a 10 % wt solid concentration. Hydrodynamic behaviors, such as solid volume fraction, velocity, pressure, and shear stress, in both the Rushton and Rushton V-cut impellers were investigated. Computational fluid dynamics (CFD) software able to understanding of hydrodynamics of stirring liquids which contain solid particles. CFD programme was to predict the mixing flow of the highly viscous system. Therefore, the present work was carried out using CFD software with the Eulerian-Eulerian approach with a turbulent k-ε model. The simulation of mixing tanks was consisted of moving and stationary zones by using moving references frame method or MFR. The results were observed that Rushton V-cut can dramatically reduce pressure up to 20% and the shear stress up to 64.38% while keeping the liquid-solid mixing at a considerable degree. Therefore, this design can reduce the power consumption.


1977 ◽  
Vol 82 (4) ◽  
pp. 705-723 ◽  
Author(s):  
T. R. Heidrick ◽  
S. Banerjee ◽  
R. S. Azad

This paper is the second of a pair describing two-point velocity measurements in fully developed pipe flow. A method of processing hot-film anemometer signals to identify intervals of high energy production (‘bursts’) in wall turbulence is presented. The method uses filtered cross-stream spatial derivatives of the axial velocity fluctuations. It is demonstrated to be more sensitive to ‘bursts’ than several other methods of indentification. The bursts identified in this manner are shown to have similar characteristics to those observed in visual studies.The technique has been applied to the wall region of turbulent pipe flow. Mean burst rates have been obtained at various distances from the wall for three Reynolds numbers. It is shown that the mean burst rate cannot be reliably obtained from a previously used technique based on the autocorrelation of the axial velocity fluctuations.On the basis of our experiments, the mean burst rate and the turbulent shear stress have been found to vary similarly with distance from the wall. In the region near the wall where the shear stress is constant the mean burst rate is independent of the kinematic viscosity.Some characteristics of the velocity fluctuations during burst intervals have been studied. All the bursts began with a relative minimum in the axial velocity fluctuations followed by a peak in the cross-stream spatial derivative. A second peak always occurred midway through the burst. The sequence of events is somewhat similar to that in the last stage of laminar-to-turbulent transition.


Sign in / Sign up

Export Citation Format

Share Document