Abstract

In this paper, we introduce a Crank-Nicolson split least-squares Galerkin finite element procedure for parabolic integro-differential equations, arising in the modeling of nonlocal reactive flows in porous media. By selecting the least-squares functional properly, the procedure can be split into two independent sub-procedures, one of which is for the primitive unknown and the other is for the flux. By carefully choosing projections, we get optimal order H1(Ω) and L2(Ω) norm error estimates for u and sub-optimal (L2(Ω))d norm error estimate for σ with second-order accuracy in time increment. The numerical examples are given to testify the efficiency of the introduced scheme.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call