Next Article in Journal
The Impacts of Gustiness on Air–Sea Momentum Flux
Next Article in Special Issue
Flagellar Cooperativity and Collective Motion in Sperm
Previous Article in Journal
The Coupled Volume of Fluid and Brinkman Penalization Methods for Simulation of Incompressible Multiphase Flows
Previous Article in Special Issue
Remarks on Regularized Stokeslets in Slender Body Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regularized Stokeslets Lines Suitable for Slender Bodies in Viscous Flow

1
Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, UK
2
Department of Mathematics and Statistics, Macquarie University, 192 Balaclava Rd, Macquarie Park, NSW 2113, Australia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fluids 2021, 6(9), 335; https://doi.org/10.3390/fluids6090335
Submission received: 27 August 2021 / Revised: 10 September 2021 / Accepted: 15 September 2021 / Published: 21 September 2021

Abstract

:
Slender-body approximations have been successfully used to explain many phenomena in low-Reynolds number fluid mechanics. These approximations typically use a line of singularity solutions to represent flow. These singularities can be difficult to implement numerically because they diverge at their origin. Hence, people have regularized these singularities to overcome this issue. This regularization blurs the force over a small blob and thereby removing divergent behaviour. However, it is unclear how best to regularize the singularities to minimize errors. In this paper, we investigate if a line of regularized Stokeslets can describe the flow around a slender body. This is achieved by comparing the asymptotic behaviour of the flow from the line of regularized Stokeslets with the results from slender-body theory. We find that the flow far from the body can be captured if the regularization parameter is proportional to the radius of the slender body. This is consistent with what is assumed in numerical simulations and provides a choice for the proportionality constant. However, more stringent requirements must be placed on the regularization blob to capture the near field flow outside a slender body. This inability to replicate the local behaviour indicates that many regularizations cannot satisfy the no-slip boundary conditions on the body’s surface to leading order, with one of the most commonly used blobs showing an angular dependency of velocity along any cross section. This problem can be overcome with compactly supported blobs, and we construct one such example blob, which can be effectively used to simulate the flow around a slender body.

Graphical Abstract

1. Introduction

Slender bodies immersed in fluids are frequently studied in biology, polymer mechanics and colloids. These bodies typically have arc length much larger than their radius and can be difficult to resolve numerically. Let a denote the ratio of the radius of a slender body to its length. In the small a limit, asymptotic theories have been developed for the flow around slender bodies [1,2,3,4,5,6,7,8,9,10]. These theories are called slender-body theories (SBTs), and they expand the system in powers of a or 1 / ln a . Algebraically accurate SBTs have proven to be very accurate and provide exact results for prolate spheroids [3,5,7,11,12,13,14,15,16], while logarithmically accurate SBTs are useful for analytical estimations [17,18,19,20,21,22,23,24]. Many of these theories place the singularity solutions of Stokes equations along the centreline of the body. This singularity representation diverges on the centreline of the body, and it can become difficult to implement numerically.
One approach to overcome these problems is to regularize the singularities by distributing the force over a small blob [25,26,27]. This regularization of singularities was introduced by Cortez [25] with applications to the boundary integral representation for the flow [26,28,29,30]. It has since been extended to describe the flow around slender bodies by using different regularized slender-body theories [31,32,33,34,35,36].
Two major versions of these regularized SBTs exist: one in which both regularized point forces and source dipoles are distributed over the centreline [31,32,36] and one that consists of placing a series of regularized point forces along the centreline of the body [34,35,37,38,39,40]. The former uses asymptotic methods to match the boundary conditions in SBT, while the latter assumes that the regularized singularities will produce the correct flow outside the slender body with a suitable choice of regularization parameter. This is argued from unit analysis but has not been mathematically proven. Here, we investigate the accuracy of the latter representation but note that recent analysis has looked into the accuracy of the former [41].
In this paper, we compare the asymptotic flows around a line of regularized singularities to the flow around a slender body as predicted by slender-body theory. We show that the asymptotic flow far from the body can be captured by any regularization provided that the regularization parameter satisfies a specific equation that depends on the shape of the blob. This justifies the linear relationship between ϵ and the thickness of the slender body, which is typically assumed in numerical simulations. We also show that the flow close to the slender body can only be asymptotically captured for specific types of blobs. As such, not all regularization choices are capable of satisfying the no-slip boundary conditions on the body and so cannot capture the behaviour of the slender body itself. This is an issue for a commonly used power-law decaying blob but not for compactly supported blobs. We construct one such compact blob, which could be effectively used for numerical simulations.
In Section 2, we briefly review Stokes flow with a focus on classical and regularized singularity solutions and slender-body theories. Section 3 then identifies the conditions under which the asymptotic flow far from a line of regularized Stokeslets is the same as that far from a slender body before we identify the conditions necessary to capture the near flow in Section 4. Finally, we look at how these conditions apply to typical blob types in Section 5 before we conclude the paper.

2. Background into Viscous Flows around Slender Bodies

2.1. Stokes Flow and Classical Singularity Solutions

At microscopic scales the ratio of fluid inertia to fluid viscosity, the Reynolds number tends to be very small. In these cases, the behaviour of the resultant slow viscous flow can be accurately described by the incompressible Stokes equations [42,43,44]:
p   +   μ 2 u   +   f = 0 ,
· u = 0 ,
where p is the pressure of the fluid, u is the velocity of the fluid and f is the force per unit volume on the fluid. Throughout this paper, we will set μ   =   1 without any loss of generality. Although these equations are linear and time independent, they display a strong dependence on the geometry of the system considered. As a result, analytical solutions to these equations are rare except in very specific geometries [42,43,44]. Hence, for a typical problem, they must be evaluated numerically or asymptotically. Both of these approaches often use the point force solution to the flow.
The fluid flow from a point force, f   =   F δ ( r ) where δ ( r ) is the three dimensional Dirac delta function supported at the origin, in a viscous fluid is called the Stokeslet and, in free space, generates the following flow:
u ( r )   =   S ( r ) · F = r 2 I + rr 8 π r 3 · F ,
p ( r )   =   P ( r ) · F = r 4 π r 3 · F ,
where r is the position vector, r   =   | r | is its norm, S ( r ) is the Oseen tensor, I is the identity tensor and all tensor product symbols are omitted. The Stokeslet is the Green’s function of the flow and satisfies the Stokes equations everywhere outside the origin at which point it diverges. As a result, it can be used to solve the Stokes equations in two different ways: the boundary integral representation and the representation by fundamental singularities [42,44]. In the boundary integral representation, the Green’s function nature of the Stokeslet is used to convert the equations into an integral over the boundary of the domain. These integrals can then be discretized to solve for the unknowns in the problem in a method called the boundary element method (for more information on this method see Reference [44]).
The representation by fundamental singularities method, on the other hand, seeks to place the Stokeslet and its derivatives outside of the flow region (i.e., within the body) such that the boundary conditions are satisfied [45]. If such a distribution can be found, it is guaranteed to be the solution due to the uniqueness of Stokes flow. In general, the specific singularities needed depend on the geometry of the problem, the region they are distributed over and the motion of the body. However, due to the biharmonic nature of the Stokes equations, singularity representations involving the Stokeslet are consistently found to occur in combination with another singularity solution called the source dipole. The source dipole is proportional to the Laplacian of the Stokeslet and generates the following flow:
u ( r )   =   D ( r ) · A = r 2 I 3 rr 8 π r 5 · A ,
p ( r ) = 0 ,
where A is the strength of the source dipole, and D ( r )   =   2 S ( r ) / 2 is the source dipole tensor. Similar to the Stokeslet, this fundamental singularity satisfies the Stokes equations everywhere except at the origin where it diverges.

2.2. Regularized Singularity Solutions

