a b s t r a c tThe 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.

The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Terms of UseArticle is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. In analogy to gas-dynamical detonation waves, which consist of a shock with an attached exothermic reaction zone, we consider herein nonlinear traveling wave solutions to the hyperbolic ͑"inviscid"͒ continuum traffic equations. Generic existence criteria are examined in the context of the Lax entropy conditions. Our analysis naturally precludes traveling wave solutions for which the shocks travel downstream more rapidly than individual vehicles. Consistent with recent experimental observations from a periodic roadway ͓Y. Sugiyama et al., N. J. Phys. 10, 033001 ͑2008͔͒, our numerical calculations show that nonlinear traveling waves are attracting solutions, with the time evolution of the system converging toward a wave-dominated configuration. Theoretical principles are elucidated by considering examples of traffic flow on open and closed roadways. Self-sustained nonlinear waves in traffic flow

In this paper we present a method to treat interface jump conditions for constant coefficients Poisson problems that allows the use of standard "black box" solvers, without compromising accuracy. The basic idea of the new approach is similar to the Ghost Fluid Method (GFM). The GFM relies on corrections applied on nodes located across the interface for discretization stencils that straddle the interface. If the corrections are solution-independent, they can be moved to the right-hand-side (RHS) of the equations, producing a problem with the same linear system as if there were no jumps, only with a different RHS. However, achieving high accuracy is very hard (if not impossible) with the "standard" approaches used to compute the GFM correction terms.In this paper we generalize the GFM correction terms to a correction function, defined on a band around the interface. This function is then shown to be characterized as the solution to a PDE, with appropriate boundary conditions. This PDE can, in principle, be solved to any desired order of accuracy. As an example, we apply this new method to devise a 4 th order accurate scheme for the constant coefficients Poisson equation with discontinuities in 2D. This scheme is based on (i) the standard 9-point stencil discretization of the Poisson equation, (ii) a representation of the correction function in terms of bicubics, and (iii) a solution of the correction function PDE by a least squares minimization. Several applications of the method are presented to illustrate its robustness dealing with a variety of interface geometries, its capability to capture sharp discontinuities, and its high convergence rate.

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.2000 Mathematics Subject Classification. 65M25; 65M12; 35L04. The term "jet scheme" exists in the fields of algebra and algebraic geometry, introduced by Nash [20] and popularized by Kontsevich [16], as a concept to understand singularities. Since we are here dealing with a computational scheme for a partial differential equation, we expect no possibility for confusion.

The reference map, defined as the inverse motion function, is utilized in an Eulerian-frame representation of continuum solid mechanics, leading to a simple, explicit finite-difference method for solids undergoing finite-deformations. We investigate the accuracy and applicability of the technique for a range of finite-strain elasticity laws under various geometries and loadings. Capacity to model dynamic, static, and quasi-static conditions is shown. Specifications of the approach are demonstrated for handling irregularly shaped and/or moving boundaries, as well as shock solutions. The technique is also integrated within a fluid-solid framework using a level-set to discern phases and using a standard explicit fluid solver for the fluid phases. We employ a sharp-interface method to institute the interfacial conditions, and the resulting scheme is shown to quickly capture FSI solutions in several examples.