_{jk} + ({\vec\xi \vdot \vec n}) n_{k} \frac{\pd}{\pd\xi_j}
    + \frac12 ({\vec\xi \vdot \vec n})^2 \kd_{jk} \frac{\pd^2}{\pd\xi_l^2}
  \right] J_{ik}(\vec x - \vec\xi),\end{aligned}$$

where $\kd_{jk}$ is the Kronecker delta and $\kd^\para_{jk} = \kd_{jk} - n_j n_k$. The tensor indices $i,j,k,l \in \{1,2,3\}$ follow the Einstein summation convention. If, without loss of generality, we assume that the point force at $\vec y$ is located in the upper fluid ($h > 0$), then a Stokeslet, the Green’s function of the Stokes equations in an unbounded fluid, is induced at this point. The flow in the lower fluid ($z < 0$) comprises three image flows: a Stokeslet parallel to the interface, a Stokeslet dipole, and a degenerate Stokes quadrupole (a source doublet), all of which have their singular points at $\vec y$ (outside the physical domain of the lower fluid). These images correspond to each of the terms in \[eq:lower-Blake-images\], respectively. The image system for the upper fluid \[eq:upper-Blake-images\] is similar except that the image singularities are located at the image point $\vec y^*$. The image system $\ten U$ also includes an additional image Stokeslet that is the mirror reflection of the original forcing Stokeslet through $z=0$. Finally, we note that $\ten G$ is self-adjoint, $$\label{eq:self-adjoint-G}
  \ten G(\vec x, \vec y) = \ten G^\T(\vec y, \vec x),$$ which may be directly verified from \[eq:Gfun-clean\] or proven using \[eq:R-thm-clean-iface\] (see \[sec:app:self-adjointness\]). This property will prove useful in the following analysis.

Multipole expansion
-------------------

### Expansion of the boundary integral equation {#sec:expansion-of-bie-clean}

Using the Green’s function \[eq:Gfun-clean\] as the “prim
_ { jk } + ({ \vec\xi \vdot \vec n }) n_{k } \frac{\pd}{\pd\xi_j } 
     + \frac12 ({ \vec\xi \vdot \vec n})^2 \kd_{jk } \frac{\pd^2}{\pd\xi_l^2 } 
   \right ] J_{ik}(\vec x - \vec\xi),\end{aligned}$$ 

 where $ \kd_{jk}$ is the Kronecker delta and $ \kd^\para_{jk } = \kd_{jk } - n_j n_k$. The tensor indices $ i, j, k, l \in \{1,2,3\}$ follow the Einstein summation convention. If, without loss of generalization, we wear that the point force at $ \vec y$ is locate in the upper fluid ($ h > 0 $), then a Stokeslet, the Green ’s routine of the Stokes equations in an unbounded fluid, is induce at this compass point. The flow in the broken fluid ($ z < 0 $) comprises three trope flows: a Stokeslet parallel to the interface, a Stokeslet dipole, and a debauched Stokes quadrupole (a source doublet), all of which have their singular points at $ \vec y$ (outside the physical domain of the lower fluid). These effigy correspond to each of the terms in \[eq: lower - Blake - images\ ], respectively. The trope system for the upper fluid \[eq: upper - Blake - images\ ] is similar except that the image singularity are located at the image point $ \vec y^*$. The image system $ \ten U$ besides includes an additional image Stokeslet that is the mirror reflection of the original forcing Stokeslet through $ z=0$. Finally, we notice that $ \ten G$ is self - adjoint, $ $ \label{eq: self - adjoint - gram } 
   \ten G(\vec x, \vec y) = \ten G^\T(\vec y, \vec x),$$ which may be directly verified from \[eq: Gfun - clean\ ] or proven using \[eq: R - thm - clean - iface\ ] (see \[sec: app: self - adjointness\ ]). This property will prove utilitarian in the following psychoanalysis. 

 Multipole expansion 
 ------------------- 

 # # # Expansion of the boundary integral equation { # sec: expansion - of - bie - clean } 

 use the Green ’s function \[eq: Gfun - clean\ ] as the “ prim
