The dual reciprocity boundary element solution of Helmholtz-type equations in fluid dynamics


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: 2013

Öğrenci: NAGEHAN ALSOY AKGÜN

Danışman: MÜNEVVER TEZER

Özet:

In this thesis, the two-dimensional, unsteady, laminar and incompressible fluid flow problems governed by partial differential equations are solved by using dual reciprocity boundary element method (DRBEM). First, the governing equations are transformed to the inhomogeneous modified Helmholtz equations, and then the fundamental solution of modified Helmholtz equation is used for obtaining boundary element method (BEM) formulation. Thus, all the terms in the equation except the modified Helmholtz operator are considered as inhomogeneity. All the inhomogeneity terms are approximated by using suitable radial basis functions, and corresponding particular solutions are derived by using the annihilator method. Transforming time dependent partial differential equations to the form of inhomogeneous modified Helmholtz equations in DRBEM application enables us to use more information from the original governing equation. These are the main original parts of the thesis. In order to obtain modified Helmholtz equation for the time dependent partial differential equations, the time derivatives are approximated at two time levels by using forward finite difference method. This also eliminates the need of another time integration scheme, and diminishes stability problems. Stream function-vorticity formulations are adopted in physical fluid dynamics problems in DRBEM by using constant elements. First, the procedure is applied to the lid-driven cavity flow and results are obtained for Reynolds number values up to 2000. The natural convection flow is solved for Rayleigh numbers between 10^3 to 10^6 when the energy equation is added to the Navier-Stokes equations. Then, double diffusive mixed convection flow problem defined in three different physical domains is solved by using the same procedure. Results are obtained for various values of Richardson and Reynolds numbers, and buoyancy ratios. Behind these, DRBEM is used for the solution of natural convection flow under a magnetic field by using two different radial basis functions for both vorticity transport and energy equations. The same problem is also solved with differential quadrature method using the form of Poisson type stream function and modified Helmholtz type vorticity and energy equations. DRBEM and DQM results are obtained for the values of Rayleigh and Hartmann numbers up to 10^6 and 300, respectively, and are compared in terms of accuracy and computational cost. Finally, DRBEM is used for the solution of inverse natural convection flow under a magnetic field using the results of direct problem for the missing boundary conditions.