Motion of a slender body near the free surface of a compressible fluid

An analytic solution to the problem of motion of a slender rigid body in a semi-infinite domain of a compressible fluid is obtained for the case when the body moves in parallel to the free surface at a constant velocity. This problem is similar to the problem of motion of a hydrofoil ship whose wing-like device allows it to lift its hull above the water surface and to decrease the friction and drag forces limiting the speed of usual ships. During its motion in water, a hydrofoil produces a lift force. The obtained analytic solution allows one to derive explicit expressions for the drag force and for the lift force in the limiting cases of relatively small and large depths. When depth is small, the drag force is greater than that in an infinite medium, since the wave drag is additionally evolved. When the velocity increases and approaches the sound velocity, the forces exerted on the body increase without limit, which is typical for a linear formulation of the problem.     

Drag measurements at Mach 5 using a stress wave force balance

The stress wave force balance, which has been used for measurements of drag on short models in hypersonic impulse facilities, is investigated here for its suitability for drag measurements on a longer, axisymmetric model. The sensitivity of the balance to loading distribution is investigated and results are reported for experiments on a 5° semi-angle cone, 425 mm in length and of 1.71 kg mass. Experimental drag measurements are shown to be in good agreement with theoretical levels. An investigation into the period over which the stress wave force balance can be used is addressed and, for the present model, the balance is shown to be suitable for measurements in flows of durations of one to several milliseconds with an estimated accuracy of ±10%.     

Interaction of a shock with a sphere suspended in a vertical shock tube

Shock wave interaction with a sphere is one of the benchmark tests in shock dynamics. However, unlike wind tunnel experiments, unsteady drag force on a sphere installed in a shock tube have not been measured quantitatively. This paper presents an experimental and numerical study of the unsteady drag force acting on a 80 mm diameter sphere which was vertically suspended in a 300 mm x 300 mm vertical shock tube and loaded with a planar shock wave of M _s = 1.22 in air. The drag force history on the sphere was measured by an accelerometer installed in it. Accelerometer output signals were subjected to deconvolution data processing, producing a drag history comparable to that obtained by solving numerically the Navier-Stokes equations. A good agreement was obtained between the measured and computed drag force histories. In order to interpret the interaction of shock wave over the sphere, high speed video recordings and double exposure holographic interferometric observations were also conducted. It was found that the maximum drag force appeared not at the time instant when the shock arrived at the equator of the sphere, but at some earlier time before the transition of the reflected shock wave from regular to Mach reflection took place. A negative value of the drag force was observed, even though for a very short duration of time, when the Mach stem of the transmitted shock wave relfected and focused at the rear stagnation point of the sphere.Received: 31 March 2003, Accepted: 7 July 2003, Published online: 2 September 2003     

Axisymmetric creeping motion of a slip spherical particle in a nonconcentric spherical cavity

A semianalytical study of the creeping flow caused by a spherical fluid or solid particle with a slip surface translating in a viscous fluid within a spherical cavity along the line connecting their centers is presented in the quasisteady limit of small Reynolds number. In order to solve the Stokes equations for the flow field, a general solution is constructed from the superposition of the fundamental solutions in the two spherical coordinate systems based on both the particle and cavity. The boundary conditions on the particle surface and cavity wall are satisfied by a collocation technique. Numerical results for the hydrodynamic drag force exerted on the particle are obtained with good convergence for various values of the ratio of particle-to-cavity radii, the relative distance between the centers of the particle and cavity, the relative viscosity or slip coefficient of the particle, and the slip coefficient of the cavity wall. In the limits of the motions of a spherical particle in a concentric cavity and near a cavity wall with a small curvature, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force exerted on the particle for all cases is a monotonic increasing function of the ratio of particle-to-cavity radii, and becomes infinite in the touching limit. For a specified ratio of particle-to-cavity radii, the drag force is minimal when the particle is situated at the cavity center and increases monotonically with its relative distance from the cavity center to infinity in the limit as it is located extremely away from the cavity center. The drag force acting on the particle, in general, increases with an increase in its relative viscosity or with a decrease in its slip coefficient for a given configuration, but surprisingly, there are exceptions when the ratio of particle-to-cavity radii is large.     

