High velocity flow through fractured and porous media: the role of flow non-periodicity

In the present paper we examine the evolution of the macroscopic flow law in a crenellated channel, representing an element of fractured or porous medium and in function of the Reynolds number Re. A numerical analysis based on the Navier–Stokes equations is applied. We focus on the influence of the flow periodicity or non-periodicity upon the macroscopic law. The physical explanation of the non-linear deviation from Darcy's law is still an issue, as the Ergun–Forchheimer law admitted for high Reynolds numbers comes up against some theoretical problems. In the periodic case, three non-linear flow regimes were revealed: a cubic flow with respect to velocity at low Re, an intermediate non-quadratic law, and a self-similar mode independent of Re at very high Re. The Forchheimer law is not confirmed. The case of a non-periodic flow clearly highlights the link between the flow non-periodicity and the quadratic law. The quadratic deviation becomes all the more important as the non-periodicity degree is high.     

A hybrid RANS/LES simulation of high-Reynolds-number channel flow using additional filtering at the interface

A hybrid method combining large eddy simulation (LES) with the Reynolds-averaged Navier-Stokes (RANS) equation is used to simulate a turbulent channel flow at high Reynolds number. It is known that the mean velocity profile has a mismatch between the RANS and LES regions in hybrid simulations of a channel flow. The velocity mismatch is reproduced and its dependence on the location of the RANS/LES interface and on the type of RANS model is examined in order to better understand its properties. To remove the mismatch and to obtain better velocity profiles, additional filtering is applied to the velocity components in the wall-parallel planes near the interface. The additional filtering was previously introduced to simulate a channel flow at low Reynolds number. It is shown that the filtering is effective in reducing the mismatch even at high Reynolds number. Profiles of the velocity fluctuations of runs with and without the additional filtering are examined to help understand the reason for the mismatch. Due to the additional filtering, the wall-normal velocity fluctuation increases at the bottom of the LES region. The resulting velocity field creates the grid-scale shear stress more efficiently, and an overestimate of the velocity gradient is removed. The dependence of the velocity profile on the grid point number is also investigated. It is found that the velocity gradient in the core region is underestimated in the case of a coarse grid. Attention should be paid not only to the velocity mismatch near the interface but also to the velocity profile in the core region in hybrid simulations of a channel flow at high Reynolds number. PACS47.27.Eq; 47.27.Nz; 47.60.+i     

Polynomial Filtration Laws for Low Reynolds Number Flows Through Porous Media

In this study, we use the method of homogenization to develop a filtration law in porous media that includes the effects of inertia at finite Reynolds numbers. The result is much different than the empirically observed quadratic Forchheimer equation. First, the correction to Darcy’s law is initially cubic (not quadratic) for isotropic media. This is consistent with several other authors (Mei and Auriault, J Fluid Mech 222:647–663, 1991; Wodié and Levy, CR Acad Sci Paris t.312:157–161, 1991; Couland et al. J Fluid Mech 190:393–407, 1988; Rojas and Koplik, Phys Rev 58:4776–4782, 1988) who have solved the Navier–Stokes equations analytically and numerically. Second, the resulting filtration model is an infinite series polynomial in velocity, instead of a single corrective term to Darcy’s law. Although the model is only valid up to the local Reynolds number, at the most, of order 1, the findings are important from a fundamental perspective because it shows that the often-used quadratic Forchheimer equation is not a universal law for laminar flow, but rather an empirical one that is useful in a limited range of velocities. Moreover, as stated by Mei and Auriault (J Fluid Mech 222:647–663, 1991) and Barree and Conway (SPE Annual technical conference and exhibition, 2004), even if the quadratic model were valid at moderate Reynolds numbers in the laminar flow regime, then the permeability extrapolated on a Forchheimer plot would not be the intrinsic Darcy permeability. A major contribution of this study is that the coefficients of the polynomial law can be derived a priori, by solving sequential Stokes problems. In each case, the solution to the Stokes problem is used to calculate a coefficient in the polynomial, and the velocity field is an input of the forcing function, F, to subsequent problems. While numerical solutions must be utilized to compute each coefficient in the polynomial, these problems are much simpler and robust than solving the full Navier–Stokes equations.     

