advection.h

Go to the documentation of this file.

/** @file advection.h
 */
/**
# An advection solver
 
We wish to solve the advection equations
\f[
\partial_tf_i+\mathbf{u}\cdot\nabla f_i=0
\f]
where \f$\mathbf{u}\f$ is the velocity field and \f$f_i\f$ are a list of
passive tracers.  This can be done with a flux-based advection scheme
such as the 2nd-order, unsplit, upwind scheme of [Bell-Collela-Glaz,
1989](references.bib#bell89).
 
The main time loop is defined in [run.h](). A stable timestep needs to
respect the [CFL
condition](http://en.wikipedia.org/wiki/Courant%E2%80%93Friedrichs%E2%80%93Lewy_condition). */
 
#include "run.h"
#include "timestep.h"
 
/**
We allocate the (face) velocity field. For compatibility with the
other solvers, we allocate it as `uf` and define an alias. The
`gradient` function is used to set the type of slope-limiting
required. The default is to not use any limiting (i.e. a purely
centered slope estimation). */
 
vector uf[];
#define u uf
double (* gradient) (double, double, double) = NULL;
 
/**
Here we set the gradient functions for each tracer (as defined in the
user-provided `tracers` list). */
 
extern scalar * tracers;
 
/** @brief Event: defaults (i = 0) */
void event_defaults (void) {
  for (int _f = 0; _f < 1; _f++) /* scalar in tracers */
    f.gradient = gradient;
}
 
/**
User initialisation happens here. */
 
event init (i = 0);
 
/**
The timestep is set using the velocity field and the CFL
criterion. The integration itself is performed in the events of
[tracer.h](). */
 
/** @brief Event: velocity (i++,last) */
void event_velocity (void) {
  dt = dtnext (timestep (u, DT));
}
 
#include "tracer.h"