Motion of a circular cylinder in a viscous liquid between parallel plates

The two-dimensional motion of a cylinder in a viscous fluid between two parallel walls of a vertical channel is studied. It is found that when the cylinder moves very closely along one of the channel walls, it always rotates in the direction opposite to that of contact rolling along the nearest wall. When the cylinder is away from the walls, its rotation depends on the Reynolds number of the flow. In this study two numerical methods were used. One is for the unsteady motion of a sedimenting cylinder initially released from a position close to one of the channel walls, where the Navier-Stokes equations are solved for the fluid and Newton's equations of motion are solved for the rigid cylinder. The other method is for the steady flow in which a cylinder is fixed in a uniform flow field where the channel walls are sliding past the cylinder at the speed of the approaching flow, or equivalently a cylinder is moving with a constant velocity in a quiescent fluid. The flow field, the drag, the side force (lift), and the torque experienced by the cylinder are studied in detail. The effects of the cylinder location in the channel, the size of the channel relative to the cylinder diameter, and the Reynolds number of the flow are examined. In the limit when the cylinder is translating very closely along one of the walls, the flow in the gap between the cylinder and the wall is solved analytically using lubrication theory, and the numerical solution in the other region is used to piece together the whole flow field.This research was supported by NSF DMR91-20668 through the Laboratory for Research on the Structure of Matter at the University of Pennsylvania and from the Research Foundation of the University of Pennsylvania.     

Effect of Relaxation and Retardation Time on Peristaltic Transport of the Oldroydian Viscoelastic Fluid

The influence of relaxation and retardation time on peristaltic transport of an incompressible Oldroydian viscoelastic fluid by means of an infinite train of sinusoidal waves traveling along the walls of a two-dimensional flexible channel is investigated. A perturbation solution is obtained for the case in which the amplitude ratio (wave amplitude to channel half-width) is small. The results show that the values of the mean axial velocity of an Oldroydian viscoelastic fluid is smaller than that for a Newtonian fluid. The reflux phenomena are discussed. It is found that the critical reflux pressure gradient decreases with increasing retardation time and increases with increasing relaxation time. Numerical results are reported for different values of the physical parameters of interest. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 86–95, November–December, 2005.     

Peristaltic motion of an Oldroyd‐B fluid in a channel with compliant walls

This paper is an analysis of the peristaltic flow of an Oldroyd‐B fluid in a channel with compliant walls. The flow is induced by the sinusoidal waves on the channel walls. A series solution of the resulting boundary value problem is derived under small amplitude assumption. Emphasis is placed on determining the effects of various interesting flow parameters. Copyright © 2010 John Wiley & Sons, Ltd.     

Fluid mixing induced by vibrating walls

Recent progress in micro-fluid dynamics has identified an increased demand for efficient mixing of highly viscous fluids in small channels and cavities. One way to do this is through the steady streaming generated by the vibration of solid boundaries. In this paper we investigate the mixing properties of such streaming flows in an infinite channel. A Newtonian fluid is confined within flexible walls with transverse motion in the form of standing waves of small amplitude. The velocity field is determined using a perturbation approach with the slope of the wall as a small parameter [Phys. Fluids 16 (2004) 1822]. Streaming occurs at second order with the formation of cellular flow patterns in the channel. The Lagrangian velocities were found to mimic the Eulerian except for flows at large channel half-widths and low frequencies. Most effective mixing is observed for flows at channel half-widths of similar, or lower, order than the vibratory wavelength and for sufficiently high frequencies.     

Unidirectional Heat-Gravitational Motion of Viscous Fluid in A Flat Channel with A Given Flow Rate

This paper touches upon an initial-boundary-value problem that describes the unidirectional heat-gravitational motion of fluid in a plane channel in the case of solid immobile upper and lower walls with temperature distribution thereon and in the case of a heat-insulated upper wall. The motion is caused by a joint effect of the longitudinal temperature gradient and given nonstationary flow rate. The initial-boundary-value problem is inverse relative to the pressure gradient along the channel. An exact stationary solution is obtained. A solution of the nonstationary problems in Laplace images is determined, and the results of numerical calculations are presented.     

Study of wall properties on peristaltic transport of a couple stress fluid

