Direct numerical simulation of particle-laden plane turbulent wall jet and the influence of Stokes number

The distribution and motion of inertial particles in plane turbulent wall jet are investigated using direct numerical simulation, under the assumption of one-way coupling. To our knowledge, this appears to be the first direct numerical simulation of a particle-laden plane turbulent wall jet. It is shown that, in outer part of the wall jet, the behaviour of particles closely resembles that of a free plane jet. Due to the streamwise decay of particle Stokes number, the particle streaks formed in the near wall region of the wall jet are characterized by their intensity variation, which differs significantly from those in the channel flow. The streamwise growth of the particle velocity half-width is approximately equal to that of the fluid velocity half-width and the maximum velocity of particles decays slower than that of fluid due to inertia. The outer scaling can collapse the mean particle velocity in both the inner and outer region for heavier particles. In the buffer region, the particle–fluid velocity difference can be negative or positive depending on the Stokes number since there are two competing effects, namely the memory effect and turbophoresis. In the viscous region, the larger particles are on average faster than fluid and the velocity difference is found to be self-similar depending on outer Stokes number. The near-wall distribution of velocity difference is significantly correlated with the presence of high-momentum particles which are entrained by vortical structures generated in the outer region of the wall jet. These results are useful for environmental and engineering applications.     

Modeling Fluid Flow in Fractured Porous Media with the Interfacial Conditions Between Porous Medium and Fracture

One of the most popular models that has been applied to predict the fluid velocity inside the fracture with impermeable walls is the cubic law. It highlights that the mean flux along the fracture is proportional to the cubic of fracture aperture. However, for a fractured porous medium, the normal and tangential interface conditions between the fracture and porous matrix can change the velocity profile inside the fracture. In this paper, a correction factor is introduced for flow equation along the fracture by imposing the continuity of normal and tangential components of velocity at the interface between the fracture and porous matrix. As a result, the mean velocity inside the fracture depends not only on the fracture aperture, but also on a set of non-dimensional numbers, including the matrix porosity, the ratio of intrinsic permeability of fracture to that of matrix, the wall Reynolds number, and the ratio of normal velocity on one wall to the other. Finally, the introduced correction factor is employed within the extended finite element method, which is widely used for numerical simulation of fluid flow within the fractured porous media. Several numerical results are presented for the fluid flow through a specimen containing single fracture, in order to investigate the deviation from the cubic law in different case studies.

     

Effects of chemical reaction on mixed convection flow of a polar fluid through a porous medium in the presence of internal heat generation

This paper reports a detailed numerical investigation on mixed convection flow of a polar fluid through a porous medium due to the combined effects of thermal and mass diffusion. The energy equation accounts for heat generation or absorption, while the nth order homogeneous chemical reaction between the fluid and the diffusing species is included in the mass diffusion equation. The governing equations of the linear momentum, angular momentum, energy and concentration are obtained in a non-similar form by introducing a suitable group of transformations. The final set of non-similar coupled non-linear partial differential equations is solved using an implicit finite-difference scheme in combination with quasi-linearization technique. The effects of various parameters on the velocity, angular velocity, temperature and concentration fields are investigated. Numerical results for the skin friction coefficient, wall stress of angular velocity, Nusselt number and Sherwood number are also presented.     

A numerical study for three-dimensional viscoelastic flow inspired by non-linear radiative heat flux

Present article examines the three-dimensional flow of upper-convected Maxwell (UCM) fluid over a radiative bi-directional stretching surface. Novel non-linear Rosseland formula for thermal radiation is utilized in the formulation of energy equation. The conventional transformations lead to a strongly non-linear differential system which is treated numerically through Runge–Kutta integration procedure together with the shooting approach. We found that heat transfer rate from the sheet has inverse as well as non-linear relationship with wall to ambient temperature ratio. Moreover an increase in viscoelastic fluid parameter (Deborah number) corresponds to a decrease in the fluid velocity and the boundary layer thickness.     

Exact flow of a third-grade fluid on a porous wall

