fluid

NAME

fluid

SYNOPSIS

          fluid [-f filename] [-m mu] [-r rho]

DEFAULTS

          -m 1.0 -r 1.0

DESCRIPTION

fluid specifies the physical fluid parameters of the calculation. fluid interprets the following options.

• -f filename - Initial fluid velocities stored in filename.
• -m mu - Fluid viscosity
• -r rho - Fluid density

FLUID VELOCITY FILE

If the option -f filename is used, then the fluid velocities are assumed to be stored in a file. This file must contain nx*ny lines and each line must contain two numbers corresponding to the components of the fluid velocity at a given point. The data will be read with the following FORTRAN code, so setup the file accordingly.

        open (10, file=filename, status='old')
        do j = 1, ny
          do i = 1, nx
            read(10,*) u(i,j,1), u(i,j,2)
          enddo
        enddo
        close (10)

Notice that the code reads the u component first, then the v component of the velocity field, and it reads along a complete row before the column is changed.

EXAMPLE

          fluid -m 0.6 -r 1.0

The viscosity of the fluid is 0.6 and the density of the fluid is 1.