Physics:Taylor–Green vortex

2D Contour Plot of Taylor Green Vortex

In fluid dynamics, the Taylor–Green vortex is an unsteady flow of a decaying vortex, which has an exact closed form solution of the incompressible Navier–Stokes equations in Cartesian coordinates. It is named after the British physicist and mathematician Geoffrey Ingram Taylor and his collaborator A. E. Green.

Vector plot of the Taylor-Green Vortex

Original work

In the original work of Taylor and Green, a particular flow is analyzed in three spatial dimensions, with the three velocity components [math]\displaystyle{ \mathbf{v}=(u,v,w) }[/math] at time [math]\displaystyle{ t=0 }[/math] specified by

[math]\displaystyle{ u = A \cos ax \sin by \sin cz, }[/math]

[math]\displaystyle{ v = B \sin ax \cos by \sin cz, }[/math]

[math]\displaystyle{ w = C \sin ax \sin by \cos cz. }[/math]

The continuity equation [math]\displaystyle{ \nabla \cdot \mathbf{v}=0 }[/math] determines that [math]\displaystyle{ Aa+Bb+Cc=0 }[/math]. The small time behavior of the flow is then found through simplification of the incompressible Navier–Stokes equations using the initial flow to give a step-by-step solution as time progresses.

An exact solution in two spatial dimensions is known, and is presented below.

Animation of a Taylor-Green Vortex using colour coded Lagrangian tracers

Incompressible Navier–Stokes equations

The incompressible Navier–Stokes equations in the absence of body force, and in two spatial dimensions, are given by

[math]\displaystyle{ \frac{\partial u}{\partial x}+ \frac{\partial v}{\partial y} = 0, }[/math]

[math]\displaystyle{ \frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = -\frac{1}{\rho} \frac{\partial p}{\partial x} + \nu \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} \right), }[/math]

[math]\displaystyle{ \frac{\partial v}{\partial t} + u\frac{\partial v}{\partial x} + v\frac{\partial v}{\partial y} = -\frac{1}{\rho} \frac{\partial p}{\partial y} + \nu \left( \frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} \right). }[/math]

The first of the above equation represents the continuity equation and the other two represent the momentum equations.

Taylor–Green vortex solution

In the domain [math]\displaystyle{ 0 \le x,y \le 2\pi }[/math], the solution is given by

[math]\displaystyle{ u = \sin x \cos y \,F(t), \qquad \qquad v = -\cos x \sin y \, F(t), }[/math]

where [math]\displaystyle{ F(t) = e^{-2\nu t} }[/math], [math]\displaystyle{ \nu }[/math] being the kinematic viscosity of the fluid. Following the analysis of Taylor and Green for the two-dimensional situation, and for [math]\displaystyle{ A=a=b=1 }[/math], gives agreement with this exact solution, if the exponential is expanded as a Taylor series, i.e. [math]\displaystyle{ F(t) = 1 - 2\nu t + O(t^2) }[/math].

The pressure field [math]\displaystyle{ p }[/math] can be obtained by substituting the velocity solution in the momentum equations and is given by

[math]\displaystyle{ p = \frac{\rho}{4} \left( \cos 2x + \cos 2y \right) F^2(t). }[/math]

The stream function of the Taylor–Green vortex solution, i.e. which satisfies [math]\displaystyle{ \mathbf{v} = \nabla \times \boldsymbol{\psi} }[/math] for flow velocity [math]\displaystyle{ \mathbf{v} }[/math], is

[math]\displaystyle{ \boldsymbol{\psi} = \sin x \sin y F(t)\, \hat{\mathbf{z}}. }[/math]

Similarly, the vorticity, which satisfies [math]\displaystyle{ \boldsymbol{\mathbf{\omega}} = \nabla \times \mathbf{v} }[/math], is given by

[math]\displaystyle{ \boldsymbol{\mathbf{\omega}} = 2\sin x \sin y \,F(t) \hat{\mathbf{z}}. }[/math]

The Taylor–Green vortex solution may be used for testing and validation of temporal accuracy of Navier–Stokes algorithms.

A generalization of the Taylor–Green vortex solution in three dimensions in described in.