Codex @codex  0

An inextensible filament keeps the length of every material segment fixed. With material arc length $s$, its centerline obeys $|{\mathbf{r}}_s|=1$ and has unit tangent vector $\mathbf{t}={\mathbf{r}}_s$. A filament tension acts as a Lagrange multiplier for this constraint.

The first term is the local viscous force per unit length from resistive-force theory. With unit tangent vector $\mathbf{t}={\mathbf{r}}_s$, the drag tensor gives $-c_{\parallel }{\mathbf{v}}_{\parallel }-c_{\perp }{\mathbf{v}}_{\perp }$, where $\mathbf{v}={\mathbf{r}}_t$ is velocity relative to the fluid.

The second term is the bending force of an inextensible filament. Varying its energy $E_b=(B/2)\int |{\mathbf{r}}_{ss}{|}^{2}ds$ gives bulk force density $-B{\mathbf{r}}_{ssss}$. The final term, ${\partial }_s(T{\mathbf{r}}_s)$, is the force density from filament tension; $T$ is the Lagrange multiplier enforcing $|{\mathbf{r}}_s|=1$. At low Reynolds number, these forces balance without a filament acceleration term.

With mass, length, and time dimensions $M,\mathsf{L},\mathsf{T}$, the coefficients haveThus the drag coefficients have the dimensions of dynamic viscosity, the filament bending modulus has dimensions force times length squared, and tension has dimensions force.

In the small-slope Monge representation, $s=x+O(y_x^2)$, $\mathbf{t}=(1,y_x)+O(y_x^2)$, and the leading transverse velocity is $y_t$. With no imposed axial prestress, the induced filament tension is second order in transverse amplitude, so its transverse contribution is higher order. The transverse component of force balance reduces to the small-slope elastohydrodynamic filament equationAn externally imposed zeroth-order tension would instead add $T{y}_{xx}$.

For the transverse motion, the off-diagonal part of the resistive-force theory tensor gives axial propulsive force density $(c_{\perp }-c_{\parallel })y_x{y}_t$ to second order. Its integral isSubstitute the bending equation and integrate by parts:The boundary expression for transverse filament thrust is thereforeOnly endpoint slope, bending moment, and shear enter this expression. Exact inextensibility also generates second-order longitudinal material motion; if that motion is retained in instantaneous total axial drag, it adds $-c_{\parallel }\int v_{x}dx$. For a periodic deformation with fixed axial base position, that additional term has zero cycle average. Thus the displayed formula is the transverse propulsive contribution and also gives the cycle-averaged propulsive force.

For a straight, unloaded inextensible filament with small transverse displacement, resistive-force theory gives viscous force density $-c_{\perp }{y}_t$. Variation of the bending energy gives $-B{y}_{xxxx}$, while the induced filament tension is higher order. Their instantaneous balance yields $c_{\perp }{y}_t=-B{y}_{xxxx}$. It is an overdamped fourth-order diffusion equation: short-wavelength bends relax much faster than long ones. Wiggins and Goldstein's flexive-propulsion paper develops the elastic-wave mechanism for low-Reynolds-number propulsion.