Contour advection

Contour advection is a Lagrangian method
of simulating the evolution of one or more contours or isolines of
a tracer as it is stirred by a moving fluid.
Consider a blob of dye injected into a river or stream: to first order it could be modelled by tracking only the motion of its outlines.
It is an excellent method for studying chaotic mixing:
even when advected by smooth or finitely-resolved velocity fields,
through a continuous process of stretching and folding,
these contours often develop into intricate fractals.
The tracer is typically passive as in

but may also be active as in,
representing a dynamical property of the fluid such as vorticity.
At present, advection of contours is limited to two dimensions,
but generalizations to three dimensions are possible.

Method

First we need a set of points that accurately define the contour.
These points are advected forward using a trajectory
integration technique.
To maintain its integrity,
points must be added to or removed from the curve
at regular intervals based on some criterion or metric.
The most obvious criterion is to maintain the distance between adjacent points
within a certain interval.
A better method is to use curvature since fewer points are required for
the same level of precision.
The curvature of a two-dimensional, Cartesian curve is given as:

:$\frac = \sqrt$

where $r$ is the radius of curvature
and $s$ is the path.
We need to keep the fraction of arc traced out between two adjacent points,
$r \Delta s$, where $\Delta s$ is the path difference between them,
roughly constant

In,{{Cite journal
, author = Peter Mills
, title = Isoline retrieval: An optimal method for validation of advected contours
, journal = Computers & Geosciences
, volume = 35
, number = 11
, pages = 2020–2031
, year = 2009
, doi = 10.1016/j.cageo.2008.12.015
, url = http://peteysoft.users.sourceforge.net/Mills2009.pdf
, arxiv = 1202.5659
, s2cid = 1637949

cubic spline fitting is used both to calculate the curvature
and interpolate new points into the contour.
The spline, which
is fitted parametrically,
returns a set of second-order derivatives.

Surgery

A powerful refinement to the technique involves cutting out filaments that have become too
narrow to be significant.
If the distance method of adding/removing points is used,
then it is relatively straight forward
to check the distances between all combinations of points.
If a distance between non-adjacent points is too small,
then the two points are separated from their neighbours,
joined together and their neighbours joined also.
Points may then be removed if necessary.
Once we allow surgery, we allow multiply connected domains inside the same contour.
A piece of the contour only one point in length would be removed from the simulation.
The most challenging part of the exercise is keeping track of all the points in order
to reduce the number of distance calculations---see nearest neighbour search.
If the curvature method is used,
then it may be difficult to recognize when two sections of the contour
are close enough to apply the surgery because of differing spacing
in strongly curved versus relatively straight sections.

Validation

Advected contours, e.g. of trace gases (such as ozone) in the stratosphere,
can be validated with satellite remote sensing instruments using a method called isoline retrieval.

External links

ctraj
A library for Lagrangian advection simulations.

References

Fluid dynamics
Meteorological concepts
Numerical climate and weather models