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.

Nonlinear hyperbolic systems with relaxations may encounter different scales of relaxation time, which is a prototype multiscale phenomenon that arises in many applications. In such a problem the relaxation time is of O(1) in part of the domain and very small in the remaining domain in which the solution can be approximated by the zero relaxation limit which can be solved numerically much more efficiently. For the Jin–Xin relaxation system in such a two-scale setting, we establish its wellposedness and singular limit as the (smaller) relaxation time goes to zero. The limit is a multiscale coupling problem which couples the original Jin–Xin system on the domain when the relaxation time is O(1) with its relaxation limit in the other domain through interface conditions which can be derived by matched interface layer analysis.As a result, we also establish the well-posedness and regularity (such as boundedness in sup norm with bounded total variation and L 1-contraction) of the coupling problem, thus providing a rigorous mathematical foundation, in the general nonlinear setting, to the multiscale domain decomposition method for this two-scale problem originally proposed in Jin et al. in Math. Comp. 82, 749–779, 2013.

We show how to construct the Perron–Bremermann function by using proper analytic disks. We apply this result to the polynomial hull of a compact set$K$defined on the boundary of the unit ball.

In this paper, we consider the initial boundary value problem in a cylindrical domain to the three dimensional primitive equations with full eddy viscosity in the momentum equations but with only horizontal eddy diffusivity in the temperature equation. Global well-posedness of $z$-weak solution is established for any such initial datum that itself and its vertical derivative belong to $L^2$. This not only extends the results in \cite{Cao5} from the spatially periodic case to general cylindrical domains but also weakens the regularity assumptions on the initial data which are required to be $H^2$ there.

In this work, we give a perturbation theorem for strong polynomial solutions to the zero surface tension Hele-Shaw equation driven by injection or suction, the so called Polubarinova–Galin equation. This theorem enables us to explore properties of solutions with initial functions close to polynomials. Applications of this theorem are given in the suction and injection cases. In the former case, we show that if the initial domain is close to a disk, most of the fluid will be sucked before the strong solution blows up. In the latter case, we obtain precise large-time rescaling behaviors for large data to Hele-Shaw flows in terms of invariant Richardson complex moments. This rescaling behavior result generalizes a recent result regarding large-time rescaling behavior for small data in terms of moments. As a byproduct of a theorem in this paper, a short proof of existence and uniqueness of strong solutions to the Polubarinova–Galin equation is given.