圆柱管道内运动的小球受到的黏滞阻力的数值分析

为了更好地运用落球法测量研究流体的黏滞系数,研究小球在黏性流体中下落的受力情况,本文对小球在充满黏性流体的圆柱管道的下落过程进行分析。利用COMSOL4.4仿真模拟,建立了合理的仿真模型,并分析了小球受到的黏滞阻力与小球的大小、下落位置的关系。结果表明:选择速度项二阶近似、压强项一阶近似的离散化方法,可以得到和理论值非常相符的仿真结果;当下落过程中小球球心始终在圆柱轴线上时,小球受到的黏滞阻力相对于Stokes力的修正系数,是小球半径与圆柱管道半径的比例函数,本文得到了更大范围的符合理论解的修正系数;当下落过程中小球的球心偏离圆柱轴线时,对于同样大小的小球,黏滞阻力、压强力、黏性力均随着球心到轴线的距离先减小后增大,且具有不同的极小值点。     

Drag of a Plate Planing on the Shallow Water with Formation of Waves

The plane problem of the plate planing at a constant velocity on the surface of a heavy, ideal, incompressible, finite-depth fluid is considered. The approximate, depth-independent expression for the force acting on the plate is derived from the linear distribution of the fluid velocity along the plate and the height of the flow stagnation point, without regard for jet formation near the leading edge. In this approximate formulation the plate drag depends on its velocity and the trailing edge immersion and does not depend on the planing angle. Experiments and numerical calculations in the exact formulation are performed in the near-critical flow regimes. It is shown that the wave patterns in the experiments and numerical calculations coincide, the formula for the drag being in agreement with the numerical experiments. An approximate criterion of the formation of waves going away from the plate in the forward direction is proposed.     

A numerical simulation of suspension flow using a constitutive model based on anisotropic interparticle interactions

We report a Brownian configuration field implementation of a recent constitutive equation for suspensions, reported by Phan-Thien et al. 1999. The numerical method is a hybrid technique, which combines a modification of the Brownian configuration field method described by Hulsen et al. 1997 and the adaptive viscosity split stress formulation proposed by Sun et al. 1996. The implementation is used to examine the flow past a sphere in a tube. The relative viscosity derived from the drag force/sedimentation velocity agrees well with a well-known empiricism. In addition, the ratio of the pressure force to the drag on the sphere seems to be weakly dependent on the volume fraction, and is somewhat higher than Brenner's results of 1962, which were derived for Newtonian fluids. Received: 5 April 1999/Accepted: 27 September 1999     

Constitutive relations for the interaction force in multicomponent particulate flows

In the mechanics of multiphase (or multicomponent) mixtures, one of the outstanding issues is the formulation of constitutive relations for the interaction force. In this paper, we give a brief review of the various relations proposed for this interaction force. The review is tilted toward presenting the works of those who have used the mixture theory (or the theory of interacting continua) to derive or to propose a relationship for the interaction (or diffusive) force. We propose a constitutive relation which is general and frame-indifferent and thus suitable for use in many flow conditions. At the end, we provide an alternative approach for finding the drag force on a particle in a particulate mixture. This approach has been used in the non-Newtonian fluid mechanics to find the drag force on surfaces.     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Slow motion of a slip spheroid along its axis of revolution

A combined analytical and numerical study of the Stokes flow caused by a rigid spheroidal particle translating along its axis of revolution in a viscous fluid is presented. The fluid is allowed to slip at the surface of the particle. The general solution for the stream function in prolate and oblate spheroidal coordinates can be expressed in an infinite-series form of semi-separation of variables. The slip boundary condition incorporating the shear stress at the particle surface is applied to this general solution to determine its unknown coefficients of the leading orders. The solution of these coefficients can be either numerical results obtained from a boundary-collocation method or explicit formulas derived analytically. The drag force exerted on the spheroidal particle by the fluid is evaluated with good convergence behavior for various values of the slip parameter and aspect ratio of the particle. The agreement between our hydrodynamic drag results and the relevant numerical solutions obtained previously using a singularity method is excellent. Although the drag force acting on the translating spheroid normalized by that on a corresponding sphere with equal equatorial radius increases monotonically with an increase in the axial-to-radial aspect ratio for a no-slip spheroid, it decreases monotonically as this aspect ratio increases for a perfect-slip spheroid. The normalized drag force exerted on a spheroid with a given surface slip coefficient in between the no-slip and perfect-slip limits is not a monotonic function of its aspect ratio. For a spheroid with a fixed aspect ratio, its drag force is a monotonically decreasing function of the slip coefficient of the particle.     

