Next Article in Journal
Anisotropic RANS Turbulence Modeling for Wakes in an Active Ocean Environment
Next Article in Special Issue
Modeling Heavy-Gas Dispersion in Air with Two-Layer Shallow Water Equations
Previous Article in Journal
Computational Aerodynamics Analysis of Non-Symmetric Multi-Element Wing in Ground Effect with Humpback Whale Flipper Tubercles
Previous Article in Special Issue
Lightning Solvers for Potential Flows
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stokes Equation in a Semi-Infinite Region: Generalization of the Lamb Solution and Applications to Marangoni Flows

Laboratoire Ondes et Matière d’Aquitaine (UMR 5798), CNRS, University of Bordeaux, F-33400 Talence, France
*
Author to whom correspondence should be addressed.
Fluids 2020, 5(4), 249; https://doi.org/10.3390/fluids5040249
Submission received: 1 December 2020 / Revised: 12 December 2020 / Accepted: 15 December 2020 / Published: 18 December 2020
(This article belongs to the Special Issue Recent Advances in Fluid Mechanics: Feature Papers)

Abstract

:
We consider the creeping flow of a Newtonian fluid in a hemispherical region. In a domain with spherical or nearly spherical geometry, the solution of the Stokes equation can be expressed as a series of spherical harmonics. However, the original Lamb solution is not complete when the flow is restricted to a semi-infinite space. The general solution in hemispherical geometry is then constructed explicitly. As an application, we discuss the solutions of Marangoni flows due to a local source at the liquid–air interface.

1. Introduction

At small length scales, the creeping flow of a Newtonian fluid is described by the linear Stokes equation. Since the latter is directly related to the Laplace equation, it admits a wide variety of analytical solutions that are well documented and constitute the subject of classical textbooks [1,2]. Applications are found in everyday life, ranging from the actuation of micro-organisms [3] to the draining of liquid foams [4] or the wetting dynamics of thin films [5], to name only a few.
Despite these facts, there are still numerous flow problems involving the Stokes equation that remain unresolved. In particular, flows driven by interfacial stresses, such as Marangoni flows, can be especially tenacious [6,7]. The first reason for that comes from the presence of the interface itself, which breaks the symmetry in the direction perpendicular to it. The second reason lies in the fact that Marangoni flows are usually coupled to transport phenomena, leading to intricate relations even in the linear regime [8,9]. It was also suggested in recent experimental studies that the presence of surface-active impurities at the interface might induce hydrodynamic instabilities at vanishing Reynolds numbers [10,11]. There is thus a strong need for a more thorough description of Marangoni flows.
Here, we focus more specifically on the Stokes equation in spherical geometry. The general solution, which was first derived by Lamb [12], can be written as a series of Legendre functions. The situation gets more involved when the flow is restricted to a semi-infinite region. In this case, the solution of Lamb is not complete anymore, and several authors have pointed out that additional terms have to be considered [13,14]. The aim of this work was thus to derive in a systematic manner the solution of the Stokes equation in hemispherical geometry. The paper is organized as follows. In Section 2, we first remind the reader of the main properties of the Lamb solution. We then discuss in Section 3 the thermocapillary flow due to a point-like heat source at the interface. This simple example shows unambiguously that the Lamb solution happens to be incomplete when considering a hemispherical domain. The general solution is then constructed in Section 4, and some applications to thermocapillary flows are given in Section 5. The main results are finally summarized in Section 6.

2. The Stokes Equation in Spherical Coordinates

2.1. Lamb Solution

We consider the viscous flow of an incompressible fluid in a domain with spherical, or nearly spherical, geometry. We note r the radial coordinate, θ the polar angle ( 0 θ π ) and φ the azimuthal angle ( 0 φ < 2 π ). In the limit of the vanishing Reynolds number, the velocity and pressure fields are governed by the Stokes equation
η 2 v = p ,
together with the continuity equation
· v = 0 .
Due to the linearity of Equations (1) and (2), the general solution can be written as the sum of a homogeneous and a particular solution, v = v H + v P .
  • The homogeneous solution v H satisfies the equations 2 v H = 0 and · v H = 0 . It can thus be expressed as the linear combination of a potential and a toroidal field:
    v H = Φ + × r χ ,
    where Φ and χ are both harmonic functions, 2 Φ = 2 χ = 0 .
  • The particular solution v P is related to the theory of harmonic functions as well. By taking the divergence of the equation 2 v p = p together with · v p = 0 , it is straightforward to show that the pressure field satisfies the Laplace equation, 2 p = 0 .
In spherical coordinates, the solution of the Laplace equation 2 ψ = 0 can be expressed as a series of spherical harmonics. It can be written quite generally as ψ ( r , θ , φ ) = l = ψ l ( r , θ , φ ) , with
ψ l ( r , θ , φ ) = m = l l A ( l m ) r l Y l m ( θ , φ ) .
The spherical harmonics of degree l and order m are defined (up to a numerical prefactor) by
Y l m ( θ , φ ) = P l m ( cos θ ) e i m φ ,
where the associated Legendre polynomials can be expressed, for instance, by Rodrigues formula:
P l m ( x ) = ( 1 ) l + m 2 l l ! 1 x 2 | m | / 2 d l + | m | d x l + | m | 1 x 2 l for l 0 .
The latter relation states in particular that Y l m = 0 if | m | > l . The definition of Y l m can be extended to negative degrees thanks to the equality Y l m = Y l 1 m . An essential property of spherical harmonics is that they form an orthogonal basis for scalar functions on the sphere surface.
Coming back to the Stokes equation, we can implement the same decomposition as in Equation (4) for the potential, toroidal and pressure fields, respectively. It is then possible to show that
v ( r ) = l = ( l + 3 ) r 2 p l 2 η ( l + 1 ) ( 2 l + 3 ) r l p l η ( l + 1 ) ( 2 l + 3 ) + Φ l + × r χ l .
This expression is usually named after Lamb [12]. It is therefore a series in powers (for l > 0 ) or inverse powers (for l < 0 ) of the distance r. As a matter of fact, the convention l l 1 is often assumed for the exterior flows in order to avoid negative indices. The components of the velocity v = v r e r + v θ e θ + v φ e φ are finally obtained by applying the differential operators in Equation (7).

