next up previous contents
Next: Cylinder Flow Up: Incompressible Viscous Magnetohydrodynamics Previous: Magnetic Pearson's Vortex

Orszag-Tang Vortex

The Orszag-Tang vortex is an initial value problem in a square periodic domain (length L). It has been investigated by [66] and [64] amongst others . It demonstrates that turbulent scales can result from a coherent initial condition with only two spatial frequencies. The initial conditions are:

An appropriate choice of initial magnetic stream potential is:

\begin{displaymath}
\phi = \frac{L}{2 \pi} cos(\frac{2 \pi y}{L}) - \frac{L}{4 \pi} cos(\frac{4 \pi x}{L})\end{displaymath}


 
Table 6.3: Simulation parameters for the incompressible MHD Orszag-Tang vortex simulation.
Parameter Value
Dimension 2d
Sv 100
Sr 100
$\Delta t$ 0.001
N-Range 10 to 25
KQuad 64
Method Galerkin

As the solution evolves in time we see that the initial vortex splits into two vortices. The primitive variable fields are shown in figure 6.4 at t=1. We see that sharp gradients build up and it took 64 quadrilateral elements at N=25 to resolve these features well. We see in figure 6.4 that the vorticity is smooth and well resolved, and that the velocity divergence is small. The vorticity exhibits complex small scale structures. Also in figure 6.4 we see that the divergence free conditions for the velocity field is met to 10-4 and the magnetic field is divergence free to 10-8. A summary of the simulation parameters is given in table 6.3.


\begin{figure}
\centerline{
\psfig {file=/crunch/crunch7/tcew/Thesis/Figures1/Ep...
 ...h7/tcew/Thesis/Figures1/Eps/tang.inc.yvel.t1.n25.eps,width=3.00in}
}\end{figure}

\begin{figure}
\centerline{
\psfig {file=/crunch/crunch7/tcew/Thesis/Figures1/Ep...
 ...h7/tcew/Thesis/Figures1/Eps/tang.inc.xmag.t1.n25.eps,width=3.00in}
}\end{figure}

  
Figure 6.4: Incompressible Orszag-Tang vortex (t=1, instantaneous fields). Top Left: x component of velocity, Top Right: y component of velocity, Middle: Pressure, Bottom Left: x component of magnetic field, Bottom Right: y component of magnetic field.
\begin{figure}
\centerline{
\psfig {file=/crunch/crunch7/tcew/Thesis/Figures1/Eps/tang.inc.ymag.t1.n25.eps,width=3.00in}
}\end{figure}


\begin{figure}
\centerline{
\psfig {file=/crunch/crunch7/tcew/Thesis/Figures1/Ep...
 .../tcew/Thesis/Figures1/Eps/tang.inc.magstreams_t1.eps,width=3.00in}
}\end{figure}

\begin{figure}
\centerline{
\psfig {file=/crunch/crunch7/tcew/Thesis/Figures1/Ep...
 .../tcew/Thesis/Figures1/Eps/tang.inc.divvel.t1.n25.eps,width=3.00in}
}\end{figure}

  
Figure 6.5: Incompressible Orszag-Tang vortex (t=1, instantaneous fields). Top Left: Velocity streamlines, Top Right: Magnetic streamlines, Middle Left: Vorticity, Middle right: Divergence of velocity, Bottom Left: Curl of magnetic field, Bottom Right: Divergence of magnetic field.
\begin{figure}
\centerline{
\psfig {file=/crunch/crunch7/tcew/Thesis/Figures1/Ep...
 .../tcew/Thesis/Figures1/Eps/tang.inc.divmag.t1.n25.eps,width=3.00in}
}\end{figure}


next up previous contents
Next: Cylinder Flow Up: Incompressible Viscous Magnetohydrodynamics Previous: Magnetic Pearson's Vortex
T. Warburton
10/24/1998