On contour crossings in contour-advective simulations – part 1 – algorithm for detection and quantification

• Published in: Journal of computational physics Vol. 231; no. 2; pp. 465 - 480
• Main Authors: Schaerf, T.M., Macaskill, C.
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Inc 2012; Elsevier
• Subjects: Algorithms; CASL; Computational techniques; Computer simulation; Contour crossings; Errors; Exact sciences and technology; Fluid flow; Mathematical analysis; Mathematical methods in physics; Numerical error; Physics; Quasigeostrophic turbulence; Representations; Shape; Turbulent flow; CASL; Numerical error; Contour crossings; Quasigeostrophic turbulence; Turbulence; Algorithms; Lagrangian; Calculation; Polygon; Calculation methods
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2011.09.015

This is the first of two papers devoted to the analysis of contour crossing errors that occur in contour-advective simulations of fluid motion. Here an algorithm is presented for quantifying the error due to contour crossings. The first step is to determine the relative proximity of all possible pairs of contours. A digital representation of each contour is produced to aid in the corresponding calculation. Simple analysis of functions is then used to find any crossings between contours deemed close to each other by the digital representation method. Next, the area in error of a pair of crossing contours is calculated by identifying the polygon or polygons that approximately bound the erroneous region. Finally, some preliminary results of analysis of contour crossings that occur in contour-advective semi-lagrangian (CASL) simulations of single layer quasigeostrophic turbulence are presented. It is shown that the error due to contour crossings is small in the simulations considered here.