The divergent nature of the Stokeslets requires special care when integrating over or determining the flow close to the origin. However this integration is necessary for the boundary integral representation, and the nearby flow is often required when constructing representations using fundamental singularities. Cortez proposed the regularization of these singularities to overcome this issue [25]. This was achieved by replacing the Dirac delta function in the point force representation by a suitable mollifier (hereafter called a blob), f ϵ ( r )   =   ϵ 3 f ( r / ϵ ) where f ϵ ( r ) d r   =   1 , that depends on some regularization parameter ϵ . As ϵ 0 , these blobs become the Dirac delta function; thus, the flow from a Stokeslet is recovered. The flow from a specific blob can be found by solving the forced Stokes equations or represented as a convolution with the relevant singularity. For example, the regularized Oseen tensor can be generally expressed as:
S ϵ ( r )   =   S f ϵ   =   S ( r r ) f ϵ ( r ) d r
but a popular regularization takes the following form:
S ϵ ( r )   =   I ( r 2 + 2 ϵ 2 ) + r r 8 π r ϵ 3 , r ϵ   =   r 2 + ϵ 2
when the blob is f ϵ ( r )   =   15 ϵ 4 / ( 8 π r ϵ 7 ) . The flow from these tensors define the regularized Stokeslet.
The substitution of this regularized Stokeslet into the boundary integral representation develops the regularized boundary element method. This method overcomes any singularity issues caused by the Stokeslet but introduces an error related to the regularization parameter ϵ [26]. Careful analysis of this error showed that, for spherically symmetric blobs, the additional error in the boundary integral representation is proportional to ϵ 2 [30].
This error can be understood through the additional flows generated by the regularization process [27]. We previously showed that when r ϵ , the regularized Stokeslet S ϵ can be written as an infinite sum of singularity solutions, and when r ϵ it becomes isotropic. Furthermore, when f is spherically symmetric, the regularized Oseen tensor can be expressed exactly [27] as:
S ϵ ( r )   =   S ( r ) 0 r / ϵ 4 π t 2 f ( t ) d t   +   ϵ 2 D ( r ) 0 r / ϵ 4 π 3 t 4 f ( t ) d t   +   2 I 3 ϵ r / ϵ t f ( t ) d t ,
where we have now expressed it in terms of fundamental singularity solutions to the Stokes equations. The above representation shows that the flow far from a regularized Stokeslet can always be expressed as that of a Stokeslet and a source dipole proportional to ϵ 2 regardless of the choice of blob. Hence, a spherically symmetric blob generally produces additional flows at order ϵ 2 .
This ϵ 2 error can be controlled by keeping ϵ much smaller than all other length scales within a given problem, thereby allowing the benefits of regularization to outweigh any drawbacks. As such, regularized boundary element simulations have become an increasingly popular tool to explore the behaviour of micro-swimmers [28,34,37,38,39,46,47] and have been extended beyond free space singularities [46,48,49]. This process has proved very popular and has been extended to create regularised slender-body theories [31,32,33,35,36,40].

2.3. Classical Slender-Body Theory

The slow viscous hydrodynamics of long thin bodies is important in many systems [2]. However, this slender-body hydrodynamics can be hard to simulate thanks to the very large aspect ratio 1 / a (major/minor lengths) of the body. This large aspect ratio causes boundary element methods to typically require high surface resolutions. Hence, asymptotic methods, called slender-body theories (SBTs), have been developed to overcome these issues [1,2,3,4,5,6,7,10]. These theories seek to expand the governing equations in terms of the large aspect ratio 1 / a to produce simplified models for the hydrodynamics. Broadly, SBTs have been divided into two types: those that expand the system in powers of 1 / ln ( a ) [1,3] and those that expand the system in powers of a [4,5,6,7]. The former is often called resistive force theory and provides useful analytical estimates [17,18,19,20,21,22,23,24], while the latter has been shown to provide highly accurate results in numerous situations [11,12,13,14,15,16].
Probably the most commonly used algebraically accurate SBT was first developed by Keller and Rubinow [4]. This version was developed by matching the flow around an infinite cylinder to that of a line of Stokeslets. These equations were later rederived and extended by Johnson [5] and Gotz [6] by using the representation of fundamental singularities and matched asymptotic expansions. They achieved this by parametrizing the surface of the slender body, x S ( s , θ ) , as:
x S ( s , θ )   =   x 0 ( s )   +   a e r ( s , θ ) ρ 0 ( s ) ,
where s [ 0 , 1 ] is the arclength parameter of the centreline, x 0 ( s ) is the centreline of the body, 0     a ρ 0 ( s )     a is the body radius at s, e t   =   x 0 ˙ / | x 0 ˙ | is the unit tangent vector along the centreline and e r ( s , θ ) is the local radial vector perpendicular to e t (Figure 1). The radial vector, e r , is often represented through the Frenet–Serret coordinates as e r ( s , θ )   =   cos ( θ ) n ^   +   sin ( θ ) b ^ , where n ^ and b ^ are the normal and bi-normal vectors relative to the centreline, respectively. Johnson and Gotz assumed that the leading-order flow around such a body could be captured by a system of Stokeslets with strength α S ( s ) and source dipoles with strength β S   =   1 2 a 2 ρ 0 2 α S placed along the body’s centreline, x 0 ( s ) . This choice of singularity representation comes from the exact solutions known for a prolate spheroid [45] and generates the following flow:
u ( x )   =   0 1 S ( x x 0 ( t ) ) · α S ( t )   +   D ( x x 0 ( t ) ) · β S ( t ) d t .
The uniqueness of Stokes flow tells us that this representation must be the solution if the fluid velocity u ( x ) matches the surface velocity of the body when x   =   x S . They then assumed that the surface velocity, U ( s ) , is uniform along each cross-section (i.e., it only depends on the position along the centreline x 0 ) and performed a matched asymptotic expansion for u ( x S ) in the limit a 0 while keeping s, and θ fixed.
This matched asymptotic expansion divided the system into an outer and an inner region, and it expanded the behaviour in each region in the limit a 0 . The different regions are then matched together by using Van Dyke’s matching rule, and a composite representation formed [50]. We will provide a summary of the major results of this expansion here and refer readers to the thesis of Gotz [6] for the full derivation. In the inner region, ξ   =   ( t     s ) / a is held constant, and the integral kernels expand to:
S ( i ) ( x S x 0 ( t ) ) = I 8 π a δ 0   +   ρ 0 2 e r e r ρ 0 ξ ( e r e t + e t e r ) + ξ 2 e t e t 8 π a δ 0 3   +   O ( 1 ) ,
D ( i ) ( x S x 0 ( t ) ) = I 8 π a 3 δ 0 3     3 ρ 0 2 e r e r ρ 0 ξ ( e r e t + e t e r ) + ξ 2 e t e t 8 π a 3 δ 0 5   +   O ( a 2 ) ,
α S ( t ) = α S ( s )   +   O ( a ) ,
where the superscript ( i ) denotes the inner region expanded form and δ 0 ( t , s , a )   = ξ 2 + ρ 0 2 ( s ) . Similarly, the outer limit fixes t and produces the expanded forms:
S ( o ) ( x S x 0 ( t ) ) · α S ( t ) = 1 8 π α S ( t ) R 0   +   1 8 π R 0 R 0 α S ( t ) R 0 3   +   O ( a ) ,
D ( o ) ( x S x 0 ( t ) ) · β S ( t ) = 0   +   O ( a ) ,
where the superscript ( o ) denotes the outer region expanded term, R 0 ( t , s )   =   x 0 ( s )     x 0 ( t ) , and its norm is R 0   =   | R 0 | . The common part of these two asymptotic forms is found by expanding the outer form in terms of the inner regions variables or equivalently the inner form in terms of the outer regions variables [50], and it is given by:
S ( o ) ( i ) ( x S x 0 ( t ) ) · α S ( t )   =   S ( i ) ( o ) ( x S x 0 ( t ) ) · α S ( t ) = α S ( s ) 8 π | s t |   +   e t e t · α S ( s ) 8 π | s t |   +   O ( a ) ,
D ( o ) ( i ) ( x S x 0 ( t ) ) · β S ( t )   =   D ( i ) ( o ) ( x S x 0 ( t ) ) · β S ( t ) = 0   +   O ( a ) ,
where the superscript ( o ) ( i ) denotes the outer asymptotic representations, Equations (15) and (16), expanded in the inner region variables ξ , and ( i ) ( o ) denotes the inner asymptotic representations, Equations (12) and (13), expanded in the outer region variable t.
Since these regions match a composite representation, the integral can be formed by adding the inner and outer limits, subtracting the common part and then integrating the representation over the entire length. By integrating this composite representation, we find that the flow on the surface of the body is asymptotically given by:
u ( x S )   =   U ( s ) = 0 1 S ( x S x 0 ( t ) ) · α S ( t ) + D ( x S x 0 ( t ) ) · β S ( t ) d t = 0 1 S ( o ) ( x S x 0 ( t ) ) · α S ( t ) + D ( o ) ( x S x 0 ( t ) ) · β S ( t ) d t +   0 1 S ( i ) ( x S x 0 ( t ) ) · α S ( s ) + D ( i ) ( x S x 0 ( t ) ) · β S ( s ) d t   0 1 S ( o ) ( i ) ( x S x 0 ( t ) ) · α S ( s ) + D ( o ) ( i ) ( x S x 0 ( t ) ) · β S ( s ) d t +   O ( a ) = 1 8 π 0 1 α S ( t ) R 0 + R 0 R 0 · α S ( t ) R 0 3 I + e t e t | s t | · α S ( s ) d t +   0 1 S ( i ) ( x S x 0 ( t ) ) · α S ( s ) + D ( i ) ( x S x 0 ( t ) ) · β S ( s ) d t +   O ( a ) ,
where we have explicitly substituted the forms of the outer and common parts of the expansion. The remaining integrals over the inner region expansion can be asymptotically evaluated [6,7] in order to find the following:
8 π U ( s ) = 0 1 α S ( t ) R 0 + R 0 R 0 · α S ( t ) R 0 3 I + e t e t | s t | · α S ( s ) d t +   [ L I + e t e t + I 3 e t e t ] · α S ( s )   +   O ( a ) ,
where L   =   ln [ 4 s ( 1 s ) / a 2 ρ 0 2 ] . The above equation is the classical slender-body theory integral equation for the unknown Stokeslet strengths α S ( t ) . These equations have been used successfully in the modeling of several systems [12,13,14,15,16,51,52,53], and have been rederived from the boundary integral representation [7] and extended to treat bodies of different shapes [54,55]. Recent studies have also studied the effectiveness of these equation in different situations and bounded error on the flow from a given force to O ( a ln a ) [56,57,58].
Once the force per unit length along the slender body is found from Equation (20), the resultant flow can be determined with Equation (11). This predicted flow has been shown to compare well with experiments of rotating helices [11] and should generally be accurate in regions close to the slender body or far away from the slender body. This is because it must asymptotically satisfy the no-slip boundary condition on the surface of the body, and the flow can be represented by a multipole expansion with Stokeslets and source dipoles far away.
Although the classic SBT equations, Equations (11) and (20), have been very successful, the components of the integrand within Equation (20) are individually divergent. Held together, these divergent factors asymptotically cancel, but this makes them difficult to treat numerically. Furthermore, the operator is known to encounter numerical instabilities related to eigenvalues changing sign. Hence, people have created regularized slender-body theories to overcome these issues.

