Data Sets


The following list contains a number of numerical data sets that we released for any use. If you use the data in your publications, please acknowledge the contributors by using the provided citation next to the download links. Most data sets are provided in NetCDF format, the Amira format and the VTK format. The NetCDF format is highly recommended, as it contains many useful properties, such as the physical units and the parameters of the simulation (obstacles, viscosity, Reynolds number if available).


Numerical Simulations


2D Unsteady Fluid Simulation Ensemble for Machine Learning

Image Example flows Regular grid: Amira (1.95 GB per file)
DownloadCitation

Ensemble simulation of 8,000 unsteady 2D flows with spatially-periodic boundary conditions. The simulation was done with Gerris flow solver and is carried out on a regular grid. The 8,000 flows are ordered from laminar configurations to turbulent configurations. The last column and row are left out, since they are the same as the first column and row. A detailed description of the data set can be found in Jakob et al.

Regular grid resolution (X x Y x T): 512 x 512 x 1001
Simulation domain: [0, 1] x [0, 1] x [0, 10]
Reynolds Number range: [1, 4096]
Kinematic viscosity range: [1e-5, 1e-4]

2D Unsteady Cylinder Flow with von Karman Vortex Street

Image Line integral convolution Image Vorticity Image Finite-time Lyapunov exponent (backward) Regular grid: NetCDF (440 MB), Amira (460 MB), VTK (448 MB)
Unstructured grid: VTK (4.83 GB)
DownloadCitation

Simulation of a viscous 2D flow around a cylinder. The fluid was injected to the left of a channel bounded by solid walls with a slip boundary condition. The simulation was done with Gerris flow solver and was resampled onto a regular grid. In the original simulation, the unstructured grid was adaptively discretized based on the vorticity. Over the course of the simulation, the characteristic von-Karman vortex street is forming. The image on the side shows a later time step, in which the street is fully formed. The vortices move with almost constant speed, except directly in the wake of the obstacle, where they accelerate.

Regular grid resolution (X x Y x T): 640 x 80 x 1501
Simulation domain: [-0.5, 7.5] x [-0.5, 0.5] x [0, 15]
Reynolds Number: 160
Kinematic viscosity: 0.00078125
Obstacle at (0,0) with radius: 0.0625

2D Unsteady Heated Cylinder with Boussinesq Approximation

Image LIC Image Vorticity Image FTLE
Regular grid: NetCDF (903 MB), Amira (922 MB), VTK (908 MB)
Unstructured grid: VTK (1.68 GB)
DownloadCitation

Simulation of a 2D flow generated by a heated cylinder. To solve the buoyancy problem, the Boussinesq approximation is used. The simulation was done with Gerris flow solver and was resampled onto a regular grid. The turbulent plume contains numerous small vortices that in part rotate around each other.

Grid resolution (X x Y x T): 150 x 450 x 2001
Simulation domain: [-0.5, 0.5] x [-0.5, 2.5] x [0, 20]
Obstacle at (0,-0.15) with radius: 0.0625

2D Unsteady Cylinder Flow Around Corners

Image Vorticity Image Finite-time Lyapunov exponent (backward) Regular grid: NetCDF (409 MB), Amira (414 MB), VTK (412 MB)
Unstructured grid: VTK (518 MB)
DownloadCitation

Simulation of a viscous 2D flow around two cylinders. The fluid was injected to the left of a channel bounded by solid walls with a slip boundary condition. The simulation was done with Gerris flow solver and was resampled onto a regular grid. In the original simulation, the unstructured grid was adaptively discretized based on the vorticity. Initially, a vortex street forms behind the first obstacle, which then flows around two corners. Behind each corner, a standing vortex forms. The latter one blocks half of the flow to the second obstacle, creating a one-sided vortex street.

Regular grid resolution (X x Y x T): 450 x 150 x 1501
Simulation domain: [-0.5, 5.5] x [-0.5, 1.5] x [0, 15]
Reynolds Number: 160
Kinematic viscosity: 0.00078125
Obstacles at (0,0) and (3,1) both with radius: 0.0625

3D Steady Cloud-topped Boundary Layer

Image Vorticity Regular grid: NetCDF (201 MB), Amira (202 MB), VTK (202 MB)
DownloadCitation

