Third Year Project

This page contains the informations for the third year students working under the supervision of Dr Revell working on their research project.

Taylor Green Vortices

Taylor Green vortices is a classical numerical test case in order to check the energy conservation of a code. The initial solution has the following form
 \begin{equation} \left\{ \begin{array}{c} u_{1}=-\sin(kx_{1})\cdot\cos(kx_{2})\\ u_{2}=\cos(kx_{1})\cdot\sin(kx_{2})\end{array}\right.\label{eq: T-G vel field} \end{equation}
and the pressure is expressed by:
 \begin{equation} p=-\frac{1}{4}\left[\cos(2kx_{1})+\cos(2kx_{2})\right]\label{eq: p field T-G}\end{equation}
where $k=1,2,...$ is an arbitrary constant. The number of vortices in each direction is equal to $2\times k$. In this analysis the constant is set equal to 1.

Steps to do to run the simulation:

  1. Set the calculation as transient
  2. Laminar case (no turbulence model)
  3. inviscid case (viscosity $\mu = 0$)
  4. periodic boundary conditions:
    • periodic on the 4 edges of the square $ x = 0 \, \& \, x= 2 \pi $ and $ y = 0 \, \& \, y= 2 \pi $
    • symmetry on the front and rear
  5. Set a time step for which the CFL max is well below 1 (an example could be $\delta t = 1e-3$
  6. Set the convective scheme to a full central difference (CD).
  7. Set the residual for a variable resolution to 1e-8 (this is not the convergence criteria, it is the residual for the resolution of the linear system).
  8. Run the simulation for 30000 time steps.
  9. At the end of every time step compute the total kinetic enery as $ K=0.5(u_{1}^2+u_{2}^2)$

The conformal 60x60 mesh for the case can be download from here.

  • Code_Saturne: refer to this page to have the fortran to run the case
  • STAR-CD: here you can download the setup for the case for version 4.02. The archive cointain:
    • .MAC: macro that set up the case, from meshing generation to boundary conditions, to all numerical parameters
    • ufile: folder with the user fortran (NB in order to use the fortran the installation of Absoft fortran could be required, have a try and let me know if you have problem)
      • posdat: subroutine for the post processing (call at the end of every time step, the analogue of usproj under saturne)
      • initfi: subroutine to inizialize the field
  • FLUENT: try to follow the same step, if you have any problem to set up the case let me know

Current Tags:
create new tag
, view all tags
Topic attachments
I Attachment Action Size Date Who Comment
zipzip manage 137.0 K 2009-10-15 - 12:57 StefanoRolfo T&G vortices STAR-CD V4.02
elseunv tay60_conf.unv manage 2090.7 K 2009-10-15 - 12:39 StefanoRolfo T&G vortices mesh conformal
Edit | Attach | Print version | History: r5 < r4 < r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r3 - 2009-11-13 - 15:03:12 - StefanoRolfo
Main Web
22 Mar 2019


Manchester CfdTm

Ongoing Projects


Previous Projects


Useful Links:

User Directory
Photo Wall
Upcoming Events
Add Event

Computational Fluid Dynamics and Turbulence Mechanics
@ the University of Manchester
Copyright © by the contributing authors. Unless noted otherwise, all material on this web site is the property of the contributing authors.