The flow of a third-grade fluid occupying the space over a wall is studied. At the surface of the wall suction or blowing velocity is applied. By introducing a velocity field, the governing equations are reduced to a non-linear partial differential equation. The resulting equation is analysed analytically using Lie group methods.     

A plane turbulent wall jet on a fully rough surface

The mean velocity field and skin friction characteristics of a plane turbulent wall jet on a smooth and a fully rough surface were studied using Particle Image Velocimetry. The Reynolds number based on the slot height and the exit velocity of the jet was Re = 13,400 and the nominal size of the roughness was k = 0.44 mm. For this Reynolds number and size of roughness element, the flow was in the fully rough regime. The surface roughness results in a distinct change in the shape of the mean velocity profile when scaled in outer coordinates, i.e. using the maximum velocity and outer half-width as the relevant velocity and length scales, respectively. Using inner coordinates, the mean velocity in the lower region of the inner layer was consistent with a logarithmic profile which characterizes the overlap region of a turbulent boundary layer; for the rough wall case, the velocity profile was shifted downward due to the enhanced wall shear stress. For the fully rough flow, the decay rate of the maximum velocity of the wall jet is increased, and the skin friction coefficient is much larger than for the smooth wall case. The inner layer is also thicker for the rough wall case. The effects of surface roughness were observed to penetrate into the outer layer and slightly enhance the spread rate for the outer half-width, which was not observed in most other studies of transitionally rough wall jet flows.     

Turbulent Channel Flow with an Artificial Two-Dimensional Wall Layer

Turbulent flow of an incompressible fluid in a plane channel with parallel walls is considered. The three-dimensional time-dependent Navier-Stokes equations are solved numerically using the spectral finite-difference method. An artificial force which completely suppresses lateral oscillations of the velocity is introduced in the near-wall zone (10 % of the channel half-width in the neighborhood of each wall). Thus, the three-dimensional flow zone, in which turbulent oscillations can develop, is separated from the wall by a fluid layer. It is found that the elimination of three-dimensionality in the neighborhood of the walls leads to a significant reduction in the drag. However, complete laminarization does not occur. The flow in the stream core remains turbulent and can be interpreted as a turbulent flow in a channel with walls located on the boundary of the two-dimensional layer and traveling at the local mean-flow velocity. The oscillations developing inside the two-dimensional layer, which have significant amplitude, distort the flow only in the adjacent zone. Beyond this zone the distributions of the mean characteristics and the structure of instantaneous fields completely correspond to ordinary turbulent flow in a channel with rigid walls. The results obtained confirm the hypothesis of the unimportance of the no-slip boundary conditions for the fluctuating velocity component in the mechanism of onset and self-maintenance of turbulence in wall flows.     

Modeling of Solute Transport in a 3D Rough-Walled Fracture–Matrix System

Fluid flow and solute transport in a 3D rough-walled fracture–matrix system were simulated by directly solving the Navier–Stokes equations for fracture flow and solving the transport equation for the whole domain of fracture and matrix with considering matrix diffusion. The rough-walled fracture–matrix model was built from laser-scanned surface tomography of a real rock sample, by considering realistic features of surfaces roughness and asperity contacts. The numerical modeling results were compared with both analytical solutions based on simplified fracture surface geometry and numerical results by particle tracking based on the Reynolds equation. The aim is to investigate impacts of surface roughness on solute transport in natural fracture–matrix systems and to quantify the uncertainties in application of simplified models. The results show that fracture surface roughness significantly increases heterogeneity of velocity field in the rough-walled fractures, which consequently cause complex transport behavior, especially the dispersive distributions of solute concentration in the fracture and complex concentration profiles in the matrix. Such complex transport behaviors caused by surface roughness are important sources of uncertainty that needs to be considered for modeling of solute transport processes in fractured rocks. The presented direct numerical simulations of fluid flow and solute transport serve as efficient numerical experiments that provide reliable results for the analysis of effective transmissivity as well as effective dispersion coefficient in rough-walled fracture–matrix systems. Such analysis is helpful in model verifications, uncertainty quantifications and design of laboratorial experiments.     

OPERATOR–SPLITTING COMPUTATION OF TURBULENT FLOW IN AN AXISYMMETRIC 180° NARROWING BEND USING SEVERAL k—ϵ MODELS AND WALL FUNCTIONS

