# Stokes problem

In fluid dynamics, Stokes problem also known as Stokes second problem or sometimes referred to as Stokes boundary layer or Oscillating boundary layer is a problem of determining the flow created by an oscillating solid surface, named after Sir George Stokes. This is considered as one of the simplest unsteady problem that have exact solution for the Navier-Stokes equations[1][2]. In turbulent flow, this is still named a Stokes boundary layer, but now one has to rely on experiments, numerical simulations or approximate methods in order to obtain useful information on the flow.

## Flow description[3][4]

Consider an infinitely long plate which is oscillating with a velocity ${\displaystyle U\cos \omega t}$ in the ${\displaystyle x}$ direction, which is located at ${\displaystyle y=0}$ in an infinite domain of fluid, where ${\displaystyle \omega }$ is the frequency of the oscillations. The incompressible Navier-Stokes equations reduce to

${\displaystyle {\frac {\partial u}{\partial t}}=\nu {\frac {\partial ^{2}u}{\partial y^{2}}}}$

where ${\displaystyle \nu }$ is the kinematic viscosity. The pressure gradient does not enter into the problem. The initial, no-slip condition on the wall is

${\displaystyle u(0,t)=U\cos \omega t,\quad u(\infty ,t)=0,}$

and the second boundary condition is due to the fact that the motion at ${\displaystyle y=0}$ is not felt at infinity. The flow is only due to the motion of the plate, there is no imposed pressure gradient.

### Solution[5][6]

The initial condition is not required because of periodicity. Since both the equation and the boundary conditions are linear, the velocity can be written as the real part of some complex function

${\displaystyle u=U\Re \left[e^{i\omega t}f(y)\right]}$

because ${\displaystyle \cos \omega t=\Re e^{i\omega t}}$.

Substituting this into the partial differential equation reduces it to ordinary differential equation

${\displaystyle f''-{\frac {i\omega }{\nu }}f=0}$

with boundary conditions

${\displaystyle f(0)=1,\quad f(\infty )=0}$

The solution to the above problem is

${\displaystyle f(y)=\exp \left[-{\frac {1+i}{\sqrt {2}}}{\sqrt {\frac {\omega }{\nu }}}y\right]}$
${\displaystyle u(y,t)=Ue^{-{\sqrt {\frac {\omega }{2\nu }}}y}\cos \left(\omega t-{\sqrt {\frac {\omega }{2\nu }}}y\right)}$

The disturbance created by the oscillating plate travels as the transverse wave through the fluid, but it is highly damped by the exponential factor. The depth of penetration ${\displaystyle \delta ={\sqrt {2\nu /\omega }}}$ of this wave decreases with the frequency of the oscillation, but increases with the kinematic viscosity of the fluid.

The force per unit area exerted on the plate by the fluid is

${\displaystyle F=\mu \left({\frac {\partial u}{\partial y}}\right)_{y=0}={\sqrt {\rho \omega \mu }}U\cos \left(\omega t-{\frac {\pi }{4}}\right)}$

There is a phase shift between the oscillation of the plate and the force created.

### Vorticity oscillations near the boundary

An important observation from Stokes' solution for the oscillating Stokes flow is that vorticity oscillations are confined to a thin boundary layer and damp exponentially when moving away from the wall.[7] This observation is also valid for the case of a turbulent boundary layer. Outside the Stokes boundary layer – which is often the bulk of the fluid volume – the vorticity oscillations may be neglected. To good approximation, the flow velocity oscillations are irrotational outside the boundary layer, and potential flow theory can be applied to the oscillatory part of the motion. This significantly simplifies the solution of these flow problems, and is often applied in the irrotational flow regions of sound waves and water waves.

### Fluid bounded by an upper wall

If the fluid domain is bounded by an upper, stationary wall, located at a height ${\displaystyle y=h}$, the flow velocity is given by

