scholarly journals Turbulent Bubble-Laden Channel Flow of Power-Law Fluids: A Direct Numerical Simulation Study

Fluids ◽  
2021 ◽  
Vol 6 (1) ◽  
pp. 40
Author(s):  
Felix Bräuer ◽  
Elias Trautner ◽  
Josef Hasslberger ◽  
Paolo Cifani ◽  
Markus Klein

The influence of non-Newtonian fluid behavior on the flow statistics of turbulent bubble-laden downflow in a vertical channel is investigated. A Direct Numerical Simulation (DNS) study is conducted for power-law fluids with power-law indexes of 0.7 (shear-thinning), 1 (Newtonian) and 1.3 (shear-thickening) in the liquid phase at a gas volume fraction of 6%. The flow is driven downward by a constant volumetric flow rate corresponding to a friction Reynolds number of Reτ≈127.3. The Eötvös number is varied between Eo=0.3125 and Eo=3.75 in order to investigate the influence of quasi-spherical as well as wobbling bubbles and thus the interplay of the bubble deformability with the power-law behavior of the liquid bulk. The resulting first- and second-order fluid statistics, i.e., the gas fraction, mean velocity and velocity fluctuation profiles across the channel, show clear trends in reply to varying power-law indexes. In addition, it was observed that the bubble oscillations increase with decreasing power-law index. In the channel core, the bubbles significantly increase the dissipation rate, which, in contrast to its behavior at the wall, shows similar orders of magnitude for all power-law indexes.

2017 ◽  
Vol 826 ◽  
pp. 918-941 ◽  
Author(s):  
A. Bougouin ◽  
L. Lacaze ◽  
T. Bonometti

Experiments on the collapse of non-colloidal and neutrally buoyant particles suspended in a Newtonian fluid column are presented, in which the initial volume fraction of the suspension $\unicode[STIX]{x1D719}$, the viscosity of the interstitial fluid $\unicode[STIX]{x1D707}_{f}$, the diameter of the particles $d$ and the mixing protocol, i.e. the initial preparation of the suspension, are varied. The temporal evolution of the slumping current highlights two main regimes: (i) an inertial-dominated regime followed by (ii) a viscous-dominated regime. The inertial regime is characterized by a constant-speed slumping which is shown to scale as in the case of a classical inertial dam-break. The viscous-dominated regime is observed as a decreasing-speed phase of the front evolution. Lubrication models for Newtonian and power-law fluids describe most of situations encountered in this regime, which strongly depends on the suspension parameters. The temporal evolution of the propagating front is used to extract the rheological parameters of the fluid models. At the early stages of the viscous-dominated regime, a constant effective shear viscosity, referred to as an apparent Newtonian viscous regime, is found to depend only on $\unicode[STIX]{x1D719}$ and $\unicode[STIX]{x1D707}_{f}$ for each mixing protocol. The obtained values are shown to be well fitted by the Krieger–Dougherty model whose parameters involved, say a critical volume fraction $\unicode[STIX]{x1D719}_{m}$ and the exponent of divergence, depend on the mixing protocol, i.e. the microscale interaction between particles. On a longer time scale which depends on $\unicode[STIX]{x1D719}$, the front evolution is shown to slightly deviate from the apparent Newtonian model. In this apparent non-Newtonian viscous regime, the power-law model, indicating both shear-thinning and shear-thickening behaviours, is shown to be more appropriate to describe the front evolution. The present experiments indicate that the mixing protocol plays a crucial role in the selection of a shear-thinning or shear-thickening type of collapse, while the particle diameter $d$ and volume fraction $\unicode[STIX]{x1D719}$ play a significant role in the shear-thickening case. In all cases, the normalized effective consistency of the power-law fluid model is found to be a unique function of $\unicode[STIX]{x1D719}$. Finally, an apparent viscoplastic regime, characterized by a finite length spreading reached at finite time, is observed at high $\unicode[STIX]{x1D719}$. This regime is mostly observed for volume fractions larger than $\unicode[STIX]{x1D719}_{m}$ and up to a volume fraction $\unicode[STIX]{x1D719}_{M}$ close to the random close packing fraction at which the initial column remains undeformed on opening the gate.


2018 ◽  
Vol 857 ◽  
pp. 270-290 ◽  
Author(s):  
Josef Hasslberger ◽  
Markus Klein ◽  
Nilanjan Chakraborty

This paper presents a detailed investigation of flow topologies in bubble-induced two-phase turbulence. Two freely moving and deforming air bubbles that have been suspended in liquid water under counterflow conditions have been considered for this analysis. The direct numerical simulation data considered here are based on the one-fluid formulation of the two-phase flow governing equations. To study the development of coherent structures, a local flow topology analysis is performed. Using the invariants of the velocity gradient tensor, all possible small-scale flow structures can be categorized into two nodal and two focal topologies for incompressible turbulent flows. The volume fraction of focal topologies in the gaseous phase is consistently higher than in the surrounding liquid phase. This observation has been argued to be linked to a strong vorticity production at the regions of simultaneous high fluid velocity and high interface curvature. Depending on the regime (steady/laminar or unsteady/turbulent), additional effects related to the density and viscosity jump at the interface influence the behaviour. The analysis also points to a specific term of the vorticity transport equation as being responsible for the induction of vortical motion at the interface. Besides the known mechanisms, this term, related to surface tension and gradients of interface curvature, represents another potential source of turbulence production that lends itself to further investigation.


