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.

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.

. Let $u_t+f(u)_x=0$ be a strictly hyperbolic $n\times n$ system of conservation laws, each characteristic field being linearly degenerate or genuinely nonlinear. In this paper we explicitly define a functional $\Phi=\Phi(u,v)$, equivalent to the $\L^1$ distance, which is almost decreasing i.e., $$ \Phi\big( u(t),~v(t)\big)-\Phi\big( u(s),~v(s)\big)\leq \O(\ve)\cdot (t-s)\quad\hbox{for all}~~t&gt;s\geq 0,$$ for every pair of <i></i>-approximate solutions <i>u</i>, <i>v</i> with small total variation, generated by a wave front tracking algorithm. The small parameter <i></i> here controls the errors in the wave speeds, the maximum size of rarefaction fronts and the total strength of all non-physical waves in <i>u</i> and in <i>v</i>. From the above estimate, it follows that front-tracking approximations converge to a unique limit solution, depending Lipschitz continuously on the initial data, in the ${\vec L}^1$ norm. This provides a new proof of the existence of the

No abstract uploaded!

No abstract uploaded!