Boundary element methods for transient convective diffusion. Part II: 2D implementation Original Research Article

A boundary element method for transient convective diffusion phenomena presented in Part I of the paper is extended to two dimensional problems. We introduce a series representation for the transient convective kernel and perform a time integration for the double integrals to evaluate coefficients of the time-discrete boundary integral equation. The time-integrated kernels are evaluated for the linear, quadratic and quartic time interpolation functions utilized in the paper. Then, linear, quadratic and quartic boundary elements as well as bi-linear, bi-quadratic and bi-quartic volume cells are introduced to ensure proper resolution in space for the two-dimensional formulation. Due to the singular nature of the transient convective diffusion kernels, integration of the kernels over the boundary elements and volume cells requires a considerable effort to maintain a desired level of accuracy. We define influence domains due to time-integrated and time-delayed kernels arising for the surface and volume integrals, respectively. Note that the kernel influences are extremely localized due to the convective nature of the kernels, thus, the surface and volume integrations are performed only within these domains of influence. The localization of the kernels becomes more prominent as the Peclet number of the flow increases. Due to increasing sparsity of the global matrix, iterative solvers become the primary choice for the convective diffusion problems.