Acoustic radiation force on a rigid cylinder in an off-axis Gaussian beam near an impedance boundary

The exact equations of the axial and transverse acoustic radiation force functions of a Gaussian beam arbitrarily incident on an infinite rigid cylinder close to an impedance boundary and immersed in an ideal fluid are deduced by expressing the incident wave, the scattering wave and the boundary reflected wave in terms of the cylindrical wave function. The effects of the beam waist, the sound reflection coefficient, the cylinder position and the distance from the impedance boundary on the acoustic radiation force are studied using numerical simulations. The simulation results show that the amplitude of the acoustic radiation force function increases with beam width. Moreover, the values of the acoustic radiation force in both the axial and transverse directions reach those of a plane wave when the beam width is considerably larger than the wavelength of the Gaussian beam. The properties of the impedance boundary and the position of the cylinder in the Gaussian beam have a considerable effect on the magnitude and direction of the force. The simulation results, particularly in the case of a transverse force, indicate the presence of a negative acoustic radiation force that is related to the nondimensional frequency and position of the cylinder in the Gaussian beam.     

Acoustic radiation force on an elastic cylinder in a Gaussian beam near an impedance boundary

Acoustic radiation force (ARF) is studied by considering an infinite elastic cylinder near an impedance boundary when the cylinder is illuminated by a Gaussian beam. The surrounding fluid is an ideal fluid. Using the method of images and the translation-addition theorem for the cylindrical Bessel function, the resulting sound field including the incident wave, its reflection from the boundary, the scattered wave from the elastic cylinder, and its image are expressed in terms of the cylindrical wave function. Then, we deduce the exact equations of the axial and transverse ARFs. The solutions depend on the cylinder position, cylinder material, beam waist, reflection coefficient, distance from the impedance boundary, and absorption in the cylinder. To analyze the effects of the various factors intuitively, we simulate the radiation force for non-absorbing elastic cylinders made of stainless steel, gold, and beryllium as well as for an absorbing elastic cylinder made of polyethylene, which is a well-known biomedical polymer. The results show that the impedance boundary, cylinder material, absorption in the cylinder, and cylinder position in the Gaussian beam significantly affect the magnitude and direction of the force. Both stable and unstable equilibrium regions are found. Moreover, a larger beam waist broadens the beam domain, corresponding to non-zero axial and transverse ARFs. More importantly, negative ARFs are produced depending on the choice of the various factors. These results are particularly important for designing acoustic manipulation devices operating with Gaussian beams.     

Single Bessel tractor-beam tweezers

The tractor behavior of a zero-order Bessel acoustic beam acting on a fluid sphere, and emanating from a finite circular aperture (as opposed to waves of infinite extent) is demonstrated theoretically. Conditions for an attractive force acting in opposite direction of the radiating waves, determined by the choice of the beam’s half-cone angle, the size of the radiator, and its distance from a fluid sphere, are established and discussed. Numerical predictions for the radiation force function, which is the radiation force per unit energy density and cross-sectional surface, are provided using a partial-wave expansion method stemming from the acoustic scattering. The results suggest a simple and reliable analysis for the design of Bessel beam acoustical tweezers and tractor beam devices.     

Influence of an acoustic wave on a system of two spheres in a viscous fluid

In this article we formulate and solve the problem of the influence of radiation forces (forces created by the radiation pressure) on two spheres in a viscous fluid during the transmission of an acoustic wave. On the basis of these forces we investigate the nature of the interaction between the spheres as determined by the mutual disturbance of the flow fields around them as a result of interference between the primary and secondary waves reflected from the spheres. A previously proposed [2] approach is used in the investigations. The radiation force acting on one of the spheres is filtered by averaging the convolution of the stress tensor in the fluid with the unit normal to the surface of the sphere over a time interval and over the surface of the sphere. The stresses in the fluid are represented, to within second-order quantities in the parameters of the wave field, in terms of the velocity potentials obtained from the solution of the linear problem of the diffraction of the primary wave by the free spheres. The diffraction problem is formulated and solved within the framework of the theory of linear viscoelastic solids [6]. The case of an ideal fluid has been studied previously [3–5, 7]. Radiation forces are one of the causes of the relative drift of solid particles situated in a fluid in an acoustic field.S. P. Timoshenko Institute of Mechanics, Academy of Sciences of Ukraine, Kiev. Translated from Prikladnaya Mekhanika, Vol. 30, No. 2, pp. 33–40, February, 1994.     

Axial acoustic radiation force on a rigid cylinder near an impedance boundary for on-axis Gaussian beam