2.4. Regularized Slender-Body Theories

The numerical issues in implementing the classical SBT equations have prompted people to create classical SBT like models using regularized Stokeslets. Two methods have been developed for this process: one where the flow is represented by a line of regularized Stokeslets and regularized source dipoles [31,32,36] and one where the flow is represented by a line of regularized Stokeslets alone [35]. Although both are based on regularized singularities, these methods treat the no-slip boundary condition and the regularization parameter very differently.
The former method, with regularized Stokeslets and regularized source dipoles, seeks to choose the source dipole strength such that the no-slip boundary condition asymptotically applies at the surface of the body. This can be shown with a similar asymptotic expansion, as is performed in the classical case. These methods can leave the regularization parameter free [32] or use it to help with the expansion [31]. In both cases, however, the regularization parameter, ϵ , helps remove the numerical issue with the implementation. The flow predicted by this method, therefore, satisfies the boundary conditions at the surface of the body and far from it but has been shown to exhibit an additional O ( 1 ) difference between it and the classical SBT that cannot be fixed with the choice of regularization parameter [41]. In practice, this difference may be small and likely stems from the fact that the flow outside the body is not described by the unforced incompressible Stokes equations.
The latter method, with a line of regularized Stokeslets alone, tries to use the regularization parameter itself to capture the thickness of the slender body [35]. As such, the regularization parameter is typically chosen to be ϵ   =   k a ρ 0 ( s ) , where k is some arbitrary constant. The strength of the Stokeslets is then set by applying the velocity condition on an effective line in space. Often this line is chosen to be the centreline of the body, although it could be applied on a line along the surface as well. This writes the problem as a line integral similar to that of classical SBT. Similarly to the previous representation, the use of regularized Stokeslets means that the flow predicted by this model will not satisfy the unforced incompressible Stokes equations outside the body for all regularizations. However, the difference between the classical SBT prediction and that from this representation has not been determined.
In the remaining sections of this paper, we will investigate how well a line distribution of regularized Stokeslets replicates the flow around a slender body for a general spherically symmetric blob and regularization parameter. This will provide us with a set of conditions that must be met in order for the far and near field flows to asymptotically reflect that around a slender body. Importantly we show that one common power-law blob cannot satisfy the no-slip condition on the body’s surface, while compactly supported blobs do.

3. Regularizations for the Flow Far from a Slender Body

Far away from the slender body, the flow, predicted by SBT, is well approximated by Equation (11). This is the because it is a leading order multipole representation [11]. Similarly, the flow far from a line of regularized Stokeslets placed along the centreline of the body can be approximated by expanding S ϵ ( x x 0 ( t ) ) in the limit | x | 1 ,   ϵ and integrating over the length. This expansion provides the following:
u ϵ far ( x ) = 0 1 S ( x x 0 ( t ) ) · γ ( t ) d t   +   O 1 | x | 2 +   ϵ 2 0 4 π 3 u 4 f ( u ) d u 0 1 D ( x x 0 ( t ) ) · γ ( t ) d t   +   O ϵ 2 | x | 4 ,
where u ϵ far ( x ) is the far field flow from the line of regularized Stokeslets, and we have used Equation (9). The above approximation assumes that the functions are well behaved enough to interchange limits and integrals. We note that the source dipole term is subdominate to the Stokeslet in the far-field limit. However, we retain both contributions here because the Stokeslet and source dipole generate different flows, and the far-field contributions of both singularities, from the motion of a slender body, are asymptotically determined by Equation (11). If we compare the Stokeslet and source dipole flows, we observe that the above equation approximates Equation (11) if the following holds:
γ ( t ) = α ( t ) ,
2 h ( s ) 2 0 4 π 3 u 4 f ( u ) d u = ρ 0 ( s ) 2 ,
where we have set ϵ   =   a h ( s )   +   O ( a 2 ) . The first condition tells us that, in order for the flow around a line of regularized Stokeslets to replicate that of a slender body, the strength of the regularized Stokeslets must be equal to the force per unit length along the body. This is unsurprising. The second equation, however, provides a relationship between the regularization parameter, ϵ   =   a h ( s ) , the blob used for the regularization, f ( v ) , and the cross section of the body, ρ 0 ( s ) . It shows that, as assumed in the existing literature, the regularization parameter is proportional to the radius of the cross-section and this proportionality constant depends on the fourth moment of the regularizing blob. Provided this moment exists, the far field contributions from any regularized Stokeslet can be made to approximate the Stokeslet and source dipole flow contributions far from a slender body.