This paper describes a finite element implementation of an operator-splitting algorithm for solving transient/steady turbulent flows and presents solutions for the turbulent flow in an axisymmetric 180° narrowing bend, a benchmark problem dealt with at the 1994 WUA-CFD annual meeting. Three k–ϵ based models are used: the standard linear k–ϵ model, a non-linear k–ϵ model and an RNG k–ϵ model. Flow separation after the bend, as observed in the experiment, is predicted by the RNG model and by both the linear and non-linear k–ϵε models with van Driest mixing length wall functions. Good agreement with experimental data of pressure distribution on bending walls is obtained by the present numerical simulation. Results show that there is very little difference between the linear and non-linear k–ϵε models in terms of predicted velocity fields and that the non-linearities mainly affect the distribution of turbulent normal stress and pressure, in analogy to the effect of second-order viscoelastic fluid models on laminar flow. Both the linear and non-linear k–ϵε models fail to predict any flow separation if logarithmic wall functions are used.     

Peristaltic motion of a power-law fluid in an asymmetric channel

The flow of a power-law fluid is investigated in an asymmetric channel caused by the movement of peristaltic waves with the same speed but with different amplitudes and phases on the flexible walls of the channel. The differential equation governing the flow is non-linear and can admit non-unique solutions. There exist two different physically meaningful solutions one satisfying the boundary conditions at the upper wall and the other at the lower wall. The effects of the power-law nature of the fluid on the pumping characteristics and axial velocity are studied in detail.     

Flow and heat transfer over a generalized stretching/shrinking wall problem—Exact solutions of the Navier-Stokes equations

In this paper, we investigate the steady momentum and heat transfer of a viscous fluid flow over a stretching/shrinking sheet. Exact solutions are presented for the Navier-Stokes equations. The new solutions provide a more general formulation including the linearly stretching and shrinking wall problems as well as the asymptotic suction velocity profiles over a moving plate. Interesting non-linear phenomena are observed in the current results including both exponentially decaying solution and algebraically decaying solution, multiple solutions with infinite number of solutions for the flow field, and velocity overshoot. The energy equation ignoring viscous dissipation is solved exactly and the effects of the mass transfer parameter, the Prandtl number, and the wall stretching/shrinking strength on the temperature profiles and wall heat flux are also presented and discussed. The exact solution of this general flow configuration is a rare case for the Navier-Stokes equation.     

On the modelling conditions of mass transfer in porous media presenting capacitance effects by a dispersion-convection equation for the mobile fluid and a diffusion equation for the stagnant fluid

The modelling of mass transfer in porous media presenting capacitance effects by a dispersion-convection equation for the mobile fluid and a diffusion equation for the stagnant fluid has been shown (Piquemal, 1992) to be erroneous in the general case, because it is assumed that the mean concentration of the flowing fluid equals the point concentration at the boundary of the stagnant fluid. This boundary condition cannot be realized. This paper gives the conditions that allows the use of this kind of model with an acceptable approximation. The problem has been solved in the case of two schematic structures: the first is a cylindrical tube with stagnant pockets in its wall, the second a stratified medium. The characteristic lengths of the mobile and immobile domain must obey a criterion in which the porous medium characteristics and the flow velocity appear.     

CALCULATION OF STEADY AND OSCILLATING FLOWS IN TUBES USING A VORTICITY TRANSPORT ALGORITHM

Steady and oscillating axisymmetric tube flows are modelled using a vorticity transport algorithm. The axisymmetric convective –diffusive Navier–Stokes equations are solved using a splitting technique. Axisymmetric ring vortex filaments are introduced on the walls and subsequently convected and diffused throughout the flow field. An axisymmetric equation similar to the Oseen diffusion equation is used to diffuse the ring vortex filaments. Vorticity is reflected from the tube walls using two techniques. Results are presented for the developing Poiseuille flow and for the developed flow in the form of the entrance length and the axial velocity and vorticity profiles. Good agreement is achieved with a finite difference method in the developing region of Poiseuille flow. The developed flow results are compared with the analytical solutions. The developed profiles of velocity and vorticity have errors of less than 0ċ3 per cent for both methods of dealing with reflection of diffusion at the bounding surfaces and similar accuracy is obtained for the velocity profiles in oscillating flow except at the wall. Oscillating flow is produced with a discretized sinusoidal piston motion. Velocity profiles, boundary layer thickness and entrance length are presented for oscillating flow. Good agreement is achieved for low-Womersley-number non-dimensional frequency. At higher values of this parameter, flows are inaccurately simulated, because the number of piston positions used to discretize the piston motion is inversely proportional to the non-dimensional frequency.     

