A gradient-augmented level set method with an optimally local, coherent advection scheme

The level set approach represents surfaces implicitly, and advects them by evolving a level set function, which is numerically defined on an Eulerian grid. Here we present an approach that augments the level set function values by gradient information, and evolves both quantities in a fully coupled fashion. This maintains the coherence between function values and derivatives, while exploiting the extra information carried by the derivatives. The method is of comparable quality to WENO schemes, but with optimally local stencils (performing updates in time by using information from only a single adjacent grid cell). In addition, structures smaller than the grid size can be located and tracked, and the extra derivative information can be employed to obtain simple and accurate approximations to the curvature. We analyze the accuracy and the stability of the new scheme, and perform benchmark tests.Comment: 28 pages, 14 figure

The Semi Implicit Gradient Augmented Level Set Method

Here a semi-implicit formulation of the gradient augmented level set method is presented. By tracking both the level set and it's gradient accurate subgrid information is provided,leading to highly accurate descriptions of a moving interface. The result is a hybrid Lagrangian-Eulerian method that may be easily applied in two or three dimensions. The new approach allows for the investigation of interfaces evolving by mean curvature and by the intrinsic Laplacian of the curvature. In this work the algorithm, convergence and accuracy results are presented. Several numerical experiments in both two and three dimensions demonstrate the stability of the scheme.Comment: 19 Pages, 14 Figure

Level Set Jet Schemes for Stiff Advection Equations: The SemiJet Method

Many interfacial phenomena in physical and biological systems are dominated by high order geometric quantities such as curvature. Here a semi-implicit method is combined with a level set jet scheme to handle stiff nonlinear advection problems. The new method offers an improvement over the semi-implicit gradient augmented level set method previously introduced by requiring only one smoothing step when updating the level set jet function while still preserving the underlying methods higher accuracy. Sample results demonstrate that accuracy is not sacrificed while strict time step restrictions can be avoided

A General, Mass-Preserving Navier-Stokes Projection Method

The conservation of mass is common issue with multiphase fluid simulations. In this work a novel projection method is presented which conserves mass both locally and globally. The fluid pressure is augmented with a time-varying component which accounts for any global mass change. The resulting system of equations is solved using an efficient Schur-complement method. Using the proposed method four numerical examples are performed: the evolution of a static bubble, the rise of a bubble, the breakup of a thin fluid thread, and the extension of a droplet in shear flow. The method is capable of conserving the mass even in situations with morphological changes such as droplet breakup.Comment: Submitted to Computer Physics Communication

A comparative study of the efficiency of jet schemes

We present two versions of third order accurate jet schemes, which achieve high order accuracy by tracking derivative information of the solution along characteristic curves. For a benchmark linear advection problem, the efficiency of jet schemes is compared with WENO and Discontinuous Galerkin methods of the same order. Moreover, the performance of various schemes in tracking solution contours is investigated. It is demonstrated that jet schemes possess the simplicity and speed of WENO schemes, while showing several of the advantages as well as the accuracy of DG methods.Comment: 12 pages, 6 figures, presented at the conference Mathematical Modeling and Applications to Industrial Problems 201

A Characteristic Mapping Method for Transport on the Sphere

A semi-Lagrangian method for the solution of the transport equation on a sphere is presented. The method evolves the inverse flow-map using the Characteristic Mapping (CM) [1] and Gradient-Augmented Level Set (GALS) [2] frameworks. Transport of the advected quantity is then computed by composition with the numerically approximated inverse flow-map. This framework allows for the advection of complicated sets and multiple quantities with arbitrarily fine-features using a coarse computational grid. We discuss the CM method for linear transport on the sphere and its computational implementation. Standard test cases for solid body rotation, deformational and divergent flows, and numerical mixing are presented. The unique features of the method are demonstrated by the transport of a multi-scale function and a fractal set in a complex flow environment

Jet schemes for advection problems

We present a systematic methodology to develop high order accurate numerical approaches for linear advection problems. These methods are based on evolving parts of the jet of the solution in time, and are thus called jet schemes. Through the tracking of characteristics and the use of suitable Hermite interpolations, high order is achieved in an optimally local fashion, i.e. the update for the data at any grid point uses information from a single grid cell only. We show that jet schemes can be interpreted as advect-and-project processes in function spaces, where the projection step minimizes a stability functional. Furthermore, this function space framework makes it possible to systematically inherit update rules for the higher derivatives from the ODE solver for the characteristics. Jet schemes of orders up to five are applied in numerical benchmark tests, and systematically compared with classical WENO finite difference schemes. It is observed that jet schemes tend to possess a higher accuracy than WENO schemes of the same order.Comment: 26 pages, 6 figures; based on a poster presentation at the conference Fluid Dynamics, Analysis, and Numerics 201

Computation and optimal perturbation of finite-time coherent sets for aperiodic flows without trajectory integration

Understanding the macroscopic behavior of dynamical systems is an important tool to unravel transport mechanisms in complex flows. A decomposition of the state space into coherent sets is a popular way to reveal this essential macroscopic evolution. To compute coherent sets from an aperiodic time-dependent dynamical system we consider the relevant transfer operators and their infinitesimal generators on an augmented space-time manifold. This space-time generator approach avoids trajectory integration, and creates a convenient linearization of the aperiodic evolution. This linearization can be further exploited to create a simple and effective spectral optimization methodology for diminishing or enhancing coherence. We obtain explicit solutions for these optimization problems using Lagrange multipliers and illustrate this technique by increasing and decreasing mixing of spatial regions through small velocity field perturbations