4. Regularization’s for the Flow Near a Slender Body

The flow from a line of regularized Stokeslets near the boundary of the slender body can be approximated by using a similar matched asymptotic expansion as that in slender-body theory. This involves expanding the regularised Stokeslet kernel in the limit of small a for an inner region, where ξ   =   ( t s ) / a , s and θ are fixed, and outer region, where t is fixed, and then matching the solutions. The composite representation of the integral is then found by adding the outer and inner regions and subtracting the common part, similarly to before, to find the following:
0 1 S ϵ ( x x 0 ( t ) ) · α ( t ) d t = 0 1 S ϵ ( o ) ( x x 0 ( t ) ) · α ( t ) d t +   0 1 S ϵ ( i ) ( x x 0 ( t ) ) · α ( s ) d t   0 1 S ϵ ( o ) ( i ) ( x x 0 ( t ) ) · α ( s ) d t   +   O ( a ) ,
where the superscripts ( o ) , ( i ) and ( o ) ( i ) again denote the expansion in the outer, inner and the outer subsequently expanded in the inner (common), respectively. In the above, x   =   x 0 ( s )   +   e r a ρ is a point close to the surface of the body, and ρ is a radial distance from the centreline that is of the same order as ρ 0 ( s ) . Setting ρ   =   ρ 0 ( s ) corresponds to the flow on the boundary, while ρ   >   ρ 0 ( s ) represents the flow close to the boundary.
In the previous section, we established that far from the line of regularised Stokeslets, the kernel, S ϵ , behaves as a line distribution of point forces and source dipoles (refer to Equation (9)). As a consequence, the outer region expansion of the regularised kernel is asymptotically equivalent to that of a Stokeslet and a source dipole in the same limits. Hence, the kernel can be expressed as:
S ϵ ( o ) ( x x 0 ( t ) ) · α ( t ) = S ( o ) ( x x 0 ( t ) ) · α ( t ) +   ϵ 2 0 4 π 3 u 4 f ( u ) d u D ( o ) ( x x 0 ( t ) ) · α ( t ) +   O ( a ) = 1 8 π α S ( t ) R 0   +   1 8 π R 0 R 0 · α S ( t ) R 0 3   +   O ( a ) ,
which is identical to the outer expanded kernel in classical SBT. In the above, we have assumed that ϵ is linearly related to a, as suggested by the analysis of the previous section. Since this outer kernel is the same as that found in classical SBT, the expansion of it in terms of the inner region variables will also be the same. The common behaviour can, therefore, be expressed as:
S ϵ ( o ) ( i ) ( x x 0 ( t ) ) · α S ( s ) = S ( o ) ( i ) ( x x 0 ( t ) ) · α S ( s )   +   O ( a ) = α S ( s ) 8 π | s t |   +   e t e t · α S ( s ) 8 π | s t |   +   O ( a ) .
These results can be substituted back into our composite representation for the flow near the surface of the body to find the following:
0 1 S ϵ ( x x 0 ( t ) ) · α ( t ) d t = 1 8 π 0 1 α ( t ) R 0 + R 0 R 0 · α ( t ) R 0 3 I + e t e t | s t | · α ( s ) d t +   0 1 S ϵ ( i ) ( x x 0 ( t ) ) · α ( s ) d t   +   O ( a ) .
The remaining inner region expansion can be deduced from Equation (9), and it takes the following form:
S ϵ ( i ) ( x x 0 ( t ) ) = S ( i ) ( x x 0 ( t ) ) 0 δ / h ( s ) 4 π r 2 f ( r ) d r +   a 2 h 2 ( s ) D ( i ) ( x x 0 ( t ) ) 0 δ / h ( s ) 4 π 3 r 4 f ( r ) d r +   2 I 3 a h ( s ) δ / h ( s ) r f ( r ) d r   +   O ( a ) ,
where ϵ   =   a h ( s ) , δ ( t , s , a )   =   ξ 2 + ρ 2 ( s ) , and we noted that r / ϵ   =   a δ ( t , s , a ) / ϵ   +   O ( a )   =   O ( 1 ) . If we consider a point on the surface δ   =   δ 0 . The above relation expresses the inner region representation of the regularized Stokeslet in terms of the inner region representation of the Stokeslet and source dipole. Similar inner region expansions were stated in the summary of SBT for flow on the surface of the body, Equations (12) and (13). The expanded form of the Stokeslet and source dipole near the surface are structurally equivalent to these but with δ 0 ,   ρ 0 replaced with δ ,   ρ . Explicitly, this provides the following:
S ( i ) ( x x 0 ( t ) ) = I 8 π a δ   +   ρ 2 e r e r ρ ξ ( e r e t + e t e r ) + ξ 2 e t e t 8 π a δ 3 ,
D ( i ) ( x x 0 ( t ) ) = I 8 π a 3 δ 3 3 ρ 2 e r e r ρ ξ ( e r e t + e t e r ) + ξ 2 e t e t 8 π a 3 δ 5 .
These equations can be combined with Equations (28) and (27) to asymptotically determine the flow close to the surface of the slender body from a line of regularised Stokeslets.
The difference in the near body flow from a line of regularized Stokeslets and the flow around a slender body asymptotically predicted by SBT can be determined by subtracting Equation (27) from Equation (19). Since the outer and common limits are the same in each case, the difference between the leading order contributions at x   =   x 0 ( s )   +   a e r ( s ) ρ ( s ) is:
0 1 Δ SBT ϵ d t = 0 1 S ( i ) ( x x 0 ( t ) ) + 1 2 a 2 ρ 0 2 ( t ) D ( i ) ( x x 0 ( t ) ) S ϵ ( i ) ( x x 0 ( t ) ) d t · α ( s ) = 0 1 S ( i ) ( x x 0 ( t ) ) 1 0 δ / h ( s ) 4 π r 2 f ( r ) d r d t · α ( s ) +   1 2 a 2 ρ 0 2 ( s ) 0 1 D ( i ) ( x x 0 ( t ) ) 1 2 h 2 ρ 0 2 0 δ / h ( s ) 4 π 3 r 4 f ( r ) d r d t · α ( s )   0 1 2 I 3 a h ( s ) δ / h ( s ) r f ( r ) d r d t · α ( s ) = 0 1 S ( i ) ( x x 0 ( t ) ) δ / h ( s ) 4 π r 2 f ( r ) d r d t · α ( s ) 2 I 3 a h ( s ) 0 1 δ / h ( s ) r f ( r ) d r d t · α ( s ) + a 2 h 2 ( s ) 0 1 D ( i ) ( x x 0 ( t ) ) δ / h ( s ) 4 π 3 r 4 f ( r ) d r d t · α ( s ) +   a 2 1 2 ρ 0 2 ( s ) h 2 0 4 π 3 r 4 f ( r ) d r 0 1 D ( i ) ( x x 0 ( t ) ) d t · α ( s )
where 0 1 Δ SBT ϵ d t is the difference between the leading order contributions relative to the near field flows from a line of regularized Stokeslets and classical SBT. This difference needs to be O(a) for the flow around a line of regularized Stokeslets to asymptotically replicate the flow predicted by SBT to the same order. Since α ( s ) is typically unknown, this condition must hold for arbitrary α (s). The integrand of the above integral can be expressed as:
Δ SBT ϵ   =   E 1 I   +   E 2 e r e r   +   E 3 e t e t   +   E 4 ( e r e t + e t e r ) ,
where
E 1 = δ 2 a h 3 1 2 3 ϱ 1 2 ϱ 2 1 6 ϱ 4 f ϱ δ h d ϱ   +   1 8 π a δ 3 h ( s ) 2 0 4 π 3 u 4 f ( u ) d u 1 2 ρ 0 ( s ) 2 ,
E 2 = ρ 2 2 a h 3 1 ( ϱ 4 ϱ 2 ) f ϱ δ h d ϱ     3 ρ 2 8 π a δ 5 h ( s ) 2 0 4 π 3 u 4 f ( u ) d u 1 2 ρ 0 ( s ) 2 ,
E 3 = ξ 2 2 a h 3 1 ( ϱ 4 ϱ 2 ) f ϱ δ h d ϱ     3 ξ 2 8 π a δ 5 h ( s ) 2 0 4 π 3 u 4 f ( u ) d u 1 2 ρ 0 ( s ) 2 ,
E 4 = ρ ξ 2 a h 3 1 ( ϱ 2 ϱ 4 ) f ϱ δ h d ϱ   +   3 ρ ξ 8 π a δ 5 h ( s ) 2 0 4 π 3 u 4 f ( u ) d u 1 2 ρ 0 ( s ) 2 .
In the above we used the integral substitution r   =   ϱ h ( s ) / δ to change the lower bounds of the first integrals to one and Equations (29) and (30). These E i values represent the relative error between the flow predicted by SBT and that from the line of regularized Stokeslets in each direction. We require all four terms to integrate to O ( a ) for the flow to match. The integral of E 4 can be shown to be O ( a ) by parity arguments. This is because E 4 is odd in ξ , and in the asymptotic limit a 0 the integral bounds become ± . However, the other three terms provide non-trivial conditions on f ( u ) .
If we assume that Equation (23) is satisfied, then the conditions E 1 , 2 , 3 d t   =   O ( a ) simplify to:
3 h 3 2 ρ 3 0 1 E 1 + 1 3 E 2 + E 3 d t = s / ( a ρ ) ( 1 s ) / ( a ρ ) ( ζ 2 + 1 ) 1 ϱ ϱ 2 f ϱ ρ 1 + ζ 2 h d ϱ d ζ   =   O ( a ) ,
2 h 3 ρ 3 0 1 E 2 d t = s / ( a ρ ) ( 1 s ) / ( a ρ ) 1 ( ϱ 4 ϱ 2 ) f ϱ ρ 1 + ζ 2 h d ϱ d ζ   =   O ( a ) ,
2 h 3 ρ 3 0 1 E 3 d t = s / ( a ρ ) ( 1 s ) / ( a ρ ) ζ 2 1 ( ϱ 4 ϱ 2 ) f ϱ ρ 1 + ζ 2 h d ϱ d ζ   =   O ( a ) ,
where ζ   =   ξ / ρ   =   ( t s ) / ( a ρ ) , and the equality needs to hold for ρ     ρ 0 . The above conditions provide a set of restrictions on the types of blobs, f ( v ) , that can replicate the near field flow around a slender body. These conditions are necessary and sufficient for the line regularized Stokeslets to agree with classical SBT to O(a) near the boundary as the integrals of E i represent the difference between the two representations. These conditions, therefore, specify the types of regularizations that can be used in numerical simulations to asymptotically capture the flow near a slender boy. These conditions can be viewed as an integral transform into ρ , which depends on the tails of the blobs. Although the range over which these equations apply can change from blob to blob, these equations must be true when ρ   =   ρ 0 . If this condition is not satisfied, the no-slip boundary condition on the surface of the body cannot be satisfied. In these cases the solutions found with the line of regularized Stokeslets cannot be assumed to relate to the real solution a priori.