In this paper, we investigate the peristaltic transport of a couple stress fluid in a channel with compliant walls. Perturbation method has been used to get the solution. The flow is induced by sinusoidal traveling waves along the channel walls. The effects of wall damping, wall elastance, wall tension and couple stress parameter on the flow are investigated using the equations of fluid as well as deformable boundaries. It is found that the mean velocity at boundaries decreases with increasing couple-stress parameter and wall damping and increases with increasing wall tension and wall elastance, while the mean axial velocity increases with increasing wall tension and wall elastance and decreases with couple-stress parameter and wall damping.     

Darcy Flow in a Wavy Channel Filled with a Porous Medium

Flow in channels bounded by wavy or corrugated walls is of interest in both technological and geological contexts. This paper presents an analytical solution for the steady Darcy flow of an incompressible fluid through a homogeneous, isotropic porous medium filling a channel bounded by symmetric wavy walls. This packed channel may represent an idealized packed fracture, a situation which is of interest as a potential pathway for the leakage of carbon dioxide from a geological sequestration site. The channel walls change from parallel planes, to small amplitude sine waves, to large amplitude nonsinusoidal waves as certain parameters are increased. The direction of gravity is arbitrary. A plot of piezometric head against distance in the direction of mean flow changes from a straight line for parallel planes to a series of steeply sloping sections in the reaches of small aperture alternating with nearly constant sections in the large aperture bulges. Expressions are given for the stream function, specific discharge, piezometric head, and pressure.     

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.     

Unsteady Rarefied Gas Flow with Shock Wave in a Channel

The two-dimensional time-dependent problem of rarefied gas flow in a plane channel, formed by parallel plates of finite length and closed at one end, is solved on the basis of the kinetic S-model. The flow develops as a result of rupture of a diaphragm which separates the gas at rest in the channel and the gas at rest in a reservoir of infinite volume. The effect of gas deceleration at the channel walls under the conditions of diffuse molecular reflection from the channel walls and end face is studied. Decay of a shock wave and disappearance of a homogeneous flow zone behind the shock wave is traced for three variants of conditions at the channel inlet: (1) gas enters the channel from a reservoir of infinite length and width (as the basic variant), the simultaneous motion in the reservoir and channel being studied; (2) the high-pressure reservoir represents a usual channel section; and (3) the motion of the gas in the reservoir is not considered at all, instead of this, the boundary conditions of the evaporation-condensation type under the conditions of gas at rest in the reservoir are imposed in the inlet cross-section.     

Three-dimensionality of grooved channel flows at intermediate Reynolds numbers

 This paper describes the three-dimensional flow structure in grooved channels with different cavity lengths at intermediate Reynolds numbers. For steady flow, the three-dimensional effects are dominant near the side walls of the channel. However, after the onset of self-sustained oscillatory flow due to Tollmien–Schlichting waves as the primary instability, a secondary instability produces a three-dimensional flow with Taylor–Geortler-like vortical structure, at the bottom of the groove. This trend becomes more significant as the cavity length increases. Furthermore, the reason for three-dimensional flow is discussed using additional numerical analysis, and it is confirmed that the source of three-dimensional instability is the groove vortices due to the presence of side walls, rather than the channel traveling wave. Received: 7 September 1999/Accepted: 11 November 2000     

Random Oscillations of an Antiphase-Excited Aeroelastic System with Synchronization

We examine synchronization of the oscillatory motion of thin elastic cylindrical plates forming the walls of a channel filled by a gas. The gas motion in the channel is described by a system of Navier–Stokes equations solved by the MacCormack method of second-order accuracy. The motion of the channel walls is described by a system of dynamic, geometrically nonlinear equations of the thin-shell theory; this system is solved by the finite difference method. Kinematic and dynamic contact conditions are set at the interface between the media. By means of a numerical experiment, possible scenarios of the transition of the aeroelastic system to in-phase oscillations were identified, and the regime of random oscillations in the system with synchronization under antiphase external excitation was found.     

Particle manipulation via integration of electroosmotic flow of power-law fluids with standing surface acoustic waves (SSAW)

