Application of the boundary element method to parabolic type equations

Download
2010
Bozkaya, Nuray
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.

Suggestions

Oscillation of even order nonlinear delay dynamic equations on time scales
Erbe, Lynn; Mert, Raziye; Peterson, Allan; Zafer, Ağacık (Institute of Mathematics, Czech Academy of Sciences, 2013-03-01)
One of the important methods for studying the oscillation of higher order differential equations is to make a comparison with second order differential equations. The method involves using Taylor's Formula. In this paper we show how such a method can be used for a class of even order delay dynamic equations on time scales via comparison with second order dynamic inequalities. In particular, it is shown that nonexistence of an eventually positive solution of a certain second order delay dynamic inequality is...
Least-squares differential quadrature time integration scheme in the dual reciprocity boundary element method solution of diffusive-convective problems
Bozkaya, Canan (Elsevier BV, 2007-01-01)
Least-squares differential quadrature method (DQM) is used for solving the ordinary differential equations in time, obtained from the application of dual reciprocity boundary element method (DRBEM) for the spatial partial derivatives in diffusive-convective type problems with variable coefficients. The DRBEM enables us to use the fundamental solution of Laplace equation, which is easy to implement computation ally. The terms except the Laplacian are considered as the nonhomogeneity in the equation, which ar...
Hyperbolic conservation laws on manifolds with limited regularity
Lefloch, Philippe G.; Okutmuştur, Baver (Elsevier BV, 2008-05-01)
We introduce a formulation of the initial and boundary value problem for nonlinear hyperbolic conservation laws posed on a differential manifold endowed with a volume form, possibly with a boundary; in particular, this includes the important case of Lorentzian manifolds. Only limited regularity is assumed on the geometry of the manifold. For this problem, we establish the existence and uniqueness of an L-1 semi-group of weak solutions satisfying suitable entropy and boundary conditions.
Forced oscillation of second-order nonlinear differential equations with positive and negative coefficients
ÖZBEKLER, ABDULLAH; Wong, J. S. W.; Zafer, Ağacık (Elsevier BV, 2011-07-01)
In this paper we give new oscillation criteria for forced super- and sub-linear differential equations by means of nonprincipal solutions.
Accurate numerical bounds for the spectral points of singular Sturm-Liouville problems over 0 < x < infinity
Taşeli, Hasan (Elsevier BV, 2004-03-01)
The eigenvalues of singular Sturm-Liouville problems defined over the semi-infinite positive real axis are examined on a truncated interval 0<x<l as functions of the boundary point l. As a basic theoretical result, it is shown that the eigenvalues of the truncated interval problems satisfying Dirichlet and Neumann boundary conditions provide, respectively, upper and lower bounds to the eigenvalues of the original problem. Moreover, the unperturbed system in a perturbation problem, where l remains sufficient...
Citation Formats
N. Bozkaya, “Application of the boundary element method to parabolic type equations,” Ph.D. - Doctoral Program, Middle East Technical University, 2010.