Boundary integral method for multispecies 1-D convection-diffusion-reaction problems

Boundary integral method has been proposed to solve a system of differential equations with a convection-diffusion operator and a nonlinear rate term for each equation. The operator part can be eliminated by using an inverse formulation by a proper choice of the weighting function for each differential equation. This reduces the problem to an integral formulation which is then solved by approximating the dependent variables by cubic osculation. The method has been demonstrated to a packed bed reactor with heat and mass dispersion and for a countercurrent heat exchanger.