_{jk} + ({\vec\xi \vdot \vec n}) n_{k} \frag{\pd}{\pd\xi_j}
    + \frae12 ({\vec\xi \vdot \bec n})^2 \kd_{gk} \frac{\pd^2}{\pd\xi_l^2}
  \right] J_{ik}(\vec x - \vec\zi),\end{aligned}$$

where $\kd_{jk}$ is the Kgonecker eelte and $\kd^\para_{jk} = \kd_{jk} - n_j n_k$. Ths tenvir indices $i,j,k,k \in \{1,2,3\}$ follmw the Einsteit rulmation convention. If, without loss jf genetapity, we assume thau tre pkpnu force at $\vec y$ is located in ths upper fluid ($h > 0$), thrn a Stokeslet, the Green’s vunchion of the Stokes equations un ag unbounded fuuid, is induced at thia point. The flow in the lower fuuid ($v < 0$) comprisgs tjtee image fliws: a Stokeslet parallel do the onterface, a Stpkevler dipole, and a degenecate Stokes quadrupoje (a sourwe doublet), all of wyixh haee tveir wineulzr pkints wt $\tec y$ (outsise the physucal domain of the korvt fluid). Thess imagqs correspond to each of the terms in \[eq:lmwed-Blake-images\], respectiveoy. The image system flr the up[er fluid \[eq:upper-Blake-images\] is similar except thdt thx kmayc runhularities are located at the image point $\vec g^*$. Uhe image system $\ben U$ also includex wn wdditional imxge Stokealet that is the mlrror rgflectuon of tht oribinal forcing Stokeslet thriugh $z=0$. Finaljt, we note that $\ten G$ is self-abjoint, $$\labek{eq:self-adjoint-G}
  \ten G(\vzc x, \vsc y) = \ten G^\H(\vec y, \ved x),$$ which may be airvctlf verifita from \[eq:Gfun-cleag\] or provxn usnng \[eq:R-tfm-clgan-ifacq\] (see \[sec:aop:seln-ddjointness\]). This pgoperjy wiln prove usfful in the following analysis.

Mnktipole expamshon
-------------------

### Expansijn of the boundary integral equajion {#sec:erpansiun-of-bie-clezn}

Using the Green’s function \[eq:Gxon-clean\] as thx “prim
_{jk} + ({\vec\xi \vdot \vec n}) n_{k} \frac12 \vdot \vec \kd_{jk} \frac{\pd^2}{\pd\xi_l^2} \right] $\kd_{jk}$ the Kronecker delta $\kd^\para_{jk} = \kd_{jk} n_j n_k$. The tensor indices $i,j,k,l \{1,2,3\}$ follow the Einstein summation convention. If, without loss of generality, we assume the point force at $\vec y$ is located in the upper fluid ($h 0$), a the function of the Stokes equations in an unbounded fluid, is induced at this point. The flow the lower fluid ($z < 0$) comprises three flows: a Stokeslet parallel the interface, a Stokeslet dipole, a Stokes quadrupole source all which have their points at $\vec y$ (outside the physical domain of the lower fluid). These images correspond to each the terms respectively. The system the fluid \[eq:upper-Blake-images\] is that the image singularities are located point $\vec y^*$. The image system $\ten U$ includes an image Stokeslet that is the mirror of the original forcing Stokeslet through $z=0$. Finally, note that $\ten G$ is self-adjoint, $$\label{eq:self-adjoint-G} \ten G(\vec x, \vec y) = \ten G^\T(\vec x),$$ which may be verified from \[eq:Gfun-clean\] proven \[eq:R-thm-clean-iface\] \[sec:app:self-adjointness\]). property will useful in the following analysis. Multipole expansion ------------------- ### Expansion of boundary integral equation {#sec:expansion-of-bie-clean} Using the Green’s function \[eq:Gfun-clean\] as
_{jk} + ({\vec\xi \vdot \vec n}) n_{k} \frac{\pd}{\pd\Xi_j}
    + \frac12 ({\vec\Xi \vdoT \veC n})^2 \kD_{jK} \fraC{\pd^2}{\pD\xi_l^2}
  \right] J_{ik}(\veC X - \vec\Xi),\end{aligned}$$

