“School of Mathematics”
Back to Papers HomeBack to Papers of School of Mathematics
Paper IPM / M / 7374  


Abstract:  
Numerical techniques based on finite difference schemes leading to
parallel algorithms have been developed for obtaining approximate
solutions to an initialboundary value problem for the
twodimensional parabolic partial differential equation (PDE) with
a nonlinear boundary condition. This class of parabolic PDEs
plays a very important role in many branches of science and
engineering. The nonlinear condition is in the form of a double
integral giving the specification of mass in the solution domain.
Not only the problem has both Neumann and Dirichlet boundary
conditions but the Dirichlet boundary condition is in a
nonstandard form. While sharing some common features with the
onedimensional models, the solution of twodimensional equations
are substantially more difficult, thus some considerations are
taken to be able to extend some ideas of onedimensional case. Due
to the structure of the boundary conditions a reduced
onedimensional parabolic equation for the new unknown
v(y,t)=∫^{1}_{0} u(x,y,t)dx can be formulated. The resulting
problem has a nonlocal boundary condition. The new
twodimensional parabolic PDE with Neumann's boundary conditions
will be solved numerically by using the method of lines
semidiscretization approach. The space derivatives in the PDE are
approximated by finite difference replacements. The solution of
the resulting system of firstorder linear ordinary differential
equations satisfies a recurrence relation which involves a matrix
exponential function. The accuracy in time is controlled by
choosing several subdiagonal Pade approximants to replace this
matrix exponential term. Numerical techniques are developed to
compute the required solution using a splitting method, leading to
algorithms for sequential and parallel implementation. The
algorithms are tested on a model problem from the literature. The
article concludes with the results of some numerical experiments.
Download TeX format 

back to top 