Showing papers on "Hele-Shaw flow published in 2015"

Mechanism of frequency lock-in in vortex-induced vibrations at low Reynolds numbers

Abstract: In this study, a CFD-based linear dynamics model combined with the direct Computational Fluid Dynamics/Computational Structural Dynamics (CFD/CSD) simulation method is utilized to study the physical mechanisms underlying frequency lock-in in vortex-induced vibrations (VIVs). An identification method is employed to construct the reduced-order models (ROMs) of unsteady aerodynamics for the incompressible flow past a vibrating cylinder at low Reynolds numbers ( ). Reduced-order-model-based fluid–structure interaction models for VIV are also constructed by coupling ROMs and structural motion equations. The effects of the natural frequency of the cylinder, mass ratio and structural damping coefficient on the dynamics of the coupled system at are investigated. The results show that the frequency lock-in phenomenon at low Reynolds numbers can be divided into two patterns according to different induced mechanisms. The two patterns are ‘resonance-induced lock-in’ and ‘flutter-induced lock-in’. When the natural frequency of the cylinder is in the vicinity of the eigenfrequency of the uncoupled wake mode (WM), only the WM is unstable. The dynamics of the coupled system is dominated by resonance. Meanwhile, for relatively high natural frequencies (i.e. greater than the eigenfrequency of the uncoupled WM), the structure mode becomes unstable, and the coupling between the two unstable modes eventually leads to flutter. Flutter is the root cause of frequency lock-in and the higher vibration amplitude of the cylinder than that of the resonance region. This result provides evidence for the finding of De Langre (J. Fluids Struct., vol. 22, 2006, pp. 783–791) that frequency lock-in is caused by coupled-mode flutter. The linear model exactly predicts the onset reduced velocity of frequency lock-in compared with that of direct numerical simulations. In addition, the transition frequency predicted by the linear model is in close coincidence with the amplitude of the lift coefficient of a fixed cylinder for high mass ratios. Therefore, it confirms that linear models can capture a significant part of the inherent physics of the frequency lock-in phenomenon.

Effects of cylinder Reynolds number on the turbulent horseshoe vortex system and near wake of a surface-mounted circular cylinder

Abstract: The turbulent horseshoe vortex (HV) system and the near-wake flow past a circular cylinder mounted on a flat bed in an open channel are investigated based on the results of eddy-resolving simulations and supporting flow visualizations. Of particular interest are the changes in the mean flow and turbulence statistics within the HV region as the necklace vortices wrap around the cylinder’s base and the variation of the mean flow and turbulence statistics in the near wake, in between the channel bed and the free surface. While it is well known that the drag crisis induces important changes in the flow past infinitely long circular cylinders, the changes are less understood and more complex for the case of flow past a surface-mounted cylinder. This is because even at very high cylinder Reynolds numbers, ReD, the flow regime remains subcritical in the vicinity of the bed surface due to the reduction of the incoming flow velocity within the bottom boundary layer. The paper provides a detailed discussion of the changes in the flow physics between cylinder Reynolds numbers at which the flow in the upstream part of the separated shear layers (SSLs) is laminar (ReD = 16 000, subcritical flow regime) and Reynolds numbers at which the transition occurs inside the attached boundary layers away from the bed and the flow within the SSLs is turbulent (ReD = 5 ∗ 105, supercritical flow regime). The changes between the two regimes in the dynamics and level of coherence of the large-scale coherent structures (necklace vortices, vortex tubes shed in the SSLs and roller vortices shed in the wake) and their capacity to induce high-magnitude bed friction velocities in the mean and instantaneous flow fields and to amplify the near-bed turbulence are analyzed. Being able to quantitatively and qualitatively describe these changes is critical to understand Reynolds-number-induced scale effects on sediment erosion mechanisms around cylinders mounted on a loose bed, which is a problem of great practical relevance (e.g., for pier scour studies).

Direct numerical simulation of moderate-Reynolds-number flow past arrays of rotating spheres