Nonlinear evolution of instabilities behind spheres and disks

We have performed precise and systematic experiments with PIV in order to measure the velocity field in the wake of a solid sphere and of a disk in a water channel, in the range of intermediate Reynolds number in which stationary and oscillatory instabilities appear, including the hairpin shedding regime. From these experimental data, we study the modal decomposition of the streamwise vorticity in an unsteady case and we describe the full nonlinear evolution of the bifurcating branches. We compare these results with recent theoretical and numerical studies on instability in the vortex shedding process at these intermediate Reynolds numbers.     

Instabilities of a gas-liquid flow between two inclined plates analyzed using the Navier–Stokes equations

The paper is devoted to a theoretical analysis of a counter-current gas-liquid flow between two inclined plates. We linearized the Navier–Stokes equations and carried out a stability analysis of the basic steady-state solution over a wide variation of the liquid Reynolds number and the gas superficial velocity. As a result, we found two modes of the unstable disturbances and computed the wavelength and phase velocity of their neutral disturbances varying the liquid and gas Reynolds number. The first mode is a “surface mode” that corresponds to the Kapitza's waves at small values of the gas superficial velocity. We found that the dependence of the neutral disturbance wavelength on the liquid Reynolds number strongly depends on the gas superficial velocity, the distance between the plates and the channel inclination angle for this mode. The second mode of the unstable disturbances corresponds to the transition to a turbulent flow in the gas phase and there is a critical value of the gas Reynolds number for this mode. We obtained that this critical Reynolds number weakly depends on both the channel inclination angle, the distance between the plates and the liquid flow parameters for the conditions considered in the paper. Despite a thorough search, we did not find the unstable modes that may correspond to the instability in frame of the viscous (or inviscid) Kelvin–Helmholtz heuristic analysis.     

Intense Reynolds-stress events in turbulent ducts

The aim of the present work is to investigate the role of intense Reynolds shear-stress events in the generation of the secondary flow in turbulent ducts. We consider the connected regions of flow where the product of the instantaneous fluctuations of two velocity components is higher than a threshold based on the long-time turbulence statistics, in the spirit of the three-dimensional quadrant analysis proposed by Lozano-Durán et al. (J. Fluid Mech., vol. 694, 2012, pp. 100–130). We examine both the geometrical properties of these structures and their contribution to the mean in-plane velocity components, and we perfom a comparison with turbulent channel flow at similar Reynolds number. The contribution to a certain mean quantity is defined as the ensemble average over the detected coherent structures, weighted with their own occupied volume fraction. In the core region of the duct, the contribution of intense events to the wall-normal component of the mean velocity is in very good agreement with that in the channel, despite the presence of the secondary flow in the former. Additionally, the shapes of the three-dimensional objects do not differ significantly in both flows. In the corner region of the duct, the proximity of the walls affects both the geometrical properties of the coherent structures and the contribution to the mean component of the vertical velocity. However, such contribution is less relevant than that of the complementary portion of the flow not included in such objects. Our results show that strong Reynolds shear-stress events are affected by the presence of a corner but, despite the important role of these structures in the dynamics of wall-bounded turbulent flows, their contribution to the secondary flow is relatively low, both in the core and in the corner.     

Linear stability of two-dimensional steady flow in wavy-walled channels

Linear stability of fully developed two-dimensional periodic steady flows in sinusoidal wavy-walled channels is investigated numerically. Two types of channels are considered: the geometry of wavy walls is identical and the location of the crest of the lower and upper walls coincides (symmetric channel) or the crest of the lower wall corresponds to the furrow of the upper wall (sinuous channel). It is found that the critical Reynolds number is substantially lower than that for plane channel flow and that when the non-dimensionalized wall variation amplitude is smaller than a critical value (about 0.26 for symmetric channel, 0.28 for sinuous channel), critical modes are three-dimensional stationary and for larger , two-dimensional oscillatory instabilities set in. Critical Reynolds numbers of sinuous channel flows are smaller for three-dimensional disturbances and larger for two-dimensional disturbances than those of symmetric channel flows. The disturbance velocity distribution obtained by the linear stability analysis suggests that the three-dimensional stationary instability is mainly caused by local concavity of basic flows near the reattachment point, while the critical two-dimensional mode resembles closely the Tollmien–Schlichting wave for plane Poiseuille flow.     

