Constricted channel flow with different cross-section shapes

Pressure driven steady flow through a uniform circular channel containing a constricted portion is a common problem considering physiological flows such as underlying human speech sound production. The influence of the constriction’s cross-section shape (circle, ellipse, circular sector) on the flow within and downstream from the constriction is experimentally quantified. An analytical boundary layer flow model is proposed which takes into account the hydraulic diameter of the cross-section shape. Comparison of the model outcome with experimental and three-dimensional numerically simulated flow data shows that the pressure distribution within the constriction can be modeled accurately so that the model is of interest for analytical models of fluid–structure interaction without the assumption of two-dimensional flow

Conforming Chebyshev spectral collocation methods for the solution of laminar flow in a constricted channel

The numerical simulation of steady planar two-dimensional, laminar flow of an incompressible fluid through an abruptly contracting channel using spectral domain decomposition methods is described. The key features of the method are the decomposition of the flow region into a number of rectangular subregions and spectral approximations which are pointwise C(1) continuous across subregion interfaces. Spectral approximations to the solution are obtained for Reynolds numbers in the range 0 to 500. The size of the salient corner vortex decreases as the Reynolds number increases from 0 to around 45. As the Reynolds number is increased further the vortex grows slowly. A vortex is detected downstream of the contraction at a Reynolds number of around 175 that continues to grow as the Reynolds number is increased further

Discrete element study of liquid-solid slurry flows through constricted channels

Discrete element model is used to simulate the flow of liquid-granule mixtures in an inclined channel containing a linear contraction. All the relevant particle/particle and particle/fluid interactions are included in the numerical model. The presence of the contraction induces different steady morphologies of the solid phase or the mixture depending on whether closed or open channels are used. These flows behave quite differently depending on the upstream Froude number and the contraction size ratio. The model is first validated by comparing with the existing results for dry granular (glass particles) chute flows (Vreman et al., 2007). Then simulations of a chute of glass particles in water flowing in a closed channel are compared to the dry granular case. With the same solid flux at the inlet, the hydrodynamic forces in the liquid-solid mixture induce higher particle solid volume fractions in the part of the flow containing the solid phase. The streamwise particle velocity (resp. depth of the solid phase) has the same evolution along the channel with smaller (larger) values than in the dry granular flow case

Critical point network for drainage between rough surfaces

In this paper, we present a network method for computing two-phase flows between two rough surfaces with significant contact areas. Low-capillary number drainage is investigated here since one-phase flows have been previously investigated in other contributions. An invasion percolation algorithm is presented for modeling slow displacement of a wetting fluid by a non wetting one between two rough surfaces. Short-correlated Gaussian process is used to model random rough surfaces.The algorithm is based on a network description of the fracture aperture field. The network is constructed from the identification of critical points (saddles and maxima) of the aperture field. The invasion potential is determined from examining drainage process in a flat mini-channel. A direct comparison between numerical prediction and experimental visualizations on an identical geometry has been performed for one realization of an artificial fracture with a moderate fractional contact area of about 0.3. A good agreement is found between predictions and observations

Comparison of computations of asymptotic flow models in a constricted channel

International audienceWe aim at comparing computations with asymptotic models issued from incom- pressible Navier-Stokes at high Reynolds number: the Reduced Navier-Stokes/Prandtl (RNS/P) equations and the Double Deck (DD) equations. We treat the case of the steady two dimensional flow in a constricted pipe. In particular, finite differences and finite element solvers are compared for the RNS/P equations. It results from this study that the two codes compare well. Numerical examples also illustrate the interest of these asymptotic models as well as the flexibility of the finite element solver

Large-eddy simulation of physiological pulsatile flow through a channel with double constriction

Pulsatile flow in a 3D model of arterial double stenoses is investigated using a large eddy simulation (LES) technique. The computational domain that has been chosen is a simple channel with a biological-type stenosis formed eccentrically on the top wall. The pulsation was generated at the inlet using the first four harmonics of the Fourier series of the pressure pulse. The flow Reynolds numbers, which are typically suitable for a large human artery, are chosen in the present work. In LES, a top-hat spatial grid-filter is applied to the Navier–Stokes equations of motion to separate the large-scale flows from the sub-grid scale (SGS). The large-scale flows are then resolved fully while the unresolved SGS motions are modelled using a localized dynamic model. It is found that the narrowing of the channel causes the pulsatile flow to undergo a transition to a turbulent condition in the downstream region; as a consequence, a severe level of turbulent fluctuations is achieved in these zones. Transitions to turbulent of the pulsatile flow in the post stenosis are examined through the various numerical results, such as velocity, streamlines, wall pressure, shear stresses and root mean square turbulent fluctuations

Vortex dynamics under pulsatile flow in axisymmetric constricted tubes

An improved understanding of how vortices develop and propagate under pulsatile flow can shed important light on the mixing and transport processes including the transition to turbulent regime occurring in such systems. For example, the characterization of pulsatile flows in obstructed artery models serves to encourage research into flow-induced phenomena associated with changes in morphology, blood viscosity, wall elasticity and flow rate. In this work, an axisymmetric rigid model was used to study the behaviour of the flow pattern with varying constriction degree ($d_0$), mean Reynolds number ($\bar{Re}$) and Womersley number ($\alpha$). Velocity fields were acquired experimentally using Digital Particle Image Velocimetry and generated numerically. For the acquisition of data, $\bar{Re}$ was varied from 385 to 2044, $d_0$ was 1.0 cm and 1.6 cm, and $\alpha$ was varied from 17 to 33 in the experiments and from 24 to 50 in the numerical simulations. Results for the considered Reynolds number, showed that the flow pattern consisted of two main structures: a central jet around the tube axis and a recirculation zone adjacent to the inner wall of the tube, where vortices shed. Using the vorticity fields, the trajectory of vortices was tracked and their displacement over their lifetime calculated. The analysis led to a scaling law equation for the maximum vortex displacement as a function of a dimensionless variable dependent on the system parameters Re and $\alpha$

Taylor-Goertler instabilities of Tollmien-Schlichting waves and other flows governed by the interactive boundary layer equations

The Taylor-Gortler vortex instability equations are formulated for steady and unsteady interacting boundary layer flows of the type which arise in triple-deck theory. The effective Gortler number is shown to be a function of the all shape in the boundary layer and the possibility of both steady and unsteady Taylor-Gortler modes exists. As an example the steady flow in a symmetrically constricted channel is considered and it is shown that unstable Gortler vortices exist before the boundary layers at the wall develop the Goldstein singularity. As an example of an unsteady spatially varying basic state the instability of high frequency large amplitude Tollmien-Schlichting waves in a curved channel were considered. It is shown that they are unstable in the first Stokes layer stage of the hierarchy of nonlinear states. The Tollmien-Schlichting waves are shown to be unstable in the presence of both convex and concave curvature