Numerical study of particulate suspension flow through wavy‐walled channels

The particulate suspension flow in a channel whose walls describe a travelling wave motion is examined numerically. A perturbation method is employed and the primitive variables are expanded in a series with the wall amplitude as the perturbation parameter. The boundary conditions are applied at the mean surface of the channel and the first‐order perturbation quantities are numerically determined by solving the governing system of ordinary differential equations by shooting technique. The present approach does not impose any restriction on the Reynolds number of the flow and the wave number and frequency of the wavy‐walled channel, although it is limited by the linear analysis. The wall shear stress and the positions of flow separation and reattachment points are computed and the influence of the volume fraction density of the particles is examined. The variations of velocity and pressure of the particulate suspension flow with frequency of excitation are also presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytic Approximate Solution for a Flow of a Second-Grade Viscoelastic Fluid in a Converging Porous Channel

The problem of a two-dimensional steady flow of a second-grade fluid in a converging porous channel is considered. It is assumed that the fluid is injected into the channel through one wall and sucked from the channel through the other wall at the same velocity, which is inversely proportional to the distance along the wall from the channel origin. The equations governing the flow are reduced to ordinary differential equations. The boundary-value problem described by the latter equations is solved by the homotopy perturbation method. The effects of the Reynolds and crossflow Reynolds number on the flow characteristics are examined.     

On the validity of the perturbation approach for the flow inside weakly modulated channels

The equations governing the flow of a viscous fluid in a two‐dimensional channel with weakly modulated walls have been solved using a perturbation approach, coupled to a variable‐step finite‐difference scheme. The solution is assumed to be a superposition of a mean and perturbed field. The perturbation results were compared to similar results from a classical finite‐volume approach to quantify the error. The influence of the wall geometry and flow Reynolds number have extensively been investigated. It was found that an explicit relation exists between the critical Reynolds number, at which the wall flow separates, and the dimensionless amplitude and wavelength of the wall modulation. Comparison of the flow shows that the perturbation method requires much less computational effort without sacrificing accuracy. The differences in predicted flow is kept well around the order of the square of the dimensionless amplitude, the order to which the regular perturbation expansion of the flow variables is carried out. Copyright © 2002 John Wiley & Sons, Ltd.     

Magnetohydrodynamic mixed convection in a vertical channel

The problem of combined free and forced convective magnetohydrodynamic flow in a vertical channel is analysed by taking into account the effect of viscous and ohmic dissipations. The channel walls are maintained at equal or at different constant temperatures. The velocity field and the temperature field are obtained analytically by perturbation series method and numerically by finite difference technique. The results are presented for various values of the Brinkman number and the ratio of Grashof number to the Reynolds number for both equal and different wall temperatures. Nusselt number at the walls is determined. It is found that the viscous dissipation enhances the flow reversal in the case of downward flow while it counters the flow in the case of upward flow. It is also found that the analytical and numerical solutions agree very well for small values of ε.     

