Transactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
|
|
- Ethel Barber
- 6 years ago
- Views:
Transcription
1 Variable separation in acoustic radiation problems using Chebyshev polynomials Z. Wozniak* and C. Purcell** *ORTECH Corporation, Mississauga, Ontario, Canada, L5K1B3 **Defence Research Establishment Atlantic, Dartmouth, Nova SWz'a, Canada, B27 JZ7 Abstract A new method of computation of Boundary Element matrices for acoustic radiation problems has been developed. The method utilizes Chebyshev polynomials for separation of the frequency from space variables in the free-space Green's function. The BEM matrices resulting from discretization of the Helmholtz equation are static and can be synthesized for each required frequency using an analytical formula. Significant savings in computer time have been achieved. 1 Introduction The radiation and scattering of sound waves is a vital part of many engineering problems such as fluid-structure interactions, acoustic radiation, earthquake and vibrating foundation problems. As the Boundary Element Method (BEM) has been found to be an effective tool for solving this kind of problem, especially for unbounded domains, it is also recognized that the most time consuming part of the computational process is the formulation of the BEM matrices. Due to the oscillatory nature of the Green's
2 420 Boundary Elements function involved, the order of integration over the element surface has to increase with frequency and often the mesh itself has to be refined to achieve satisfactory accuracy. In the past few years, considerable effort has been expended in developing methods for speeding up the BEM computations. The two different approaches have been explored; series approximation and interpolation. An example of the series approximation is given by Koopman and Benner^. For small values of kr they used a truncated McLaurin series expansion for exp(ikr}. Alternatively, for larger values of kr the farfield approximation was used. The interpolation method was introduced by Schenck and Benthien^. The rapid variation of the Green's function integrands was eliminated by factoring out the exp(-ikrmn) quantity, where rmn is a distance between two reference points on the radiating surface. Modified this way integrands are slowly varying with frequency, thus allowing linear interpolation in frequency domain. Review of the range of methods used for speeding up the BEM solution is given by Kirkup and Henwood ^. They proposed an efficient interpolation method in combination with the collocation method for the Burton and Miller* formulation of BIE. The method presented in this work separates the frequency from space variables using an expansion of the Green's function in terms of the Chebyshev polynomials. This alternate formulation of the fundamental solution to the Helmholtz equation has a significant computational advantage over the standard approach since the integration over the space variables need to be performed only once and can be used subsequently for any frequency required. A similar method** has been successfully developed by the author for modeling of water surface gravity waves. Since the Chebyshev polynomials are bounded in the [-1,1] interval, a better numerical performance of the expansion has been achieved. 2 Theory The fundamental solution to the Helmholtz equation^ satisfying the Sommerfeld radiation condition is the free-space Green's function <9<1> ^K^T+^O] (2) G(P,Q) = ^- (3) where <& is the velocity flow potential, r = \P Q\ is the distance between any two points P and Q in the acoustic medium or on the body surface,
3 Boundary Elements 421 and k is the wavenumber. In order to separate the wavenumber from the space variable r in Equation (3), the following Chebyshev expansion of the exponential function is examined where Jn(p) are Bessel's functions of thefirstkind and Tn(q) are the Chebyshev polynomials. The above series converges uniformly for all p > 0 and for q within the range [-1,1]. Since Tn(q) are bounded within [-!,+!], the convergence of the series is governed mainly by the Bessel's functions Jn(p)- The error of the partial sums 5oc - &, of the Equation (4) shown in Figure (1) was computed assuming 7^ 1 for all n. n=l (4) kr/2 = 0.1 kr/2=1.0 kr/2 = 5.0 kr/2=10 kr/2= Number of terms Figure 1: Error of approximation \S^ Sn using Chebyshev expansion It can be seen that for n» p the error is rapidly decreasing due to decrease of the Bessel's terms Jn(p) as the index n increases. The magnitude of the terms truncated can be estimated using the asymptotic estimate^ applies /-(*) ~( (5) which indicates that terms of the expansion for n» k are monotonically decreasing as shown in Figure (1) The argument of the Chebyshev polynomials requires mapping of the space variable r on the [-1,1] interval. Such a mapping is always possible since practically all sound radiating bodies have a limited maximum dimension (6)
4 422 Boundary Elements for both internal and external problems. Thus a new space variable p may be defined as follows P = 2^-l (7) and consequently the Green's function can now be rewritten as Substituting expansion (4) for the second term in the above formula (8) yields the alternate equivalent form of the Green's function 1 00 G(P,Q) = e-< '*-[Jo(kR/2] - 2 i"jn(kr/2)t»(p)] (9) with the wave number k separated from space variable r. It can be seen that the convergence of the above formula is governed by the argument of the Bessel's function kr/2. A better convergence exists for smaller wavenumbers as indicated by the estimate (5). The Chebyshev polynomials can be generated using well known recurrence relation To(p) = 1 Ti(p) = P (10) Tk+i(p) = 3 Numerical Implementation Following the standard procedure in BEM\ the boundary integral equation is based on the Helmholtz integral equation (ii) with boundary conditions specified on the body surface F = FI -f where the coefficient c(p) = 2?r for a node on the smooth boundary or ;(V (13) for an edge node or a corner node. The I\ and 1^ are the parts of the body surface with specified potential and the normal gradient of the potential, respectively.
5 Boundary Elements 423 In the course of discretization, the surface integral (11) is replaced by a sum of integrals over a number of surface elements nodes m, C,;$,; = (Aifr - B^) J = l elem,. -/-. AH = Is^N^dS^ (14) m=l Bij = *E' Is^NjGjdSm m=l where i and j are the node numbers. Here, the mapping of the global-to-local coordinates is accomplished by using the shape functions Nj(r, s) X(r,s) = ^XjNj(r,s) (15) j For convenience, the same shape functions are being used for interpolation of thefluidpotential 0 together with the gradient d$/dn on the boundary where j is the surface element node number. The discretized version of the surface integral (11) is given in a form J 00 (A]{*} = [B}{^} (17) where complex matrices A and B are assembled row- by-row for each frequency from surface element integrals. The advantage of the proposed series formulation of the Green's function is that it makes possible to synthesize these matrices using formulas derived from the Chebyshev expansion (4). Substituting Equation (4) into the formulas for element vectors A^ and B^ (14), and grouping terms, one can obtain alternate expressions for the element vectors oc"=0 (18) B% = E#k&uk 6=0 where the coefficients ak(kr/2) depend on frequency only
6 424 Boundary Elements and the coefficients a^, CLijkt &^d 6^-& depend on the space variable only, and may be obtained by integration over the element surfaces Sm elem. at* = Js^NftgdSn ra=l at* = "izfs^n&sidsn (20) 771=1 elem. bijk = E Js^NftdSm 771=1 Since integrands of the coefficients a^ and a^ (20) are well behaved in r, they can be integrated by simple Gaussian quadrature rule. For the coefficients bijk however, the problem of an integrable singularity arises when the "root" node is in the element under consideration. All the Chebyshev polynomials have a value +1 or -1 at the root node and consequently all the integrals b^k (20) are singular. The integration method used in this work is a adaptation of the method given by Lachat and Watson^. The quadrilateral element is divided into two triangles with one vertex located at the point of singularity. As shown elsewhere^, that singularity can be effectively removed by replacing triangles with quadrilaterals collapsed at the "root" node. Assembling of the BEM matrices A and B of Equation (17) for any desired frequency proceeds in a usual manner (Equation 14) with an addition of one inner loop computing element vectors (18) from the Chebyshev expansion. The process of assembling is very fast even if the coefficients of expansion are stored on disk. The most time consuming part is a computation of the coefficients of expansion, but it occurs only once at the initial stage of computation. The time gains of fast assembling of the matrices A and B clearly outpace this initial expense especially if the restart option is implemented and stored coefficients can be used over and over again for analysis involving different frequency ranges. Another important advantage of this method is that the matrices A and B are the continuous functions of the frequency. This feature effectively renders an interpolation unnecessary when A and B are to be evaluated only at certain frequencies. 4 Examples The proposed method of solution to the Helmholtz equation has been implemented into a Boundary Element analysis program. The purpose of the examples presented is to examine performance of the method and compare the results to known analytical and numerical solutions. Both internal and external problems are considered. The first example of an internal problem is an acoustic wave inside a square tube. The piston at x 0 vibrates with uniform velocity V(0) = V exp(-iujt), and at the other end of the tube, at x=l, there is a prescribed
7 Boundary Elements 425 velocity potential $(!/) = 0. The analytical solution for this simple internal problem is $W = V^-^ kcos(kl) (21) The numerical solution to this problem was obtained using a mesh with 88 quadrilateral elements and 116 nodes. Tube dimensions were 1 m x 1 m x 5 m and the velocity V = 1 m/s. _O CD O Q_ BEM tf exact <f> ct d<f/dn Wave number k Figure 2: Transmission of the acoustic wave inside the tube Summary of the results is shown in Figure (2). A second example, of an external problem is calculation of the acoustic radiation from a pulsating sphere. The sphere is centered at r=0 and has radius r=r. There is a prescribed uniform radial velocity V(R) = Ve-xp(-iujt), and therefore complete spherical symmetry exists. The velocity potential for this case is given by Junger and Feit l + ikj? On the surface of the sphere the potential is r (22) (23) A sphere of radius R = 1 m, with the velocity V = 1 m/s was discretized with 384 quadrilateral elements and 386 nodes.
8 426 Boundary Elements Wave number k Figure 3: Velocity potential on the surface of the pulsating sphere The resulting velocity potential on the surface of the sphere is shown in Figure (3). The noticeable discrepancies between BEM and theoretical values at the wave numbers k = TT, 2?r,... are due to nonuniqueness of the BEM solution^ at the frequencies coinciding with the eigenfrequencies of the corresponding inner problem. For both problems, the results obtained using the proposed method based on the Chebyshev expansion are practically undistinguishable (accuracy e 10~^ was imposed) from those obtained using standard BEM approach. This should be expected because both methods are using the same Green's function and the same order of integration was applied during computation. The only difference is in the way the Green's function is represented. It is interesting to compare in Table (1), the execution time for these problems Table 1: Comparison of the execution time, sec Problem Tube Sphere Chebyshev expansion coefficients assembling Standard method assembling All computations were done on an AST Bravo 486/25 PC with a hard disk of 19 ms average access time. The compiler used was Z or tech C+H- version 3.0.
9 Boundary Elements Conclusions A highly efficient method of computation of acoustic problems has been developed which is suitable for both internal and external problems. When using the traditional BEM the formulation of the matrices takes n^tr CPU time, where n^ is the number of wave numbers and i? is the time taken for the computation at each wave number. The new method takes ts + n^n where {3 is the set-up time and t^ is the time required to assemble the BEM matrices for each individual wave number. For the two examples shown, the new method is about 4 to 6 times faster in assembling than the standard BEM procedure. The method is particularly attractive for problems which require computation of a large number of frequency steps or extraction of the natural frequencies of complex fluid-solid systems. Bibliography 1. Burton, A.J., Miller, G.F. The Application of Integral Equation Methods to the Numerical Solution of Some Exterior Boundary Value Problems, Proc. o/zae Ao?/a/ 5ocze%/ Z/oWon, Vol. A323 (1971), Ciskowski, R.D., Brebbia, C. A., Boundary Element Methods in Acoustics, Comp. Mech. Publ. and Elsevier Appl. Science (1992). 3. Jeans, R.A., Mathews, I. C., Solution of fluid-structure interaction problems using coupled finite element and variational boundary element oc./w. 88, (November 1990), Koopman, G.H., Benner, H. Method for computing the sound power of machines based on the Helmholtz integral, J.Acoust.Soc.Am. 71(1) (Jan. 1992), Kirkup, S.M., Kenwood, D.J. Methods for Speeding Up the Boundary Element Solution of Acoustic Radiation Problems, Journal of Vibraand /Icows^cs, 114 (1992), Lachat, J.C., Watson, J.O., Effective numerical treatment of boundary integral equations, Int. J.N urn. Methods Eng. 10, (1976), Lamb, ^..Hydrodynamics, Dover, New York, Paszkowski, S., Numerical Applications of the Chebyshew Polynomials and Chebyshev Series, Pans twowe Wydawnictwo Naukowe, Warsaw (1975).
10 428 Boundary Elements 9. Press, W.H, at al,numerical Recipes in C, Cambridge University Press (1988). 10. Schenck, H.A., Benthien, G.W. The Application of a Coupled Finite- Element Technique to Large-Scale Structural Acoustic Problems, Advances in Boundary Elements, Proceedings of the Eleventh International Conference on BEM, Cambridge, Massachusetts, Brebbia, C.A., and Connor, J.J.,eds., Vol Wozniak, Z.Efficient Algorithm for Computation of the Green's Function in Frequency Domain, Proc. llth Int'l Conf. of the Canadian Applied Mathematics Society, (1990).
An adaptive fast multipole boundary element method for the Helmholtz equation
An adaptive fast multipole boundary element method for the Helmholtz equation Vincenzo Mallardo 1, Claudio Alessandri 1, Ferri M.H. Aliabadi 2 1 Department of Architecture, University of Ferrara, Italy
More informationNew Developments of Frequency Domain Acoustic Methods in LS-DYNA
11 th International LS-DYNA Users Conference Simulation (2) New Developments of Frequency Domain Acoustic Methods in LS-DYNA Yun Huang 1, Mhamed Souli 2, Rongfeng Liu 3 1 Livermore Software Technology
More informationOn the use of multipole methods for domain integration in the BEM
On the use of multipole methods for domain integration in the BEM A.A. Mammoii, M.S. Ingber & M.J. Brown Department of Mechanical Engineering, University of New Mexico, USA Abstract Often, a single application
More informationThe use of exact values at quadrature points in the boundary element method
The use of exact values at quadrature points in the boundary element method N.D. Stringfellow, R.N.L. Smith Applied Mathematics & OR Group, Cranfield University, Shrivenham, Swindon, U.K., SN6 SLA E-mail:
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 HIGH FREQUENCY ACOUSTIC SIMULATIONS VIA FMM ACCELERATED BEM
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 HIGH FREQUENCY ACOUSTIC SIMULATIONS VIA FMM ACCELERATED BEM PACS: 43.20.Fn Gumerov, Nail A.; Duraiswami, Ramani; Fantalgo LLC, 7496
More informationSimulation of acoustic and vibroacoustic problems in LS-DYNA using boundary element method ABSTRACT:
Simulation of acoustic and vibroacoustic problems in LS-DYNA using boundary element method Yun Hang, Mhamed Souli, Rogelio Perez Livermore Software Technology Corporation USA & University of Lille Laboratoire
More informationBEM Methods for Acoustic and Vibroacoustic Problems in LS-DYNA
12 th International LS-DYNA Users Conference Simulation(2) BEM Methods for Acoustic and Vibroacoustic Problems in LS-DYNA Mhamed Souli, Yun Huang, Rongfeng Liu Livermore Software Technology Corporation
More informationFastBEM Acoustics. Verification Manual , Advanced CAE Research, LLC (ACR) Cincinnati, Ohio, USA All Rights Reserved
FastBEM Acoustics Verification Manual 2007-2017, Advanced CAE Research, LLC (ACR) Cincinnati, Ohio, USA All Rights Reserved www.fastbem.com Copyright 2007-2017, Advanced CAE Research, LLC, All Rights Reserved
More informationSimulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method
10 th International LS-DYNA Users Conference Simulation Technolog (2) Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundar Element Method Yun Huang Livermore Software Technolog Corporation
More informationLearning Acoustics through the Boundary Element Method: An Inexpensive Graphical Interface and Associated Tutorials
Learning Acoustics through the Boundary Element Method: An Inexpensive Graphical Interface and Associated Tutorials ABSTRACT Laura A. Brooks, Rick C. Morgans, Colin H. Hansen School of Mechanical Engineering,
More informationAbstract INTRODUCTION
Iterative local minimum search for eigenvalue determination of the Helmholtz equation by boundary element formulation N. Kamiya,E. Andoh, K. Nogae Department ofinformatics and Natural Science, School of
More informationComputing the Acoustic Field of a Radiating Cavity by the Boundary Element - Rayleigh Integral Method (BERIM)
Computing the Acoustic Field of a Radiating Cavity by the Boundary Element - Rayleigh Integral Method (BERIM) Stephen Kirkup and Ambrose Thompson Abstract This paper describes the Fortran subroutine BERIM3
More informationOn spherical-wave scattering by a spherical scatterer and related near-field inverse problems
IMA Journal of Applied Mathematics (2001) 66, 539 549 On spherical-wave scattering by a spherical scatterer and related near-field inverse problems C. ATHANASIADIS Department of Mathematics, University
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Field reconstruction by acoustic holography technique based on BEM Y.A. He, Z.Y. He Acoustic Engineering Institute, Harbin Engineering University, Harbin, Heilongjiang 150001, China Abstract This paper
More informationAcoustic radiation by means of an acoustic dynamic stiffness matrix in spherical coordinates
Acoustic radiation by means of an acoustic dynamic stiffness matrix in spherical coordinates Kauê Werner and Júlio A. Cordioli. Department of Mechanical Engineering Federal University of Santa Catarina
More informationSolution of Kirchhoff's Time-Domain Integral Equation in Acoustics
Solution of Kirchhoff's Time-Domain Integral Equation in Acoustics Dr. George Benthien and Dr. Stephen Hobbs March 29, 2005 E-mail: george@gbenthien.net In this paper we will loo at the time-domain equivalent
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationStructural Acoustics Applications of the BEM and the FEM
Structural Acoustics Applications of the BEM and the FEM A. F. Seybert, T. W. Wu and W. L. Li Department of Mechanical Engineering, University of Kentucky Lexington, KY 40506-0046 U.S.A. SUMMARY In this
More informationSTRUCTURAL-ACOUSTIC COUPLING G.W. Benthien and H.A. Schenck Naval Ocean Systems Center, San Diego, CA , USA
STRUCTURAL-ACOUSTIC COUPLING G.W. Benthien and H.A. Schenck Naval Ocean Systems Center, San Diego, CA 92152-5000, USA 1 Introduction In recent years there has been much interest in problems involving acoustic
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 9, No. 3, pp. 07 6 c 007 Society for Industrial and Applied Mathematics NUMERICAL QUADRATURES FOR SINGULAR AND HYPERSINGULAR INTEGRALS IN BOUNDARY ELEMENT METHODS MICHAEL CARLEY
More informationFirst- and second-order sensitivity analysis schemes by collocation-type Trefftz method
Computer Assisted Mechanics and Engineering Sciences, 4: 477 489, 1997 Copyright c 1997 by Polska Akademia Nauk First- and second-order sensitivity analysis schemes by collocation-type Trefftz method Eisuke
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationSolutions of the Schrödinger equation for an attractive 1/r 6 potential
PHYSICAL REVIEW A VOLUME 58, NUMBER 3 SEPTEMBER 998 Solutions of the Schrödinger equation for an attractive /r 6 potential Bo Gao Department of Physics and Astronomy, University of Toledo, Toledo, Ohio
More informationA comparative study of the direct boundary element method and the dual reciprocity boundary element method in solving the Helmholtz equation
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 2007 A comparative study of the direct boundary element method and the
More informationVerification of Sound Absorption Characteristics Constituted Porous Structure
Verification of Sound Absorption Characteristics Constituted Porous Structure Toru Yoshimachi 1, Ryo Ishii 1, Kuniharu Ushijima 2, Naoki Masuda 2, Takao Yamaguchi 3, Yun Huang 4, Zhe Cui 4 1 JSOL Corporation,
More informationFrank J. Rizzo Department of Aerospace Engineering and Engineering Mechanics, Iowa State University, Ames, Iowa 50011
cattering of elastic waves from thin shapes in three dimensions using the composite boundary integral equation formulation Yijun Liu Department of Mechanical, Industrial and Nuclear Engineering, University
More informationarxiv: v1 [math.na] 14 Feb 2019
Closed-form evaluation of potential integrals in the Boundary Element Method Michael Carley (m.j.carley@bath.ac.uk) arxiv:192.551v1 [math.na] 14 Feb 219 February 15, 219 Abstract A method is presented
More informationA Novel Single-Source Surface Integral Method to Compute Scattering from Dielectric Objects
SUBMITTED TO IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS ON NOVEMBER 18, 2016 1 A Novel Single-Source Surface Integral Method to Compute Scattering from Dielectric Objects Utkarsh R. Patel, Student
More informationTransactions on Modelling and Simulation vol 2, 1993 WIT Press, ISSN X
A boundary integral equation method for water wave-structure interaction problems S. Liapis Aerospace and Ocean Engineering, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
More information4 The matrix approach
4 The matrix approach 4. Introduction It is not a coincidence that the finite element method blossomed at the time when computers were being developed, because it calls for a vast amount of computation
More informationCalculation of sound radiation in infinite domain using a meshless method
PROCEEDINGS of the 22 nd International Congress on Acoustics Structural Acoustics and Vibration: Paper ICA2016-41 Calculation of sound radiation in infinite domain using a meshless method Shaowei Wu (a),
More informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationEfficient calculation for evaluating vast amounts of quadrupole sources in BEM using fast multipole method
PROCEEDINGS of the 22 nd International Congress on Acoustics Boundary Element and Meshless Methods on Acoustics and Vibrations: Paper ICA2016-309 Efficient calculation for evaluating vast amounts of quadrupole
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Analysis of singular and near-singular integrals by the continuation approach D. Rosen & D.E. Cormack Department of Chemical Engineering and Applied Chemistry, University of Toronto, Toronto M5S 1A4, Canada
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More informationarxiv: v1 [math.na] 19 Nov 2018
A Bivariate Spline Solution to the Exterior Helmholtz Equation and its Applications Shelvean Kapita Ming-Jun Lai arxiv:1811.07833v1 [math.na] 19 Nov 2018 November 20, 2018 Abstract We explain how to use
More informationin the CBIE when it is applied to domains surrounding INTRODUCTION
A unified boundary element method for the analysis of sound and shell-like structure interactions. I. Formulation and verification haohai Chen and Yijun Liu a) Department of Mechanical, Industrial and
More information1 Introduction. 2 Boundary Integral Equations
Analysis of sound transmission through a thin elastic plate by using boundary integral equations T. Terai*, Y. Kawai* "Department ofarchitecture, Faculty ofengineering, Kinki University, 1 Umenobe Takaya,
More informationFast multipole boundary element method for the analysis of plates with many holes
Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational
More informationAnalysis of thin acoustic barriers over an
Analysis of thin acoustic barriers over an impedance plane using a dual BEM formulation L.A. de Lacerda*\ L.C. Wrobel*, WJ. Mansur^ "Wessex Institute of Technology, Ashurst, Southampton, UK ^COPPE/URRJ,
More informationMuffler Transmission Loss Simple Expansion Chamber
Muffler Transmission Loss Simple Expansion Chamber 1 Introduction The main objectives of this Demo Model are Demonstrate the ability of Coustyx to model a muffler using Indirect model and solve the acoustics
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationSchur decomposition in the scaled boundary finite element method in elastostatics
IOP Conference Series: Materials Science and Engineering Schur decomposition in the scaled boundary finite element method in elastostatics o cite this article: M Li et al 010 IOP Conf. Ser.: Mater. Sci.
More informationInverse Source Identification based on Acoustic Particle Velocity Measurements. Abstract. 1. Introduction
The 2002 International Congress and Exposition on Noise Control Engineering Dearborn, MI, USA. August 19-21, 2002 Inverse Source Identification based on Acoustic Particle Velocity Measurements R. Visser
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More informationNumerical techniques. Chapter Difference equations
Chapter 6 Numerical techniques The differential equations (4.61), (4.62) and (4.64), in combination with boundary conditions such as equations (4.65) (4.68), constitute a two point boundary value problem.
More informationA Coupled FEM-BEM Formulation in Structural Acoustics for Imaging a Material Inclusion
A Coupled FEM-BEM Formulation in Structural Acoustics for Imaging a Material Inclusion Abhiit Sarkar, Venkata R. Sonti and Rudra Pratap Facility for Research in Technical Acoustics, Department of Mechanical
More informationfree space (vacuum) permittivity [ F/m]
Electrostatic Fields Electrostatic fields are static (time-invariant) electric fields produced by static (stationary) charge distributions. The mathematical definition of the electrostatic field is derived
More informationCh. 03 Numerical Quadrature. Andrea Mignone Physics Department, University of Torino AA
Ch. 03 Numerical Quadrature Andrea Mignone Physics Department, University of Torino AA 2017-2018 Numerical Quadrature In numerical analysis quadrature refers to the computation of definite integrals. y
More informationTransactions on Modelling and Simulation vol 3, 1993 WIT Press, ISSN X
Boundary element method in the development of vehicle body structures for better interior acoustics S. Kopuz, Y.S. Unliisoy, M. Qali kan Mechanical Engineering Department, Middle East Technical University,
More informationFEM/FMBEM coupling for acoustic structure interaction and acoustic design sensitivity analysis with sound-absorbing materials
Boundary Elements and Other Mesh Reduction Methods XXXVIII 113 FEM/FMBEM coupling for acoustic structure interaction and acoustic design sensitivity analysis with sound-absorbing materials Y. M. Xu, H.
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationSPECIAL FUNCTIONS OF MATHEMATICS FOR ENGINEERS
SPECIAL FUNCTIONS OF MATHEMATICS FOR ENGINEERS Second Edition LARRY C. ANDREWS OXFORD UNIVERSITY PRESS OXFORD TOKYO MELBOURNE SPIE OPTICAL ENGINEERING PRESS A Publication of SPIE The International Society
More informationBoundary element methods for acoustics
Boundary element methods for acoustics 1. Solve the integral equation Problem Sheet Part One y(s) = 1 + 1 0 s 2 t 2 y(t)dt by the method which was used in Section 2.1 of the notes to obtain the solution
More informationAEROELASTIC ANALYSIS OF SPHERICAL SHELLS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationLIMITATIONS AND ERROR ANALYSIS OF SPHERICAL MICROPHONE ARRAYS. Thushara D. Abhayapala 1, Michael C. T. Chan 2
ICSV4 Cairns Australia 9-2 July, 2007 LIMITATIONS AND ERROR ANALYSIS OF SPHERICAL MICROPHONE ARRAYS Thushara D. Abhayapala, Michael C. T. Chan 2 Research School of Information Sciences and Engineering
More informationUNCLASSIFIED UNCLASSIFIED DEFENSE DOCUMENTATION CENTER SCIENTIFIC AND TECHNICAL INFORMATION FOR CAMERON STATION ALEXANDRIA VIRGINIA
UNCLASSIFIED AD DEFENSE DOCUMENTATION CENTER FOR SCIENTIFIC AND TECHNICAL INFORMATION CAMERON STATION ALEXANDRIA VIRGINIA DOWHGRADED AT 3 TOAR INTERVALS: DECLASSIFIED ATTER 12 YEARS DOD DIR 520010 UNCLASSIFIED
More informationJournal of Computational Science and Technology
Science and Technology Explicit Evaluation of Hypersingular Boundary Integral Equation for 3-D Helmholtz Equation Discretized with Constant Triangular Element Toshiro MATSUMOTO, Changjun ZHENG,, Shinya
More informationCONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION
Journal of Computational Acoustics, Vol. 8, No. 1 (2) 139 156 c IMACS CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION MURTHY N. GUDDATI Department of Civil Engineering, North Carolina
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More informationH 2 -matrices with adaptive bases
1 H 2 -matrices with adaptive bases Steffen Börm MPI für Mathematik in den Naturwissenschaften Inselstraße 22 26, 04103 Leipzig http://www.mis.mpg.de/ Problem 2 Goal: Treat certain large dense matrices
More informationMETHODS OF THEORETICAL PHYSICS
METHODS OF THEORETICAL PHYSICS Philip M. Morse PROFESSOR OF PHYSICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Herman Feshbach PROFESSOR OF PHYSICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY PART II: CHAPTERS 9
More informationHolographic measuring techniques in acoustics
Holographic measuring techniques in acoustics Ferenc Márki and Fülöp Augusztinovicz Budapest University of Technology and Economics, Dept. of Telecommunications H-1111 Budapest, Sztoczek u., ferko@hit.bme.hu
More informationA Plane Wave Expansion of Spherical Wave Functions for Modal Analysis of Guided Wave Structures and Scatterers
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 51, NO. 10, OCTOBER 2003 2801 A Plane Wave Expansion of Spherical Wave Functions for Modal Analysis of Guided Wave Structures and Scatterers Robert H.
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationON THE COUPLING OF BEM AND FEM FOR EXTERIOR PROBLEMS FOR THE HELMHOLTZ EQUATION
MATHEMATICS OF COMPUTATION Volume 68, Number 227, Pages 945 953 S 0025-5718(99)01064-9 Article electronically published on February 15, 1999 ON THE COUPLING OF BEM AND FEM FOR EXTERIOR PROBLEMS FOR THE
More informationNumerical comparison of two boundary meshless methods for water wave problems
Boundary Elements and Other Mesh Reduction Methods XXXVI 115 umerical comparison of two boundary meshless methods for water wave problems Zatianina Razafizana 1,2, Wen Chen 2 & Zhuo-Jia Fu 2 1 College
More informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationMassachusetts Institute of Technology Instrumentation Laboratory Cambridge, Massachusetts
Massachusetts Institute of Technology Instrumentation Laboratory Cambridge, Massachusetts Space Guidance Analysis Memo #9 To: From: SGA Distribution James E. Potter Date: November 0, 96 Subject: Error
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationFast computation of the Green function for steady ship wave resistance B. Ponizy", M. Ba", M. Gmlbaud\ F. Noblesse"
Fast computation of the Green function for steady ship wave resistance B. Ponizy", M. Ba", M. Gmlbaud\ F. Noblesse" Guillaume VII 86684 Poitiers, France ^C.E.A.T., Universite de Poitiers, 4 3 Route de
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationAXISYMMETRIC ELECTRIC FIELD CALCULATION WITH ZONAL HARMONIC EXPANSION
Progress In Electromagnetics Research B, Vol. 32, 319 350, 2011 AXISYMMETRIC ELECTRIC FIELD CALCULATION WITH ZONAL HARMONIC EXPANSION F. Glück1, 2, * 1 Karlsruhe Institute of Technology, IEKP, POB 3640,
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationA Fast Regularized Boundary Integral Method for Practical Acoustic Problems
Copyright 2013 Tech Science Press CMES, vol.91, no.6, pp.463-484, 2013 A Fast Regularized Boundary Integral Method for Practical Acoustic Problems Z.Y. Qian, Z.D. Han 1, and S.N. Atluri 1,2 Abstract: To
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationMulti Acoustic Prediction Program (MAPP tm ) Recent Results Perrin S. Meyer and John D. Meyer
Multi Acoustic Prediction Program (MAPP tm ) Recent Results Perrin S. Meyer and John D. Meyer Meyer Sound Laboratories Inc., Berkeley, California, USA Presented at the Institute of Acoustics (UK), Reproduced
More informationSome negative results on the use of Helmholtz integral equations for rough-surface scattering
In: Mathematical Methods in Scattering Theory and Biomedical Technology (ed. G. Dassios, D. I. Fotiadis, K. Kiriaki and C. V. Massalas), Pitman Research Notes in Mathematics 390, Addison Wesley Longman,
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationME equations. Cylindrical symmetry. Bessel functions 1 kind Bessel functions 2 kind Modifies Bessel functions 1 kind Modifies Bessel functions 2 kind
Δϕ=0 ME equations ( 2 ) Δ + k E = 0 Quasi static approximation Dynamic approximation Cylindrical symmetry Metallic nano wires Nano holes in metals Bessel functions 1 kind Bessel functions 2 kind Modifies
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationSound Transmission in an Extended Tube Resonator
2016 Published in 4th International Symposium on Innovative Technologies in Engineering and Science 3-5 November 2016 (ISITES2016 Alanya/Antalya - Turkey) Sound Transmission in an Extended Tube Resonator
More informationSound radiation from the open end of pipes and ducts in the presence of mean flow
Sound radiation from the open end of pipes and ducts in the presence of mean flow Ray Kirby (1), Wenbo Duan (2) (1) Centre for Audio, Acoustics and Vibration, University of Technology Sydney, Sydney, Australia
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationCoupling the BEM/TBEM and the MFS for the numerical simulation of acoustic wave propagation and transient conduction heat transfer
Boundary Elements and Other Mesh Reduction Methods XXXII 261 Coupling the BEM/TBEM and the MFS for the numerical simulation of acoustic wave propagation and transient conduction heat transfer A. Tadeu
More informationVibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
More informationApplication of pseudo-symmetric technique in dynamic analysis of concrete gravity dams
Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams V. Lotfi Department of Civil and Environmental Engineering, Amirkabir University, Iran Abstract A new approach is
More informationThaddeus V. Samulski Department of Radiation Oncology, Duke University Medical Center, Durham, North Carolina 27710
An efficient grid sectoring method for calculations of the near-field pressure generated by a circular piston Robert J. McGough a) Department of Electrical and Computer Engineering, Michigan State University,
More informationAEROELASTIC ANALYSIS OF COMBINED CONICAL - CYLINDRICAL SHELLS
Proceedings of the 7th International Conference on Mechanics and Materials in Design Albufeira/Portugal 11-15 June 2017. Editors J.F. Silva Gomes and S.A. Meguid. Publ. INEGI/FEUP (2017) PAPER REF: 6642
More informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
More informationRayleigh wave correction for the BEM analysis of two-dimensional elastodynamic problems in a half-space
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2004; 60:1 16 [Version: 2002/09/18 v2.02] Rayleigh wave correction for the BEM analysis of two-dimensional elastodynamic
More informationA Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions
mathematics of computation volume 39, number 159 july 1982, pages 201-206 A Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions By John P. Boyd Abstract. The theorem proved here extends
More information