In this paper, we use the weak Galerkin (WG) finite element method to solve the mixed form linear elasticity problem. In the mixed form, we get the discrete of proximation of the stress tensor and the displacement field. For the WG methods, we define the weak function and the weak differential operator in an optimal polynomial approximation spaces. The optimal error estimates are given and numerical results are presented to demonstrate the efficiency and the accuracy of the weak Galerkin finite element method.

In this article, two block-centered finite difference schemes are introduced and analyzed to solve the parabolic integro-differential equation arising in modeling non-Fickian flow in porous media. One scheme is Euler backward scheme with first order accuracy in time increment while the other is Crank-Nicolson scheme with second order accuracy in time increment. Stability analysis and second-order error estimates in spatial mesh size for both pressure and velocity in discrete L² norms are established on non-uniform rectangular grid. Numerical experiments using the schemes show that the convergence rates are in agreement with the theoretical analysis.

Recently, a kind of high order hybrid methods based on Runge-Kutta discontinuous Galerkin (RKDG) method and weighted essentially non-oscillatory finite difference (WENO-FD) scheme was proposed. Those methods are computationally efficient, however, stable problems might sometimes be encountered in practical applications. In this work, we first analyze the linear stabilities of those methods based on the Heuristic theory. We find that the conservative hybrid method is linearly unstable if the numerical flux at the coupling interface is chosen to be 'downstream'. Then we introduce two ways of healing this defect. One is to choose the numerical flux at the coupling interface to be 'upstream'. The other is to employ a slope limiter function to enforce the hybrid method satisfying the local total variation diminishing (TVD) condition. In the end, numerical experiments are provided to validate the effectiveness of the proposed methods.

A quasi-Newton waveform relaxation (WR) algorithm for semi-linear reaction-diffusion equations is presented at first in this paper. Using the idea of energy estimate, a general proof method for convergence of the continuous case and the discrete case of quasi-Newton WR is given, which appears to be the superlinear rate. The semi-linear wave equation and semi-linear coupled equations can similarly be solved by quasi-Newton WR algorithm and be proved as convergent with the energy inequalities. Finally several parallel numerical experiments are implemented to confirm the effectiveness of the above theories.

This paper presents an investigation on the anomalous diffusion in finite length fingers comb frame, the time and space Riesz fractional Cattaneo-Christov flux is introduced with the Oldroyds' upper convective derivative and the effect of Poiseuille flow is also taken into account. Formulated governing equation possesses the coexisting characteristics of parabolicity and hyperbolicity. Numerical solution is obtained by the L1-scheme and shifted Grünwald formulae, which is verified by introducing a source item to construct an exact solution. The effects, such as time and space fractional parameters, relaxation parameter and the ratio of the pressure gradient and viscosity coefficient, on the spatial and temporal evolution of particles distribution and dynamic characteristics are shown graphically and analyzed in detail.

The preconditioner for parameterized inexact Uzawa methods has been used to solve some indefinite saddle point problems. Firstly, we modify the preconditioner by making it more generalized, then we use theoretical analyses to show that the iteration method converges under certain conditions. Moreover, we discuss the optimal parameter and matrices based on these conditions. Finally, we propose two improved methods. Numerical experiments are provided to show the effectiveness of the modified preconditioner. All methods have fantastic convergence rates by choosing the optimal parameter and matrices.

In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrödinger-KdV (CNLS-KdV) equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing some canonical momenta. To simulate the problem efficiently, the CNLS-KdV equations are approximated by a high order compact method in space which preserves $N$ semi-discrete multisymplectic conservation laws. We then discretize the semi-discrete system by using a symplectic midpoint scheme in time. Thus, a full-discrete multisymplectic scheme is obtained for the CNLS-KdV equations. The conservation laws of the full-discrete scheme are analyzed. Some numerical experiments are presented to further verify the convergence and conservation laws of the new scheme.

We present a new Dirichlet boundary condition for the rate-type non-Newtonian diffusive constitutive models. The newly proposed boundary condition is compared with two such well-known and popularly used boundary conditions as the pure Neumann condition [1] and the Dirichlet condition by Sureshkumar and Beris [2]. Our condition is demonstrated to be more stable and robust in a number of numerical test cases. A new Dirichlet boundary condition is implemented in the framework of the finite difference Marker and Cell (MAC) method. In this paper, we also present an energy-stable finite difference MAC scheme that preserves the positivity for the conformation tensor and show how the addition of the diffusion helps the energy-stability in a finite difference MAC scheme-setting.