This data set contains a cloud resolving boundary layer simulation that was simulated with the UCLA‑LES model. In the domain, you will find a number of rising cumulus clouds. The data set is kindly provided by the German Climate Computing Center (DKRZ) and the Max Planck Institute for Meteorology (MPI‑M), and was simulated in the HD(CP)² project. The simulation uses double‑periodic boundary conditions, a homogeneous surface forcing, and takes its large‑scale information from the COSMO‑DE simulation model.

Regular grid resolution (X x Y x Z): 384 x 384 x 130
Simulation domain: [0, 10] x [0, 10] x [0, 3.2]

3D Unsteady Half Cylinder Ensemble

Image Vorticity for Reynolds number: 160 Image Vorticity for Reynolds number: 320 Image Vorticity for Reynolds number: 640 Image Vorticity for Reynolds number: 6400
Re = 160
 Regular grid: NetCDF (17.7 GB), Amira (18.7 GB), VTK (17.9 GB)
 Unstructured grid: VTK (1.2 GB)
Re = 320
 Regular grid: NetCDF (17.6 GB), Amira (18.6 GB), VTK (17.8 GB)
 Unstructured grid: VTK (1.4 GB)
Re = 640
 Regular grid: NetCDF (17.6 GB), Amira (18.6 GB), VTK (17.9 GB)
 Unstructured grid: VTK (1.5 GB)
Re = 6400
 Regular grid: NetCDF (17.6 GB), Amira (18.6 GB), VTK (17.9 GB)
 Unstructured grid: VTK (1.9 GB)
DownloadCitation

Small ensemble of numerical simulations of an incompressible 3D flow around a half cylinder. Each ensemble member was simulated with a different Reynolds number (Re). The simulations were done with Gerris flow solver and were resampled onto a regular grid. In the original simulation, the unstructured grid was adaptively discretized based on the vorticity. This flow might be useful for ensemble visualization techniques or for tests on flows with varying degree of turbulence.

Regular grid resolution (X x Y x Z x T): 640 x 240 x 80 x 151
Simulation domain: [-0.5, 7.5] x [-1.5, 1.5] x [-0.5, 0.5] x [0, 15]
Reynolds Numbers: 160, 320, 640, 6400

3D Unsteady Research Vessel Tangaroa

Image Vorticity in wake of the ship Regular grid: NetCDF (12.5 GB), Amira (13.1 GB), VTK (12.7 GB)
Unstructured grid: VTK (2.1 GB)
DownloadCitation

Simulation of an incompressible 3D flow around a CAD model of the research vessel Tangaroa. The simulation was done with Gerris flow solver and a region of interest was resampled onto a regular grid. This is one of the example simulations of the Gerris flow solver. In the original simulation, the unstructured grid was adaptively discretized based on the vorticity.

Regular grid resolution (X x Y x Z x T): 300 x 180 x 120 x 201
Simulation domain: [-0.35, 0.65] x [-0.3, 0.3] x [-0.5, -0.3] x [0, 2]



Analytic Vector Fields


It can come in very handy to test algorithms on analytic vector fields. For each analytic data set, the formular of the vector field itself and the first-order partial derivatives are provided in C++, Matlab and Python. For ease of use and testing, the vector fields are sampled onto regular grids, as well.

2D Unsteady Double Gyre

Image LIC with vorticity Image Finite-time Lyapunov exponent Regular grid: NetCDF (108 MB), Amira (113 MB), VTK (117 MB)
Code:C++,  Matlab,  Python
DownloadCitation

The double gyre is a periodic time‑dependent vector field, in which a separating boundary oscillates horizontally between two oppositely rotating vortices. The flow was introduced by Shadden et al. and became the prime benchmark for finite‑time Lyapunov exponents. The analytic formula has several parameters that steer magnitude (A), oscillation frequency (omega), and oscillation amplitude (eps). The resampled versions use the standard parameters listed below. The flow is also spatially periodic, which allows variations that contain saddles in the interior of the domain, such as the quad gyre with [0,2] × [-1,1] × [0,10].

Grid resolution (X × Y × T): 256 × 128 × 512
Space‑time domain: [0, 2] × [0, 1] × [0, 10]
Standard parameters: A = 0.1, omega = π/5, eps = 0.25

2D Unsteady Beads Problem

Image Vortex corelines in the Beads problem, image from here Regular grid: NetCDF (775 KB), Amira (13.8 MB), VTK (837 KB)
Code:C++,  Matlab,  Python
DownloadCitation

