CFD codes list - free software

[pub] UIFlow
[pub] Foil 1.0.2
[pub] PPM -- Piecewise Parabolic Method (DE and LR versions)
[pub] General Relativity NCSA group (black hole evolution)
[pub] NACA airfoils
[pub] College of Marine Studies (sci.geo.fluids models)
[pub] Ocean models (MOM, POM, POP, ...)
[pub] CLAWPACK (library of codes by Randall LeVeque)
[pub] HEATING (a multdimensional heating conduction code)
[pub] PHI3D (3D FEM Navier-Stokes code)
[pub] HYDRO (2D Lagrangian hydrodynamics code)
[pub] Smoothed Particle Hydrodynamics
[pub] Cullimore and Ring Technologies, Inc. (SINDA/FLUINT, SINAPS)
[pub] STSWM (Spectral Transform Shallow Water Model)
[pub] FCT (Flux-Corrected Transport)
[pub] NSC2KE (N-S finite-volume Galerkin 2D code)
[pub] FEMLAB (2D FEM with automatic error control)
[pub] Riemann Problem Package
[pub] CHAMMP repository (s-w equations in spherical geometry)
[pub] NACHOS (NACHOS - FEM for incompressible flows)
[pub] FMS-1D (Fluid Modeling System)
[pub] NPARC Flow Solver
[pub] QUICK 'n SIMPLE (2-D, Macintosh)
[pub] ViewProf (inviscid 2-d multiblock Oellers method panel solver)
[pub] Numerical Models at IMCS (ocean circulation/bottom boundary layer models)
[pub] Synergium's CFD Tools (free version of AVIA)
[pub] ALLSPD-3D Combustor Code
[pub] Unicom Technology Systems (VORSTAB-PC)
[pub] FORTRAN codes for computing the discrete Helmholtz integral operators
[pub] NaSt2D (2D Navier-Stokes solver)
[pub] WinPipeD (pipe hydraulics)
[pub] DROP
[pub] Test Problem Code Archive
[pub] USGS water resources applications software
[pub] The Center for Computational Sciences and Engineering, LBNL
[pub] Pab3d
[pub] CFX User Subroutine Archive
[pub] MicroTunnel
[pub] Mouse
[pub] Mathtools


[index] [up] [down]

link to UIFLOW

From Richard Hilson /
Subject: UIFlow and Post3D
Date: 29 Apr 94 12:23:14 -0800
Article: 177 of sci.physics.computational.fluid-dynamics

