A sphere in a second degree polynomial creeping flow parallel to a wall

Comprehensive results are provided for the creeping flow arounda spherical particle in a viscous fluid close to a plane wall,when the external velocity is parallel to the wall and variesas a second degree polynomial in the coordinates. By linearityof Stokes equations, the solution is a sum of flows for typicalunperturbed flows: a pure shear flow, a ‘modulated shearflow’, for which the rate of shear varies linearly inthe direction normal to the wall, and a quadratic flow. Solutionsconsidered here use the bipolar coordinates technique. Theycomplement the accurate results of Chaoui and Feuillebois (2003)for the pure shear flow. The solution of Goren and O'Neill (1971)for the quadratic flow is reconsidered and a new analyticalsolution is derived for the ambient modulated shear flow. Theperturbed flow fields for these two cases are presented in detailand discussed. Results for the force and torque friction factorsare provided with a 5 x 10^–17 accuracy as a reference.For the quadratic flow, there is a force and a torque on a fixedsphere. A minimum value of the torque is found for a gap ofabout 0·18a, where a is the sphere radius. This minimumis interpreted in term of the corresponding flow structure.For the modulated shear flow, there is only a torque. The freemotion of a sphere in an ambient quadratic flow is also determined.     

The motion of a submerged sphere in a liquid under an ice sheet

The solution of the linear steady problem of the flow of an inviscid, incompressible and infinitely deep liquid around a sphere under an ice sheet, which is modelled by a thin elastic stressed plate of constant thickness is constructed. Special cases of this problem are the motion of a submerged sphere under broken ice, a membrane, and also under the free surface both in the presence and absence of capillary effects. The method of multipole expansions is used in the framework of the linear potential wave theory. The hydrodynamic loads (the wave drag and the buoyancy) acting on the body and also the distribution of the deflections of the ice sheet are calculated as a function of the body velocity, the ice thickness and the value of the compressing or stretching forces. It is shown that all the flow characteristics depend considerably on the ratio of the body velocity and the critical velocity of flexural-gravitational waves.     

Slip at the surface of a sphere translating perpendicular to a plane wall in micropolar fluid

The Stokes axisymmetrical flow caused by a sphere translating in a micropolar fluid perpendicular to a plane wall at an arbitrary position from the wall is presented using a combined analytical-numerical method. A linear slip, Basset type, boundary condition on the surface of the sphere has been used. To solve the Stokes equations for the fluid velocity field and the microrotation vector, a general solution is constructed from fundamental solutions in both cylindrical, and spherical coordinate systems. Boundary conditions are satisfied first at the plane wall by the Fourier transforms and then on the sphere surface by the collocation method. The drag acting on the sphere is evaluated with good convergence. Numerical results for the hydrodynamic drag force and wall effect with respect to the micropolarity, slip parameters and the separation distance parameter between the sphere and the wall are presented both in tabular and graphical forms. Comparisons are made between the classical fluid and micropolar fluid.      

Stokes flow for the axisymmetric motion of several spherical particles perpendicular to a plane wall

The creeping flow around several spherical particles moving on a line perpendicular to a plane wall is calculated numerically using the boundary integral method. The locations of the point forces on the surfaces of the spheres are chosen so as to describe precisely the lubrication regions when the surfaces are close to one another. Earlier results are recovered for the cases of a single sphere and a wall and of two equal spheres far from a wall. New results are presented for two (equal or unequal) spheres close to a plane wall and several equal spheres far from a wall.     

The lift on a small sphere in a linear shear flow near the interface of two immiscible fluids

Analytical expressions have been derived which predict, to lowest order, the inertial lift and the lateral migration velocity of a rigid sphere translating and rotating in a linear shear flow field near the flat interface of two immiscible fluids. This asymptotic analysis is primarily based on the assumption that the two Reynolds numbers defined by the gap width between the interface and the sphere, the shear rate and the translational slip velocity with which the spherical particle moves parallel to the interface are small. Furthermore, the radius of the sphere is assumed to be small compared to the gap width. To leading order in this creeping flow regime, the linear Stokes equations are obtained and a symmetry argument can be used to show that the Stokes solution does not predict any lift force. The transverse force experienced by the sphere and its migration velocity are due to the small but finite inertial terms in the Navier-Stokes equations, which can be studied by perturbation techniques. By applying a Green's function approach and matched asymptotic methods, which also incorporate the effects of the outer Oseen-like flow regime, the three components comprising the lift velocity have been calculated in closed form: the one induced by the shear rate only, the purely slip induced one and the one due to the interaction of the slip velocity with the shear flow field. The thus obtained expressions for the case of two immiscible fluids with arbitrary density and viscosity ratios extend the results that already exist in the literature for other flow configurations, such as an unbounded shear flow field [1] or a wall-bounded one, where the wall lies either within the leading order Stokes region [2] or in the outer Oseen region [3]. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Retracted: Unsteady flow and heat transfer on a rotating disk in a porous medium

