In [H. Guo, Q. Zhang, J. Wang, Applied Mathematics and Computation, 259 (2015), 88-105], a nonconservative local discontinuous Galerkin (LDG) method for both flow and transport equations was introduced for the one-dimensional coupled system of compressible miscible displacement problem. In this paper, we will continue our effort and develop a conservative LDG method for the problem in two space dimensions. Optimal error estimates in $L^{\infty}(0, T; L^{2})$ norm for not only the solution itself but also the auxiliary variables will be derived. The main difficulty is how to treat the inter-element discontinuities of two independent solution variables (one from the flow equation and the other from the transport equation) at cell interfaces. Numerical experiments will be given to confirm the accuracy and efficiency of the scheme.

We present two uniformly accurate numerical methods for discretizing the Zakharovsystem (ZS) with a dimensionless parameter 0< ε ≤ 1, which is inversely proportional to theacoustic speed. In the subsonic limit regime, i.e., 0< ε << 1, the solution of ZS propagates waves with O(ε)- andO(1)-wavelengths in time and space, respectively, and/or rapid outgoing initial layerswith speed O(1/ε) in space due to the singular perturbation of the wave operator in ZS and/or theincompatibility of the initial data. By adopting an asymptotic consistent formulation of ZS, wepresent a time-splitting exponential wave integrator (TS-EWI) method via applying a time-splittingtechnique and an exponential wave integrator for temporal derivatives in the nonlinear Schr ̈odingerequation and wave-type equation, respectively. By introducing a multiscale decomposition of ZS, wepropose a time-splitting multiscale time integrator (TS-MTI) method. Both methods are explicitand convergent exponentially in space for all kinds of initial data, which is uniformly for ε ∈ (0,1].The TS-EWI method is simpler to be implemented and it is only uniformly and optimally accuratein time for well-prepared initial data, while the TS-MTI method is uniformly and optimally accuratein time for any kind of initial data. Extensive numerical results are reported to show their efficiencyand accuracy, especially in the subsonic limit regime. Finally, the TS-MTI method is applied tostudy numerically convergence rates of ZS to its limiting models when ε→0+.

In this paper, we review some recent progress made in [4, 5, 6] on finite difference schemes for viscous incompressible flows using vorticity formulation. The main purpose of this series of papers [4, 5, 6] is to resurrect the idea of using local vorticity boundary condition for unsteady calculation. The emphasis is on simplicity of the methods. Three main issues will be discussed: efficient time-stepping procedures and cell Reynolds number constraints, efficient methods in 3D on non-staggered grids and efficient high order methods using compact differencing.

In this paper, we illustrate a new concept regarding unitary elements defined on Lorentz cone, and establish some basic properties under the so-called unitary transformation associated with Lorentz cone. As an application of unitary transformation, we achieve a weaker version of the triangle inequality and several (weak) majorizations defined on Lorentz cone.

In this paper, we propose an unfitted Nitsche's method to compute the band structures of phononic crystal with impurities of general geometry. The proposed method does not require the background mesh to fit the interfaces of impurities, and thus avoids the expensive cost of generating body-fitted meshes and simplifies the inclusion of interface conditions in the formulation. The quasi-periodic boundary conditions are handled by the Floquet-Bloch transform, which converts the computation of band structures into an eigenvalue problem with periodic boundary conditions. More importantly, we show the well-posedness of the proposed method using a delicate argument based on the trace inequality, and further prove the convergence by the Babuška-Osborn theory. We achieve the optimal convergence rate at the presence of the impurities of general geometry. We confirm the theoretical results by two numerical examples, and show the capability of the proposed methods for computing the band structures without fitting the interfaces of impurities.