5. Testing the Conditions on Common Blob Types

In the previous sections, we identified one condition on the regularization parameter, Equation (23), required to approximately match the far field flow and three conditions on the blob type, Equations (37)–(39), in order to asymptotically match the flow near the surface of the body. These conditions apply in limit that the body is long and slender and were found by comparing the asymptotic behaviour around a line of regularised Stokeslets to that predicted by SBT. Furthermore, although these latter conditions can hold for a range of ρ , we know they must hold when ρ   =   ρ 0 if the resultant flow is to approximate that of a slender body. In this section, we apply these conditions to three common regularization blob types: power-law blobs, compact blobs and exponential blobs. In each case, we identify if a line of them can be used to represent the flow from a slender body.

5.1. Power-Law Blobs

Probably the most common blob type used within the literature is the power-law type blob. These blobs are given by:
f p ( v )   =   Γ ( n / 2 ) π 3 / 2 Γ ( ( n 3 ) / 2 ) 1 ( v 2 + 1 ) n / 2 ,
where n is a integer greater then five, and Γ ( x ) is the gamma function. We remind the reader that the above blob relates to the regularizing function through f ϵ   =   ϵ 3 f ( r / ϵ ) . The above representation reduces to one of the most popular examples, Equation (8), when n   =   7 .
By inserting this blob into the far field flow condition, Equation (23), we find the following:
h ( s )   =   ρ 0 ( s ) n 5 2 .
This, therefore, sets the regularization the blobs must use to satisfy the far field condition. Significantly, we observe that when n   =   7 , h ( s )   =   ρ 0 ( s ) . This indicates that the regularisation parameter ϵ should be taken as half the thickness of the slender body, as is typically assumed in numerical simulations. Similarly, the near field conditions, Equations (37)–(39), asymptotically evaluate to the following:
s / ( a ρ ) ( 1 s ) / ( a ρ ) ( ζ 2 + 1 ) 1 ϱ ϱ 2 f p ϱ ρ 1 + ζ 2 h d ϱ d ζ = ( n 5 ) n ( n 2 ) ( ( n 5 ) 2 + 4 ρ 2 ) 3 n 2 16 π ( n 3 ) ρ 3 ( n 5 ) n ( ( n 10 ) n 25 ) F 2 ( ρ ) 2 n + 3 ρ n + 2 π ( n 1 ) + ( n 5 ) n F 1 ( ρ ) 2 n ρ n π   +   ( n 5 ) n + 2 F 3 ( ρ ) 2 n + 2 ρ n + 2 π n +   O ( a ) ,
s / ( a ρ ) ( 1 s ) / ( a ρ ) 1 ( ϱ 4 ϱ 2 ) f p r ρ 1 + ζ 2 h d ϱ d ζ = ( n 5 ) n 1 ( ( n 5 ) 2 + 4 ρ 2 ) 5 n 2 32 π ρ 5 + O ( a ) ,
s / ( a ρ ) ( 1 s ) / ( a ρ ) ζ 2 1 ( ϱ 4 ϱ 2 ) f p ϱ ρ 1 + ζ 2 h d ϱ d ζ = ( n 5 ) n 1 ( ( n 5 ) 2 + 4 ρ 2 ) 5 n 2 64 π ρ 5 n 5 n F 4 ( ρ ) 2 n + 1 ρ n π ( n 3 ) + O ( a ) ,
where ρ   =   ρ / ρ 0 , 2 F 1 ( a , b ; c , x ) is the ordinary hypergeometric function, and
F 1 ( ρ ) = 2 F 1 n 2 2 , n 1 2 ; n 2 ; ( n 5 ) 2 4 ρ 2 ,
F 2 ( ρ ) = 2 F 1 n 1 2 , n 1 2 ; n + 1 2 ; ( n 5 ) 2 4 ρ 2 ,
F 3 ( ρ ) = 2 F 1 n 1 2 , n 2 ; n + 2 2 ; ( n 5 ) 2 4 ρ 2 ,
F 4 ( ρ ) = 2 F 1 n 3 2 , n 3 2 ; n 1 2 ; ( n 5 ) 2 4 ρ 2 .
Clearly the above forms are not O ( a ) to leading order; thus, the power law blobs do not formally capture the flow. Furthermore, since the conditions are also not satisfied when ρ   =   1 , the power-law blobs are incapable of capturing the no-slip boundary conditions. However, when these functions are plotted against ρ for different n (Figure 2), we find that they are small for small n and increase as n increases. This suggests that although these representations cannot replicate the boundary conditions directly, when n   =   6 or n   =   7 , the error induced by these terms may be small. This may explain why n   =   7 blobs have seemingly been successfully used in many problems, although in practice further validation is needed to be certain that they actually represent real flow.
Finally, we note that the same conclusion can be drawn by completing the inner region expansion of the resultant regularized Stokeslet. For example, when n   =   7 the regularized Stokeslet from this blob has the form of Equation (8). If we expand this in the inner limit, we find the following:
S ϵ ( i ) ( x x 0 ( t ) )   =   I ( δ 2 + 2 h 2 ) + ρ 2 e r e r ρ ξ ( e r e t + e t e r ) + ξ 2 e t e t 8 π a δ 2 + h 2 3   +   O ( a ) ,
which can then be integrated to give:
8 π 0 1 S ϵ ( i ) ( x x 0 ( t ) ) d t = I L 2 + 2 ρ 2 + 4 h 2 ( s ) ρ 2 + h 2 ( s ) + 2 ρ 2 ρ 2 + h 2 ( s ) e r e r + L 2 e t e t ,
where L 2   =   ln [ 4 s ( 1 s ) / ( a 2 ( ρ 2 + h 2 ) ) ] . The inner structure clearly has a non-zero e r e r term that cannot be eliminated by any choice of h. Hence, this blob cannot remove this angular dependence and cannot satisfy the boundary conditions on the surface of the slender body when ρ   =   ρ 0 . Therefore, it cannot be assumed a priori that the flow from a line of these regularized Stokeslets approximates the flow around a slender body. The presence of this e r e r term is directly related to Equation (43) containing an O(1) contribution. This is because the condition used in Equation (43) comes from the E 2 terms, which contain the difference in the leading order terms from a line of regularised Stokeslets and SBT in the e r e r direction.