Abstract: Direct numerical simulations with an immersed boundary-lattice Boltzmann method are used to investigate the effects of particle rotation on flows past random arrays of mono-disperse spheres at moderate particle Reynolds numbers. This study is an extension of a previous study of the authors [Q. Zhou and L.-S. Fan, “Direct numerical simulation of low-Reynolds-number flow past arrays of rotating spheres,” J. Fluid Mech. 765, 396–423 (2015)] that explored the effects of particle rotation at low particle Reynolds numbers. The results of this study indicate that as the particle Reynolds number increases, the normalized Magnus lift force decreases rapidly when the particle Reynolds number is in the range lower than 50. For the particle Reynolds number greater than 50, the normalized Magnus lift force approaches a constant value that is invariant with solid volume fractions. The proportional dependence of the Magnus lift force on the rotational Reynolds number (based on the angular velocity and the diameter of the spheres) observed at low particle Reynolds numbers does not change in the present study, making the Magnus lift force another possible factor that can significantly affect the overall dynamics of fluid-particle flows other than the drag force. Moreover, it is found that both the normalized drag force and the normalized torque increase with the increase of the particle Reynolds number and the solid volume fraction. Finally, correlations for the drag force, the Magnus lift force, and the torque in random arrays of rotating spheres at arbitrary solids volume fractions, rotational Reynolds numbers, and particle Reynolds numbers are formulated.

Purely elastic flow instabilities in microscale cross-slot devices

Abstract: We present an experimental investigation of viscoelastic fluid flow in a cross-slot microgeometry under low Reynolds number flow conditions. By using several viscoelastic fluids, we investigate the effects of the microchannel bounding walls and the polymer solution concentration on the flow patterns. We demonstrate that for concentrated polymer solutions, the flow undergoes a bifurcation above a critical Weissenberg number (Wi) at which the flow becomes asymmetric but remains steady. The appearance of this elastic instability depends on the channel aspect ratio, defined as the ratio between the depth and the width of the channels. At high aspect ratios, when bounding wall effects are reduced, two types of elastic instabilities were observed, one in which the flow becomes asymmetric and steady, followed by a second instability at higher Wi, in which the flow becomes time-dependent. When the aspect ratio decreases, the bounding walls have a stabilizing effect, preventing the occurrence of steady asymmetric flow and postponing the transition to unsteady flow to higher Wi. For less concentrated solutions, the first elastic instability to steady asymmetric flow is absent and only the time-dependent flow instability is observed.

Time-resolved flow dynamics and Reynolds number effects at a wall-cylinder junction

Abstract: This study investigated the physics of separated turbulent flows near the vertical intersection of a flat wall with a cylindrical obstacle. The geometry imposes an adverse pressure gradient on the incoming boundary layer. As a result, flow separates from the wall and reorganizes to a system of characteristic flow patterns known as the horseshoe vortex. We studied the time-averaged and instantaneous behaviour of the turbulent horseshoe vortex using planar time-resolved particle image velocimetry (TRPIV). In particular, we focused on the effect of Reynolds number based on the diameter of the obstacle and the bulk approach velocity, . Experiments were carried out at : , and . Data analysis emphasized time-averaged and turbulence quantities, time-resolved flow dynamics and the statistics of coherent flow patterns. It is demonstrated that two large-scale vortical structures dominate the junction flow topology in a time-averaged sense. The number of additional vortices with intermittent presence does not vary substantially with . In addition, the increase of turbulence kinetic energy (TKE), momentum and vorticity content of the flow at higher is documented. The distinctive behaviour of the primary horseshoe vortex for the case is manifested by episodes of rapid advection of the vortex to the upstream, higher spatio-temporal variability of its trajectory, and violent eruptions of near-wall fluid. Differences between this experimental run and those at lower Reynolds numbers were also identified with respect to the spatial extents of the bimodal behaviour of the horseshoe vortex, which is a well-known characteristic of turbulent junction flows. Our findings suggest a modified mechanism for the aperiodic switching between the dominant flow modes. Without disregarding the limitations of this work, we argue that Reynolds number effects need to be considered in any effort to control the dynamics of junction flows characterized by the same (or reasonably similar) configurations.

