Discussion for Home page

Neil Vaytet — Sun, 09 Feb 2014 19:28:38 -0000

               MINYDRO: Minimal Cell-Based Hydrodynamics
                   Written by Neil Vaytet (ENS Lyon)
                        Version 1.0 - July 2013

MINYDRO is a very simple code for 3D hydrodynamical simulations written in
Fortran 90. It is less than 700 lines long, and has a pretty neat structure!
Each cell is treated as a object which has a density, pressure, etc... and
interacts with other cells (i.e. its neighbours).

Being very short, it can't do as many things as other codes used in research
(such as MORPHEUS) can do. Indeed, it only has a cartesian grid, and currently
only has free-flow and periodic boundary conditions, but the simplicity of the
code makes it ideal for teaching purposes.

The code can be compiled with:
gfortran -o minydro minydro.f90 -ffree-line-length-none
If you wish to use OpenMP parallelisation, you can add the "-fopenmp" flag at
the end of your compilation command line.

The default setup is a Sedov explosion inside a 2D rectangular box. Some
parameters such as the simulation end time or the output frequency can be
changed at the beginning of the "minydro" program (lines 101 to 103).
To change the simulation setup, one needs to modify the "setup" subroutine. The
number of dimensions, the grid resolution, the simulation box size and the
boundary conditions can be modified at the start of the subroutine, while the
initial conditions of the gas conservative variables are modified starting from
line 220.

A small plotting program is also supplied in this archive. It reads the outputs
from the 2D simulation and writes density maps to ppm images. It can be
compiled simply by entering:
gfortran -o read_minydro_data read_minydro_data.f90
It can then be run using:
./read_minydro_data noutput_start [noutput_end increment]
where noutput_start is the output number from which to start plotting, and the
other two optional arguments are the output number at which plotting should
stop and the increment to be used for stepping through the outputs. If not
supplied, the default values are noutput_end = noutput_start and increment = 1.

This program is free to use by anyone and for any purpose, there is absolutely
no restrictions on use or distribution. It was designed mainly for teaching
purposes and it is made freely available in the hope that someone may find it
useful.

Home modified by Neil Vaytet

Neil Vaytet — Sun, 09 Feb 2014 17:21:10 -0000

Welcome to your wiki!

This is the default page, edit it as you see fit. To add a new page simply reference it within brackets, e.g.: [SamplePage].

The wiki uses Markdown syntax.

Project Members:

Neil Vaytet (admin)

Recent changes to Home

Discussion for Home page

Home modified by Neil Vaytet

Project Members: