This study is concerned with the numerical approximation of the extended Fisher-Kolmogorov equation with a modified boundary integral method. A key aspect of this formulation is that it relaxes the domain-driven approach of a typical boundary element (BEM) technique. While its discretization keeps faith with the second order accurate BEM formulation, its implementation is element-based. This leads to a local solution of all integral equation and their final assembly into a slender and banded coefficient matrix which is far easier to manipulate numerically. This outcome is much better than working with BEM’s fully populated coefficient matrices resulting from a numerical encounter with the problem domain especially for nonlinear, transient, and heterogeneous problems. Faithful results of high accuracy are achieved when the results obtained herein are compared with those available in literature.
The major thrust of the paper is a numerical study of a nonlinear fourth order partial differential equation in a hybrid boundary integral setting.
Many problems in engineering and applied sciences can be represented by fourth order differential equations. A good understanding of the processes that lead to their derivation can only be obtained via numerical solutions. Several examples of these can be found in scientific literature for example, the Euler-Bernoulli beam theory equation [
Earlier work by Danumjaya and Pani [
However a review of current literature clearly demonstrates paucity of information on the application of the boundary integral method to the numerical solution of the EFK equation. The reason is not far-fetched and primarily arises from several factors among which are its transient nature, its fourth order, its dimensionality and nonlinearity. These are factors which pose serious numerical challenge to a typical BEM formulation. The hybrid boundary integral technique applied herein is able to deal satisfactorily with these problems [
This approach is initiated by transforming the governing differential equation into its integral analog via the Green’s second identity followed by a FEM-like element-by-element implementation of the resulting equations throughout the problem domain. Unlike other BEM techniques, no extraneous mathematical techniques are applied to do all integration on the boundary of the problem domain. The whole process is amply simplified by the presence of a source point inside a generic element which facilitates a “local support” based numerical integration akin to the FEM implementation. The gain resulting from this approach is immense because the resulting integral equations appear in the form of local element equations whose contributions add up to produce slender and banded coefficient matrices [
A lot of effort has been expended in dealing with domain integrals arising from BEM implementation (Portapila and Power [
The EFK equation is given by
where
and boundary conditions
We apply a splitting technique by converting Equation (1) into a two coupled system of differential equations.
and Equation (1) becomes
with the following change in one of the boundary conditions
Numerical approximation of Equations (5) and (6) in line with a hybrid boundary integral implementation involves the following steps:
1) We seek the solution of a prescribed auxiliary differential equation of a governing differential equation. This solution (fundamental solution or the free space Green function) constitutes a key element for obtaining the integral analog of the governing differential equation via the Green’s second identity.
2) Classical BEM approach requires that any part of the governing differential equation that does not contribute to the auxiliary differential equation is treated as a boundary integral in the “integralization” procedure. This is the origin of most of the numerical difficulties encountered in BEM numerical implementation. The implication here is that various components of a differential equation whose physics dictates direct interaction with the problem domain are forced to the boundary. We refer to body force terms, heterogeneity, nonlinearity, temporal terms etc. The hybrid method employed herein does not consider domain integration as a disadvantage because it is domain-element-driven. As a result the ensuing computations are relatively straightforward.
3) The finite element component of the hybrid formulation demands that dicretization be carried out throughout the problem domain, with initial and boundary conditions satisfied as well as the enforcement of inter element continuity. Finally all element equations are assembled to form a global coefficient matrix. These are slender and banded because the procedure takes advantage of the “local support” inherent in the finite element formulation.
4) The resulting coefficient matrix together with the right hand side vector of known values (boundary and initial conditions and specified problem parameters) is solved to yield the problem dependent variables. We hasten to comment that boundary element formulation of the hybrid algorithm guarantees a second order accuracy in the computed results.
Numerical discretization of the coupled system is sought by specifying a suitable complementary differential equation whose solution straightforward to obtain. A good example of this follows from previous work [
where G is the Green’s function and
Substitution of the complementary equation, its fundamental solution, and its derivative into the Green’s second identity provides a platform for transforming the governing differential equation into its integro-differ- rential analogs. Without any loss of generality these can be represented by:
where
where
The nonlinearity in Equation (12) can be resolved by the Newton-Raphson or Picard schemes [
In this section we apply the method developed herein to numerically solve the extended Fisher-Kolmogorov equation and compare the displayed results with those available in literature.
We consider the EFK Equation (1) with the following boundary and initial conditions:
and initial condition
The numerical solution of the governing differential equation has been found for
The problem is specified to read:
The initial condition as well as the problem parameters is left the same.
This problem is taken from Danumjaya and Pani [
and an initial condition given by
Changing the value of gamma from
The exact solutions of the EFK equation for the given initial and boundary conditions are not known. In order for us to have meaningful tests, the exact solution is replaced with the EFK numerical solution for 160 partitions and the numerical solutions at different intervals are compared with these. To this end, the following parameters are applied.
No. of Partitions | Order of Convergence | Order of Convergence | ||
---|---|---|---|---|
20 | 0.01123876 | - | 0.0077612 | |
40 | 0.01003659 | 1.86659 | 0.0033576 | 1.86766 |
80 | 0.00326537 | 1.92254 | 0.0011254 | 1.91856 |
160 | 0.00011265 | 2.00312 | 0.0005632 | 2.00136 |
where
We have demonstrated a straightforward hybrid BEM approach for solving a fourth order nonlinear equation of considerable significance in mathematical physics. The approach demonstrated here is simple and straightforward and involves the splitting of the governing differential equations into two second order components followed by an application of BEM-FEM numerical procedure. The utility of the present formulation rests on the fact that all integrations are performed locally and a slender coefficient matrix is created by assembling all local solutions. This approach facilitates the handling of some of the notorious disadvantages that plague the direct BEM formulation and pave the way for more challenging computations.
Okey OselokaOnyejekwe, (2015) A Direct Implementation of a Modified Boundary Integral Formulation for the Extended Fisher-Kolmogorov Equation. Journal of Applied Mathematics and Physics,03,1262-1269. doi: 10.4236/jamp.2015.310155