---
abstract: 'We investigate the motion of a thin rigid body in Stokes flow and the corresponding slender body approximation used to model sedimenting fibers. In particular, we derive a rigorous error bound comparing the rigid slender body approximation to the classical PDE for rigid motion in the case of a closed loop with constant radius. Our main tool is the slender body PDE framework established by the authors and D. Spirn in [@closed_loop; @free_ends], which we adapt to the rigid setting.'
author:
- |
    Yoichiro Mori, Laurel Ohm [^1]\
    *School of Mathematics, University of Minnesota, Minneapolis, MN 55455*
bibliography:
- 'rigid\_bib.bib'
title: 'An error bound for the slender body approximation of a thin, rigid fiber sedimenting in Stokes flow'
---

Introduction
============

Determining the motion of a three-dimensional rigid body sedimenting in a Stokesian fluid is an important problem in both theoretical and computational fluid mechanics. This motion is described by a classical PDE [@corona2017integral; @galdi1999steady; @weinberger1972variational], which we write below in the case of a thin rigid body. We use ${\mathcal{E}}({\bm{u}}) = \frac{1}{2}(\nabla{\bm{u}}+(\nabla{\bm{u}})^{\rm T})$ to denote the symmetric gradient, and $\bm{\sigma}=\bm{\sigma}({\bm{u}},p) = 2{\mathcal{E}}({\bm{u}})-p{\bf I}$ to denote the stress tensor. Let $\Sigma_\epsilon$ denote a closed loop slender body of radius $\epsilon>0$ (to be made precise in Section \[geometry\]) and let $\Omega_\epsilon ={\mathbb{R}}^3 \setminus \overline{\Sigma_\epsilon}$ and $\Gamma_\epsilon={\partial}\Sigma_\epsilon$ (see Figure \[fig:coord\_sys\]). The full PDE description of a slender body undergoing a rigid motion in Stokes flow may be written as follows: $$\label{rigid}
\begin{aligned}
-\Delta {\bm{u}}^{\rm r} +\nabla p^{\rm r} &=0   \hspace{2.6cm} \text{ in } \Omega_\epsilon \\
{{\rm{div}\,}}{\bm{u}}^{\rm r} &= 0  \hspace{2.6cm} \text{ in } \Omega_\epsilon \\
{\bm{u}}^{\rm r}({\bm{x}}) &= {\bm{v}}^{\rm r} + \bm{\omega}^{\rm r}\times {\bm{x}}, \qquad {\bm{x}}\in \Gamma_\epsilon \\
{\bm{u}}^{\rm r}({\bm{x}}) &\to 0  \hspace{2.6cm} \text{as }{\left\lvert {\bm{x}}\right\rvert}\to \infty 
\end{aligned}$$ and $$\begin{aligned}
\int_{\Gamma_\epsilon} \bm{\sigma}^{\rm r}\bm{n} \; dS &= \bm{F}, \quad \int_{\Gamma_\epsilon} {\bm{x}}\times (\bm{\sigma}^{\rm r}\bm{n}) \; dS = \bm{T} .\end{aligned}$$ Here the total force $\bm{F}\in {\mathbb{R}}^3$ and torque $\bm{T}\in {\mathbb{R}}^3$ are given and we aim to solve for the linear velocity ${\bm{v}}^{\rm r}\in {\mathbb{R}}^3$ and angular velocity $\bm{\omega}^{\rm r}\in{\mathbb{R}}^3$ of the body. Note that the boundary value problem is in fact valid for rigid bodies of arbitrary shape, but for the purposes of this paper we specifically consider here a slender closed loop. Using the variational framework of [@galdi1999steady; @gonzalez2004dynamics; @weinberger1972variational], it can be shown that is a well-posed PDE.\
On the computational side, there has been much recent interest in numerical simulations of rigid particle sedimentation [@guazzelli2006sedimentation; @guazzelli2011fluctuations], and various tools have been developed to facilitate these simulations [@corona2017integral; @jung2006periodic; @mitchell2015sedimentation].\

![The geometry of the rigid fiber may be parameterized with respect to the orthogonal frame ${\bm{e}}_t(s)$, ${\bm{e}}_{n_1}(s)$, ${\bm{e}}_{n_2}(s)$ defined in Section \[geometry\].[]{data-label="fig:coord_sys"}](SB_geometry_rigid "fig:")\

For a thin rigid body, a commonly-used tool for simplifying simulations is slender body theory, which exploits the thin geometry of the body by approximating the filament as a one-dimensional force density distributed along the fiber centerline. Slender body theory is a popular method for modeling sedimentation of thin fibers, both rigid [@butler2002dynamic; @park2010cloud; @saintillan2005smooth; @shin2009structure] and semi-flexible [@li2013sedimentation; @manikantan2014instability]. Here we will specifically consider the slender body theory established by Keller and Rubinow [@keller1976slender] and further developed in [@gotz2000interactions; @johnson1980improved; @tornberg2004simulating].\
Let ${\bm{X}}: {\mathbb{T}}\equiv {\mathbb{R}}/ {\mathbb{Z}}\to {\mathbb{R}}^3$ denote the coordinates of the slender body centerline, parameterized by arclength $s$ and defined more precisely in Section \[geometry\]. Given a line force density $\bm{f}^{\rm s}(s)$, $s\in {\mathbb{T}}$, the slender body approximation yields a direct expression approximating the velocity of the fiber, given by [@shelley2000stokesian]: $$\label{SBT_expr}
\begin{aligned}
{\bm{u}}^{\rm s}_{\rm C}(s) &= \bm{\Lambda}[\bm{f}^{\rm s}](s) + \bm{K}[\bm{f}^{\rm s}](s),\\
\bm{\Lambda}[\bm{f}](s) &:=  \frac{1}{8\pi}\big[({\bf I}- 3{\bm{e}}_t{\bm{e}}_t^{\rm T})-2({\bf I}+{\bm{e}}_t{\bm{e}}_t^{\rm T}) \log(\pi\epsilon/4) \big]{\bm f}(s) \\
\bm{K}[\bm{f}](s) &:= \frac{1}{8\pi}\int_{{\mathbb{T}}} \left[ \left(\frac{{\bf I}}{|\bm{R}_0|}+ \frac{\bm{R}_0\bm{R}_0^{\rm T}}{|\bm{R}_0|^3}\right){\bm f}(s') - \frac{{\bf I}+{\bm{e}}_t(s){\bm{e}}_t(s)^{\rm T} }{|\sin (\pi(s-s'))/\pi|} {\bm f}(s)\right] \, ds'.
\end{aligned}$$ Here ${\bm{e}}_t(s)$ is the unit tangent vector to ${\bm{X}}(s)$ and $\bm{R}_0(s,s') = {\bm{X}}(s) - {\bm{X}}(s')$. The slender body approximation generally allows for bending and flexing of the filament along its centerline and requires specifying the one-dimensional force density over the length of the fiber centerline. If the fiber is constrained to be fully rigid, only the total force $\bm{F}$ and torque $\bm{T}$ must be specified, where $$\label{SBT_cond1}
 \int_{{\mathbb{T}}} \bm{f}^{\rm s}(s) \, ds = \bm{F}, \quad \int_{{\mathbb{T}}} {\bm{X}}(s)\times \bm{f}^{\rm s}(s) \, ds = \bm{T}.$$ Additionally, we constrain the motion of the fiber centerline to be rigid, i.e. $$\label{SBT_cond2}
{\bm{u}}^{\rm s}_{\rm C}(s) = {\bm{v}}^{\rm s}+\bm{\omega}^{\rm s}\times {\bm{X}}(s)$$ for constant vectors ${\bm{v}}^{\rm s}$, $\bm{\omega}^{\rm s}$. These constraints give rise to a system of integral equations which must be solved to obtain the line force density along the slender body (see [@gustavsson2009gravity; @tornberg2006numerical]).\
The aim of this paper is to establish a rigorous error bound between the slender body approximation for rigid motion - and the classical PDE describing the sedimentation of a rigid fiber immersed in Stokes flow. We show the following theorem.

