Hybrid boundary element formulation in acoustics
|
|
- Jasmin Greene
- 6 years ago
- Views:
Transcription
1 Hybrid boundary element formulation in acoustics L. Gaul, M. Wagner Institute A of Mechanics, University of Stuttgart, Germany. Abstract A symmetric hybrid boundary element method in the frequency domain is introduced for the computation of acoustic radiation and scattering in closed and infinite domains. The Hybrid Stress Boundary Element Method (HSBEM) in a frequency domain formulation is based on the dynamical Hellinger-Reissner potential and leads to a Hermitian, frequency-dependent stiffness equation. As compared to previous results published by the authors, new considerations concerning the interpretation of singular contributions in the stiffness matrix are communicated. The field variables are separated into boundary variables, which are approximated by piecewise polynomial functions, and domain variables, which are approximated by a superposition of singular fundamental solutions. This approximation cancels the domain integral over the equation of motion in the hybrid principle and leads to a boundary integral formulation, incorporating singular integrals. This approach results in a linear system of equations with a symmetric stiffness and mass matrix. 1 Introduction The rationale for the search for alternative boundary element methods in dynamics is guided by the demand for simple symmetric formulations, the simplification of field point evaluations and a consistent structure of the spatially discretized equations of motion in boundary element and finite element formulations to ease the treatment of fluid-structure interaction problems For these reasons the authors present a new approach for calculating boundary as well as domain data for dynamic field problems in frequency and time domain.
2 In contrast to domain discretizations, the suggested procedure is a boundary discretization approach which takes advantage of the reduction of the problem dimension by one. It is not derived from a boundary integral equation (BIE) such as the conventional direct boundary element method (CREM) but based on multifield variational principles which are the underlying principles for hybrid finite elements together with the Hellinger- Reissner and the Hu-Washizu principles. The choice of weighted fundamental solutions as domain approximations selected independently from polynomial approximations of boundary variables, allows us to map domain integrals to boundary integrals. The authors generalized known static formulations to dynamic hybrid boundary element formulations by the use of Hamilton's principle in frequency and time domain for field problems of linear acoustics and elastodynamics [4, 5, 8, 91. Symmetric system matrices are obtained by construction whereas the collocation or Galerkin approach to CBEM lead to nonsymmetric matrices. Double integration of the shape functions in the symmetric Bubnov- Galerkin formulation is avoided and replaced by single integration and weighted superposition of the domain approximation. As the bybrid boundary element method (HBEM) in frequency domain is based on the field approximation with weighted dynamic fundamental solutions, the Sommerfeld radiation condition for the exterior domain is fulfilled. This is why radiating irregular shaped surfaces can be discretized by fluid boundary elements. 2 Hybrid boundary element method Other than BEM based on a weighted residual statement, the hybrid boundary element method is derived from variational principles. In this section the HSBEM [l] is applied to acoustical problems in the frequency domain. Considering a compressible, inviscid fluid with linear changes of variables such as velocity and density, wave propagation in the frequency domain in a fluid domain fl with boundary I' can be described by the Helmholtz equation (qi + /G2p = 0 (1) with the velocity potential 4 as field variable, which is defined by +,i = vi where vi is the particle velocity. Domain sources are neglected in the presented formulation. The wave number is defined as K = wlc where c is the speed of sound. Along the boundary, the outward flux 11, is descri'bed by the normal derivative of the velocity potential 2 = =,ini 11, in which ni is the outward normal vector. Mixed boundary conditions are prescribed on the boundaries I'+ and I'+, respectively, as r$(x) = $(t) on I'+ and $(X) = 4(t) on I'+. (2) The HSBEM is derived from the Hellinger-Reissner potential, which is a complementary energy functional. Thy the derived field variab1.e is
3 considered as an independent field. In elastodynamics this is the stress tensor in the field and in acoustics it is the gradient of the velocity potential in the domain and the flux 1C, on the boundary. This gradient generates a potential field $v, which is an exclusive function of the gradient and is a-priori independent of the real potential field Variational Principle In the following, only steady state time-harmonic motion is considered, thus the field variables can be separated into space and time harmonic functions, where only the real part of the complex representation yields the physical result. Complex variables are denoted by 0. The conjugate complex variable is denoted by ()*. The Hellinger-Reissner potential for acoustics in frequency domain reads with the four independent fields, the potential gradient 4,;, the two independent potential fields $ and $V, and a Lagrangian multiplier introducing the static part of the Helmholtz equation. Additionally, essential boundary conditions have to be considered Note that the potential field 4 is assumed only on the boundary, whereas $,i exists on the boundary and in the domain. Applying the fundamental lemma to the vanishing first variation of the potential gives yields an interpretation of the Lagrangian multiplier as Moreover, the governing field equation Eq. (1) arises. Additionally, natural boundary conditions, as well as compatibility conditions of the two potential fields arise. Considering that the variations on the adjoint boundaries vanish and introducing the interpretation of the Lagrangian multiplier in Eq. (6)) yields
4 as a starting point for the derivation of the HSBEM. If the Dirichlet boundary condition is enforced as subsidiary condition, Eq. (7) completely ~describes the acoustical field problem. 2.2 Approximation Functions The numerical solution requires approximation functions for the field variables. The approximation of the boundary variable $ is carried out with a shape function vector N and a nodal vector 4 The approximations of the potential $V and gradient $,i in the domain R, as well as the flux 4 on the boundary, are given as a superposition of n fundamental W), obtained as the solution of where W) is the Dirac-Distribution evaluated at the loadpoint $k) that is collocated with the boundary node and T(~) = lx - <ck)l is the Euclidean distance to the loadpoint ~(~1. The fundamental solutions are weighted by unknown parameters yi, i = l..n. Hence, the approximation of the potential and flux fields are obtained as By modifying the domain R such that small spheres with radii E, centered at the load points 6 where the fundamental solutions are singular, are subtracted, the domain R' with boundary I" is introduced. The properties of the Dirac distribution acting at the points now outside of the considered domain let the domain integral in Eq. (7) vanish in the limit R' -+ 0, if the domain approximations in Eq. (10) are introduced 2.3 Matrix Formulation Inserting the approximations given byeqs. (8) and (10) yields a discretized formulation of the vanishing first variation
5 where OH denotes the complex conjugate transpose. In the absence of domain sources, only boundary integrals remain. The vector f of equivalent nodal fluxes contains the transformations of the given Neumann data to the nodes of the discretization. The matrix H is known as the double-layer potential and will be interpreted as a kinematic transformation matrix in section The frequency-dependent matrix can be expanded in a static part H. and a frequency-dependent part H, For the following, the feature of a frequency-dependent fundamental solution, that it tends to the behavior of the adjoint static fundamental solution in the limit, when the angular frequency or the distance between field and load point vanish is recalled [6]. Due to this the main diagonal of the static part H. contains a strong singularity. The Hermitian matrix F contains the product of the weakly and strongly singular fundamental solutions at the load points i and j. This matrix can also be expanded in the same way A main diagonal entry Fo,, of the static part is made up of a hypersingular integral over the boundary extension FE, similar to the integral-free term in the conventional BEM, and the remaining boundary P. The second integral can be evalutated as finite-part integral, whereas the entries over FE cannot be evaluated directly, since the transformation to spherical coordinates does not regularize a hypersingularity. The computation of these terms will be discussed in detail in section 2.5. Again applying the fundamental lemma yields a matrix equation for the unknowns 7 and 4. Solving for c$ one obtains ~ c $ = f, with the Hermitian stiffness matrix 2.4 Computation of Field Point Data By solving the boundary-value problem in Eq. (14), the primary unknows of the field problem, the potential $ and the generalized flux 7, become known. The potential in the domain can be determined by Eq. (10) only up to a constant. Thus, a complete description of the potential field can be written =(@i+@cci)ri+rq5c, in R. (17)
6 where r is a scaling factor of the additional constant potential This field is arbitrary and can be defined in any appropriate way. Eq. (17) can be looked at as establishing two orthogonal subspaces of the space of potential~ W". The subspace W' contains the potentials that purely produce equilibrated fluxes in the domain, since only these are physically meaningful. In this sense, potentials belonging to W' are called admissible herein. The subspace W contains the constant potentials, that yield no flux. The subspaces are complementary, hence W" = W' + W. It cannot be assured beforehand that the generalized fluxes 7 do not contain a certain amount of unequilibrated fluxes. Thus, a vector ci is introduced that maps the constant potential to the generalized fluxes, compensating for contributions of constant potential. Thus, the first term of Eq. (17) describes only potentials generating equilibrated fluxes and the second term describes constant potentials. As these vectors are in complementary subspaces, the following orthogonality condition for arbitrary yi can be established Eq. (18) yields a way to determine the unknown vector ci. For more detail on the computation of this quantity refer to [2]. One requirement on the field description in Eq. (17) is that it must coincide with the nodal values 4 for compatibility. In a discretized form an orthonormal basis for the N X 1 vector subspace of constant nodal potentials shall be denoted by W. For the nodal values, Eq. (17) reads where 9 is a symmetric N X N matrix of fundamental solutions aij - [(j)l,w). Generally, the computation of field point data involves no additional integration, hence field point data can be gained with much less effort than in the direct boundary element method [l]. This feature can be combined with the conventional boundary element method to obtain an effective way to compute field point data, if symmetry of the system matrices is not required [7]. 2.5 Orthogonality Properties of the HBEM For the computation of the matrix F the hypersingular main diagonal of the static part F. must be evaluated, where the load points of the Dirichlet and Neumann fundamental solutions coincide. For this reason orthogonality properties of this matrix [l] are used to obtain the main diagonal entries without any further integrations Mechanical Coordinate Systems A physical interpretation starts from the virtual work expression for the
7 static matrix F. lim [614*)~dr = 67H~07 = 67 H-v 4. --to 7 is interpreted as generalized flux and ;pv as generalized potential. This expression yields ;pv = F07, (=I which defines F. as a flexibility matrix [l]. For the further development, two sets of coordinates have to be considered. The external coordinate system consists of the potential field 4 on the boundary, approximated by Eq. (8). The potentials cause a nodal flux f. The internal coordinate system uses the generalized flux 7, introduced in the approximation of the domain, Eq. (10). The adjoint potential field ;pv in the internal system is defined via Eq. (21). Note that this internal system does not physically exist but is an energetically equivalent formulation arising from the Hellinger-Reissner potential Consideration of Interior Domains Since constant potentials generate zero flux, it follows for the static stiffness equation that Kow = 0. The matrix KO is singular and W is an orthonormal basis of the nullspace of KO. Considering the virtual work principle, it can be deduced that constant potentials do not produce virtual work, thus the fundamental orthogonality condition holds. In order to solve Eq. (14), relations between the two coordinate systems are necessary. It follows from the first line of Eq. (14), also considering Eq. (21), that -v 4 =Ho~. (23) Thus, the matrix H. transforms potentials from the external into the internal system. Analogously, the second line of Eq. (14) states a transformation relation between the internal and external fluxes. Using this and Eq. (22) shows that the basis of constant potentials also forms the nullspace of the matrix H. How = 0, (24) meaning that constant potentials are not transformed into the internal system, substantiating the fact that the internal system incorporates only equilibrated states. There also exists the left nullspace v of Ho, vtho = 0. This implies that unequilibrated fluxes in the internal system are not transformed into the external system. The conclusion is that if the internal flux belongs to the left nullspace of Ho, 7 E V, it transforms the external flux into the null element of the range of Ho, f E 0. In this case it follows
8 from Eq. (24) that the external potential 4 belongs to the nullspace of No, 4 E W and therefore the internal potential field transforms into the null element, 4v E 0. Finally, it can be deduced from Eq. (21) that the left and right nullspace of F. are identical to the left nullspace of H. This equation is uniquely solvable in a potential problem for the unkcswn main diagonal entries of Fo Orthogonality Properties of i@o In order to apply Eq. (19), the matrix + must be computed. The offdiagonal terms are analytic, but the main diagonal cannot be evaluated due to the singular behaviour of the fundamental solution. As before this singularity is exclusively contained in the static part +o. To compute the main diagonal an orthogonality condition for this matrix is exploited. By construction, the first term in Eq. (19) is set up in such a way that the generalized flux 7 yields only equilibrated potential fields, orthogonal to constant potentials. Then, it can be deduced that unequilibrated generalized fluxes must be orthogonal to the matrix (i@o + wct). Since these are in the subspace V, it follows that as an orthogonality condition. This equation can be used to uniquely compute the yet undetermined main diagonal entries of the static part of the fundamental solution Consideration of Exterior Domains Finally, taking into account Eq. (26), the second case of an unbounded exterior domain fl can be considered. Only recently orthogonality properties for this case have been published [2, 31. Differently from the interior problem, unequilibrated fluxes yield a completely determined solution for the field problem, due to Sommerfelds radiation condition. Since constant potentials are a valid solution, the field in the domain can be computed by evaluating $V E fl contrary to Eq. (17), where an overbar denotes exterior variables. This must also hold for the nodal description, thus 4 = i@ 9. Note that i@ is the same matrix as in the interior case, so that Eq. (26) can be used to compute the singular main-diagonal elements. To obtain and also the boundary solution for the potentials, Eq. (14) for the exterior domain must be solved. The exterior matrices H and P are related to the interior matrices via H+H=I and P+F=9. (27) The first term is a well-known property of the double-layer potential in the direct BEM and the second term can be found by integrating the involved
9 Boundaq Element Techtlology XIV 12 1 matrices. Still, the main diagonal of the static part is undetermined. For this, the orthogonality properties in Eqs. (25) and (26) are exploited. By adding wct and multiplying v to the second term it follows that as a means to determine the hypersingular main diagonal of PO. Note that neither B. nor F. are singular in the exterior case, as it must be expected, since unbalanced fluxes are a valid solution of the field problem. 3 Numerical Example As a numerical example the 2D model of a car compartment is shown here. The boundary is discretized with 105 linear elements and 1063 field-points are used. The domain is excited by a velocity boundary condition as depicted in Fig. l. This models the vibrations of the engine coupled to the passenger compartment. The remaining boundary carries a homogeneous Neumann condition. Fig. 1 shows the solution of the program hybem de- (4 h$em (b) SYSNOISE Figure 1: Distribution of pressure amplitude in car compartment. Excitation frequency is f = 200Hz. veloped by the authors in comparison with the commercial acoustics code SYSNOISE. The distribution of the pressure amplitude in the car compartment at an excitation frequency of f = 200Hz is displayed. The color scales on both pictures are the same. It can be seen that the results agree very well. Note, that the solution close to the boundary of SYSNOISE are worse. This is due to the fact that only every second evaluated field node on the boundary coincides with a boundary node. Obviously the evaluation close to the singular points with the Somigliana identity as done in SYSNOISE generates more error than the evaluation with the HBEM. 4 Conclusion In the present paper a method for the calculation of interior and exterior acoustical problems with the HSBEM has been presented. The HSBEM
10 is developed from the dynamical Hellinger-Reissner potential using fundamental solutions for the approximation of the domain variable, giving rise to matrix F with a hyper-singular main diagonal. The main diagonal cannot be evaluated mathematically, thus a criterion for the evaluation of the main diagonal for an exterior domain has been applied, starting from physical interpretations. The presented example of wave propagation shows that the proposed method yields accurate results. Acknowledgements Support by the DFG of project B5 in SFB 404 "multi-field problems in continuum mechanics" at Stuttgart University is gratefully acknowledged. References [l] N. A. Dumont. The hybrid boundary element method. Applied Mechanics Review, 42(11):54-63, [2] N. A. Dumont and R. A. P. Chaves. The simplified hybrid boundary element method. In CILAMCE XX, Sgo Paulo, Brasilien, [3] N. A. Dumont, L. Gaul, M. Noronha, and M. Wagner. Hybrid boundary element method in acoustics. In P. Gon~alves et al., editors, 6th PACAM, pages , Rio de Janeiro, Brasilien, [4] L. Gaul and M. Wagner. A new hybrid symmetric boundary element method in elastodynamics. In C. Constanda et. al., editors, Integral Methods in Science and Engineering, Pitman Research Notes, pages 3-11, Harlow, UK, [5] L. Gaul and M. Wagner. Beam response derived from 3d hybrid boundary integral method in elastodynamics. Mechanical Systems and Signal Processing, 11(2): , [6] L. Gaul and M. Wagner. Zur Dynamik der Fluid-Struktur-Interaktion. Zeitschrift f. angew. Mathematik und Mechanik, 79(S2): , (71 L. Gaul, M. Wagner, and W. Wenzel. Efficient field point evaluation by combined direct and hybrid boundary element methods. Engineering Analysis with Boundary Element Methods, 21: , [8] L. Gaul, M. Wagner, and W. Wenzel. Hybrid boundary element methods in frequency and time domain. In 0, von Estorff, editor, Boundary Elements in Acoustics, Southampton, WIT Press. [g] M. Wagner. Die hybride Randelementmethode in der Akustik und zur Struktur-Fluid-Interaktion. Bericht aus dem Institut A fiir Mechanik , Universitat Stuttgart, 2000.
Structural 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 informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
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 informationConvergence of modes in exterior acoustics problems with infinite elements
INTER-NOISE 216 Convergence of modes in exterior acoustics problems with infinite elements Lennart Moheit 1 ; Steffen Marburg Chair of Vibroacoustics of Vehicles and Machines Technische Universität München,
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
More informationThe sound power output of a monopole source in a cylindrical pipe containing area discontinuities
The sound power output of a monopole source in a cylindrical pipe containing area discontinuities Wenbo Duan, Ray Kirby To cite this version: Wenbo Duan, Ray Kirby. The sound power output of a monopole
More informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
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 informationPrediction of the radiated sound power from a fluid-loaded finite cylinder using the surface contribution method
Prediction of the radiated sound power from a fluid-loaded finite cylinder using the surface contribution method Daipei LIU 1 ; Herwig PETERS 1 ; Nicole KESSISSOGLOU 1 ; Steffen MARBURG 2 ; 1 School of
More informationFinal Exam Solution Dynamics :45 12:15. Problem 1 Bateau
Final Exam Solution Dynamics 2 191157140 31-01-2013 8:45 12:15 Problem 1 Bateau Bateau is a trapeze act by Cirque du Soleil in which artists perform aerial maneuvers on a boat shaped structure. The boat
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 informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
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 informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationSPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES
2010/2 PAGES 1 8 RECEIVED 21. 9. 2009 ACCEPTED 20. 1. 2010 Y. KOLEKOVÁ, M. PETRONIJEVIĆ, G. SCHMID SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES ABSTRACT
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationA study on regularization parameter choice in Near-field Acoustical Holography
Acoustics 8 Paris A study on regularization parameter choice in Near-field Acoustical Holography J. Gomes a and P.C. Hansen b a Brüel & Kjær Sound and Vibration Measurement A/S, Skodsborgvej 37, DK-285
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
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 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 informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
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 informationON THE MODELING AND SIMULATION OF AN ACOUSTIC CLOAK
5 th International Conference Advanced Composite Materials Engineering COMAT 014 16-17 October 014, Braşov, Romania ON THE MODELING AND SIMULATION OF AN ACOUSTIC CLOAK Veturia Chiroiu 1, Rodica Ioan 1,,
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
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 informationAn accelerated predictor-corrector scheme for 3D crack growth simulations
An accelerated predictor-corrector scheme for 3D crack growth simulations W. Weber 1 and G. Kuhn 2 1,2 1 Institute of Applied Mechanics, University of Erlangen-Nuremberg Egerlandstraße 5, 91058 Erlangen,
More informationBOUNDARY ELEMENT METHOD IN REFRACTIVE MEDIA
BOUDARY ELEMET METHOD I REFRACTIVE MEDIA David R. Bergman Exact Solution Scientific Consulting LLC, Morristown J, USA email: davidrbergman@essc-llc.com Boundary Element Method (BEM), an application of
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationVirtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing example
AMAS Workshop on Smart Materials and Structures SMART 03 (pp.313 324) Jadwisin, September 2-5, 2003 Virtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing
More informationTransactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,
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 informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationInstructions Please answer the five problems on your own paper. These are essay questions: you should write in complete sentences.
Instructions Please answer the five problems on your own paper. These are essay questions: you should write in complete sentences.. Recall that P 3 denotes the vector space of polynomials of degree less
More informationIdentification of Noise Sources by Means of Inverse Finite Element Method
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Identification of Noise Sources by Means of Inverse Finite Element Method M. Weber *,1, T. Kletschkowski 2 and B. Samtleben 3 1,2 Helmut-Schmidt-University
More informationFreund, Jouni; Salonen, Eero-Matti Sensitized principle of virtual work and the single-element strain energy test
Powered by CPDF (www.tcpdf.org) his is an electronic reprint of the original article. his reprint may differ from the original in pagination and typographic detail. Freund, Jouni; Salonen, Eero-Matti Sensitized
More informationLECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES
LECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES Figure 3.47 a. Two-mass, linear vibration system with spring connections. b. Free-body diagrams. c. Alternative free-body
More informationA NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 212) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 1-14, 212 A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE
More informationComparison of Galerkin and collocation Trefftz formulations for plane elasticity
Comparison of Galerkin and collocation Trefftz formulations for plane elasticity V.M.A. Leitão December 8, 2000 Abstract The purpose of this work is to compare and assess, more in terms of computational
More informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationICSV14 Cairns Australia 9-12 July, 2007
ICSV14 Cairns Australia 9-1 July, 7 STUDY ON THE CHARACTERISTICS OF VIBRATION AND ACOUSTIC RADIATION OF DAMAGED STIFFENED PANELS Hong Ming 1, Guo Xin-yi, Liu Lian-hai 3 and Li Gen-tian 4 1 Department of
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 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 informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationMatrix Representation
Matrix Representation Matrix Rep. Same basics as introduced already. Convenient method of working with vectors. Superposition Complete set of vectors can be used to express any other vector. Complete set
More informationStructural Dynamics A Graduate Course in Aerospace Engineering
Structural Dynamics A Graduate Course in Aerospace Engineering By: H. Ahmadian ahmadian@iust.ac.ir The Science and Art of Structural Dynamics What do all the followings have in common? > A sport-utility
More informationSub-structuring of mechanical systems based on the path concept
Sub-structuring of mechanical systems based on the path concept Francesc Xavier MAGRANS 1 ; Jordi POBLET-PUIG 2 ; Antonio RODRÍGUEZ-FERRAN 3 1 Ingeniería para el Control del Ruido S.L., Spain 1,2,3 Universitat
More informationTable of Contents. Preface... 13
Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...
More informationTechnische Universität Graz
Technische Universität Graz Robust boundary element domain decomposition solvers in acoustics O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2009/9 Technische Universität
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationAnalytical Solution for a Fluid-Structure Interaction Problem in Comparison with Finite Element Solution
Analytical Solution for a Fluid-Structure Interaction Problem in Comparison with Finite Element Solution Amirhossein Keivani & Ahmad Shooshtari Ferdowsi University of Mashhad, Mashhad, Iran. Ahmad Aftabi
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationStatic & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering
Static & Dynamic Analysis of Structures A Physical Approach With Emphasis on Earthquake Engineering Edward LWilson Professor Emeritus of Civil Engineering University of California, Berkeley Fourth Edition
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationMathematics Department Stanford University Math 61CM/DM Inner products
Mathematics Department Stanford University Math 61CM/DM Inner products Recall the definition of an inner product space; see Appendix A.8 of the textbook. Definition 1 An inner product space V is a vector
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationTRANSMISSION LOSS OF EXTRUDED ALUMINIUM PANELS WITH ORTHOTROPIC CORES
TRANSMISSION LOSS OF EXTRUDED ALUMINIUM PANELS WITH ORTHOTROPIC CORES PACS REFERENCE: 43.40-Rj RADIATION FROM VIBRATING STRUCTURES INTO FLUID MEDIA Names of the authors: Kohrs, Torsten; Petersson, Björn
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
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 informationTyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin
Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhauser Boston Basel Berlin Preface to the Fourth Edition Preface to the Third Edition
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationEfficient boundary element analysis of periodic sound scatterers
Boundary Element and Meshless Methods in Acoustics and Vibrations: Paper ICA2016-418 Efficient boundary element analysis of periodic sound scatterers M. Karimi, P. Croaker, N. Kessissoglou 1 School of
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
More informationBoundary. DIFFERENTIAL EQUATIONS with Fourier Series and. Value Problems APPLIED PARTIAL. Fifth Edition. Richard Haberman PEARSON
APPLIED PARTIAL DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fifth Edition Richard Haberman Southern Methodist University PEARSON Boston Columbus Indianapolis New York San Francisco
More informationEigenvalues and Eigenvectors A =
Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationA simple Galerkin boundary element method for three-dimensional crack problems in functionally graded materials
Materials Science Forum Vols. 492-493 (2005) pp 367-372 Online available since 2005/Aug/15 at www.scientific.net (2005) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/msf.492-493.367
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More information1 Introduction. Green s function notes 2018
Green s function notes 8 Introduction Back in the "formal" notes, we derived the potential in terms of the Green s function. Dirichlet problem: Equation (7) in "formal" notes is Φ () Z ( ) ( ) 3 Z Φ (
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationA geometric proof of the spectral theorem for real symmetric matrices
0 0 0 A geometric proof of the spectral theorem for real symmetric matrices Robert Sachs Department of Mathematical Sciences George Mason University Fairfax, Virginia 22030 rsachs@gmu.edu January 6, 2011
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationAcoustic Radiation Modes of a Tire on a Reflecting Surface
Purdue University Purdue e-pubs Publications of the Ray W. Herrick Laboratories School of Mechanical Engineering 3-2005 Acoustic Radiation Modes of a Tire on a Reflecting Surface Kiho Yum Purdue University
More informationBoundary element methods in the prediction of the acoustic damping of ship whipping vibrations
ANZIAM J. 45 (E) ppc845 C856, 2004 C845 Boundary element methods in the prediction of the acoustic damping of ship whipping vibrations D. S. Holloway G. A. Thomas M. R. Davis (Received 8 August 2003) Abstract
More informationTopic 1: Matrix diagonalization
Topic : Matrix diagonalization Review of Matrices and Determinants Definition A matrix is a rectangular array of real numbers a a a m a A = a a m a n a n a nm The matrix is said to be of order n m if it
More informationarxiv: v1 [math.ap] 21 Dec 2018
Uniqueness to Inverse Acoustic and Electromagnetic Scattering From Locally Perturbed Rough Surfaces Yu Zhao, Guanghui Hu, Baoqiang Yan arxiv:1812.09009v1 [math.ap] 21 Dec 2018 Abstract In this paper, we
More informationwhich implies that we can take solutions which are simultaneous eigen functions of
Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationAn 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 informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More information