NUMERICAL SIMULATION OF STEADY FLOW IN A COMPLIANT TUBE OR CHANNEL WITH TAPERED WALL THICKNESS

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a compliant channel and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-sectional area, indicative of wave-speed flow limitation or “choking” (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, the flow rate decreases. Cross-sectional narrowing is significant but localized. For the particular wall and fluid properties used in these simulations, the area throat is located near the downstream end when the ratio of downstream-to-upstream wall thickness is 2; as wall taper is increased to 3, the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.     

Propagation of a Tollmien-Schlichting wave over the junction between rigid and compliant surfaces

The propagation of an instability wave over the junction region between rigid and compliant panels is studied theoretically. The problem is investigated using three different methods with reference to flow in a plane channel containing sections with elastic walls. Within the framework of the first approach, using the solution of the problem of flow receptivity to local wall vibration, the problem considered is reduced to the solution of an integro-differential equation for the complex wall oscillation amplitude. It is shown that at the junction of rigid and elastic channel walls the instability-wave amplitude changes stepwise. For calculating the step value, another, analytical, method of investigating the perturbation propagation process, based on representing the solution as a superposition of modes of the locally homogeneous problem, is proposed. This method is also applied to calculating the flow stability characteristics in channels containing one or more elastic sections or consisting of periodically alternating rigid and compliant sections. The third method represents the unknown solution as the sum of a local forced solution and a superposition of orthogonal modes of flow in a channel with rigid walls. The latter method can be used for calculating the eigenvalues and eigenfunctions of the stability problem for flow in a channel with uniformly compliant walls.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 31–48. Original Russian Text Copyright © 2004 by Manuilovich.     

Buoyancy-Opposed Darcy’s Flow in a Vertical Circular Duct with Uniform Wall Heat Flux: A Stability Analysis

An analytical and numerical study is presented to show that buoyancy-opposed mixed convection in a vertical porous duct with circular cross-section is unstable. The duct wall is assumed to be impermeable and subject to a uniform heat flux. A stationary and parallel Darcy’s flow with a non-uniform radial velocity profile is taken as a basic state. Stability to small-amplitude perturbations is investigated by adopting the method of normal modes. It is proved that buoyancy-opposed mixed convection is linearly unstable, for every value of the Darcy–Rayleigh number, associated with the wall heat flux, and for every mass flow rate parametrised by the Péclet number. Axially invariant perturbation modes and general three-dimensional modes are investigated. The stability analysis of the former modes is carried out analytically, while general three-dimensional modes are studied numerically. An asymptotic analytical solution is found, suitable for three-dimensional modes with sufficiently small wave number and/or Péclet number. The general conclusion is that the onset of instability selects the axially invariant modes. Among them, the radially invariant and azimuthally invariant mode turns out to be the most unstable for all possible buoyancy-opposed flows.     

Experiments on the linear instability of flow in a wavy channel

The effects of wall corrugation on the stability of wall-bounded shear flows have been examined experimentally in plane channel flows. One of the channel walls has been modified by introduction of the wavy wall model with the amplitude of 4% of the channel half height and the wave number of 1.02. The experiment is focused on the two-dimensional travelling wave instability and the results are compared with the theory [J.M. Floryan, Two-dimensional instability of flow in a rough channel, Phys. Fluids 17 (2005) 044101 (also: Rept. ESFD-1/2003, Dept. of Mechanical and Materials Engineering, The University of Western Ontario, London, Ontario, Canada, 2003)]. It is shown that the flow is destabilized by the wall corrugation at subcritical Reynolds numbers below 5772, as predicted by the theory. For the present corrugation geometry, the critical Reynolds number is decreased down to about 4000. The spatial growth rates, the disturbance wave numbers and the distribution of disturbance amplitude measured over such wavy wall also agree well with the theoretical results.     

Mixed Convection in a Vertical Porous Channel

A numerical study of mixed convection in a vertical channel filled with a porous medium including the effect of inertial forces is studied by taking into account the effect of viscous and Darcy dissipations. The flow is modeled using the Brinkman–Forchheimer-extended Darcy equations. The two boundaries are considered as isothermal–isothermal, isoflux–isothermal and isothermal–isoflux for the left and right walls of the channel and kept either at equal or at different temperatures. The governing equations are solved numerically by finite difference method with Southwell–Over–Relaxation technique for extended Darcy model and analytically using perturbation series method for Darcian model. The velocity and temperature fields are obtained for various porous parameter, inertia effect, product of Brinkman number and Grashof number and the ratio of Grashof number and Reynolds number for equal and different wall temperatures. Nusselt number at the walls is also determined for three types of thermal boundary conditions. The viscous dissipation enhances the flow reversal in the case of downward flow while it counters the flow in the case of upward flow. The Darcy and inertial drag terms suppress the flow. It is found that analytical and numerical solutions agree very well for the Darcian model. An erratum to this article is available at .     

Vibrations and stability of a periodically supported rectangular plate immersed in axial flow

Vibrations and stability of a thin rectangular plate, infinitely long and wide, periodically supported in both directions (so that it is composed by an infinite number of supported rectangular plates with slope continuity at the edges) and immersed in axial liquid flow on its upper side is studied theoretically. The flow is bounded by a rigid wall and the model is based on potential flow theory. The Galerkin method is applied to determine the expression of the flow perturbation potential. Then the Rayleigh–Ritz method is used to discretize the system. The stability of the coupled system is analyzed by solving the eigenvalue problem as a function of the flow velocity; divergence instability is detected. The convergence analysis is presented to determine the accuracy of the computed eigenfrequencies and stability limits. Finally, the effects of the plate aspect ratio and of the channel height ratio on the critical velocity giving divergence instability and vibration frequencies are investigated.     

Peristaltic flow of a Johnson-Segalman fluid through a deformable tube

To understand theoretically the flow properties of physiological fluids we have considered as a model the peristaltic motion of a Johnson–Segalman fluid in a tube with a sinusoidal wave traveling down its wall. The perturbation solution for the stream function is obtained for large wavelength and small Weissenberg number. The expressions for the axial velocity, pressure gradient, and pressure rise per wavelength are also constructed. The general solution of the governing nonlinear partial differential equation is given using a transformation method. The numerical solution is also obtained and is compared with the perturbation solution. Numerical results are demonstrated for various values of the physical parameters of interest.      

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

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

Numerical calculation of two-dimensional flow in the channel with a partially compliant wall

The present paper is concerned with the flow in a two-dimensional channel whose wall is partially compliant. The flow field is calculated by the finite-difference method. Results are as follows: (1) When the upstream condition is given by steady flow (Reynolds number Re = 50), a compliant part of the wall oscillates with a frequency nearly equal to the characteristic frequency of the elastic wall. Absolute values of the pressure drop across the compliant part become small compared with those of the plane Poiseuille flow with wholly rigid walls. This ensures under physiological conditions that the blood can be transported more easily toward distal parts due to the compliance of vessel walls. (2) When the upstream condition is given by a pulsatile flow (Womersley number α = 8), interaction arises between characteristic frequency of the wall and basic frequency of the main stream near the compliant wall. As the basic frequency of pulsatile flow decreases, absolute values of mean pressure, which drop across the compliant wall, also become small compared with those of pulsatile flow between wholly rigid walls.     

The effect of buoyancy forces on laminar convection with viscous dissipation in a vertical channel

The interaction of the viscous dissipation effect with the presence of buoyancy forces is investigated for laminar-flow heat transfer in a parallel-plate vertical channel. One of the channel walls is considered as isothermal with a prescribed temperature, while the other wall is considered as insulated. The velocity field is assumed to be parallel. The velocity field, the temperature field and the Nusselt number are obtained by a perturbation series method which employs the ratio between the Grashof number and the Reynolds number as the perturbation parameter. The radius of convergence of the perturbation series is estimated. Received on 10 December 1997     

Rayleigh–Benard convection subject to time dependent wall temperature/gravity in a fluid-saturated anisotropic porous medium

 The effect of time-periodic temperature/gravity modulation at the onset of convection in a Boussinesq fluid-saturated anisotropic porous medium is investigated by making a linear stability analysis. Brinkman flow model with effective viscosity larger than the viscosity of the fluid is considered to give a more general theoretical result. The perturbation method is applied for computing the critical Rayleigh and wave numbers for small amplitude temperature/gravity modulation. The shift in the critical Rayleigh number is calculated as a function of frequency of the modulation, viscosity ratio, anisotropy parameter and porous parameter. We have shown that it is possible to advance or delay the onset of convection by time-periodic modulation of the wall temperature and to advance convection by gravity modulation. It is also shown that the small anisotropy parameter has a strong influence on the stability of the system. The effect of viscosity ratio, anisotropy parameter, the porous parameter and the Prandtl number is discussed. Received on 28 July 2000 / Published online: 29 November 2001     

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.     

Effect of Thermal Modulation on the Onset of Convection in a Porous Medium Layer Saturated by a Nanofluid

The effect of time-periodic temperature modulation at the onset of convection in a Boussinesq porous medium saturated by a nanofluid is studied analytically. The model used for the nanofluid incorporates the effects of Brownian motion. Three types of boundary temperature modulations are considered namely, symmetric, asymmetric, and only the lower wall temperature is modulated while the upper wall is held at constant temperature. The perturbation method is applied for computing the critical Rayleigh and wave numbers for small amplitude temperature modulation. The shift in the critical Rayleigh number is calculated as a function of frequency of modulation, concentration Rayleigh number, porosity, Lewis number, and thermal capacity ratio. It has been shown that it is possible to advance or delay the onset of convection by time-periodic modulation of the wall temperature. The nanofluid is found to have more stabilizing effect when compared to regular fluid. Low frequency is destabilizing, while high frequency is always stabilizing for symmetric modulation. Asymmetric modulation and only lower wall temperature modulation is stabilizing for all frequencies when concentration Rayleigh number is greater than one.     

An interpretation of the stability of a turbulent shear flow near a rough wall

A conventional, small perturbation, stability analysis has been applied to a fully developed turbulent shear flow in a parallel duct with rough walls. This is an attempt to detect the inherent state of flow stability to quasi-regular disturbances emanating from the surface roughness elements. The surface roughness is represented by the usual roughness Reynolds number; it is fed into the analysis through an assumed mean velocity profile valid between the viscous sublayer and the inner (wall) region. An eddy viscosity model is used to secure the equation closure and the final equation for the perturbation amplitude has been solved numerically using the techniques developed for the Orr-Sommerfeld equation.Within the domain of realistic flow conditions, and for a range of surface roughness amplitudes, a local minimum of stability in terms of the longitudinal wave number has been found. However, it is not implied that it is a minimum minomorum, as only a limited range of surface roughnesses has been tried.     

Laminar flow through a channel with uniformly porous walls of different permeability

Summary The steady laminar flow of a viscous incompressible fluid through a two-dimensional channel, having fluid sucked or injected with different velocities through its uniformly porous parallel walls is considered. A solution for small suction Reynolds number has been given by the authors in a previous paper. The purpose of this paper is to present a solution valid for large Reynolds numbers for the cases of (i) suction at both walls, and (ii) suction at one wall and injection at the other. A technique of matching outer and inner expansions is used to obtain an asymptotic solution for both of these cases. Further a perturbation solution for the case of suction at one wall and injection at the other is obtained by choosing the difference between two wall velocities as the perturbation parameter. Both asymptotic and perturbation solutions are confirmed by exact numerical solutions. As expected, the resulting solutions show the presence of the usual suction boundary layers in both types of flow considered in this paper.     

Strength characteristics of the self-sustained wave in grooved channels with different groove length

The self-sustained oscillations arising in a series of grooved channels are investigated experimentally. Pressure drop, time-averaged and time-various local pressure in the grooved channels with six kinds of groove length are measured with the differential transducer and the pressure sensor, respectively, and the flow structures are visualized using the aluminum dust method. The local pressure signal shows that the self-sustained wave appears in the first or second frequency, and the Strouhal number, based on the nature frequency of the self-sustained wave, is almost equivalent for the first or second frequency in the same channel. Meanwhile, the Strouhal number for each channel decreases monotonously with the groove length. Furthermore, it is found that increasing pressure will result in higher amplitude of the self-sustained wave, this behavior is significant for the efficient heat transfer in practical engineering.