\[rigid\_theorem\] Let $\Sigma_\epsilon$ be a slender body as defined in Section \[geometry\]. Suppose the total force $\bm{F}\in {\mathbb{R}}^3$ and torque $\bm{T}\in {\mathbb{R}}^3$ are given, and assume that rigid slender body approximation - is satisfied by some $\bm{f}^{\rm s}\in C^1({\mathbb{T}})$. Then the difference ${\bm{v}}^{\rm r}-{\bm{v}}^{\rm s}$, $\bm{\omega}^{\rm r}-\bm{\omega}^{\rm s}$ between the linear and angular velocities of true rigid motion and the slender body approximation - satisfies $$\label{rigid_err}
{\left\lvert {\bm{v}}^{\rm r}-{\bm{v}}^{\rm s} \right\rvert} + {\left\lvert \bm{\omega}^{\rm r}-\bm{\omega}^{\rm s} \right\rvert} \le C\sqrt{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{3/2} \big({\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert \bm{f}^{\rm s} \right\rVert}_{C^1({\mathbb{T}})} + {\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert} \big)$$ for $C$ depending on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.

The constants $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$ have to do only with the shape of the fiber centerline and are defined in Section \[geometry\]. Notice that we must assume that we can find $\bm{f}^{\rm s}\in C^1({\mathbb{T}})$ satisfying -, and that this $\bm{f}^{\rm s}$ then appears in the final error bound. In this sense, our bound should be considered as an [*a posteriori*]{} error estimate, similar to the type of estimates commonly used in finite element analysis [@ainsworth2011posteriori]. To obtain an [*a priori*]{} bound, we would need to able to say that such an $\bm{f}^{\rm s}$ is then bounded by the given $\bm{F}$ and $\bm{T}$, but the proof of such a bound is hindered by the likely ill-posedness of the rigid slender body approximation itself. Before discussing this issue in greater detail, we note how it arises in the error analysis.\
In order to compare the classical PDE with the rigid slender body approximation, we introduce an intermediary PDE which we will call the [*slender body PDE for rigid motion*]{}. The idea follows from the notion of slender body PDE proposed by the authors and D. Spirn in [@closed_loop] and [@free_ends] as a framework for analyzing the error introduced by the Keller-Rubinow slender body approximation for closed-loop and open-ended fibers, respectively. To construct the rigid slender body PDE, we impose that the velocity of the slender body is uniform over each cross section $s$ of the fiber. In particular, we approximate ${\bm{x}}\in \Gamma_\epsilon$ as its $L^2$ projection onto the fiber centerline ${\bm{X}}(s)$, thereby ignoring slight differences in torque across the slender body. Note that the slender body geometry is defined in Section \[geometry\] such that this projection onto the fiber centerline is unique; i.e. the notion of “fiber cross section" is well-defined. We define the slender body PDE for rigid motion as follows: $$\label{SB_PDE}
\begin{aligned}
-\Delta {\bm{u}}^{\rm p} +\nabla p^{\rm p} &=0  \hspace{3.27cm} \text{in } \Omega_\epsilon  \\
 {{\rm{div}\,}}{\bm{u}}^{\rm p} &=0  \hspace{3.27cm} \text{in } \Omega_\epsilon  \\
{\bm{u}}^{\rm p}({\bm{x}}) &= {\bm{v}}^{\rm p} + \bm{\omega}^{\rm p}\times {\bm{X}}(s) \qquad \text{on } \Gamma_\epsilon  \\
{\bm{u}}^{\rm p}({\bm{x}}) &\to 0 \hspace{3.23cm} \text{as }{\left\lvert {\bm{x}}\right\rvert}\to \infty
\end{aligned}$$ and $$\begin{aligned}
\int_{\Gamma_\epsilon} \bm{\sigma}^{\rm p}\bm{n} \; {\mathcal{J}}_\epsilon(s,\theta) \, d\theta \, ds &= \bm{F}, \qquad
\int_{{\mathbb{T}}} {\bm{X}}(s)\times \bigg(\int_0^{2\pi} \bm{\sigma}^{\rm p}\bm{n} \, {\mathcal{J}}_\epsilon(s,\theta) d\theta \bigg) ds = \bm{T}.\end{aligned}$$ Here we have written $dS = {\mathcal{J}}_\epsilon(s,\theta)\, d\theta\, ds$, where ${\mathcal{J}}_\epsilon$ is the Jacobian factor on the slender body surface, which we parameterize as a tube about ${\bm{X}}(s)$ using surface angle $\theta$ (see Section \[geometry\] and expression ). We show that for a closed filament, the rigid slender body PDE is in fact close to the classical PDE for rigid motion [@galdi1999steady; @weinberger1972variational] – in particular, the variation in torque over any cross section of the slender body is higher order in $\epsilon$.\
In the case of a flexible filament with a prescribed force density per unit length along the centerline, the slender body PDE of [@closed_loop; @free_ends] is well-posed, and the difference between the slender body approximation and the PDE solution can be estimated in terms of the slender body radius and the given line force density. We would like to use the existing error analysis in [@closed_loop] to bound the difference between the rigid slender body approximation and the rigid slender body PDE solution. The rigid case is complicated by the fact that the existing error bound relies on knowledge of line force density along the filament, and therefore on the solvability of the integral equation for this force density arising from -.\
This raises the issue of the lack of a general solution theory for this particular type of integral equation. Specifically, a detailed spectral analysis by Götz [@gotz2000interactions] in the case of a straight slender body centerline shows that the slender body operator $(\bm{\Lambda}+\bm{K})$ in is not necessarily invertible for all values of $\epsilon$ and all centerline velocities ${\bm{u}}^{\rm s}_{\rm C}$. A similar result for a perfectly circular, planar centerline was shown by Shelley-Ueda in [@shelley2000stokesian]. For fibers with more general centerline curvature, a spectral analysis of the slender body integral operator is complicated and the invertibility properties remain unclear.\
In rigid slender body theory, this same integral operator arises in the integral equations for the line force density $\bm{f}^{\rm s}$. Therefore, to carry out an error analysis for the rigid slender body problem we must assume from the start that we are considering a slender body approximation that gives rise to $\bm{f}^{\rm s}\in C^1({\mathbb{T}})$ satisfying -. This $\bm{f}^{\rm s}$ must be included in the final error bound along with the total force $\bm{F}$ and torque $\bm{T}$, giving rise to a type of [*a posteriori*]{} error estimate. Furthermore, the $\sqrt{\epsilon}$ bound is the best we can do using our current techniques, but it is possible that a more complete solution theory for the slender body approximation could improve this bound.\
In practice, various regularizations of the slender body integral operator are used to combat this ill-posedness, although rigorous justification for these regularizations is still needed. Given that the slender body PDE framework of [@closed_loop; @free_ends] [*is*]{} well-posed, we should be able to use this framework to come up with the best regularization for the ill-posed slender body approximation. However, this is truly a deeper issue that we plan to explore in future work.

Outline of the proof of Theorem \[rigid\_theorem\]
--------------------------------------------------

The strategy for proving Theorem \[rigid\_theorem\] is to show that, given $\bm{F}$ and $\bm{T}$, the solution to the rigid slender body PDE is close to both the classical rigid PDE solution and the rigid slender body approximation - .\
First, we must show that the rigid slender body PDE is well-posed. Using Definition \[rigid\_weak\_p\] of a weak solution to the rigid slender body PDE , we show the following.

\[SB\_PDE\_well\] Let $\Sigma_\epsilon$ be a slender body as defined in Section \[geometry\]. Given $\bm{F}$ and $\bm{T}\in{\mathbb{R}}^3$, there exists a unique weak solution $({\bm{u}}^{\rm p},p^{\rm p})\in {\mathcal{R}}_\epsilon^{{{\rm{div}\,}}} \times L^2(\Omega_\epsilon)$ to the slender body PDE for rigid motion satisfying the estimate $$\label{totalPest}
{\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} + {\left\lVert p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2} ( {\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert} )$$ for $C$ depending on $c_\Gamma$ and $\kappa_{\max}$.

Theorem \[SB\_PDE\_well\] can be established using many of the same tools from the well-posedness theory in [@closed_loop]. In addition, we will make use of the following bound along the slender body centerline ${\bm{X}}(s)$:

\[v\_omega\_bd\] Let ${\bm{X}}$ be as in Section \[geometry\] and consider constant vectors ${\bm{v}}$, $\bm{\omega} \in {\mathbb{R}}^3$. Then $${\left\lvert {\bm{v}}\right\rvert}+{\left\lvert \bm{\omega} \right\rvert} \le C{\left\lVert {\bm{v}}+ \bm{\omega}\times {\bm{X}}\right\rVert}_{L^2({\mathbb{T}})}$$ for $C$ depending only on $c_\Gamma$ and $\kappa_{\max}$.

We will first prove Lemma \[v\_omega\_bd\] in Section \[2ndLemma\]; then Theorem \[SB\_PDE\_well\] quickly follows using some of the key inequalities collected in Section \[ineq\].\
With the variational framework for , comparing to is relatively straightforward. Using Lemma \[v\_omega\_bd\], we show that the difference between the true rigid motion and the slender body PDE description satisfies the following lemma.

\[true\_vs\_SB\] Let ${\bm{X}}$ be as in Section \[geometry\]. Given $\bm{F}$ and $\bm{T}\in {\mathbb{R}}^3$, let $({\bm{v}}^{\rm r},\bm{\omega}^{\rm r})$ be the corresponding boundary values satisfying and let $({\bm{v}}^{\rm p},\bm{\omega}^{\rm p})$ be the boundary values satisfying . Then $$\label{true_vs_SB_bd}
{\left\lvert \bm{\omega}^{\rm r}-\bm{\omega}^{\rm p} \right\rvert}+ {\left\lvert {\bm{v}}^{\rm r}-{\bm{v}}^{\rm p} \right\rvert} \le \epsilon{\left\lvert \log\epsilon \right\rvert}C( {\left\lvert \bm{T} \right\rvert}+ {\left\lvert \bm{F} \right\rvert})$$ where $C$ depends only on $c_\Gamma$ and $\kappa_{\max}$.

The main difficulties in proving Theorem \[rigid\_theorem\] arise in comparing to - . As discussed, the solvability issue for the rigid slender body approximation is addressed by assuming that we are only considering a rigid slender body approximation - that gives rise to a force density $\bm{f}^{\rm s}\in C^1({\mathbb{T}})$. An additional difficulty arises in that in order to use the error analysis framework of [@closed_loop], the line force density along the slender body must be the same for both the slender body approximation and the slender body PDE. Therefore we need to define yet another intermediary PDE.\
Given $\bm{f}^{\rm s}\in C^1({\mathbb{T}})$ satisfying - for given $\bm{F}$ and $\bm{T}\in {\mathbb{R}}^3$, we define ${\bm{u}}^{\rm p,s}$ as the solution to the PDE: $$\label{SB_PDE_ps}
\begin{aligned}
-\Delta {\bm{u}}^{\rm p,s} +\nabla p^{\rm p,s} &=0  \hspace{2.75cm} \text{in } \Omega_\epsilon  \\
 {{\rm{div}\,}}{\bm{u}}^{\rm p,s} &=0 \hspace{2.75cm} \text{in } \Omega_\epsilon  \\
\int_0^{2\pi} (\bm{\sigma}^{\rm p,s}\bm{n}) \; {\mathcal{J}}_\epsilon(s,\theta) \, d\theta &= \bm{f}^{\rm s}(s)  \hspace{2.05cm} \text{on }\Gamma_\epsilon \\
{\bm{u}}^{\rm p,s}\big|_{\Gamma_\epsilon} &= {\rm Tr}({\bm{u}}^{\rm p,s})(s), \qquad \text{ unknown but independent of } \theta \\
{\bm{u}}^{\rm p,s}&\to 0  \hspace{2.65cm} \text{as } {\left\lvert {\bm{x}}\right\rvert}\to\infty.
\end{aligned}$$ Here ${\rm Tr}({\bm{u}}^{\rm p,s})(s)$ denotes the trace of ${\bm{u}}^{\rm p,s}$ on $\Gamma_\epsilon$, and $\theta$ refers to the parameterization of $\Gamma_\epsilon$ as a tube about ${\bm{X}}(s)$ (see Section \[geometry\]). By [@closed_loop], we know that a (weak) solution $({\bm{u}}^{\rm p,s},p^{\rm p,s})$ exists and is unique. Now, ${\rm Tr}({\bm{u}}^{\rm p,s})(s)$ may not be precisely a rigid motion, but we can show that it is close. In particular, by Theorem 1.3 in [@closed_loop], we may bound the difference between ${\rm Tr}({\bm{u}}^{\rm p,s})(s)$ and ${\bm{u}}^{\rm s}_{\rm C}(s) =  {\bm{v}}^{\rm s}+\bm{\omega}^{\rm s}\times {\bm{X}}(s)$ by $$\label{centerline_err} 
{\left\lVert {\rm Tr}( {\bm{u}}^{\rm p,s}) - ({\bm{v}}^{\rm s}+\bm{\omega}^{\rm s} \times {\bm{X}})  \right\rVert}_{L^2({\mathbb{T}})} \le C\epsilon {\left\lvert \log\epsilon \right\rvert}^{3/2}{\left\lVert \bm{f}^{\rm s} \right\rVert}_{C^1({\mathbb{T}})}$$ for $C$ depending only on $c_\Gamma$ and $\kappa_{\max}$.\
A further technical issue arises in comparing to . In order to obtain a useful estimate of the difference between $({\bm{u}}^{\rm p,s}-{\bm{u}}^{\rm p},p^{\rm p,s}-p^{\rm p})$ in terms of only $\bm{F}$, $\bm{T}$, and $\bm{f}^{\rm s}(s)$, we will need a careful characterization of the $\epsilon$-dependence in a higher regularity estimate for solutions to (see Lemma \[high\_reg\]). Note that for a (sufficiently smooth) sedimenting rigid body, once well-posedness of the PDE has been established, higher regularity of the solution follows by standard arguments for a Stokes Dirichlet boundary value problem. In our case, the novelty is determining how the higher regularity bound scales with $\epsilon$. Our proof (see Appendix \[reg\_lem\]) makes use of the local coordinate system valid near the slender body. We obtain commutator estimates for the tangential derivatives along the slender body surface and use an integration by parts argument, along with the form of the Stokes equations in local coordinates, to show that the bound for an additional derivative of the rigid slender body PDE solution scales like $1/\epsilon$, up to logarithmic corrections.\
Now, using the variational framework for the slender body PDE along with this higher regularity lemma, we can show the following estimate.

\[ups\_up\_err\] Let ${\bm{u}}^{\rm p,s}$ satisfy and let $({\bm{v}}^{\rm p},\bm{\omega}^{\rm p})$ denote the rigid slender body PDE boundary values satisfying . Then $${\left\lVert {\rm Tr}( {\bm{u}}^{\rm p,s}) - ( {\bm{v}}^{\rm p} + \bm{\omega}^{\rm p}\times {\bm{X}}) \right\rVert}_{L^2({\mathbb{T}})} \le C\sqrt{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{3/2} \big({\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert \bm{f}^{\rm s} \right\rVert}_{C^1({\mathbb{T}})} + {\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert} \big)$$ for $C$ depending on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.

Combining estimate with Lemma \[ups\_up\_err\] and using Lemma \[v\_omega\_bd\] with ${\bm{v}}^{\rm p}-{\bm{v}}^s$ and $\bm{\omega}^{\rm p}-\bm{\omega}^{\rm s}$ in place of ${\bm{v}}$ and $\bm{\omega}$, we obtain the following bound for the difference between the slender body approximation - and the slender body PDE : $$\label{PDE_vs_SBT}
{\left\lvert {\bm{v}}^{\rm p}-{\bm{v}}^{\rm s} \right\rvert} + {\left\lvert \bm{\omega}^{\rm p}-\bm{\omega}^{\rm s} \right\rvert} \le C\sqrt{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{3/2} \big({\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert \bm{f}^{\rm s} \right\rVert}_{C^1({\mathbb{T}})} + {\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert} \big).$$

Finally, combining the estimate with Lemma \[true\_vs\_SB\] yields Theorem \[rigid\_err\]. The remainder of this paper is thus devoted to showing Lemmas \[v\_omega\_bd\] - \[ups\_up\_err\]. We will begin by introducing the variational framework for and noting some key inequalities in Section \[variational0\]. In Section \[2ndLemma\], we show Lemma \[v\_omega\_bd\] and use it to derive estimates for $({\bm{u}}^{\rm p},p^{\rm p},{\bm{v}}^{\rm p},\bm{\omega}^{\rm p})$ satisfying . These estimates can then be used to show Theorem \[SB\_PDE\_well\]. In Section \[3rdLemma\], we use the variational framework for the rigid slender body PDE to prove Lemma \[true\_vs\_SB\]. Finally, in Section \[1stLemma\], we prove Lemma \[ups\_up\_err\] to complete the proof of Theorem \[rigid\_err\].

Geometry and variational framework {#variational0}
==================================

We begin in Section \[geometry\] with a precise definition of the slender body geometry. In Section \[variational\], we introduce the variational form of the slender body PDE for rigid motion , which, along with the variational form of , will provide the framework for obtaining Theorem \[rigid\_theorem\]. Finally, in Section \[ineq\], we make note of some key inequalities that will be used throughout the remainder of this paper.

Slender body geometry {#geometry}
---------------------

As in [@closed_loop], we let ${\bm{X}}: {\mathbb{T}}\equiv {\mathbb{R}}/ {\mathbb{Z}}\to {\mathbb{R}}^3$ denote the coordinates of a closed, non-self-intersecting $C^3$ curve in ${\mathbb{R}}^3$, parameterized by arclength $s$. We require that $$\label{cgamma}
 \inf_{s\neq s'}\frac{{\left\lvert {\bm{X}}(s)-{\bm{X}}(s') \right\rvert}}{{\left\lvert s-s' \right\rvert}} \ge c_\Gamma$$ for some constant $c_\Gamma>0$.\
Along ${\bm{X}}(s)$ we consider the orthonormal frame $({\bm{e}}_t(s),{\bm{e}}_{n_1}(s),{\bm{e}}_{n_2}(s))$ defined in [@closed_loop]. Here ${\bm{e}}_t(s) = \frac{d{\bm{X}}}{ds}$ is the unit tangent vector to ${\bm{X}}(s)$ and $({\bm{e}}_{n_1}(s),{\bm{e}}_{n_2}(s))$ span the plane normal to ${\bm{e}}_t(s)$. The frame satisfies the ODEs $$\label{C2_frame}
\frac{d}{ds}{\bm{e}}_t = \kappa_1{\bm{e}}_{n_1}+\kappa_2{\bm{e}}_{n_2}, \quad \frac{d}{ds}{\bm{e}}_{n_1} = -\kappa_1{\bm{e}}_t+\kappa_3{\bm{e}}_{n_2}, \quad  \frac{d}{ds}{\bm{e}}_{n_2} = -\kappa_2{\bm{e}}_t-\kappa_3{\bm{e}}_{n_1}$$ where $\kappa_1^2(s) +\kappa_2^2(s) = \kappa^2(s)$, the fiber curvature, and $\kappa_3$ is a constant satisfying ${\left\lvert \kappa_3 \right\rvert}\le \pi$. We require the orthonormal frame to be $C^2$ and denote $$\kappa_{\max} := \max_{s\in{\mathbb{T}}} {\left\lvert \kappa(s) \right\rvert}, \quad \xi_{\max} = \max_{s\in{\mathbb{T}}} {\left\lvert \frac{{\partial}^3{\bm{X}}}{{\partial}s^3} \right\rvert} .$$ Note that ${\left\lvert {\partial}\kappa_1 /{\partial}s \right\rvert} +{\left\lvert {\partial}\kappa_2/{\partial}s \right\rvert} \le \xi_{\max}+2(\kappa_{\max}+\pi)$.\
We define $${\bm{e}}_\rho(s,\theta) := \cos\theta{\bm{e}}_{n_1}(s)+\sin\theta{\bm{e}}_{n_2}(s)$$ and, for some $r_{\max}=r_{\max}(c_\Gamma,\kappa_{\max})\le \frac{1}{2\kappa_{\max}}$, we can uniquely parameterize points ${\bm{x}}$ within a neighborhood ${\rm dist}({\bm{x}},{\bm{X}})< r_{\max}$ of ${\bm{X}}(s)$ as $${\bm{x}}= {\bm{X}}(s) + \rho{\bm{e}}_\rho(s,\theta), \quad 0\le \rho <r_{\max}.$$

For $\epsilon<r_{\max}/4$, we may then define a slender body of uniform radius $\epsilon$ as $$\label{SB_def} 
\Sigma_\epsilon:= \big\{{\bm{x}}\in {\mathbb{R}}^3 \; : \; {\bm{x}}= {\bm{X}}(s) + \rho{\bm{e}}_\rho(s,\theta), \; \rho <\epsilon \big\}.$$ We parameterize the slender body surface $\Gamma_\epsilon = {\partial}\Sigma_\epsilon$ as $$\Gamma_\epsilon =  {\bm{X}}(s) + \epsilon{\bm{e}}_\rho(s,\theta).$$ In addition, we may parameterize the Jacobian factor ${\mathcal{J}}_\epsilon(s,\theta)$ on the slender body surface as $$\label{jac_fac}
{\mathcal{J}}_\epsilon(s,\theta) = \epsilon\big(1-\epsilon(\kappa_1(s)\cos\theta+\kappa_2(s)\sin\theta) \big).$$

Variational form of  {#variational}
--------------------

Letting $\Omega_\epsilon={\mathbb{R}}^3\backslash\overline{\Sigma_\epsilon}$ for $\Sigma_\epsilon$ as in Section \[geometry\], we recall the following function spaces, used in [@closed_loop] to study a slender body PDE of the form . We use $D^{1,2}(\Omega_\epsilon)$ to denote the homogeneous Sobolev space $$\label{D12}
D^{1,2}(\Omega_\epsilon) = \big\{ {\bm{u}}\in L^6(\Omega_\epsilon) \; : \; \nabla {\bm{u}}\in L^2(\Omega_\epsilon) \big\},$$ which, due to the Sobolev inequality in $\Omega_\epsilon \subset {\mathbb{R}}^3$ (see Lemma \[sobolev\]), is a Hilbert space with norm ${\left\lVert \nabla {\bm{u}}\right\rVert}_{L^2(\Omega_\epsilon)}$. We define $D^{1,2}_0(\Omega_\epsilon)$ as the closure of $C_0^\infty(\Omega_\epsilon)$ (smooth, compactly supported test functions) in $D^{1,2}(\Omega_\epsilon)$.\
We also recall the space ${\mathcal{A}}_\epsilon$, the subspace of $D^{1,2}(\Omega_\epsilon)$ with $\theta$-independent boundary values: $$\label{Aepsilon}
{\mathcal{A}}_\epsilon = \big\{ {\bm{u}}\in D^{1,2}(\Omega_\epsilon) \; : \; {\bm{u}}\big|_{\Gamma_\epsilon} = {\bm{u}}(s) \big\}.$$ Here the boundary value ${\bm{u}}\big|_{\Gamma_\epsilon} = {\bm{u}}(s)$ is not directly specified but is required to be independent of the surface angle $\theta$. We define ${\mathcal{A}}_\epsilon^{{\rm{div}\,}}$ to be the divergence-free subspace of ${\mathcal{A}}_\epsilon$.\
We also recall the variational form of , examined in detail in [@closed_loop].

A weak solution ${\bm{u}}^{\rm p,s}\in {\mathcal{A}}_\epsilon^{{{\rm{div}\,}}}$ to satisfies $$\label{weak_no_p}
\int_{\Omega_\epsilon} 2 {\mathcal{E}}({\bm{u}}^{\rm p,s}): {\mathcal{E}}({\bm{v}}) \, d{\bm{x}}= \int_{\mathbb{T}}{\bm{v}}(s)\cdot\bm{f}^{\rm s}(s) \, ds$$ for any ${\bm{v}}\in{\mathcal{A}}_\epsilon^{{{\rm{div}\,}}}$. In addition, for ${\bm{u}}^{\rm p,s}$ satisfying , there exists a unique pressure $p^{\rm p,s}\in L^2(\Omega_\epsilon)$ satisfying $$\label{weak_with_p}
\int_{\Omega_\epsilon} \big( 2 {\mathcal{E}}({\bm{u}}^{\rm p,s}): {\mathcal{E}}({\bm{v}}) - p^{\rm p,s}\, {{\rm{div}\,}}{\bm{v}}\big)\, d{\bm{x}}= \int_{\mathbb{T}}{\bm{v}}(s)\cdot\bm{f}^{\rm s}(s) \, ds$$ for any ${\bm{v}}\in {\mathcal{A}}_\epsilon$.

To study , we define the following subspace of ${\mathcal{A}}_\epsilon$, where we further restrict the boundary value to be a rigid motion: $$\label{Repsilon}
{\mathcal{R}}_\epsilon =  \big\{ {\bm{u}}\in D^{1,2}(\Omega_\epsilon) \; : \; {\bm{u}}\big|_{\Gamma_\epsilon} = {\bm{v}}+ \bm{\omega}\times{\bm{X}}(s) \text{ for }{\bm{v}},\, \bm{\omega}\in {\mathbb{R}}^3  \big\}.$$ Again, ${\bm{v}}$ and $\bm{\omega}$ are not directly specified but are required to be constant vectors in ${\mathbb{R}}^3$. We let ${\mathcal{R}}_\epsilon^{{\rm{div}\,}}$ denote the divergence-free subspace of ${\mathcal{R}}_\epsilon$.\
We then define a weak solution to the rigid motion slender body PDE as follows.

\[rigid\_weak\] A weak solution ${\bm{u}}^{\rm p}\in {\mathcal{R}}_\epsilon^{{\rm{div}\,}}$ to satisfies $$\int_{\Omega_\epsilon} 2\, {\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}(\bm{\varphi}) \, d{\bm{x}}= {\bm{v}}_\varphi \cdot\bm{F} + \bm{\omega}_\varphi \cdot\bm{T}$$ for any $\bm{\varphi}\in {\mathcal{R}}_\epsilon^{{\rm{div}\,}}$, where we denote $\bm{\varphi}\big|_{\Gamma_\epsilon} =  {\bm{v}}_\varphi + \bm{\omega}_\varphi \times {\bm{X}}(s)$.

Given the existence and uniqueness of ${\bm{u}}^{\rm p}$ satisfying Definition \[rigid\_weak\], using an essentially identical proof to that in Section 2.2 of [@closed_loop], we can establish an equivalent notion of weak solution that includes a corresponding weak pressure $p^{\rm p}\in L^2(\Omega_\epsilon)$ and removes the divergence-free restriction on test functions $\bm{\varphi}$.

\[rigid\_weak\_p\] Given ${\bm{u}}^{\rm p}\in {\mathcal{R}}_\epsilon^{{\rm{div}\,}}$ satisfying Definition \[rigid\_weak\], there exists a unique $p^{\rm p}\in L^2(\Omega_\epsilon)$ satisfying $$\int_{\Omega_\epsilon} \big(2\, {\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}(\bm{\varphi}) - p\, {{\rm{div}\,}}\bm{\varphi} \big) \, d{\bm{x}}= {\bm{v}}_\varphi \cdot\bm{F} + \bm{\omega}_\varphi\cdot \bm{T}$$ for any $\bm{\varphi}\in {\mathcal{R}}_\epsilon$. Here we again denote $\bm{\varphi}\big|_{\Gamma_\epsilon} =  {\bm{v}}_\varphi + \bm{\omega}_\varphi \times {\bm{X}}(s)$.

Important inequalities {#ineq}
----------------------

In addition to the definitions of Section \[variational\], we collect the statements of various inequalities that are used throughout the paper, keeping track of the $\epsilon$-dependence in any constants that arise. The proofs of these inequalities are mostly contained in [@closed_loop], with the exception of Lemma \[trace2\], which appears in Appendix \[appendix\].\
First, we note the following ${\mathcal{A}}_\epsilon$ trace inequality, which holds for functions ${\bm{u}}\in {\mathcal{A}}_\epsilon$ due to $\theta$-independence on $\Gamma_\epsilon$. As a slight abuse of notation, the trace operator ${\rm Tr}$, when applied to ${\mathcal{A}}_\epsilon$ functions, will be considered as both a function on $\Gamma_\epsilon$ and on ${\mathbb{T}}$.

*($L^2({\mathbb{T}})$ trace inequality)*\[trace1\] Let $\Omega_{\epsilon}={\mathbb{R}}^3 \backslash \overline{\Sigma_{\epsilon}}$ be as in Section \[geometry\]. Then any ${\bm{u}}\in {\mathcal{A}}_\epsilon$ satisfies $$\|{\rm Tr}({\bm{u}})\|_{L^2({\mathbb{T}})} \le C |\log\epsilon|^{1/2} \| \nabla {\bm{u}}\|_{L^2(\Omega_{\epsilon})},$$ where the constant $C$ depends on $\kappa_{\max}$ and $c_{\Gamma}$ but is independent of $\epsilon$.

The proof of this lemma appears in Appendix A.2.1 of [@closed_loop].\
On the other hand, for general $D^{1,2}(\Omega)$ functions, the following trace inequality holds over the surface $\Gamma_\epsilon$:

*($L^2(\Gamma_\epsilon)$ trace inequality)*\[trace2\] Let $\Omega_{\epsilon}={\mathbb{R}}^3 \backslash \overline{\Sigma_{\epsilon}}$ be as in Section \[geometry\]. Then any ${\bm{u}}\in D^{1,2}(\Omega_\epsilon)$ satisfies $$\|{\rm Tr}({\bm{u}})\|_{L^2(\Gamma_\epsilon)} \le C \sqrt{\epsilon}|\log\epsilon|^{1/2} \| \nabla {\bm{u}}\|_{L^2(\Omega_{\epsilon})},$$ where the constant $C$ depends on $\kappa_{\max}$ and $c_{\Gamma}$ but is independent of $\epsilon$.

The proof of Lemma \[trace2\] appears in Appendix \[appendix\].\
We will also need the following Korn inequality.

*(Korn inequality)*\[korn\] Let $\Omega_{\epsilon}={\mathbb{R}}^3 \backslash \overline{\Sigma_{\epsilon}}$ be as in Section \[geometry\]. Then any ${\bm{u}}\in D^{1,2}(\Omega_{\epsilon})$ satisfies $$\|\nabla {\bm{u}}\|_{L^2(\Omega_{\epsilon})} \le C\|{\mathcal{E}}({\bm{u}})\|_{L^2(\Omega_{\epsilon})},$$ where the constant $C$ depends only on $\kappa_{\max}$ and $c_{\Gamma}$.

The proof of $\epsilon$-independence in the Korn constant is given in Appendix A.2.2 - A.2.3 in [@closed_loop].\
Finally, we make use of the following pressure estimate.

\[pressure\] For $({\bm{u}},p)$ satisfying the Stokes equations in $\Omega_\epsilon$, we have $${\left\lVert p \right\rVert}_{L^2(\Omega_\epsilon)} \le C {\left\lVert {\mathcal{E}}({\bm{u}}) \right\rVert}_{L^2(\Omega_\epsilon)}$$ for $C$ independent of $\epsilon$.

The proof of this lemma exactly follows the proof of estimate (2.17) in [@closed_loop].

Proof of Lemma \[v\_omega\_bd\] and a corollary {#2ndLemma}
===============================================

Here we prove Lemma \[v\_omega\_bd\] and make note of a corollary which allows us to obtain a useful bound for functions in $\rm{R}_\epsilon$. This corollary, along with the Korn inequality (Lemma \[korn\]) and pressure estimate (Lemma \[pressure\]), then allows us to prove Theorem \[SB\_PDE\_well\].

Note that Lemma \[v\_omega\_bd\] is obviously true when ${\bm{v}}=\bm{\omega}=0$; thus we can assume that at least one of ${\bm{v}},\bm{\omega}$ is nonzero. Suppose that Lemma \[v\_omega\_bd\] does not hold. Then we may choose a sequence of triples $({\bm{v}}_k,\bm{\omega}_k,{\bm{X}}_k(s))$ such that the following properties hold for each $k=1,2,3,\dots$. First, ${\bm{v}}_k,\bm{\omega}_k\in{\mathbb{R}}^3$ satisfy ${\left\lvert {\bm{v}}_k \right\rvert}^2+{\left\lvert \bm{\omega}_k \right\rvert}^2=1$, and ${\bm{X}}_k(s)$ is a closed curve satisfying the geometric constraints of Section \[geometry\] – in particular, ${\left\lvert {\bm{X}}_k'' \right\rvert}\le \kappa_{\max}$. For simplicity, we also take $\int_{{\mathbb{T}}}{\bm{X}}_k(s) \, ds=0$. In addition, $$1 = {\left\lvert {\bm{v}}_k \right\rvert}^2+{\left\lvert \bm{\omega}_k \right\rvert}^2 > k^2\int_{{\mathbb{T}}}{\left\lvert {\bm{v}}_k+\bm{\omega}_k\times {\bm{X}}_k(s) \right\rvert}^2 \, ds.$$ Then $$\int_{{\mathbb{T}}}{\left\lvert {\bm{v}}_k+\bm{\omega}_k\times {\bm{X}}_k(s) \right\rvert}^2 \, ds < \frac{1}{k^2} \to 0$$ as $k\to\infty$. Since ${\bm{v}}_k,\bm{\omega}_k$ are just vectors in ${\mathbb{R}}^3$, some limit ${\bm{v}}_\infty,\bm{\omega}_\infty$ exists. Furthermore, since each ${\bm{X}}_k$ is controlled in $C^2$ by $\kappa_{\max}$, we have that (passing to a subsequence) ${\bm{X}}_k\to {\bm{X}}_\infty$ in $C^1$ for some closed, unit length curve ${\bm{X}}_\infty(s)$. Thus $$\int_{{\mathbb{T}}}{\left\lvert {\bm{v}}_\infty +\bm{\omega}_\infty \times {\bm{X}}_\infty(s) \right\rvert}^2 \, ds =0,$$ and therefore $\bm{\omega}_\infty\times {\bm{X}}_\infty(s)\equiv -{\bm{v}}_\infty$. But $\bm{\omega}_\infty$ and ${\bm{v}}_\infty$ are both constant vectors with ${\left\lvert {\bm{v}}_\infty \right\rvert}^2+{\left\lvert \bm{\omega}_\infty \right\rvert}^2=1$, while ${\bm{X}}_\infty(s)$ necessarily has nonzero curvature. Thus $\bm{\omega}_\infty\times {\bm{X}}_\infty(s)$ cannot identically equal the constant vector $-{\bm{v}}_\infty$. Furthermore, because ${\bm{X}}_k$ was allowed to vary among curves satisfying the constraints of Section \[geometry\], the constant $C$ arising in Lemma \[v\_omega\_bd\] depends only on $c_\Gamma$ and $\kappa_{\max}$.

Given $\bm{F}$ and $\bm{T}$, as an immediate corollary to Lemma \[v\_omega\_bd\] we obtain the following useful bound for any function $\bm{\varphi}\in{\mathcal{R}}_\epsilon$.

\[omegaP\] Consider $\bm{\varphi}\in {\mathcal{R}}_\epsilon$ with boundary value denoted by ${\bm{v}}_\varphi+\bm{\omega}_\varphi\times {\bm{X}}(s)$. Then $$\label{curved_est}
|{\bm{v}}_\varphi|+|\bm{\omega}_\varphi| \le C |\log\epsilon|^{1/2} {\left\lVert {\mathcal{E}}(\bm{\varphi}) \right\rVert}_{L^2(\Omega_\epsilon)}$$ where $C$ depends only on $c_\Gamma$ and $\kappa_{\max}$.

Using Lemma \[v\_omega\_bd\] along with the slender body trace estimate (Lemma \[trace1\]) and Korn inequality (Lemma \[korn\]), we immediately obtain $$\begin{aligned}
{\left\lvert {\bm{v}}_\varphi \right\rvert}+ {\left\lvert \bm{\omega}_\varphi \right\rvert} &\le {\left\lVert {\bm{v}}_\varphi+\bm{\omega}_\varphi\times{\bm{X}}\right\rVert}_{L^2({\mathbb{T}})} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert \nabla \bm{\varphi} \right\rVert}_{L^2(\Omega_\epsilon)} \\
&\le C{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert {\mathcal{E}}(\bm{\varphi}) \right\rVert}_{L^2(\Omega_\epsilon)}\end{aligned}$$ for $C$ depending only on $c_\Gamma$ and $\kappa_{\max}$.

Using Corollary \[omegaP\] and the variational formulation of , we may now prove Theorem \[SB\_PDE\_well\].

We first show the existence of a weak solution ${\bm{u}}^{\rm p}\in {\mathcal{R}}_\epsilon^{{{\rm{div}\,}}}$ satisfying Definition \[rigid\_weak\]. Note that the bilinear form appearing on the left hand side of Definition \[rigid\_weak\] is bounded on ${\mathcal{R}}_\epsilon^{{{\rm{div}\,}}}$, as $${\left\lvert \int_{\Omega_\epsilon} 2\, {\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}(\bm{\varphi}) \, d{\bm{x}}\right\rvert} \le 2{\left\lVert  {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2(\Omega_\epsilon)}{\left\lVert {\mathcal{E}}(\bm{\varphi}) \right\rVert}_{L^2(\Omega_\epsilon)} \le 2{\left\lVert  \nabla{\bm{u}}^{\rm p}  \right\rVert}_{L^2(\Omega_\epsilon)}{\left\lVert \nabla\bm{\varphi} \right\rVert}_{L^2(\Omega_\epsilon)}.$$ Coercivity of the bilinear form also follows by the Korn inequality (Lemma \[korn\]). Furthermore, using Corollary \[omegaP\], the linear functional on the right hand side of Definition \[rigid\_weak\] is bounded for $\bm{\varphi}\in {\mathcal{R}}_\epsilon^{{{\rm{div}\,}}}$, as $${\left\lvert  {\bm{v}}_\varphi \cdot\bm{F} + \bm{\omega}_\varphi\cdot \bm{T} \right\rvert} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert {\mathcal{E}}(\bm{\varphi}) \right\rVert}_{L^2(\Omega_\epsilon)}({\left\lvert \bm{F} \right\rvert} +{\left\lvert \bm{T} \right\rvert}) \le C{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert \nabla\bm{\varphi} \right\rVert}_{L^2(\Omega_\epsilon)}({\left\lvert \bm{F} \right\rvert} +{\left\lvert \bm{T} \right\rvert}).$$ Then, by the Lax-Milgram theorem, there exists a unique weak solution ${\bm{u}}^{\rm p}\in {\mathcal{R}}_\epsilon^{{{\rm{div}\,}}}$ to .\
In addition, using the variational form of along with Corollary \[omegaP\], we have that ${\bm{u}}^{\rm p}$ satisfies $$\begin{aligned}
\int_{\Omega_\epsilon} {\left\lvert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rvert}^2 \, d{\bm{x}}&= \int_{\Gamma_\epsilon} ({\bm{v}}^{\rm p}+\bm{\omega}^{\rm p}\times {\bm{X}}(s))\cdot(\bm{\sigma}^{\rm p}\bm{n}) \, dS \\
&= {\bm{v}}^{\rm p}\cdot\int_{\Gamma_\epsilon} \bm{\sigma}^{\rm p}\bm{n} \, dS + \bm{\omega}^{\rm p}\cdot \int_{{\mathbb{T}}} {\bm{X}}(s)\times \bigg(\int_0^{2\pi} (\bm{\sigma}^{\rm p}\bm{n}) \, {\mathcal{J}}_\epsilon(s,\theta) d \theta \bigg) ds \\
&\le {\left\lvert {\bm{v}}^{\rm p} \right\rvert}{\left\lvert \bm{F} \right\rvert} + {\left\lvert \bm{\omega}^{\rm p} \right\rvert}{\left\lvert \bm{T} \right\rvert} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2(\Omega_\epsilon)}({\left\lvert \bm{F} \right\rvert} +{\left\lvert \bm{T} \right\rvert}) \\
&\le \frac{1}{2}{\left\lVert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2(\Omega_\epsilon)}^2 + C{\left\lvert \log\epsilon \right\rvert}({\left\lvert \bm{F} \right\rvert}^2 +{\left\lvert \bm{T} \right\rvert}^2),\end{aligned}$$ where we have used Young’s inequality in the last line. We thus obtain the estimate $$\label{Eup_est}
{\left\lVert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2(\Omega_\epsilon)} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2} ( {\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert} ).$$

As noted after Definition \[rigid\_weak\], the existence of a unique corresponding weak pressure $p^{\rm p}\in L^2(\Omega_\epsilon)$ satisfying Definition \[rigid\_weak\_p\] as well as Lemma \[pressure\] follows by an essentially identical proof to that appearing in Section 2.2 of [@closed_loop].\
Combining with Lemmas \[korn\] and \[pressure\] then yields the bound .

Classical versus slender body PDE description of rigid motion {#3rdLemma}
=============================================================

Using the variational framework of Section \[variational0\] along with Lemma \[v\_omega\_bd\], we prove Lemma \[true\_vs\_SB\] comparing the classical PDE and slender body PDE descriptions of rigid slender body motion.

The difference $\overline{\bm{u}}= {\bm{u}}^{\rm r}- {\bm{u}}^{\rm p}$, $\overline p= p^{\rm r}-  p^{\rm p}$, $\overline{\bm{\sigma}}=\bm{\sigma}^{\rm r}- \bm{\sigma}^{\rm p}$, $\overline{\bm{\omega}}= \bm{\omega}^{\rm r}- \bm{\omega}^{\rm p}$, $\overline{\bm{v}}= {\bm{v}}^{\rm r}- {\bm{v}}^{\rm p}$ satisfies the PDE $$\label{diff_PDE}
\begin{aligned}
-\Delta \overline{\bm{u}}+\nabla \overline p &=0, \quad {{\rm{div}\,}}\overline{\bm{u}}=0 \hspace{2.3cm} \text{in } \Omega_\epsilon \\
\overline{\bm{u}}({\bm{x}}) &= \overline{\bm{v}}+ \overline{\bm{\omega}}\times {\bm{x}}+ \epsilon \bm{\omega}^{\rm p} \times {\bm{e}}_\rho, \qquad {\bm{x}}\in \Gamma_\epsilon \\
\overline{\bm{u}}({\bm{x}}) &\to 0 \hspace{4.3cm} \text{as }{\left\lvert {\bm{x}}\right\rvert}\to \infty \\
\int_{\Gamma_\epsilon} \overline{\bm{\sigma}}\bm{n} \, dS &=0, \quad \int_{\Gamma_\epsilon} {\bm{x}}\times (\overline{\bm{\sigma}}\bm{n}) \, dS = - \epsilon \int_{\Gamma_\epsilon} {\bm{e}}_\rho \times (\bm{\sigma}^{\rm p}\bm{n}) \; dS.
\end{aligned}$$

Then, (formally) multiplying by $\overline{\bm{u}}$ and integrating by parts, we have that $\overline{\bm{u}}$ satisfies $$\label{diff_bound0}
\begin{aligned}
\int_{\Omega_\epsilon} 2|{\mathcal{E}}(\overline{\bm{u}})|^2 \, d{\bm{x}}&= \int_{\Gamma_\epsilon} \big( \overline{\bm{v}}+ \overline{\bm{\omega}}\times {\bm{x}}+ \epsilon\bm{\omega}^{\rm p}\times{\bm{e}}_\rho \big) \cdot(\overline{\bm{\sigma}}\bm{n}) \, dS \\
&= \overline{\bm{v}}\cdot\int_{\Gamma_\epsilon}\overline{\bm{\sigma}}\bm{n} \, dS + \overline{\bm{\omega}}\cdot\int_{\Gamma_\epsilon} {\bm{x}}\times(\overline{\bm{\sigma}}\bm{n}) \, dS + \epsilon\bm{\omega}^{\rm p}\cdot\int_{\Gamma_\epsilon} {\bm{e}}_\rho\times(\overline{\bm{\sigma}}\bm{n}) \, dS \\
&= -\epsilon \overline{\bm{\omega}}\cdot\int_{\Gamma_\epsilon} {\bm{e}}_\rho \times (\bm{\sigma}^{\rm p}\bm{n}) \, dS + \epsilon \bm{\omega}^{\rm p}\cdot\int_{\Gamma_\epsilon} {\bm{e}}_\rho \times(\overline{\bm{\sigma}}\bm{n}) \, dS.
\end{aligned}$$

To estimate the right hand side of , we first need to define a smooth cutoff function $\phi(\rho)$ satisfying $$\label{cutoff_def}
\phi(\rho) = \begin{cases}
1, & \rho < 2 \\
0, & \rho >4
\end{cases}$$ with smooth decay satisfying $$\label{phi_decay}
{\left\lvert \frac{d\phi}{d\rho} \right\rvert} \le c_\phi.$$ Then for ${\bm{x}}={\bm{X}}(s)+\rho{\bm{e}}_\rho(\theta,s)$ in a neighborhood of $\Gamma_\epsilon$, we define $\phi_\epsilon(\rho):= \phi(\rho/\epsilon)$.\
We estimate the second term on the right hand side first, noting that the estimation of the first term will be essentially identical. Using index notation (the subscript $\cdot_{,j}$ signifies $\frac{{\partial}\cdot}{{\partial}x_j}$; sum over repeated indices) along with the divergence theorem, we may write $$\label{cross_est0}
\begin{aligned} 
\int_{\Gamma_\epsilon} \big({\bm{e}}_\rho \times(\overline{\bm{\sigma}}\bm{n}) \big)_i \, dS &=  \int_{\Gamma_\epsilon} \varepsilon_{ijk} (e_\rho)_j \overline\sigma_{k\ell} n_\ell \, dS = \int_{\Omega_\epsilon} (\phi \, \varepsilon_{ijk}(e_\rho)_j \overline\sigma_{k\ell})_{,\ell} \, d{\bm{x}}\\
&= \int_{\Omega_\epsilon} \varepsilon_{ijk}\big( \phi_{,\ell}(e_\rho)_j \overline\sigma_{k\ell} + \phi \,(e_\rho)_{j,\ell} \overline\sigma_{k\ell} \big)\, d{\bm{x}}.
\end{aligned}$$ Here $\varepsilon_{ijk}$ is the alternating symbol $$\varepsilon_{ijk} = \begin{cases}
1, & \text{ for even permutations of }i,j,k \\
-1, & \text{ for odd permutations of }i,j,k \\
0, & \text{ if } i=j,j=k, \text{ or }k=i, \end{cases}$$ and we have used that $\overline{\bm{\sigma}}$ is divergence-free.

Although $\overline{\bm{\sigma}}$ only belongs to $L^2(\Omega_\epsilon)$ [*a priori*]{}, we may justify equation because ${{\rm{div}\,}}\overline{\bm{\sigma}}=0\in L^2(\Omega_\epsilon)$, and therefore $\overline{\bm{\sigma}}$ belongs to the space $H_{{{\rm{div}\,}}}(\Omega_\epsilon):= \{\bm{\varphi}\in L^2(\Omega_\epsilon) \, : \, {{\rm{div}\,}}\bm{\varphi}\in L^2(\Omega_\epsilon) \}$. Then we may use the divergence theorem for $H_{{{\rm{div}\,}}}$ functions given in Lemma IV.3.3 of [@boyer2012mathematical] to obtain .

Now, due to the cutoff $\phi_\epsilon$, the integrand on the right hand side of is supported only within the region $${\mathcal{O}}_\epsilon := \big\{{\bm{X}}(s)+\rho {\bm{e}}_\rho(s,\theta) \; : \; s\in{\mathbb{T}}, \; \epsilon \le \rho\le4\epsilon, \; 0\le \theta<2\pi \big\}$$ with ${\left\lvert {\mathcal{O}}_\epsilon \right\rvert}= C\epsilon^2$ for some $C$ depending only on $c_\Gamma$ and $\kappa_{\max}$.\
Within ${\mathcal{O}}_\epsilon$, defining ${\widehat{\kappa}}(s,\theta):= \kappa_1(s)\cos\theta+\kappa_2(s)\sin\theta$, we have $$\label{grad_rho}
\begin{aligned}
{\left\lvert \nabla {\bm{e}}_\rho(s,\theta) \right\rvert} &=  {\left\lvert \frac{1}{\rho}\frac{{\partial}{\bm{e}}_\rho}{{\partial}\theta}{\bm{e}}_\theta^{\rm T}+ \frac{1}{1-\rho{\widehat{\kappa}}}\bigg(\frac{{\partial}{\bm{e}}_\rho}{{\partial}s} - \kappa_3\frac{{\partial}{\bm{e}}_\rho}{{\partial}\theta} \bigg){\bm{e}}_t^{\rm T}  \right\rvert} \\
&= {\left\lvert \frac{1}{\rho}{\bm{e}}_\theta{\bm{e}}_\theta^{\rm T} - \frac{1}{1-\rho{\widehat{\kappa}}}(\kappa_1\cos\theta+\kappa_2\sin\theta) {\bm{e}}_t{\bm{e}}_t^{\rm T}  \right\rvert} \le \frac{1}{\epsilon}+ 4\kappa_{\max},
 \end{aligned}$$ where the final $\kappa_{\max}$ bound is shown in Appendix \[reg\_lem\].\
Using , , and Cauchy-Schwarz, we may estimate as $$\label{cross_bound}
\begin{aligned}
{\left\lvert \int_{\Gamma_\epsilon} {\bm{e}}_\rho \times(\overline{\bm{\sigma}}\bm{n}) \, dS \right\rvert} &\le \int_{{\mathcal{O}}_\epsilon} \big({\left\lvert \phi_\epsilon\nabla{\bm{e}}_\rho \right\rvert}+{\left\lvert \nabla\phi_\epsilon \right\rvert}){\left\lvert \overline{\bm{\sigma}} \right\rvert} \, d{\bm{x}}\\
&\le {\left\lvert {\mathcal{O}}_\epsilon \right\rvert}^{1/2} \bigg(\frac{1}{\epsilon} + 4\kappa_{\max} + \frac{c_\phi}{\epsilon} \bigg)\bigg(\int_{\Omega_\epsilon} {\left\lvert \overline{\bm{\sigma}} \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2} \\
& \le C \bigg(\int_{\Omega_\epsilon} {\left\lvert \overline{\bm{\sigma}} \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2}
\end{aligned}$$ where $C$ depends only on the shape of ${\bm{X}}$ – in particular, $c_\Gamma$ and $\kappa_{\max}$. Finally, using Lemma \[pressure\], we obtain $$\label{cross_bound2}
{\left\lvert \int_{\Gamma_\epsilon} {\bm{e}}_\rho \times(\overline{\bm{\sigma}}\bm{n}) \, dS \right\rvert} \le C \bigg(\int_{\Omega_\epsilon} \big({\left\lvert {\mathcal{E}}(\overline{\bm{u}}) \right\rvert}^2 + \overline p^2 \big) \, d{\bm{x}}\bigg)^{1/2} \le C \bigg(\int_{\Omega_\epsilon} {\left\lvert {\mathcal{E}}(\overline{\bm{u}}) \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2}.$$

Following exactly the same procedure, we can also show $$\label{cross_bound3}
{\left\lvert \int_{\Gamma_\epsilon} {\bm{e}}_\rho \times(\bm{\sigma}^{\rm p}\bm{n}) \, dS \right\rvert} \le C \bigg(\int_{\Omega_\epsilon} {\left\lvert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2}.$$

Furthermore, in the same way as in Lemma \[omegaP\], it can be shown that $$\label{omegav_bound}
{\left\lvert \overline{\bm{\omega}} \right\rvert}+ {\left\lvert \overline{\bm{v}}\right\rvert} \le C{\left\lVert \overline{\bm{v}}+\overline{\bm{\omega}}\times{\bm{x}}+ \epsilon\bm{\omega}^{\rm p}\times {\bm{e}}_\rho \right\rVert}_{L^2(\Gamma_\epsilon)} \le C\sqrt{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert {\mathcal{E}}(\overline{\bm{u}}) \right\rVert}_{L^2(\Omega_\epsilon)},$$ where we have used the $L^2(\Gamma_\epsilon)$ trace estimate (Lemma \[trace2\]).\
Then, using and in along with Lemma \[omegaP\] and , we have $$\label{diff_bound1}
\begin{aligned}
\int_{\Omega_\epsilon} 2|{\mathcal{E}}(\overline{\bm{u}})|^2 \, d{\bm{x}}&\le \epsilon C{\left\lvert \overline{\bm{\omega}} \right\rvert}\bigg(\int_{\Omega_\epsilon}{\left\lvert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2} + \epsilon C{\left\lvert \bm{\omega}^{\rm p} \right\rvert} \bigg(\int_{\Omega_\epsilon}{\left\lvert {\mathcal{E}}(\overline{\bm{u}}) \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2} \\
&\le \epsilon{\left\lvert \log\epsilon \right\rvert}^{1/2}C\bigg(\int_{\Omega_\epsilon}{\left\lvert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2}\bigg(\int_{\Omega_\epsilon}{\left\lvert {\mathcal{E}}(\overline{\bm{u}}) \right\rvert}^2 \, d{\bm{x}}\bigg)^{1/2} \\
&\le \epsilon^2{\left\lvert \log\epsilon \right\rvert}C\int_{\Omega_\epsilon}{\left\lvert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rvert}^2 \, d{\bm{x}}+ \int_{\Omega_\epsilon}{\left\lvert {\mathcal{E}}(\overline{\bm{u}}) \right\rvert}^2 \, d{\bm{x}},
\end{aligned}$$ where we have used Young’s inequality in the last line. Then, using , we obtain $${\left\lVert {\mathcal{E}}(\overline{\bm{u}}) \right\rVert}_{L^2(\Omega_\epsilon)} \le \epsilon{\left\lvert \log\epsilon \right\rvert}C( {\left\lvert \bm{T} \right\rvert}+ {\left\lvert \bm{F} \right\rvert}).$$

Finally, using again, we obtain Lemma \[true\_vs\_SB\].

Proof of Lemma \[ups\_up\_err\]  {#1stLemma}
================================

Finally, we prove Lemma \[ups\_up\_err\] comparing the rigid slender body PDE to the intermediary slender body PDE .\
We begin by defining $$\label{fp_def}
\bm{f}^{\rm p}(s):=\int_0^{2\pi} (\bm{\sigma}^{\rm p}\bm{n}) \, {\mathcal{J}}_\epsilon(s,\theta) d\theta$$ for $\bm{\sigma}^{\rm p}$ as in , and establish the following:

\[fp\_FT\] Suppose the slender body $\Sigma_\epsilon$ is as in Section \[geometry\] – in particular, ${\bm{X}}\in C^3({\mathbb{T}})$. Let the total force $\bm{F}$ and torque $\bm{T}$ be given, and let $\bm{f}^{\rm p}$ be as defined in . Then $${\left\lVert \bm{f}^{\rm p} \right\rVert}_{L^2({\mathbb{T}})} \le C{\left\lvert \log\epsilon \right\rvert}^{3/2}({\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert})$$ for $C$ depending only on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.

The proof of this lemma relies on a higher regularity estimate for $\bm{\sigma}^{\rm p}$. Note that once Theorem \[SB\_PDE\_well\] has been established, we immediately obtain that ${\bm{u}}^{\rm p}\big|_{\Gamma_\epsilon} = {\bm{v}}^{\rm p}+\bm{\omega}^{\rm p}\times {\bm{X}}(s)$ is in $C^3(\Omega_\epsilon)$, since ${\bm{v}}^{\rm p}$ and $\bm{\omega}^{\rm p}$ are just constants in ${\mathbb{R}}^3$ and the fiber centerline ${\bm{X}}$ is in $C^3({\mathbb{T}})$. Given this $C^3$ Dirichlet data, $\bm{\sigma}^{\rm p}\in H^1(\Omega_\epsilon)$ follows by standard higher regularity arguments for the exterior Stokes Dirichlet boundary value problem (see the proof of Lemma V.4.3 in [@galdi2011introduction] or Theorem IV.5.8 in [@boyer2012mathematical]). Note that since ${\bm{X}}\in C^3({\mathbb{T}})$, $\bm{\sigma}^{\rm p}$ should in fact be even more regular, but the method we use to show Lemma \[high\_reg\] only allows us to quantify the $\epsilon$-dependence in the estimate for ${\left\lVert \nabla\bm{\sigma}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)}$. In particular, we can show the following bound on $\nabla \bm{\sigma}^{\rm p}$.

\[high\_reg\] Given $\Omega_\epsilon$ as in Section \[geometry\], the solution $\bm{\sigma}^{\rm p}$ to belongs to $H^1(\Omega_\epsilon)$ and satisfies $$\label{high_reg_eqn}
{\left\lVert \nabla\bm{\sigma}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \le {\left\lVert \nabla^2{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} + {\left\lVert \nabla p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \le \frac{C}{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{1/2}\big({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} + {\left\lVert p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big),$$ where $C$ depends on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.

The proof of the $\epsilon$-dependence in Lemma \[high\_reg\] is given in Appendix \[reg\_lem\].\
Using Lemma \[high\_reg\] and Corollary \[omegaP\], we have the higher regularity estimate $$\label{high_reg2}
\begin{aligned}
{\left\lVert \nabla\bm{\sigma}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} &\le \frac{C}{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{1/2}\bigg({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} + {\left\lVert p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \bigg) \le \frac{C}{\epsilon}{\left\lvert \log\epsilon \right\rvert}({\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert}).
\end{aligned}$$

Now, using that ${\mathcal{J}}_\epsilon(s,\theta)>0$ for each $(s,\theta)\in \Gamma_\epsilon$ and the surface measure ${\left\lvert \Gamma_\epsilon \right\rvert} = \int_{\mathbb{T}}\int_0^{2\pi} {\mathcal{J}}_\epsilon(s,\theta) d\theta ds =  \epsilon$, we have $$\begin{aligned}
{\left\lVert \bm{f}^{\rm p} \right\rVert}_{L^2({\mathbb{T}})}^2 &= \int_{\mathbb{T}}{\left\lvert \int_0^{2\pi} \bm{\sigma}^{\rm p}\bm{n} \, {\mathcal{J}}_\epsilon(s,\theta) d\theta  \right\rvert}^2 ds \le {\left\lvert \Gamma_\epsilon \right\rvert} \int_{\mathbb{T}}\int_0^{2\pi} {\left\lvert {\rm Tr}(\bm{\sigma}^{\rm p}) \right\rvert}^2 \, {\mathcal{J}}_\epsilon(s,\theta) d\theta \, ds \\
&\le C\epsilon^2{\left\lvert \log\epsilon \right\rvert}{\left\lVert \nabla\bm{\sigma}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)}^2 \le C{\left\lvert \log\epsilon \right\rvert}^3({\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert})^2.\end{aligned}$$ Here we have applied both the $L^2(\Gamma_\epsilon)$ trace inequality (Lemma \[trace2\]) and the higher regularity estimate in the last line.

With Lemma \[fp\_FT\], we are now equipped to show Lemma \[ups\_up\_err\].

The proof relies on estimates for the PDE satisfied by the difference between solutions to and . Letting ${\widetilde{{\bm{u}}}} = {\bm{u}}^{\rm p,s}-{\bm{u}}^{\rm p}$, ${\widetilde{p}} = p^{\rm p,s}- p^{\rm p} $, ${\widetilde{{\bm{v}}}} = {\bm{v}}^{\rm s}-{\bm{v}}^{\rm p}$, ${\widetilde{\bm{\omega}}} = \bm{\omega}^{\rm s}-\bm{\omega}^{\rm p}$, ${\widetilde{\bm{\sigma}}}=\bm{\sigma}^{\rm p,s} - \bm{\sigma}^{\rm p}$, we consider the following boundary value problem: $$\label{wtbu_eqn}
\begin{aligned}
-\Delta {\widetilde{{\bm{u}}}} +\nabla {\widetilde{p}} &=0, \quad {{\rm{div}\,}}{\widetilde{{\bm{u}}}} =0 \hspace{2.25cm} \text{in } \Omega_\epsilon \\
{\widetilde{{\bm{u}}}}({\bm{x}}) &= {\widetilde{{\bm{v}}}} + {\widetilde{\bm{\omega}}}\times {\bm{X}}(s) + \bm{R}(s), \qquad {\bm{x}}\in \Gamma_\epsilon \\
{\widetilde{{\bm{u}}}}({\bm{x}}) &\to 0 \hspace{4.3cm} \text{as }{\left\lvert {\bm{x}}\right\rvert}\to \infty \\
\int_{\Gamma_\epsilon} {\widetilde{\bm{\sigma}}}\bm{n} \; dS &= 0, \quad \int_{{\mathbb{T}}} {\bm{X}}(s)\times \bigg(\int_0^{2\pi}  {\widetilde{\bm{\sigma}}}\bm{n} \, {\mathcal{J}}_\epsilon(s,\theta) d\theta \bigg) ds = 0,
\end{aligned}$$ where $\bm{R}(s):={\rm Tr}({\bm{u}}^{\rm p,s})(s) - \big( {\bm{v}}^{\rm s} + \bm{\omega}^{\rm s}\times {\bm{X}}(s)\big)$ satisfies $$\label{Rest}
{\left\lVert \bm{R} \right\rVert}_{L^2({\mathbb{T}})} \le C\epsilon{\left\lvert \log\epsilon \right\rvert}^{3/2}{\left\lVert \bm{f}^{\rm s} \right\rVert}_{C^1({\mathbb{T}})},$$ by . We consider the variational form of : formally, multiplying by by ${\widetilde{{\bm{u}}}}$ and integrating by parts, we have $$\begin{aligned}
\int_{\Omega_\epsilon} 2|{\mathcal{E}}({\widetilde{{\bm{u}}}})|^2 \, d{\bm{x}}&= \int_{\Gamma_\epsilon} \big({\widetilde{{\bm{v}}}} + {\widetilde{\bm{\omega}}}\times {\bm{X}}(s) +\bm{R}(s) \big) \cdot({\widetilde{\bm{\sigma}}}\bm{n}) \; dS \\
&= {\widetilde{{\bm{v}}}} \cdot \int_{\Gamma_\epsilon}({\widetilde{\bm{\sigma}}}\bm{n}) \; dS +  {\widetilde{\bm{\omega}}}\cdot\int_{{\mathbb{T}}} {\bm{X}}(s)\times \bigg(\int_0^{2\pi} {\widetilde{\bm{\sigma}}}\bm{n}\, {\mathcal{J}}_\epsilon(s,\theta) d\theta \bigg) ds \\
&\qquad + \int_{\Gamma_\epsilon} \big({\widetilde{{\bm{v}}}} + {\widetilde{\bm{\omega}}}\times {\bm{X}}(s) +\bm{R}(s) \big) \cdot({\widetilde{\bm{\sigma}}}\bm{n}) \; dS \\
&=\int_{\mathbb{T}}\bm{R}(s) \cdot \big( \bm{f}^{\rm s} -\bm{f}^{\rm p} \big) \; ds \le {\left\lVert \bm{R} \right\rVert}_{L^2({\mathbb{T}})} \big( {\left\lVert \bm{f}^{\rm s} \right\rVert}_{L^2({\mathbb{T}})} + {\left\lVert \bm{f}^{\rm p} \right\rVert}_{L^2({\mathbb{T}})} \big) \\
&\le C\epsilon{\left\lvert \log\epsilon \right\rvert}^{3/2}{\left\lVert \bm{f}^{\rm s} \right\rVert}_{C^1({\mathbb{T}})}\big({\left\lVert \bm{f}^{\rm s} \right\rVert}_{L^2({\mathbb{T}})} + {\left\lvert \log\epsilon \right\rvert}^{3/2}({\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert})  \big). \end{aligned}$$ Here we have used and Lemma \[fp\_FT\] in the final line.

It would seem to make sense to try to bound the difference $\bm{f}^{\rm s} -\bm{f}^{\rm p}$ appearing in the second-to-last line by ${\left\lVert {\mathcal{E}}({\widetilde{{\bm{u}}}}) \right\rVert}_{L^2(\Omega_\epsilon)}$, or try to use an extension $\overline{\bm{R}}({\bm{x}})\in D^{1,2}(\Omega_\epsilon)$ with $\overline{\bm{R}}\big|_{\Gamma_\epsilon} = \bm{R}(s)$ and instead take ${\widetilde{{\bm{u}}}} - \overline{\bm{R}}$ as a test function in the above variational estimate to get rid of the boundary term. In either case, we run into difficulties in that we only have an $L^2({\mathbb{T}})$ estimate for $\bm{R}(s)$, when at least an $H^{1/2}({\mathbb{T}})$ estimate would be needed. However, as noted in Lemma \[high\_reg\], bounding the gradient of a function on $\Omega_\epsilon$ incurs an additional factor of $1/\epsilon$. By scaling, an $H^{1/2}({\mathbb{T}})$ estimate for $\bm{R}(s)$ would likely yield the same $\sqrt{\epsilon}$ factor appearing in Lemma \[ups\_up\_err\].

Now, using the $L^2({\mathbb{T}})$ trace inequality (Lemma \[trace1\]), the Korn inequality (Lemma \[korn\]), and Young’s inequality, along the with above ${\left\lVert {\mathcal{E}}({\widetilde{{\bm{u}}}}) \right\rVert}_{L^2(\Omega_\epsilon)}$ estimate, we have $$\begin{aligned}
{\left\lVert {\rm Tr}({\bm{u}}^{\rm p,s}) - ({\bm{v}}^{\rm p}+\bm{\omega}^{\rm p}\times{\bm{X}}) \right\rVert}_{L^2({\mathbb{T}})} &\le C{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert \nabla{\widetilde{{\bm{u}}}} \right\rVert}_{L^2(\Omega_\epsilon)} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert {\mathcal{E}}({\widetilde{{\bm{u}}}}) \right\rVert}_{L^2(\Omega_\epsilon)} \\
&\le C\sqrt{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{3/2} \big({\left\lvert \log\epsilon \right\rvert}^{1/2}{\left\lVert \bm{f}^{\rm s} \right\rVert}_{C^1({\mathbb{T}})} + {\left\lvert \bm{F} \right\rvert}+{\left\lvert \bm{T} \right\rvert} \big),\end{aligned}$$ yielding Lemma \[ups\_up\_err\].

Appendix
========

Here we provide proofs for the $L^2(\Gamma_\epsilon)$ trace inequality (Lemma \[trace2\]) and the higher regularity estimate (Lemma \[high\_reg\]).\
We first recall the following lemma, which will be used throughout the appendix.

*(Sobolev inequality)*\[sobolev\] Let $\Omega_{\epsilon}={\mathbb{R}}^3\backslash\overline{\Sigma_{\epsilon}}$ be as in Section \[geometry\]. For any ${\bm{u}}\in D^{1,2}(\Omega_{\epsilon})$, we have $$\label{sobolev_const}
\| {\bm{u}}\|_{L^6(\Omega_{\epsilon})} \le C\|\nabla{\bm{u}}\|_{L^2(\Omega_{\epsilon})}$$ where $C$ depends only on $c_\Gamma$ and $\kappa_{\max}$.

The proof of $\epsilon$-independence of $C$ appears in Appendix A.2.4 of [@closed_loop].

Proof of Lemma \[trace2\]
-------------------------

The proof of the $L^2(\Gamma_\epsilon)$ trace inequality follows the same outline as the proof of Lemma \[trace1\], contained in Appendix A.2.1 of [@closed_loop]. In particular, using the $\epsilon$-independent $C^2$-diffeomorphisms $\psi_j$ (defined in Appendix A.2.1, [@closed_loop]) which map segments of the curved slender body $\Sigma_\epsilon$ to a straight cylinder, it suffices to show the $\sqrt{\epsilon}{\left\lvert \log\epsilon \right\rvert}$ dependence of the trace constant for a straight cylinder.\
Accordingly, let ${\mathcal{D}}_\rho\subset{\mathbb{R}}^2$ denote the open disk of radius $\rho$ in ${\mathbb{R}}^2$, centered at the origin, and, for some $a<\infty$, define the cylindrical surface $\Gamma_{\epsilon,a}={\partial}{\mathcal{D}}_\epsilon\times [-a,a]$ and the cylindrical shell ${\mathcal{C}}_{\epsilon,a}= ({\mathcal{D}}_1\backslash\overline{{\mathcal{D}}_\epsilon}) \times [-a,a]$. Consider the function space $$D^{1,2}_\Gamma({\mathcal{C}}_{\epsilon,a})= \big\{ {\bm{u}}\in D^{1,2}({\mathcal{C}}_{\epsilon,a}) \; : \; {\bm{u}}\big|_{{\partial}{\mathcal{C}}_{\epsilon,a}\backslash\Gamma_{\epsilon,a}} = 0 \big\}.$$ As in the proof of Lemma \[trace1\], it suffices to show the $\sqrt{\epsilon}{\left\lvert \log\epsilon \right\rvert}$ dependence of the $L^2(\Gamma_{\epsilon,a})$ trace constant for functions belonging to $D^{1,2}_\Gamma({\mathcal{C}}_{\epsilon,a})$.\
By estimate (A.4) in [@closed_loop], any ${\bm{u}}\in C^1({\mathcal{C}}_{\epsilon,a}) \cap C^0(\overline{{\mathcal{C}}_{\epsilon,a}})\cap D^{1,2}_\Gamma({\mathcal{C}}_{\epsilon,a})$ satisfies $$\begin{aligned}
{\left\lvert {\rm Tr}({\bm{u}}) \right\rvert}^2 \le {\left\lvert \log\epsilon \right\rvert} \int_\epsilon^1 {\left\lvert \frac{{\partial}{\bm{u}}}{{\partial}\rho} \right\rvert}^2 \rho \, d\rho.\end{aligned}$$ Then, noting that the surface element on $\Gamma_{\epsilon,a}$ is simply $\epsilon$, we have $$\begin{aligned}
{\left\lVert {\rm Tr}({\bm{u}}) \right\rVert}_{L^2(\Gamma_{\epsilon,a})}^2 &= \int_{-a}^a \int_0^{2\pi} {\left\lvert {\rm Tr}({\bm{u}}) \right\rvert}^2 \epsilon \, d\theta\, ds \\
&\le \epsilon{\left\lvert \log\epsilon \right\rvert}  \int_{-a}^a \int_0^{2\pi}\int_\epsilon^1 {\left\lvert \frac{{\partial}{\bm{u}}}{{\partial}\rho} \right\rvert}^2 \rho \, d\rho\, d\theta\, ds  
\le \epsilon{\left\lvert \log\epsilon \right\rvert}{\left\lVert \nabla{\bm{u}}\right\rVert}_{L^2({\mathcal{C}}_{\epsilon,a})}^2.\end{aligned}$$ The same result for ${\bm{u}}\in D^{1,2}_\Gamma({\mathcal{C}}_{\epsilon,a})$ follows by density.

Proof of Lemma \[high\_reg\] {#reg_lem}
----------------------------

To determine the $\epsilon$-dependence of the constant in , it suffices to work locally near the slender body surface and show that Lemma \[high\_reg\] holds within an $\epsilon$-independent region about the slender body centerline. We define the region $$\label{mcO}
{\mathcal{O}} = \big\{{\bm{x}}\in \Omega_{\epsilon} \; : \; {\bm{x}}= {\bm{X}}(s) + \rho {\bm{e}}_{\rho}(s,\theta), \quad \epsilon < \rho<r_{\max} \big\},$$ where $r_{\max}$ is as in Section \[geometry\]. Within ${\mathcal{O}}$, we can use the orthonormal frame . We will use the notation ${\partial}_s,{\partial}_\theta,{\partial}_\rho$ to denote derivatives ${\partial}/{\partial}s$, ${\partial}/{\partial}\theta$, ${\partial}/{\partial}\rho$ with respect to the variables $s,\theta,\rho$, defined with respect to the orthonormal frame. We verify the $\epsilon$-dependence in the bound for $\nabla^2{\bm{u}}^{\rm p}$ and $\nabla p^{\rm p}$ in two parts: we first show an $L^2$ bound for derivatives $\nabla({\partial}_s{\bm{u}}^{\rm p})$, $\nabla({\partial}_\theta{\bm{u}}^{\rm p})$, ${\partial}_s p^{\rm p}$, and ${\partial}_\theta p^{\rm p}$ in directions tangent to the slender body surface $\Gamma_\epsilon$, and then use these bounds to estimate the derivatives $\nabla({\partial}_\rho{\bm{u}}^{\rm p})$, ${\partial}_\rho p^{\rm p}$ normal to $\Gamma_\epsilon$.\
We begin by estimating the tangential derivatives $\nabla({\partial}_s{\bm{u}}^{\rm p})$ and $\nabla({\partial}_\theta{\bm{u}}^{\rm p})$. Since the derivatives ${\partial}_s$ and ${\partial}_\theta$ with respect to the orthonormal frame do not commute with the “straight” differential operators $\nabla$ and ${{\rm{div}\,}}$, we will need to make use of the following commutator bounds.

\[comm\_ests\] For any function ${\bm{u}}\in D^{1,2}_0({\mathcal{O}})$ and for each of the differential operators $D={{\rm{div}\,}},\, \nabla,\, {\mathcal{E}}(\cdot)$, the following commutator estimates hold: $$\begin{aligned}
{\left\lVert [D,{\partial}_\theta]{\bm{u}}\right\rVert}_{L^2({\mathcal{O}})} &\le C{\left\lVert \nabla{\bm{u}}\right\rVert}_{L^2({\mathcal{O}})}, \quad {\left\lVert [D,{\partial}_s]{\bm{u}}\right\rVert}_{L^2({\mathcal{O}})} \le C{\left\lVert \nabla{\bm{u}}\right\rVert}_{L^2({\mathcal{O}})}, \end{aligned}$$ where the constant $C$ depends only on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.

We begin by denoting $$\begin{aligned}
 {\bm{e}}_\theta(s,\theta) &= -\sin\theta{\bm{e}}_{n_1}(s) + \cos\theta {\bm{e}}_{n_2}(s),\\
 u_\rho &= {\bm{u}}\cdot{\bm{e}}_\rho, \; u_\theta={\bm{u}}\cdot{\bm{e}}_\theta, \; u_s = {\bm{u}}\cdot{\bm{e}}_t.
 \end{aligned}$$

Then, with respect to the orthonormal frame , the divergence and gradient are given by $$\begin{aligned}
{{\rm{div}\,}}{\bm{u}}&= \frac{1}{1-\rho{\widehat{\kappa}}} \bigg( \frac{1}{\rho} \frac{{\partial}(\rho(1-\rho{\widehat{\kappa}})u_\rho)}{{\partial}\rho} + \frac{1}{\rho}\frac{{\partial}((1-\rho{\widehat{\kappa}})u_\theta)}{{\partial}\theta} + \frac{{\partial}u_s}{{\partial}s} \bigg) \\
\nabla {\bm{u}}&= {\bm{e}}_\rho(s,\theta)\frac{{\partial}{\bm{u}}}{{\partial}\rho}^{\rm T} + {\bm{e}}_\theta(s,\theta)\frac{1}{\rho}\frac{{\partial}{\bm{u}}}{{\partial}\theta}^{\rm T} + {\bm{e}}_t(s)\frac{1}{1-\rho{\widehat{\kappa}}} \bigg(\frac{{\partial}{\bm{u}}}{{\partial}s} -\kappa_3 \frac{{\partial}{\bm{u}}}{{\partial}\theta} \bigg)^{\rm T},\end{aligned}$$ where $$\label{kappahat}
{\widehat{\kappa}}(s,\theta) = \kappa_1(s)\cos\theta + \kappa_2(s)\sin\theta.$$

Direct computation of the commutators yields $$\begin{aligned}
[{{\rm{div}\,}},{\partial}_\theta]{\bm{u}}&= \frac{({\partial}_\theta{\widehat{\kappa}})}{1-\rho{\widehat{\kappa}}}\bigg( \rho\,{{\rm{div}\,}}{\bm{u}}- \frac{1}{\rho}\frac{{\partial}}{{\partial}\rho} \big(\rho^2 u_\rho \big)  - \frac{{\partial}u_\theta}{{\partial}\theta} \bigg) - \frac{({\partial}_\theta^2{\widehat{\kappa}})}{1-\rho{\widehat{\kappa}}}u_\theta  \\
[{{\rm{div}\,}},{\partial}_s]{\bm{u}}&= \frac{({\partial}_s{\widehat{\kappa}})}{1-\rho{\widehat{\kappa}}}\bigg( \rho\,{{\rm{div}\,}}{\bm{u}}- \frac{1}{\rho}\frac{{\partial}}{{\partial}\rho} \big(\rho^2 u_\rho \big)  - \frac{{\partial}u_\theta}{{\partial}\theta} \bigg) - \frac{({\partial}_\theta{\partial}_s{\widehat{\kappa}})}{1-\rho{\widehat{\kappa}}}u_\theta  \\
[\nabla,{\partial}_\theta]{\bm{u}}&= {\bm{e}}_\theta\frac{{\partial}{\bm{u}}}{{\partial}\rho}^{\rm T} - {\bm{e}}_\rho\frac{1}{\rho}\frac{{\partial}{\bm{u}}}{{\partial}\theta}^{\rm T} + {\bm{e}}_t\frac{\rho({\partial}_\theta{\widehat{\kappa}})}{(1-\rho{\widehat{\kappa}})^2} \bigg(\frac{{\partial}{\bm{u}}}{{\partial}s} -\kappa_3 \frac{{\partial}{\bm{u}}}{{\partial}\theta} \bigg)^{\rm T} \\
[\nabla,{\partial}_s]{\bm{u}}&= ({\partial}_s{\bm{e}}_\rho)\frac{{\partial}{\bm{u}}}{{\partial}\rho}^{\rm T} + ({\partial}_s{\bm{e}}_\theta)\frac{1}{\rho}\frac{{\partial}{\bm{u}}}{{\partial}\theta}^{\rm T} + \bigg({\bm{e}}_t \frac{\rho({\partial}_s{\widehat{\kappa}})}{1-\rho{\widehat{\kappa}}}+({\partial}_s{\bm{e}}_t) \bigg)\frac{1}{1-\rho{\widehat{\kappa}}} \bigg(\frac{{\partial}{\bm{u}}}{{\partial}s} -\kappa_3 \frac{{\partial}{\bm{u}}}{{\partial}\theta} \bigg)^{\rm T} \end{aligned}$$

Using and the orthonormal frame ODEs , we have $$\begin{aligned}
{\left\lvert {\partial}_\theta {\widehat{\kappa}} \right\rvert} &= {\left\lvert -\kappa_1\sin\theta+  \kappa_2\cos\theta \right\rvert}\le \kappa_{\max}, \quad {\left\lvert {\partial}_s {\widehat{\kappa}} \right\rvert} = {\left\lvert \kappa_1'\cos\theta+\kappa_2'\sin\theta \right\rvert}\le \xi_{\max} + 2(\kappa_{\max}+\pi), \\
{\left\lvert {\partial}_\theta^2{\widehat{\kappa}} \right\rvert} &= {\left\lvert -{\widehat{\kappa}} \right\rvert}\le \kappa_{\max}, \quad {\left\lvert {\partial}_\theta{\partial}_s {\widehat{\kappa}} \right\rvert} = {\left\lvert -\kappa_1'\sin\theta+  \kappa_2'\cos\theta \right\rvert}\le \xi_{\max}+ 2(\kappa_{\max}+\pi),  \\
{\left\lvert {\partial}_s {\bm{e}}_\rho \right\rvert} &= {\left\lvert -{\widehat{\kappa}}{\bm{e}}_t + \kappa_3{\bm{e}}_\theta \right\rvert}\le \kappa_{\max} +\pi, \quad {\left\lvert {\partial}_s{\bm{e}}_\theta \right\rvert} = {\left\lvert -({\partial}_\theta {\widehat{\kappa}}){\bm{e}}_t -\kappa_3{\bm{e}}_t \right\rvert}\le \kappa_{\max} +\pi, \\
{\left\lvert \frac{1}{1-\rho{\widehat{\kappa}}}  \right\rvert} &\le \frac{1}{1- r_{\max}\kappa_{\max}(\cos\theta+\sin\theta) } \le \frac{1}{1- \frac{1}{2\kappa_{\max}}\kappa_{\max}\sqrt{2} } \le 4 .\end{aligned}$$

Finally, noting that, by Lemma \[sobolev\], $${\left\lVert u_\theta \right\rVert}_{L^2({\mathcal{O}})} \le {\left\lvert {\mathcal{O}} \right\rvert}^{1/3} {\left\lVert {\bm{u}}\right\rVert}_{L^6({\mathcal{O}})} \le C{\left\lVert \nabla{\bm{u}}\right\rVert}_{L^2({\mathcal{O}})},$$ the desired $L^2(\Omega)$ bounds follow for each of $D={{\rm{div}\,}},\nabla$. The estimate for the symmetric gradient ${\mathcal{E}}({\bm{u}})$ then follows from the gradient commutator bound.

Now, to derive an estimate for $\nabla({\partial}_s{\bm{u}}^{\rm p})$, we will make use of Definition \[rigid\_weak\_p\] with a particular test function $\bm{\varphi}$, which we will construct here. First, we want our test function to be supported only within ${\mathcal{O}}$. We define a smooth cutoff function $$\label{Ocutoff}
\psi(\rho) = \begin{cases}
1, & \rho < r_{\max}/4 \\
0, & \rho > r_{\max}/2, 
\end{cases} \quad {\left\lvert \frac{{\partial}\psi}{{\partial}\rho}  \right\rvert} \le C,$$ where $C$ depends only on $r_{\max}$. Note that $\psi(\rho)$ commutes with both ${\partial}_\theta$ and ${\partial}_s$.\
We would like to use ${\partial}_s^2(\psi{\bm{u}}^{\rm p})$ as a test function in Definition \[rigid\_weak\_p\], but it will be more convenient to work with a function which vanishes on $\Gamma_\epsilon$. We therefore construct a correction ${\bm{g}}\in C^2(\Omega_\epsilon)$ supported only in ${\mathcal{O}}$ and satisfying $$\label{g_correction}
{\bm{g}}\big|_{\Gamma_\epsilon} = ({\partial}_s {\bm{u}}^{\rm p})\big|_{\Gamma_\epsilon} = \bm{\omega}^{\rm p}\times{\bm{e}}_t(s), \quad {\left\lVert \nabla {\bm{g}}\right\rVert}_{L^2({\mathcal{O}})} \le C{\left\lvert \bm{\omega} \right\rvert},$$ where $C$ depends on $c_\Gamma$ and $\kappa_{\max}$. To build ${\bm{g}}$, we follow a similar construction used in Section 4.1 of [@closed_loop]. We define $${\bm{g}}_0(\rho,\theta,s) = \begin{cases}
\bm{\omega}^{\rm p}\times {\bm{e}}_t(s) & \text{if } \rho<4\epsilon \\
0 & \text{otherwise}
\end{cases}$$ and take $${\bm{g}}(\rho,\theta,s):= \phi_\epsilon(\rho){\bm{g}}_0(\rho,\theta,s),$$ where $\phi_\epsilon(\rho)$ is the smooth cutoff defined in -. Note that ${\bm{g}}\in C^2$ and is supported within the region $${\mathcal{O}}_\epsilon := \big\{ {\bm{X}}(s) + \rho {\bm{e}}_\rho(s,\theta) \; : \; s\in {\mathbb{T}}, \, \epsilon \le \rho \le 4\epsilon, \, 0\le \theta <2\pi \big\},$$ where ${\left\lvert {\mathcal{O}}_\epsilon \right\rvert} \le C\epsilon^2$. Then, using and , we have $$\begin{aligned}
{\left\lVert \nabla {\bm{g}}\right\rVert}_{L^2({\mathcal{O}})} &\le \sqrt{{\left\lvert {\mathcal{O}}_\epsilon \right\rvert}}{\left\lVert \nabla {\bm{g}}\right\rVert}_{C({\mathcal{O}}_\epsilon)} \\
& \le \sqrt{{\left\lvert {\mathcal{O}}_\epsilon \right\rvert}}\bigg( {\left\lVert \frac{{\partial}\phi_\epsilon}{{\partial}\rho} \right\rVert}_{C({\mathcal{O}}_\epsilon)}{\left\lVert {\bm{g}}_0 \right\rVert}_{C({\mathcal{O}}_\epsilon)}+ {\left\lVert \frac{1}{1-\rho{\widehat{\kappa}}}\frac{{\partial}{\bm{g}}_0}{{\partial}s} \right\rVert}_{C({\mathcal{O}}_\epsilon)} \bigg) \le C{\left\lvert \bm{\omega}^{\rm p} \right\rvert}.\end{aligned}$$

Now, we could just use ${\partial}_s({\partial}_s(\psi{\bm{u}}^{\rm p}) -{\bm{g}})$ as a test function in Definition \[rigid\_weak\_p\], but it will actually be useful to include a second correction term in the following way. We consider ${\bm{z}}\in D^{1,2}_0({\mathcal{O}})$ satisfying $$\label{zee_def}
\begin{aligned}
{{\rm{div}\,}}{\bm{z}}&= {{\rm{div}\,}}(\psi{\partial}_s{\bm{u}}^{\rm p} -{\bm{g}}) \quad \text{in } {\mathcal{O}} \\
{\left\lVert \nabla {\bm{z}}\right\rVert}_{L^2({\mathcal{O}})} &\le C{\left\lVert {{\rm{div}\,}}(\psi{\partial}_s{\bm{u}}^{\rm p} -{\bm{g}}) \right\rVert}_{L^2({\mathcal{O}})} 
\end{aligned}$$ for $C$ depending only on $c_\Gamma$ and $\kappa_{\max}$. We know that such a ${\bm{z}}$ exists due to [@galdi2011introduction], Section III.3, and the constant $C$ is independence of $\epsilon$ due to Appendix A.2.5 of [@closed_loop]. Furthermore, since ${{\rm{div}\,}}{\bm{u}}^{\rm p}=0$, by Proposition \[comm\_ests\] we have $$\begin{aligned}
{\left\lVert {{\rm{div}\,}}(\psi{\partial}_s{\bm{u}}^{\rm p}-{\bm{g}}) \right\rVert}_{L^2({\mathcal{O}})} &\le {\left\lVert {{\rm{div}\,}}({\partial}_\theta{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})}+ C{\left\lVert {\partial}_\theta{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}+{\left\lVert \nabla {\bm{g}}\right\rVert}_{L^2({\mathcal{O}})} \\
 &\le {\left\lVert [{{\rm{div}\,}},{\partial}_\theta]{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} + C{\left\lVert {\partial}_\theta{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} + C{\left\lvert \bm{\omega}^{\rm p} \right\rvert} \\
 &\le C{\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} + C{\left\lvert \bm{\omega}^{\rm p} \right\rvert}.\end{aligned}$$ Here we have also used that ${\left\lVert {\partial}_\theta{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} \le {\left\lVert \rho \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} \le r_{\max}{\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}$. In particular, ${\bm{z}}$ satisfying also satisfies $$\label{zestimate}
{\left\lVert \nabla{\bm{z}}\right\rVert}_{L^2({\mathcal{O}})} \le C{\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}  + C{\left\lvert \bm{\omega}^{\rm p} \right\rvert}.$$

Using extension by zero to consider ${\bm{z}}$ as a function over all $\Omega_\epsilon$, we can now construct our desired test function for use in Definition \[rigid\_weak\_p\]. In particular, we will use the function ${\partial}_s({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}})$ in place of $\bm{\varphi}$ in Definition \[rigid\_weak\_p\]. Note that by definition of ${\bm{z}}$, this function may only belong to $L^2(\Omega_\epsilon)$. In this case, we can make sense of the following integration-by-parts argument using finite differences rather than full derivatives (see [@boyer2012mathematical], Section III.2.7 for construction of finite difference operators along a curved boundary). Thus we really only need ${\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}\in D^{1,2}(\Omega_\epsilon)$ to make sense of the following result. Note that in integrating by parts, we will also need to make use of the fact that, for $i=s,\theta$, $$\label{jacfac_i}
 {\partial}_i (d{\bm{x}}) = -\frac{\rho{\partial}_i{\widehat{\kappa}}}{1-\rho{\widehat{\kappa}}} d{\bm{x}}:= {\mathcal{J}}_i \, d{\bm{x}}, \quad {\left\lvert {\mathcal{J}}_i \right\rvert} \le C; \quad i=s,\theta ,$$ where $C$ depends on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.\
Then, using ${\partial}_s({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}})$ in Definition \[rigid\_weak\_p\], we have $$\begin{aligned}
0 &= \int_{{\mathcal{O}}} \bigg(2{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}\big({\partial}_s({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \big) - p^{\rm p} \,{{\rm{div}\,}}({\partial}_s({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}})) \bigg) \, d{\bm{x}}\\
&= \int_{{\mathcal{O}}} 2{\mathcal{E}}({\bm{u}}^{\rm p}): {\partial}_s{\mathcal{E}}({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, d{\bm{x}}+ \int_{{\mathcal{O}}} 2{\mathcal{E}}({\bm{u}}^{\rm p}): [{\mathcal{E}}(\cdot),{\partial}_s]({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, d{\bm{x}}\\
&\qquad - \int_{{\mathcal{O}}} p^{\rm p} \,{\partial}_s({{\rm{div}\,}}({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}})) \, d{\bm{x}}- \int_{{\mathcal{O}}} p^{\rm p} \,[{{\rm{div}\,}},{\partial}_s]({\partial}_s(\psi{\bm{u}}^{\rm p})- {\bm{g}}-{\bm{z}}) \, d{\bm{x}}\\
&= -\int_{{\mathcal{O}}} 2{\partial}_s{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, d{\bm{x}}-\int_{{\mathcal{O}}} 2{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, {\mathcal{J}}_s \, d{\bm{x}}\\
&\qquad + \int_{{\mathcal{O}}} 2{\mathcal{E}}({\bm{u}}^{\rm p}): [{\mathcal{E}}(\cdot),{\partial}_s]({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, d{\bm{x}}- \int_{{\mathcal{O}}} p^{\rm p} \,[{{\rm{div}\,}},{\partial}_s]({\partial}_s(\psi{\bm{u}}^{\rm p})- {\bm{g}}-{\bm{z}}) \, d{\bm{x}}\\
&= -\int_{{\mathcal{O}}} 2{\mathcal{E}}({\partial}_s{\bm{u}}^{\rm p}):{\mathcal{E}}({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, d{\bm{x}}- \int_{{\mathcal{O}}} 2{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, {\mathcal{J}}_s \, d{\bm{x}}\\
&\qquad +\int_{{\mathcal{O}}} 2[{\mathcal{E}}(\cdot),{\partial}_s]{\bm{u}}^{\rm p}: {\mathcal{E}}({\partial}_s(\psi{\bm{u}}^{\rm p})- {\bm{g}}-{\bm{z}}) \, d{\bm{x}}+ \int_{{\mathcal{O}}} 2{\mathcal{E}}({\bm{u}}^{\rm p}): [{\mathcal{E}}(\cdot),{\partial}_s]({\partial}_s(\psi{\bm{u}}^{\rm p})-{\bm{g}}-{\bm{z}}) \, d{\bm{x}}\\
&\qquad  - \int_{{\mathcal{O}}} p^{\rm p} \,[{{\rm{div}\,}},{\partial}_s]({\partial}_s(\psi{\bm{u}}^{\rm p})- {\bm{g}}-{\bm{z}}) \, d{\bm{x}}.\end{aligned}$$ Note that the first integral in the third line vanishes due to the definition of ${\bm{z}}$. In this way we we can avoid having to deal with a ${\partial}_s p^{\rm p}$ term in the resulting estimate.\
Then, using Proposition \[comm\_ests\], estimates and , and Lemma \[korn\], we have $$\begin{aligned}
{\left\lVert {\mathcal{E}}({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})}^2 &\le  C{\left\lVert {\mathcal{E}}({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} \big( {\left\lVert {\partial}_s {\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert {\mathcal{E}}({\bm{z}}) \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert {\mathcal{E}}({\bm{g}}) \right\rVert}_{L^2({\mathcal{O}})} \big) \\
&\quad + C{\left\lVert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} \big({\left\lVert {\mathcal{E}}(\psi{\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} +{\left\lVert {\mathcal{E}}({\bm{g}}) \right\rVert}_{L^2({\mathcal{O}})} +{\left\lVert {\mathcal{E}}({\bm{z}}) \right\rVert}_{L^2({\mathcal{O}})} \big) \\
&\quad + 2{\left\lVert [{\mathcal{E}}(\cdot),{\partial}_s]{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}\big( {\left\lVert {\mathcal{E}}(\psi{\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert {\mathcal{E}}({\bm{z}}) \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert {\mathcal{E}}({\bm{g}}) \right\rVert}_{L^2({\mathcal{O}})} \big) \\
&\quad +2{\left\lVert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} \big( {\left\lVert [{\mathcal{E}}(\cdot),{\partial}_s](\psi{\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})}  + {\left\lVert [{\mathcal{E}}(\cdot),{\partial}_s]({\bm{z}}) \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert [{\mathcal{E}}(\cdot),{\partial}_s]({\bm{g}}) \right\rVert}_{L^2({\mathcal{O}})} \big) \\
&\quad + {\left\lVert  p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} \big( {\left\lVert [{{\rm{div}\,}},{\partial}_s](\psi{\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert [{{\rm{div}\,}},{\partial}_s]({\bm{z}}) \right\rVert}_{L^2({\mathcal{O}})} +{\left\lVert [{{\rm{div}\,}},{\partial}_s]({\bm{g}}) \right\rVert}_{L^2({\mathcal{O}})} \big)  \\
&\le C({\left\lVert {\mathcal{E}}({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} +{\left\lvert \bm{\omega} \right\rvert})\big({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}+{\left\lVert p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} +{\left\lvert \bm{\omega}^{\rm p} \right\rvert} \big) \\
&\le \delta{\left\lVert {\mathcal{E}}({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})}^2 + C(\delta) \big({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}^2+{\left\lVert p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}^2 + {\left\lvert \bm{\omega}^{\rm p} \right\rvert}^2\big)\end{aligned}$$ for any $0<\delta\in{\mathbb{R}}$, by Young’s inequality. Taking $\delta=\frac{1}{2}$ and using Lemma \[korn\], we obtain $$\label{ps_est}
\begin{aligned}
{\left\lVert \nabla({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} &\le {\left\lVert {\mathcal{E}}({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} \le C \big({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}+{\left\lVert p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} + {\left\lvert \bm{\omega}^{\rm p} \right\rvert} \big) \\
& \le C{\left\lvert \log\epsilon \right\rvert}^{1/2} \big({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)}+{\left\lVert p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big),
\end{aligned}$$ where we have used Corollary \[omegaP\] to bound ${\left\lvert \bm{\omega}^{\rm p} \right\rvert}$. Here $C$ depends only on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.\
We may estimate ${\partial}_\theta{\bm{u}}^{\rm p}$ in a similar way. In fact, the construction of the analogous test function is simpler since $({\partial}_\theta{\bm{u}}^{\rm p})\big|_{\Gamma_\epsilon} = {\partial}_\theta({\bm{v}}+\bm{\omega}\times {\bm{X}}(s))=0$ and thus we do not need to correct for a nonzero boundary value. Following the same steps used to estimate ${\partial}_s{\bm{u}}^{\rm p}$, we obtain $$\label{ptheta_est}
{\left\lVert \nabla({\partial}_\theta{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} \le C \big({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)}+{\left\lVert p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big),$$ where $C$ depends only on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.\
In addition to the estimates and , we need bounds for the tangential derivatives ${\partial}_s p^{\rm p}$ and ${\partial}_\theta p^{\rm p}$ of the pressure. We begin by estimating ${\partial}_s p^{\rm p}$; the bound for ${\partial}_\theta p^{\rm p}$ is similar. Since we already know that ${\partial}_s p^{\rm p}\in L^2(\Omega_\epsilon)$, we may consider ${\widetilde{{\bm{z}}}} \in D^{1,2}_0({\mathcal{O}})$ satisfying $$\label{wtzee_def}
\begin{aligned}
{{\rm{div}\,}}{\widetilde{{\bm{z}}}} &= \psi {\partial}_s p^{\rm p} \quad \text{in }{\mathcal{O}}, \\
{\left\lVert \nabla{\widetilde{{\bm{z}}}} \right\rVert}_{L^2({\mathcal{O}})} &\le C{\left\lVert \psi {\partial}_s p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})},
\end{aligned}$$ where $\psi$ is as in . Again, we know that such a ${\widetilde{{\bm{z}}}}$ exists due to [@galdi2011introduction], Section III.3 and [@closed_loop], Appendix A.2.5.\
Using ${\partial}_s{\widetilde{{\bm{z}}}}$ as a test function in Definition \[rigid\_weak\_p\] (again, we can make sense of the following computation using finite differences, and thus only require ${\widetilde{{\bm{z}}}}\in D^{1,2}({\mathcal{O}})$), we have $$\begin{aligned}
0&= \int_{{\mathcal{O}}} \bigg(2{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\partial}_s {\widetilde{{\bm{z}}}}) - p^{\rm p}\, {{\rm{div}\,}}({\partial}_s{\widetilde{{\bm{z}}}})\bigg) \, d{\bm{x}}= \int_{{\mathcal{O}}}2{\mathcal{E}}({\bm{u}}^{\rm p}): {\partial}_s{\mathcal{E}}({\widetilde{{\bm{z}}}}) \, d{\bm{x}}\\
&\quad +  \int_{{\mathcal{O}}}2{\mathcal{E}}({\bm{u}}^{\rm p}): [{\mathcal{E}}(\cdot),{\partial}_s]{\widetilde{{\bm{z}}}} \, d{\bm{x}}- \int_{{\mathcal{O}}} p^{\rm p} \,{\partial}_s{{\rm{div}\,}}{\widetilde{{\bm{z}}}} \, d{\bm{x}}- \int_{{\mathcal{O}}} p^{\rm p} \, [{{\rm{div}\,}},{\partial}_s]{\widetilde{{\bm{z}}}} \, d{\bm{x}}\\
&= -\int_{{\mathcal{O}}}2{\partial}_s{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\widetilde{{\bm{z}}}}) \, d{\bm{x}}- \int_{{\mathcal{O}}}2{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\widetilde{{\bm{z}}}}) \, {\mathcal{J}}_s \, d{\bm{x}}+  \int_{{\mathcal{O}}}2{\mathcal{E}}({\bm{u}}^{\rm p}): [{\mathcal{E}}(\cdot),{\partial}_s]{\widetilde{{\bm{z}}}} \, d{\bm{x}}\\
&\quad - \int_{{\mathcal{O}}} p^{\rm p} \, [{{\rm{div}\,}},{\partial}_s]{\widetilde{{\bm{z}}}} \, d{\bm{x}}+ \int_{{\mathcal{O}}}({\partial}_s p){{\rm{div}\,}}{\widetilde{{\bm{z}}}} \, d{\bm{x}}+ \int_{{\mathcal{O}}} p^{\rm p} \,{{\rm{div}\,}}{\widetilde{{\bm{z}}}} \, {\mathcal{J}}_s \, d{\bm{x}}\\
&= \int_{{\mathcal{O}}}\psi({\partial}_s p)^2 \, d{\bm{x}}-\int_{{\mathcal{O}}}2{\partial}_s{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\widetilde{{\bm{z}}}}) \, d{\bm{x}}- \int_{{\mathcal{O}}}2{\mathcal{E}}({\bm{u}}^{\rm p}): {\mathcal{E}}({\widetilde{{\bm{z}}}}) \, {\mathcal{J}}_s \, d{\bm{x}}\\
&\quad  +  \int_{{\mathcal{O}}}2{\mathcal{E}}({\bm{u}}^{\rm p}): [{\mathcal{E}}(\cdot),{\partial}_s]{\widetilde{{\bm{z}}}} \, d{\bm{x}}- \int_{{\mathcal{O}}} p^{\rm p} \, [{{\rm{div}\,}},{\partial}_s]{\widetilde{{\bm{z}}}} \, d{\bm{x}}+ \int_{{\mathcal{O}}} p^{\rm p} \,{{\rm{div}\,}}{\widetilde{{\bm{z}}}} \, {\mathcal{J}}_s \, d{\bm{x}},\end{aligned}$$ where ${\mathcal{J}}_s \, d{\bm{x}}$ is as in and we have used . Then, using that $\psi^2\le \psi$, we have $$\begin{aligned}
{\left\lVert \psi{\partial}_s p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}^2 &\le 2{\left\lVert {\mathcal{E}}({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} {\left\lVert {\mathcal{E}}({\widetilde{{\bm{z}}}}) \right\rVert}_{L^2({\mathcal{O}})} + 2{\left\lVert [{\mathcal{E}}(\cdot),{\partial}_s]{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} {\left\lVert {\mathcal{E}}({\widetilde{{\bm{z}}}}) \right\rVert}_{L^2({\mathcal{O}})} \\
&\quad + C{\left\lVert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} {\left\lVert {\mathcal{E}}({\widetilde{{\bm{z}}}}) \right\rVert}_{L^2({\mathcal{O}})} + 2{\left\lVert {\mathcal{E}}({\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})}{\left\lVert  [{\mathcal{E}}(\cdot),{\partial}_s]{\widetilde{{\bm{z}}}} \right\rVert}_{L^2({\mathcal{O}})} \\
&\quad +{\left\lVert  p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} {\left\lVert  [{{\rm{div}\,}},{\partial}_s]{\widetilde{{\bm{z}}}} \right\rVert}_{L^2({\mathcal{O}})} + C{\left\lVert p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} {\left\lVert {{\rm{div}\,}}{\widetilde{{\bm{z}}}} \right\rVert}_{L^2({\mathcal{O}})}\\
&\le C\big({\left\lVert \nabla({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})} + {\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} +{\left\lVert  p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} \big) {\left\lVert  \psi{\partial}_s p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} \\
&\le \delta {\left\lVert  \psi{\partial}_s p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}^2 + C(\delta)\big({\left\lVert \nabla({\partial}_s{\bm{u}}^{\rm p}) \right\rVert}_{L^2({\mathcal{O}})}^2 + {\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}^2 +{\left\lVert  p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})}^2 \big)\end{aligned}$$ for $0<\delta\in {\mathbb{R}}$. Here we have used , , Proposition \[comm\_ests\], and Young’s inequality. Taking $\delta=\frac{1}{2}$ and using , we obtain $${\left\lVert \psi{\partial}_s p^{\rm p} \right\rVert}_{L^2({\mathcal{O}})} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2}\big( {\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} +{\left\lVert  p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big).$$ Then, using , within the region $${\mathcal{O}}' =  \bigg\{{\bm{x}}\in \Omega_{\epsilon} \; : \; {\bm{x}}= {\bm{X}}(s) + \rho {\bm{e}}_{\rho}(s,\theta), \quad \epsilon < \rho<\frac{r_{\max}}{4} \bigg\},$$ we have $$\label{psp_est}
{\left\lVert {\partial}_s p^{\rm p} \right\rVert}_{L^2({\mathcal{O}}')} \le C{\left\lvert \log\epsilon \right\rvert}^{1/2}\big( {\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} +{\left\lVert  p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big)$$ for $C$ depending only on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$.\
We can similarly use to show $$\label{pthetap_est}
{\left\lVert {\partial}_\theta p^{\rm p} \right\rVert}_{L^2({\mathcal{O}}')} \le C\big( {\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} +{\left\lVert  p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big).$$

Now we can use the tangential bounds , , , and to obtain an estimate for derivatives $\nabla({\partial}_\rho{\bm{u}}^{\rm p})$ normal to $\Gamma_\epsilon$. For this, we will use the full Stokes equations , written with respect to the orthonormal frame ${\bm{e}}_t$, ${\bm{e}}_\rho$, ${\bm{e}}_\theta$ in ${\mathcal{O}}$ as $$\begin{aligned}
- \Delta {\bm{u}}^{\rm p} +\nabla p^{\rm p} &= -\Delta {\bm{u}}^{\rm p} + \frac{{\partial}p^{\rm p}}{{\partial}\rho}{\bm{e}}_{\rho} + \frac{1}{\rho}\frac{{\partial}p^{\rm p}}{{\partial}\theta}{\bm{e}}_{\theta} +
\frac{1}{1-\rho{\widehat{\kappa}}}\bigg(\frac{{\partial}p^{\rm p}}{{\partial}s}-\kappa_3\frac{{\partial}p^{\rm p}}{{\partial}\theta}\bigg){\bm{e}}_t =0 \\
{{\rm{div}\,}}{\bm{u}}^{\rm p} &= \frac{1}{1-\rho{\widehat{\kappa}}}\bigg(\frac{1}{\rho}\frac{{\partial}(\rho (1-\rho{\widehat{\kappa}}) u_{\rho})}{{\partial}\rho}+\frac{1}{\rho}\frac{{\partial}((1-\rho{\widehat{\kappa}}) u_{\theta})}{{\partial}\theta} + \frac{{\partial}u_s}{{\partial}s} \bigg) = 0.\end{aligned}$$ Here ${\widehat{\kappa}}$ is as in and we recall the notation $u_\rho = {\bm{u}}^{\rm p}\cdot{\bm{e}}_\rho$, $u_\theta={\bm{u}}^{\rm p}\cdot{\bm{e}}_\theta$, $u_s={\bm{u}}^{\rm p}\cdot{\bm{e}}_t$.\
From the divergence-free condition on ${\bm{u}}^{\rm p}$, after multiplying through by $\rho (1-\rho{\widehat{\kappa}})$ and differentiating once with respect to $\rho$, we obtain $$\begin{aligned}
{\left\lVert  \frac{{\partial}^2 u_{\rho}}{{\partial}^2 \rho} \right\rVert}_{L^2(\mathcal{O})} &\le C\bigg({\left\lVert \frac{1}{\rho} \nabla {\bm{u}}^{\rm p}  \right\rVert}_{L^2(\mathcal{O})} + \bigg\|\frac{1}{\rho}\bigg\|_{L^{\infty}(\mathcal{O})}|\mathcal{O}|^{1/3}\big\|{\bm{u}}^{\rm p} \big\|_{L^6(\mathcal{O})} \\
&\hspace{2cm}+ \bigg\|\frac{1}{\rho}\frac{{\partial}}{{\partial}\rho}\bigg(\frac{{\partial}u_{\theta}}{{\partial}\theta}\bigg)\bigg\|_{L^2(\mathcal{O})} +\bigg\|\frac{{\partial}}{{\partial}\rho}\bigg(\frac{{\partial}u_s}{{\partial}s}\bigg)\bigg\|_{L^2(\mathcal{O})} \bigg)\\
&\le \frac{C}{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{1/2}\big( \|\nabla {\bm{u}}^{\rm p}\|_{L^2(\Omega_\epsilon)} +{\left\lVert  p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big),\end{aligned}$$ where we have used and along with the Sobolev inequality on $\Omega_\epsilon$.\
Furthermore, using the ${\bm{e}}_\rho$ component of $-\Delta {\bm{u}}^{\rm p} +\nabla p=0$, we have $$\begin{aligned}
\frac{{\partial}p^{\rm p}}{{\partial}\rho} &= (\Delta{\bm{u}}^{\rm p}) \cdot{\bm{e}}_{\rho} \\
&= \frac{1}{\rho (1-\rho{\widehat{\kappa}})}\frac{{\partial}}{{\partial}\rho}\left(\rho (1-\rho{\widehat{\kappa}})\frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}\rho}\right)\cdot{\bm{e}}_{\rho} +\frac{1}{\rho^2(1-\rho{\widehat{\kappa}})}\frac{{\partial}}{{\partial}\theta}\bigg((1-\rho{\widehat{\kappa}})\frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}\theta}\bigg)\cdot{\bm{e}}_{\rho} \\
&\hspace{4cm} +\frac{1}{1-\rho{\widehat{\kappa}}} \frac{{\partial}}{{\partial}s}\bigg( \frac{1}{1-\rho{\widehat{\kappa}}}\bigg[ \frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}s}- \kappa_3\frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}\theta} \bigg] \bigg)\cdot{\bm{e}}_{\rho} \\
&= \frac{1}{\rho (1-\rho{\widehat{\kappa}})}\frac{{\partial}}{{\partial}\rho}\left(\rho (1-\rho{\widehat{\kappa}})\frac{{\partial}u_{\rho}}{{\partial}\rho}\right) +\frac{1}{\rho^2(1-\rho{\widehat{\kappa}})}\frac{{\partial}}{{\partial}\theta}\bigg((1-\rho{\widehat{\kappa}})\frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}\theta}\bigg)\cdot{\bm{e}}_{\rho} \\
&\hspace{4cm} +\frac{1}{1-\rho{\widehat{\kappa}}} \frac{{\partial}}{{\partial}s}\bigg( \frac{1}{1-\rho{\widehat{\kappa}}}\bigg[ \frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}s}- \kappa_3\frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}\theta} \bigg] \bigg)\cdot{\bm{e}}_{\rho}, \end{aligned}$$ since ${\bm{e}}_{\rho}(s,\theta)$ does not vary with $\rho$. Therefore, using , , , and , along with the the bound on $\frac{{\partial}^2 u_{\rho}}{{\partial}\rho^2}$, we have $$\|\nabla p^{\rm p}\|_{L^2(\mathcal{O}')} \le \frac{C}{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{1/2}\big( \|\nabla {\bm{u}}^{\rm p}\|_{L^2(\Omega_\epsilon)} +{\left\lVert  p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big).$$

Finally, to estimate $\frac{{\partial}^2 u_j}{{\partial}\rho^2}$, $j=\theta,s$, we again use that $$\begin{aligned}
 \nabla p^{\rm p}\cdot{\bm{e}}_j &= (\Delta {\bm{u}}^{\rm p})\cdot {\bm{e}}_j(s,\theta) \\
&= \frac{1}{\rho (1-\rho{\widehat{\kappa}})}\frac{{\partial}}{{\partial}\rho}\bigg(\rho (1-\rho{\widehat{\kappa}})\frac{{\partial}u_j}{{\partial}\rho}\bigg) +\frac{1}{\rho^2(1-\rho{\widehat{\kappa}})}\frac{{\partial}}{{\partial}\theta}\bigg((1-\rho{\widehat{\kappa}})\frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}\theta}\bigg)\cdot{\bm{e}}_j  \\
&\hspace{4cm}+\frac{1}{1-\rho{\widehat{\kappa}}} \frac{{\partial}}{{\partial}s}\bigg( \frac{1}{1-\rho{\widehat{\kappa}}}\bigg[ \frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}s}- \kappa_3\frac{{\partial}{\bm{u}}^{\rm p}}{{\partial}\theta} \bigg] \bigg)\cdot{\bm{e}}_j , \quad j=\theta,s,\end{aligned}$$ since each of ${\bm{e}}_t(s)$, ${\bm{e}}_\rho(s,\theta)$ and ${\bm{e}}_\theta(s,\theta)$ are independent of $\rho$. Then we have $$\begin{aligned}
 \bigg\|\frac{{\partial}^2 u_j}{{\partial}\rho^2}\bigg\|_{L^2(\mathcal{O}')} &\le C \bigg( \bigg\|\frac{1}{\rho}\bigg\|_{L^{\infty}(\mathcal{O}')} \| \nabla {\bm{u}}^{\rm p}\|_{L^2(\mathcal{O}')} + \bigg\|\frac{{\partial}^2{\bm{u}}^{\rm p}}{{\partial}s^2}\bigg\|_{L^2(\mathcal{O}')} \\
 &\hspace{2cm} + \bigg\|\frac{{\partial}^2{\bm{u}}^{\rm p}}{{\partial}s{\partial}\theta}\bigg\|_{L^2(\mathcal{O}')} +\bigg\|\frac{{\partial}^2{\bm{u}}^{\rm p}}{{\partial}\theta^2}\bigg\|_{L^2(\mathcal{O}')}+\|\nabla p^{\rm p}\|_{L^2(\mathcal{O}')} \bigg) \\
&\le  \frac{C}{\epsilon}{\left\lvert \log\epsilon \right\rvert}^{1/2} \big({\left\lVert \nabla{\bm{u}}^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} + {\left\lVert p^{\rm p} \right\rVert}_{L^2(\Omega_\epsilon)} \big), \quad j=\theta, s, 
 \end{aligned}$$ where $C$ depends only on $c_\Gamma$, $\kappa_{\max}$, and $\xi_{\max}$. Altogether, we obtain Lemma \[high\_reg\]. $\square$

We note that the factor of $\frac{1}{\epsilon}$ in Lemma \[high\_reg\] is necessary. As a heuristic, we consider an infinite straight cylinder of radius $\epsilon$ and take ${\bm{u}}= (\frac{1}{\rho}-\frac{1}{\epsilon}) {\bm{e}}_{\theta}$, where ${\bm{e}}_{\theta}$ is now the (constant) angular vector in straight cylindrical coordinates, and $p\equiv$ constant. Ignoring decay conditions toward infinity along the cylinder, $({\bm{u}},p)$ solves the Stokes equations with ${\bm{u}}=0$ on the cylinder surface. Then $$\begin{aligned}
|\nabla^2{\bm{u}}| =\bigg| \frac{{\partial}^2}{{\partial}\rho^2} \frac{1}{\rho}\bigg| = \bigg|\frac{2}{\rho^3}\bigg| = \frac{2}{\rho}\big| \nabla {\bm{u}}\big|,\end{aligned}$$ and within the region $\epsilon < \rho \le 2 \epsilon$, we have $|\nabla^2{\bm{u}}| \ge \frac{1}{\epsilon}|\nabla {\bm{u}}|$.

[^1]: This research was supported in part by NSF grant DMS-1620316 and DMS-1516978, awarded to Y.M., and a Torske Klubben Fellowship, awarded to L.O. L.O. also thanks Dallas Albritton for helpful discussion.