${\displaystyle u(y,t)={\frac {U}{2(\cosh 2\lambda h-\cos 2\lambda h)}}[e^{-\lambda (y-2h)}\cos(\omega t-\lambda y)+e^{\lambda (y-2h)}\cos(\omega t+\lambda y)-e^{-\lambda y}\cos(\omega t-\lambda y+2\lambda h)-e^{\lambda y}\cos(\omega t+\lambda y-2\lambda h)]}$

where ${\displaystyle \lambda ={\sqrt {\omega /(2\nu )}}}$.

### Flow due to an oscillating pressure gradient near a plane rigid plate

Stokes boundary layer due to the sinusoidal oscillation of the far-field flow velocity. The horizontal velocity is the blue line, and the corresponding horizontal particle excursions are the red dots.

The case for an oscillating far-field flow, with the plate held at rest, can easily be constructed from the previous solution for an oscillating plate by using linear superposition of solutions. Consider a uniform velocity oscillation ${\displaystyle u(\infty ,t)=U_{\infty }\cos \omega t}$ far away from the plate and a vanishing velocity at the plate ${\displaystyle u(0,t)=0}$. Unlike the stationary fluid in the original problem, the pressure gradient here at infinity must be a harmonic function of time. The solution is then given by

${\displaystyle u(y,t)=U_{\infty }\left[\,\cos \omega t-{\text{e}}^{-{\sqrt {\frac {\omega }{2\nu }}}y}\,\cos \left(\omega t-{\sqrt {\frac {\omega }{2\nu }}}y\right)\right],}$

which is zero at the wall z = 0, corresponding with the no-slip condition for a wall at rest. This situation is often encountered in sound waves near a solid wall, or for the fluid motion near the sea bed in water waves. The vorticity, for the oscillating flow near a wall at rest, is equal to the vorticity in case of an oscillating plate but of opposite sign.

## Stokes problem in cylindrical geometry

### Torsional oscillation

Consider an infinitely long cylinder of radius ${\displaystyle a}$ exhibiting torsional oscillation with angular velocity ${\displaystyle \Omega \cos \omega t}$ where ${\displaystyle \omega }$ is the frequency. Then the velocity for the steady state (i.e. neglecting the transient time) is given by[8]

${\displaystyle v_{\theta }=a\Omega \ \Re \left[{\frac {K_{1}(r{\sqrt {i\omega /\nu }})}{K_{1}(a{\sqrt {i\omega /\nu }})}}e^{i\omega t}\right]}$

where ${\displaystyle K_{1}}$ is the modified Bessel function of the second kind.

This solution can be expressed with real argument[9] as:

{\displaystyle {\begin{aligned}v_{\theta }\left(r,t\right)&=\Psi \left\lbrace \left[{\textrm {kei}}_{1}\left({\sqrt {R_{\omega }}}\right){\textrm {kei}}_{1}\left({\sqrt {R_{\omega }}}r\right)+{\textrm {ker}}_{1}\left({\sqrt {R_{\omega }}}\right){\textrm {ker}}_{1}\left({\sqrt {R_{\omega }}}r\right)\right]\cos \left(t\right)\right.\\&+\left.\left[{\textrm {kei}}_{1}\left({\sqrt {R_{\omega }}}\right){\textrm {ker}}_{1}\left({\sqrt {R_{\omega }}}r\right)-{\textrm {ker}}_{1}\left({\sqrt {R_{\omega }}}\right){\textrm {kei}}_{1}\left({\sqrt {R_{\omega }}}r\right)\right]\sin \left(t\right)\right\rbrace \\\end{aligned}}}

where

${\displaystyle \Psi =\left[{\textrm {kei}}_{1}^{2}\left({\sqrt {R_{\omega }}}\right)+{\textrm {ker}}_{1}^{2}\left({\sqrt {R_{\omega }}}\right)\right]^{-1}}$

and ${\displaystyle R_{\omega }}$is to the dimensionless oscillatory Reynolds number defined as ${\displaystyle R_{\omega }=\omega a^{2}/\nu }$, being ${\displaystyle \nu }$ the kinematic viscosity.

### Axial oscillation