Transport and buckling dynamics of an elastic fibre in a viscous cellular flow

Abstract: We study, using both experiment and theory, the coupling of transport and shape dynamics for elastomeric fibres moving through an inhomogeneous flow. The cellular flow, created electromagnetically in our experiment, comprises many identical cells of counter-rotating vortices, with a global flow geometry characterized by a backbone of stable and unstable manifolds connecting hyperbolic stagnation points. Our mathematical model is based upon slender-body theory for the Stokes equations, with the fibres modelled as inextensible elastica. Above a certain threshold of the control parameter, the elasto-viscous number, transport of fibres is mediated by their episodic buckling by compressive stagnation point flows, lending an effectively chaotic component to their dynamics. We use simulations of the model to construct phase diagrams of the fibre state (buckled or not) near stagnation points in terms of two variables that arise in characterizing the transport dynamics. We show that this reduced statistical description quantitatively captures our experimental observations. By carefully reproducing the experimental protocols and time scales of observation within our numerical simulations, we also quantitatively explain features of the measured buckling probability curve as a function of the effective flow forcing. Finally, we show within both experiment and simulation the existence of short and long time scales in the evolution of fibre conformation.

Thermodynamics analysis of hydromagnetic third grade fluid flow through a channel filled with porous medium

Abstract: In this work, analysis has been carried out to study the entropy generation rate in the flow and heat transfer of hydromagnetic third grade fluid between horizontal parallel plates saturated with porous materials. The flow is induced by a constant pressure gradient applied in the flow direction and also influenced by a uniform magnetic field that is applied across the flow channel. The equations governing the fluid flow are modeled, non-dimensionalized and solved analytically using regular perturbation method. The effect of various flow parameters on the fluid flow is presented graphically and discussed.

Transition of transient channel flow after a change in Reynolds number

Abstract: It has previously been shown that the transient flow in a channel following a step increase of Reynolds number from 2800 to 7400 (based on channel half-height and bulk velocity) is effectively a laminar-turbulent bypass transition even though the initial flow is turbulent (He & Seddighi, J. Fluid Mech., vol. 715, 2013, pp. 60-102). In this paper, it is shown that the transient flow structures exhibit strong contrasting characteristics in large and small flow perturbation scenarios. When the increase of Reynolds number is large, the flow is characterized by strong elongated streaks during the initial period, followed by the occurrence and spreading of isolated turbulent spots, as shown before. By contrast, the flow appears to evolve progressively and the turbulence regeneration process remains largely unchanged during the flow transient when the Reynolds number ratio is low, and streaks do not appear to play a significant role. Despite the major apparent differences in flow structures, the transient flow under all conditions considered is unambiguously characterized by laminar-turbulent transition, which exhibits itself clearly in various flow statistics. During the pre-transition period, the time-developing boundary layers in all the cases show a strong similarity to each other and follow closely the Stokes solution for a transient laminar boundary layer. The streamwise fluctuating velocity also shows good similarity in the various cases, irrespective of the appearance of elongated streaks or not, and the maximum energy growth exhibits a linear rate similar to that in a spatially developing boundary layer. The onset of transition is clearly definable in all cases using the minimum friction factor, and the critical time thus defined is strongly correlated with the free-stream turbulence in a power-law form.

A self-consistent model for the saturation dynamics of the vortex shedding around the mean flow in the unstable cylinder wake