5.2. Compact Blobs

Another class of blobs that can be used for regularized singularities are compact blobs. These blobs are zero beyond a certain radius and can significantly simplify the integration. Although such compact blobs regularize the singularity, they often generate discontinuities in the flow; thus, they are not used as often. For example, consider the following blob:
f c ( v )   =   15 32 π ( 21 v 2 11 ) v   <   1 0 v     1 .
This blob isolates the regularization to within the slender body and is chosen such that the total dipole strength is the same as in classical slender-body theory. Given these conditions we expect that the flow outside the body should solve the unforced Stokes equations and identically match the slender-body flow. This can be verified directly with the far field condition, Equation (23), which sets the regularization parameter as
h   =   ρ 0 ,
and the near field conditions become
s / ( a ρ ) ( 1 s ) / ( a ρ ) ( ζ 2 + 1 ) 1 ϱ ϱ 2 f c ϱ ρ 1 + ζ 2 h d ϱ d ζ = 0 ,
s / ( a ρ ) ( 1 s ) / ( a ρ ) 1 ( ϱ 4 ϱ 2 ) f c ϱ ρ 1 + ζ 2 h d ϱ d ζ = 0 ,
s / ( a ρ ) ( 1 s ) / ( a ρ ) ζ 2 1 ( ϱ 4 ϱ 2 ) f c ϱ ρ 1 + ζ 2 h d ϱ d ζ = 0 ,
because r ρ 1 + ζ 2     1 ; thus, f c r ρ 1 + ζ 2   =   0 . Again, we observe that the regularization parameter ϵ equals half the body thickness, similar to those used in simulations; however, this representation also replicates all the near field conditions identically. Hence, a line of these compactly supported blobs could be used to asymptotically determine the flow around a slender body.

5.3. Gaussian Blobs

Finally we consider Gaussian blobs. These blobs decay to zero rapidly but are continuous and so can be thought of similar to an intermediate between the power law and the compact blobs. Gaussian blobs have the following form:
f g ( v )   =   1 π 3 / 2 e v 2 .
The far field condition for these blobs becomes
h   =   ρ 0 ,
similarly to the compact blob and power-law blob with n   =   7 solution. However, in this case the near surface conditions become
s / ( a ρ ) ( 1 s ) / ( a ρ ) ( ζ 2 + 1 ) 1 ϱ ϱ 2 f g ϱ ρ 1 + ζ 2 h d ϱ d ζ = 1 4 π ρ 3 ρ 2 e u u d u   +   O ( a ) ,
s / ( a ρ ) ( 1 s ) / ( a ρ ) 1 ( ϱ 4 ϱ 2 ) f g ϱ ρ 1 + ζ 2 h d ϱ d ζ = e ρ 2 2 π ρ 5   +   O ( a ) ,
s / ( a ρ ) ( 1 s ) / ( a ρ ) ζ 2 1 ( ϱ 4 ϱ 2 ) f g ϱ ρ 1 + ζ 2 h d ϱ d ζ = 1 4 π ρ 5 e ρ 2 ρ 2 ρ 2 e u u d u   +   O ( a ) .
Again, the above evaluation shows that the exponential blobs do not satisfy the conditions; however, this error is exponentially small in ρ 2 with the largest term being 1 / 2 π when ρ   =   1 . This again suggests that the asymptotic error from this regularization could be small; thus, it may still work for simulating the flow around the slender body. However, it would be important to validate the simulation sufficiently beforehand. We note that the difference between the SBT model and the line of Gaussian regularized Stokeslets at the surface can be constructed from the above conditions by using them to reconstruct the E i when ρ   =   1 .

6. Conclusions

In this paper, we compared the asymptotic flow from a line of regularized Stokeslets to that from a slender body. This was possible by comparing SBT with a matched asymptotic expansion of the flow from a generalised line of regularized Stokeslets. We focused on the flow far from the body and near the body surface as the SBT asymptotic representations are expected to be the most valid in these regions. We found that any line of regularized Stokeslets can approximate the flow far from the body provided that the regularization parameter, ϵ , is equal to the thickness of the body, a ρ 0 , multiplied by a constant that depends on the regularization blob, ( 8 π 0 u 4 f ( u ) d u / 3 ) 1 / 2 , Equation (23). This constant can be used for the unknown proportionality constant used in numerical simulations. The asymptotic flow near the surface of a slender body, however, could only be captured by specific regularisation blobs.
The regularized Stokeslets which are capable of asymptotically describing the flow near the surface of a slender body are found to satisfy three conditions (Equations (37)–(39)). These conditions depend on the tails of the regularization blob and are required to apply over a region. Importantly, these conditions must be met at the surface of the slender body if the drag from the line of regularized Stokeslets is to correctly predict the same force and torque on the body from a given motion. If it does not, then the representation cannot be assumed to describe the flow around a slender body. This restricts the types of regularization blobs that could be used when simulating a slender body with a line of regularised Stokeslets.
Finally, we tested these conditions on three blob types: power-law blobs, compact blobs and Gaussian blobs. While all could satisfy the far field condition without issue, we found that power-law blobs and Gaussian blobs could not satisfy the near field conditions anywhere. As such, these blobs cannot match the correct no-slip condition on the bodies surface and so models using these regularization’s may not describe the desired system. Plots of the error suggested that these differences are small; thus, they may not make a large difference practically. The compact blobs investigated were found to satisfy all the conditions without issue. This means that, out of the three blobs considered, the compact blob is the best suited for modelling the behaviour of slender bodies with a line of regularised Stokeslets.
In this article, we restricted ourselves to spherically symmetric regularizations and focused on slender bodies with circular cross sections. It would be of interest to extend the results to non-spherically symmetric regularizations and slender bodies with non-circular cross sections to observe how the conditions change. Furthermore, slender-body theory has also been used in potential flow and diffusion problems; thus, the method could be extended to investigate the accuracy of a line of singularities for a slender body in those systems. Finally, we note that the only regularization we found that satisfied all the conditions was a compact blob; thus, it would be of interest to investigate any non-compact examples that can work.

Author Contributions

Conceptualization, L.K.; methodology, L.K.; formal analysis, B.Z.; investigation, B.Z.; writing—original draft preparation, B.Z.; writing—review and editing, L.K.; supervision, L.K.; project administration, L.K.; funding acquisition, L.K. Both authors have read and agreed to the published version of the manuscript.