Film thickness measurements in liquid–liquid slug flow regimes

At present there is significant interest in the development of small scale medical diagnostic equipment. These devices offer faster processing times and require smaller sample volumes than equivalent macro scale systems. Although significant attention has been focused upon their outputs, little attention has been devoted to the detailed fluid mechanics that govern the flow mechanisms within these devices. Conventionally, the samples in these small scale devices are segmented into distinct discrete droplets or slugs which are suspended in an organic carrier phase. Separating these slugs from the channel wall is a very thin film of the organic carrier phase.The magnitude of this film is the focus of the present study and the effects of sample slug length and carrier phase fluidic properties on the film are examined over a range of Capillary numbers. A non-intrusive optical technique was used to capture images of the flow from which the magnitude of the film was determined.The experimental results show that the film is not constant along the length of the slug; however above a threshold value for slug length, a region of constant film thickness exists. When compared with existing correlations in the literature, the experimental data showed reasonable agreement with the Bretherton model when the Capillary number was calculated based on the mean two phase flow velocity. However, significant differences were observed when the Capillary number was redefined to account for the mean velocity at the liquid interface, i.e., the mean slug velocity.Analysis of the experimental data revealed that it fell into two distinct flow regimes; a visco-capillary regime and a visco-inertial regime. A modified Taylor expression is presented to estimate the magnitude of the film for flows in the visco-capillary regime while a new model is put forward, based on Capillary and Weber numbers, for flows in the visco-inertial regime. Overall, this study provides some novel insights into parameters, such as aqueous slug length and carrier phase fluidic properties, that affect the thickness of the film in liquid–liquid slug flow regimes.     

Simulation of proppant transport at intersection of hydraulic fracture and natural fracture of wellbores using CFD-DEM

Proppants transport is an advanced technique to improve the hydraulic fracture phenomenon, in order to promote the versatility of gas/oil reservoirs. A numerical simulation of proppants transport at both hydraulic fracture (HF) and natural fracture (NF) intersection is performed to provide a better understanding of key factors which cause, or contribute to proppants transport in HF–NF intersection. Computational fluid dynamics (CFD) in association with discrete element method (DEM) is used to model the complex interactions between proppant particles, host fluid medium and fractured walls. The effect of non-spherical geometry of particles is considered in this model, using the multi-sphere method. All interaction forces between fluid flow and particles are considered in the computational model. Moreover, the interactions of particle–particle and particle–wall are taken into account via Hertz–Mindlin model. The results of the CFD-DEM simulations are compared to the experimental data. It is found that the CFD-DEM simulation is capable of predicting proppant transport and deposition quality at intersections which are in agreement with experimental data. The results indicate that the HF–NF intersection type, fluid velocity and NF aperture affect the quality of blockage occurrence, presenting a new index, called the blockage coefficient which indicates the severity of the blockage.     

Effect of free-stream turbulence on conical vortex dynamics

The influence of Free-Stream Turbulence (FST) on the space–time dynamics of a conical vortex developing along a A-pillar is studied experimentally. Measurements of unsteady wall pressure and velocity by High Speed-Stereo PIV highlight the important effects of turbulence on the mean and instantaneous properties of the vortex. Very significant increases in Reynolds stresses into the vortex region and in wall fluctuating pressure are observed in the presence of FST. In smooth flow, the frequency content of the pressure and velocity fields is very rich with low and high frequency contributions due to the meandering of the vortex and instabilities in the vortex core. Meandering shows, for the different integral length scales and intensities of turbulence tested, a great receptivity to the presence of a FST and we observe a global motion of the vortex structure at low frequency. This frequency is modulated by the value of the integral length scale of the FST. We show that the mean conical structure is a wave guide for the perturbations of the core but that, with FST, the spatio-temporal evolution of the envelope overwhelms the intrinsic instability of the vortex core observed in smooth flow.     

