Solitary waves bifurcated from edges of Bloch bands in two-dimensional periodic media are determined both analytically and numerically in the context of a two-dimensional nonlinear Schrödinger equation with a periodic potential. Using multiscale perturbation methods, the envelope equations of solitary waves near Bloch bands are analytically derived. These envelope equations reveal that solitary waves can bifurcate from edges of Bloch bands under either focusing or defocusing nonlinearity, depending on the signs of the second-order dispersion coefficients at the edge points. Interestingly, at edge points with two linearly independent Bloch modes, the envelope equations lead to a host of solitary wave structures, including reduced-symmetry solitons, dipole-array solitons, vortex-cell solitons, and so on—many of which have not been reported before to our knowledge. It is also shown analytically that the centers of envelope solutions can be positioned at only four possible locations at or between potential peaks. Numerically, families of these solitary waves are directly computed both near and far away from the band edges. Near the band edges, the numerical solutions spread over many lattice sites, and they fully agree with the analytical solutions obtained from the envelope equations. Far away from the band edges, solitary waves are strongly localized, with intensity and phase profiles characteristic of individual families.

This paper is devoted to developing the Il'in-Allen-Southwell (IAS) parameter-uniform difference scheme on uniform meshes for solving strongly coupled systems of singularly perturbed convection-diffusion equations whose solutions may display boundary and/or interior layers, where strong coupling means that the solution components in the system are coupled together mainly through their first derivatives. By decomposing the coefficient matrix of convection term into the Jordan canonical form, we first construct the IAS scheme for 1-D systems and then extend the scheme to 2-D systems by employing an alternating direction technique. The robustness of the developed IAS scheme is illustrated through a series of numerical examples, including the magnetohydrodynamic duct flow problem with a high Hartmann number. Numerical evidence indicates that the IAS scheme is formally second order accurate in the sense that it is second order convergent when the perturbation parameter $\varepsilon$ is not too small and when the $\varepsilon$ is sufficiently small, the scheme is first order convergent in the discrete maximum norm uniformly in $\varepsilon$.

In this paper, we discuss high order finite difference weighted essentially non-oscillatory (WENO) schemes, coupled with total variation diminishing (TVD) Runge-Kutta (RK) temporal integration, for solving the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology. Since the solutions to this system are non-negative, we discuss a positivity-preserving limiter without compromising accuracy. Analysis is performed to justify the maintanance of third order spatial / temporal accuracy when the limiters are applied to a third order finite difference scheme and third order TVD-RK time discretization for solving this model. Numerical results are also provided to demonstrate these methods up to fifth order accuracy.

Global geometrical optics method is a new semi-classical approach for the high frequency linear waves proposed by the author in [13]. In this paper, we rederive it in a more concise way. It is shown that the right candidate of solution ansatz for the high frequency wave equations is the extended WKB function, other than the WKB function used in the classical geometrical optics approximation. A new and main contribution of this paper is an interface analysis for the Helmholtz equation when the incident wave is of extended WKB-type. We derive asymptotic expressions for the reflected and/or transmitted propagating waves in the general case. These expressions are valid even when the incident rays include caustic points.

In this paper, we develop an efficient and energy stable fully-discrete local discontinuous Galerkin (LDG) method for the Cahn–Hilliard–Hele–Shaw (CHHS) system. The semi-discrete energy stability of the LDG method is proved firstly. Due to the strict time step restriction ( t = O( x4)) of explicit time discretization methods for stability, we introduce a semi- implicit time integration scheme which is based on a convex splitting of the discrete Cahn– Hilliard energy. The unconditional energy stability has been proved for this fully-discrete LDG scheme. The fully-discrete equations at the implicit time level are nonlinear. Thus, the nonlinear Full Approximation Scheme (FAS) multigrid method has been applied to solve this system of algebraic equations, which has been shown the nearly optimal complexity numerically. Numerical results are also given to illustrate the accuracy and capability of the LDG method coupled with the multigrid solver.