On solving plane problems of non-stationary heat conduction by the collocation boundary element method

D. Yu. Ivanov

Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mehanika, no. 50 (2017), pp. 9-29

Résumé

In this paper, we propose a fully justified collocation boundary element method allowing one to obtain numerical solutions of internal and external initial-boundary value problems (IBVPs) with boundary conditions of the first, second, and third kind for the equation $\partial_tu=a^2\Delta_2u-pu$ with constants $a, p>0$ in a plane spatial domain $\Omega$ (in a bounded one $\Omega^+$ or in its exterior $\Omega^-$) on a finite time interval $I_t\equiv[0, T]$ at a zero initial condition. The solutions are found in the form of the double-layer potential for the Dirichlet IBVP and in the form of the simple layer potential for the Neumann–Robin IBVP with unknown density functions determined from the boundary integral equations (BIEs) of the second kind. In this paper, instead of the usual piecewise-polynomial interpolation of the density function on time variable $\tau$, the BIEs are approximated by the piecewise-quadratic interpolation (PQI) of the $C_0$-semigroup of right shifts on time. Also, on the basis of the PQI, the approximation of the multiplier $e^{-p\tau}$ in kernels of the integral operators is carried out. In addition, the PQI of density functions is performed: for the BIE, only on arc-length $s$; for the potentials, on both variables $s$ and $\tau$. Then, the integration with respect to the variable $\tau$ on the boundary elements (BEs) is performed exactly. The integration with respect to the variable $s$ on the BE for the potentials is performed approximately by using the Gaussian quadrature with $\gamma\geqslant 2$ points. For the BIE, the integration with respect to the arc-length s is carried out in two ways. On singular BEs and on nearby singular BEs, adjacent to a singular BE in some fixed arc-length region, an exact integration with respect to the variable r is carried out ($r$ is the distance from the boundary point at which the integral is calculated as a function of parameter to the current boundary point of the integration). In this integration, functions of the variable $r$ are taken as the weighting functions. The functions of $r$ are generated by the fundamental solution of the heat equation and the rest of the integrand is approximated by quadratic interpolation on $r$. The integrals with respect to $s$ on the remaining BEs are calculated using the Gaussian quadrature with $\gamma$ points. The cubic convergence of approximate solutions of the IBVP at any point of the set $\Omega\times I_t$ is proved under conditions $\partial\Omega\in C^5\cap C^{2\gamma}$ and $w\in C_{1,3}^3(\partial \Omega)$. It is also proved that such solutions are resistant to perturbations of the boundary function $w$ in the norm of the space $C_1^0(\partial\Omega)$. Here, $C_{m,n}^k(\partial\Omega)\equiv C_m^k(\partial\Omega)\cap C_{m+n}^0(\partial\Omega)$ and $C_m^k(\partial\Omega)$ is the Banach space of $k$ times continuously differentiable on $\partial\Omega$ vector functions with values in Sobolev's space which is the domain of definition of the operator $\mathbf{B^m}$ ($(\mathbf{B f})(t)=f'(t)$, $f(t=0)=0$). In conclusion, results of the numerical experiments are presented. They confirm the cubic convergence of approximate solutions for all three IBVPs in a circular domain.

Détail

Keywords: boundary integral equation, boundary element method, singular boundary elements, non-stationary heat conduction, collocation, operator, approximation, stability.

Bibliographie

[1] Brebbia C. A., Telles J. C. F., Wrobel L. C., Boundary element techniques, Springer-Verlag, New York, 1984, 464 pp. | MR | Zbl

[2] Onishi K., “Convergence in the boundary element method for heat equation”, Teaching for Robust Understanding of Mathematics, 17 (1981), 213–225 | MR | Zbl

[3] Costabel M., Onishi K., Wendland W. L., “A boundary element collocation method for the Neumann problem of the heat equation”, Inverse and Ill-Posed Problems, eds. H.W. Engl, C.W. Groetsch, Academic Press, Boston, 1987, 369–384 | DOI | MR

[4] Hongtao Y., “On the convergence of boundary element methods for initial-Neumann problems for the heat equation”, Mathematics of Computation, 68:226 (1999), 547–557 | DOI | MR | Zbl

[5] Hamina M., Saranen J., “On the spline collocation method for the single layer heat operator equation”, Mathematics of Computation, 62:205 (1994), 41–64 | DOI | MR | Zbl

[6] Ivanov D. Y., “Solution of two-dimensional boundary-value problems corresponding to initial-boundary value problems of diffusion on a right cylinder”, Differential Equations, 46:8 (2010), 1104–1113 | DOI | MR | Zbl

[7] Ivanov D. Yu., Dzerzhinskiy R. I., “Solution of the Robin problems for two-dimensional differential-operator equations describing the thermal conductivity in a straight cylinder”, Nauchno-tekhnicheskiy vestnik Povolzh'ya – Scientific and Technical Journal of the Volga Region, 2016, no. 1, 15–17

[8] Ivanov D. Yu., “Stable solvability in spaces of differentiable functions of some two-dimensional integral equations of heat conduction with an operator-semigroup kernel”, Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics, 2015, no. 6(38), 33–45 | DOI

[9] Berezin I. S., Zhidkov N. P., Methods of computations, v. 2, GIFML, M., 1962 | MR

[10] Fihktengol'ts G. M., A course of differential and integral calculus, v. 2, Fizmatlit, M., 2001

[11] Ivanov D. Yu., “Calculation of operators solving problems of heat conduction in straight cylinders using semigroup symmetry”, Izvestiya Moskovskogo gosudarstvennogo tekhnicheskogo universiteta MAMI – Proceedings of Moscow State Technical University MAMI, 4:4(22) (2014), 26–38