Abstract: The supercritical instability leading to the Benard-von Karman vortex street in a cylinder wake is a well known example of supercritical Hopf bifurcation: the steady solution becomes linearly unstable and saturates into a periodic limit cycle. Nonetheless, a simplified physical formulation accurately predicting the transition dynamics of the saturation process is lacking. Building upon our previous work, we present here a simple self-consistent model that provides a clear description of the saturation mechanism in a quasi-steady manner by means of coupling the instantaneous mean flow with its most unstable eigenmode and its instantaneous amplitude through the Reynolds stress. The system is coupled for different oscillation amplitudes, providing an instantaneous mean flow as function of an equivalent time. A transient physical picture is described, wherein a harmonic perturbation grows and changes in amplitude, frequency, and structure due to the modification of the mean flow by the Reynolds stress forcing, saturating when the flow is marginally stable. Comparisons with direct numerical simulations show an accurate prediction of the instantaneous amplitude, frequency, and growth rate, as well as the saturated mean flow, the oscillation amplitude, frequency, and the resulting mean Reynolds stresses.

The effects of a submerged non-erodible triangular obstacle on bottom propagating gravity currents

Abstract: The flow induced by a compositional (e.g., temperature or salinity driven) gravity current propagating over a fixed non-erodible triangular bottom-mounted obstacle is investigated based on 3-D large eddy simulations. The paper discusses how the flow physics (e.g., type and characteristics of the reflected bore, dynamics of flow instabilities affecting mixing, and turbulence structure) and main flow variables (e.g., the proportion of the flow advected over the obstacle, the height of the reflected flow and speed of the reflected bore, the height of the lower layer and front speed of the current downstream of the obstacle, and the drag force) change as a function of the incoming gravity current type (lock-exchange vs constant-flux), relative obstacle height, and Reynolds number. A particular focus is on the flow structure during the two possible quasi-steady regimes that can occur in such flows. The predictive capabilities of shallow flow theory models to estimate the main flow parameters during these two regimes are investigated. An analytical model is proposed to estimate the mean streamwise drag force on the obstacle during the two regimes. Finally, the bed friction velocity distributions are used to identify regions where significant erosion will occur in the case of a loose surface.

Dynamics of isolated confined air bubbles in liquid flows through circular microchannels: an experimental and numerical study

Abstract: Experimental and numerical studies are performed to characterise the dynamics of isolated confined air bubbles in laminar fully developed liquid flows within channels of diameters d = 0.5 mm and d = 1 mm. Water and glycerol are used as the continuous liquid phase, and therefore, a large range of flow capillary numbers 10−4 < Ca < 10−1 and Reynolds numbers 10−3 < Re < 103 are covered. An extensive investigation is performed on the effect of bubble size and flow capillary number on different flow parameters, such as the shape and velocity of bubbles, thickness of the liquid film formed between the bubbles and the channel wall, and the development lengths in front and at the back of the bubbles. The micro-particle shadow velocimetry technique (μPSV) is employed in the experimental measurements allowing simultaneous quantification of important flow parameters using a single sequence of high-speed greyscale images recorded at each test condition. Bubble volume and flow rate of the continuous liquid phase are precisely determined in the post-processing stage using the μPSV images. These parameters are then used as initial and boundary conditions to set up CFD simulations reproducing the corresponding two-phase flow. Simulations based on the volume of fluid technique with the aim of capturing the interface dynamics are performed with both ANSYS Fluent v. 14.5, here augmented by implementing self-defined functions to improve the accuracy of the surface tension force estimation, and ESI OpenFOAM v. 2.1.1. The present approach not only results in valuable findings on the underlying physics involved in the problem of interest but also allows us to directly compare and validate results that are currently obtained by the experimental and computational methods. It is believed that similar methodology can be employed to rigorously investigate more complex two-phase flow regimes in micro-geometries.

Effect of Péclet number on miscible rectilinear displacement in a Hele-Shaw cell.