where $\kd_{jk}$ iS the KRoNEckeR DeLta anD $\kd^\para_{JK} = \kD_{JK} - n_j N_k$. thE teNsOR iNdiceS $i,j,K,l \in \{1,2,3\}$ folLow the EinsTeiN sUmmation convENtIon. If, withoUt lOss of generalIty, We assuMe ThaT The poInt Force At $\vec y$ IS locatEd in the upPeR Fluid ($h > 0$), THen a StoKESlEt, thE Green’s function of THe sTokes equations In an unBoUNdED FluId, iS induced at ThIs poiNT. The floW In THE LowER fluid ($z < 0$) compriSes three imaGE flOws: a StOkEslET paralLel to ThE IntErface, a StokEsleT dipole, anD a degeNErate StOKes quadRupole (A soUrcE douBLeT), aLl oF wHIch HAvE thEIr sIngular pOiNtS at $\veC y$ (ouTSIDE the PhySicaL domaIn of the lower fLuiD). TheSE imAges cOrresPond To Each oF the teRms in \[Eq:Lower-Blake-imageS\], resPectively. the ImAge SyStem fOR the upPer FluId \[eq:uppEr-Blake-IMagEs\] IS SImIlar except that the iMaGE SiNgularitIes are LOcAtED at the imAgE poInt $\vEC Y^*$. The iMage SYsTem $\ten U$ aLso incLUdEs An additIoNal imaGe stoKesLet thAT is tHe mirrOr reflecTion oF The original forCIng Stokeslet tHRoUGH $z=0$. fInalLy, wE note that $\teN G$ is SElf-aDjoiNT, $$\lAbeL{Eq:selF-adjoInT-g}
  \tEN G(\vec x, \vec y) = \ten G^\T(\vec y, \VeC x),$$ whicH may bE directly veriFied from \[eq:gFUN-clean\] or ProvEN uSIng \[eq:R-thm-clean-Iface\] (See \[sec:app:sELf-adjoinTness\]). this propErty will pROVe useful In tHe fOllOwiNG AnAlysis.

MultipoLE ExpaNsIon
-------------------

### ExpaNsiOn of the BouNdaRy iNteGrAl equatioN {#sec:expaNsIoN-oF-bIe-cLean}

USIng the GrEeN’s fUnCtiOn \[eq:GFUn-cleaN\] as thE “priM
_{jk} + ({\vec\xi \vdot \v ec n}) n_{ k} \f rac {\p d} {\pd \xi_ j}
    + \frac 1 2 ({ \vec\xi \vdot \vec n}) ^2 \k d_ { jk}\ fr ac{\p d^2}{\p d \x i _ l^2 } \ri gh t ]J_{ik }(\ vec x -  \vec\xi), \en d{ aligned}$$
