Application of the boundary element method to parabolic type equations


Tezin Türü: Doktora

Tezin Yürütüldüğü Kurum: Orta Doğu Teknik Üniversitesi, Fen Bilimleri Enstitüsü, Fen Bilimleri Enstitüsü, Türkiye

Tezin Onay Tarihi: 2010

Öğrenci: NURAY BOZKAYA

Danışman: MÜNEVVER TEZER

Özet:

In this thesis, the two-dimensional initial and boundary value problems governed by unsteady partial differential equations are solved by making use of boundary element techniques. The boundary element method (BEM) with time-dependent fundamental solution is presented as an efficient procedure for the solution of diffusion, wave and convection-diffusion equations. It interpenetrates the equations in such a way that the boundary solution is advanced to all time levels, simultaneously. The solution at a required interior point can then be obtained by using the computed boundary solution. Then, the coupled system of nonlinear reaction-diffusion equations and the magnetohydrodynamic (MHD) flow equations in a duct are solved by using the time-domain BEM. The numerical approach is based on the iteration between the equations of the system. The advantage of time-domain BEM are still made use of utilizing large time increments. Mainly, MHD flow equations in a duct having variable wall conductivities are solved successfully for large values of Hartmann number. Variable conductivity on the walls produces coupled boundary conditions which causes difficulties in numerical treatment of the problem by the usual BEM. Thus, a new time-domain BEM approach is derived in order to solve these equations as a whole despite the coupled boundary conditions, which is one of the main contributions of this thesis. Further, the full MHD equations in stream function-vorticity-magnetic induction-current density form are solved. The dual reciprocity boundary element method (DRBEM), producing only boundary integrals, is used due to the nonlinear convection terms in the equations. In addition, the missing boundary conditions for vorticity and current density are derived with the help of coordinate functions in DRBEM. The resulting ordinary differential equations are discretized in time by using unconditionally stable Gear's scheme so that large time increments can be used. The Navier-Stokes equations are solved in a square cavity up to Reynolds number 2000. Then, the solution of full MHD flow in a lid-driven cavity and a backward facing step is obtained for different values of Reynolds, magnetic Reynolds and Hartmann numbers. The solution procedure is quite efficient to capture the well known characteristics of MHD flow.