The sine-Gordon (SG) equation and perturbed nonlinear Schrdinger (NLS) equations are studied numerically for modeling the propagation of two space dimensional (2D) localized pulses (the so-called <i>light bullets</i>) in nonlinear dispersive optical media. We begin with the (2+1) SG equation obtained as an asymptotic reduction in the two level dissipationless MaxwellBloch system, followed by the review on the perturbed NLS equation in 2D for SG pulse envelopes, which is globally well posed and has all the relevant higher order terms to regularize the collapse of standard critical (cubic focusing) NLS. The perturbed NLS is approximated by truncating the nonlinearity into finite higher order terms undergoing focusingdefocusing cycles. Efficient semi-implicit sine pseudospectral discretizations for SG and perturbed NLS are proposed with rigorous error estimates. Numerical comparison results between light bullet

We introduce two exponential-type integrators for the “good” Boussinesq equation. They are of orders one and two, respectively, and they require lower spatial regularity of the solution compared to classical exponential integrators. For the first integrator, we prove first-order convergence in H^r for solutions in H^{r+1} with r>1/2. This new integrator even converges (with lower order) in H^r for solutions in H^r, i.e., without any additional smoothness assumptions. For the second integrator, we prove second-order convergence in H^r for solutions in H^{r+3} with r>1/2 and convergence in L^2 for solutions in H^3. Numerical results are reported to illustrate the established error estimates. The experiments clearly demonstrate that the new exponential-type integrators are favorable over classical exponential integrators for initial data with low regularity.

We consider an elliptic pseudo-differential equation with a highly oscillating linear potential modeled as a stationary ergodic random field. The random field is a function composed with a centered long-range correlated Gaussian process. In the limiting of vanishing correlation length, the heterogeneous solution converges to a deterministic solution obtained by averaging the random potential. We characterize the deterministic and stochastic correctors. With proper rescaling, the mean-zero stochastic corrector converges to a Gaussian random process in probability and weakly in the spatial variables. In addition, for two prototype equations involving the Laplacian and the fractional Laplacian operators, we prove that the limit holds in distribution in some Hilbert spaces. We also determine the size of the deterministic corrector when it is larger than the stochastic corrector. Depending on the correlation structure of the random field and on the singularities of the Green’s function, we show that either the deterministic or the random part of the corrector dominates.

No abstract uploaded!

In this paper, we consider the Cauchy problem to the heat conductive compressible Navier-Stokes equations in the presence of vacuum and with vacuum far field. Global well-posedness of strong solutions is established under the assumption, among some other regularity and compatibility conditions, that the scaling invariant quantity $\|\rho_0\|_\infty(\|\rho_0\|_3+\|\rho_0\|_\infty^2\|\sqrt{\rho_0}u_0\|_2^2)(\|\nabla u_0\|_2^2+ \|\rho_0\|_\infty\|\sqrt{\rho_0}E_0\|_2^2)$ is sufficiently small, with the smallness depending only on the parameters $R, \gamma, \mu, \lambda,$ and $\kappa$ in the system. Notably, the smallness assumption is imposed on the above scaling invariant quantity exclusively, and it is independent of any norms of the initial data, which is different from the existing papers, see, e.g., \cite{HLX12,HUANGLI11,WENZHU17}. The total mass can be either finite or infinite. A new equation for the density, more precisely for the cubic of the density, derived from combing the continuity and momentum equations, is employed to get the $L^\infty_t(L^3)$ type estimate of the density.