Comparison of different notation for equations of motion of a body in a medium flow

In [1–6], a model of a nonstationary action of a medium flow on a body moving in this flow was constructed in the form of an associated dynamical system of second order. In the literature, the representation of the aerodynamic force in integral form with a Duhamel type integral is often used (e.g., see [7, 8]). In the present paper, we pay attention to the fact that a system of ODE is equivalent not to a single integro-differential equation but to a family of such equations. Therefore, it is necessary to discuss the problem of the correspondence between their solutions. The integro-differential representation of the aerodynamic force is reduced to a form convenient to realize the procedure of separation of motions. In this case, we single out the first two approximations with respect to a small parameter. It turns out that in the case of actual airfoils one can speak of “detached” rather than “attached” mass. In the problem on the forced drag of an airfoil in a flow, it is shown that for a sufficiently large acceleration the aerodynamic force can change its direction and turn from a drag force into an “accelerating” force for some time. At the same time, in the case of free drag of a sufficiently light plate, the “acceleration” effect is not observed, but in the course of deceleration the plate moves from it original position in the direction opposite to the initial direction of motion.     

A general formula for the drag on a sphere placed in a creeping unsteady micropolar fluid flow

In the present work, we investigate the creeping unsteady motion of an infinite micropolar fluid flow past a fixed sphere. The technique of Laplace transform is used. The drag formula is obtained in the physical domain analytically by using the complex inversion formula of the Laplace transform. The well known formula of Basset for the drag on a sphere placed in an unsteady viscous fluid flow and that of Ramkissoon and Majumdar for steady motion in the case of micropolar fluids are recovered as special cases. The obtained formula is employed to calculate the drag force for some micropolar fluid flows. Numerical results are obtained and represented graphically.     

Numerical study of the turbulent channel flow under space-dependent electromagnetic force control at different Reynolds numbers

In this paper, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with the direct numerical simulation(DNS) methods for different Reynolds numbers. A formulation is derived to express the relation between the drag and the Reynolds shear stress. With the application of optimal electromagnetic force, the in-depth relations among characteristic structures in the flow field, mean Reynolds shear stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The regular quasi-streamwise vortex structures, which appear in the flow field, have the same period with that of the electromagnetic force.These structures suppress the random velocity fluctuations, which leads to the absolute value of mean Reynolds shear stress decreasing and the distribution of that moving away from the wall. Moreover, the wave number of optimal electromagnetic force increases,and the scale of the regular quasi-streamwise vortex structures decreases as the Reynolds number increases. Therefore, the rate of drag reduction decreases with the increase in the Reynolds number since the scale of the regular quasi-streamwise vortex structures decreases.     

Slow motion of a slip spherical particle perpendicular to two plane walls

The problem of the quasisteady motion of a spherical fluid or solid particle with a slip-flow surface in a viscous fluid perpendicular to two parallel plane walls at an arbitrary position between them is investigated theoretically in the limit of small Reynolds number. To solve the axisymmetric Stokes equation for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both circular cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the plane walls by the Hankel transform and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force exerted on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle and of the relative separation distances between the particle and the confining walls. For the motions of a spherical particle normal to a single plane wall and of a no-slip sphere perpendicular to two plane walls, our drag results are in good agreement with the available solutions in the literature for all relative particle-to-wall spacings. The boundary-corrected drag force acting on the particle in general increases with an increase in its relative viscosity or with a decrease in its slip coefficient for a given geometry, but there are exceptions. For a specified wall-to-wall spacing, the drag force is minimal when the particle is situated midway between the two plane walls and increases monotonically when it approaches either of the walls. The boundary effect on the particle motion normal to two plane walls is found to be significant and much stronger than that parallel to them.     