2.2. Explicit Expression for Exterior Flows

The Lamb solution is relevant to describe the motion of the fluid inside or outside a region with spherical symmetry. In this work, we focus specifically on exterior flows and assume that all fields vanish when r . The pressure can then be conveniently expressed as
p = η l 1 m π ( l m ) r l + 1 2 ( 2 l 1 ) ( l + 1 ) Y l m ,
with l m l . The components of the velocity v = v r e r + v θ e θ + v φ e φ are given by [1,2]
(9) v r = l 1 m π ( l m ) r l Y l m + l 3 m ρ ( l 2 , m ) r l Y l 2 m , (10) v θ = l 1 m π ( l m ) r l ( l 2 ) l ( l + 1 ) s c Y l m + l 3 m ρ ( l 2 , m ) r l 1 ( l 1 ) s c Y l 2 m + l 2 m σ ( l 1 , m ) r l i m s Y l 1 m , (11) v φ = l 1 m π ( l m ) r l ( l 2 ) l ( l + 1 ) i m s Y l m l 3 m ρ ( l 2 , m ) r l 1 ( l 1 ) i m s Y l 2 m + l 2 m σ ( l 1 , m ) r l s c Y l 1 m ,
where we have introduced c = cos θ and s = sin θ . Note that our notation differs slightly from that usually found in the literature, since we anticipate the generalization to a semi-infinite region.
As an illustration, let us consider the axisymmetric solutions that decay from the origin as the inverse of the distance. The radial and polar components of velocity field then read v r ( r , θ ) = Φ r ( 0 ) ( θ ) / r and v θ ( r , θ ) = Φ θ ( 0 ) ( θ ) / r , where the functions Φ r ( 0 ) and Φ θ ( 0 ) are given by
Φ r ( 0 ) ( θ ) = π ( 10 ) P 1 0 , and Φ θ ( 0 ) ( θ ) = π ( 10 ) 2 s c P 1 0 ,
with P 1 0 = cos θ and s c P 1 0 = sin θ . The integration constant π ( 10 ) has then to be set by enforcing the relevant boundary conditions.

3. A First Encounter with Hemispheric Legendre Functions

So far, we have made no restriction regarding the polar angle θ . The issue gets more involved, however, when it is limited to part of the interval [ 0 , π ] . Such a situation typically occurs when the fluid is bounded by a flat surface or interface. This encompasses a large number of physical realizations, ranging from a droplet placed on a solid surface to the motion of a camphor scrap at the water–air interface. In a semi-infinite domain, Lamb’s solution is not complete anymore, and several authors have recently derived additional terms to complement Equation (7). To illustrate our point, we first discuss a classic problem of fluid mechanics: the steady-state Marangoni flow due to a local heat source at the liquid–air interface [13].

3.1. Thermocapillary Flow Due to a Point Source

We study the stationary heat flow induced by a source located at z = 0 . This can be achieved, for instance, with a colloidal particle heated by a laser beam [15]. Assuming that the size of the source is small compared to the other relevant length-scales, heat may be regarded as emerging from a point source q ( r ) = Q δ ( r ) , with Q the injected power ( [ Q ] = W ), in a semi-infinite liquid bounded by a flat interface. The z-axis is perpendicular to the interface and oriented downward; the liquid phase corresponds to z > 0 . The system is schematically shown in Figure 1. We assume that the temperature variations are small enough so that both the viscosity and the fluid density can be regarded as constant.
The discussion is restricted to the linearized version of the transport equations; i.e., we focus on the limits of vanishing Reynolds and Péclet numbers. The temperature field is then solution of the heat equation
ρ c t T = κ 2 T + Q δ ( r ) ,
with κ being the thermal conductivity ( [ κ ] = W · m 1 · K 1 ), ρ the mass density ( [ ρ ] = kg · m 3 ) and c the specific heat capacity ( [ c ] = J · K 1 · kg 1 ) of the fluid. Here, δ is the Dirac delta function ( [ δ ( r ) ] = m 3 in 3D). We focus in this work on the stationary version of the heat equation
2 T = Q κ δ ( r ) .
The results presented hereafter are thus valid in the long time limit r D t , with D = κ / ( ρ c ) being the thermal diffusion coefficient. The velocity field satisfies the incompressible Stokes equation (Equations (1) and (2)). Regarding the boundary conditions, it is assumed that there is neither mass nor heat flux across the interface
v · e z = 0 and T · e z = 0 at z = 0 ,
whereas both the velocity and the temperature return to their equilibrium values far away from the disturbance
lim r v ( r ) = 0 and lim r T ( r ) = T 0 .
Finally, since the surface tension is usually a decreasing function of the temperature, the actuation of the fluid arises from stresses at the interface. This is expressed by the Marangoni boundary condition, which relates the shear stress and the gradient of interfacial tension
1 e z e z · σ · e z γ | z = 0 = 0 ,
with σ the hydrodynamic stress tensor (whose Cartesian components are σ i j = p δ i j + η i v j + j v i ). For moderate deviations with respect to the equilibrium temperature, one can assume a linearized relationship for the surface tension: γ ( T ) = γ 0 γ T ( T T 0 ) , with γ 0 = γ ( T 0 ) . The coefficient γ T = | γ / T | characterizes the rate of change of surface tension with respect to temperature.

3.2. Solution of the Transport Equations