Abstract: The influence of fluid dispersion on the Saffman-Taylor instability in miscible fluids has been investigated in both the linear and the nonlinear regimes. The convective characteristic scales are used for the dimensionless formulation that incorporates the Peclet number (Pe) into the governing equations as a measure for the fluid dispersion. A linear stability analysis (LSA) is performed in a similarity transformation domain using the quasi-steady-state approximation. LSA results confirm that a flow with a large Pe has a higher growth rate than a flow with a small Pe. The critical Peclet number (Pec) for the onset of instability for all possible wave numbers and also a power-law relation of the onset time and most unstable wave number with Pe are observed. Unlike the radial source flow, Pec is found to vary with t0. A Fourier spectral method is used for direct numerical simulations (DNS) of the fully nonlinear system. The power-law dependence of the onset of instability ton on Pe is obtained from the DNS and found to be inversely proportional to Pe and it is in good agreement with that obtained from the LSA. The influence of the anisotropic dispersion is analyzed in both the linear and the nonlinear regimes. The results obtained confirm that for a weak transverse dispersion merging happens slowly and hence the small wave perturbations become unstable. We also observ that the onset of instability sets in early when the transverse dispersion is weak and varies with the anisotropic dispersion coefficient, e, as ∼√[e], in compliance with the LSA. The combined effect of the Korteweg stress and Pe in the linear regime is pursued. It is observed that, depending on various flow parameters, a fluid system with a larger Pe exhibits a lower instantaneous growth rate than a system with a smaller Pe, which is contrary to the results when such stresses are absent.

A microscopic particle image velocimetry method for studying the dynamics of immiscible liquid–liquid interactions in a porous micromodel

Abstract: The development of an experimental protocol to investigate the flow field produced by the interaction of two immiscible liquids flowing through a porous network is reported. The experimental protocol allows simultaneous quantification of the velocity distribution in a multi-liquid system based on the microscopic particle image velocimetry technique. The experimental challenges associated with this unique application are discussed, including two-liquid imaging and interface tracking, and solutions that couple refractive index matching and fluorescent signal separation are described. The technique was applied to both single- and two-liquid flows in a two-dimensional pore network comprising a staggered array of circular pillars wherein the flow was driven by a steady pressure gradient. Both drainage and imbibition were considered herein with a focus on fluid–fluid front migration and effects owing to the passage of the interface. The velocity distribution obtained for these two-liquid-phase flow scenarios revealed several peculiarities when compared to the reference case of single-liquid-phase flow. In particular, the instabilities associated with the interfacial processes propagate downstream and perturb the flow field, resulting in dramatic differences from the regular and periodic flow paths typical of steady-state, single-phase flow. Additionally, the passage of the interface does not restore previous flow patterns, but instead yields complex preferential flow paths that mutually interact with residual trapped pockets of fluid. Such dynamical events must be quantified in order to properly model the pore-scale physics central to fully understanding the wealth of practical applications represented by this model flow system.

Vortex dynamics in a pipe T-junction: Recirculation and sensitivity

Abstract: In the last few years, many researchers have noted that regions of recirculating flow often exhibit particularly high sensitivity to spatially localized feedback. We explore the flow through a T-shaped pipe bifurcation—a simple and ubiquitous, but generally poorly understood flow configuration—and provide a complex example of the relation between recirculation and sensitivity. When Re ≥ 320, a phenomenon resembling vortex breakdown occurs in four locations in the junction, with internal stagnation points appearing on vortex axes and causing flow reversal. The structure of the recirculation is similar to the traditional bubble-type breakdown. These recirculation regions are highly sensitive to spatially localized feedback in the linearized Navier–Stokes operator. The flow separation at the corners of the “T,” however, does not exhibit this kind of sensitivity. We focus our analysis on the Reynolds number of 560, near the first Hopf bifurcation of the flow.

Simulation of Secondary and Separated Flow in Diffusing S Ducts

Abstract: The focus of this paper is on the numerical simulation of compressible flow in diffusing S-ducts; this flow is characterized by secondary flow as well as regions of boundary-layer separation. The S-duct geometry produces streamline curvature and an adverse pressure gradient resulting in these flow characteristics. Two S-duct geometries are employed in this investigation: one was used in an experimental study conducted at NASA John H. Glenn Research Center at Lewis Field in the early 1990s, and the other is a benchmark configuration proposed by AIAA Propulsion Aerodynamics Workshop to assess the accuracy and best practices of computational fluid dynamics solvers. The computational fluid dynamics flow solver ANSYS-FLUENT is employed in the investigation of compressible turbulent flow through the S duct. A second-order accurate, steady, density-based solver is employed in a finite-volume framework. The three-dimensional Reynolds-averaged Navier–Stokes equations are solved on a structured mesh with a number o...

