Abstract:This article describes a numerical method based on the dual reciprocity boundary elements method (DRBEM) for solving some well-known nonlinear parabolic partial differential equations (PDEs). The equations include the classic and generalized Fisher's equations, Allen-Cahn equation, Newell-Whithead equation, Fitz-HughNagumo equation and generalized Fitz-HughNagumo equation with time-dependent coefficients. The concept of the dual reciprocity is used to convert the domain integral to the boundary that leads to an integration free method. We employ the time stepping scheme to approximate the time derivative, and the linear radial basis functions (RBFs) are used as approximate functions in presented method. The nonlinear terms are treated iteratively within each time step. The developed formulation is verified in some numerical test examples. The results of numerical experiments are compared with analytical solution to confirm the accuracy and efficiency of the presented scheme.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to provide numerical solutions for some well - known nonlinear parabolic partial differential equations (PDEs). Specifically, the author proposes a numerical method based on the Dual Reciprocity Boundary Element Method (DRBEM) to solve one - dimensional nonlinear parabolic partial differential equations. These equations include the classical and generalized Fisher equations, Allen - Cahn equation, Newell - Whithead equation, FitzHugh - Nagumo equation, and the generalized FitzHugh - Nagumo equation with time - dependent coefficients.
### Specific Problem Description
The forms of the equations mentioned in the paper are as follows:
\[ u_t+\nu(t)u_x-\mu(t)u_{xx}-\eta(t)F(u) = 0,\quad(x,t)\in[a,b]\times[0,T], \]
where the initial condition and the boundary conditions are respectively:
\[ u(x,0)=f(x),\quad x\in[a,b], \]
\[ u(a,t)=g_1(t),\quad u(b,t)=g_2(t),\quad t\in[0,T]. \]
Here, \( F(u) \) can be divided into linear and nonlinear parts \( F(u)=F_l(u)+F_n(u) \), and \( \nu(t) \), \( \mu(t) \), and \( \eta(t) \) are arbitrary real - valued functions.
### Solution
In order to eliminate the domain integral and convert it into a boundary integral, the author uses the Dual Reciprocity Boundary Element Method (DRBEM). The specific steps include:
1. **Time Discretization**: Use a time - step scheme to approximate the time derivative.
2. **Radial Basis Functions (RBFs)**: Use linear radial basis functions as approximation functions.
3. **Iterative Processing of Nonlinear Terms**: Iteratively process the nonlinear terms within each time step.
Through this method, the author transforms the original partial differential equation into a system of linear algebraic equations and uses a simple predictor - corrector algorithm to handle the nonlinear terms.
### Verification and Application
The author verifies the effectiveness and accuracy of the proposed method through some numerical test examples and compares the numerical experimental results with the analytical solutions to confirm the accuracy and efficiency of the method.
### Summary
The main contribution of this paper is to provide an effective numerical method for solving a class of important nonlinear parabolic partial differential equations. By introducing the Dual Reciprocity Boundary Element Method, the complex domain integral problem is successfully simplified into a boundary integral problem, thus improving the computational efficiency and accuracy.