The mathematical problem defined by Equations (14)–(17) is solved in the semi-infinite region z > 0 . First, the temperature field can be derived by analogy with electrostatics. Indeed, Equation (14) is equivalent to the Poisson equation for the electrostatic potential due to a point charge. One therefore gets
T ( r ) = T 0 + Q 2 π κ r .
Note that the numerical coefficient differs from the usual value ( 2 π vs. 4 π ) since heat diffusion is restricted to a half-space. Note also that this solution is only valid down to a cut-off distance set by the size of the physical heat source.
Next, Equation (17) suggests that the velocity fields follows the same power law as the temperature. We can therefore search for solutions of the Stokes equations (Equations (1) and (2)) in the form v r ( r , θ ) = Φ r ( θ ) / r and v θ ( r , θ ) = Φ θ ( θ ) / r . After some algebra, one obtains
Φ r ( θ ) = Q γ T 4 π κ η ( 1 2 cos θ ) , and Φ θ ( θ ) = Q γ T 4 π κ η cos θ sin θ 1 + cos θ .
At first sight, the connection between the θ -dependence in Equation (19) and the generic Lamb solution (Equation (12)) is not striking. Still, it can be noticed that the function Φ r ( θ ) can also be written as
Φ r ( θ ) = α P 1 0 ( c ) α 2 ,
with P 1 0 ( c ) = c = cos θ and α = Q γ T / ( 2 π κ η ) . Similarly, Φ θ ( θ ) can be re-expressed as
Φ θ ( θ ) = α 2 s c P 1 0 + α 2 1 c 1 + c ,
with s c P 1 0 = s = sin θ . In both expressions, the first term involves an associated Legendre function (or its derivative) with l = 1 and m = 0 , as expected from Lamb’s solution (12). The second term, which is not present in Equation (12), is more puzzling. Let us focus, for instance, on the expression (21) for Φ θ : although the second term is well defined in the upper half-space 0 c 1 , it is obviously singular in the limit c 1 . Since the original Lamb solution only involves functions that are regular in the full range 1 c 1 , this contribution is necessarily missing from the general solution. As we shall see in the following, the additional terms in Equations (20) and (21) can actually be defined as hemispheric Legendre functions.

4. Laplace and Stokes Equations in a Semi-Infinite Region

As argued in Section 2, the solution of the Stokes equation is closely related to the theory of harmonic functions. In this section, we first focus on Laplace equation in order to introduce generalized spherical harmonics. The derivation follows closely that of reference [16]. We then discuss how to extend Lamb’s solution in a semi-infinite domain.

4.1. Method of Separation of Variables Applied to Laplace Equation

Let us consider the Laplace equation 2 ψ ( r ) = 0 satisfied by a real scalar function ψ . The position vector r = ( r , θ , φ ) is restricted to the upper half-space z 0 , which corresponds to the interval 0 θ π / 2 . Results pertaining to the lower half-space z 0 can be deduced in a straightforward manner thanks to the symmetry θ π θ .
In spherical polar coordinates, the Laplace equation reads
2 ψ = 1 r 2 r r 2 ψ r + 1 r 2 sin θ θ sin θ ψ θ + 1 r 2 sin 2 θ 2 ψ φ 2 = 0 .
This equation can be solved using the method of separation of variables. To that end, we seek a general solution of the form ψ ( r , θ , φ ) = f ( r ) g ( θ ) h ( φ ) , where f, g and h are functions of a single variable. Equation (22) then leads to a set of ordinary differential equations
(23) d d r r 2 d f d r λ f = 0 , (24) 1 sin θ d d θ sin θ d g d θ + λ μ sin 2 θ g = 0 , (25) d 2 h d φ 2 + μ h = 0 ,
where λ and μ are two (yet unknown) constants.
We first consider the radial Equation (23). Without loss of generality, we can set λ = l ( l + 1 ) . Searching for power-law solutions, one readily gets
f ( r ) = A r l + B r l 1 ,
where A and B are two integration constants. In the particular case l = 1 / 2 , the general solution (26) has to be replaced by f ( r ) = A r 1 / 2 + B r 1 / 2 ln r .
We next focus on the azimuthal Equation (25). On physical grounds, one requires the function h to be single-valued. Since ( r , θ , φ ) and ( r , θ , φ + 2 π ) represent the same location in space, the condition h ( φ + 2 π ) = h ( φ ) demands that μ = m 2 , with m being an integer. We write the solution of Equation (25) as
h ( φ ) = e i m φ , with m Z .

4.2. Hemispheric Legendre Functions

Regarding Equation (24), it can be recast in a more usual form by setting c = cos θ . Since we focus on the semi-infinite region 0 θ π / 2 , the variable c is always positive: 0 c 1 . Equation (24) is then equivalent to Legendre differential equation
( 1 c 2 ) d 2 g d c 2 2 c d g d c + l ( l + 1 ) m 2 1 c 2 g ( c ) = 0 .
In the situations commonly encountered in physics, e.g., in quantum mechanics or in fluid mechanics, one is interested in solutions of Equation (28) that are regular over the entire interval c [ 1 , 1 ] . These standard solutions are the well-known associated Legendre functions P l m ( c ) defined in Section 2. In this work, however, the regularity condition is only required in the semi-infinite region I + = [ 0 , 1 ] . Additional solutions, which can be singular on the complementary interval I = [ 1 , 0 ] , are therefore permitted. The essential point is that the extra singularities that may exist in I are physically irrelevant, since we exclusively focus on the interval I + . Although this issue is fairly well established in classical physics, e.g., in the context of the sharp point effect in electrostatics [17], we consider it worth it to remind the reader of the main steps of the derivation—see Appendix A. The general solution of Equation (28) is therefore g ( c ) = P l m ( c ) , where the hemispheric Legendre functions are defined as
P l m ( c ) = 1 c 1 + c | m | / 2 F l , l + 1 ; 1 + | m | ; 1 c 2 , c [ 0 , 1 ] ,
where the hypergeometric function F is given by the series
F ( α , β ; γ ; x ) = n = 0 ( α ) n ( β ) n ( γ ) n x n n ! ,
with ( q ) n being the Pochhammer symbol: ( q ) n = q ( q + 1 ) ( q + n 1 ) for n > 0 , and ( q ) 0 = 1 . The main properties of the hypergeometric function are listed in Appendix A.
In Equation (29), the degree l can assume any real value. Working in a semi-infinite region thus implies a rather unconventional definition of the associated Legendre functions. For 0 | m | l , they coincide (up to a numerical prefactor) with the “usual” functions P l m ( c ) defined in (6)—see Table 1. However, since the series (30) actually converges for all c [ 0 , 1 ] , there is no restriction regarding the order m anymore. Solutions with l < | m | do exist as well, although this is normally forbidden in R 3 . This assertion is of paramount importance, since these additional terms play a central role in accounting for interfacial phenomena such as the Marangoni effect. It can be noticed that the generalized Legendre functions with | m | > l are proportional to ( 1 + c ) | m | / 2 . These functions are indeed regular in the semi-infinite region 0 c 1 , but would be otherwise singular over the whole interval 1 c 1 . To illustrate our point, the first hemispheric Legendre functions are listed in Table 1.
Finally, the general solution of Laplace equation 2 ψ ( r ) = 0 can be written in spherical coordinates as
ψ ( r , θ , φ ) = d l m = A ( l m ) r l Y l m ( θ , φ ) .
Here, we define the hemispheric spherical harmonics
Y l m ( θ , φ ) = P l m ( c ) e i m φ ,
with P l m ( c ) being given in Equation (29). The first difference between the solution (4) in R 3 and Equation (31) in the semi-infinite region 0 c 1 is that the sum now extends to any value of m, be it smaller or larger than the degree l. The second difference is that the degree l varies continuously, the sum being replaced by an integral. Note also that the functions Y l m (or P l m ) are in general not orthogonal on the half-sphere 0 c 1 (although orthogonality is preserved regarding the azimuthal angle φ ).