蜿蜒边界下平面流的线性流动稳定性

为便于数值分析,蜿蜒河流水动力和演变模型中一般隐性假设二次时均流-二次涡的关系与明渠流时均流-明渠湍流的关系相同,但由于高雷诺数下的DNS算力限制和实验尺度限制,这种隐含假设是否成立目前尚无相关湍流研究来支撑.文章试图通过分析明渠湍流和二次湍流发展初期的研究,侧面揭示其湍流结构的异同.通过对曲线正交坐标系下的平面二维NS方程使用双参数摄动的方法,建立了一种求解蜿蜒边界弱非线性层流的摄动解法,并推导得出一个适用于蜿蜒边界的EOS方程以及其特征值问题的解法.蜿蜒边界下弱非线性层流解为一系列蜿蜒谐波分量的叠加,其中线性部分使得两壁产生流速差,非线性部分随着雷诺数增大呈指数增长.水流的扰动增长率特征谱的第一模态与直道流相似,由3条曲线、4个波段合成,但其长波段和短波段的扰动流场与直道流不同,所有短波段的扰动流速近似于KH涡.蜿蜒边界对内部水流扰动有一定的选择性.偏角幅值越大扰动增长越快;蜿蜒波数的影响则为先增后减,有一个使扰动增长最快的蜿蜒波数.扰动流场由一个典型的TS波和一对波包形式的二次涡叠加而成,波包只有纵向流速分量,包络线由蜿蜒波数控制,波包内是与直道扰动波参数相同的TS波.     

PIV with volume lighting in a narrow cell: An efficient method to measure large velocity fields of rapidly varying flows

In this work we test a methodology for PIV measurements when a large field of view is required in planar confined geometries. Using a depth of field larger than the channel width, we intend to measure the in-plane variations of the velocity of the fluid averaged through the width of the channel, and we examine in which operating conditions this becomes possible. Measurements of the flow through a narrow channel by PIV are challenging because of the strong velocity gradients that develop between the walls. In particular, all techniques that use small particles as tracers have to deal with the possible migration of the tracers in the direction perpendicular to the walls. Among the complex mechanisms for migration, we focus on the so called Segré-Silberberg effect which can lead to transverse migration of neutrally buoyant tracers of finite size. We report experimental PIV measurements in a Hele-Shaw cell of 1 mm gap, which have been carried out by using neutrally buoyant tracers of size around 10 μm. By considering steady flows, we have observed, in particular flow regimes, the effect of an accumulation of the tracers at a certain distance to the wall due to the so called Segré-Silberberg effect. The particle migration is expected to occur at any Reynolds numbers but the migration velocity depends on the Reynolds number. A significant migration therefore takes place each time the observation duration is large enough compared to the migration time. For a given observation duration, the tracers remain uniformly distributed at low Reynolds numbers whereas they all accumulate at the equilibrium position at large ones. When using volume lighting, the PIV algorithm provides the average velocity of the flow through the gap at low Reynolds number, while it leads to the velocity of the flow at the equilibrium position of the tracers at large Reynolds numbers. By considering unsteady flows, we have observed that the migration does not occur if the timescale of flow variation is short compared to the time required for the parabolic flow to develop across the gap. In this case, there is no transverse velocity gradient and the PIV algorithm provides the fluid velocity. Altogether, these results allow us to propose guidelines for the interpretation of PIV measurements in confined flow, which are based on the theoretical predictions of the tracer migration derived by Asmolov [1].     

