This article is intended to fill in the blank of the numerical schemes with second-order convergence accuracy in time for nonlinear Stokes’ first problem for a heated generalized second grade fluid with fractional derivative. A linearized difference scheme is proposed. The time fractional-order derivative is discretized by second-order shifted and weighted Gr¨unwald-Letnikov difference operator. The convergence accuracy in space is improved by performing the average operator. The presented numerical method is unconditionally stable with the global convergence order of in maximum norm, where τ and h are the step sizes in time and space, respectively. Finally, numerical examples are carried out to verify the theoretical results, showing that our scheme is efficient indeed.

In the paper, we develop and analyze a new mass-preserving splitting domain decomposition method over multiple sub-domains for solving parabolic equations. The domain is divided into non-overlapping multi-bock sub-domains. On the interfaces of sub-domains, the interface fluxes are computed by the semi-implicit (explicit) flux scheme. The solutions and fluxes in the interiors of sub-domains are computed by the splitting one-dimensional implicit solution-flux coupled scheme. The important feature is that the proposed scheme is mass conservative over multiple non-overlapping sub-domains. Analyzing the mass-preserving S-DDM scheme is difficult over non-overlapping multi-block sub-domains due to the combination of the splitting technique and the domain decomposition at each time step. We prove theoretically that our scheme satisfies conservation of mass over multi-block non-overlapping sub-domains and it is unconditionally stable. We further prove the convergence and obtain the error estimate in L2-norm. Numerical experiments confirm theoretical results.

The paper is concerned with the unconditional stability and error estimates of fully discrete Galerkin-Galerkin FEMs for the equations of incompressible miscible flows in porous media. We prove that the optimal L2 error estimates hold without any time-step (convergence) conditions, while all previous works require certain time-step restrictions. Theoretical analysis is based on a splitting of the error into two parts: the error from the time discretization of the PDEs and the error from the finite element discretization of the corresponding time-discrete PDEs, which was proposed in our previous work [26, 27]. Numerical results for both two and three-dimensional flow models are presented to confirm our theoretical analysis.

The symmetric energy-conserved splitting FDTD scheme developed in is a very new and efficient scheme for computing the Maxwell’s equations. It is based on splitting the whole Maxwell’s equations and matching the x-direction and y-direction electric fields associated to the magnetic field symmetrically. In this paper, we make further study on the scheme for the 2D Maxwell’s equations with the PEC boundary condition. Two new energy-conserved identities of the symmetric EC-S-FDTD scheme in the discrete H1-norm are derived. It is then proved that the scheme is uncondi-tionally stable in the discrete H1-norm. By the new energy-conserved identities, the super-convergence of the symmetric EC-S-FDTD scheme is further proved that it is of second order convergence in both time and space steps in the discrete H1-norm. Numerical experiments are carried out and confirm our theoretical results.