A Strum-Liouville integral transform technique is novelly applied to solve system of coupled nonlinear boundary-value problems approximately. The systems of differential equations consist of a linear differential operator and a nonlinear function of the dependent variables. To illustrate the potential of this technique we consider an example which comes from the modeling of diffusion and nonlinear chemical reaction systems in chemical engineering. The approximate solutions obtained by our technique agree surprising well with the numerically exact solutions obtained by the orthogonal collocation technique. To improve the approximation an iteration scheme in transform space is also defined. Scope—Today, mathematical modeling of physical phenomena often produces (single or coupled) nonlinear differential equations. The true physical situation can, in many cases, be more closely described if the differential equations are allowed to be nonlinear. However, nonlinear differential equations are generally too difficult to be solved analytically apart from a few “tricks” or substitutions which apply only to a handful of equations [1]. An alternative approach is to look for a method which will reduce the problem, via analytical techniques, to a point where a “simple” computer program can solve the rest of the problem. The method introduced in this paper belongs to this class of solution techniques. The method, which in this paper is applied to solving coupled nonlinear boundary-value problems, is a generalization of an idea in a paper by Do and Bailey [2] who apply it to a single nonlinear differential equation of boundary-value type. The equations, to which the technique is applied, arise from Fick's law diffusion into a porous solid and nonlinear reaction within the solid. The solution method employs a Strum-Loiuville integral transform and to account for the nonlinear part an approximation is introduced. An iteration scheme is defined to improved the accuracy of the solution. The system of coupled nonlinear differential equations is reduced to a system of coupled nonlinear algebraic equations which is solved using a Newton-Raphson process. Finally, the solution is expressed as an infinite series, which is summed using a computer. In response to papers by Do and Bailey [3] and Do and Weiland [4], Jerri [5] has tried to put this method on a more mathematical footing, and he shows that this method is a special case of a more general technique he has devised. Jerri uses the idea of Fourier transforms and convolution products to justify his method. The results for the example he considered are good, but he did not state how many iterations he required to obtain the solutions reported. Conclusions and Significance—This paper has presented a very powerful method of solving boundary-value problems with linear operators and a nonlinear function of the dependent variable. The method works well for a single equation or coupled equations and can handle any kind of nonlinear function. We have shown through extensive numerical calculation the accuracy of this solution method, where the accuracy is measured in terms of a ratio of norms. In most cases an error of 4% can be achieved with just one iteration (Tables 2 and 3). Even though the present method has been applied to problems which have arisen from the modeling of chemical engineering problems, it would also be applicable to differential equations arising in other areas, provided they are of the same form.
Read full abstract