This work presents a theoretical model to calculate the acoustic radiation force on a rigid cylindrical particle immersed in an ideal fluid near a boundary for an on-axis Gaussian beam. An exact solution of the axial acoustic radiation force function is derived for a cylindrical particle by applying the translation addition theorem of cylindrical Bessel function. We analyzed the effects of the impedance boundary on acoustic radiation force of a rigid cylinder immersed in water near an impedance boundary with particular emphasis on the radius of the rigid cylinder and the distance from the cylinder center to impedance boundary. Simulation results reveal that the existence of particle trapping behavior depends on the choice of nondimensional frequency as well as the offset distance from the impedance boundary. The value of the radiation force function varies when the cylinder lies at the different position of the on-axis Gaussian beam. For the particle with different radius, the acoustic radiation force functions vary significantly with frequency. This study provides a theoretical basis for acoustic manipulation, which may benefit to the improvement and development of the acoustic control technology.     

New expressions for the radiation force function of spherical targets in stationary and quasi-stationary waves

A new expression for the radiation force function – which is the radiation force per unit energy density and unit cross-sectional surface area – for spheres in a stationary (or standing) and quasi-stationary wave is obtained based on the far-field acoustic scattering field. The radiation force function formulation has been simplified mathematically and improved into a more general form. Numerical results are presented for rigid and elastic spheres, air bubbles in water as well as liquid drops in air to illustrate the theory. It is demonstrated that expressions for the radiation force functions obtained from the far-field derivation approach are equivalent to those obtained from the near-field-based derivation.     

Acoustic radiation force of a solid elastic sphere immersed in a cylindrical cavity filled with ideal fluid

An expression for the acoustic radiation force function on a solid elastic spherical particle placed in an infinite rigid cylindrical cavity filled with an ideal fluid is deduced when the incident wave is a plane progressive wave propagated along the cylindrical axis. The acoustic radiation force of the spherical particle with different materials was computed to validate the theory. The simulation results demonstrate that the acoustic radiation force changes demonstrably because of the influence of the reflective acoustic wave from the cylindrical cavity. The sharp resonance peaks, which result from the resonance of the fluid-filled cylindrical cavity, appear at the same positions in the acoustic radiation force curve for the spherical particle with different radii and materials. Relative radius, which is the ratio of the sphere radius and the cylindrical cavity radius, has more influence on acoustic radiation force. Moreover, the negative radiation forces, which are opposite to the progressive directions of the plane wave, are observed at certain frequencies.     

Dynamic stability of pre-twisted beams with non-constant spin rates under axial random forces

This paper investigates the dynamic stability of a pre-twisted cantilever beam spinning along its longitudinal axis with a periodically varying speed and acted upon by an axial random force at the free end. The spin rate of the beam is characterized as a small periodic perturbation superimposed on a constant speed, and the axial force is assumed as the sum of a static force and a weakly stationary random process with a zero mean. Both the periodically varying spin rate and the axial random force may lead to parametric instability of the beam. In this work, the finite element method is applied first to get rid of the dependence on the spatial coordinate. The method of stochastic averaging is then adopted to obtain Ito’s equations for the system response under different resonant frequency combinations. Finally, the first-moment and the second-moment stability conditions of the beam are derived explicitly. Numerical results are presented for a simple harmonic speed perturbation and a Gaussian white noise axial force.     

Development and validation of a reduced order history force model

The history force model accounts for temporal development in fluid gradients in the viscous region surrounding a particle in point particle methods. The calculation of the history force typically requires storing and using relative velocity information during the life time of the particle. For a large number of particles integrated over large times, history force calculation can become prohibitively expensive. The current work presents a new modeling approach to calculate the history force in which a decay function is applied to a stored cumulative value of the history force. The proposed formulation is equivalent to applying the same function obtained from a constant acceleration assumption to a running average of the acceleration within the memory time of the particle. The new force model is validated with experimental measurements of settling spheres at Reynolds numbers ranging from around one to a few hundreds and at density ratios from 1.2 to about 9.32. More validation work was carried-out with experimental measurements of oscillating spheres at different frequencies and amplitudes, as well as bouncing spheres at different Reynolds numbers and density ratios. The model shows very good agreement with the experiments of settling spheres and reasonable/good agreement with oscillating and bouncing sphere experiments. The proposed model significantly reduces the computational resources required to calculate the history force especially when large number of particles need to be integrated over long times.     

Drag forces of interacting spheres in power-law fluids

