A system of singularly perturbed convection-diffusion equations with Robin boundary conditions is considered on the interval \([0,1]\). It is shown that any solution of such a problem can be expressed to a system of first-order singularly perturbed initial value problem, which is discretized by the backward Euler formula on an arbitrary nonuniform mesh. An a posteriori error estimation in maximum norm is derived to design an adaptive grid generation algorithm. Besides, in order to establish the initial values of the original problems, we construct a nonlinear optimization problem, which is solved by the Nelder–Mead simplex method. Numerical results are given to demonstrate the performance of the presented method.
In this paper, an adaptive grid method is proposed to solve one-dimensional unsteady singularly perturbed Burger-Huxley equation with appropriate initial and boundary conditions. Firstly, we use the classical backward-Euler scheme on a uniform mesh to approximate time derivative. The resulting nonlinear singularly perturbed semi-discrete problem is linearized by using Newton-Raphson-Kantorovich approximation method which is quadratically convergent. Then, an upwind finite difference scheme on an adaptive nonuniform grid is used for space derivative. The nonuniform grid is generated by equidistribution of a positive monitor function, which is similar to the arc-length function. It is shown that the presented adaptive grid method is first order uniform convergent in the time and spatial directions, respectively. Finally, numerical results are given to validate the theoretical results.
1 January 2024 in Networks and Heterogeneous Media
In this paper, a second finite difference method on a graded grid is proposed for a Volterra integro-differential equation with a weakly singular kernel. The proposed scheme is obtained by using the two-step backward differentiation formula (BDF2) to discretize the first derivative term and the first-order interpolation scheme to approximate the integral term. The analysis of stability is proved and used to prove the convergence of our presented numerical method in the discrete maximum norm. Finally, Numerical experiments are given to verify the theoretical results.
A nonlinear fractional differential equation with a Caputo derivative of order is studied. This problem is discretized by using the L1 scheme on an arbitrary nonuniform mesh. By utilizing the Taylor expansion with integral remainder term, an optimal local truncation error estimation of L1 scheme is proved. Based on this truncation error estimation and the mesh equidistribution principle, a new monitor function is constructed to construct an adaptive grid generation algorithm. Numerical experiments are performed to confirm the accuracy of our new adaptive grid algorithm.
1 January 2025 in Networks and Heterogeneous Media
In this paper, a higher order nonsymmetric interior penalty Galerkin (NIPG) method on a Bakhvalov mesh is developed for a weakly coupled system of singularly perturbed reaction-diffusion equations. At first, by selecting special penalty parameters at different mesh points, the supercloseness of the $ k+\frac{1}{2} $ order of our proposed method is derived, where $ k $ is the degree of polynomials space. Then, an optimal order of uniform convergence analysis in a balanced norm is performed. Finally, some numerical experiments are given to support our theoretical findings.