On the interplay between inertial and viscoelastic effects for the flow in weakly modulated channels

The flow inside a spatially modulated channel is examined for viscoelastic fluids of the Oldroyd‐B type. The lower wall is flat and the upper wall is sinusoidally modulated. The modulation amplitude is assumed to be small. Thus, a regular perturbation expansion of the flow field coupled to a variable‐step finite‐difference scheme is used to solve the problem. Convergence and accuracy assessment against earlier experimental results indicate that there is a significant range of validity of the perturbation approach. The influences of wall geometry, inertia and viscoelasticity on the flow kinematics and stresses are investigated systematically. In particular, the interplay between the flow and fluid parameters effects on the conditions for the onset of backflow, number of vortices, their size and location is revealed. The distance between the flow separation and reattachment locations identifies the vortex size. Non‐monotonic dependence of the vortex size on elasticity is reported. The critical conditions for the onset of negative elasticity effects on vortex size are identified. The critical Reynolds number for the onset of backflow initially decreases then levels off or even increases as elasticity increases. For highly elastic fluid and large enough Reynolds number, more than one vortex appear near the lower wall. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Hydrodynamic force and torque models for a particle moving near a wall at finite particle Reynolds numbers

This research work is aimed at proposing models for the hydrodynamic force and torque experienced by a spherical particle moving near a solid wall in a viscous fluid at finite particle Reynolds numbers. Conventional lubrication theory was developed based on the theory of Stokes flow around the particle at vanishing particle Reynolds number. In order to account for the effects of finite particle Reynolds number on the models for hydrodynamic force and torque near a wall, we use four types of simple motions at different particle Reynolds numbers. Using the lattice Boltzmann method and considering the moving boundary conditions, we fully resolve the flow field near the particle and obtain the models for hydrodynamic force and torque as functions of particle Reynolds number and the dimensionless gap between the particle and the wall. The resolution is up to 50 grids per particle diameter. After comparing numerical results of the coefficients with conventional results based on Stokes flow, we propose new models for hydrodynamic force and torque at different particle Reynolds numbers. It is shown that the particle Reynolds number has a significant impact on the models for hydrodynamic force and torque. Furthermore, the models are validated against general motions of a particle and available modeling results from literature. The proposed models could be used as sub-grid scale models where the flows between particle and wall can not be fully resolved, or be used in Lagrangian simulations of particle-laden flows when particles are close to a wall instead of the currently used models for an isolated particle.     

A numerical study of flow through wavy-walled channels

A numerical procedure is developed for the analysis of flow in a channel whose walls describe a travelling wave motion. Following a perturbation method, 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 calculated using the pseudospectral collocation method. Although limited by the linear analysis, the present approach is not restricted by the Reynolds number of the flow and the wave number and frequency of the wavy-walled channel. Using the computed wall shear stresses, the positions of flow separation and reattachment are determined. The variations in velocity and pressure with frequency of excitation are also presented. © 1998 John Wiley & Sons, Ltd.     

饱和蒸发和热毛细力作用下低雷诺数液膜在波纹竖壁上的流动

考虑表面蒸发压力和热毛细力作用情况下,对饱和蒸发状态下低雷诺数自由降落液膜在小波幅正弦型波纹壁面上的流动进行理论分析。对控制微分方程及边界条件进行量纲一化并引入流函数,对微分方程及边界条件进行摄动展开,得到了这种情况下液膜流动的简化分析模型,求出了近似解析解。讨论了壁面波纹、表面张力、蒸发压力、热毛细力对液膜流动的影响。研究表明:液膜的波动幅度随蒸发强度和热毛细力的增大而增大;液膜波动与壁面波纹的相位差随蒸发强度增大而增大,随热毛细力增大而减小。     

Numerical simulation of flow of micropolar fluids in a channel with a porous wall

Two‐dimensional steady, laminar, and incompressible flow of a micropolar fluid in a channel with no‐slip at one wall and constant uniform injection through the other wall is considered for different values of the Reynolds number R. The main flow stream is superimposed by constant injection velocity at the porous wall. The micropolar model introduced by Eringen is used to describe the working fluid. An extension of Berman's similarity transformations is used to reduce governing equations to a set of nonlinear coupled ordinary differential equations (ODEs) in dimensionless form. An algorithm based on finite difference method is employed to solve these ODEs and Richardson's extrapolation is used to obtain higher order accuracy. It has been found that the magnitude of shear stress increases strictly at the impermeable wall whereas it decreases steadily at the permeable wall, by increasing the injection velocity. The maximum value of streamwise velocity and that of the microrotation both increase with increasing the magnitude of R. Copyright © 2010 John Wiley & Sons, Ltd.     

Shear‐thinning flow in weakly modulated channels