Pre-asymptotic Transport Upscaling in Inertial and Unsteady Flows Through Porous Media

Abstract: In most classical formulations of flow and transport through porous media Reynolds numbers are assumed to be small ( $$\mathrm{Re}<1$$ ), meaning that the role of inertia is considered negligible. However, many examples of practical relevance exist where this is not the case and inertial effects can be important leading to changes in flow structure and even giving rise to unsteady and turbulent flows as Reynolds numbers become larger. This change in flow structure can have a profound impact on how solutes are transported through the porous medium, influencing how effective large-scale transport should be modeled. Here we simulate, using high-resolution numerical models, flow and transport through an idealized porous medium for flow conditions over a range of Reynolds numbers, including steady and unsteady flows. For all these conditions we propose and test three upscaled models for transport—an advection dispersion equation, an uncorrelated spatial model (USM) and a spatial Markov model (SMM). The USM and SMM fall into the wider and more general family of continuous time random walk models. We test these models by their ability to reproduce pre-asymptotic and asymptotic plume second centered moments and breakthrough curves. We demonstrate that for steady flows where inertial effects are strong, the spatial Markov model outperforms the other two, faithfully capturing many of the non-Fickian features of transport, while for unsteady flows the uncorrelated spatial model performs best, due to the fact that unsteadiness in the flow field dampens the role of correlation on large scale transport. We conclude that correlation must be accounted for to properly upscale transport in steady flows, while it can be neglected in unsteady flows.

The gas flow diode effect: theoretical and experimental analysis of moderately rarefied gas flows through a microchannel with varying cross section

Abstract: Moderately rarefied gas flows are clearly distinguished from viscous flow in the continuum regime and from free molecular flow at high rarefaction. Being of relevance for various technical applications, the understanding of such flow processes is crucial for considerable enhancement in micro electromechanical systems (MEMS) and vacuum techniques. In this work, we focus on the isothermal rarefied gas flow through long channels with longitudinally varying cross section. We apply two approaches, an analytical one and a numerical one that is based on the solution of the linearized S-model, both allowing us to predict the mass flow rate in diverging and converging flow directions for arbitrary pressure gradients. Both approaches are validated by CO2, N2 and Ar permeation experiments on tapered microchannels manufactured by means of micromilling. The local Knudsen numbers ranged from 0.0471 to 0.2263. All the numerical and analytical results are in good agreement to the experimental data and show that the mass flow rate is significantly higher when the duct is perfused in converging direction. The understanding of the physical phenomenon of this gas flow diode effect might pave the way for novel components in MEMS such as static one-way valves.

Turbulent states in plane Couette flow with rotation

Abstract: Shearing and rotational forces in fluids can significantly alter the transport of momentum. A numerical investigation was undertaken to study the role of these forces using plane Couette flow subject to rotation about an axis perpendicular to both wall-normal and streamwise directions. Using a set of progressively higher Reynolds numbers up to Re = 5200, we find that the momentum flux, measured by the wall shear stress, for a given Re is a non-monotonic function of rotation number, Ro. For low-to-moderate Reynolds numbers, we find a maximum that is associated with flow fields that are dominated by downstream vortices and calculations of 2D vortices capture the maximum also quantitatively. For higher Reynolds numbers, a second stronger maximum emerges at smaller rotation numbers, closer to non-rotating plane Couette flow. It is carried by flows with a markedly 3D structure and cannot be captured by 2D vortex studies. As the Reynolds number increases, this maximum becomes stronger and eventually overtakes the one associated with the 2D flow state.

Laboratory Study on 3D Flow Structures Induced by Zero-Height Side Weir and Implications for 1D Modeling