Heat transfer in laminar entrance region of a flat channel for the temperature boundary condition of the third kind

To explore the influence of the developing flow in a flat channel on the laminar forced convection heat transfer, the non-linear momentum and the linear energy equation are solved successively by employing the Galerkin-Kantorowich method of variational calculus. Assuming constant fluid properties, negligible axial diffusion and temperature boundary condition of the third kind, semi analytic solution for velocity and temperature is derived. The local Nusselt numbers are tabulated for various values of Biot and Prandtl number.     

One-way coupled fluid–structure interaction of gas–liquid slug flow in a horizontal pipe: Experiments and simulations

Pipelines conveying a multiphase mixture must withstand the cyclic induced stresses that occur due to the alternating motion of gas pockets and liquid slugs. Few previous studies have considered gas–liquid slug flow and the associated fluid–structure interaction problems. In this study, experimental and numerical techniques were adopted to simulate and analyze the two-phase slug flow and the associated stresses in the pipe structure. In the numerical simulation, a one-way coupled fluid–structure framework was developed to explore the slug flow interaction with a horizontal pipe assembly under various superficial gas and liquid velocities. A modified Volume of Fluid and finite element methods were utilized to model the fluid and structure domains. The file-based coupling technique was adopted to execute the coupling mechanism. By contrast, slug characteristics were measured experimentally, while Bi-axial strain gauges were used to capture time-varying strain signals. Excellent agreements between the predicted and measured stress results were achieved with a maximum error of 10.2 %. It was found that at constant superficial liquid velocity, the maximum induced stresses on the pipe wall increased with increasing the slug length and slug velocity. While for the slug frequency, the maximum principal stresses decreased with increasing the slug frequency.     

Wave propagation in an elastic tube: a numerical study

In this paper, a fluid–wall interaction model, called the elastic tube model, is introduced to investigate wave propagation in an elastic tube and the effects of different parameters. The unsteady flow was assumed to be laminar, Newtonian and incompressible, and the vessel wall to be linear-elastic, isotropic and incompressible. A fluid–wall interaction scheme is constructed using a finite element method. The results demonstrate that the elastic tube plays an important role in wave propagation. It is shown that there is a time delay between the velocity waveforms at two different locations and that the peak velocity increases while the low velocity decreases in the elastic tube model, contrary to the rigid tube model where velocity waveforms overlap each other. Compared with the elastic tube model, the increase of the wall thickness makes wave propagation faster and the time delay cannot be observed clearly, however, the velocity amplitude is reduced slightly due to the decrease of the internal radius. The fluid–wall interaction model simulates wave propagation successfully and can be extended to study other mechanical properties considering complicated geometrical and material factors.     

Hydrodynamics of viscous fluid through porous slit with linear absorption

The exact solutions for the viscous fluid through a porous slit with linear ab-sorption are obtained. The Stokes equation with non-homogeneous boundary conditions is solved to get the expressions for the velocity components, pressure distribution, wall shear stress, fractional absorption, and leakage flux. The volume flow rate and mean flow rate are found to be useful in obtaining a convenient form of the longitudinal velocity component and pressure difference. The points of the maximum velocity components for a fixed axial distance are identified. The value of the linear absorption parameter is ran-domly chosen, and the rest available data of the rat kidney to the tabulate pressure drop and fractional absorption are incorporated. The effects of the linear absorption, uniform absorption, and flow rate parameters on the flow properties are discussed by graphs. It is found that forward flow occurs only if the volume flux per unit width is greater than the absorption velocity throughout the length of the slit, otherwise back flow may occur. The leakage flux increases with the increase in the linear absorption parameter. Streamlines are drawn to help the analysis of the flow behaviors during the absorption of the fluid flow through the renal tubule and purification of blood through an artificial kidney.