Some boundary element methods for heat conduction problems 

Author:  Hamina, Martti 
Organizations:  University of Oulu, Faculty of Technology, Mathematics Division 
Format:  eBook 
Online Access:  PDF Full Text (PDF, 0.8 MB) 
Persistent link:  http://urn.fi/urn:isbn:951425614X 
Language:  English 
Published: 
2000

Publish Date:  20000412 
Thesis type:  Doctoral Dissertation 
Defence Note:  Academic Dissertation to be presented with the assent of the Faculty of Science, University of Oulu, for public discussion in Raahensali (Auditorium L 10), Linnanmaa, on June 21st, 2000, at 12 noon. 
Reviewer: 
Professor Arvet Pedas Professor Rolf Stenberg 
Description: 
AbstractThis thesis summarizes certain boundary element methods applied to some initial and boundary value problems. Our model problem is the twodimensional homogeneous heat conduction problem with vanishing initial data. We use the heat potential representation of the solution. The given boundary conditions, as well as the choice of the representation formula, yield various boundary integral equations. For the sake of simplicity, we use the direct boundary integral approach, where the unknown boundary density appearing in the boundary integral equation is a quantity of physical meaning. We consider two different sets of boundary conditions, the Dirichlet problem, where the boundary temperature is given and the Neumann problem, where the heat flux across the boundary is given. Even a nonlinear Neumann condition satisfying certain monotonicity and growth conditions is possible. The approach yields a nonlinear boundary integral equation of the second kind. In the stationary case, the model problem reduces to a potential problem with a nonlinear Neumann condition. We use the spaces of smoothest splines as trial functions. The nonlinearity is approximated by using the L^{2}orthogonal projection. The resulting collocation scheme retains the optimal L^{2}convergence. Numerical experiments are in agreement with this result. This approach generalizes to the time dependent case. The trial functions are tensor products of piecewise linear and piecewise constant splines. The proposed projection method uses interpolation with respect to the space variable and the orthogonal projection with respect to the time variable. Compared to the Galerkin method, this approach simplifies the realization of the discrete matrix equations. In addition, the rate of the convergence is of optimal order. On the other hand, the Dirichlet problem, where the boundary temperature is given, leads to a single layer heat operator equation of the first kind. In the first approach, we use tensor products of piecewise linear splines as trial functions with collocation at the nodal points. Stability and suboptimal L^{2}convergence of the method were proved in the case of a circular domain. Numerical experiments indicate the expected quadratic L^{2}convergence. Later, a PetrovGalerkin approach was proposed, where the trial functions were tensor products of piecewise linear and piecewise constant splines. The resulting approximative scheme is stable and convergent. The analysis has been carried out in the cases of the single layer heat operator and the hypersingular heat operator. The rate of the convergence with respect to the L^{2}norm is also here of suboptimal order. 
Series: 
Acta Universitatis Ouluensis. A, Scientiae rerum naturalium 
ISSNE:  1796220X 
ISBN:  951425614X 
ISBN Print:  9514256131 
Issue:  344 
Subjects:  
Copyright information:  This publication is copyrighted. You may download, display and print it for your own personal use. Commercial use is prohibited. 