Numerical Approximations of a Class of Nonlinear Second-Order Boundary Value Problems using Galerkin-Compact Finite Difference Method

Shovan Sourav Datta Pranta,Md. Shafiqul Islam

2023-06-17

Abstract:In this study, we examine numerical approximations for 2nd-order linear-nonlinear differential equations with diverse boundary conditions, followed by the residual corrections of the first approximations. We first obtain numerical results using the Galerkin weighted residual approach with Bernstein polynomials. The generation of residuals is brought on by the fact that our first approximation is computed using numerical methods. To minimize these residuals, we use the compact finite difference scheme of 4th-order convergence to solve the error differential equations in accordance with the error boundary conditions. We also introduce the formulation of the compact finite difference method of fourth-order convergence for the nonlinear BVPs. The improved approximations are produced by adding the error values derived from the approximations of the error differential equation to the weighted residual values. Numerical results are compared to the exact solutions and to the solutions available in the published literature to validate the proposed scheme, and high accuracy is achieved in all cases

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to improve the accuracy of numerical solutions for second - order linear and nonlinear boundary value problems (BVPs). Specifically, the author proposes a method that combines the Galerkin weighted residual method and the fourth - order compact finite difference method to solve such problems, and further improves the accuracy of the solution by modifying the residual.

### Specific problem description

1. **Research object**:
- The author focuses on second - order linear - nonlinear boundary value problems in the following form:
\[
\frac{d^2f}{dx^2}=g(x, f, \frac{df}{dx}), \quad a \leq x \leq b
\]
- The boundary conditions can be of Dirichlet, Neumann or Robin types.

2. **Limitations of existing methods**:
- Existing numerical methods (such as the Galerkin method) can provide preliminary approximate solutions, but these solutions usually have certain errors.
- An effective method is needed to reduce these errors in order to obtain more accurate solutions.

3. **Proposed new method**:
- **First step**: Use the Galerkin weighted residual method combined with Bernstein polynomials to obtain an initial approximate solution.
- **Second step**: Use the fourth - order compact finite difference method to solve the error differential equation, thereby correcting the errors in the initial approximate solution.
- **Final result**: By adding the error correction value to the initial approximate solution, an improved high - accuracy solution is obtained.

### Advantages of the method

- **High accuracy**: By introducing an error correction mechanism, the accuracy of the numerical solution is significantly improved.
- **High efficiency**: Compared with other methods, a high accuracy can be achieved with fewer basis functions (such as Bernstein polynomials), reducing the computation time.
- **Generality**: It is applicable to multiple types of boundary conditions (Dirichlet, Neumann, Robin), and is effective for both linear and nonlinear problems.

### Experimental verification

The author verifies the effectiveness of the proposed method through several specific examples (such as linear problems, nonlinear problems and the strongly nonlinear Bratu problem), and shows its convergence and stability under different grid sizes.

In summary, this paper aims to propose an efficient and high - accuracy numerical method for solving second - order linear and nonlinear boundary value problems by combining the Galerkin weighted residual method and the fourth - order compact finite difference method, and verifies the effectiveness of this method through experiments.