Contour advection

Last updated
Contour advection simulation Wv cadv ecmwf small.gif
Contour advection simulation

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 [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.

Contents

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:

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

In, [3] 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. [2]

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. [3]

Related Research Articles

A centripetal force is a force that makes a body follow a curved path. The direction of the centripetal force is always orthogonal to the motion of the body and towards the fixed point of the instantaneous center of curvature of the path. Isaac Newton described it as "a force by which bodies are drawn or impelled, or in any way tend, towards a point as to a centre". In the theory of Newtonian mechanics, gravity provides the centripetal force causing astronomical orbits.

<span class="mw-page-title-main">Horizon</span> Apparent curve that separates earth from sky

The horizon is the apparent curve that separates the surface of a celestial body from its sky when viewed from the perspective of an observer on or near the surface of the relevant body. This curve divides all viewing directions based on whether it intersects the relevant body's surface or not.

<span class="mw-page-title-main">Gaussian curvature</span> Product of the principal curvatures of a surface

In differential geometry, the Gaussian curvature or Gauss curvatureΚ of a smooth surface in three-dimensional space at a point is the product of the principal curvatures, κ1 and κ2, at the given point:

<span class="mw-page-title-main">Non-uniform rational B-spline</span> Method of representing curves and surfaces in computer graphics

Non-uniform rational basis spline (NURBS) is a mathematical model using basis splines (B-splines) that is commonly used in computer graphics for representing curves and surfaces. It offers great flexibility and precision for handling both analytic and modeled shapes. It is a type of curve modeling, as opposed to polygonal modeling or digital sculpting. NURBS curves are commonly used in computer-aided design (CAD), manufacturing (CAM), and engineering (CAE). They are part of numerous industry-wide standards, such as IGES, STEP, ACIS, and PHIGS. Tools for creating and editing NURBS surfaces are found in various 3D graphics and animation software packages.

<span class="mw-page-title-main">Scalar potential</span> When potential energy difference depends only on displacement

In mathematical physics, scalar potential, simply stated, describes the situation where the difference in the potential energies of an object in two different positions depends only on the positions, not upon the path taken by the object in traveling from one position to the other. It is a scalar field in three-space: a directionless value (scalar) that depends only on its location. A familiar example is potential energy due to gravity.

<span class="mw-page-title-main">Evolute</span> Centers of curvature of a curve

In the differential geometry of curves, the evolute of a curve is the locus of all its centers of curvature. That is to say that when the center of curvature of each point on a curve is drawn, the resultant shape will be the evolute of that curve. The evolute of a circle is therefore a single point at its center. Equivalently, an evolute is the envelope of the normals to a curve.

<span class="mw-page-title-main">Einstein–Hilbert action</span>

The Einstein–Hilbert action in general relativity is the action that yields the Einstein field equations through the stationary-action principle. With the (− + + +) metric signature, the gravitational part of the action is given as

In physics, the Brans–Dicke theory of gravitation is a competitor to Einstein's general theory of relativity. It is an example of a scalar–tensor theory, a gravitational theory in which the gravitational interaction is mediated by a scalar field as well as the tensor field of general relativity. The gravitational constant is not presumed to be constant but instead is replaced by a scalar field which can vary from place to place and with time.

<span class="mw-page-title-main">Level-set method</span>

Level-set methods (LSM) are a conceptual framework for using level sets as a tool for numerical analysis of surfaces and shapes. The advantage of the level-set model is that one can perform numerical computations involving curves and surfaces on a fixed Cartesian grid without having to parameterize these objects. Also, the level-set method makes it very easy to follow shapes that change topology, for example, when a shape splits in two, develops holes, or the reverse of these operations. All these make the level-set method a great tool for modeling time-varying objects, like inflation of an airbag, or a drop of oil floating in water.

The Clausius–Clapeyron relation, named after Rudolf Clausius and Benoît Paul Émile Clapeyron, specifies the temperature dependence of pressure, most importantly vapor pressure, at a discontinuous phase transition between two phases of matter of a single constituent. Its relevance to meteorology and climatology is the increase of the water-holding capacity of the atmosphere by about 7% for every 1 °C (1.8 °F) rise in temperature.

<span class="mw-page-title-main">Born coordinates</span> Coordinates to capture characteristics of rotating frames of reference

In relativistic physics, the Born coordinate chart is a coordinate chart for Minkowski spacetime, the flat spacetime of special relativity. It is often used to analyze the physical experience of observers who ride on a ring or disk rigidly rotating at relativistic speeds, so called Langevin observers. This chart is often attributed to Max Born, due to his 1909 work on the relativistic physics of a rotating body. For overview of the application of accelerations in flat spacetime, see Acceleration and proper reference frame.

Shape context is a feature descriptor used in object recognition. Serge Belongie and Jitendra Malik proposed the term in their paper "Matching with Shape Contexts" in 2000.

<span class="mw-page-title-main">Geographical distance</span> Distance measured along the surface of the earth

Geographical distance or geodetic distance is the distance measured along the surface of the earth. The formulae in this article calculate distances between points which are defined by geographical coordinates in terms of latitude and longitude. This distance is an element in solving the second (inverse) geodetic problem.

The Semi-Lagrangian scheme (SLS) is a numerical method that is widely used in numerical weather prediction models for the integration of the equations governing atmospheric motion. A Lagrangian description of a system focuses on following individual air parcels along their trajectories as opposed to the Eulerian description, which considers the rate of change of system variables fixed at a particular point in space. A semi-Lagrangian scheme uses Eulerian framework but the discrete equations come from the Lagrangian perspective.

<span class="mw-page-title-main">Lagrangian mechanics</span> Formulation of classical mechanics

In physics, Lagrangian mechanics is a formulation of classical mechanics founded on the stationary-action principle. It was introduced by the Italian-French mathematician and astronomer Joseph-Louis Lagrange in his 1788 work, Mécanique analytique.

<span class="mw-page-title-main">Chaotic mixing</span>

In chaos theory and fluid dynamics, chaotic mixing is a process by which flow tracers develop into complex fractals under the action of a fluid flow. The flow is characterized by an exponential growth of fluid filaments. Even very simple flows, such as the blinking vortex, or finitely resolved wind fields can generate exceptionally complex patterns from initially simple tracer fields.

In fluid mechanics, meteorology and oceanography, a trajectory traces the motion of a single point, often called a parcel, in the flow.

Isoline retrieval is a remote sensing inverse method that retrieves one or more isolines of a trace atmospheric constituent or variable. When used to validate another contour, it is the most accurate method possible for the task. When used to retrieve a whole field, it is a general, nonlinear inverse method and a robust estimator.

In computational fluid dynamics, the Stochastic Eulerian Lagrangian Method (SELM) is an approach to capture essential features of fluid-structure interactions subject to thermal fluctuations while introducing approximations which facilitate analysis and the development of tractable numerical methods. SELM is a hybrid approach utilizing an Eulerian description for the continuum hydrodynamic fields and a Lagrangian description for elastic structures. Thermal fluctuations are introduced through stochastic driving fields. Approaches also are introduced for the stochastic fields of the SPDEs to obtain numerical methods taking into account the numerical discretization artifacts to maintain statistical principles, such as fluctuation-dissipation balance and other properties in statistical mechanics.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

References

  1. D. W. Waugh; R. A. Plumb (1994). "Contour advection with surgery: a technique for investigating the fine scale structure in tracer transport". Journal of the Atmospheric Sciences . 51 (4): 415–422. doi:10.1175/1520-0469(1994)051<0530:CAWSAT>2.0.CO;2.
  2. 1 2 D. G. Dritschel (1988). "Contour surgery: A topological reconnection scheme". Journal of Computational Physics . 77: 240–266. doi:10.1016/0021-9991(88)90165-9.
  3. 1 2 Peter Mills (2009). "Isoline retrieval: An optimal method for validation of advected contours" (PDF). Computers & Geosciences. 35 (11): 2020–2031. arXiv: 1202.5659 . doi:10.1016/j.cageo.2008.12.015. S2CID   1637949.