Stratified flow under an ice keel—a numerical approach

A numerical investigation of the flow of two immiscible stratified fluids under an isolated keel has been undertaken. The investigation utilized the two-dimensional Euler equations for incompressible flow, and the solution of these equations has been obtained by using the well-known finite volume marker and cell approach. Experimental drag-force measurements are also presented for a family of two-dimensional topographic models of fixed height with increasing surface slopes in a two-layer density system. The range of flow speeds explored covers the Froude number range from subcritical to fully supercritical. The drag force measurements are augmented by detailed observations of the interface distortion. The results clearly show large drag increases arising from the internal wave systems generated in the stratified flow. Very good agreement has been found between the experimental and numerical results for both the interface shape between the two fluids and the drag force on a variety of keels.     

Slow motions of a circular cylinder experiencing slip near a plane wall

A theoretical study is presented for the two-dimensional creeping flow caused by a long circular cylindrical particle translating and rotating in a viscous fluid near a large plane wall parallel to its axis. The fluid is allowed to slip at the surface of the particle. The Stokes equations for the fluid velocity field are solved in the quasi-steady limit using cylindrical bipolar coordinates. Semi-analytical solutions for the drag force and torque acting on the particle by the fluid are obtained for various values of the slip coefficient associated with the particle surface and of the relative separation distance between the particle and the wall. The results indicate that the translation and rotation of the confined cylinder are not coupled with each other. For the motion of a no-slip cylinder near a plane wall, our hydrodynamic drag force and torque results reduce to the closed-form solutions available in the literature. The boundary-corrected drag force and torque acting on the particle decrease with an increase in the slip coefficient for an otherwise specified condition. The plane wall exerts the greatest drag on the particle when its migration occurs normal to it, and the least in the case of motion parallel to it. The enhancement in the hydrodynamic drag force and torque on a translating and rotating particle caused by a nearby plane wall is much more significant for a cylinder than for a sphere.     

Numerical investigation of turbulent channel flow controlled by spatially oscillating spanwise Lorentz force

A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation (DNS) of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodically well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.     

Motion of a permeable shell in a spherical container filled with non-Newtonian fluid

This paper presents an analytical study of creeping motion of a permeable sphere in a spherical container filled with a micro-polar fluid. The drag experienced by the permeable sphere when it passes through the center of the spherical container is studied.Stream function solutions for the flow fields are obtained in terms of modified Bessel functions and Gegenbauer functions. The pressure fields, the micro-rotation components,the drag experienced by a permeable sphere, the wall correction factor, and the flow rate through the permeable surface are obtained for the frictionless impermeable spherical container and the zero shear stress at the impermeable spherical container. Variations of the drag force and the wall correction factor with respect to different fluid parameters are studied. It is observed that the drag force, the wall correction factor, and the flow rate are greater for the frictionless impermeable spherical container than the zero shear stress at the impermeable spherical container. Several cases of interest are deduced from the present analysis.     

Three-dimensional hydrodynamic modelling of tidal flows interacting with aquaculture fish cages

Three-dimensional, non-linear, non-hydrostatic simulations of rotating tidal flows interacting with aquaculture cages (represented as drag elements) at the geophysical scale are performed using an adaptive, finite volume fluid code “Gerris”. Exploiting the Gerris grid structure, sub-metre scale resolution can be obtained even for the farm scale experiments, enabling examination of the impact of the cage on the imposed tidal flows. Passive tracers are used to try to quantify these cage impacts, representing either feed or faecal matter (with specified fall speeds), or other biogeochemical markers such as dissolved oxygen. Using a relatively simple drag formulation, we show that the model is able to reproduce laboratory observations. The farm scale simulations can also be “tuned” in a similar fashion, for example by comparison with observations of total drag force on such structures, or with field measurements of flow retardation by cages. Single and multi-scale cage experiments are then examined to explore the potential impacts of perturbed horizontal and vertical flows on material redistribution through and within the cages. Even with the relatively smooth forcing and drag formulation the experiments reveal a surprising level of complexity in terms of the perturbed flows and their impact on transporting and diffusing passive material.