Drag forces of interacting particles suspended in power-law fluid flows were investigated in this study. The drag forces of interacting spheres were directly measured by using a micro-force measuring system. The tested particles include a pair of interacting spheres in tandem and individual spheres in a cubic matrix of multi-sphere in flows with the particle Reynolds number from 0.7 to 23. Aqueous carboxymethycellulose (CMC) solutions and glycerin solutions were used as the fluid media in which the interacting spheres were suspended. The range of power-law index varied from 0.6 to 1.0. In conjunction to the drag force measurements, the flow patterns and velocity fields of power-law flows over a pair of interacting spheres were also obtained from the laser assisted flow visualization and numerical simulation.

Both experimental and computational results suggest that, while the drag force of an isolated sphere depends on the power-index, the drag coefficient ratio of an interacting sphere is independent from the power-law index but strongly depends on the separation distance and the particle Reynolds number. Our study also shows that the drag force of a particle in an assemblage is strongly positions dependent, with a maximum difference up to 38%.     

Optical measurements of pore geometry and fluid velocity in a bed of irregularly packed spheres

Imaging methods are proposed for the characterisation of liquid flows through transparent porous media of matched refractive index. The methods are based on the analysis of laser-illuminated slices, and specialized for the case in which the porous medium is composed of irregularly packed spheres. They include algorithms for the reconstruction of the three-dimensional (3D) sphere arrangement based on a laser scan of the packed bed, particle tracking velocimetry applied to the motions of micro-tracers in a laser-illuminated plane, and techniques for the co-registration of geometry and velocity measurements acquired from different slices. The methods are applied to a cylindrical flow cell filled with mono-sized spheres and operated at Reynolds number Re = 28. The data produced include the full 3D geometry of the packed spheres assembly, the 2D fluid velocity field in the axial centre-plane of the flow cell, and the corresponding porosity and velocity distributions.     

Viscous fluid flow induced by rotational-oscillatory motion of a porous sphere

Viscous fluid flow induced by rotational-oscillatorymotion of a porous sphere submerged in the fluid is determined. The Darcy formula for the viscous medium drag is supplementedwith a term that allows for the medium motion. The medium motion is also included in the boundary conditions. Exact analytical solutions are obtained for the time-dependent Brinkman equation in the region inside the sphere and for the Navier–Stokes equations outside the body. The existence of internal transverse waves in the fluid is shown; in these waves the velocity is perpendicular to the wave propagation direction. The waves are standing inside the sphere and traveling outside of it. The particular cases of low and high oscillation frequencies are considered.     

CFD-DEM simulation of the hole cleaning process in a deviated well drilling: The effects of particle shape

We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics (CFD) with the discrete element method (DEM). This numerical method allows us to incorporate the fluid–particle interactions (drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force) using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle−particle, particle−wall, and particle−drill pipe are taken into account with the Hertz–Mindlin model. We compare the transport of spheres with non-spherical particles (non-smooth sphere, disc, and cubic) constructed via the multi-sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid–particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid–solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.     

Interaction between two spheres placed in tandem arrangement in steady and pulsating flow

The interaction among two spheres in tandem formation are studied for a Reynolds number of 300 using both steady and pulsating inflow conditions. The purpose is to further investigate the force characteristics as well as the shedding patterns of the two spheres as the separation distance is changed from 1.5 to 12 sphere diameters. The method used for the simulations is the volume of solid (VOS) method, an approach based on the volume of fluid (VOF) method. Comparisons with other computational methods have shown VOS to accurately resolve the flow field around solid spheres. The results show that the separation distance plays a significant role in changing the flow patterns and shedding frequencies at moderate separation distances, whereas effect on drag is observed even at a separation distance of 12 diameters.     

Drift of spheres in a rotating fluid

The drift of spheres in a rotating fluid is investigated. The problem is studied experimentally and numerically using the Galerkin method. It is shown that for small angular velocities of the fluid Ω the drift velocity of the spheres is almost independent of Ω, but once a certain threshold value Ω_* is attained the drift velocity rapidly decreases. The experimental dependence of the translational velocity of the sphere on the fluid angular velocity is explained on the basis of a theoretical analysis.     

The hydrodynamic force and torque on a bounded sphere in Poiseuille flow