The steady flow inside a spatially modulated channel is examined for shear‐thinning and shear‐thickening fluids. The flow is induced by the translation of the lower plate. The modulation amplitude is assumed to be small. A regular perturbation expansion of the flow field is used, coupled to a variable‐step finite‐difference scheme, to solve the problem. Convergence and accuracy assessment against finite‐volume calculations indicates that there is a significant range of validity of the perturbation approach. The influence of the wall geometry, inertia and non‐Newtonian effects are investigated systematically. In particular, the influence of the flow and fluid parameters is examined on the conditions for the onset of separation, vortex size and location. Copyright © 2005 John Wiley & Sons, Ltd.     

Axisymmetric stagnation flow obliquely impinging on a moving circular cylinder with uniform transpiration

Laminar stagnation flow, axisymmetrically yet obliquely impinging on a moving circular cylinder, is formulated as an exact solution of the Navier–Stokes equations. Axial velocity is time‐dependent, whereas the surface transpiration is uniform and steady. The impinging free stream is steady with a strain rate k?. The governing parameters are the stagnation‐flow Reynolds number Re=k?a²/2ν, and the dimensionless transpiration S=U₀/k?a. An exact solution is obtained by reducing the Navier–Stokes equations to a system of differential equations governed by Reynolds number and the dimensionless wall transpiration rate, S. The system of Boundary Value Problems is then solved by the shooting method and by deploying a finite difference scheme as a semi‐similar solution. The results are presented for velocity similarity functions, axial shear stress and stream functions for a variety of cases. Shear stresses in all cases increase with the increase in Reynolds number and suction rate. The effect of different parameters on the deflection of viscous stagnation circle is also determined. Copyright © 2010 John Wiley & Sons, Ltd.     

Die swell measurements of second‐order fluids: numerical experiments

An analysis of the flow of a second‐order fluid is presented. Reference values for some variables are defined, and with these a non‐dimensional formulation of the governing equations. From this formulation, three dimensionless numbers appear; one is the Reynolds number, and two numbers that are called the first‐ and second‐dimensionless normal stress (NSD) coefficients. The equations of motion are solved by a finite element method using a commercially available program (Fidap), and the steady state converged solution was used to measure the die swell. The factors that influence die swell and that are studied in this work include: the die geometry for circular cross sectional dies, including tubular, converging, diverging, half‐converging/half‐tubular shapes; fluid characteristics such as Reynolds number and first‐ and second‐DNS coefficients (both positive and negative values); and flow rates, as determined by the maximum velocity in a parabolic velocity profile at the entrance to the die. The results suggest that shear and deformation histories of the fluid directly influence not only swell characteristics, but also convergence characteristics of the numerical simulation. © 1999 John Wiley & Sons, Ltd.     

Three-dimensional numerical simulation of sedimenting drops inside a vertical channel

Numerical simulation of sedimenting deformable drops inside a vertical channel has been performed at finite Reynolds numbers. The channel is confined by two vertical walls in x-direction and is periodic in y- and z-directions. Results are obtained using a finite difference/front-tracking method. The main dimensionless parameters are the Reynolds number, the Bond number and ratio of the length of channel to the diameter of drops. The effect of these parameters on lateral migration of drops is investigated. It is found that the wall repulsion is the main mechanism of the lateral migration of the drop, and drop migrates toward the channel axis. When the Reynolds number is relatively low, two different lateral migration regimes are observed: migration with monotonic approach and migration with damped oscillations. These regimes are affected by the dimensionless parameters. When the Bond number increases, the oscillations of drop around the centerline of channel are stronger and drop reaches the channel centerline in a larger period. Results of lateral migration of one drop are consistent with perturbation theory, and two-dimensional numerical simulations performed by Feng et al. (1994). The drag coefficient has also been calculated, and effect of various parameters has been discussed. Two drops interaction is similar to that observed by Feng et al. (1994) for two-dimensional circular cylinders. Results are consistent with experiments performed by Wu and Manasseh (1998). Simulations of four sedimenting drops show that depending on the relative size of drops, they either fall in two rows or they form a single horizontal layer and settle with a unique velocity.     

Regimes of vortex shedding from an in-line oscillating circular cylinder in the uniform flow

A finite volume method for the time dependent viscous incompressible flow around an in-line oscillating circular cylinder at Reynolds number of 200, 855 is presented in this paper. The Navier-Stokes equations in a finite volume form are solved with a moving grid system, based on a time dependent coordinate transformation. To investigate the vortex-shedding characteristics behind the circular cylinder and the effects of Reynolds number and other non-dimensional parameters such as reduced amplitude and reduced frequency， several numerical schemes have been tested with different amplitude and frequency close to Sto and a harmonic at each Reynolds number. Present numerical results indicate several types of vortex shedding mode which is known mainly depending on the reduced frequency and also the reduced amplitude, which is called synchronization or lock-on.     