Funding

L.K. was funded by Australian Research Council (ARC) under the Discovery Early Career Research Award scheme (grant agreement DE200100168).

Acknowledgments

The authors also thank Eric Lauga for useful discussions and advice.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cox, R.G. The motion of long slender bodies in a viscous fluid Part 1. General theory. J. Fluid Mech. 1970, 44, 791. [Google Scholar] [CrossRef]
  2. Lighthill, J. Flagellar Hydrodynamics: The John von Neumann Lecture, 1975. SIAM Rev. 1976, 18, 161–230. [Google Scholar] [CrossRef]
  3. Batchelor, G.K. Slender-body theory for particles of arbitrary cross-section in Stokes flow. J. Fluid Mech. 1970, 44, 419–440. [Google Scholar] [CrossRef]
  4. Keller, J.B.; Rubinow, S.I. Slender-body theory for slow viscous flow. J. Fluid Mech. 1976, 75, 705–714. [Google Scholar] [CrossRef] [Green Version]
  5. Johnson, R.E. An improved slender-body theory for Stokes flow. J. Fluid Mech. 1979, 99, 411–431. [Google Scholar] [CrossRef]
  6. Götz, T. Interactions of Fibers and Flow: Asymptotics, Theory and Numerics. Ph.D. Thesis, University of Kaiserslautern, Kaiserslautern, Germany, 2000. [Google Scholar]
  7. Koens, L.; Lauga, E. The boundary integral formulation of Stokes flows includes slender-body theory. J. Fluid Mech. 2018, 850, R1. [Google Scholar] [CrossRef] [Green Version]
  8. Andersson, H.I.; Celledoni, E.; Ohm, L.; Owren, B.; Tapley, B.K. An integral model based on slender body theory, with applications to curved rigid fibers. Phys. Fluids 2021, 33, 041904. [Google Scholar] [CrossRef]
  9. Barta, E.; Liron, N. Slender Body Interactions for Low Reynolds Numbers—Part I: Body-Wall Interactions. SIAM J. Appl. Math. 1988, 48, 992–1008. [Google Scholar] [CrossRef]
  10. Koens, L.; Montenegro-Johnson, T.D. Local drag of a slender rod parallel to a plane wall in a viscous fluid. Phys. Rev. Fluids 2021, 6, 064101. [Google Scholar] [CrossRef]
  11. Kim, M.J.; Kim, M.J.; Bird, J.C.; Park, J.; Powers, T.R.; Breuer, K.S. Particle image velocimetry experiments on a macro-scale model for bacterial flagellar bundling. Exp. Fluids 2004, 37, 782–788. [Google Scholar] [CrossRef]
  12. Das, D.; Lauga, E. Computing the motor torque of Escherichia coli. Soft Matter 2018, 14, 5955–5967. [Google Scholar] [CrossRef] [Green Version]
  13. Higdon, J.J.L. A hydrodynamic analysis of flagellar propulsion. J. Fluid Mech. 1979, 90, 685. [Google Scholar] [CrossRef]
  14. Barta, E.; Weihs, D. Creeping flow around a finite row of slender bodies in close proximity. J. Fluid Mech. 2006, 551, 1–17. [Google Scholar] [CrossRef]
  15. Cummins, C.; Seale, M.; Macente, A.; Certini, D.; Mastropaolo, E.; Viola, I.M.; Nakayama, N. A separated vortex ring underlies the flight of the dandelion. Nature 2018, 562, 414–418. [Google Scholar] [CrossRef]
  16. Rodenborn, B.; Chen, C.H.; Swinney, H.L.; Liu, B.; Zhang, H.P. Propulsion of microorganisms by a helical flagellum. Proc. Natl. Acad. Sci. USA 2013, 110, E338–E347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Gray, J.; Hancock, G.J. The Propulsion of Sea-Urchin Spermatozoa. J. Exp. Biol. 1955, 32, 802–814. [Google Scholar] [CrossRef]
  18. Chattopadhyay, S.; Moldovan, R.; Yeung, C.; Wu, X. Swimming efficiency of bacterium Escherichia coli. Proc. Natl. Acad. Sci. USA 2006, 103, 13712–13717. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Koens, L.; Zhang, H.; Moeller, M.; Mourran, A.; Lauga, E. The swimming of a deforming helix. Eur. Phys. J. E 2018, 41, 119. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, X.; Chan, F.K.; Parthasarathy, T.; Gazzola, M. Modeling and simulation of complex dynamic musculoskeletal architectures. Nat. Commun. 2019, 10, 4825. [Google Scholar] [CrossRef] [Green Version]
  21. Becker, L.E.; Koehler, S.A.; Stone, H.A. On self-propulsion of micro-machines at low Reynolds number: Purcell’s three-link swimmer. J. Fluid Mech. 2003, 490, 15–35. [Google Scholar] [CrossRef] [Green Version]
  22. Tătulea-Codrean, M.; Lauga, E. Asymptotic theory of hydrodynamic interactions between slender filaments. Phys. Rev. Fluids 2021, 6, 074103. [Google Scholar] [CrossRef]
  23. Waszkiewicz, R.; Szymczak, P.; Lisicki, M. Stability of sedimenting flexible loops. J. Fluid Mech. 2021, 919, A14. [Google Scholar] [CrossRef]
  24. Man, Y.; Page, W.; Poole, R.J.; Lauga, E. Bundling of elastic filaments induced by hydrodynamic interactions. Phys. Rev. Fluids 2017, 2, 123101. [Google Scholar] [CrossRef] [Green Version]
  25. Cortez, R. The Method of Regularized Stokeslets. SIAM J. Sci. Comput. 2001, 23, 1204–1225. [Google Scholar] [CrossRef]
  26. Cortez, R.; Fauci, L.; Medovikov, A. The method of regularized Stokeslets in three dimensions: Analysis, validation, and application to helical swimming. Phys. Fluids 2005, 17, 031504. [Google Scholar] [CrossRef]
  27. Zhao, B.; Lauga, E.; Koens, L. Method of regularized stokeslets: Flow analysis and improvement of convergence. Phys. Rev. Fluids 2019, 4, 084104. [Google Scholar] [CrossRef]
  28. Smith, D.J. A boundary element regularized Stokeslet method applied to cilia- and flagella-driven flow. Proc. R. Soc. A Math. Phys. Eng. Sci. 2009, 465, 3605–3626. [Google Scholar] [CrossRef] [Green Version]
  29. Gallagher, M.T.; Smith, D.J. The art of coarse Stokes: Richardson extrapolation improves the accuracy and efficiency of the method of regularized stokeslets. R. Soc. Open Sci. 2021, 8, 210108. [Google Scholar] [CrossRef] [PubMed]
  30. Nguyen, H.N.; Cortez, R. Reduction of the Regularization Error of the Method of Regularized Stokeslets for a Rigid Object Immersed in a Three-Dimensional Stokes Flow. Commun. Comput. Phys. 2014, 15, 126–152. [Google Scholar] [CrossRef]
  31. Walker, B.J.; Curtis, M.P.; Ishimoto, K.; Gaffney, E.A. A regularised slender-body theory of non-uniform filaments. J. Fluid Mech. 2020, 899, A3. [Google Scholar] [CrossRef]
  32. Cortez, R.; Nicholas, M. Slender body theory for Stokes flows with regularized forces. Commun. Appl. Math. Comput. Sci. 2012, 7, 33–62. [Google Scholar] [CrossRef] [Green Version]
  33. Cortez, R. Regularized Stokeslet segments. J. Comput. Phys. 2018, 375, 783–796. [Google Scholar] [CrossRef] [Green Version]
  34. Buchmann, A.; Fauci, L.J.; Leiderman, K.; Strawbridge, E.; Zhao, L. Mixing and pumping by pairs of helices in a viscous fluid. Phys. Rev. E 2018, 97, 023101. [Google Scholar] [CrossRef] [Green Version]
  35. Martindale, J.D.; Jabbarzadeh, M.; Fu, H.C. Choice of computational method for swimming and pumping with nonslender helical filaments at low Reynolds number. Phys. Fluids 2016, 28, 021901. [Google Scholar] [CrossRef] [Green Version]
  36. Nguyen, H.; Cortez, R.; Fauci, L. Computing Flows Around Microorganisms: Slender-Body Theory and Beyond. Am. Math. Mon. 2014, 121, 810. [Google Scholar] [CrossRef] [Green Version]
  37. Nguyen, H.; Koehl, M.A.R.; Oakes, C.; Bustamante, G.; Fauci, L. Effects of cell morphology and attachment to a surface on the hydrodynamic performance of unicellular choanoflagellates. J. R. Soc. Interface 2019, 16, 20180736. [Google Scholar] [CrossRef] [Green Version]
  38. Bouzarth, E.L.; Hutson, K.R.; Miller, Z.L.; Saine, M.E. Using Regularized Singularities to Model Stokes Flow: A Study of Fluid Dynamics Induced by Metachronal Ciliary Waves. In An Introduction to Undergraduate Research in Computational and Mathematical Biology. Foundations for Undergraduate Research in Mathematics; Birkhäuser: Cham, Switzerland, 2020; pp. 443–469. [Google Scholar] [CrossRef]
  39. Olson, S.D.; Suarez, S.S.; Fauci, L.J. Coupling biochemistry and hydrodynamics captures hyperactivated sperm motility in a simple flagellar model. J. Theor. Biol. 2011, 283, 203–216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Montenegro-Johnson, T.D.; Koens, L.; Lauga, E. Microscale flow dynamics of ribbons and sheets. Soft Matter 2017, 13, 546–553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Ohm, L. Remarks on Regularized Stokeslets in Slender Body Theory. Fluids 2021, 6, 283. [Google Scholar] [CrossRef]
  42. Kim, S.; Karrila, S.J. Microhydrodynamics: Principles and Selected Applications; Courier Corporation: Boston, MA, USA, 2005; p. 507. [Google Scholar]
  43. Lamb, H. Hydrodynamics, 6th ed.; Cambridge University Press: Cambridge, UK, 1932; p. 605. [Google Scholar]
  44. Pozrikidis, C. Boundary Integral and Singularity Methods for Linearized Viscous Flow; Cambridge University Press: Cambridge, UK, 1992; p. 259. [Google Scholar]
  45. Chwang, A.T.; Wu, T.Y. Hydromechanics of low-Reynolds-number flow. Part 2. Singularity method for Stokes flows. J. Fluid Mech. 1975, 67, 787–815. [Google Scholar] [CrossRef] [Green Version]
  46. Montenegro-Johnson, T.; Smith, A.; Smith, D.; Loghin, D.; Blake, J. Modelling the fluid mechanics of cilia and flagella in reproduction and development. Eur. Phys. J. E Soft Matter 2012, 35, 111. [Google Scholar] [CrossRef] [Green Version]
  47. Montenegro-Johnson, T.D. Microtransformers: Controlled microscale navigation with flexible robots. Phys. Rev. Fluids 2018, 3, 062201. [Google Scholar] [CrossRef] [Green Version]
  48. Ainley, J.; Durkin, S.; Embid, R.; Boindala, P.; Cortez, R. The method of images for regularized Stokeslets. J. Comput. Phys. 2008, 227, 4600–4616. [Google Scholar] [CrossRef]
  49. Cortez, R.; Varela, D. A general system of images for regularized Stokeslets and other elements near a plane wall. J. Comput. Phys. 2015, 285, 41–54. [Google Scholar] [CrossRef]
  50. Hinch, E.J. Perturbation Methods; Cambridge University Press: Cambridge, UK, 1991; p. 160. [Google Scholar]
  51. Lauga, E. Bacterial Hydrodynamics. Annu. Rev. Fluid Mech. 2016, 48, 105–130. [Google Scholar] [CrossRef] [Green Version]
  52. Gaffney, E.A.; Gadêlha, H.; Smith, D.J.; Blake, J.R.; Kirkman-Brown, J.C. Mammalian Sperm Motility: Observation and Theory. Annu. Rev. Fluid Mech. 2011, 43, 501–528. [Google Scholar] [CrossRef] [Green Version]
  53. Katsamba, P.; Michelin, S.; Montenegro-Johnson, T.D. Slender Phoretic Theory of chemically active filaments. J. Fluid Mech. 2020, 898, A24. [Google Scholar] [CrossRef]
  54. Borker, N.S.; Koch, D.L. Slender body theory for particles with non-circular cross-sections with application to particle dynamics in shear flows. J. Fluid Mech. 2019, 877, 1098–1133. [Google Scholar] [CrossRef] [Green Version]
  55. Koens, L.; Lauga, E. Slender-ribbon theory. Phys. Fluids 2016, 28, 013101. [Google Scholar] [CrossRef] [Green Version]
  56. Mori, Y.; Ohm, L. Accuracy of slender body theory in approximating force exerted by thin fiber on viscous fluid. Stud. Appl. Math. 2021, 147, 127–179. [Google Scholar] [CrossRef]
  57. Mori, Y.; Ohm, L.; Spirn, D. Theoretical Justification and Error Analysis for Slender Body Theory. Commun. Pure Appl. Math. 2020, 73, 1245–1314. [Google Scholar] [CrossRef] [Green Version]
  58. Mori, Y.; Ohm, L.; Spirn, D. Theoretical Justification and Error Analysis for Slender Body Theory with Free Ends. Arch. Ration. Mech. Anal. 2020, 235, 1905–1978. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Example geometry of a slender body. x 0 ( s ) is the centreline, x S is the surface, e t ( s ) is the tangent to the centreline, e r ( s ) is a radial vector perpendicular to the tangent and 2 a ρ 0 ( s ) represents the thickness at s.