The lattice Boltzmann (LB) method is used to study the hydrodynamic force and torque acting on a sphere held stationary between parallel plates in pressure‐driven flow. This and associated flow configurations are explored in this paper. LB results are in excellent agreement with existing theory and numerical results for simple pressure‐driven flow between parallel plates, for flow through a periodic medium of spheres [Zick AA, Homsy GM. Stokes flow through periodic arrays of spheres. Journal of Fluid Mechanics 1982; 115: 13], and for the force and torque acting on a sphere held fixed at the quarter vertical position in a pressure‐driven flow between parallel plates. In the latter case, LB calculations reveal a screening effect caused by neighboring periodic images of the test sphere. It is shown that the test sphere is hydrodynamically decoupled from its periodic images when separated by approximately 30 sphere radii. LB results for force and torque as a function of sphere height and flow cell height are also reported. Copyright © 2001 John Wiley & Sons, Ltd.     

Efficient field computation using Gaussian beams for both transmission and reception

An exact representation is presented for the field inside a sphere (the observation sphere) due to primary sources enclosed by a second sphere (the source sphere). The regions bounded by the two spheres have no common points. The field of the primary sources is expressed in terms of Gaussian beams whose branch-cut disks are all centered at the origin of the source sphere. The expansion coefficients for the standing spherical waves in the observation sphere are expressed in terms of the output of Gaussian-beam receivers, whose branch-cut disks are all centered at the origin of the observation sphere. In this configuration the patterns of the transmitting and receiving beams “multiply” to produce a higher directivity than is usually seen with Gaussian beams. The areas on the unit sphere, which must be covered by the transmitting and receiving disk normals to achieve a given accuracy, diminish as 1/(ka) for ka → ∞, where a is the disk radius and k is the wavenumber. This 1/(ka) behavior leads to a single-level method with O(N^3/2) complexity for computing matrix-vector multiplications in scattering calculations (N is the number of unknowns).     

Acoustic radiation from a submerged hollow FGM sphere

An exact study based on the linear theory of elasticity is presented for the steady-state sound radiation characteristics of an arbitrarily thick radially inhomogeneous elastic isotropic hollow sphere, immersed in and filled with ideal compressible fluids, and subjected to an arbitrary axisymmetric time-harmonic driving force at its internal surface. A modal state equation with variable coefficients is set up in terms of appropriate displacement and stress functions and their spherical harmonics by means of the laminated approximation approach. Taylor’s expansion theorem is subsequently employed to solve the modal state equation, ultimately calculating a global transfer matrix. Numerical results are presented for a water-submerged/air-filled steel/zirconia FGM hollow sphere under an axisymmetric distributed internal pressure force. The effects of shell wall thickness, the material compositional gradient, frequency, and subtended polar angle of the internal pressure force on the far-field radiated pressure directivity patterns as well as the total radiated power are examined. It is demonstrated that the material gradient can significantly change the acoustical characteristics of hollow inhomogeneous sphere, especially for thick shells at high excitation frequencies. Limiting cases are considered and good agreements with available results as well as with the computations made by using a finite element package are obtained.     

Non-linear dynamic response of a rotating radial Timoshenko beam with periodic pulse loading at the free-end

Consideration is given to the dynamic response of a Timoshenko beam under repeated pulse loading. Starting with the basic dynamical equations for a rotating radial cantilever Timoshenko beam clamped at the hub in a centrifugal force field, a system of equations are derived for coupled axial and lateral motions which includes the transverse shear and rotary inertia effects, as well. The hyperbolic wave equation governing the axial motion is coupled with the flexural wave equation governing the lateral motion of the beam through the velocity-dependent skew-symmetric Coriolis force terms. In the analytical formulation, Rayleigh-Ritz method with a set of sinusoidal displacement shape functions is used to determine stiffness, mass and gyroscopic matrices of the system. The tip of the rotating beam is subjected to a periodic pulse load due to local rubbing against the outer case introducing Coulomb friction in the system. Transient response of the beam with the tip deforming due to rub is discussed in terms of the frequency shift and non-linear dynamic response of the rotating beam. Numerical results are presented for this vibro-impact problem of hard rub with varying coefficients of friction and the contact-load time. The effects of beam tip rub forces transmitted through the system are considered to analyze the conditions for dynamic stability of a rotating blade with intermittent rub.     

Modeling of the random packing of a loose layer of polydisperse spherical particles

A method for calculating the loose packing structure of polydisperse spherical particles with a predetermined size distribution function is proposed. The coordinates of the particle centers in the loose layer are determined as the result of random fall of single spheres on a substrate under the action of gravity, assuming the inelastic collision of the spheres and considering the force of their adhesive interaction, and also assuming that the motion of one sphere on the surface of the other is pure slip. Numerical simulation is used to obtain the pattern of arrangement of polydisperse spherical particles in the loose powder layer, whose porosity depends on the particle size distribution function. The results are compared with experimental data.