Contour advection is a Lagrangian method of simulating the evolution of one or more contours or isolines ofa 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 [1] but may also be active as in,[2] 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.
First we need a set of points that accurately define the contour.These points are advected forward using a trajectoryintegration 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 pointswithin a certain interval.A better method is to use curvature since fewer points are required forthe same level of precision.The curvature of a two-dimensional, Cartesian curve is given as:
| 1 | |
| r |
=\sqrt{\left(
| d2x | |
| ds2 |
\right)2+\left(
| d2y | |
| ds2 |
\right)2}
where
r
s
r\Deltas
\Deltas
In,[3] cubic spline fitting is used both to calculate the curvatureand interpolate new points into the contour.The spline, whichis fitted parametrically,returns a set of second-order derivatives.
A powerful refinement to the technique involves cutting out filaments that have become toonarrow 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 orderto 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 contourare close enough to apply the surgery because of differing spacingin strongly curved versus relatively straight sections.[2]
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.[3]