In this paper, we initiate the study of the instability of naked singularities without symmetries. In a series of papers, Christodoulou proved that naked singularities are not stable in the context of the spherically symmetric Einstein equations coupled with a massless scalar field. We study in this paper the next simplest case: a characteristic initial value problem of this coupled system with the initial data given on two intersecting null cones, the incoming one of which is assumed to be spherically symmetric and singular at its vertex, and the outgoing one of which has no symmetries. It is shown that, arbitrarily fixing the initial scalar field, the set of the initial conformal metrics on the outgoing null cone such that the maximal future development does not have any sequences of closed trapped surfaces approaching the singularity, is of first category in the whole space in which the shear tensors are continuous. Such a set can then be viewed as exceptional, although the exceptionality is weaker than the at least 1 co-dimensionality in spherical symmetry. Almost equivalently, it is also proved that, arbitrarily fixing an incoming null cone C––ε to the future of the initial incoming null cone, the set of the initial conformal metrics such that the maximal future development has at least one closed trapped surface before C––ε, contains an open and dense subset of the whole space. Since the initial scalar field can be chosen such that the singularity is naked if the initial shear is set to be zero, we may say that the spherical naked singularities of a self-gravitating scalar field are not stable under gravitational perturbations. This in particular gives new families of non-spherically symmetric gravitational perturbations different from the original spherically symmetric scalar perturbations given by Christodoulou.

A new augmented method is proposed for elliptic interface problems with a piecewise variable coefficient that has a finite jump across a smooth interface. The main motivation is to get not only a second order accurate solution but also a second order accurate gradient from each side of the interface. Key to the new method is introducing the jump in the normal derivative of the solution as an augmented variable and rewriting the interface problem as a new PDE that consists of a leading Laplacian operator plus lower order derivative terms near the interface. In this way, the leading second order derivative jump relations are independent of the jump in the coefficient that appears only in the lower order terms after the scaling. An upwind type discretization is used for the finite difference discretization at the irregular grid points on or near the interface so that the resulting coefficient matrix is an M-matrix. A multigrid solver is used to solve the linear system of equations, and the GMRES iterative method is used to solve the augmented variable. Second order convergence for the solution and the gradient from each side of the interface is proved in this paper. Numerical examples for general elliptic interface problems confirm the theoretical analysis and efficiency of the new method.

In this work, we consider the numerical solution of the nonlinear Schrödinger equation with a highly oscillatory potential (NLSE-OP). The NLSE-OP is a model problem which frequently occurs in recent studies of some multiscale dynamical systems, where the potential introduces wide temporal oscillations to the solution and causes numerical difficulties. We aim to analyze rigorously the error bounds of the splitting schemes for solving the NLSE-OP to a fixed time. Our theoretical results show that the Lie–Trotter splitting scheme is uniformly and optimally accurate at the first order provided that the oscillatory potential is integrated exactly, while the Strang splitting scheme is not. Our results apply to general dispersive or wave equations with an oscillatory potential. The error estimates are confirmed by numerical results.

As V. I. Arnold observed in the 1960s, the Euler equations of incompressible fluid flow correspond formally to geodesic equations in a group of volume-preserving diffeomorphisms. Working in an Eulerian framework, we study incompressible flows of shapes as critical paths for action (kinetic energy) along transport paths constrained to be shape densities (characteristic functions). The formal geodesic equations for this problem are Euler equations for incompressible, inviscid potential flow of fluid with zero pressure and surface tension on the free boundary. The problem of minimizing this action exhibits an instability associated with microdroplet formation, with the following outcomes: Any two shapes of equal volume can be approximately connected by an Euler spray---a countable superposition of ellipsoidal geodesics. The infimum of the action is the Wasserstein distance squared, and is almost never attained except in dimension 1. Every Wasserstein geodesic between bounded densities of compact support provides a solution of the (compressible) pressureless Euler system that is a weak limit of (incompressible) Euler sprays. Each such Wasserstein geodesic is also the unique minimizer of a relaxed least-action principle for a two-fluid mixture theory corresponding to incompressible fluid mixed with vacuum.

The general volume of a star body, a notion that includes the usual volume, the $q$th dual volumes, and many previous types of dual mixed volumes, is introduced. A corresponding new general dual Orlicz curvature measure is defined that specializes to the $(p,q)$-dual curvature measures introduced recently by Lutwak, Yang, and Zhang. General variational formulas are established for the general volume of two types of Orlicz linear combinations. One of these is applied to the Minkowski problem for the new general dual Orlicz curvature measure, giving in particular a solution to the Minkowski problem posed by Lutwak, Yang, and Zhang for the $(p,q)$-dual curvature measures when $p>0$ and $q<0$. A dual Orlicz-Brunn-Minkowski inequality for general volumes is obtained, as well as dual Orlicz-Minkowski-type inequalities and uniqueness results for star bodies. Finally, a very general Minkowski-type inequality, involving two Orlicz functions, two convex bodies, and a star body, is proved, that includes as special cases several others in the literature, in particular one due to Lutwak, Yang, and Zhang for the $(p,q)$-mixed volume.