The onset of turbulence in physiological pulsatile flow in a straight tube

 An empirical correlation for the onset of turbulence in physiological pulsatile flow is presented. We pumped three different test fluids of kinematic viscosity 0.008–0.035 cm²/s through four straight tubes 0.4–3.0 cm in diameter. A Scotch yoke mechanism provided an oscillatory sine wave flow component of known stroke volume and frequency. We adjusted the mean flow independently until we detected signal instabilities from hot film wall shear stress probes. The critical peak Reynolds number was found to correlate with the Womersley parameter and the Strouhal number as a power law function with a root-mean-square (rms) error of 15.2%. Experimental measurements of the laminar velocity profile are compared to theoretical predictions from Poiseuille’s law and Womersley’s solution. Received: 30 October 1995/Accepted: 7 April 1997     

Flow over periodic hills: an experimental study

Two-dimensional flow over periodically arranged hills was investigated experimentally in a water channel. Two-dimensional particle image velocimetry (PIV) and one-dimensional laser Doppler anemometry (LDA) measurements were undertaken at four Reynolds numbers ( \text5,600 ￡ Re ￡ \text37,000\text{5,600} \le Re \le \text{37,000}). Two-dimensional PIV field measurements were thoroughly validated by means of point-by-point 1D LDA measurements at certain positions of the flow. A detailed study of the periodicity and the homogeneity was undertaken, which demonstrates that the flow can be regarded as two-dimensional and periodic for Re 3 \text10,000Re \ge \text{10,000}. We found a decreasing reattachment length with increasing Reynolds number. This is connected to a higher momentum in the near-wall zone close to flow separation which comes from the velocity speed up above the obstacle. This leads to a velocity overshoot directly above the hill crest which increases with Reynolds number as the inner layer depth decreases. The flow speed up above that layer is independent of the Reynolds number which supports the assumption of inviscid flow disturbance in the outer layer usually made in asymptotic theory for flow over small hills.     

Direct numerical simulation of supersonic pipe flow at moderate Reynolds number

We study compressible turbulent flow in a circular pipe at computationally high Reynolds number. Classical related issues are addressed and discussed in light of the DNS data, including validity of compressibility transformations, velocity/temperature relations, passive scalar statistics, and size of turbulent eddies. Regarding velocity statistics, we find that Huang’s transformation yields excellent universality of the scaled Reynolds stresses distributions, whereas the transformation proposed by Trettel and Larsson (2016) yields better representation of the effects of strong variation of density and viscosity occurring in the buffer layer on the mean velocity distribution. A clear logarithmic layer is recovered in terms of transformed velocity and wall distance coordinates at the higher Reynolds number under scrutiny (Re_τ ≈ 1000), whereas the core part of the flow is found to be characterized by a universal parabolic velocity profile. Based on formal similarity between the streamwise velocity and the passive scalar transport equations, we further propose an extension of the above compressibility transformations to also achieve universality of passive scalar statistics. Analysis of the velocity/temperature relationship provides evidence for quadratic dependence which is very well approximated by the thermal analogy proposed by Zhang et al. (2014). The azimuthal velocity and scalar spectra show an organization very similar to canonical incompressible flow, with a bump-shaped distribution across the flow scales, whose peak increases with the wall distance. We find that the size growth effect is well accounted for through an effective length scale accounting for the local friction velocity and for the local mean shear.     

Study of fracture permeability using Lattice Gas Automata