Figure 1. Example geometry of a slender body. x 0 ( s ) is the centreline, x S is the surface, e t ( s ) is the tangent to the centreline, e r ( s ) is a radial vector perpendicular to the tangent and 2 a ρ 0 ( s ) represents the thickness at s.
Fluids 06 00335 g001
Figure 2. The value of near flow conditions for power-law blobs of n   =   6 (blue), n   =   7 (orange) and n   =   8 (green). (a) Displays the result from condition Equation (37) (Equation (42)), (b) displays the result from condition Equation (38) (Equation (43)) and (c) displays the result from condition Equation (39) (Equation (44)).
Figure 2. The value of near flow conditions for power-law blobs of n   =   6 (blue), n   =   7 (orange) and n   =   8 (green). (a) Displays the result from condition Equation (37) (Equation (42)), (b) displays the result from condition Equation (38) (Equation (43)) and (c) displays the result from condition Equation (39) (Equation (44)).
Fluids 06 00335 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, B.; Koens, L. Regularized Stokeslets Lines Suitable for Slender Bodies in Viscous Flow. Fluids 2021, 6, 335. https://doi.org/10.3390/fluids6090335

AMA Style

Zhao B, Koens L. Regularized Stokeslets Lines Suitable for Slender Bodies in Viscous Flow. Fluids. 2021; 6(9):335. https://doi.org/10.3390/fluids6090335

Chicago/Turabian Style

Zhao, Boan, and Lyndon Koens. 2021. "Regularized Stokeslets Lines Suitable for Slender Bodies in Viscous Flow" Fluids 6, no. 9: 335. https://doi.org/10.3390/fluids6090335

APA Style

Zhao, B., & Koens, L. (2021). Regularized Stokeslets Lines Suitable for Slender Bodies in Viscous Flow. Fluids, 6(9), 335. https://doi.org/10.3390/fluids6090335

Article Metrics

Back to TopTop