Program title: DNSLab
Catalogue identifier: AEZX_v1_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEZX_v1_0.html
Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 227177
No. of bytes in distributed program, including test data, etc.: 40553359
Distribution format: tar.gz
Programming language: Matlab (2014&2015) Compatible with Octave (tested on v-3.x.x and v-4.0.0).
Computer: Any work-station or laptop computer where Matlab is running.
Operating system: Linux and Windows, and Mac OS X.
Has the code been vectorized or parallelized?: Vectorization in Matlab is utilized.
RAM: The peak memory consumption for the distributed default test cases is about 600 Mb for Channel3dLab (723 grid) and about 150 Mb for NS2dLab (2562 grid). Matlab as such requires about 350 Mb of memory without GUI.
Classification: 12, 4.3, 4.6.
Nature of problem: Numerical solution of the Navier–Stokes equations in turbulent state is demonstrated in Matlab environment for two test problems: turbulent 3d channel flow and 2d periodic array of vortices. The high-level, interpreted language Matlab enables the solution of turbulent flows using compact and short code syntax. Both of the problems are of high relevance in numerical test phases of research, and in education and numerical simulation of turbulence.
Solution method: The two solvers of DNSLab are based on the fractional step projection methods utilizing finite differences in 3d, and Fourier pseudo-spectral method in 2d. The time integrator in both solvers is the classical fourth order Runge–Kutta scheme.
Restrictions: As such, the implemented codes are limited to either periodic or simple channel flow configurations. In general, stability of fluid flow solvers is dependent on the Reynolds number, Courant number, Courant–Friedrichs–Lewy number, the initial data, as well as linear solver settings.
Additional comments: !!!!! The distribution file for this program is over 40 Mbytes and therefore is not delivered directly when download or Email is requested. Instead a html file giving details of how the program can be obtained is sent. !!!!!
Running time: Results for test cases can be produced in tens of seconds in 2d and in a few minutes in 3d. Runtime grid size dependence is investigated in the manuscript.