Wiebel et al. studied particle motion in a rotating petri‑dish, which proved to become a challenging benchmark for vortex coreline extraction, since the vortex center is moving on a circular path, which is not covered by Galilean invariant vortex extractors. An analytic approximation to this flow was for instance given by Weinkauf and Theisel, which is provided here. The image on the side shows pathlines that are attracted to the true coreline (blue), compared to a Galilean invariant coreline extraction (green).

Grid resolution (X × Y × T): 128 × 128 × 512
Simulation domain: [-2, 2] × [-2, 2] × [0, 2π]

2D Unsteady Four Rotating Centers

Image LIC with vorticity Image FTLE
Regular grid: NetCDF (54 MB), Amira (56.3 MB), VTK (57.5 MB)
Code:C++,  Matlab,  Python
DownloadCitation

This analytic data set contains four vortices. The flow is made unsteady by performing a uniform reference frame rotation. User parameters are a scaling factor of the magnitude (scale) and the speed of the reference frame rotation (al_t). The vortex centers are positioned at t = 0 at ± 2−1/2. The construction of the data set is described here in Section 7.2. This data set is a good benchmark for reference‑frame‑invariant flow feature extraction, since simple Galilean invariance is not enough.

Grid resolution (X × Y × T): 128 × 128 × 512
Space‑time domain: [-2, 2] × [-2, 2] × [0, 2π]
Standard parameters: scale = 1, al_t = 1

2D Unsteady Forced‑Damped Duffing Oscillator

Image LIC Image FTLE
Regular grid: NetCDF (30.1 MB), Amira (31 MB), VTK (29.9 MB)
Code:C++,  Matlab,  Python
DownloadCitation

This analytic vector field describes the phase space of a Duffing oscillator with forcing and damping. For the standard parameters, the vector field contains a saddle‑type periodic orbit.

Grid resolution (X × Y × T): 128 × 128 × 512
Space‑time domain: [-2, 2] × [-2, 2] × [0, 4]
Standard parameters: alpha = −0.25, beta = 0.4

2D Unsteady Cylinder Flow (Synthetic)

Image Line integral convolution Regular grid: NetCDF (311 MB), Amira (317 MB), VTK (313 MB)
DownloadCitation

This synthetic vector field represents a simple model of a von‑Kármán vortex street generation and was constructed by Jung, Tel and Ziemniak as co‑gradient to a stream function. The obstacle is positioned at (0,0) and has a radius = 1. In the LIC image on the side, the flow in the interior of the obstacle has not been set to zero. Note that only two vortices are present at the same time. In the files above, we sampled four periods onto a regular grid.

Grid resolution (X × Y × T): 450 × 200 × 500
Spatial domain: [-3, 7] × [-2, 2] × [1.107, 5.535]

3D Steady Tornado (Synthetic)

Image Streamlines from here Image Stream surfaces from here
Regular grid: NetCDF (16.9 MB), Amira (17.7 MB), VTK (17.0 MB)
DownloadCitation

This synthetic model of a tornado was created by Roger Crawfis. We scaled the flow to a larger domain and sampled it onto a regular grid. If higher resolutions are required, simply run the C code. Note, however, that we have scaled the domain.

Grid resolution (X × Y × Z): 128 × 128 × 128
Spatial domain: [-10, 10] × [-10, 10] × [-10, 10]



Data from Publications


The following collection lists derived data sets that have been used in publications.

3D Steady Opacity Optimization Data

Image Decoupled opacity optimization Points, lines, scalar fields: Amira (483 MB)
DownloadCitation

This folder contains some of the data sets used in the decoupled/Fourier opacity optimization papers. All files are stored in the Amira format, containing either points, lines, or a scalar field. The scalar field defines the importance everywhere in the domain and can be sampled onto the points/lines.



More Scientific Data Sets


Have not found the right data set yet? The following list contains a number of pointers.

Point Data

Scalar Fields

Vector Fields

How Can I Read the Data Sets?

In general, NetCDF is the recommend format. Readers exist for almost all languages (C++, Matlab, Python, ...). If you write C++ code and do not want to depend on other libraries, have a look at the Amira reader of Tino Weinkauf here. If you are compiling with VTK anyway, you can use the VTK readers. ParaView should have no problems to open the VTK files.

Let me know if you need help to open the data. We can learn from it and provide more information on this website to also help the next person.

Contribute to This List?

Send an email to ✉ tobias.guenther@fau.de if you would like a link to your data set here or are interested in sharing your data in our repository.