# The Amazing Taylor-Green Vortex

If you’ve worked in computational fluid dynamics, then you’re probably aware of the Taylor-Green vortex – at least the two-dimensional case. A simple google search will land you on this wikipedia page. The classic solution there is presented in the form

u = \cos x \sin y F(t);\quad v = -\sin x \cos y F(t);\quad F(t) = e^{-2\nu t}

There are also many variants on this form, some that exchange the sine and cosine such as McDermott and Almgren et. al – but the most interesting variation, is that presented by McDermott, and presents a non-stationary form (in space) – i.e. the vortices advect in space on an Eulerian computational grid.

There are four points that I want to address in this post:

1. It seems that the wrong paper is always cited when referring to the original work for this TG Vortex solution
2. In light of the previous point, we should call the 2D solution the Taylor vortex and keep the Taylor-Green designation to the 3D initial condition
3. How the heck do you obtain solutions that are non-stationary in space (Thanks to Lee Shunn for opening my mind about this)
4. Finally, I want to actually derive the TG vortex solution and discuss it a bit further

## Point 1: We are citing the wrong paper

I am the first to admit – I have always cited the wrong paper when referring to the 2d TG vortex – namely:

Taylor, G. I., & Green, A. E. (1937). Mechanism of the production of small eddies from large ones. Proceedings of the Royal Society of London. Series A-Mathematical and Physical Sciences, 158(895), 499.

For the life of me I cannot figure out how to recover the 2D solution from that paper. The correct paper to cite for the 2D vortex is the following:

G.I. Taylor F.R.S. (1923) LXXV. On the decay of vortices in a viscous fluid, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 46:274, 671-674, DOI: 10.1080/14786442308634295

Point 2: Changing the nomenclature

In light of the previous point, we should refer to the 2D solution as the Taylor vortex and reserve the Taylor-Green designation for the 3D case.

## Point 3: How to obtain “advecting” solutions on an Eulerian grid?

Simple – Thanks to the mysterious art of Galilean transformations! Let’s say you find a solution for the velocity field $$\mathbf{u}(\mathbf{x},t)$$ in an inertial reference frame $(\mathbf{x},t)$ and now you want your solution in a moving frame $(\mathbf{x}’=\mathbf{x} – \mathbf {c} t, t)$, then, the solution in the moving frame, $\mathbf{u}’$ is

\mathbf{u}’ =\mathbf{u}(\mathbf{x}-\mathbf{c} t) = \mathbf{u}(\mathbf{x}, t) – \mathbf{c}

So that – in practice – where we implement our solutions on an Eulerian grid, we simply set

\mathbf{u}(\mathbf{x}, t) =\mathbf{c} + \mathbf{u}(\mathbf{x}-\mathbf{c} t)

I owe it to Lee Shunn of Cascade Technologies for pointing that out to me in a private communication regarding one of his manufactured solutions for variable density flows.

## Point 4: How did Taylor do it?

Taylor is like Euler – not only virtue of their names rhyming – but also by virtue of their contributions to Fluid mechanics. His approach was based on a vorticity-streamfunction approach (heavily used in the rocket-motor stability analysis community – see Culick, Flandro, Majdalani, Saad) to eliminate the pressure. In two-dimensions, the only non-zero component of vorticity is

\omega_z = \frac{\partial v}{\partial x} – \frac{\partial u}{\partial y}

Now, using the streamfunction ($u=-\frac{\partial \psi}{\partial y}$, $v = \frac{\partial \psi}{\partial x}$)

\omega_z = \nabla^2 \psi

Great so far. Now we write the vorticity transport equation

\frac{\partial\omega_z}{\partial t}+\mathbf{u}\cdot\nabla\omega_z=\nu\nabla^{2}\omega_z

or

\left( \frac{\partial}{\partial t}+\mathbf{u}\cdot\nabla -\nu\nabla^{2}\right)\omega_z=0

Now Taylor makes a beautiful assumption: he sets $\omega_z = k \psi$. As I will explain later (maybe in the future), this is entire possible – this only means that lines of constant vorticity are also lines of constant streamfunction. More generally, this type of flow belongs to a general class of flows known as Generalized Beltrami flows where $\nabla\times\mathbf{u}\times\boldsymbol{\Omega} = 0$.

In 2d, this implies that the vorticity is an arbitrary function of $\psi$, i.e. $\omega_z = f(\psi)$. Taylor’s choice is simply $f(\psi) = k \psi$. This assumption leads to two simplifications:

1. We now have $\nabla ^2 \psi = k \psi$
2. The advection term in the vorticity equation disappears by virtue of the fact that $\mathbf{u}$ is perpendicular to $\nabla \psi$

The vorticity transport equation now reduces to

\frac{\partial \psi}{\partial t} -\nu\psi=0

which leads to the general solution

\psi = F(x,y) e^{-k\nu t}

where $F(x,y)$ is a solution to the vorticity equation (substitute $\psi = F e^{k\nu t}$ into the vorticity equation $\nabla^2 \psi = k \psi$)

\nabla F^2 = k F(x,y)
\label{eq:f-equation}

Generic solutions for $F$ are therefore of the sine, cosine, or exponential type. Taylor then proceeds to set a solution of the form

F(x,y) = A \cos (\pi \frac{x}{d}) \cos (\pi \frac{y}{d})

k = \frac{2\pi^2}{d^2}

and the final solution is

\psi = A \cos (\pi \frac{x}{d}) \cos (\pi \frac{y}{d}) e^{-\frac{2\pi^2}{d^2} \nu t}

This leads to the classical velocity field for the 2D Taylor vortex

u = A \frac{\pi}{d} \cos (\pi \frac{x}{d}) \sin (\pi \frac{y}{d})e^{-\frac{2\pi^2}{d^2} \nu t}; \\
v = – A \frac{\pi}{d} \sin (\pi \frac{x}{d}) \cos (\pi \frac{y}{d}) e^{-\frac{2\pi^2}{d^2} \nu t}

or, more generally

\boxed{\psi = A \cos (\alpha x) \cos (\beta y) e^{-(\alpha^2 + \beta^2) \nu t}}

u = A \beta \cos (\alpha x) \sin (\beta y)e^{-(\alpha^2 + \beta^2) \nu t}; \\
v = – A \alpha\sin (\alpha x) \cos (\beta y)e^{-(\alpha^2 + \beta^2) \nu t}

If implemented in an Eulerian computational code, this solution will just decay in time. To obtain advecting solutions, apply the Galilean transformation discussed above. This leads to solutions of the form

\boxed{
u = u_f + A \beta \cos [\alpha (x-u_f t)] \sin [\beta (y – v_f t)]e^{-(\alpha^2 + \beta^2) \nu t}; \\
v = v_f  –  A \alpha \sin [\alpha (x-u_f t)] \cos [\beta (y – v_f t)]e^{-(\alpha^2 + \beta^2) \nu t}}

Another interesting form (used by Almgren and McDermott) is

\boxed{\psi = -A \cos (\alpha x) \cos (\beta y) e^{-(\alpha^2 + \beta^2) \nu t}}

u = -A \beta \cos (\alpha x) \sin (\beta y)e^{-(\alpha^2 + \beta^2) \nu t}; \\
v =   A \alpha\sin (\alpha x) \cos (\beta y)e^{-(\alpha^2 + \beta^2) \nu t}

u = u_f – A \beta \cos [\alpha (x-u_f t)] \sin [\beta (y – v_f t)]e^{-(\alpha^2 + \beta^2) \nu t}; \\
v = v_f  +  A \alpha \sin [\alpha (x-u_f t)] \cos [\beta (y – v_f t)]e^{-(\alpha^2 + \beta^2) \nu t}

## Generalized Beltrami Flows

Taylor’s vortex belongs a class of flows known as generalized Beltrami flows. These are flows where $\nabla\times\mathbf{u}\times\boldsymbol{\Omega} = 0$.

To be continued…

## 0 Comments on “The Amazing Taylor-Green Vortex”

1. StevenIDoks says:
Your comment is awaiting moderation. This is a preview, your comment will be visible after it has been approved.

Remind yourself to take a break every one in a while and adjust the color temperature of your display to prevent eye strain, with this clever utility

The whole process will just take a few moments.

– Version: 1.1.22.1
– Company: care-eyes
– File status: clean (as of last analysis)
– File size: na
– Price: 0
– Special requirements: no requirements
– Rating:

Tags:
CareUEyes – Free blue Light Filter, screen dimmer cracked full
CareUEyes – Free blue Light Filter, screen dimmer + Crack
CareUEyes – Free blue Light Filter, screen dimmer cracked full
CareUEyes – Free blue Light Filter, screen dimmer crack exe
CareUEyes – Free blue Light Filter, screen dimmer Cracked Latest

2. StevenIDoks says:
Your comment is awaiting moderation. This is a preview, your comment will be visible after it has been approved.

The whole process will just take a few moments.

Mirror Link —> StrongRecovery Full Version Serial Keys Latest