4.3. Generalized Solution for Exterior Flows

We now have all the elements to extend the Lamb solution to a hemispheric region. Two strategies are possible. The first is to directly transpose the solution given by Equations (8)–(11), provided the correspondence with Y l m Y l m and
l l 0 m = l l 0 d l m = ,
with l 0 = 1 , 2 or 3. As we shall see, this approach works for all but a few terms of the series. The second strategy consists of solving explicitly the Stokes equation. This gives of course the same results, but with the advantage that singular contributions can be handled consistently. Details regarding the algebra are given in Appendix B. Focusing as previously on the exterior solution, we find that the pressure and the components of the velocity can be written as
p ( r , θ , φ ) = 0 d l m = 1 r l + 1 p l m ( c ) e i m φ ,
v r ( r , θ , φ ) = 0 d l m = 1 r l u l m ( c ) e i m φ ,
v θ ( r , θ , φ ) = 0 d l m = 1 r l v l m ( c ) e i m φ ,
v φ ( r , θ , φ ) = 0 d l m = 1 r l w l m ( c ) e i m φ ,
where we get
(38) p l m ( c ) = π ( l m ) 2 ( 2 l 1 ) ( l + 1 ) P l m . (39) u l m ( c ) = π ( l m ) P l m + ρ ( l 2 , m ) P l 2 m , (40) v l m ( c ) = π ( l m ) ( l 2 ) l ( l + 1 ) s c P l m + ρ ( l 2 , m ) 1 ( l 1 ) s c P l 2 m + σ ( l 1 , m ) i m s P l 1 m , (41) w l m ( c ) = π ( l m ) ( l 2 ) l ( l + 1 ) i m s P l m ρ ( l 2 , m ) 1 ( l 1 ) i m s P l 2 m + σ ( l 1 , m ) s c P l 1 m ,
Equations (34)–(41) are the main outcomes of this work. They express the exterior solution of a Stokes flow as a straightforward generalization of the original Lamb solution. Still, there are a few terms in the series that are singular and should therefore be handled separately.

4.3.1. Polar and Azimuthal Components of the Velocity for l = 1 and m 0

It can be noticed in Equations (40) and (41) that both the polar and azimuthal components of the velocity are not properly defined for l = 1 . To regularize the expressions, we need to solve explicitly the incompressible Stokes equations. This can be achieved by appealing the recurrence relations of Legendre functions. The main steps of the derivation are given in Appendix B, so that we directly get the results
v 1 m ( c ) = π ( 1 m ) 2 s c P 1 m + ( 1 + | m | ) ρ ( 1 , m ) P 1 m s + σ ( 0 , m ) i m s P 0 m ,
w 1 m ( c ) = π ( 1 m ) 2 φ P l m s + i s m ρ ( 1 , m ) P 0 m ( 1 + | m | ) c P 1 m + σ ( 0 , m ) s c P 0 m .

4.3.2. Polar Component of the Velocity for l = 1 and m = 0

As a matter of fact, Equation (42) is still singular when the order m is equal to zero. To remove this very last singularity, we proceed as previously and solve directly the differential equations—see Appendix B. We obtain the explicit expression
v 10 ( c ) = π ( 10 ) 2 s ρ ( 1 , 0 ) 1 c 1 + c ,
Regarding the azimuthal component for m = 0 , one can deduce from Equation (41) that w l 0 ( θ , φ ) = σ ( l 1 , 0 ) s c Y l 1 0 , so that finally w 10 ( θ , φ ) = 0 .

5. Application of the Generalized Lamb Solution to Marangoni Flows

Bringing everything together, we finally reach the conclusion that the generalized exterior solution for the Stokes flow is given by Equations (34)–(41), supplemented with Equations (42)–(44) for l = 1 . We now apply these results to thermocapillary flows.

5.1. Back to the First Encounter

In Section 3, we discussed the thermocapillary flow due to a point source at the liquid–air interface. At the vanishing Péclet number, the temperature field is found to decay as r 1 ; see Equation (18). The Marangoni boundary condition then suggests that the velocity field assumes the same functional form: v r ( r , θ ) = Φ r ( θ ) / r and v θ ( r , θ ) = Φ θ ( θ ) / r . Keeping only the terms proportional to r 1 in the generalized exterior solution, we can directly write
(45) Φ r ( θ ) = π ( 10 ) c + ρ ( 1 , 0 ) , (46) Φ θ ( θ ) = π ( 10 ) 2 s ρ ( 1 , 0 ) 1 c 1 + c .
If we enforce the condition of vanishing normal velocity at the interface, i.e., v θ ( r , π 2 ) = 0 , then the integration constants are related through π ( 10 ) = 2 ρ ( 1 , 0 ) . The solution then reads
(47) Φ r ( θ ) = ρ ( 1 , 0 ) ( 1 2 c ) , (48) Φ θ ( θ ) = ρ ( 1 , 0 ) c s 1 + c ,
which is precisely the result given in Equation (19). The remaining constant can be set by enforcing the stress continuity condition (17). We get ρ ( 1 , 0 ) = Q γ T / ( 4 π κ η ) , as expected.