I got a copy of the 2D CFD code *UIFlow* from NCSA for the Mac. It seems
to work. For those who are interested, UIFlow can be found in ( in the file Mac/UIFlow.

Foil 1.0.2

[index] [up] [down]

links to Foil: Stanford, France, ICEM CFD Engineering

From Robin B. Lake /
Subject: Airfoil Program for Mac Available
Date: 3 May 1994 00:34:46 GMT
Article: 193 of sci.physics.computational.fluid-dynamics

I happened across a program in the info-mac archives at today.  To my surprise, it turns out to be:
"This is the ReadMe file for Foil 1.0.2, a Macintosh airfoil
generation and display program.  Foil 1.0.2 is System 6/7 happy,
meaning it should run equally well under both systems.  Foil is free,
but not public domain.  I retain all rights, and it may not be
distributed in modified form, or without any of the files.  Since it's
free, please read the file "A Note About the Author".

PPM -- Piecewise Parabolic Method

[index] [up] [down]

Colella, P., Woodward, P.R., 1984, J. Comput. Phys., 54, 174.

PPM is high-order (3/2), shock-capturing Godunov-type method. Before using PPM take a look at some papers, and read infomation about LR version.

Direct Eulerian version

Tomasz Plewa
Date: 10 Jun 94
link to DE version
DE version is generally accessible through WWW ftp link.
Please, let me know if you have no access to WWW server.

Lagrangian Remap version (1)

Greg Lindahl
University of Virginia
Date: 29 Jun 94
link to LR version

Lagrangian Remap version (2)

Virginia Hydrodynamics-1
Theoretical Astrophysics, Dept. of Physics
Date: 03 Apr 95

Lagrangian Remap version (3)

The sPPM benchmark code
Paul Woodward
Laboratory for Computational Science and Engineering, Dept. of Astronomy
Twin Cities Campus, University of Minnesota
Date: 21 Oct 96

General Relativity NCSA group (black hole evolution)

[index] [up] [down]

Relativity Group
National Center for Supercomputing Applications
University of Illinois at Urbana-Champaign

Home Page and summary.

ONEDBH (no source here?)
ONEBH (overview and source)
3D codes front end.

NACA airfoils

[index] [up] [down]

Collection of articles and links related to NACA airfoils is here.

College of Marine Studies (sci.geo.fluids models)

[index] [up] [down]

README file and ftp link.

Ocean models

[index] [up] [down]

Geophysical Fluid Dynamics Laboratory

GFDL's home page and MOM page (also ftp).

Princeton Program in Atmospheric and Oceanic Sciences (AOS)
thanks to Arkady Terzhevik

AOS' home page and POM home page.

The Parallel Ocean Program (POP)

POP's home page, Fluid Dynamics Group (T-3), LANL.

CLAWPACK (library of codes by Randall LeVeque)

[index] [up] [down]

Randall J. LeVeque
Applied Mathematics Dept., FS-20
University of Washington
Seattle, WA 98195

Version 2.0 of clawpack is now available from netlib, in pdes/claw.
See also the homepage at URL
There are several new features and applications, and expanded documentation.

        CLAWPACK  --  Conservation LAWs software PACKage

Version 1.0 of CLAWPACK is now available from netlib.  This is a
package of Fortran subroutines for solving hyperbolic systems of
conservation laws in one and two space dimensions. (Future versions
should include 3D.)  High resolution flux-limiter methods are used,
based on solving one-dimensional Riemann problems.  In two dimensions,
multi-dimensional wave-propagation methods are used.  The current
version supports only uniform Cartesian grids in rectangular domains,
and is intended primarily as a research and teaching tool.  The
modular form should make it easy to modify and experiment with other
methods, as well as to apply it to new problems.

The code requires that the user supply a Riemann solver for the
problem being solved and also a subroutine that implements the
boundary conditions.  Source terms can also be handled (via Strang
splitting) in which case an ODE solver for the source terms must also
be supplied.  Several examples of different Riemann solvers are
included with the package, including e.g., advection equation,
Burgers' equation, Euler equations, isothermal equations, shallow
water equations.  Various different boundary conditions are also
demonstrated, such as periodic, inflow, nonreflecting, and solid
walls.  Numerous example drivers are included to demonstrate the use
of these subroutines.  Matlab m-files are included to graphically
display the output, or they can be used as a model for writing
graphics routines in other languages.  Documentation is included in
the package, including the postscript file for an introductory paper
on the package.

The package can be obtained from netlib, where it resides in the library
pdes/claw.  This directory currently contains:
   index      the index
   clawpack   the basic package, with 1d, 2d routines and many examples
   doc        documentation (postscript files of papers, slides)
   advection  applications of clawpack to 2d advection equations, with 
	      examples from a recent paper on this subject (in doc).
   nozzle     application of clawpack to the quasi-1d nozzle problem.
   vBurgers   viscous Burgers' equation, with diffusion equation as source.
   L-domain   L-shaped domains, such as the forward facing step problem.

My hope is that a library of more sophisticated applications will be
gradually built up, and contributions from other users will be
gratefully accepted.
The package may be obtained by anonymous ftp from, where it resides in netlib/pdes/claw. The file netlib/pdes/claw/clawpack.shar is a tar file of the entire directory. Files can also be obtained by sending e-mail to Send the message "help" to this address for more information. Alternatively, you can obtain files by ftp through Mosaic using the URL . It is also possible to browse through the library using Mosaic from .

HEATING (a multdimensional heating conduction code)

[index] [up] [down]

Kenneth. W. Childs /
Heat Transfer and Fluid Flow Group
Oak Ridge National Laboratory

HEATING is a general-purpose, conduction, heat transfer program
written and maintained by the Heat Transfer and Fluid Flow Group
More about HEATING.

From Jonas T. Holdeman, Jr. /
Date: Tue, 07 Feb 1995 15:09:52 -0500

I understand both workstation and PC versions of the HEATING code are 
available freely from the authors, Gary Giles ( and Ken Childs 
( [...] Please check with the authors for latest 

From Shelly Spencer /
Date: Tue, 06 Aug 1996 13:07:58 -0400

On your net page regarding free software, the following entry is free
only if your funding comes from the US Dept of Energy, although they
don't tell you this until they ask you to fill out and return several
pages of forms!  The fee is $400 for the HEATING code.

I just wanted to let you know so other folks don't lose valuable time
filling out the forms for a code that is not free to most people.

Date: Thu, 9 Jan 1997 17:24:44 -0500

HEATING was developed for the U.S. Department of Energy (DOE). DOE
requires that their codes be distributed through one on their code
centers. In the United States HEATING is distributed by the Radiation
Safety Information Computational Center (RSICC) in Oak Ridge, Tennessee.
RSICC can be contacted via e-mail at PDC@ORNL.GOV or by phone at
423-574-6176. RSICC has a WWW page at
HEATING is code number PSR-199 in the RSICC code collection.
RSICC charges a fee to cover their cost of distributing codes. This fee
varies based on the affiliation of the requester. I don't know how much
the fee currently is.
Documentation for HEATING is available on line as a pdf (Portable Document
Format) file at

PHI3D (3D FEM Navier-Stokes code)

[index] [up] [down]

Paul T. Williams /
Heat Transfer and Fluid Flow Group
Oak Ridge National Laboratory

As a part of validation studies for the locally developed finte-element
incompressible Navier-Stokes Code PHI3D [...]
More about PHI3D.

From Jonas T. Holdeman, Jr. /
Date: Tue, 07 Feb 1995 15:09:52 -0500

PHI3D was developed by Paul Williams ( for his Ph.D.
dissertation. He says it needs to be cleaned up before he will release
it.  If someone were willing to give him some funding, this process
could be speeded up, as he has no funding to work on the code at

HYDRO (Two-dimensional Lagrangian hydrodynamics code)

[index] [up] [down]

From Hao Chen /

ftp link

Los Alamos National Laboratory
Benchmarking Section of Group C-3
Dennis V. Brockway, Fred Gama-Lobo, Karl Wallick, John Romero, & Steve W. White
May 1, 1984


ABSTRACT:   This is a two-dimensional Lagrangian hydrodynamics code based on
            an algorithm of W.D. Schultz.
NUMERICAL:  According to the report below, HYDRO is representative of a large
            class of codes used in the Los Alamos Laboratory and is 100% 

NOTES:      HYDRO is very similar to SIMPLE - Almost identical methods.
            However, HYDRO can handle General Geometry, while SIMPLE uses Fixed
            Geometry (Equation of state table 2X2) 

AUTHORS:    produced under U.S. government contract (w-7405-eng-36)
            by Los Alamos National Laboratory
            by D. Brockway, F. Gama-Lobo, K. Wallick, J. Romero, S. White
            Contact benchmarking section of group c-3 LANL (505)667-7028
REFERENCES: Report: Los Alamos Laboratory Computer Benchmarking 1988
            H.J. Wasserman
            William D. Schulz
            Methods of Computational Physics
            Volume 3, 1964

STSWM (Spectral Transform Shallow Water Model)

[index] [up] [down]

From Hao Chen / ftp link


Copyright (C) 1992
University Corporation for Atmospheric Research
All Rights Reserved

       Ruediger Jakob
National Center for Atmospheric Research
    Boulder, CO 80307-3000

         August 1992


1. Software Distribution Conditions
2. Description of Software
3. Directory of Files
4. Corrections and Changes

Another version of this code (or postprocessor?) is called REFSOL. REFSOL uses netcdf package.

PSTSWM: a message-passing version is also available.

NA Digest   Sunday, May 28, 1995   Volume 95 : Issue 22

From: Pat Worley 
Date: Wed, 24 May 1995 12:48:20 -0400
Subject: PICL and PSTSWM software updates



PSTSWM v4.0 is now available from
(It has been available for awhile, but the user guide has just been

PSTSWM is a message-passing benchmark code and parallel algorithm
testbed that solves the nonlinear shallow water equations on a
rotating sphere using the spectral transform method. PSTSWM was
developed by to evaluate parallel algorithms for the spectral
transform method as it is used in global atmospheric circulation
models. Multiple parallel algorithms are embedded in the code and can
be selected at run-time, as can the problem size, number of
processors, and data decomposition. This flexibility allows the code
to be tuned on a parallel platform before benchmarking, thus
evaluating the multiprocessor on its ability to solve the numerical
problem rather than it ability to execute a given fixed parallel
implementation. Six different problem test cases are also supported,
each with associated reference solutions and solution and error
analysis options.

As of 4/1/95, PSTSWM runs on the Cray Research T3D, the IBM SP-1 and
SP-2, the Intel iPSC/2, iPSC/860, DELTA, and Paragon (on both GP and
MP nodes and using either the NX or SUNMOS operating systems), the
nCUBE/2 and nCUBE/2S, across a network of SUN and IBM workstations,
and on a Cray vector machine (as a serial application). Message
passing is implemented using MPI, PICL, PVM, and/or native message
passing libraries, with the choice being made at compile time.  In
principle, it should also run on any other platform on which MPI,
PICL, or PVM is available.

To aid in tuning and in understanding the parallel performance, PSTSWM
has been instrumented for the collection of performance data using the
PICL trace and profile collection interface. The PICL implementation
of the code must be used in order to collect performance data on
interprocessor communication but a mixed PICL/native implementation is
also provided that can be used to collect data on events not related
to message passing. In the mixed implementation, the performance
sensitive message passing uses native commands and PICL is only used
in the collection of the performance data.

For more information or for alternative ways of acquiring source code,

  Pat Worley

FCT (Flux-Corrected Transport)

[index] [up] [down]

link to FCT page


Flux-Corrected Transport (FCT) is a conservative, monotone technique
for integrating generalized continuity and hydromagnetic
equations. [...]

FCT was the first of the class of high-order, monotone schemes for
solving generalized continuity equations (e.g., the equations of
Eulerian hydrodynamics). An FCT bibliography of basic papers is posted

A suite of FCT modules with test programs is being developed for
posting to the HPCC/ESS Software Exchange. [...]

The software for solving 2D hydrodynamical problems will be found in
the package LCPFCT2. For 2D magnetohydrodynamical (MHD) applications,
the package MHDFCT2 also will be needed. Similar software packages
LCPFCT3 and MHDFCT3 eventually will be posted for solving 3D problems.


C. Richard DeVore 
Laboratory for Computational Physics & Fluid Dynamics, Naval Research

NSC2KE (N-S finite-volume Galerkin 2D code)

[index] [up] [down]

Bijan Mohammadi
Domaine de Voluceau
BP 105
78153 Le Chesnay, Rocquencourt

      2D and AXI Euler and Navier-Stokes equations solver 
   o  explicit multi-steps time integration process
   o  upwind schemes and linear interpolation method for 
      the computation of the convective fluxes using a finite 
      volume formulation.
   o  classical central galerkin p1-finite element method 
      for the computation of the diffusive fluxes
   o  k-epsilon turbulence model with two-layer approach or wall laws
      Copyright(C) 1994 Bijan Mohammadi-Stephane Lanteri

Abstract and user guide.


[index] [up] [down]

Public ftp archive of the CFD Community Club at the Rutherford Appleton Laboratory.


[index] [up] [down]

FEMLAB group
Dep. of Mathematics
Chalmers University of Technology

FEMLAB is an interactive program for the numerical solution
of ordinary and partial differential equations based on the
Finite Element Method in adaptive form with automatic error
Femlab is distributed freely to anyone interested. You may read more about FEMLAB and find what do you need to get FEMLAB running. Finally, you may download FEMLAB.

Riemann Problem Package

[index] [up] [down]

Eli L. Isaacson, Dan Marchesin, and Bradley J. Plohr
SUNY Department of Applied Mathematics and Statistics

The executable file "rp" contained in the present directory has been compiled on a SPARCstation 2. It requires OpenWindows 2.0 and the Sun Fortran library to work.

The program "rp" is an interactive graphics program for
exploring systems of two quadratic conservation laws.
It includes the ability to construct:

        integral curves
        Hugoniot curves
        wave curves (comprising rarefaction, shock,
                and composite segments)
        wave speed diagrams
        bifurcation loci (including coincidence, inflection,
                secondary bifurcation, and double sonic loci)
        solutions of Riemann problems

NOTICE:  You are free to use this program for studying quadratic models
and formulating conjectures about their solutions.  If this program
provides useful insight, please acknowledge its use in any relevant


CHAMMP repository

[index] [up] [down]

Access via ftp. Here is an excerpt from README file:

Downloaded from; 2/24/94

#  chammp
#  ======
#  This directory contains items relating to the numerical solution 
#  of the shallow water equations in spherical geometry.
#  The shallow water equations are used as a kernel for both
#  oceanic and atomospheric general circulation models and
#  are of interest in evaluating numerical methods for weather
#  forecasting and climate modeling.  The DOE Computer Hardware, 
#  Advanced Mathematics and Model Physics (CHAMMP) program 
#  is interested in the development of new mathematical methods
#  for these problems.  To promote this developmet, a set of 
#  test cases has been proposed and example software and reference
#  solutions are provided.

NACHOS - FEM for incompressible flows

[index] [up] [down]

From Jonas T. Holdeman, Jr. /
Date: Tue, 07 Feb 1995 15:09:52 -0500

I am surprised you don't have any mention of the NACHOS II code by D. K. 
Gartling of Sandia National Laboratory, Albuquerque, New Mexico 97185, USA. 
I had no trouble obtaining a free copy. The code is documented in two 
reports: "NACHOS II - A Finite Element Computer Program for Incompressible 
Flow Problems, Part I - Theoretical Background (SAND86-1816 UC-32) and Part 
II - User's Manual (SAND86-1817 UC-32). The code is written in Fortran. It 
features triangular and quadrilateral elements, primative variables, with 
post-processing for stream-function, fluid stress and flux computation. 

I heard that Gartling was working on a 3D version of Nachos, but I don't 
know for sure. There should be several other codes from Sandia. 

FMS-1D (Fluid Modeling System)

[index] [up] [down]

Barry Merriman / bmerriman@UCSD.EDU
UCSD Fusion Energy Research Center
UCLA Dept. of Math

FMS (Fluid Modeling System) is a system for constructing and
numerically  solving  the  kinds  of equations that arise in
general  fluid  flow  problems, i.e. systems of conservation
laws.  The  current  version  handles  flows  in  one  space
dimension,  plus  time. FMS is intended to be a flexible and
powerful  general  purpose  tool  for this type of modeling.

The  user  supplies  FMS  with  a description of the spatial
domain,  and the system of equations to be solved, including
any  boundary  conditions  or  initial  conditions,  and FMS
attempts  to  compute the the corresponding solution (steady
state or time dependant).

NPARC Flow Solver

[index] [up] [down]

The WIND code is a structured, multi-zone, compressible flow solver with flexible turbulence and chemistry models. It is being developed jointly by the NPARC Alliance (the Air Force Arnold Engineering Development Center and NASA Glenn Research Center) and the Boeing Company.

Please see for more information.


[index] [up] [down]

From Gerald Recktenwald /
PSU Mechanical Engineering
P.O. Box 751
Portland, OR 97207
(503) 725-4296

We have developed a free, 2-D CFD program for the Macintosh
called QUICK 'n SIMPLE, or "QnS" for short.

The program combines a graphical user interface, for specification
of the problem, with an iterative control-volume finite-difference
solution algorithm to solve two-dimensional convection with
diffusion of a passive scalar.  It uses either the SIMPLE or
SIMPLER coupling algorithms and allows the user to select from
three popular upwinding schemes: pure upwinding, power-law
differencing scheme (PLDS), and the quadratic upstream
interpolation for convective kinematics (QUICK) scheme.

As the name implies, the program is designed to be both quick and
simple to use.

You can check it out, and download a copy if you like, from our web
server at
If you don't have a web browser, then you can also get a
copy via anonymous ftp site at,
but its not nearly as fun, and you don't get to see our cool
examples and tutorial.

We have three versions available: 1) 68K, 2) 68K w/ fpu, and 3)
PowerMac.  The only limit on the size of the problem you can solve
is the memory of your Mac.  The 68k version (without fpu) can't
solve large problems in reasonable amounts of time, but QnS
really sings on a PowerMac. The PowerMac version runs over 200
times faster (not a typo) than the 68K version, so a problem
which may take 30 seconds on a PowerMac would take nearly two hours
on a 68K machine.  The fpu version is somewhere in between, but
closer in speed to the 68K version.

Let us know what you think.

Scott J. Forbes
Department of Mechanical Engineering
University of Akron, Akron, OH  44325-3903

Gerald W. Recktenwald
Department of Mechanical Engineering
Portland State University, Portland, OR  97207-0751


[index] [up] [down]

From Matti E Hy|tyniemi /
Helsinki University of Technology
Date: Tue, 26 Sep 1995 19:59:01 +0200

ViewProf is inviscid 2-d multiblock Oellers method panel solver for
Windows 3.1.

It is public domain and available via Web at,
or by FTP at

Numerical models at IMCS

[index] [up] [down]

Kate Hedstrom /
Institute of Marine and Coastal Sciences
Rutgers, The State University of New Jersey

For more information about SPEM/SCRUM/SEOM/BBL see

Synergium's CFD Tools

[index] [up] [down]

From Alex Portalatin /
Date: Tue, 14 Nov 95 12:37:37 -0800

We have now made available a FREE evaluation version of the AVIA system 
at our site

See also: Synergium's CFD Tools (commercial).

ALLSPD-3D Combustor Code

[index] [up] [down]

From David Fricker
Date: 10 Feb 1996 07:31:57 GMT
Organization: NASA Lewis Research Center, Cleveland, OH Article: 3966 of sci.physics.computational.fluid-dynamics

The ALLSPD Team at NASA Lewis Research Center is pleased to announce
their new web page for the ALLSPD-3D CFD code for gas turbine

The ALLSPD-3D Combustor Code is a numerical tool developed by the
Internal Fluid Mechanics Division at the NASA Lewis Research Center
for simulating chemically reacting flows in aerospace propulsion
systems.  It provides the designer of advanced engines an analysis
tool that employs state-of-the-art computational technology. The code
can simulate multi-phase, swirling flows over a wide Mach-number range
in combustors of complex geometry.

This site provides details on ALLSPD-3D, the user manual, and
procedures for obtaining the ALLSPD-3D code.  (NOTE: Distribution of
documentation is open to anyone, but distribution of the code is
restricted to U.S. residents.)

Check it out at:

Unicom Technology Systems

[index] [up] [down]

From Alan Hsing /
Unicom Technology Systems
2721 Bishop Street
Lawrence, KS 66046

For anyone interested in a handy aerodynamic design code, please check
out VORSTAB-PC 3.1.   You can ftp to '', and

   - login as 'ftp',
   - use your email id as password,
   - change directory to /pub/MSDOS_UPLOADS/aero,
   - get VORSTB31.EXE in binary mode and VORSTB31.TXT in ascii mode,
   - read VORSTB31.TXT for instructions on installation.

Email me if you have any problems doing that.  Note: the ftp site may get
highly loaded during busy hours. In that case, just be patient and try it

The homepage for VORSTAB-PC 3.1 has now been setup.  You can download
the free demo from

Acoustic Boundary Elements with Fortran source codes

[index] [up] [down]

From S.M.Kirkup /
Integrated Sound Software
25 Smithwell Lane
Hebden Bridge
West Yorkshire HX7 7NX

Tel: 01422 844607

These are Fortran subroutines for the evaluation of the discrete form
of the Helmholtz integral operators for 2D, 3D or 3D axisymmetric
problems.  The subroutines are useful for the solution of the
Helmholtz equation by boundary element and related methods.

The codes are available free from the author.


[index] [up] [down]

From Helge Tielboerger /
Lehrstuhl A fuer Thermodynamik
Technische Universitaet Muenchen
Arcisstrasse 21
80333 Muenchen

From Tilman Neunhoeffer /
Institut fuer Informatik
TU Muenchen
Gabelsbergerstr. 39, 310,
D-80290 Muenchen
Tel 089/2105-5375

Here is the place where you can get a public domain solver for the

It is a rather simple Finite-Volume-code but it can solve free-surface-flows.
Additionally there is a book available which describes step by step the
developement of the code. 
So the code together with the book is an excellent
introduction into CFD and a good basis to develop more enhanced code.
The code is available via ftp-server:
Solutions of some problems are available from

The ReadMe-file from that directory:

|                                                                             |
|              The 2D CFD Program NaSt2D associated to the book               |
|                                                                             |
|            "Numerische Simulation in der Stroemungsmechanik"                |
|            by M. Griebel, Th. Dornseifer, and T. Neunhoeffer                |
|                              Vieweg 1995                                    |
|                                                                             |
|              The english translation is published by SIAM                   |
|               Title: "Numerical Simulation in Fluid Dynamics"               |
|                                                                             |


- The program is a 2D solver for the incompressible, transient Navier-Stokes
  equations including the temperature equation and free boundary problems.
  It uses finite differences for discretization on a structured equidistant
  staggered grid, central and upwind (donor-cell) discretization of the 
  convective parts and an explicit time stepping scheme (Chorin's projection 
  method). The free boundary value problems are treated with the MAC technique.

- The program was developed for education at a Computer Science Institute.
  Thus, it is easy to understand and easy to implement but it is not a
  "state of the art" program. 

- The files are compressed in
  Please uncompress with


  Then, the files are located in the directory PRG.

  Furthermore, the files are available in the directory UNCOMPRESSED.

- The program is written in the C programming language.

- Compilation of the program with "make".

- Used C-compiler is "gcc" (see "makefile")

- Run the program by
        "run "    e.g."run dcavity.par".

- In the inputfiles "*.par" (examples can be found in the sub-directory
  "PARAMETER"), the parameters for the program are given, such as
  the parameters for the discretization, for the fluid, for the
  initial data,... The inputfile given when starting "run" is read
  at the beginning of the program.

- The following problems are implemented and inputfiles are available
  in the subdirectory "PARAMETER"
  * dcavity:            Driven Cavity
  * backstep:           Flow over a backward facing step
  * balken:             Flow past an inclined plate
  * kreis:              Flow past a circular obstacle
  * klappen:            Flow through a T-junction
  * damm:               Breaking dam (free boundary value problem FBVP)
  * tropfen:            Splash of a liquid drop (FBVP)
  * spritzguss:         Injection moulding (FBVP)
  * wafa:               Flow over a backward facing step with free surface
  * convection:         Buoyancy flow with heated side walls
  * fluidtrap:          Buoyancy flow with obstacles
  * rayleigh:           Buoyancy flow with heated upper and lower wall
                        (Rayleigh-Benard Flow)

- For some problems (e.g. backstep, balken, kreis) the pressure correction
  does not converge in the first time steps. This is due to the fact that 
  the initial velocity field does not satisfy the continuity equation.

- It is possible to store intermediate data in "outfile" specified in the 
  inputfile. Later on, you can use this intermediate data as initial data
  for continuing your simulation if you specify this file as input-data 

- The comments in the code are not very detailed. For a complete
  description of the code, please see the book.

- For the implementation of new problems, it is necessary to write
  a new inputfile, to modify the boundary data in the function
  "SETSPECBCOND" in the file "boundary.c" and the definition of the
  obstacles in "INIT_FLAG" ("init.c").

- Parallelization and turbulence simulation described in the book are
  not included in this package.

- For visualization, we used the graphics tool IDL. Using other
  visualization software, it might be necessary to adapt the functions 
  "OUTPUTVEC_bin" and "streaklines_bin" in visual.c.


ATTENTION: The program was developed on a HP workstation. 
           There might appear some problems under DOS or WINDOWS95


If you have questions, please contact

Stephan Knapek
Universität Bonn
Institut für Angewandte Mathematik
Abteilung Wissenschaftliches Rechnen und Numerische Simulation
Wegelerstr. 6
 D-53115 Bonn


[index] [up] [down]

From David C. Shaw /
Subject: Re: CFD Software for Educational Purposes
Date: Fri, 26 Jul 1996 11:38:07
Organization: NeoSoft, Inc.
Article: 5027 of sci.physics.computational.fluid-dynamics

WinPipeD (b), the basic version, is available for non-commercial use,
including education, via Internet.  Please visit us at:

The source is also available, which makes a great research tool.


[index] [up] [down]

From U. Fritsching /
Subject: CFD codes list, new PD software
Date: Fri, 15 Nov 1996 15:09:49 +0100

The program DROP describes the behaviour of individual spherical
particles like drops, bubbles or solid particles in multiphase flow
situation of a moving fluid. Starting with a force balance on an
individual particle, taking into account the relevant fluid-particle
interactions, the equation of particle movement in a Lagrangian
reference frame is derived. The resulting ODE-system is solved based
on a Runge-Kutta algorithm. As a result the velocity versus time or
position versus time function can be plotted or data can be written
into a file.  Boundary conditions regarding the properties of the
particle and the fluid can be changed to describe different areas of
applications. The program is dialog driven and is running on PC´s
with Windows 3.1. The new version DROP V1.1 allows for analysis of
accelerating fluid flow in two dimensions.

The software can be downloaded from:

Dr.-Ing. U. Fritsching
Institut für Werkstofftechnik
Badgasteiner Str. 3
28359 Bremen

Tel.: +49-421-218-3663
Fax: +49-421-218-5378


[index] [up] [down]

From U. Fritsching /
Subject: CFD codes list, new PD software
Date: Fri, 15 Nov 1996 15:09:49 +0100

The program DROP describes the behaviour of individual spherical
particles like drops, bubbles or solid particles in multiphase flow
situation of a moving fluid. Starting with a force balance on an
individual particle, taking into account the relevant fluid-particle
interactions, the equation of particle movement in a Lagrangian
reference frame is derived. The resulting ODE-system is solved based
on a Runge-Kutta algorithm. As a result the velocity versus time or
position versus time function can be plotted or data can be written
into a file.  Boundary conditions regarding the properties of the
particle and the fluid can be changed to describe different areas of
applications. The program is dialog driven and is running on PC´s
with Windows 3.1. The new version DROP V1.1 allows for analysis of
accelerating fluid flow in two dimensions.

The software can be downloaded from:

Dr.-Ing. U. Fritsching
Institut für Werkstofftechnik
Badgasteiner Str. 3
28359 Bremen

Tel.: +49-421-218-3663
Fax: +49-421-218-5378

Test Problem Code Archive

[index] [up] [down]

From Robert Fisher / bobf@physics.Berkeley.EDU
Subject: Test Codes
Date: Thu, 12 Dec 1996 19:57:23 -0800

   I have written two short programs to calculate the exact Sod shock
tube problem solution and the exact Taylor-Sedov blast wave solution.
Even though these codes are very short and straightforward, I could not
find links to similar programs under the CFD Codes List. When you have
time, could you please provide a link to these codes? You can cut and
paste the following description :

Test Problem Code Archive
sod.f -- A f77 program to calculate the exact solution to the
         shock tube problem. Initial configuration can be easily 
         altered to handle similar shock tube problems.
ts.c  -- A C program to calculate the exact solution to the Taylor-     
                 Sedov blastwave problem.

The site is located at 

USGS Water Resources Applications Software

[index] [up] [down]

From Michele Y Crouse /
Contact: Steve Regan 703-648-5896,

United States Department of the Interior
U.S. Geological Survey
Office of Water Information
Reston, Virginia 20192


A suite of 48 software packages and related materials, used by the U.S.
Geological Survey for hydrologic analysis and modeling, is now available for
electronic retrieval through an on-line repository on the World Wide Web.
The repository is accessible via the WWW from the USGS Water Resources Home
Page at:

by selecting the announcement for "USGS Water Resources Applications 
Software" or directly at:

In addition, the repository is available via anonymous File Transfer
Protocol (FTP) from the USGS Water Resources Information server:

 or (path: pub/software)

Note: a mirror site ( is available 
at Friedrich-Schiller University, Department of Geoinformatics, Jena, 
Germany that some users may find more convenient.  FTP access is also 
available at (path: pub/software).

Each software distribution package consists of compiled or source code,
test data sets and documentation files.  All of the USGS water resources
hydrologic analysis and modeling software applications available publicly on
the WWW are documented by published USGS reports.  The software distribution
packages have been prepared primarily for the Data General AViiON DG/UX
platform and for compilation on other UNIX-based computers.  The USGS 
continues to prepare software packages for its own use on UNIX-based and 
other computer platforms such as DOS-based personal computers.  As these 
packages are available, they will be added to the repository.

The software available in the repository is grouped into the following
categories: geochemical, ground water, surface water, water quality and
general.  Examples of the software now available through the on-line
repository include:

     PHREEQC--A program for aqueous geochemical calculations;

     MOC--Two-dimensional method of characteristics ground-water flow and
          transport model;

     MODFLOW--A modular three-dimensional finite-difference ground-water flow

     SUTRA--Saturated and/or unsaturated, constant or variable-density fluid
            flow, and solute or energy transport (two-dimensional finite 
            element code);

     ANNIE--Interactive hydrologic (time series) data management;

     BRANCH--Branch-network one-dimensional dynamic flow model;

     HSPEXP--Expert system for calibration of the Hydrological Simulation
             Program--Fortran (HSPF).

Since June 1996, when these programs and associated documentation became
publicly available on the WWW, the number of accesses by users outside the
USGS has increased rapidly.  Currently, the site is "visited" about 2500
times per month by users outside the USGS.

For general information on USGS water resources applications software,
please contact:  USGS, Hydrologic Analysis Software Support Team, 437
National Center, Reston, VA 20192 (e-mail:

The Center for Computational Sciences and Engineering, LBNL

[index] [up] [down]

From John B. Bell /
Date: Mon, 03 Feb 1997 10:44:28 -0800

Center for Computational Sciences and Engineering
Lawrence Berkeley National Laboratory

The Center for Computational Sciences and Engineering at the Lawrence
Berkeley National Laboratory has recently released some CFD software
that would be appropriate for the CFD codes list. The software is
available from our web site at:

Included are two applications codes:

Gas Dynamics AMR - HAMR

An adaptive mesh refinement algorithm for solving the compressible
Euler equations in two and three dimensions.  The code uses an
operator-split second-order Godunov algorithm for advancing the flow

Single Grid Incompressible Flow Solver - VarDen

An incompressible Euler and Navier-Stokes algorithm based on a
projection method formulation.  There are two and three dimensional
versions of the algorithm both of which can treat large-amplitude
density variations.

Also available are support software libraries

BoxLib C++ library

C++ class library to support block-structured adaptive algorithms.  The
applications codes discussed above are based on BoxLib

Amrvis visualization tool

A visualization tool for viewing data generated by block-structured
adaptive algorithms.  The above applications generate plot files that
can be visualized with AmrVis.  AmrVis is also built using BoxLib.


[index] [up] [down]

A 3-D Navier-Stokes CFD System

PAB3D solves the three-dimesional Reynolds-averaged Navier-Stokes equations with a finite-volume formulation on structured multi-block grids. Initially developed to predict simple internal and external performance of nozzle-afterbody configurations, the code has evolved to predict complex flow in mixer-ejector nozzles, thrust vectoring, multi-species jet mixing, and propulsion installation with a variety of aerodynamic configurations. The flow solver has a robust (production quality) two-equation k-epsilon turbulence models with anisotropic algebraic Reynolds stress models. PAB3D also has the capability of modeling standard air real-gas effects.

Code Developer: Dr. Khaled S. Abdol-Hamid (Analytical Services and Materials)

Availible by request from NASA Langley Research Center See

CFX User Subroutine Archive

[index] [up] [down]

I have recently created an archive of user fortran subroutines for CFX users. The site is located at

Currently the site contains subroutines that I wrote for my thesis research at the University of Illinois. I am hoping to get submissions from other CFX users, and make the site a useful resource.

The motivation for this site was the difficulty that I had in learing to write subroutines for CFX - there are very few useful examples in the documentation. Hopefully the site will prevent others from having the same problem

David Creech
Department of Mechanical Engineering
University of Illinois, Urbana-Champaign


[index] [up] [down]

MicroTunnel analyzes 2-D subsonic, transonic, and supersonic flow (Mach 0.1 to 3.0) over multi-element configurations on a raster type Cartesian mesh. Flow visualization is provided through color contour plots of local Mach number, density, pressure, and temperature, including streamline plots. Shock waves are easily captured!

MicroTunnel runs out of Microsoft Windows through simple "point and click" operations. While the program is computing, other applications can still be run with minimal sacrifice in performance. Any screen output can be sent to a color printer.

MicroTunnel is based on a finite volume description of the unsteady Euler equations using Van Leer type flux vector splitting. A solution is marched in time from an instantaneous acceleration of the flow. Flow resolutions of 200x150 take approximately six hours to converge to steady state on a system with minimum requirements.

For more information see


[index] [up] [down]

This is the very first public version of MOUSE, a finite volume library. At the Institute of Combustion and Gas Dynamics at the University of Duisburg in Germany, we use this library for CFD applications. Since this is the first public version and it has been restructured quite a bit from the last internal one, it is likely that you discover some bugs if you use it. Nevertheless we think it is in a state where it can be useful. If you discover bugs, please do not hesitate to contact us.

MOUSE is an object oriented framework for finite volume computations on unstructured grids. It serves as a development library, as well as a ready to use application. For more information see


[index] [up] [down] - Free scientific portal for MATLAB/MIDEVA m-files and toolboxes, and Excel/Java/Fortran/C++ resources and links. - Complementary products for MATLAB, like MIDEVA (fast MATLAB replacement), MATCOM (Compiler for MATLAB), Visual MATCOM (integrate m-files into Visual C++) and others, all available for download.

Wayne Christopher /