scholarly journals Slender body theory for particles with non-circular cross-sections with application to particle dynamics in shear flows

2019 ◽  
Vol 877 ◽  
pp. 1098-1133 ◽  
Author(s):  
Neeraj S. Borker ◽  
Donald L. Koch

This paper presents a theory to obtain the force per unit length acting on a slender filament with a non-circular cross-section moving in a fluid at low Reynolds number. Using a regular perturbation of the inner solution, we show that the force per unit length has $O(1/\ln (2A))+O(\unicode[STIX]{x1D6FC}/\ln ^{2}(2A))$ contributions driven by the relative motion of the particle and the local fluid velocity and an $O(\unicode[STIX]{x1D6FC}/(\ln (2A)A))$ contribution driven by the gradient in the imposed fluid velocity. Here, the aspect ratio ($A=l/a_{0}$) is defined as the ratio of the particle size ($l$) to the cross-sectional dimension ($a_{0}$) and $\unicode[STIX]{x1D6FC}$ is the amplitude of the non-circular perturbation. Using thought experiments, we show that two-lobed and three-lobed cross-sections affect the response to relative motion and velocity gradients, respectively. A two-dimensional Stokes flow calculation is used to extend the perturbation analysis to cross-sections that deviate significantly from a circle (i.e. $\unicode[STIX]{x1D6FC}\sim O(1)$). We demonstrate the ability of our method to accurately compute the resistance to translation and rotation of a slender triaxial ellipsoid. Furthermore, we illustrate novel dynamics of straight rods in a simple shear flow that translate and rotate quasi-periodically if they have two-lobed cross-section, and rotate chaotically and translate diffusively if they have a combination of two- and three-lobed cross-sections. Finally, we show the remarkable ability of our theory to accurately predict the motion of rings, retaining great accuracy for moderate aspect ratios (${\sim}10$) and cross-sections that deviate significantly from a circle, thereby making our theory a computationally inexpensive alternative to other Stokes flow solvers.

1968 ◽  
Vol 35 (2) ◽  
pp. 285-288 ◽  
Author(s):  
N. C. Huang ◽  
C. Y. Sheu

This paper treats the optimal design of a vertical column that is built-in at the lower end. In addition to its own weight, the column is to carry an axial compressive load at its unsupported upper end. The column is to be designed as a thin-walled tube. The median line is to be the same for all cross sections; the wall thickness, though constant along the median line of any cross section, is allowed to vary along the length of the tube. Accordingly, the weight per unit length of the tube is proportional to the bending stiffness. For given length and total weight, the variation of the wall thickness along the column is to be determined to maximize the critical value of the compressive load at the upper end. The influence of a maximum allowable compressive stress on the design is also investigated.


We consider the problem of designing the section of a cylinder to minimize the drag per unit length it experiences when placed perpendicular to a uniform stream at low Reynolds number; we suppose the area of the cross-section to be given, and the flow to be two-dimensional. The relevant properties of a cylinder of general cross-section in a particular orientation can conveniently be expressed in terms of its equivalent radius; when the drag and flow at infinity are parallel, this equivalent radius is the radius of the circular cylinder giving rise to the same drag per unit length. We obtain a variational formula for this equivalent radius when the surface of the cylinder is perturbed; this shows that the optimum profile we seek must be such that the flow past it has a vorticity of constant magnitude at its surface, and this fact enables the optimum to be determined analytically. The efficacy of a particular section may be measured by its effective radius, this being the equivalent radius when the length scale is chosen to give the section an area π ; thus a circular cylinder has an effective radius of 1. The minimum possible effective radius, achieved by the optimum profile, is 0.88876. To illustrate some of the arguments we exploit in a more familiar setting, we also obtain a variational formula for the drag on a three-dimensional body in Stokes flow when its surface is perturbed.


2018 ◽  
Vol 840 ◽  
pp. 379-407 ◽  
Author(s):  
Neeraj S. Borker ◽  
Abraham D. Stroock ◽  
Donald L. Koch

A ring with a cross-section that has a blunt inner and sharper outer edge can attain an equilibrium orientation in a Newtonian fluid subject to a low Reynolds number simple shear flow. This may be contrasted with the continuous rotation exhibited by most rigid bodies. Such rings align along an orientation when the rotation due to fluid vorticity balances the counter-rotation due to the extensional component of the simple shear flow. While the viscous stress on the particle tries to rotate it, the pressure can generate a counter-vorticity torque that aligns the particle. Using boundary integral computations, we demonstrate ways to effectively control this pressure by altering the geometry of the ring cross-section, thus leading to alignment at moderate particle aspect ratios. Aligning rings that lack fore–aft symmetry can migrate indefinitely along the gradient direction. This differs from the periodic spatial trajectories of fore–aft asymmetric axisymmetric particles that rotate in periodic orbits. The mechanism for migration of aligned rings along the gradient direction is elucidated in this work. The migration speed can be controlled by varying the cross-sectional shape and size of the ring. Our results provide new insights into controlling motion of individual particles and thereby open new pathways towards manipulating macroscopic properties of a suspension.


Author(s):  
Mahmoud Alidadi ◽  
Sander Calisal

A numerical study is conducted to calculate the squat for a wigley hull. An approach based on slender body theory is used to convert the three dimensional ship problem into a series of two dimensional problems in cross sections from bow to stern (solved sequentially in time). A boundary element method is used to compute the flow potential at every cross section. The ship squat is calculated from the pressure integration over the hull. Numerical results for the Wigley hull is presented and compared with the experimental results.


2015 ◽  
Vol 785 ◽  
pp. 372-400 ◽  
Author(s):  
Michael Zabarankin ◽  
Olga M. Lavrenteva ◽  
Avinoam Nir

The deformation of an immiscible toroidal drop embedded in axisymmetric compressional Stokes flow is analysed via the boundary integral formulation in the case of equal viscosity. Numerical simulations are performed for the drop having initially the shape of a torus with circular cross-section. The quasi-stationary dynamic simulations reveal that, when the viscous forces, proportional to the intensity of the flow, are relatively weak compared with the surface tension (the ratio of these forces is characterized by the capillary number,$Ca$), three different scenarios of drop evolution are possible: indefinite expansion of the liquid torus, contraction to the centre and a stationary toroidal shape. When the intensity of the flow is low, the stationary shapes are shown to be close to circular tori. Once the outer flow strengthens, the cross-section of the stationary torus assumes first an elliptic and then an egg-like shape. For the capillary number greater than a critical value,$Ca_{cr}$, toroidal stationary shapes were not found. Remarkably,$Ca_{cr}$is close to the critical capillary number found previously for a simply connected drop flattened in compressional flow. Thus, a new example of non-uniqueness of stationary drop shape in viscous flow is obtained. Approximate stationary solutions in the form of tori with circular and elliptic cross-sections are obtained by minimizing the normal velocity over the drop interface. They are shown to be in good agreement with the stationary shapes from quasi-dynamic simulations for the corresponding intervals of the capillary number.


1970 ◽  
Vol 44 (3) ◽  
pp. 419-440 ◽  
Author(s):  
G. K. Batchelor

A rigid body whose length (2l) is large compared with its breadth (represented by R0) is straight but is otherwise of arbitrary shape. It is immersed in fluid whose undisturbed velocity, at the position of the body and relative to it, may be either uniform, corresponding to translational motion of the body, parallel or perpendicular to the body length, or a linear function of distance along the body length, corresponding to an ambient pure straining motion or to rotational motion of the body. Inertia forces are negligible. It is possible to represent the body approximately by a distribution of Stokeslets over a line enclosed by the body; and then the resultant force required to sustain translational motion, the net stresslet strength in a straining motion, and the resultant couple required to sustain rotational motion, can all be calculated. In the first approximation the Stokeslet strength density F(x) is independent of the body shape and is of order μUε, where U is a measure of the undisturbed velocity and ε = (log 2l/R0)−1. In higher approximations, F(x) depends on both the body cross-section and the way in which it varies along the length. From an investigation of the ‘inner’ flow field near one section of the body, and the condition that it should join smoothly with the ‘outer’ flow which is determined by the body as a whole, it is found that a given shape and size of the local cross-section is equivalent, in all cases of longitudinal relative motion, to a circle of certain radius, and, in all cases of transverse relative motion, to an ellipse of certain dimensions and orientation. The equivalent circle and the equivalent ellipse may be found from certain boundary-value problems for the harmonic and biharmonic equations respectively. The perimeter usually provides a better measure of the magnitude of the effect of a non-circular shape of a cross-section than its area. Explicit expressions for the various integral force parameters correct to the order of ε2 are presented, together with iterative relations which allow their determination to the order of any power of ε. For a body which is ‘longitudinally elliptic’ and has uniform cross-sectional shape, the force parameters are given explicitly to the order of any power of ε, and, for a cylindrical body, to the order of ε3.


1974 ◽  
Vol 64 (4) ◽  
pp. 817-826 ◽  
Author(s):  
N. J. De Mestre ◽  
D. F. Katz