Standing surface acoustic wave (SSAW) based microfluidic devices have shown great promise toward fluid and particle manipulation applications in medicine, chemistry, and biotechnology. In this article, we present an analytical model for investigating continuous manipulation of particles (both synthetic and biological) within electroosmotic flow of non-Newtonian bio-fluids in a microfluidic channel under the influence of standing surface acoustic waves (SSAW). The particles are injected along the center of channel into the electroosmotically driven flow of power-law fluids, wherein their transport through the SSAW region is dictated by the hydrodynamic, electrophoretic, and acoustic forces. We first present a mathematical model to analyze the characteristics of electroosmotic flow of non-Newtonian power-law fluids in a hydrophobic slit microchannel. Next, we investigate the trajectories of particles in the flow field due to the combined effect of electroosmotic, electrophoretic, and acoustophoretic forcing mechanisms. The effect of key parameters such as particle size, their physical properties, input power, flow rate, and flow behavior index on the particle trajectories is examined while including the effect of the channel walls. The presented model delineates the methodologies of improving SSAW-based particle separation technology by considering the fluid rheology as well as the surface properties of the channel walls. Therefore, we believe that this model can serve as an efficient tool for device design and quick optimizations to explore novel applications concerning the integration of electroosmotic flows with acoustofluidic technologies.     

The PELskin project: part II—investigating the physical coupling between flexible filaments in an oscillating flow

The fluid-structure interaction mechanisms of a coating composed of flexible flaps immersed in a periodically oscillating channel flow is here studied by means of numerical simulation, employing the Euler-Bernoulli equations to account for the flexibility of the structures. A set of passively actuated flaps have previously been demonstrated to deliver favourable aerodynamic impact when attached to a bluff body undergoing periodic vortex shedding. As such, the present configuration is identified to provide a useful test-bed to better understand this mechanism, thought to be linked to experimentally observed travelling waves. Having previously validated and elucidated the flow mechanism in Paper 1 of this series, we hereby undertake a more detailed analysis of spectra obtained for different natural frequency of structures and different configurations, in order to better characterize the mechanisms involved in the organized motion of the structures. Herein, this wave-like behaviour, observed at the tips of flexible structures via interaction with the fluid flow, is characterized by examining the time history of the filaments motion and the corresponding effects on the fluid flow, in terms of dynamics and frequency of the fluid velocity. Results indicate that the wave motion behaviour is associated with the formation of vortices in the gaps between the flaps, which itself are a function of the structural resistance to the cross flow. In addition, formation of vortices upstream of the leading and downstream of the trailing flap is seen, which interact with the formation of the shear-layer on top of the row. This leads to a phase shift in the wave-type motion along the row that resembles the observation in the cylinder case.     

Plotnikov——Toland板模型中水弹性孤立波的迎撞

通过奇异摄动方法研究了在薄冰层覆盖的不可压缩理想流体表面上传播的两个水弹性孤立波之间的迎面碰撞.借助特殊的 Cosserat 超弹性壳 理论以及Kirchhoff--Love 板理论,冰层由 Plotnikov--Toland板模型描述.流体运动采用浅水假设和Boussinesq 近似. 应用Poincaré--Lighthill--Kuo 方法进行坐标变形,进而渐近求解控制方程及边界条件, 给出了三阶解的显式表达. 可以观察到碰撞后的孤立波不会改变它们的形状和振幅. 波浪轮廓在碰撞之前是对称的, 而在碰撞之后变成不对称的并且在波传播方向上向后倾斜. 弹性板和流体表面张力减小了波幅. 图示比 较了本文与已有结果可知线性板模型可作为本文的一个特例.      

The PELskin project—part I: fluid–structure interaction for a row of flexible flaps: a reference study in oscillating channel flow

Previous studies of flexible flaps attached to the aft part of a cylinder have demonstrated a favourable effect on the drag and lift force fluctuation. This observation is thought to be linked to the excitation of travelling waves along the flaps and as a consequence of that, periodic shedding of the von Kármán vortices is altered in phase. A more general case of such interaction is studied herein for a limited row of flaps in an oscillating flow; representative of the cylinder case since the transversal flow in the wake-region shows oscillating character. This reference case is chosen to qualify recently developed numerical methods for the simulation of fluid–structure interaction in the context of the EU funded ‘PELskin’ project. The simulation of the two-way coupled dynamics of the flexible elements is achieved via a structure model for the flap motion, which was implemented and coupled to two different fluid solvers via the immersed boundary method. The results show the waving behaviour observed at the tips of the flexible elements in interaction with the fluid flow and the formation of vortices in the gaps between the flaps. In addition, formation of vortices upstream of the leading and downstream of the trailing flap is seen, which interact with the formation of the shear-layer on top of the row. This leads to a phase shift in the wave-type motion along the row that resembles the observation in the cylinder case.     