5.2. Thermocapillary Flow with Dipolar Symmetry

As a second example, we consider the self-propulsion of a heated particle at the water–air interface. This question has recently attracted much attention, both theoretically and experimentally [14,15,18]. The particle acts as a point-like heat source. It is assumed that the temperature profile in the liquid phase has a dipolar symmetry
T ( r ) = T 0 + Q 2 π κ 1 r + b · r r 3 ,
with b = b e x . The resulting thermocapillary flow is solution of the Stokes equations (Equations (1) and (2)) together with the Marangoni boundary condition (17). The explicit expression for the velocity field was obtained in reference [14] after tedious calculations. Here, we can get it directly from the generalized Lamb solution. According to the superposition principle, the velocity field can be written as
v ( r ) = v 10 ( r ) + b v 21 ( r ) ;
the axisymmetric term v 10 has been discussed in the previous paragraph. It is then straightforward to get the dipolar contribution
(51) u 2 , 1 ( c ) = π ( 2 , 1 ) P 2 1 ( c ) + ρ ( 0 , 1 ) P 0 1 ( c ) = 1 2 s 2 ρ ( 0 , 1 ) + π ( 2 , 1 ) 2 ρ ( 0 , 1 ) c π ( 2 , 1 ) c 3 , (52) v 2 , 1 ( c ) = s ρ ( 0 , 1 ) c P 0 1 + i s σ ( 1 , 1 ) P 1 1 ( c ) = i σ ( 1 , 1 ) 2 ρ ( 0 , 1 ) 1 1 + c , (53) w 2 , 1 ( c ) = i s ρ ( 0 , 1 ) P 0 1 ( c ) + s σ ( 1 , 1 ) c P 1 1 = i ρ ( 0 , 1 ) 1 1 + c i σ ( 1 , 1 ) 2 c .
This solution matches exactly with the one that was calculated previously, see Equation (4.6) in reference [14].

6. Conclusions

To summarize, we have derived a generalization of Lamb’s solution for the Stokes flow in a semi-infinite domain. Equations (38)–(44) are the main results of this work. This new solution is relevant, for instance, when the liquid phase is bounded by a flat interface, or in the case of a hemispherical liquid drop moving on an inhomogeneous surface. The range of applications is therefore relatively wide.
We emphasize that, although our results look quite similar to the original Lamb solution, they present several features that are rather unconventional. First of all, the degree l is not limited to integers but is in general a real number. Second, the order m is not limited to the interval l m l but can take values larger than the degree l. These peculiarities are directly related to the fact that the flow is restricted to a semi-infinite region.
The expression of the velocity components (38)–(44) can be used directly for any problem with hemispherical geometry, exactly as one would proceed with Lamb’s solution. As shown in Section 5, it is especially well suited for interfacial phenomena involving the Marangoni effect. We thus expect that this formalism might shed a new light on several important issues in interfacial science, including the spreading of surfactant from a point source or the actuation of Marangoni surfers.

Author Contributions

Conceptualization, methodology, validation, writing—original draft preparation, writing—review and editing, project administration and funding acquisition: T.B. Investigation and formal analysis: G.K. and T.B. All authors have read and agreed to the published version of the manuscript.

Funding

We gratefully acknowledge the University of Bordeaux and its “Initiative d’Excellence” program for financial support.

Acknowledgments

The authors wish to thank S. Michelin, F. Nadal and A. Würger for many useful and fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. From Legendre to Gauss Hypergeometric Equation

In this appendix, we follow reference [16] to construct the general solutions of the associated Legendre equation in a semi-infinite region
( 1 c 2 ) d 2 g d c 2 2 c d g d c + l ( l + 1 ) m 2 1 c 2 g ( c ) = 0 .
First, it should be noticed that Equation (A1) admits three regular singular points at c = 1 , c = 1 and c = + , respectively. As a consequence, it is always possible to recast Equation (A1) into Gauss hypergeometric equation [16]
x ( 1 x ) y ( x ) + γ ( α + β + 1 ) x y ( x ) α β y ( x ) = 0 ,
where α , β and γ are three real (yet unknown) parameters. A fundamental solution of Equation (A2) is given by the hypergeometric function F ( α , β ; γ ; x ) , which is defined as
F ( α , β ; γ ; x ) = n = 0 ( α ) n ( β ) n ( γ ) n x n n ! = Γ ( γ ) Γ ( α ) Γ ( β ) n = 0 Γ ( α + n ) Γ ( β + n ) Γ ( γ + n ) x n n ! ,
with ( q ) n = q ( q + 1 ) ( q + n 1 ) the Pochhammer symbol [ ( q ) 0 = 1 ], and Γ ( x ) = 0 t x 1 e t d t the gamma function. The hypergeometric function verifies the symmetry relation:
F ( α , β ; γ ; x ) = F ( β , α ; γ ; x ) .
The radius of convergence of the series is R = 1 . The series reduces to a polynomial of degree l N when either α or β is equal to l .
The solution of Equation (A2) being known, the issue is to find the mathematical transformation that relates Equations (A1) and (A2). To this end, we make the substitution g ( c ) = χ ( x ) y ( x ) , with χ ( x ) being an auxiliary function such that y ( x ) is solution of Equation (A2). We also define the new variable x = ( 1 c ) / 2 . Inserting in Equation (A1), one gets an intermediate differential equation for y
x ( 1 x ) χ y + 2 x ( 1 x ) χ + ( 1 2 x ) χ y + x ( 1 x ) χ + ( 1 2 x ) χ + l ( l + 1 ) m 2 4 x ( 1 x ) χ y = 0 .
On comparing the second and first-order terms with those in Equation (A2), one can infer that the term in square brackets must be proportional to χ . Without loss of generality, one can always absorb any overall constant factor in the definition of y. We thus obtain a first equation for χ
2 x ( 1 x ) χ + ( 1 2 x ) χ = γ ( α + β + 1 ) x χ ,
whose solution reads
χ ( x ) = x ( γ 1 ) / 2 ( 1 x ) ( α + β γ ) / 2 .
Then, comparison of the zeroth-order terms in Equations (A2) and (A5) imposes that the term in curly brackets must be proportional to χ as well
x ( 1 x ) χ + ( 1 2 x ) χ + l ( l + 1 ) m 2 4 x ( 1 x ) χ = α β χ .
Inserting the solution (A7) into Equation (A8), we can deduce that
γ 1 = ϵ 1 | m | , α + β γ = ϵ 2 | m | , and α β = l ( l + 1 ) .
Here we define ϵ i = ± 1 . The problem therefore admits four independent solutions. For the discussion, it is convenient to switch back to the original variable c.
  • If ϵ 1 = ϵ 2 = + 1 , then the auxiliary function reads χ ( c ) = 1 4 ( 1 c 2 ) | m | / 2 . This is the canonical definition usually adopted in the literature (see, e.g., reference [16]) in order to account for the regularity condition over the whole interval c [ 1 , 1 ] .
  • If ϵ 1 = + 1 and ϵ 2 = 1 , then we get
    χ ( c ) = 1 c 1 + c | m | / 2 .
    This function has a singularity at the south pole c = 1 , but is regular otherwise. It is thus relevant in the upper half-space.
  • If ϵ 1 = 1 and ϵ 2 = + 1 , we obtain χ ( c ) = ( 1 + c ) / ( 1 c ) | m | / 2 . The singularity is now located at the north pole c = 1 . This configuration being complementary to the previous one, it is suitable in the lower half-space.
  • If ϵ 1 = ϵ 2 = 1 , the auxiliary function χ ( c ) = 4 ( 1 c 2 ) | m | / 2 is singular at both poles. It is therefore inappropriate with regard to the semi-infinite domains discussed in this work.
