Recent Advances and Emerging Applications of the Boundary Element Method Y

Total Page:16

File Type:pdf, Size:1020Kb

Recent Advances and Emerging Applications of the Boundary Element Method Y Recent Advances and Emerging Applications of the Boundary Element Method Y. J. Liu, S. Mukherjee, N. Nishimura, M. Schanz, W. Ye, A. Sutradhar, E. Pan, A. Dumont, A. Frangi, A. Saez To cite this version: Y. J. Liu, S. Mukherjee, N. Nishimura, M. Schanz, W. Ye, et al.. Recent Advances and Emerging Applications of the Boundary Element Method. Applied Mechanics Reviews, American Society of Mechanical Engineers, 2011, 64 (3), 10.1115/1.4005491. hal-01401752 HAL Id: hal-01401752 https://hal.archives-ouvertes.fr/hal-01401752 Submitted on 23 Nov 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution| 4.0 International License 3FDFOU "EWBODFT BOE &NFSHJOH "QQMJDBUJPOT PG UIF #PVOEBSZ &MFNFOU .FUIPE 4QPOTPSFE CZ UIF 64 /BUJPOBM 4DJFODF 'PVOEBUJPO B XPSLTIPQ PO UIF CPVOEBSZ FMFNFOU NFUIPE #&. XBT IFME PO UIF DBNQVT PG UIF 6OJWFSTJUZ PG "LSPO EVSJOH 4FQUFNCFS o /4' i8PSLTIPQ PO UIF &NFSHJOH "QQMJDBUJPOT BOE 'VUVSF %JSFDUJPOT PG UIF #PVOEBSZ &MFNFOU .FUIPE w 6OJWFSTJUZ PG "LSPO 0IJP 4FQUFNCFS o 5IJT QBQFS XBT QSFQBSFE BGUFS UIJT XPSLTIPQ CZ UIF PSHBOJ[FST BOE QBSUJDJQBOUT CBTFE PO UIF QSFTFOUBUJPOT BOE EJTDVTTJPOT BU UIF XPSLTIPQ 5IF QBQFS BJNT UP SFWJFX UIF NBKPS SFTFBSDI BDIJFWFNFOUT JO UIF MBTU EFDBEF UIF DVSSFOU TUBUVT BOE UIF GVUVSF EJSFDUJPOT PG UIF #&. JO UIF OFYU EFDBEF 5IF SFWJFX TUBSUT XJUI B CSJFG JOUSPEVDUJPO UP UIF #&. 5IFO OFX EFWFMPQNFOUT JO (SFFOT GVODUJPOT TZNNFUSJD (BMFSLJO GPSNVMBUJPOT CPVOEBSZ NFTIGSFF NFUIPET BOE WBSJBUJPOBMMZ CBTFE #&. GPSNVMBUJPOT BSF SFWJFXFE /FYU GBTU TPMVUJPO NFUIPET GPS FGmDJFOUMZ TPMWJOH UIF #&. TZTUFNT PG FRVBUJPOT OBNFMZ UIF GBTU NVMUJQPMF NFUIPE UIF QSFDPSSFDUFE GBTU 'PVSJFS USBOTGPSNBUJPO NFUIPE BOE UIF BEBQUJWF DSPTT BQQSPYJNBUJPO NFUIPE BSF QSFTFOUFE &NFSHJOH BQQMJDBUJPOT PG UIF #&. JO TPMWJOH NJDSPFMFDUSPNFDIBOJDBM TZTUFNT DPNQPTJUFT GVODUJPOBMMZ HSBEFE NBUFSJBMT GSBDUVSF NFDIBOJDT BDPVTUJD FMBTUJD BOE FMFDUSPNBHOFUJD XBWFT UJNFEPNBJO QSPCMFNT BOE DPVQMFE NFUIPET BSF SFWJFXFE 'JOBMMZ GVUVSF EJSFDUJPOT PG UIF #&. BT FOWJTJPOFE CZ UIF BVUIPST GPS UIF OFYU mWF UP UFO ZFBST BSF EJTDVTTFE 5IJT QBQFS JT JOUFOEFE GPS TUV EFOUT SFTFBSDIFST BOE FOHJOFFST XIP BSF OFX JO #&. SFTFBSDI BOE XJTI UP IBWF BO PWFSWJFX PG UIF mFME 5FDIOJDBM EFUBJMT PG UIF #&. BOE SFMBUFE BQQSPBDIFT EJTDVTTFE JO UIF SFWJFX DBO CF GPVOE JO UIF 3FGFSFODF TFDUJPO XJUI NPSF UIBO QBQFST DJUFE JO UIJT SFWJFX Y. J. Liu Mechanical Engineering, University of Cincinnati, 598 Rhodes Hall, 2600 Clifton Avenue, Cincinnati, OH 45221-0072 e-mail: [email protected] S. Mukherjee Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY 14853 e-mail: [email protected] N. Nishimura Department of Applied Analysis and Complex Dynamical Systems, Graduate School of Informatics, Kyoto University, Kyoto 606-8501, Japan e-mail: [email protected] M. Schanz Institute of Applied Mechanics, Graz University of Technology, Technikerstr. 4, A-8010 Graz, Austria e-mail: [email protected] W. Ye Department of Mechanical Engineering, Hong Kong University of Science and Technology, Clearwater Bay, Kowloon, Hong Kong, China e-mail: [email protected] A. Sutradhar Department of Surgery and Department of Mechanical Engineering, The Ohio State University, 915 Olentangy River Road, Suite 2100, Columbus, OH 43212 e-mail: [email protected] E. Pan Department of Civil Engineering, The University of Akron, ASEC Room 534, Akron, OH 44325-3905 e-mail: [email protected] N. A. Dumont Pontifical Catholic University of Rio de Janeiro, 22451-900 Brazil e-mail: [email protected] A. Frangi Department of Structural Engineering, Politecnico of Milano, P.za L. da Vinci 32, 20133 Milano, Italy e-mail: [email protected] A. Saez Department of Continuum Mechanics and Structural Analysis, University of Seville, Camino de los Descubrimientos s/n, Seville, E-41092, Spain e-mail: [email protected] 1 Table 1 A Comparison of the FEM and BEM Finite element method Boundary element method Unknown DOFs everywhere Unknown DOFs only on boundary in domain for linear problems Banded symmetric matrices Smaller full nonsymmetric matrices Very efficient Less efficient in standard form (fast with new forms) Exterior problems somewhat Exterior problems in infinite domains difficult easy Nonhomogeneous problems Needs Green’s functions for maxi fine mum efficiency Moving boundary problems Efficient for moving boundary can be difficult problems Fig. 1 A 3D domain V with boundary S Needs compatible elements Incompatible=nonconforming ele ments fine Inter element discontinuity of Continuous internal derivatives derivatives 1.2 A Review of the BEM Formulation. We use the 3D potential theory problem as an example to show the key steps and main results in the BEM formulation. For a domain V with bound 1 Introduction ary S (Fig. 1), we have the following boundary value problem: The boundary element method (BEM) is a numerical method for r2/ðxÞ¼0; x 2 V (1) solving boundary value or initial value problems formulated in boundary integral equations (BIEs). In some literature, it is also where / is the potential field. The boundary conditions (BCs) are called the boundary integral equation method because of this rela tionship. The most important advantage of the BEM as compared /ðxÞ¼/ðxÞ; x 2 S/ ðDirichlet BCÞ (2) with other domain based numerical methods is the reduction of the @ dimensions of the problems to be solved, from 3D to 2D, or from 2D / qðxÞ ðxÞ¼qðxÞ; x 2 Sq ðNeumann BCÞ (3) to 1D, leading to much easier mesh generation as compared with @n domain based methods. The second advantage is the accuracy the BEM offers, due to the nature of integrals used in the formulations. in which the overbar indicates the prescribed value, S/ [ Sq ¼ S, The third advantage of the BEM is its capability of accurately model and n is the outward normal of the boundary S (Fig. 1). ing infinite domain problems without introducing additional, often The fundamental solution (also known as full space Green’s cumbersome, conditions at infinity. Table 1 is a comparison of the function) Gðx; yÞ for 3D potential problems is given by BEM with the finite element method (FEM), which may change in 1 the future based on new developments in either one or both methods. Gðx; yÞ¼ (4) 4pr 1.1 A Brief History. The direct BIE formulations (in which where r is the distance between source point x and field point y. the unknown functions have clear physical meanings) and their The normal derivative of G is modern numerical solutions with boundary elements originated more than 40 years ago during the 1960s. The 2D potential prob @Gðx; yÞ 1 Fðx; yÞ ¼ r; n ðyÞ (5) lem was first formulated in terms of a direct BIE and solved @nðyÞ 4pr2 k k numerically by Jaswon, et al. [1 3]. This work was later extended to the 2D elastostatic case by Rizzo in the early 1960s as his dis with r;k ¼ @r=@yk ¼ðyk xkÞ=r. The fundamental solution sertation work at the University of Illinois, which was later pub Gðx; yÞ represents the response (potential) at y due to a unit source lished as a journal article in 1967 [4]. Following these early at x (See Sec. 2.1 for the definition of a Green’s function). works, extensive research efforts were carried out in the BIE for If we apply the second Green’s identity with / and Gðx; yÞ in mulations of many problems in applied mechanics and in the nu domain V with boundary S, we obtain the following representation merical solutions during the 1960s and 1970s [5 20]. The name integral: boundary element method appeared in the middle of the 1970s in ð an attempt to make an analogy with the FEM [21 23]. /ðxÞ¼ ½Gðx; yÞqðyÞ Fðx; yÞ/ðyÞ dSðyÞ; x 2 V (6) Some personal accounts of the historical developments of the S BIE and BEM in various places of the U.S., Europe, and China Equation (6) is an expression of the solution / inside the domain were given by Rizzo [24], Cruse [25], Watson [26], Shippy [27], V for Eq. (1). Once the boundary values of both / and q are Mukherjee [28], Telles [29], and Yao and Du [30]inaspecialissue known on the entire boundary S, Eq. (6) can be applied to calcu of the Electronic Journal of Boundary Elements in 2003 [31]. A late / everywhere in V, if needed. To solve the unknown bound comprehensive review of the heritage and early history of the BIE ary values of / and q on S, we let x tend to the boundary S to and BEM was given by Cheng and Cheng in 2005 [32]. This review obtain the following BIE: provides vivid descriptions of the long history and major contribu ð tions to the mathematical foundations of the integral methods by cðxÞ/ðxÞ¼ ½Gðx; yÞqðyÞ Fðx; yÞ/ðyÞ dSðyÞ; x 2 S (7) the pioneers from the 18th to the first half of the 20th centuries and S the early days of the BEM until the late of 1970s. Another compre = hensive review of the BEM research up to the early of 1980s was in which cðxÞ is a coefficient and cðxÞ¼1 2ifS is smooth around given by Tanaka in 1983 [33]. A special review was given by x. Equation (7) is called a conventional or singular BIE for solving Tanaka, Sladek and Sladek in 1994 [34] on the research efforts in potential problems.
Recommended publications
  • Coupling Finite and Boundary Element Methods for Static and Dynamic
    Comput. Methods Appl. Mech. Engrg. 198 (2008) 449–458 Contents lists available at ScienceDirect Comput. Methods Appl. Mech. Engrg. journal homepage: www.elsevier.com/locate/cma Coupling finite and boundary element methods for static and dynamic elastic problems with non-conforming interfaces Thomas Rüberg a, Martin Schanz b,* a Graz University of Technology, Institute for Structural Analysis, Graz, Austria b Graz University of Technology, Institute of Applied Mechanics, Technikerstr. 4, 8010 Graz, Austria article info abstract Article history: A coupling algorithm is presented, which allows for the flexible use of finite and boundary element meth- Received 1 February 2008 ods as local discretization methods. On the subdomain level, Dirichlet-to-Neumann maps are realized by Received in revised form 4 August 2008 means of each discretization method. Such maps are common for the treatment of static problems and Accepted 26 August 2008 are here transferred to dynamic problems. This is realized based on the similarity of the structure of Available online 5 September 2008 the systems of equations obtained after discretization in space and time. The global set of equations is then established by incorporating the interface conditions in a weighted sense by means of Lagrange Keywords: multipliers. Therefore, the interface continuity condition is relaxed and the interface meshes can be FEM–BEM coupling non-conforming. The field of application are problems from elastostatics and elastodynamics. Linear elastodynamics Ó Non-conforming interfaces 2008 Elsevier B.V. All rights reserved. FETI/BETI 1. Introduction main, this approach is often carried out only once for every time step which gives a staggering scheme.
    [Show full text]
  • On Multigrid Methods for Solving Electromagnetic Scattering Problems
    On Multigrid Methods for Solving Electromagnetic Scattering Problems Dissertation zur Erlangung des akademischen Grades eines Doktor der Ingenieurwissenschaften (Dr.-Ing.) der Technischen Fakultat¨ der Christian-Albrechts-Universitat¨ zu Kiel vorgelegt von Simona Gheorghe 2005 1. Gutachter: Prof. Dr.-Ing. L. Klinkenbusch 2. Gutachter: Prof. Dr. U. van Rienen Datum der mundliche¨ Prufung:¨ 20. Jan. 2006 Contents 1 Introductory remarks 3 1.1 General introduction . 3 1.2 Maxwell’s equations . 6 1.3 Boundary conditions . 7 1.3.1 Sommerfeld’s radiation condition . 9 1.4 Scattering problem (Model Problem I) . 10 1.5 Discontinuity in a parallel-plate waveguide (Model Problem II) . 11 1.6 Absorbing-boundary conditions . 12 1.6.1 Global radiation conditions . 13 1.6.2 Local radiation conditions . 18 1.7 Summary . 19 2 Coupling of FEM-BEM 21 2.1 Introduction . 21 2.2 Finite element formulation . 21 2.2.1 Discretization . 26 2.3 Boundary-element formulation . 28 3 4 CONTENTS 2.4 Coupling . 32 3 Iterative solvers for sparse matrices 35 3.1 Introduction . 35 3.2 Classical iterative methods . 36 3.3 Krylov subspace methods . 37 3.3.1 General projection methods . 37 3.3.2 Krylov subspace methods . 39 3.4 Preconditioning . 40 3.4.1 Matrix-based preconditioners . 41 3.4.2 Operator-based preconditioners . 42 3.5 Multigrid . 43 3.5.1 Full Multigrid . 47 4 Numerical results 49 4.1 Coupling between FEM and local/global boundary conditions . 49 4.1.1 Model problem I . 50 4.1.2 Model problem II . 63 4.2 Multigrid . 64 4.2.1 Theoretical considerations regarding the classical multi- grid behavior in the case of an indefinite problem .
    [Show full text]
  • Coupling of the Finite Volume Element Method and the Boundary Element Method – an a Priori Convergence Result
    COUPLING OF THE FINITE VOLUME ELEMENT METHOD AND THE BOUNDARY ELEMENT METHOD – AN A PRIORI CONVERGENCE RESULT CHRISTOPH ERATH∗ Abstract. The coupling of the finite volume element method and the boundary element method is an interesting approach to simulate a coupled system of a diffusion convection reaction process in an interior domain and a diffusion process in the corresponding unbounded exterior domain. This discrete system maintains naturally local conservation and a possible weighted upwind scheme guarantees the stability of the discrete system also for convection dominated problems. We show existence and uniqueness of the continuous system with appropriate transmission conditions on the coupling boundary, provide a convergence and an a priori analysis in an energy (semi-) norm, and an existence and an uniqueness result for the discrete system. All results are also valid for the upwind version. Numerical experiments show that our coupling is an efficient method for the numerical treatment of transmission problems, which can be also convection dominated. Key words. finite volume element method, boundary element method, coupling, existence and uniqueness, convergence, a priori estimate 1. Introduction. The transport of a concentration in a fluid or the heat propa- gation in an interior domain often follows a diffusion convection reaction process. The influence of such a system to the corresponding unbounded exterior domain can be described by a homogeneous diffusion process. Such a coupled system is usually solved by a numerical scheme and therefore, a method which ensures local conservation and stability with respect to the convection term is preferable. In the literature one can find various discretization schemes to solve similar prob- lems.
    [Show full text]
  • Uncertainties in the Solutions to Boundary Element
    UNCERTAINTIES IN THE SOLUTIONS TO BOUNDARY ELEMENT METHOD: AN INTERVAL APPROACH by BARTLOMIEJ FRANCISZEK ZALEWSKI Submitted in partial fulfillment of the requirements For the degree of Doctor of Philosophy Dissertation Adviser: Dr. Robert L. Mullen Department of Civil Engineering CASE WESTERN RESERVE UNIVERSITY August, 2008 CASE WESTERN RESERVE UNIVERSITY SCHOOL OF GRADUATE STUDIES We hereby approve the thesis/dissertation of ______________Bartlomiej Franciszek Zalewski______________ candidate for the Doctor of Philosophy degree *. (signed)_________________Robert L. Mullen________________ (chair of the committee) ______________Arthur A. Huckelbridge_______________ ______________Xiangwu (David) Zeng_______________ _________________Daniela Calvetti__________________ _________________Shad M. Sargand_________________ (date) _______4-30-2008________ *We also certify that written approval has been obtained for any proprietary material contained therein. With love I dedicate my Ph.D. dissertation to my parents. TABLE OF CONTENTS LIST OF TABLES ………………………………………………………………………. 5 LIST OF FIGURES ……………………………………………………………………… 6 ACKNOWLEDGEMENTS ………………………………………………………..…… 10 LIST OF SYMBOLS ……………………………………………………….……………11 ABSTRACT …………………………………………………………………………...... 13 CHAPTER I. INTRODUCTION ……………………………………………………………… 15 1.1 Background ………………………………………………...………… 15 1.2 Overview ……………………………………………..…………….… 20 1.3 Historical Background ………………………………………….……. 21 1.3.1 Historical Background of Boundary Element Method ……. 21 1.3.2 Historical Background of Interval
    [Show full text]
  • The Boundary Element Method in Acoustics
    The Boundary Element Method in Acoustics by Stephen Kirkup 1998/2007 . The Boundary Element Method in Acoustics by S. M. Kirkup First edition published in 1998 by Integrated Sound Software and is published in 2007 in electronic format (corrections and minor amendments on the 1998 print). Information on this and related publications and software (including the second edi- tion of this work, when complete) can be obtained from the web site http://www.boundary-element-method.com The author’s publications can generally be found through the web page http://www.kirkup.info/papers The nine subroutines and the corresponding nine example test programs are available in Fortran 77 on CD ABEMFULL. The original core routines can be downloaded from the web site. The learning package of subroutines ABEM2D is can be downloaded from the web site. See the web site for a full price list and alternative methods of obtaining the software. ISBN 0 953 4031 06 The work is subject to copyright. All rights are reserved, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, reproduction on microfilm or in any other way, and storage in data banks. Duplication of this publication or parts thereof is subject to the permission of the author. c Stephen Kirkup 1998-2007. Preface to 2007 Print Having run out of hard copies of the book a number of years ago, the author is pleased to publish an on-line version of the book in PDF format. A number of minor corrections have been made to the original, following feedback from a number of readers.
    [Show full text]
  • A Boundary-Spheropolygon Element Method for Modelling Sub-Particle Stress and Particle Breakage
    A Boundary-Spheropolygon Element Method for Modelling Sub-particle Stress and Particle Breakage YUPENG JIANG1, HANS J. HERRMANN2, FERNANDO ALONSO-MARROQUIN*1, 1. School of Civil Engineering, The University of Sydney, Australia 2. Computational Physics for Engineering Materials, ETH Zurich, Switzerland Abstract: We present a boundary-spheropolygon element method that combines the boundary element method (BEM) and the spheropolygon-based discrete element method (SDEM). The interaction between particles is simulated via the SDEM, and the sub-particle stress is calculated by BEM. The framework of BSEM is presented. Then, the accuracy and efficiency of the method are analysed by comparison with both analytical solutions and a well-established finite element method. The results demonstrate that BSEM can efficiently provide instant sub-particle stress for particles with an optimized compromise between computational time and accuracy, thus overcoming most of the disadvantages of existing numerical methods. Keywords: Boundary-Spheropolygon element method; Particle stress state; Particle breakage; Numerical simulation 1. Introduction Particle breakage is ubiquitous for granular materials. Breakage is caused by many reasons, such as strong compressive loading, large shear displacement, extreme temperature, dehydration or chemical effects. Based on laboratory studies, it is concluded that breakage causes changes in particle size distribution [1-4] and particle shape [5,6], which in turn change the state variables, such as density and compressibility, that govern the mechanical behaviour of the bulk granular material. These studies have provided a framework for studying particle breakage. However, due to the limitations in experimental observations, laboratory studies are not able to offer more detailed information about the breakage process.
    [Show full text]
  • Optimal Additive Schwarz Preconditioning for Adaptive 2D IGA
    OPTIMAL ADDITIVE SCHWARZ PRECONDITIONING FOR ADAPTIVE 2D IGA BOUNDARY ELEMENT METHODS THOMAS FUHRER,¨ GREGOR GANTNER, DIRK PRAETORIUS, AND STEFAN SCHIMANKO Abstract. We define and analyze (local) multilevel diagonal preconditioners for isogeo- metric boundary elements on locally refined meshes in two dimensions. Hypersingular and weakly-singular integral equations are considered. We prove that the condition number of the preconditioned systems of linear equations is independent of the mesh-size and the re- finement level. Therefore, the computational complexity, when using appropriate iterative solvers, is optimal. Our analysis is carried out for closed and open boundaries and numerical examples confirm our theoretical results. 1. Introduction In the last decade, the isogeometric analysis (IGA) had a strong impact on the field of scientific computing and numerical analysis. We refer, e.g., to the pioneering work [HCB05] and to [CHB09, BdVBSV14] for an introduction to the field. The basic idea is to utilize the same ansatz functions for approximations as are used for the description of the geometry by some computer aided design (CAD) program. Here, we consider the case, where the geometry is represented by rational splines. For certain problems, where the fundamental solution is known, the boundary element method (BEM) is attractive since CAD programs usually only provide a parametrization of the boundary ∂Ω and not of the volume Ω itself. Isogeometric BEM (IGABEM) has first been considered for 2D BEM in [PGK+09] and for 3D BEM in [SSE+13]. We refer to [SBTR12, PTC13, SBLT13, NZW+17] for numerical experiments, to [HR10, TM12, MZBF15, DHP16, DHK+18, DKSW18] for fast IGABEM based on wavelets, fast multipole, -matrices resp.
    [Show full text]
  • Numerical Analysis of Heat Conduction and Potential Flow Problems Over a Sphere and a 3- Dimensional Prolate Spheroidal Body
    ISTANBUL TECHNICAL UNIVERSITY INSTITUTE OF SCIENCE AND TECHNOLOGY NUMERICAL ANALYSIS OF HEAT CONDUCTION AND POTENTIAL FLOW PROBLEMS OVER A SPHERE AND A 3- DIMENSIONAL PROLATE SPHEROIDAL BODY M.Sc. Thesis by Onur ÖLMEZ, B.Sc. Department: Ocean Engineering Programme: Ocean Engineering Supervisor : Assist. Prof. Şafak Nur Ertürk JUNE 2008 ISTANBUL TECHNICAL UNIVERSITY INSTITUTE OF SCIENCE AND TECHNOLOGY NUMERICAL ANALYSIS OF HEAT CONDUCTION AND POTENTIAL FLOW PROBLEMS OVER A SPHERE AND A 3- DIMENSIONAL PROLATE SPHEROIDAL BODY M.Sc. Thesis by Onur ÖLMEZ, B.Sc. (508051107) Date of submission : 21 April 2008 Date of defence examination: 10 June 2008 Supervisor (Chairman): Assist. Prof. Şafak Nur ERTÜRK (I.T.U) Members of the Examining Committee: Prof. Dr. Ömer GÖREN (I.T.U) Prof. Dr. Serdar BEJ İ (I.T.U) JUNE 2008 İSTANBUL TEKN İK ÜN İVERS İTES İ FEN B İLİMLER İ ENSTİTÜSÜ KÜRESEL VE 3 BOYUTLU PROLAT KÜRESEL GÖVDE ÇEVRES İNDE ISI TRANSFER İ VE POTANS İYEL AKI Ş PROBLEMLER İNİN NUMER İK ANAL İZİ YÜKSEK L İSANS TEZ İ Müh. Onur ÖLMEZ (508051107) Tezin Enstitüye Verildi ği Tarih : 21 Nisan 2008 Tezin Savunuldu ğu Tarih : 10 Haziran 2008 Tez Danı şmanı : Yrd. Doç. Dr. Şafak Nur ERTÜRK ( İ.T.Ü.) Di ğer Jüri Üyeleri: Prof. Dr. Ömer GÖREN ( İ.T.Ü.) Prof. Dr. Serdar BEJ İ ( İ.T.Ü.) HAZ İRAN 2008 ACKNOWLEDGEMENTS I would like to thank to my advisor Assist. Prof. Şafak Nur ERTÜRK who supported me to complete this study and illuminated me whenever I encountered with problems. I would also wish to thank to my friends for sophisticated discussions.
    [Show full text]
  • A Coupled Level Set–Boundary Integral Method for Moving Boundary Simulations
    Interfaces and Free Boundaries 7 (2005), 1–26 A coupled level set–boundary integral method for moving boundary simulations M. GARZON† Department of Applied Mathematics, University of Oviedo, Oviedo, Spain D. ADALSTEINSSON‡ Department of Mathematics, University of North Carolina, Chapel Hill, NC, USA L. GRAY§ Computer Science and Mathematics Division, Oak Ridge National Laboratory AND J. A. SETHIAN¶ Department of Mathematics, University of California, Berkeley, and Mathematics Department, Lawrence Berkeley National Laboratory, Berkeley, CA, USA [Received 27 September 2004 and in revised form 12 February 2005] A numerical method for moving boundary problems based upon level set and boundary integral formulations is presented. The interface velocity is obtained from the boundary integral solution using a Galerkin technique for post-processing function gradients on the interface. We introduce a new level set technique for propagating free boundary values in time, and couple this to a narrow band level set method. Together, they allow us to both update the function values and the location of the interface. The methods are discussed in the context of the well-studied two-dimensional nonlinear potential flow model of breaking waves over a sloping beach. The numerical results show wave breaking and rollup, and the algorithm is verified by means of convergence studies and comparisons with previous techniques. 1. Introduction and overview In this paper, we develop an algorithm for solving moving boundary problems. Our approach is based on a combined level set method and boundary integral method, allowing us to handle topological changes while maintaining the highly accurate aspects of boundary integral formulations. There are many scientific and engineering areas where this capability is expected to be important, such as free surface flows [12] and electrostatically driven flows [24, 36].
    [Show full text]
  • A Non-Standard Finite Element Method Using Boundary Integral Operators
    JOHANNES KEPLER UNIVERSITAT¨ LINZ JKU Technisch-Naturwissenschaftliche Fakult¨at A Non-standard Finite Element Method using Boundary Integral Operators DISSERTATION zur Erlangung des akademischen Grades Doktor im Doktoratsstudium der Technischen Wissenschaften Eingereicht von: DI Clemens Hofreither Angefertigt am: Doktoratskolleg Computational Mathematics Beurteilung: O.Univ.-Prof. Dipl.-Ing. Dr. Ulrich Langer (Betreuung) Prof. Dr. Sergej Rjasanow Linz, Oktober, 2012 Acknowledgments This thesis was created during my employment at the Doctoral Program “Computational Mathematics” (W1214) at the Johannes Kepler University in Linz. First of all I would like to thank my advisor, Prof. Ulrich Langer, for interesting me in the research topic treated in this thesis, for giving me the chance to work on it within the framework of the Doctoral Program, and for his guidance and advice during my work on the thesis. Even before, during my graduate studies and supervising my diploma thesis, he has taught me a lot and shaped my academic career in a significant way. Furthermore, I am very grateful to Clemens Pechstein for the countless hours he has invested in the past four years into discussions which have deepened my insight into many topics. Without him, many results of this thesis would not have been possible. I want to thank all the other past and present students at the Doctoral Program for creating a very pleasant environment in which to discuss and work, as well as for the friendships which have developed during my time at the Program, making this time not only a productive, but also a highly enjoyable one. Furthermore, the director of the Doctoral Program, Prof.
    [Show full text]
  • Application of the Boundary Element Method to Parabolic Type Equations
    1 APPLICATION OF THE BOUNDARY ELEMENT METHOD TO PARABOLIC TYPE EQUATIONS A THESIS SUBMITTED TO THE GRADUATE SCHOOL OF NATURAL AND APPLIED SCIENCES OF MIDDLE EAST TECHNICAL UNIVERSITY BY NURAY BOZKAYA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY IN MATHEMATICS JUNE 2010 Approval of the thesis: APPLICATION OF THE BOUNDARY ELEMENT METHOD TO PARABOLIC TYPE EQUATIONS submitted by NURAY BOZKAYA in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Mathematics Department, Middle East Technical Uni- versity by, Prof. Dr. Canan Ozgen¨ Dean, Graduate School of Natural and Applied Sciences Prof. Dr. Zafer Nurlu Head of Department, Mathematics Prof. Dr. Munevver¨ Tezer-Sezgin Supervisor, Department of Mathematics, METU Examining Committee Members: Prof. Dr. Tanıl Ergenc¸ Mathematics Department, Atılım University Prof. Dr. Munevver¨ Tezer-Sezgin Mathematics Department, METU Assoc. Prof. Dr. Ismail˙ Aydın Civil Engineering, METU Assoc. Prof. Dr. Songul¨ Kaya Merdan Mathematics Department, METU Prof. Dr. Nevzat Guneri¨ Genc¸er Electrical and Electronics Engineering, METU Date: I hereby declare that all information in this document has been obtained and presented in accordance with academic rules and ethical conduct. I also declare that, as required by these rules and conduct, I have fully cited and referenced all material and results that are not original to this work. Name, Last Name: NURAY BOZKAYA Signature : iii ABSTRACT APPLICATION OF THE BOUNDARY ELEMENT METHOD TO PARABOLIC TYPE EQUATIONS Bozkaya, Nuray Ph.D., Department of Mathematics Supervisor : Prof. Dr. Munevver¨ Tezer-Sezgin June 2010, 165 pages 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.
    [Show full text]
  • The Boundary Element Method
    2nd Reading April 23, 2013 16:39 WSPC/0219-8762 196-IJCM 1350037 International Journal of Computational Methods Vol. 10, No. 6 (2013) 1350037 (91 pages) c World Scientific Publishing Company DOI: 10.1142/S0219876213500370 THE BOUNDARY ELEMENT METHOD SUBRATA MUKHERJEE∗ Sibley School of Mechanical and Aerospace Engineering Kimball Hall, Cornell University Ithaca, NY 14853, USA [email protected] YIJUN LIU Department of Mechanical Engineering University of Cincinnati Cincinnati, OH 45221, USA Received 15 September 2010 Accepted 25 September 2011 Published 9 April 2013 The boundary element method (BEM), along with the finite element and finite difference methods, is commonly used to carry out numerical simulations in a wide variety of subjects in science and engineering. The BEM, rooted in classical mathematics of integral equations, started becoming a useful computational tool around 50 years ago. Many researchers have worked on computational aspects of this method during this time. This paper presents an overview of the BEM and related methods. It has three sec- tions. The first, relatively short section, presents the governing equations for classical applications of the BEM in potential theory, linear elasticity and acoustics. The second describes specialized applications in bodies with thin features including micro-electro- mechanical systems (MEMS). The final section addresses current research. It has three subsections that present the boundary contour, boundary node and fast multipole meth- ods (BCM, BNM and FMM), respectively. Several numerical examples are included in the second and third sections of this paper. Keywords: Boundary integral equation; boundary element method; boundary contour method; boundary node method; fast multipole method; potential theory; elasticity; acoustics.
    [Show full text]