2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ($\mathit{Re}_{{\it\tau}}$) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Kármán constant ${\it\kappa}=0.384\pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in wavenumber $(k)$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.


2019 ◽  
Vol 862 ◽  
pp. 781-815 ◽  
Author(s):  
Y. Kuwata ◽  
Y. Kawaguchi

Lattice Boltzmann direct numerical simulation of turbulent open-channel flows over randomly distributed hemispheres at $Re_{\unicode[STIX]{x1D70F}}=600$ is carried out to reveal the influence of roughness parameters related to a probability density function of rough-surface elevation on turbulence by analysing the spatial and Reynolds- (double-) averaged Navier–Stokes equation. This study specifically concentrates on the influence of the root-mean-square roughness and the skewness, and profiles of turbulence statistics are compared by introducing an effective wall-normal distance defined as a wall-normal integrated plane porosity. The effective distance can completely collapse the total shear stress outside the roughness sublayer, and thus the similarity of the streamwise mean velocity is clearer by introducing the effective distance. In order to examine the influence of the root-mean-square roughness and the skewness on dynamical effects that contribute to an increase in the skin friction coefficient, the triple-integrated double-averaged Navier–Stokes equation is analysed. The main contributors to the skin friction coefficient are found to be turbulence and drag force. The turbulence contribution increases with the root-mean-square roughness and/or the skewness. The drag force contribution, on the other hand, increases in particular with the root-mean-square roughness whereas an increase in the skewness does not increase the drag force contribution because it does not necessarily increase the surface area of the roughness elements. The contribution of the mean velocity dispersion induced by spatial inhomogeneity of the rough surfaces substantially increases with the root-mean-square roughness. A linear correlation is confirmed between the root-mean-square roughness and the equivalent roughness while the equivalent roughness monotonically increases with the skewness. A new correlation function based on the root-mean-square roughness and the skewness is developed with the available experimental and direct numerical simulation data, and it is confirmed that the developed correlation reasonably predicts the equivalent roughness of various types of real rough surfaces.


Author(s):  
Shengbiao Zhang ◽  
Zhandong Wang ◽  
Yan Yan ◽  
Guifang Sun ◽  
Zhonghua Ni

A fluid simulation was conducted on the flow states of water and drainage gas from the drainage cover in underwater local dry laser welding. The results of the simulation of gas volume fraction are obtained. Drainage cover was developed based on the drainage cover model. Experiments of drainage in this machined drainage cover was carried out, and were compared with the numerical simulation. The results show high similarity, which verifies the feasibility. Then, innovative structure featuring an uneven buffer around the exit of drainage cover are proposed to improve the performance of drainage. A group of 3D models of fluid domain related to differently characterized drainage covers was designed and the flow states were simulated. They are N-model with a narrow outlet, S-model with a smooth buffer added additionally based on N-model, U-model with an uneven buffer compared with S-model, and W-model with a wide outlet and the diameter of outlet twice than N-model. The gas volume fracture distributions were recorded and in comparison with each other. The results indicate that U-model with the uneven buffer realized the deepest drainage depth.


1993 ◽  
Vol 115 (2) ◽  
pp. 330-333 ◽  
Author(s):  
Ali J. Chamkha

Equations governing flow of a particulate suspension exhibiting finite volume fraction in non-Newtonian power-law fluids are developed and applied to the problem of unsteady flow past an infinite porous flat plate with suction. Numerical results for small volume fraction for the displacement thicknesses for both phases and the skin-friction coefficient for the fluid phase are obtained using an implicit finite difference scheme and presented graphically to elucidate interesting features of the solutions.


2014 ◽  
Vol 751 ◽  
pp. 184-215
Author(s):  
Liyan Yu ◽  
John Hinch

AbstractWe study the solitary wave solutions in a thin film of a power-law fluid coating a vertical fibre. Different behaviours are observed for shear-thickening and shear-thinning fluids. For shear-thickening fluids, the solitary waves are larger and faster when the reduced Bond number is smaller. For shear-thinning fluids, two branches of solutions exist for a certain range of the Bond number, where the solitary waves are larger and faster on one and smaller and slower on the other as the Bond number decreases. We carry out an asymptotic analysis for the large and fast-travelling solitary waves to explain how their speeds and amplitudes change with the Bond number. The analysis is then extended to examine the stability of the two branches of solutions for the shear-thinning fluids.


Sign in / Sign up

Export Citation Format

Share Document