Index of /rjl/clawpack-4.x/doc/sphinx/users/claw/apps/advection/2d/swirl

CLAWPACK Sample Code

Description

Swirling flow in a box

Edge velocities are stored in aux array (see setaux.f [.html]) with velocity specified by differencing the streamfunction psi.f

The velocities are time-dependent giving reversing flow. These velocities are computed in b4step2.f [.html]. The period is specified by the parameter tperiod in setprob.data.

In theory the solution should agree with the initial data at times t = N*tperiod/2 for all integers N, but because of the numerical diffusion this won't happen.

As a special case, if tperiod = 0, then the velocities are constant in time and b4step2 does nothing. The velocities specified in setaux.f are then used at all times.

An example using AMRCLAW is also included in directory ./amr

This test problem was used in:

@article{rjl:advect,
  author="R. J. LeVeque",
  title="High-resolution conservative algorithms for advection in
  incompressible flow",
  journal="SIAM J. Numer. Anal.",
  volume="33",
  year="1996",
  pages="627-665"
}

See faculty.washington.edu/rjl/pubs/hiresadv

See Instructions

Plots of results

Fortran files

Makefile [.html]

Determines which version of fortran files are used when compiling the code with make.

qinit.f95 [.html]

Sets initial conditions.

setprob.f [.html]

A routine by this name is called by the library routine clawpack/1d/lib/claw1ez.f [.html] and is generally used to set any values needed for the specific problem being solved. In this example, the value of the advection velocity $u$ and a parameter $\beta$ used in setting the initial conditions are read from the file setprob.data [.html].

setaux.f [.html]

Set values in the aux arrays.

psi.f [.html]

Stream function for advection velocities.

b4step2.f [.html]

Called before each time step. Used here to make advection velocity time dependent.

Riemann solvers:

rpn2ad1.f [.html]

Normal Riemann solver (normal to cell interface).

rpt2ad1.f [.html]

Transverse Riemann solver.

Data files

    $ make .output