This article has been retracted. See retraction notice DOI: 10.1002/mma.850 . An unsteady flow and heat transfer in a porous medium of a viscous incompressible fluid over a rotating disk in an otherwise ambient fluid are studied. The unsteadiness in the flow field is caused by the angular velocity of the disk which varies with time. The new self‐similar solution of the Navier–Stokes and energy equations is obtained numerically. The solution obtained here is not only the solution of the Navier–Stokes equations, but also of the boundary layer equations. Also, for a simple scaling factor, it represents the solution of the flow and heat transfer in the forward stagnation‐point region of a rotating sphere or over a rotating cone. The asymptotic behaviour of the solution for a large porosity or for a large independent variable is also examined. The surface shear stresses in the radial and tangential directions and the surface heat transfer increase as the acceleration parameter increases. Also, the surface shear stress in the radial direction and the surface heat transfer decrease with increasing porosity, but the surface shear stress in the tangential direction increases. Copyright © 2006 John Wiley & Sons, Ltd.     

Direct numerical simulation of particles in a fluid interacting with a wall

We develop a code to be applied in the context of the cleaning of wafer surfaces by hydrodynamic forces. Our goal is to study the detachment of (submicron) particles, exposed to a shear flow, from a wall by means of direct numerical simulation. The particles are treated as rigid bodies fully interacting with the fluid. To simulate moving particles in the fluid we implement an immersed boundary method with direct forcing into OpenFOAM. The particle-wall interaction is treated with a soft contact model. As first simple examples we study the elastic normal impact of a cylinder onto a wall as well as the onset of sliding of a cylinder pressed to a horizontal wall by gravity under a time-depended drag force. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stresslet resistance functions for low Reynolds number flow using deforming spheres

The stresslets of two rigid spheres in an ambient pure straining flow are obtained at low Reynolds number by defining and solving an equivalent problem of flow around deforming spheres. If the spheres are separated by a small gap, the stresslet of each sphere (the symmetric first moment of the surface stress) is a singular function of the gap width. For spheres in an ambient pure straining flow, the singularities manifest themselves as the slow convergence of numerical calculations. The methods of lubrication theory are used to calculate the singularities in the stresslets and it is shown that these new singularities can be related to singularities already found in other resistance functions. It is also shown that the singular terms can be used to improve the rate of convergence of series expressions for the stresslets. The series expressions then become valid for all separations of the spheres.     

Hydrodynamics of a particle impact on a wall

The problem of a particle impacting on a wall, a common phenomenon in particle-laden flows in the minerals and process industries, is investigated computationally using a spectral-element method with the grid adjusting to the movement of the particle towards the wall. Remeshing is required at regular intervals to avoid problems associated with mesh distortion and the constantly reducing maximum time-step associated with integration of the non-linear convective terms of the Navier–Stokes equations. Accurate interpolation between meshes is achieved using the same high-order interpolation employed by the spectral-element flow solver. This approach allows the full flow evolution to be followed from the initial approach, through impact and afterwards as the flow relaxes. The method is applied to the generic two-dimensional and three-dimensional bluff body geometries, the circular cylinder and the sphere. The principal case reported here is that of a particle colliding normally with a wall and sticking. For the circular cylinder, non-normal collisions are also considered. The impacts are studied for moderate Reynolds numbers up to approximately 1200. A cylindrical body impacting on a wall produces two vortices from its wake that convect away from the cylinder along the wall before stalling while lifting induced wall vorticity into the main flow. The situation for a sphere impact is similar, except in this case a vortex ring is formed from the wake vorticity. Again, secondary vorticity from the wall and particle plays a role. At higher Reynolds number, the secondary vorticity tends to form a semi-annular structure encircling the primary vortex core. At even higher Reynolds numbers, the secondary annular structure fragments into semi-discrete structures, which again encircle and orbit the primary core. Vorticity fields and passive tracer particles are used to characterize the interaction of the vortical structures. The evolution of the pressure and viscous drag coefficients during a collision are provided for a typical sphere impact. For a Reynolds number greater than approximately 1000 for a sphere and 400 for a cylinder, the primary vortex core produced by the impacting body undergoes a short-wavelength instability in the azimuthal/spanwise direction. Experimental visualisation using dye carried out in water is presented to validate the predictions.     

Hydrodynamic analogy for direct shear force and bending moment evaluation in a plate

A brand new interpretation of the plate bending equations is given using hydrodynamic analogy. It permits one to determine directly the shear forces and bending moments of a plate without the need of finding deflections. In engineering design of a plate it is more important to know shear forces and bending moments than the deflections. The existing numerical methods of solution of plate problems consist in determining deflections; then shear forces and bending moments are obtained by differentiating the deflection three and two times which produces great loss of accuracy. The hydrodynamic analogy method has the advantage over other numerical methods because the shear forces and bending moments are obtained directly, without the need of finding deflections and because they are obtained with better accuracy. The hydrodynamic analogy can be applied to a plate of arbitrary shape, with arbitrary boundary conditions under an arbitrary loading.     

Dynamic deformation of a stationary rubber flow

The effect of dynamic deformation on the stationary flow of a rubber composition has been experimentally investigated for comparable values of the stationary and dynamic strain rates. The dependence of the effective viscosity on the stationary shear rate is not equivalent to its dependence on the periodic shear rate amplitude. An expression is given for calculating the effective viscosity in the case of combined stationary and dynamic shear deformation. The effectiveness of the dynamic deformation, estimated in terms of the effective viscosity, depends on whether it is superimposed on the stationary flow at constant stationary shear rate or at constant stress. It is proposed to estimate the effectiveness of dynamic deformation of a stationary non-Newtonian flow in terms of the change in the power of the stationary forces. When the effective viscosity is reduced by dynamic deformation of the stationary flow, the power of the stationary forces increases at constant shear stress and falls at constant stationary shear rate.Institute of Polymer Mechanics, Academy of Sciences of the Latvian SSR, Riga. Translated from Mekhanika Polimerov, No. 3, pp. 489–496, May–June, 1971.     

Fluid flow over a thin deformable porous layer

The flow of a viscous fluid over a thin, deformable porous layer fixed to the solid wall of a channel is considered. The coupled equations for the fluid velocity and the infinitesimal deformation of the solid matrix within the porous layer are developed using binary mixture theory, Darcy's law and the assumption of linear elasticity. The case of pure shear is solved analytically for the displacement of the solid matrix, the fluid velocity both in the porous medium and the fluid above it. For a thin porous layer the boundary condition for the fluid velocity at the fluid-matrix interface is derived. This condition replaces the usual no slip condition and can be applied without solving for the flow in the porous layer.     

Latitudinally deforming rotating sphere

In this study, we consider a sphere with a surface that is fully covered by a stretchable elastic material. The radius of the sphere is fixed and it is also rotating about its radial axis. We investigate how the axisymmetric motion of a triggered fluid flow around the sphere is affected by the presence of both sphere rotation and latitudinal stretching. Considering that the deformation over the sphere commences at the pole, the problem is formulated such that the fluid flow near the pole is similar to the induced flow due to a linearly stretchable rotating disk, which has been described well in previous studies. When the rotation is omitted, the flow develops two-dimensionally under the action of pure stretching; otherwise, a three-dimensional axisymmetric fluid flow occurs, which is computed at each latitudinal angle both numerically and using a perturbation approach. The solution with wall deformation is different from the traditional character of the solution due to a solely rotating sphere. This solution is then used to compute the surface shears due to the physical drag and torque acting over the sphere. The contribution of wall stretching reduces the drag, whereas high rotation suppresses the effects of stretching to enhance the drag. More torque is required to rotate the sphere when both stretching and rotation mechanisms are in action.     

A simple and efficient implicit direct forcing immersed boundary model for simulations of complex flow

This paper presents an implementation of an implicit immersed boundary (IB) method in a flow solver based on the fractional step method and the finite volume method for complex flows involving moving boundaries and complex geometries. In this implementation, a body force caused by the immersed body is first introduced into the N-S equation to model the effect of immersed boundary. However, the body force is not pre-calculated, but implicitly determined in such a way that the velocity at the immersed boundary interpolated from the corrected velocity field accurately satisfies the no-slip and no-penetration conditions. Then, the large-eddy simulation is applied in the solver, where the subgrid-scale stress is determined by the Smagorinsky–Lilly model. Near the immersed boundaries, the subgrid-scale stress is determined by a wall model where the wall shear stress is directly calculated from the Lagrangian force(which represents the action of fluid on solid) on the immersed boundary. Such treatment makes the simulations of high Reynolds number turbulent flows feasible with the IB method. The accuracy and capability of the present method are demonstrated by simulations of a variety of both two- and three-dimensional simulations, including laminar flow past static and oscillating cylinders, rotating hydrofoil and turbulent flow around a three-dimensional circular cylinder and a sphere. It shows that the present implementation provides an easy-to-use, inexpensive and accurate technique for computational fluid dynamics in industrially relevant problems.     

On the Slow Rotation of a Sphere about a Diameter Parallel to a Nearby Plane Wall