In the upper half-space, the relevant choice is then ϵ 1 = + 1 and ϵ 2 = 1 . It follows from Equation (A9) that the three Gauss parameters are
α = l , β = l + 1 , and γ = | m | + 1 .
Finally, the solutions of Equation (A1) are the generalized Legendre functions, defined for c [ 0 , 1 ] by
P l m ( c ) = 1 c 1 + c | m | / 2 F l , l + 1 ; | m | + 1 ; 1 c 2 .
In the complementary interval c [ 1 , 0 ] , the solutions are obtained thanks to the substitution c c . It can be noted that the generalized Legendre functions satisfy the condition
P ( l + 1 ) m ( c ) = P l m ( c ) .
They also verify the following properties: (i) they are eigenfunctions of the differential operator L
L P l m ( c ) ( 1 c 2 ) c 2 2 c c m 2 ( 1 c 2 ) P l m ( c ) = l ( l + 1 ) P l m ( c ) .
Note that one also has
L c P l m = ( l 2 + l + 2 ) c P l m + 2 ( 1 c 2 ) d P l m d c .
They follow the recurrence relations
(A16) 1 c 2 d P l m d c = l | m | P l 1 m ( c ) l c P l m ( c ) , (A17) 1 c 2 d P l m d c = l + | m | + 1 P l + 1 m ( c ) + ( l + 1 ) c P l m ( c ) .

Appendix B. Solutions of the Stokes Problem in a Semi-Infinite Region

The aim of this second appendix is to derive explicitly the solution of the Stokes equations (Equations (1) and (2)) in the semi-infinite region defined by c = cos θ 0 . The parameter s = sin θ = 1 c 2 is always positive since θ [ 0 , π ] . Note also that θ = s c and d s / d c = c / s . In this representation, Equation (2) reads
1 r 2 r r 2 v r 1 r c ( s v θ ) + 1 s r φ v φ = 0 .
The projection of the Stokes equations on the spherical basis e r , e θ , e φ leads to a set of partial differential equations that are listed, for instance, in [19]. We then look for solutions using the Anzats suggested in Equations (34)–(37).

Appendix B.1. Pressure Field

The pressure p ( r ) satisfies Laplace equation Δ p = 0 , where we define the scalar Laplacian
Δ ψ = r 1 r 2 r ψ + r 2 c s 2 c ψ + ( r s ) 2 φ 2 ψ .
By inserting the Ansatz p ( r , θ , φ ) = r ( l + 1 ) p l m ( c ) e i m φ , one obtains the Legendre differential equation for each mode p l m ( c )
( 1 c 2 ) d 2 p l m d c 2 2 c d p l m d c + l ( l + 1 ) m 2 1 c 2 p l m ( c ) = 0 .
Following the discussion of Appendix A, it is straightforward to conclude that
p l m ( c ) = 2 ( 2 l 1 ) l + 1 π ( l , m ) P l m ( c ) ,
with π ( l m ) being the integration constant. Here, the prefactor 2 ( 2 l 1 ) / ( l + 1 ) is introduced for later convenience.

Appendix B.2. Radial Component of the Velocity

The radial projection of the Stokes Equation (1) reads
Δ v r 2 r 2 v r + 2 r 2 c s v θ 2 s r 2 φ v φ = 1 η r p .
Then, eliminating v θ and v φ thanks to Equation (A18), one gets
Δ v r 2 r 2 v r + 2 r 3 r r 2 v r = 1 η r p .
We then assume the functional form v r ( r , θ , φ ) = r l u l m ( c ) e i m φ , so that each mode satisfies the following equation
( 1 c 2 ) d 2 u l m d c 2 2 c d u l m d c + ( l 2 ) ( l 1 ) m 2 1 c 2 u l m ( c ) = 2 ( 2 l 1 ) p l m ( c ) .
One recognizes on the left-hand side the Legendre equation of degree l 2 : the homogeneous solution is thus u l m ( h ) ( c ) P l 2 m ( c ) . Moreover, since p l m ( c ) P l m ( c ) , a particular solution can be searched under the form u l m ( p ) ( c ) = α P l m ( c ) , with α a constant. It is then straightforward to get α = π ( l m ) , so that finally
u l m ( c ) = π ( l , m ) P l m ( c ) + ρ ( l 2 , m ) P l 2 m ( c ) .

Appendix B.3. Polar Component of the Velocity