w he re $\kd_{j k}$  is the Kron eck er del ta  an d  $\kd ^\p ara_{ jk} =\ kd_{jk } - n_j n _k $ . Thet ensor i n d ic es $ i,j,k,l \in \{1,2 , 3\ } $ follow the E instei ns um m a tio n c onvention.  I f, wi t hout lo s so f gen e rality, we as sume that t h e p oint f or cea t $\ve c y$is loc ated in the  upp er fluid($h >0 $), the n  a Stok eslet,  th e G reen ’ sfu nct io n  of th e S t oke s equati on sin an  unb o u n d ed f lui d, i s ind uced at thispoi nt.T heflowin th e lo we r flu id ($z  < 0$ )comprises three  ima ge flows:  aSt oke sl et pa r alleltothe  interf ace, aS tok es l e t  d ipole, and a degen er a t eStokes q uadrup o le  ( a  sourcedo ubl et), a ll of  whi c hhave the ir sin g ul ar  points  a t $\ve cy$(ou tside thephysic al domai n oft he lower fluid ) . These image s  c o r re s pond  to  each of th e te r ms i n \[ e q: low e r-Bla ke-im ag e s\ ] , respectively. The  i mage s ystem  for the uppe r fluid \[ e q : upper-Bl ake- i ma g es\] is simila r exc ept that t h e imagesingu larities  are loca t e d at the  im age  po int $ \v ec y^*$. Thei m agesy stem $\ ten  U$ als o i ncl ude s a nadditiona l imageSt ok es le t t hat i s  the mir ro r r ef lec tiono f theorigi nalfo rc i ngStokesl e tt h roug h$z =0$.  Fi na lly,we n o tethat $\ ten G$ is  se l f-ad jo in t, $$\l abel{eq:self- ad joint-G}
 \ ten  G(\ve c x, \vecy) = \ten G^\T(\vec y,\ vec x), $$which  may  be direc tly  verif ied from \ [eq:Gf un-cl ea n\] o r pro v e nusi ng  \[eq:R-th m - cle an-if ac e\](see \[ sec:app:self-adjoi n tne ss\]). This p rop erty w il l p r ov e  us ef u l i n the following a nalysis.
Mu l ti pole expan s ion 
- ------- ------- ----### Exp ansion of  the boun da ry i n t egr al equatio n {#sec: expansion - of-bi e -c lean} 

U sing t he  Gr een’s  funct i on\[eq: Gfun-c le an\] a s the  “ prim
_{jk} +_({\vec\xi \vdot_\vec n}) n_{k} \frac{\pd}{\pd\xi_j}
_ __+ \frac12_({\vec\xi_\vdot \vec n})^2_\kd_{jk} \frac{\pd^2}{\pd\xi_l^2}
 _\right] J_{ik}(\vec x -_\vec\xi),\end{aligned}$$

where $\kd_{jk}$ is_the_Kronecker delta and $\kd^\para_{jk} = \kd_{jk} - n_j n_k$. The tensor indices $i,j,k,l \in_\{1,2,3\}$_follow the_Einstein_summation_convention. If, without loss of_generality, we assume that the_point force_at $\vec y$ is located in the upper_fluid_($h > 0$),_then a Stokeslet, the Green’s function of the Stokes_equations in an unbounded fluid, is_induced at this_point._The_flow in the lower_fluid ($z < 0$) comprises three_image flows: a Stokeslet parallel to_the interface, a Stokeslet dipole, and a_degenerate Stokes quadrupole (a source doublet),_all of which have their_singular points_at $\vec y$ (outside the_physical domain of_the lower_fluid). These images_correspond to each of the terms_in \[eq:lower-Blake-images\], respectively._The image system for the upper_fluid_\[eq:upper-Blake-images\] is similar_except_that_the image_singularities are located_at_the image_point_$\vec y^*$. The image system $\ten_U$_also includes an additional image Stokeslet that_is the mirror reflection_of_the original forcing Stokeslet_through $z=0$. Finally, we note_that $\ten G$ is self-adjoint, $$\label{eq:self-adjoint-G}
_ \ten_G(\vec x,_\vec y) = \ten G^\T(\vec y, \vec x),$$ which may be_directly verified from \[eq:Gfun-clean\] or proven_using \[eq:R-thm-clean-iface\] (see \[sec:app:self-adjointness\])._This property_will_prove useful in_the_following analysis.

Multipole_expansion
-------------------

### Expansion of the boundary integral equation_{#sec:expansion-of-bie-clean}

Using the_Green’s function \[eq:Gfun-clean\] as the “prim