We study the problem of flow permeability of fracture joints using Lattice-Gas Automata simulations. We model the fracture as a rough channel bounded by a self-affine surface. Changing the surface roughness exponent, rough walls having different microstructures are obtained. Different relative roughnesses — defined as the height of the largest surface asperity divided by the mean aperture — are obtained pulling apart the two surfaces that constitute the rough walls of the channel. We calculate the macroscopic variables volume flow rate and pressure difference using microscopic balances. In the low Reynolds number regime the pressure difference and the flow rate are linearly related (the behavior is described by Darcy's law). In this regime, we study the effect of geometry on the permeability. We have found that permeability is independent of the surface roughness exponentH and it is fully determined in terms of the relative roughness and mean aperture of the fracture joint. For larger Reynolds numbers a transition to a regime in which pressure difference and flow rate are not longer linearly related is observed. This transition is observed in a domain of Reynolds numbers for which the behavior in a smooth channel remains linear. We discuss this transition.     

Super-stable parallel flows of multiple visco-plastic fluids

We consider the stability of a multi-layer plane Poiseuille flow of two Bingham fluids. It is shown that this two-fluid flow is frequently more stable than the equivalent flow of either fluid alone. This phenomenon of super-stability results only when the yield stress of the fluid next to the channel wall is larger than that of the fluid in the centre of the channel, which need not have a yield stress. Our result is in direct contrast to the stability of analogous flows of purely viscous generalised Newtonian fluids, for which short wavelength interfacial instabilities can be found at relatively low Reynolds numbers. The results imply the existence of parameter regimes where visco-plastic lubrication is possible, permitting transport of an inelastic generalised Newtonian fluid in the centre of a channel, lubricated at the walls by a visco-plastic fluid, travelling in a stable laminar flow at higher flow rates than would be possible for the single fluid alone.     

Characteristics of solitary waves on a falling liquid film sheared by a turbulent counter-current gas flow

We report an experimental investigation of a falling water film sheared by a turbulent counter-current air flow in an inclined rectangular channel. Film thickness and wave velocity measurements associated with visual observation are conducted to study the influence of the air flow on controlled traveling waves consisting of a large wave hump preceded by capillary ripples. First, we focus on the variation of the shape, amplitude and velocity of the waves as the gas velocity is gradually increased. We demonstrate that the amplitude of the main hump grows substantially even for moderate gas velocities, whereas modification of the wave celerity becomes significant above a specific gas velocity around 4 m/s, associated with an alteration of the capillary region. The influence of the gas flow on 3D secondary instabilities of the solitary waves detected in a previous study Kofman et al. (2014), namely rugged or scallop waves, is also investigated. We show that the capillary mode is damped while the inertial mode is enhanced by the interfacial shear. Next, the gas velocity is increased until the onset of upstream-moving patterns referred to as flooding in our experiments. At moderate inclination angles (typically < 7 ^○), flooding occurs for a gas velocity around 8 m/s and is initiated at the scallop wave crests by a backward wave-breaking phenomenon preceded by the onset of ripples on the flat residual film separating two waves. At high inclination angle, a rapid development of solitons is observed as the air velocity is increased preventing the waves to turn back. Finally, at high liquid Reynolds number, sudden and intermittent events are triggered consisting of very large amplitude waves that go back upwards very fast. These “slugs” either extend over the whole width of the channel or are very localized and can thus potentially evolve towards atomization.     

Effects of two-dimensional V-shaped grooves on turbulent channel flow

Detailed Laser Doppler velocimeter (LDV) measurements have been carried out in a turbulent rectangular channel flow with one rough wall. The roughness elements of two-dimensional spanwise 120° V-shaped grooves are periodically arranged with different depths and pitches. The Reynolds number based on the centerline velocity, and the channel half height ranges from 2,740 to 20,000. The comparisons of turbulence statistics over smooth and rough walls indicate that the present roughness leads to a significant change in the turbulence both in the inner and in the outer flow. Particularly, the distribution density of the grooves is a key parameter to evaluate the effect of roughness. The low-Reynolds-number dependence of turbulence statistics is also observed. The rough walls with the same pitch-to-depth ratio exhibit the equivalent roughness function under the corresponding Reynolds numbers. The disagreement of velocity defect profiles between smooth and rough walls challenges the defect universal law. The variations of the turbulence stresses and Reynolds shear stress decomposition in the outer layer suggest that the turbulent motions may be modified by the present grooves. The importance of sweep events for the present groove-roughened walls is reflected by the differences in relative contribution to Reynolds shear stress from each quadrant and the higher-order moments over smooth and rough walls.     

Unsteady PTV velocity field past an airfoil solved with DNS: Part 1. Algorithm of hybrid simulation and hybrid velocity field at <Emphasis Type="Italic">Re</Emphasis> ≈ 10<Superscript>3</Superscript>

We develop a hybrid unsteady-flow simulation technique combining direct numerical simulation (DNS) and particle tracking velocimetry (PTV) and demonstrate its capabilities by investigating flows past an airfoil. We rectify instantaneous PTV velocity fields in a least-squares sense so that they satisfy the equation of continuity, and feed them to the DNS by equating the computational time step with the frame rate of the time-resolved PTV system. As a result, we can reconstruct unsteady velocity fields that satisfy the governing equations based on experimental data, with the resolution comparable to numerical simulation. In addition, unsteady pressure distribution can be solved simultaneously. In this study, particle velocities are acquired on a laser-light sheet in a water tunnel, and unsteady flow fields are reconstructed with the hybrid algorithm solving the incompressible Navier–Stokes equations in two dimensions. By performing the hybrid simulation, we investigate nominally two-dimensional flows past the NACA0012 airfoil at low Reynolds numbers. In part 1, we introduce the algorithm of the proposed technique and discuss the characteristics of hybrid velocity fields. In particular, we focus on a vortex shedding phenomenon under a deep stall condition (α = 15°) at Reynolds numbers of Re = 1000 and 1300, and compare the hybrid velocity fields with those computed with two-dimensional DNS. In part 2, the extension to higher Reynolds numbers is considered. The accuracy of the hybrid simulation is evaluated by comparing with independent experimental results at various angles of attack and Reynolds numbers up to Re = 10⁴. The capabilities of the hybrid simulation are also compared with two-dimensional unsteady Reynolds-Averaged Navier–Stokes (URANS) solutions in part 2. In the first part of these twin papers, we demonstrate that the hybrid velocity field approaches the PTV velocity field over time. We find that intensive alternate vortex shedding past the airfoil, which is predicted by the two-dimensional DNS, is substantially suppressed in the hybrid simulation and the resultant flow field is similar to the PTV velocity field, which is projection of the three-dimensional velocity field on the streamwise plane. We attempt to identify the motion that originates three-dimensional flow patterns by highlighting the deviation of the PTV velocity field from the two-dimensional governing equations at each snapshot. The results indicate that the intensive spots of the deviation appear in the regions in which three-dimensional instabilities are induced in the shear layer separated from the pressure side.     

Upscaling Forchheimer law

We investigate the high velocity flow in heterogeneous porous media. The model is obtained by upscaling the flow at the heterogeneity scale where the Forchheimer law is assumed to be valid. We use the method of multiple scale expansions, which gives rigorously the macroscopic behaviour without any prerequisite on the form of the macroscopic equations. We show that Forchheimer law does not generally survive upscaling. The macroscopic flow law is strongly non-linear and anisotropic. A 2-point Padé approximation of the flow law in the form of a Forchheimer law is given. However, this approximation is generally poor. These results are illustrated in two particular cases: a layered composite porous media and a composite constituted by a square array of circular porous inclusions embedded in a porous matrix. We show that non-linearities are sources of anisotropy.     

Symmetries, Invariance and Scaling-Laws in Inhomogeneous Turbulent Shear Flows

An approach to derive turbulent scaling laws based on symmetry analysis is presented. It unifies a large set of scaling laws for the mean velocity of stationary parallel turbulent shear flows. The approach is derived from the Reynolds averaged Navier–Stokes equations, the fluctuation equations, and the velocity product equations, which are the dyad product of the velocity fluctuations with the equations for the velocity fluctuations. For the plane case the results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the centre of a Couette flow and in the centre of a rotating channel flow, and a new exponential mean velocity profile that is found in the mid-wake region of high Reynolds number flat-plate boundary layers. The algebraic scaling law is confirmed in both the centre and the near wall regions in both experimental and DNS data of turbulent channel flows. For a non-rotating and a moderately rotating pipe about its axis an algebraic law was found for the axial and the azimuthal velocity near the pipe-axis with both laws having equal scaling exponents. In case of a rapidly rotating pipe, a new logarithmic scaling law for the axial velocity is developed. The key elements of the entire analysis are two scaling symmetries and Galilean invariance. Combining the scaling symmetries leads to the variety of different scaling laws. Galilean invariance is crucial for all of them. It has been demonstrated that two-equation models such as the k– model are not consistent with most of the new turbulent scaling laws.