Approximate solution of non-Newtonian viscoelastic fluid flow on a turbine disc for cooling purposes by using Adomian decomposition method

In this work, a powerful semi-analytical method, called Adomian decomposition method (ADM) is introduced to obtain the exact solution of heat transfer equation of a non-Newtonian fluid flow in an axisymmetric channel with a porous wall which is useful for turbine cooling applications. This method is employed to obtain the expressions for velocity and temperature fields. Then, the influence of the two dimensionless numbers: the Prandtl number (Pr) and the Reynolds number (Re) for the dynamic forces have been considered. Comparisons are made between the numerical method (NM) solution, Homotopy perturbation method (HPM), Variation iteration method (VIM) and the results of Adomian decomposition method (ADM). The results reveal that this method is very effective and simple and can be applied for other nonlinear problems.     

On the influence of numerical schemes and subgrid–stress models on large eddy simulation of turbulent flow past a square cylinder

Influence of finite difference schemes and subgrid‐stress models on the large eddy simulation calculation of turbulent flow around a bluff body of square cylinder at a laboratory Reynolds number, has been examined. It is found that the type and the order of accuracy of finite‐difference schemes and the subgrid‐stress model for satisfactory results are dependent on each other, and the grid resolution and the Reynolds number. Using computational grids manageable by workstation‐level computers, with which the near‐wall region of the separating boundary layer cannot be resolved, central‐difference schemes of realistic orders of accuracy, either fully conservative or non‐conservative, suffer stability problems. The upwind‐biased schemes of third order and the Smagorinsky eddy‐viscosity subgrid model can give reasonable results resolving much of the energy‐containing turbulent eddies in the boundary layers and in the wake and representing the subgrid stresses in most parts of the flow. Noticeable improvements can be obtained by either using higher order difference schemes, increasing the grid resolution and/or by implementing a dynamic subgrid stress model, but each at a cost of increased computational time. For further improvements, the very small‐scale eddies near the upstream corners and in the laminar sublayers need to be resolved but would require a substantially larger number of grid points that are out of the range of easily accessible computers. Copyright © 2002 John Wiley & Sons, Ltd.     

Analysis of a meshless solver for high Reynolds number flow

In this paper, the extension of an upwind least‐square based meshless solver to high Reynolds number flow is explored, and the properties of the meshless solver are analyzed both theoretically and numerically. Existing works have verified the meshless solver mostly with inviscid flows and low Reynolds number flows, and in this work, we are interested in the behavior of the meshless solver for high Reynolds number flow, especially in the near‐wall region. With both theoretical and numerical analysis, the effects of two parameters on the meshless solver are identified. The first one is the misalignment effect caused by the significantly skewed supporting points, and it is found that the meshless solver still yields accurate prediction. It is a very interesting property and is opposite to the median‐dual control volume based vertex‐centered finite volume method, which is known to give degraded result with stretched triangular/tetrahedral cells in the near‐wall region. The second parameter is the curvature, and according to theoretical analysis, it is found in the region with both large aspect ratio and curvature, and the streamwise residual is less affected; however, the wall‐normal counterpart suffers from accuracy degradation. In this paper, an improved method that uses a meshless solver for the streamwise residual and finite difference for wall‐normal residual is developed. This method is proved to be less sensitive to the curvature and provides improved accuracy. This work presents an understanding of the meshless solver for high Reynolds number flow computation, and the analysis in this paper is verified with a series of numerical experiments. Copyright © 2012 John Wiley & Sons, Ltd.     

Vortex shedding from a square cylinder in presence of a moving wall

A numerical study on the flow past a square cylinder placed parallel to a wall, which is moving at the speed of the far field has been made. Flow has been investigated in the laminar Reynolds number (based on the cylinder length) range. We have studied the flow field for different values of the cylinder to wall separation length. The governing unsteady Navier–Stokes equations are discretized through the finite volume method on a staggered grid system. A SIMPLE type of algorithm has been used to compute the discretized equations iteratively. A shear layer of negative vortex generates along the surface of the wall, which influences the vortex shedding behind the cylinder. The flow‐field is distinct from the flow in presence of a stationary wall. An alternate vortex shedding occurs for all values of gap height in the unsteady regime of the flow. The strong positive vortex pushes the negative vortex upwards in the wake. The gap flow in the undersurface of the cylinder is strong and the velocity profile overshoots. The cylinder experiences a downward force for certain values of the Reynolds number and gap height. The drag and lift are higher at lower values of the Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical study of effects of pulsatile amplitude for transitional turbulent pulsatile flow in pipes with ring‐type constrictions