Stokes flow is analysed for a combination body, consisting of a sphere attached to a slender body, translating along its axis in an infinite and otherwise un-disturbed fluid. The cross-section of the after-body, or tail, is circular; the radius, while not necessarily constant, is small compared with the radius of the spherical head. The tail is represented by a distribution of Stokeslets of strength per unit length F(z), located and directed along its axis. The interactive effect of head-tail attachment is manifested by the presence of image singularities located within the sphere. The image system for a single tail Stokeslet must be such that the no-slip condition is satisfied on the surface of the sphere. It is shown that this system consists of a Stokeslet, a Stokes doublet (stresslet only) and a source doublet located a t the image point. The strength F(z) is obtained by applying the no-slip condition to the combination body. The solution follows the lines of traditional slender-body theory, an expansion being performed in ascending powers of the reciprocal of the logarithm of the aspect ratio. The integral force parameters and F(z) are obtained to second order. The interactive effect is assessed, and the results are discussed in the context of a sedimenting micro-organism, such as a spermatozoon. The drag on the combination body is shown to be less by around 10% than the sum of the drags on an isolated sphere and tail. This drag, for a sperm-shaped body, is divided approximately equally between head and tail.


The semi-classical treatment of atom-atom collisions involving electronic transitions is discussed. As is well know n difficulties occur if the classical trajectories associated with the various states of importance in a collision process differ significantly. A method designed to overcome these is described. It will be referred to as the forced-common-turning-point method. The four coupled first-order differential equations which describe the new version of the semi-classical two-state treatment for an atom-atom collision may be reduced to a pair of generalized impact parameter equations. The first Born approximation to the cross-section obtained from the straightforward semiclassical treatment differs from the corresponding cross-section obtained from the full quantal treatment mainly in that it contains an anomalous multiplying factor equal to the ratio of the initial to the final velocity of relative motion. This anomaly does not arise with the forced-common-turning-point method. A model collision process which provides a very searching test is considered. Only two states are included. The initial interaction is zero, the final interaction is Coulombic and the transition matrix element is exponential. Curve-crossing may occur. The distorted wave approximation to the excitation cross-section may be found exactly and may also be computed using the forced-common-turning-point method. There is remarkable accord between the results. Thus in a case w here the reduced m ass of the colliding systems is 2 on the chemical scale, w here the excitation energy is 3.4 eV and where the incident kinetic energy of relative motion is only 0.85 eV above this the excitation cross-sections obtained differ by as little as 0.01 %; and, moreover, the patterns of the contributions to the cross-sections from the separate partial waves are similar.


2017 ◽  
Vol 830 ◽  
pp. 369-386 ◽  
Author(s):  
Stoffel D. Janssens ◽  
Vikash Chaurasia ◽  
Eliot Fried

We perform a static analysis of a circular cylinder that forms a barrier between surfactant-laden and surfactant-free portions of a liquid–gas interface. In addition to determining the general implications of the balances for forces and torques, we quantify how the imbalance $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FE}=\unicode[STIX]{x1D6FE}_{a}-\unicode[STIX]{x1D6FE}_{b}$ between the uniform surface tension $\unicode[STIX]{x1D6FE}_{a}$ of the surfactant-free portion of the interface and the uniform surface tension $\unicode[STIX]{x1D6FE}_{b}$ of the surfactant-laden portion of the interface influences the load-bearing capacity of a hydrophobic cylinder. Moreover, we demonstrate that the difference between surface tensions on either side of a cylinder with a cross-section of arbitrary shape induces a horizontal force component $f^{h}$ equal to $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FE}$ in magnitude, when measured per unit length of the cylinder. With an energetic argument, we show that this relation also applies to a rod-like barrier with cross-sections of variable shape. In addition, we apply our analysis to amphiphilic Janus cylinders and we discuss practical implications of our findings for Marangoni propulsion and surface pressure measurements.


1980 ◽  
Vol 99 (2) ◽  
pp. 411-431 ◽  
Author(s):  
Robert E. Johnson

The present study examines the flow past slender bodies possessing finite centre-line curvature in a viscous, incompressible fluid without any appreciable inertia effects. We consider slender bodies having arbitrary centre-line configurations, circular transverse cross-sections, and longitudinal cross-sections which are approximately elliptic close to the body ends (i.e. prolate-spheroidal body ends). The no-slip boundary condition on the body surface is satisfied, using a convenient stepwise procedure, to higher orders in the slenderness parameter (ε) than has previously been possible. In fact, the boundary condition is satisfied up to an error term of O(ε2) by distributing appropriate stokeslets, potential doublets, rotlets, sources, stresslets and quadrupoles on the body centre-line. The methods used here produce an integral equation valid along the entire body length, including the ends, whose solution determines the stokeslet strength or equivalently the force per unit length up to a term of O(ε2). The O(ε2) correction to the stokeslet strength is also found. The theory is used to examine the motion of a partial torus and a helix of finite length. For helical bodies comparisons are made between the present theory and the resistive-force theory using the force coefficients of Gray & Hancock and Lighthill. For the motion considered the Gray & Hancock force coefficients generally underestimate the force per unit length, whereas Lighthill's coefficients provide good agreement except in the vicinity of the body ends.


Sign in / Sign up

Export Citation Format

Share Document