Abstract: Side weir flows induce complex three-dimensional (3D) structures on the flow field in the main channel. Extensive experiments at the laboratory scale were conducted in order to investigate the flow field for the particular case of a side weir with a zero height crest in both fixed bed and steady flow conditions. Detailed measurements of the flow surface, discharge, and velocity field were performed in the main channel upstream, alongside, and downstream of the side weir location. They allowed detection and analysis of the three-dimensional flow structures. Elaborations of the measurements show the implications of such three-dimensional nature of the flow on integral quantities such as Coriolis and Boussinesq coefficients as well as the specific energy and discharge coefficient that are important in the widely used one-dimensional (1D) models. Moreover, the results may help to better understand other free-surface problems similar to the side weir flow, such as lateral diversions, river bifurcations...

Counterflow quantum turbulence of He-II in a square channel: Numerical analysis with nonuniform flows of the normal fluid

Abstract: We perform a numerical analysis of counterflow quantum turbulence of superfluid 4He with nonuni- form flows by using the vortex filament model. In recent visualization experiments nonuniform laminar flows of the normal fluid, namely, Hagen-Poiseuille flow and tail-flattened flow, have been observed. Tail-flattened flow is a novel laminar flow in which the outer part of the Hagen-Poiseuille flow becomes flat. In our simulation, the velocity field of the normal fluid is prescribed to be two nonuniform profiles. This work addresses a square channel to obtain important physics not revealed in the preceding numerical works. In the studies of the two profiles we analyze the statistics of the physical quantities. Under Hagen-Poiseuille flow, inhomogeneous quantum turbulence appears as a statistically steady state. The vortex tangle shows a characteristic space-time oscillation. Under tail-flattened flow, the nature of the quantum turbulence depends strongly on that flatness. Vortex line density increases significantly as the profile becomes flatter, being saturated above a certain flatness. The inhomogeneity is significantly reduced in comparison to the case of Hagen-Poiseuille flow. Investigating the behavior of quantized vortices reveals that tail-flattened flow is an intermedi- ate state between Hagen-Poiseuille flow and uniform flow. In both profiles we obtain a characteristic inhomogeneity in the physical quantities, which suggests that a boundary layer of superfluid appears near a solid boundary. The vortex tangle produces a velocity field opposite to the applied superfluid flow, and, consequently, the superfluid flow becomes smaller than the applied one.

Optimization of the selective frequency damping parameters using model reduction

Abstract: In the present work, an optimization methodology to compute the best control parameters, χ and Δ, for the selective frequency damping method is presented. The optimization does not suppose any a priori knowledge of the flow physics, neither of the underlying numerical methods, and is especially suited for simulations requiring large quantity of grid elements and processors. It allows for obtaining an optimal convergence rate to a steady state of the damped Navier-Stokes system. This is achieved using the Dynamic Mode Decomposition, which is a snapshot-based method, to estimate the eigenvalues associated with global unstable dynamics. Validations test cases are presented for the numerical configurations of a laminar flow past a 2D cylinder, a separated boundary-layer over a shallow bump, and a 3D turbulent stratified-Poiseuille flow.

Low frequency vibration induced streaming in a Hele-Shaw cell

Abstract: When an acoustic wave propagates in a fluid, it can generate a second order flow whose characteristic time is much longer than the period of the wave. Within a range of frequency between ten and several hundred Hz, a relatively simple and versatile way to generate streaming flow is to put a vibrating object in the fluid. The flow develops vortices in the viscous boundary layer located in the vicinity of the source of vibrations, leading in turn to an outer irrotational streaming called Rayleigh streaming. Because the flow originates from non-linear time-irreversible terms of the Navier-Stokes equation, this phenomenon can be used to generate efficient mixing at low Reynolds number, for instance in confined geometries. Here, we report on an experimental study of such streaming flow induced by a vibrating beam in a Hele-Shaw cell of 2 mm span using long exposure flow visualization and particle-image velocimetry measurements. Our study focuses especially on the effects of forcing frequency and amplitude on flow dynamics. It is shown that some features of this flow can be predicted by simple scaling arguments and that this vibration-induced streaming facilitates the generation of vortices.