The effects of pulsatile amplitude on sinusoidal transitional turbulent flows through a rigid pipe in the vicinity of a sharp‐edged mechanical ring‐type constriction have been studied numerically. Pulsatile flows were studied for transitional turbulent flow with Reynolds number (Re) of the order of 10⁴, Womersley number (Nw) of the order of 50 with a corresponding Strouhal number (St) of the order of 0.04. The pulsatile flow considered is a sinusoidal flow with dimensionless amplitudes varying from 0.0 to 1.0. Transitional laminar and turbulent flow characteristics in an alternative manner within the pulsatile flow fields were observed and studied numerically. The flow characteristics were studied through the pulsatile contours of streamlines, vorticity, shear stress and isobars. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity, turbulent shear stress are smaller during the acceleration phase when compared with those during deceleration period. Various parametric equations within a pulsatile cycle have also been formulated through numerical experimentations with different pulsatile amplitudes. In the vicinity of constrictions, the empirical relationships were obtained for the instantaneous flow rate (Q), the pressure gradient (dp/dz), the pressure loss (P_loss), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max). Elliptic relation was observed between flow rate and pressure gradient. Quadratic relations were observed between flow rate and the pressure loss, the maximum values of shear stress, wall shear stress, turbulent kinematic energy and the turbulent viscosity. Linear relationships exist between the instantaneous flow rate and the maximum values of vorticity, wall vorticity and velocity. The time‐average axial pressure gradient and the time average pressure loss across the constriction were observed to increase linearly with the pulsatile amplitude. Copyright © 1999 John Wiley & Sons, Ltd.     

Side‐wall effects on the bifurcation of the flow through a sudden expansion

Three‐dimensional computations have been performed to study the flow through a symmetric sudden expansion with an expansion ratio of 3 at low Reynolds numbers. The aspect ratio of the flow channel is allowed to vary within a wide range to examine its influence on the flow which bifurcates from a symmetric state to an asymmetric state. The results reveal that the critical Reynolds number of the symmetry‐breaking bifurcation increases while the aspect ratio is reduced. The flow behaviour near the side walls is illustrated by using limiting streamlines. The origin of the singular points identifiable on the side wall can be traced back to the recirculating flows and the relevant reattachment/separation points in the core of the channel. It is seen that the determination of the exact critical Reynolds number is not trivial because it depends on how to define asymmetric flow. Computations have also been conducted to show that a slight asymmetry in the channel geometry causes a smooth transition from symmetric to non‐symmetric states. Copyright © 2007 John Wiley & Sons, Ltd.     

On the critical problem of F. D. pressure treatment for laminar flows confined by permeable walls

In the present work the viscous (low Reynolds) flow in plane ducts confined by permeable walls has been studied. A simple model of the filtrating walls has been used, with the normal velocity component proportional to the pressure jump across the wall, resulting in a non-standard boundary value Navier-Stokes problem. A critical analysis of the appropriate boundary condition and pressure problem has led to the conclusions of employing a simple explicit finite volume approach, and of avoiding the use of higher order finite difference schemes. In this paper a special emphasis on the structure of the involved computational matrices has been given to illustrate the chosen algorithm. The latter yields a steady state solution that is second order accurate in space, and it has an accuracy in time of order ≤ Δt (the time step), due to the explicit treatment of the velocity boundary conditions along the membrane. The model has been tested to study the effects of the inlet/outlet conditions, Reynolds number and filtrating wall constant.     

被均匀流缓慢调制的有限水深毛细重力波

非线性的存在会产生高次谐波,这些谐波又反作用于原来的低次谐波,使波幅发生缓慢变化,从而产生缓慢调制现象．这里从考虑均匀流作用下的毛细重力水波基本方程出发,在不可压缩、无旋、无黏条件假设下,使用多重尺度分析方法推导出了在均匀流影响下有限深水毛细重力波振幅所满足的非线性Schr?dinger方程(NLSE)．分析了NLSE解的调制不稳定性．给出了毛细重力波调制不稳定的条件和钟型孤立波产生的条件．分析了无量纲最大不稳定增长率随无量纲水深和表面张力的变化趋势．同时给出了无量纲不稳定增长率随无量纲微扰动波数变化的曲线,呈现出了先增大后减小的趋势．最后指出均匀顺流减小了无量纲不稳定增长率及最大增长率,逆流增大了它们．由表面张力作用产生的毛细波及重力与表面张力共同作用产生的毛细重力波,与流的相互作用可以改变海表粗糙度和海洋上层流场结构,进而影响海气界面动量、热量及水汽的交换．了解海表这些短波动力机制,对卫星遥感的精确测量、海气相互作用的研究及海气耦合模式的改进等有重要意义．      

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.