A method using a matched asymptotic expansions technique ispresented for obtaining the Stokes flow solution for a rigidspherical body of radius a rotating uniformly about a diameterparallel to a fixed plane wall when the minimum clearance ais very much smaller than a. An inner solution is constructedwhich is valid for the region in the neighbourhood of the nearestpoints of the sphere and the wall where the flow is stronglysheared with large velocity gradients and pressure; in thisregion the leading term of the asymptotic expansion of the solutionsatisfies the equations of lubrication theory. A matching outersolution is constructed which is valid in the remainder of thefluid where the flow is weakly sheared and it is possible toassume = 0. The forces and couples acting on the sphere andthe wall are shown to be of the form (₀+₁) log +ß₀+0(,where ₀, ₁ and ß₀ are constants which have been determinedexplicitly. By use of these results it is shown that the problemwhen the sphere rolls on the wall is not well posed.     

A model of a reinforced tyre with a rod protector

A tyre design consisting of a steel-cord-reinforced rigid bond with sides connected to the wheel disc and a protector(tread) in contact with the road is examined. The tread is in the form of a set of rods connected by one end to the band, with the other end either free or in contact with the road. The rod end in contact with the road is acted upon by a force applied from the road, represented by a force normal to the road plane and a shear force due to dry friction. If the modulus of the shear force does not exceed the magnitude of the normal force multiplied by the dry friction coefficient, there is no slip at the contact point. In the opposite case, the rod end will be displaced along the road by an amount sufficient to distribute the normal and shear forces. The dynamics of longitudinal and transverse strains of the rods in contact with the road is analysed using the motion separation method in the quasi-static approximation. The behaviour of the tread rods as a function of the vertical displacement of the wheel centre is investigated, the contact area is found and the conditions are determined under which the contact area is divided into parts in which either slip of the rod ends occurs or does not occur, depending on the magnitude of the longitudinal displacement of the wheel centre or its turning relative to the horizontal axis. An analogue of a continuous model of a rod-like tread is considered, and the magnitudes of the forces and moments are found as a function of the wheel disc displacements. The equations of wheel rolling are obtained, and the conditions under which steady motions exist are found.     

固壁近旁Stokes流中粘性液滴的运动和变形

发展了一种模拟固壁近旁轴对称Stokes流中粘性液滴的运动和变形及直接计算固壁上应力的边界积分方法．用此方法对不同的液滴-固壁初始相对间距、粘度比、表面张力和浮力联合参数以及环境流动参数情况进行了数值实验．数值结果显示,由于环境流动和浮力的作用,随着时间的推进,液滴在轴向压缩,在径向拉伸．当环境流动的作用弱于浮力作用时,随着时间的推移,液滴上升并向上弯,固壁上由液滴运动所引起的应力不断减小．当环境流动的作用强于浮力作用时,随着时间的推移,液滴变得越来越扁．在这种情形,当大初始间距时,壁面上的应力随液滴的演变而增大;当小初始间距时,由环境流动、浮力及壁面对流动的较强作用的联合影响,此应力随液滴的演变而减小．由于液滴运动所引起的壁面应力的有效作用仅限于对称轴附近的一个小范围内,且此范围随液滴与固壁的初始间距增大而增大．应力的大小随初始间距增大而大为减小．表面张力对液滴变形有阻止作用．液滴粘性会减小液滴的变形和位置迁移．     

Application of Adomian's approximation to blood flow through arteries in the presence of a magnetic field

The present investigation deals with the application of Adomian's decomposition method to blood flow through a constricted artery in the presence of an external transverse magnetic field which is applied uniformly. The blood flowing through the tube is assumed to be Newtonian in character. The expressions for the two-term approximation to the solution of stream function, axial velocity component and wall shear stress are obtained in this analysis. The numerical solutions of the wall shear stress for different values of Reynold number and Hartmann number are shown graphically. The solution of this theoretical result for a particular Hartmann number is compared with the integral method solution of Morgan and Young [17].     

The nonsimilar laminar wall jet along a moving wall,in a free stream and in a free stream/moving wall

The flow of a laminar wall jet along either a moving plate, or in a free stream, or in combined moving plate and free stream is considered. The plate is isothermal and its temperature is different from that of the ambient fluid. The governing boundary-layer equations are converted into non-dimensional form and are solved numerically. Velocity and temperature profiles as well as the variation of wall shear stress and wall heat transfer are presented for all cases considered. For the case of a moving plate new results have been found although this problem has been investigated in the past. For the case of the free stream and the combined moving plate/free stream all the results are new and are presented for the first time in the literature. The three cases have been tackled with a unified way.     

A solution for the problem of stokes flow past a porous sphere

The problem of a general non-axisymmetric Stokes flow of a viscous fluid past a porous sphere is considered. The expressions for the velocity and pressure, both inside and outside the sphere are given, when the flow outside satisfies the Stokes equations and the flow inside the sphere is governed by Darcy's law. The expressions for drag and torque are given. It is found that the drag is greater or smaller than the drag in the rigid case, depending on whether the undisturbed velocity is a pure biharmonic or a harmonic respectively. The torque is same as in the rigid case.