Instability of a planar liquid sheet with surrounding fluids between two parallel walls

We investigate the linear and nonlinear instability of a planar liquid sheet with surrounding fluids between two parallel plane solid walls. Linear analysis shows that the maximum temporal growth rate and unstable wave number region of disturbances increase for the dilational and sinuous modes when the gap between the sheet and the wall decreases. The walls have more influence on the instability when the density ratio of the surrounding fluid to the sheet and/or the Weber number decrease. On the other hand, nonlinear analysis is performed by means of the discrete vortex method, where double vortex rows and their mirror images are placed so as to satisfy the boundary condition on the walls. Numerical results show that the walls enhance nonlinearity, which causes deformation and distortion of the sheet, whereas the nonlinearity diminishes linear growth rates except for long dilational disturbances. In particular, as the walls are placed more closely to the sheet, local sheet thinning becomes more pronounced in the long dilational mode, while the dilational mode is more strongly induced from the sinuous mode through monotonic or periodic energy exchanges between the two modes.     

Baroclinic vortex sheet production by shocks and expansion waves

Vortex sheet production by shocks and expansion waves refracting at a density discontinuity was examined and compared using an analytical solution and numerical simulations. The analytical solution showed that with a small exception, vortex sheet strength is generally stronger in fast/slow shock refractions. In contrast, expansion waves generated a stronger vortex sheet in slow/fast refractions. This difference results in larger vorticity deposited by shocks in fast/slow refractions and by expansion waves in slow/fast refractions. Shock refractions become irregular and the analytical solution fails when either incident, transmitted or reflected shock, exceeded the angle limit for an attached shock. To investigate vortex sheet production outside the range of analytical solutions and to verify the applicability of the planar-interface analytical solution to a curved interface, shock refraction through a sinusoidal interface was numerically simulated in the shock frame of reference. It is found that variation in the local incidence angle along the curved interface creates pressure waves that affect the level of deposited vorticity. This contributes to the difference between predictions from local analysis and numerical computation. Furthermore, an interesting behavior of the shock and expansion wave-deposited vorticity in supersonic ramp flow was discovered. When the high- and low-density streams were swapped, while keeping the incident flow Mach numbers constant, a vortex sheet of equal magnitude but of opposite sign was generated.     

Wave action on a movable plate immersed into a fluid. Physical and numerical experiments

The effect of the motion of a rigid inclined plate immersed into a fluid counter incident waves is studied both experimentally and numerically. In the wave water channel experiments the velocity of a trolley with a plate running freely on rails is determined as a function of the wave parameters, the immersion depth, the angle of inclination, and the plate dimensions. The interaction between the traveling waves and the plate having a single translational degree of freedom along the horizontal axis is numerically calculated in the time-dependent, two-dimensional formulation. The dependence of the upwave motion effect on the parameters varied in the full-scale experiment is analyzed. In the numerical experiment a regime of the downwave plate motion at a constant high velocity is found to exist. The channel bottom effect is estimated and the behavior of the plate with a flap is studied.     

Forced capillary-gravity wave motion of two-layer fluid in three-dimensions

A class of problems associated with forced capillary-gravity wave motion in a channel are analyzed in the presence of surface and interfacial tensions in a two-layer fluid in both the cases of finite and infinite water depths. The two and three-dimensional Green functions associated with the capillary-gravity wave problems in the presence of surface and interfacial tensions are derived using the fundamental source potentials. Using the two-dimensional Green function along with Green’s second identity, the expansion formulae for the velocity potentials associated with the capillary-gravity wavemaker problems in two-dimensions are obtained. The two-dimensional results are generalized to derive the expansion formulae for the velocity potentials associated with the forced capillary-gravity wave motion in the presence of surface and interfacial tensions in three-dimensions. Certain characteristics of the eigen-system associated with the expansion formulae are derived. The velocity potentials associated with the free oscillation of capillary-gravity waves in a closed basin and semi-infinite open channel in the presence of surface and interfacial tensions are obtained. The utility of the forced motion in a channel is demonstrated by analyzing the capillary-gravity wave reflection by a wall in a channel in the presence of surface and interfacial tensions. Long wave equations associated with capillary-gravity wave motion in the presence of surface and interfacial tensions are derived under shallow water approximation and the associated dispersion relation are obtained. Various expansion formulae and Green functions derived in the present study will be useful for analyzing a large class of physical problems in ocean engineering and mathematical physics.