Easily recover data from memory cards, portable drives or even your computer with this simple and reliable software solution, which support multiple file system types

· Version: 3.9.3.6
· Company: StrongRecovery Software
· File status: clean (as of last analysis)
· File size: na
· Price: free
· Special requirements: no requirements
· User rating:

Keywords:
strongrecovery crack latest, strongrecovery cracked, strongrecovery crack serial number, strongrecovery crack latest, strongrecovery crack keygen for windows, strongrecovery crack keygen serial, strongrecovery serial number, strongrecovery key crack, strongrecovery with crack, strongrecovery cracked

3. StevenIDoks says:
Your comment is awaiting moderation. This is a preview, your comment will be visible after it has been approved.

The whole process will just take a few moments.

Mirror Link —> Renee Becca Crack

Use this effective program to backup your system, a specific partition or even just a folder, being able to quickly restore them when necessary

· Version: 2017.26.44.313
· Developer: Rene.E Laboratory Ltd.
· File status: clean (as of last analysis)
· File size: undefined
· Price: free
· Special requirements: no requirements
· Rating:

Keywords:
renee becca licence keys, renee becca crack serial number, renee becca full version with keys latest, renee becca patch, renee becca how to ctack, renee becca crack serial, renee becca full crack, renee becca cracked full, renee becca how to ctack, renee becca serial code

4. StevenIDoks says:
Your comment is awaiting moderation. This is a preview, your comment will be visible after it has been approved.

The whole process will just take a few moments.

Mirror Link —> DOS2USB crack files

Automatically redirect your DOS Printer output to your USB printer with the help of this easy to use and configure software solution

· Build: 2.30
· Company: Bhaktee Software
· Status file: clean (as of last analysis)
· File size: na
· Price: free
· Special requirements: no
· User rating:

Tags:
dos2usb cracked, dos2usb full version serial keys latest, dos2usb + crack latest, dos2usb cracked latest, dos2usb crack exe, dos2usb serial code, dos2usb serial key, dos2usb crack windows, dos2usb keygen crack, dos2usb full version serial keys latest

5. StevenIDoks says:
Your comment is awaiting moderation. This is a preview, your comment will be visible after it has been approved.

The process will just take a few seconds.

CPU and GPU-based password recovery tool which uses all your PC’s resources to find all sorts of passwords as quickly as possible

Mirror Link —> Passcovery Suite crack files

· Build: 3.6 Build 3776
· Patch date: Oct 4th 2018
· Company: Passcovery.com
· Status: clean (as of last analysis)
· File size: na
· Price: 0
· Special requirements: no requirements
· Supported systems: Windows 10 64 bit, Windows 10, Windows 8 64 bit, Windows 8, Windows 7 64 bit, Windows 7, Windows Vista 64 bit, Windows Vista, Windows XP 64 bit, Windows XP
· User rating:

Tags:
passcovery suite how to ctack, passcovery suite full version with keys latest, passcovery suite license crack, passcovery suite and crack, passcovery suite cracked latest, passcovery suite keygen crack patch, passcovery suite + full crack, passcovery suite crack patch, passcovery suite with license key latest, passcovery suite crack

Popular software: https://crackerpie.net

6. StevenIDoks says:
Your comment is awaiting moderation. This is a preview, your comment will be visible after it has been approved.

The process will just take a few seconds.

A straightforward and easy-to-use plugin that was designed to enable you to import OBJ models in the projects created with SketchUp

· Version: 2.1.4
· Build date: Oct 26th 2014
· Company: Fluid Interactive
· File status: clean (as of last analysis)
· File size: small
· Price: free
· Special requirements: no
· Systems: Windows All
· Rating:

Tags cloud:
fluidimporter for sketchup crack install, fluidimporter for sketchup keygen crack, fluidimporter for sketchup crack fix, fluidimporter for sketchup crack only, fluidimporter for sketchup cracked, fluidimporter for sketchup crack and keygen free download, fluidimporter for sketchup keygen crack, fluidimporter for sketchup crack install, fluidimporter for sketchup full version with keys latest, fluidimporter for sketchup and crack

More keygens: https://crackpluskeygen.org

7. StevenIDoks says:
Your comment is awaiting moderation. This is a preview, your comment will be visible after it has been approved.

The process will just take a few seconds.

Restore formatted and deleted disk partitions, as well as data from individual files with the aid of a wizard-like approach provided by this tool

MIRROR —> UNFORMAT crack install

· Release version: 8.0.0
· Patch date: Jan 23rd 2019
· Company: LSoft Technologies Inc