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
and the pressure is expressed by:
where
is an arbitrary constant. The number of vortices in each direction is equal to
. In this analysis the constant is set equal to 1.
Steps to do to run the simulation:
- Set the calculation as transient
- Laminar case (no turbulence model)
- inviscid case (viscosity )
- periodic boundary conditions:
- periodic on the 4 edges of the square and
- symmetry on the front and rear
- Set a time step for which the CFL max is well below 1 (an example could be
- Set the convective scheme to a full central difference (CD).
- 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).
- Run the simulation for 30000 time steps.
- At the end of every time step compute the total kinetic enery as
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
- New meshes with different resolution:
- : Star-CD macro to generate the unv mesh