sandbox/acastillo/readme

Andrés Castillo

My name is Andrés. I store in this location some of the work done in Basilisk. I mainly work with parametric instabilities, vortex instabilities and thermal convection.

Structure of this sandbox

This sandbox is divided in different parts:

Vortex Filament Method

This includes an implementation of a Vortex Filament Method (VFM) using Biot-Savart and some examples:

• Motion of a vortex ring
• Motion of 2 vortex rings
• Motion of 3 vortex rings
• Motion of an elliptical vortex ring

Output fields

This includes a series of routines to write the resulting fields into several different formats.

• output_vtu() - This format is compatible with the VTK XML file format for unstructured grids File formats for VTK version 4.2 which can be read using Paraview. Here, the results are written in raw binary. The unstructured grid is required to write results from quadtrees and octrees. If used in MPI, each MPI task writes its own file, which may be linked together using a *.pvtu file. An example is available here.

• output_xmf() This routine is compatible with the XDMF Model and Format which can be read using Paraview or Visit. Data is split in two categories: Light data and Heavy data. Light data is stored using eXtensible Markup Language (XML) to describe the data model and the data format. Heavy data is composed of large datasets stored using the Hierarchical Data Format HDF5. As the name implies, data is organized following a hierarchical structure. HDF5 files can be read without any prior knowledge of the stored data. The list of software capable of reading HDF5 files includes Visit, Paraview, Matlab, and Tecplot, to name a few. The type, rank, dimension and other properties of each array are stored inside the file in the form of meta-data. Additional features include support for a large number of objects, file compression, a parallel I/O implementation through the MPI-IO or MPI POSIX drivers. Using this format requires the HDF5 library, which is usually installed in most computing centers or may be installed locally through a repository. Linking is automatic but requires the environment variables HDF5_INCDIR and HDF5_LIBDIR, which are usually set when you load the module hdf5. An example is available here.

  • For testing the XDMF format try:

    sudo apt install libhdf5-dev hdf5-helpers  hdf5-tools
    export HDF5_INCDIR=/usr/include/hdf5/serial
    export HDF5_LIBDIR=/usr/lib/x86_64-linux-gnu/hdf5/serial
    make test_output4.tst

• For testing using MPI try:

    sudo apt install libhdf5-mpi-dev hdf5-helpers  hdf5-tools
    export HDF5_INCDIR=/usr/include/hdf5/openmpi
    export HDF5_LIBDIR=/usr/lib/x86_64-linux-gnu/hdf5/openmpi
    CC='mpicc -D_MPI=4' make test_output4.tst

• After the simulation, it is possible inspect the contents of the Heavy data from the terminal using:

  >> h5ls -r fields.h5
  /                        Group
  /000000                  Group
  /000000/Cell             Group
  /000000/Cell/T           Dataset {512, 1}
  /000000/Cell/pid         Dataset {512, 1}
  /000000/Cell/u.x         Dataset {512, 3}
  /000000/Geometry         Group
  /000000/Geometry/Points  Dataset {900, 3}
  /000000/Topology         Dataset {512, 8}
  ...

• output_vtkhdf() This routine is compatible with the VTKHDF File Format which can be read using Paraview or Visit. VTK HDF files start with a group called VTKHDF with two attributes: Version, an array of two integers and Type, a string showing the VTK dataset type stored in the file — in our case UnstructuredGrid. The data type for each HDF dataset is part of the dataset and it is determined at write time. In the diagram below, showing the HDF file structure for UnstructuredGrid, the rounded blue rectangles are HDF groups and the gray rectangles are HDF datasets. Each rectangle shows the name of the group or dataset in bold font and the attributes underneath with regular font. In our case, scalar and vector fields are stored as CellData. The unstructured grid is split into partitions, with a partition for each MPI rank.

  Image data VTKHDF File Format	Linear cell types found in VTK

  An example is available here and here.

Articles Submitted or in Preparation

Parametric instabilities

Articles in Publication

Vortex flows

Thermal convection