# Stokes flow: Wikis

Note: Many of our articles have direct quotes from sources you can cite, within the Wikipedia article! This article doesn't yet, but we're working on it! See more info or our list of citable articles.

# Encyclopedia

An object moving through a gas or liquid experiences a force in direction opposite to its motion. Terminal velocity is achieved when the drag force is equal in magnitude but opposite in direction to the force propelling the object. Shown is a sphere in Stokes flow, at very low Reynolds number.

Stokes flow (named after George Gabriel Stokes), also named creeping flow, is a type of fluid flow where advective inertial forces are small compared with viscous forces. The Reynolds number is low, i.e. $\textit{Re} \ll 1$. This is a typical situation in flows where the fluid velocities are very slow, the viscosities are very large, or the length-scales of the flow are very small. Creeping flow was first studied to understand lubrication. In nature this type of flow occurs in the swimming of microorganisms and sperm[1] and the flow of lava. In technology, it occurs in paint, MEMS devices, and in the flow of viscous polymers generally.

## Stokes equations

For this type of flow, the inertial forces are assumed to be negligible and the Navier–Stokes equations simplify to give the Stokes equations:

$\boldsymbol{\nabla} \cdot \mathbb{P} + \boldsymbol{f} = 0$

where $\mathbb{P}$ is the comoving stress tensor, and $\boldsymbol{f}$ an applied body force. There is also an equation for conservation of mass. In the common case of an incompressible Newtonian fluid, the Stokes equations are:

$\boldsymbol{\nabla}p = \mu \nabla^2 \boldsymbol{u} + \boldsymbol{f}$
$\boldsymbol{\nabla}\cdot\boldsymbol{u}=0$

Here $\boldsymbol{u}$ is the velocity of the fluid, $\boldsymbol{\nabla} p$ is the gradient of the pressure, and μ is the dynamic viscosity.

### Properties

The Stokes equations represent a considerable simplification of the full Navier–Stokes equations, especially in the incompressible Newtonian case.[2][3][4][5]

Instantaneity
A Stokes flow has no dependence on time other than through time-dependent boundary conditions. This means that, given the boundary conditions of a Stokes flow, the flow can be found without knowledge of the flow at any other time.
Time-reversibility
An immediate consequence of instantaneity, time-reversibility means then a time-reversed Stokes flow solves the same equations as the original Stokes flow. This property can sometimes be used (in conjunction with linearity and symmetry in the boundary conditions) to derive results about a flow without solving it fully. Time reversibility means that it is difficult to mix two fluids using creeping flow; a dramatic demonstration is possible of apparently mixing two fluids and then unmixing them by reversing the direction of the mixer.[6]

While these properties are true for incompressible Newtonian Stokes flows, the non-linear and sometimes time-dependent nature of non-Newtonian fluids means that they do not hold in the more general case.

### Methods of solution

#### By stream function

It can be shown that in planar 2-D systems, the stream function for an incompressible Newtonian Stokes flow satisfies the biharmonic equation $\nabla^4 \psi = 0$.

In the 3-D axisymmetric case, the Stokes stream function Ψ solves the equation E2Ψ = 0, where $E = {\partial^2 \over \partial r^2} + {\sin{\theta} \over r^2} {\partial \over \partial \theta} { 1 \over \sin{\theta}} {\partial \over \partial \theta}.$

#### By Papkovich–Neuber solution

The Papkovich–Neuber Solution represents the velocity and pressure fields of an incompressible Newtonian Stokes flow in terms of two harmonic potentials.

#### By boundary element method

Certain problems, such as the evolution of the shape of a bubble in a Stokes flow, are conducive to numerical solution by the Boundary Element method. This technique can be applied in both 2- and 3-dimensional flows.

#### By Green's function

The linearity of the Stokes equations in the case of an incompressible Newtonian fluid means that a Green's function for the equations can be found. The solution for the pressure p and velocity $\boldsymbol{u}$ due to a point force $\boldsymbol{F}\delta(\boldsymbol{x})$ acting at the origin with $|\boldsymbol{u}|$ and p vanishing at infinity is given by[7]

$\boldsymbol{u}(\boldsymbol{x}) = \boldsymbol{F} \cdot \mathbb{J}(\boldsymbol{x})$
$p(\boldsymbol{x}) = \frac{\boldsymbol{F}\cdot\boldsymbol{x}}{4 \pi |\boldsymbol{x}|^3}$

where

$\mathbb{J}(\boldsymbol{x}) = {1 \over 8 \pi \mu} \left( \frac{\mathbb{I}}{|\boldsymbol{x}|} + \frac{\boldsymbol{x}\boldsymbol{x}^\mathrm{T}}{|\boldsymbol{x}|^3} \right).$

is a second-rank tensor known as the Oseen tensor (after Carl Wilhelm Oseen).

The solution for a continuous force distribution (density) $\boldsymbol{f}(\boldsymbol{x})$ (again vanishing at infinity) can then be constructed by superposition:

$\boldsymbol{u}(\boldsymbol{x}) = \int \boldsymbol{f}(\boldsymbol{y}) \cdot \mathbb{J}(\boldsymbol{x} - \boldsymbol{y}) \, \mathrm{d}^3\!y$
$p(\boldsymbol{x}) = \int \frac{\boldsymbol{f}(\boldsymbol{y})\cdot(\boldsymbol{x}-\boldsymbol{y})}{4 \pi |\boldsymbol{x}-\boldsymbol{y}|^3} \, \mathrm{d}^3\!y$

## References

1. ^ Dusenbery, David B. (2009). Living at Micro Scale. Harvard University Press, Cambridge, Mass. ISBN 978-0-674-03116-6.
2. ^ Leal, L.G. (2007). Advanced Transport Phenomena:Fluid Mechanics and Convective Transport Processes.
3. ^
4. ^ Happel, J. & Brenner, H. (1981) Low Reynolds Number Hydrodynamics, Springer. ISBN 9001371159.
5. ^ Batchelor, G. (2000). Introduction to Fluid Mechanics.
6. ^ Dusenbery, David B. (2009). Living at Micro Scale, pp.46. Harvard University Press, Cambridge, Mass. ISBN 978-0-674-03116-6.
7. ^ Kim, S. & Karrila, S. J. (2005) Microhydrodynamics: Principles and Selected Applications, Dover. ISBN 0486442195.
• Ockendon, H. & Ockendon J. R. (1995) Viscous Flow, Cambridge University Press. ISBN 0521458811.