VH-1
The Virginia Numerical Bull Session ideal hydrodynamics PPMLR
Version 1.0 October 1991
This version was removed from the ``loop'' in the continued development
of VH-1 on Sep 1 1990, but has followed a similar evolution until this
release date. Please address questions and comments to
John Blondin
(Email: john_blondin@ncsu.edu).
VH-1 is a multidimensional ideal compressible hydrodynamics code
written in FORTRAN, for use on desktop workstations or supercomputers.
It is based on the Lagrangian remap version of the Piecewise Parabolic
Method (PPM) developed by Paul Woodward and Phil Collela.
You are welcome to download the code and do with it what you wish. We would
appreciate (1) your acknowledging the use of the code in publications, and (2)
sending a quick
email
to let us know who is using the code and what it is being used
for. Thank you for your cooperation.
However, keep in
mind that (1) this code does not come with any guarantees, and (2) you
should not try to use this code without first understanding how it
functions (and what computational hydrodynamics is all about)!
We strongly recommend that you read Collela & Woodward (J. Comp. Phys., 54, 174)
thoroughly before trying to use VH-1.
There are several files available at the ftp site relating to VH-1:
- FORTRAN Code
- A tar file of the NEW VH-1 code,
including all subroutines, files for all test problems,
and a makefile for UNIX compilation. This version of the code uses NCSA's HDF
utilities for data storage.
- Isothermal version of new VH-1 code.
- MPI version of F90 VH-1 code using netcdf for output.
- Linux HDF Libraries
- Test Page
- Check out our new testpage full of new movie and data
files for the VH-1 test problems.
The following files pertain to the original version of VH-1:
- Old FORTRAN Code
- A compressed tar file of the VH-1 code,
including all subroutines and
a makefile for UNIX compilation. This version of the code uses NCSA's HDF
utilities for data storage.
- Old TEST Subroutines
- A compressed tar file of subroutines preconfigured to run the
test problems described in the user's guide.
- DATA Files
- A compressed tar file of HDF output generated by the test subroutines
supplied above (to see if you did it right!).
- USER'S GUIDE
- An admittedly incomplete user's guide for VH-1 in postscript format.
- It is also available in html format.
- ISOTHERMAL VH-1
- An isothermal version of VH-1 (much faster without the need for an iterative
Riemann solver).
Getting Started:
Assuming you know hydrodynamics backwards and forwards, are well versed in
computational techniques, and have read the Collela & Woodward reference, you
can follow these simple steps to multidimensional PPM hydrodynamics...
NOTE:VH-1 in its present form uses the HDF library for writing
data files. If you do not have this library on your system, you will either
have to install it (the source is available here)
or modify prin.f and images.f to use some other form of binary output.
- Read the user's guide to VH-1.
- Download VH-1.
- Copy the code from this link,
or via anonymous ftp at wonka.physics.ncsu.edu.
- Uncompress the file uncompress VH-1_Code.tar.Z
- Untar the file tar xvf VH-1_Code.tar, creating the directory Vhone/
- Compile the code by typing Make within the directory Vhone. You may have to
change some parameters (ie, compiler flags, etc.) within the file Makefile.
PY228 students, in the Makefile you should change the location of the HDF library from
-ldf to /ncsu/py228_info/lib/libdf.a.
- Run the code by typing vhone. This will run the 1D Sod shock tube
problem. Final output is in the file sod1d.1000 in ascii format. Use
your favorite plotting utility to plot the output.
- Try your hand at some more complicated test problems.
- Download the test problem subroutines (TEST)
- Uncompress and untar VH-1_Test.tar.Z as above.
- Choose a particular test problem. The directory Tests contains subroutines
for the test problems in the User's Guide. The test problem to which each file corresponds
is identified by the first letter in the filename. Unmodified subroutines are also
included here in order to restore VH-1 to the original code.
- a_ Strong shock
- b_ Supersonic flow past a ball or cylinder
- h_ Hawley-Zabusky shock tube
- m_ Supersonic flow past a step (Mach step)
- o_ Oblique wall shock
- s_ Sod shock tube in 1, 2, or 3 dimensions
- Copy the relevant subroutines to the Vhone directory,
ie, cp Tests/a_init.f Vhone/init.f
- Modify the declared array sizes in zone.h and sweep.h to
accomodate the desired problem.
- Compile and run the code as described above.
- Compare the output with the data files provided in the
Data tarfile.