The projection of the Stokes Equation (1) on the polar direction leads to
Δ v θ 1 s 2 r 2 v θ 2 s r 2 c v r 2 c s 2 r 2 φ v φ = s η r c p .
Defining v ˜ θ = s v θ and eliminating v φ thanks to Equation (A18), it can be rewritten as
s 2 c 2 v ˜ θ 2 c c v ˜ θ + 1 s 2 φ 2 v ˜ θ + r r 2 r v ˜ θ = s 2 r η c p + 2 s 2 c v r 2 c r r r 2 v r .
We then proceed along the same lines and look for solutions of the form v θ = r l v l m ( c ) e i m φ . Setting v ˜ l m = s v l m yields
( 1 c 2 ) d 2 v ˜ l m d c 2 2 c d v ˜ l m d c + l ( l 1 ) m 2 1 c 2 v ˜ l m ( c ) = ( 1 c 2 ) d p l m d c + 2 ( 1 c 2 ) d u l m d c + 2 ( l 2 ) c u l m ( c ) .
The left-hand side (LHS) corresponds to Legendre equation of degree l 1 , so that the homogeneous solution is v ˜ l m ( h ) ( c ) P l 1 m ( c ) . Coming back to the original function we get v l m ( h ) ( c ) P l 1 m ( c ) / s , which happens to be singular at c = 1 when m = 0 . However, since one can always search for a particular solution that is regular at this point (except for the case l = 1 ; see below), the homogeneous solution has to be regular as well. To remove the singularity, one is naturally led to set
v l m ( h ) ( c ) = i m σ ( l 1 , m ) P l 1 m ( c ) / s .
Note that the homogeneous term can also be written as v l m ( h ) ( c ) e i m φ = σ ( l 1 , m ) φ Y l 1 m / s , with Y l m ( c , φ ) = P l m ( c ) e i m φ , which is the form generally used in Lamb’s solution.
We next focus on the inhomogeneous equation. According to the previous results, the right-hand side (RHS) is a linear combination of P l m and P l 2 m . Due to the linearity of Equation (A28), the particular solution can be searched under the form v ˜ l m ( p ) = v ˜ l m ( p 1 ) + v ˜ l m ( p 2 ) , with
v ˜ l m ( p 1 ) = K 1 ( 1 c 2 ) d P l m d c , and v ˜ l m ( p 2 ) = K 2 ( 1 c 2 ) d P l 2 m d c .
The two constants K 1 π ( l m ) and K 2 ρ ( l m ) have yet to be determined. Thanks to the recurrence relations of the Legendre functions Equation (A16), we can write
v ˜ l m ( p 1 ) = K 1 ( l | m | ) P l 1 m l c P l m ,
and
v ˜ l m ( p 2 ) = K 2 ( l 1 ) c P l 2 m ( l + m 1 ) P l 1 m .
Inserting into the LHS of Equation (A28), one finds LHS = LHS 1 + LHS 2 , with
(A33) LHS 1 = 2 l K 1 ( 1 c 2 ) d P l m d c ( l + 1 ) c P l m , (A34) LHS 2 = 2 ( l 1 ) K 2 ( 1 c 2 ) d P l 2 m d c + ( l 2 ) c P l 2 m .
From the solution (A21) for p l m and (A25) for u l m , the RHS of Equation (A28) can also be written as RHS = RHS 1 + RHS 2 , with
(A35) RHS 1 = 2 l 2 l + 1 π ( l m ) ( 1 c 2 ) d P l m d c ( l + 1 ) c P l m , (A36) RHS 2 = 2 ρ ( l 2 , m ) ( 1 c 2 ) d P l 2 m d c + ( l 2 ) c P l 2 m .
Comparing (A33) and (A35) on the one hand, and (A34) and (A36) on the other hand, one finally obtains
K 1 = l 2 l ( l + 1 ) π ( l m ) , and K 2 = 1 l 1 ρ ( l 2 , m ) .
We can then draw the conclusion that
v l m ( c ) = s l 2 l ( l + 1 ) π ( l m ) d P l m d c + ρ ( l 2 , m ) l 1 d P l 2 m d c + i m s σ ( l 1 , m ) P l 1 m ( c ) .

Appendix B.4. Azimuthal Component of the Velocity

Although one can proceed in the same way for the azimuthal projection, it is more convenient to deduce the component v φ directly from Equation (A18). Defining v φ = r l w l m ( c ) e i m φ , the modes w l m are given by
w l m ( c ) = i s m d s v l m d c + ( l 2 ) u l m ( c ) .
Using (A25), (A38) and the definition of Legendre equation, we get after some algebra
w l m ( c ) = i m s ( l 2 ) l ( l + 1 ) π ( l m ) P l m ( c ) + ρ ( l 2 , m ) l 1 P l 2 m ( c ) + s σ ( l 1 , m ) d P l 1 m d c .

Appendix B.5. Singular Cases

It can be noticed that the expressions of v l m and w l m are not defined for l = 1 ; see Equations (A38) and (A40). These singular cases can be handled as follows. Using the recurrence relation (A17), Equation (A36) can be rewritten as
RHS 2 = 2 ρ ( l 2 , m ) ( l | m | 2 ) P l 3 m ,
so that the corresponding particular solution can be searched under the form v ˜ l m ( p 2 ) = K 2 P l 3 m . After some algebra one gets
K 2 = l | m | 2 2 l 3 ρ ( l 2 , m ) .
This expression is now regular for l = 1 . We can thus write
v 1 m ( c ) = π ( 1 m ) 2 s c P 1 m + ( 1 + | m | ) ρ ( 1 , m ) P 1 m s + σ ( 0 , m ) i m s P 0 m , m 0 .
However, the latter relation has still to be regularized for m = 0 . To do so, we set σ ˜ ( 00 ) = i m σ ( 0 , m ) m = 0 and replace in Equation (A43) the generalized Legendre functions by their explicit expression
v 10 ( c ) = π ( 10 ) 2 s c 2 + 2 ρ ( 1 , 0 ) π ( 10 ) c + 2 σ ˜ ( 0 , 0 ) π ( 10 ) 1
For this expression to be regular, c = 1 must be a root of the term in squared brackets; i.e.,
= c 1 c + 2 ρ ( 1 , 0 ) π ( 10 ) + 1 .
This finally sets the value σ ˜ ( 0 , 0 ) = ρ ( 1 , 0 ) , so that
v 10 ( c ) = 1 2 π ( 10 ) 1 c 2 ρ ( 1 , 0 ) 1 c 1 + c .

References

  1. Happel, J.; Brenner, H. Low Reynolds Number Hydrodynamics, with Special Applications to Particulate Media; Kluver Academic Publishers: The Hague, The Netherlands, 1983. [Google Scholar]
  2. Kim, S.; Karrila, S.J. Microhydrodynamics: Principles and Selected Applications; Dover: New York, NY, USA, 2005. [Google Scholar]
  3. Lauga, E.; Powers, T.R. The hydrodynamics of swimming microorganisms. Rep. Prog. Phys. 2009, 72, 096601. [Google Scholar] [CrossRef]
  4. Cantat, I. Liquid meniscus friction on a wet plate: Bubbles, lamellae, and foams. Phys. Fluids 2013, 25, 031303. [Google Scholar] [CrossRef]
  5. Bertin, V.; Niven, J.; Stone, H.A.; Salez, T.; Raphaël, E.; Dalnoki-Veress, K. Symmetrization of Thin Freestanding Liquid Films via a Capillary-Driven Flow. Phys. Rev. Lett. 2020, 124, 184502. [Google Scholar] [CrossRef] [PubMed]
  6. Scriven, L.E. The Marangoni effects. Nature 1960, 187, 186. [Google Scholar] [CrossRef]
  7. Manikantan, H.; Squires, T.M. Surfactant dynamics: Hidden variables controlling fluid flows. J. Fluid Mech. 2020, 892, P1. [Google Scholar] [CrossRef]
  8. Bickel, T.; Loudet, J.-C.; Koleski, G.; Pouligny, B. Hydrodynamic response of a surfactant-laden interface to a radial flow. Phys. Rev. Fluids 2019, 4, 124002. [Google Scholar] [CrossRef] [Green Version]
  9. Bickel, T. Effect of surface-active contaminants on radial thermocapillary flows. Eur. Phys. J. E 2019, 42, 131. [Google Scholar] [CrossRef] [PubMed]
  10. Mizev, A. Influence of an adsorption layer on the structure and stability of surface tension driven flows. Phys. Fluids 2005, 17, 122107. [Google Scholar] [CrossRef]
  11. Koleski, G.; Vilquin, A.; Loudet, J.-C.; Bickel, T.; Pouligny, B. Azimuthal instability of the radial thermocapillary flow around a hot bead trapped at the water–air interface. Phys. Fluids 2020, 32, 092108. [Google Scholar] [CrossRef]
  12. Lamb, H. Hydrodynamics, 6th ed.; Dover: New York, NY, USA, 1932. [Google Scholar]
  13. Bratukhin, Y.K.; Maurin, L.N. Thermocapillary convection in a fluid filling a half-space. J. Appl. Math. Mech. 1967, 31, 577. [Google Scholar] [CrossRef]
  14. Würger, A. Thermally driven Marangoni surfers. J. Fluid Mech. 2014, 752, 589. [Google Scholar] [CrossRef] [Green Version]
  15. Girot, A.; Danné, N.; Würger, A.; Bickel, T.; Ren, F.; Loudet, J.-C.; Pouligny, B. Motion of optically heated spheres at the water-air interface. Langmuir 2016, 32, 2687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Seaborn, J.B. Hypergeometric Functions and Their Applications; Springer: New York, NY, USA, 1991. [Google Scholar]
  17. Jackson, J.D. Classical Electrodynamics, 3rd ed.; Wiley: New York, NY, USA, 1999. [Google Scholar]
  18. Masoud, H.; Stone, H.A. A reciprocal theorem for Marangoni propulsion. J. Fluid Mech. 2014, 741, R4. [Google Scholar] [CrossRef] [Green Version]
  19. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics, 2nd ed.; Pergamon Press: Oxford, UK, 1987. [Google Scholar]
Figure 1. Schematic representation of the system. The temperature gradient induces horizontal stresses that drive the fluid from the hot to the cold areas of the interface. The z-axis is oriented downward.
Figure 1. Schematic representation of the system. The temperature gradient induces horizontal stresses that drive the fluid from the hot to the cold areas of the interface. The z-axis is oriented downward.
Fluids 05 00249 g001
Table 1. Expressions of the first hemispheric Legendre functions P l m ( c ) for c > 0 .
Table 1. Expressions of the first hemispheric Legendre functions P l m ( c ) for c > 0 .
m = 0 m = 1 m = 2 m = 3
l = 0 1 1 c 1 + c 1 / 2 1 c 1 + c 1 c 1 + c 3 / 2
l = 1 c 1 2 1 c 2 1 / 2 1 3 ( 2 + c ) 1 c 1 + c 1 4 ( 3 + c ) 1 c 1 + c 3 / 2
l = 2 1 2 ( 3 c 2 1 ) 1 2 c 1 c 2 1 / 2 1 4 ( 1 c 2 ) 1 20 ( 8 + 9 c + 3 c 2 ) 1 c 1 + c 3 / 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Koleski, G.; Bickel, T. Stokes Equation in a Semi-Infinite Region: Generalization of the Lamb Solution and Applications to Marangoni Flows. Fluids 2020, 5, 249. https://doi.org/10.3390/fluids5040249

AMA Style

Koleski G, Bickel T. Stokes Equation in a Semi-Infinite Region: Generalization of the Lamb Solution and Applications to Marangoni Flows. Fluids. 2020; 5(4):249. https://doi.org/10.3390/fluids5040249

Chicago/Turabian Style

Koleski, Goce, and Thomas Bickel. 2020. "Stokes Equation in a Semi-Infinite Region: Generalization of the Lamb Solution and Applications to Marangoni Flows" Fluids 5, no. 4: 249. https://doi.org/10.3390/fluids5040249

APA Style

Koleski, G., & Bickel, T. (2020). Stokes Equation in a Semi-Infinite Region: Generalization of the Lamb Solution and Applications to Marangoni Flows. Fluids, 5(4), 249. https://doi.org/10.3390/fluids5040249

Article Metrics

Back to TopTop