If the cylinder oscillates in the axial direction with velocity ${\displaystyle U\cos \omega t}$, then the velocity field is

${\displaystyle u=U\ \Re \left[{\frac {K_{0}(r{\sqrt {i\omega /\nu }})}{K_{0}(a{\sqrt {i\omega /\nu }})}}e^{i\omega t}\right]}$

where ${\displaystyle K_{0}}$ is the modified Bessel function of the second kind.

## Stokes-Couette flow[10]

In the Couette flow, instead of the translational motion of one of the plate, an oscillation of one plane will be executed. If we have a bottom wall at rest at ${\displaystyle y=0}$ and the upper wall at ${\displaystyle y=h}$ is executing an oscillatory motion with velocity ${\displaystyle U\cos \omega t}$, then the velocity field is given by

${\displaystyle u=U\ \Re \left\{{\frac {\sin ky}{\sin kh}}\right\},\quad {\text{where}}\quad k={\frac {1+i}{\sqrt {2}}}{\sqrt {\frac {\omega }{\nu }}}.}$

The frictional force per unit area on the moving plane is ${\displaystyle -\mu U\Re \{k\cot kh\}}$ and on the fixed plane is ${\displaystyle \mu U\Re \{k\csc kh\}}$.

## References

1. ^ Wang, C. Y. (1991). "Exact solutions of the steady-state Navier-Stokes equations". Annual Review of Fluid Mechanics. 23: 159–177. Bibcode:1991AnRFM..23..159W. doi:10.1146/annurev.fl.23.010191.001111.
2. ^ Landau & Lifshitz (1987), pp. 83–85.
3. ^ Batchelor, George Keith. An introduction to fluid dynamics. Cambridge university press, 2000.
4. ^ Lagerstrom, Paco Axel. Laminar flow theory. Princeton University Press, 1996.
5. ^ Acheson, David J. Elementary fluid dynamics. Oxford University Press, 1990.
6. ^ Landau, Lev Davidovich, and Evgenii Mikhailovich Lifshitz. "Fluid mechanics." (1987).
7. ^ Phillips (1977), p. 46.
8. ^ Drazin, Philip G., and Norman Riley. The Navier–Stokes equations: a classification of flows and exact solutions. No. 334. Cambridge University Press, 2006.
9. ^ Rivero, M.; Garzón, F.; Núñez, J.; Figueroa, A. "Study of the flow induced by circular cylinder performing torsional oscillation". European Journal of Mechanics - B/Fluids. 78: 245–251. doi:10.1016/j.euromechflu.2019.08.002.
10. ^ Landau, L. D., & Sykes, J. B. (1987). Fluid Mechanics: Vol 6. pp. 88
Alexander Ramm

Alexander G. Ramm (born 1940 in St. Petersburg, Russia) is an American mathematician. His research focuses on differential and integral equations, operator theory, ill-posed and inverse problems, scattering theory, functional analysis, spectral theory, numerical analysis, theoretical electrical engineering, signal estimation, and tomography.

Geneviève Raugel

Geneviève Raugel (27 May 1951 – 10 May 2019) was a French mathematician working in the field of numerical analysis and dynamical systems.

Gifted (2017 film)

Gifted is a 2017 American drama film directed by Marc Webb and written by Tom Flynn. It stars Chris Evans, Mckenna Grace, Lindsay Duncan, Jenny Slate and Octavia Spencer. The plot follows an intellectually gifted 7-year-old who becomes the subject of a custody battle between her uncle and grandmother. The film was released on April 7, 2017, by Fox Searchlight Pictures, and grossed $43 million worldwide. List of fluid flows named after people This is a list of fluid flows named after people (eponymous flows). List of people from Hamilton, Ontario The following people were born in, residents of, or otherwise closely connected to the city of Hamilton, Ontario. Navier–Stokes existence and smoothness The Navier–Stokes existence and smoothness problem concerns the mathematical properties of solutions to the Navier–Stokes equations, a system of partial differential equations that describe the motion of a fluid in space. Solutions to the Navier–Stokes equations are used in many practical applications. However, theoretical understanding of the solutions to these equations is incomplete. In particular, solutions of the Navier–Stokes equations often include turbulence, which remains one of the greatest unsolved problems in physics, despite its immense importance in science and engineering. Even more basic properties of the solutions to Navier–Stokes have never been proven. For the three-dimensional system of equations, and given some initial conditions, mathematicians have not yet proved that smooth solutions always exist, or that if they do exist, they have bounded energy. This is called the Navier–Stokes existence and smoothness problem. Since understanding the Navier–Stokes equations is considered to be the first step to understanding the elusive phenomenon of turbulence, the Clay Mathematics Institute in May 2000 made this problem one of its seven Millennium Prize problems in mathematics. It offered a US$1,000,000 prize to the first person providing a solution for a specific statement of the problem:

Prove or give a counter-example of the following statement:In three space dimensions and time, given an initial velocity field, there exists a vector velocity and a scalar pressure field, which are both smooth and globally defined, that solve the Navier–Stokes equations.

Rayleigh problem

In fluid dynamics, Rayleigh problem also known as Stokes first problem is a problem of determining the flow created by a sudden movement of an infinitely long plate from rest, named after Lord Rayleigh and Sir George Stokes. This is considered as one of the simplest unsteady problem that have exact solution for the Navier-Stokes equations. The impulse movement of semi-infinite plate was studied by Keith Stewartson.

Sir George Stokes, 1st Baronet

Sir George Gabriel Stokes, 1st Baronet, (; 13 August 1819 – 1 February 1903) was an Anglo-Irish physicist and mathematician. Born in County Sligo, Ireland, Stokes spent all of his career at the University of Cambridge, where he was the Lucasian Professor of Mathematics from 1849 until his death in 1903. As a physicist, Stokes made seminal contributions to fluid mechanics, including the Navier-Stokes equations, and to physical optics, with notable works on polarization and fluorescence. As a mathematician, he popularised "Stokes' theorem" in vector calculus and contributed to the theory of asymptotic expansions. Stokes, along with Felix Hoppe-Seyler, first demonstrated the oxygen transport function of hemoglobin and showed color changes produced by aeration of hemoglobin solutions.

Stokes was made a baronet (hereditary knight) by the British monarch in 1889. In 1893 he received the Royal Society's Copley Medal, then the most prestigious scientific prize in the world, "for his researches and discoveries in physical science". He represented Cambridge University in the British House of Commons from 1887 to 1892, sitting as a Tory. Stokes also served as president of the Royal Society from 1885 to 1890 and was briefly the Master of Pembroke College, Cambridge.

Stokes approximation and artificial time

This article provides an error analysis of time discretization applied to spatially discrete approximation of the stationary and nonstationary Navier-Stokes equations. The nonlinearity of the convection term is the main problem in solving a stationary or nonstationary Navier-Stokes equation or Euler equation problems. Stoke incorporated ‘the method of artificial compressibility’ to solve these problems.

Tasneem M. Shah

Tasneem Mohammad Shah (Urdu: تسنیم محمد شاه‎), SI, TI, is a Pakistani scientist and a prominent mathematician who has made pioneering and instrumental research and contributions to the field of computational fluid dynamics (CFD) at Dr. A. Q. Khan Research Laboratories (KRL). Trained as an applied mathematician, his contributions include differential geometry, numerical analysis, information security, CFD-DEM model, hydrodynamics (of explosions), computer science, fluid mechanics, Vacuum Technology and CFD-DEM.

Tasneem Mohammad Shah was a pioneer and senior member of the Kahuta Project, and, along with Dr. G.D. Alam and Dr. Abdul Qadeer Khan, Dr. Tasneem Mohammad Shah have had worked out in the nuclear physics involved in gas centrifuge and the uranium based-device. As of today, he is a full professor and chairman of the Department of Mathematics at the Air University.

This page is based on a Wikipedia article written by authors (here).
Text is available under the CC BY-SA 3.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.