Abstract Simulations of sound pressure levels are often used as a decision basis when planning location for e.g wind turbines. In previous work, an e
|
|
- Brittney Mitchell
- 5 years ago
- Views:
Transcription
1
2 Abstract Simulations of sound pressure levels are often used as a decision basis when planning location for e.g wind turbines. In previous work, an e cient model for computing sound propagation over irregular terrain has been developed. The model simulates the propagation using the time dependent wave equation. The purpose of this study was to improve this solver by finding a way to represent atmospheric attenuation in the time dependent wave equation by adding a damping term to the equation. This was done by numerically simulating wave propagation and comparing to measured data on the atmospheric attenuation. The study did not result in finding a satisfying way to model the attenuation. The results of the simulations are analyzed and suggestions for further research are given.
3 Contents 1 Introduction Introduction Atmospheric Attenuation Project definition Analysis The wave equation The SBP-SAT method/the semi discrete model Results Simulations using damping term u t Simulations using damping term u t and r 2 u t Simulations using damping term u t r 2 u t Conclusions and discussion 11
4 1 Introduction 1.1 Introduction Calculations of sound propagation are often used as decision basis when planning the construction of for example wind turbines, airports, tra c junctions and other places where extreme sound pressure levels occur and cause problems. This makes accurate modeling of sound propagation an important subject. To predict sound pressure levels there are several factors that should be included in the model. For example, the sound pressure levels depends strongly on meteorological data such as sound speed variations and atmospheric attenuation as well as the topography of the area over which the sound propagates [7]. There are several models that can be used to simulate sound propagation, varying in computational cost and accuracy. A common model used in Sweden today is NORD2000, based on ray tracing theory. Ray tracing theory models do not take complicated topography or meteorology into account [5]. More accurate models such as the finite di erence time domain method (FDTD) are often too time consuming to be used in three dimensional calculations. A good model for sound propagation is the time-dependent acoustic wave equation: u tt = c 2 u (1) Solving the wave equation can result in high computational cost if not using an e cient numerical method. The computation time can be reduced by using the Helmholtz equation instead. Helmholtz equation is a time independent, frequency domain representation of the wave equation. Transforming the wave equation to frequency domain results in less computational cost, but modeling in frequency domain demands sources consisting of a single frequency. Wind turbines are broadband sources and require a time-dependent model for e cient simulation of sound pressure levels. In previous work, an e cient numerical solver for solving wave propagation problems on large domains has been implemented [3]. This solver uses fourth order accurate finite di erence methods and runs in parallel in MATLAB. The solver can be used to model propagation over complicated terrains, but does not take into account meteorological data such as atmospheric attenuation. 1
5 1.2 Atmospheric Attenuation Atmospheric attenuation is a damping of the sound waves caused by transport processes and molecular relaxtion losses in the atmosphere. As sound waves propagate through the atmosphere the amplitude decreases due to absorption as p t = p i e s, (2) where p t is the sound pressure amplitude at a distance s from the initial pressure p i and is the attenuation coe cient dependent on frequency, background pressure, humidity and temperaure. More detailed information about atmospheric attenuation and measured values for for varying atmospherical data is given in [2]. Figure 1: Example of table with values on the attenuation constant for varying atmospherical data as given in [2]. 1.3 Project definition In this project we have studied the possibility to use a damping term to model the atmospheric attenuation in the time dependent wave equation. This would be a way to improve the model developed in [3]. The idea was to add a physical damping term to the right hand side of (1) of the form r ru t,wherethecoe cient is a function of the atmosphere, i.e. dependent on temperature, humidity and pressure, u tt = c 2 u + r ru t. (3) 2
6 The damping term should depend on frequency correspondingly to the frequency dependence in the atmospheric attenuation. Ideally, this would lead to a model that could accurately simulate the atmospheric attenuation of broad band sources, which as discussed above is desired. The main object of the project was, thus, to study the wave equation with the added damping term with the aim to find an explicit form of the damping term. This is done by simulating the wave propagation and testing it against -values from [2] inserted in the analytic solution given by (2). To begin with, a fixed atmosphere (i.e fixed values on humidity, temperature and background pressure) was examined. The aim was to later extend the modeling to more general atmospheres. In the following sections we first describe the numerical model and how the sound propagation was simulated. We then present the results of our simulations, followed by a discussion about the results and recommendations for further research. 2 Analysis 2.1 The wave equation In order to investigate the behavior of the damping term we perform simulations of damped waves and compare them to analytical data. Sound from a point source propagates with spherical symmetry in free space. Using spherical coordinates in one dimension will therefore lead to a simplified model that will still simulate e ects relevant to propagation in three dimensions. Adding the damping term as described above to the one-dimensional wave equation in spherical coordinates and setting in (3) to be independent of r yields r 2 u tt = c r r By numerically modeling this wave and comparing it to the analytical attenuated wave equation the damping term can be adjusted to fit the analytical damping. The analytical form of the damped wave is given by applying the damping described in (2) to the analytical solution of the one dimensional spherical wave equation u d (r, t) = a r cos(!t! c r) e r (5) 2.2 The SBP-SAT method/the semi discrete model The SBP-SAT method is a finite di erence method that combines Summation By Part operators (SBP) with the simultaneous approximation terms 3
7 method (SAT) to weakly impose boundary conditions. Proving stability for SBP-SAT representations is done by analyzing the energy using the semi discrete energy method as demonstrated later in this section. In order to present the SBP-SAT method we first introduce some notation and definitions (following [4]). Discretization The domain r 0 apple r apple r N is discretized using M points given by r i = r 0 + ih, h = r N r 0, i =1, 2,...,M. M 1 The numeric solution at r i is denoted by v i and the corresponding analytic solution is denoted by u(r i ). These vectors are used when introducing the simultaneous approximation terms ê 1 =[1, 0,..., 0] T, ê N =[0,..., 0, 0] T. Definition 1 A di erence operator D 1 = H 1 Q is a first derivative SBP operator if H = H T > 0, and Q + QT = B = diag( 1, 0..., 0, 1). Definition 2 The operator D (b) 2 = H 1 ( M (b) + BS) where b(r) > 0, H is diagonal and positive definite, M (b) is symmetric and positive semidefinite, S approximates the first derivative operator at the boundaries and B = diag( b 0, 0,...,0,b N ),. Definition 3 An inner product for discrete real-valued vector functions v, w 2 R M is given by (v, w) H = v T Hw where H = H T > 0, with corresponding norm v 2 H = vt Hv. The continuous term u r, for example, would have the semi-discrete approximation D 1 v using the first derivative SBP operator introduced above. Below we apply the SBP-SAT method to the wave equation given in spherical coordinates. A nonreflecting boundary conditions for spherical coordinates as described by Grote in [1] is applied on the right boundary in order to prevent non-physical reflections. The left (inflow) boundary is implemented by strongly enforcing the sound pressure levels, sometimes referred to as injection: 4
8 8 r 2 u tt = c r r 0 apple r apple r N, t > u r + 1 c u t + u r =0, r = r N >: u(r, t) = a r cos(!t! c r), r = r 0 (6) (7) (8) 1 Analytic injection Non-reflective boundary Amplitude 0-1 Computational domain 0 r 0 r N >r N Distance Figure 2: The simulation set up. The actual computational domain is between r 0 and r N. The semi-discrete approximation of (6) using the second order SBP operator is given by Bv tt = c 2 D (b) 2 v (9) where B = diag(b 0,b 2 1,...,b2 N ) and b(r) =r2. The nonreflecting boundary condition is imposed weakly by adding a simultaneous approximation term to the right hand side of the equation Bv tt = c 2 D (b) 2 v H 1 ê N {(Sv + v t c + v r ) N} (10) The SAT method can be seen as a force pulling the solution on the boundary towards the boundary condition. To prove stability we use the semi discrete energy method, given by multiplying the equation by v T t H and adding the transpose. This yields d dt ( v t 2 HB + c 2 v T M (b) v)=c 2 v T t BSv + c 2 (v T t BSv) T 2 (v T t ) N {(Sv + v t c + v r ) N} (11) 5
9 Where we have used that v T t HBv + v T HBv t = d dt vt HBv = d dt v HB. The semi discrete energy estimate for the wave equation is given by E = v t 2 HB +c2 v T M (b) v, which gives that the left hand side in (11) corresponds to de/dt. d dt E = c2 v T t BSv + c 2 (v T t BSv) T 2 (v T t ) N {(Sv + v t c + v r ) N} (12) The right hand side is further simplified using the definition of B and S d dt E = 2c2 b 0 (vt T ) N (Sv) N + 2(c 2 b N )(vt T ) N (Sv) N 2 2 (vt T ) N v N (13) c r N To prove stability we must have de/dt apple 0. In order to cancel out the term 2c 2 b N (vt T ) N (Sv) N that could lead to non-physical growth in energy we must set the SAT coe cient to = b N c 2 = r 2 Nc 2. (14) Adding the damping term u t to the equation introduces a second SAT term to the semi discrete representation. By the same method as above the penalty term is set to 2 = rn 2, and the semi discrete approximation is given by 3 Results Bv tt = c 2 D (b) 2 v + D(b) 2 v t r 2 Nc 2 H 1 ê N {(Sv + v t c + v r ) N} r 2 NH 1 ê N {(Sv t + v tt c + v t r ) N} (15) In this section we present results of simulations for a fixed atmosphere, 50% relative humidity and 15 C (see Table 1 in [2] for exact values of ). All the testing is done against -values for this specific atmosphere. The initial intention was to further extend the study to more general atmospheric data, but this was omitted given the nature of the results obtained. The simulations began with trying to adapt the damping term u t to the atmospheric attenuation data for the given atmosphere. It showed that it was possible to mimic the attenuation for a specific frequency by choosing a suitable value on, but that the damping term and the attenuation did not follow the same frequency dependence. We then moved on to testing other damping operators and combinations of operators. Below the results from these simulations are presented. 6
10 3.1 Simulations using damping term u t The simulations were started by trying to find the correct damping for a fixed frequency. This was possible, which is shown in Figure 3 where the maximum values of the amplitude for both the analytical and the numerical solution are shown for a damped wave with frequency 1000 Hz. The numerical and analytical solution follows the same decay Amplitude for frequency 1000 Hz over distance Analytic Numeric: β(1000 Hz) 0.7 Amplitude Distance (m) Figure 3: The amplitude of a signal with a frequency of 1000 Hz with optimized for 1000 Hz coincide with the measured values of the attenuation. However, it was found that signals with di erent frequencies demand different values of. The damping operator u t a ects signals with higher frequencies much stronger than the atmospheric attenuation does, which results in that a lower value on is needed for higher frequencies in order to fit the damping. The optimal -values for frequencies ranging between 50 and 1000 Hz are shown in Figure 4. Since di erent frequencies demand different values, a signal containing several frequencies can not be modelled using just one damping term as it would only give the correct damping for one frequency. This is illustrated in Figure 5 and Figure 6. An incorrect results in an incorrect damping and an incorrect modeling of the sound. 7
11 β Frequency (Hz) Figure 4: The value of as a function of frequency. Since di erent frequencies demand di erent values on, this damping term can not be used to model signals containing multiple frequencies Amplitude for frequency 1000 Hz over distance Analytic Numeric: β(1000 Hz) Numeric: β(50 Hz) Amplitude Distance (m) Figure 5: The amplitude of a signal with a frequency of 1000 Hz with optimized for 50 Hz compared to the analytical solution. 8
12 Amplitude ratio β(50 Hz) β(250 Hz) β(1000 Hz) Distance (m) Figure 6: Amplitude ratio between the numerical and the analytical solution. A signal with frequency 250 Hz is modelled using three values of. An incorrect value of leads to a increasing error as the sound wave propagates. 3.2 Simulations using damping term u t and r 2 u t Proceeding our project, we tried a damping term of a di erent form, namely u t, resulting in the equation u tt = c 2 u u t. (16) As in the previous subsection, we started out by trying to find the correct damping for a signal of fixed frequency. This did not succeed using the damping term u t. The numerical solution did not follow the same exponential decay as the analytical solution. In Figure 7, the amplitude ratio between the numerical and the analytical solution is shown. What is desired is a straight line, meaning that the numerical solution follows the same damping relation as the analytical. In the case with the damping term u t this was not observed. We found that we needed to add a factor of r 2 to the damping term in order to obtain the correct damping, u tt = c 2 u r 2 u t (17) As seen in Figure 7, the r 2 u t -term gives a decay coinciding with the analytical solution. However, this damping did not give better results than the 9
13 u t -term when simulating varying frequencies. Here too it was found that di erent frequencies demand di erent values of (see Figure 7 and Figure 8) γ Undamped γ r 2 Amplitude ratio Distance (m) Figure 7: The amplitude ratio between the numerical and analytical solution for a 250 Hz signal. Using a damping term u t, the numerical solution did not follow the same damping as the analytical solution. Using r 2 u t it was possible to obtain a correct damping. The undamped solution is included for comparison γ Frequency (Hz) Figure 8: The value of as a function of frequency. As in the previous simulations it was found that di erent frequencies demanded di erent values on. 10
14 3.3 Simulations using damping term u t r 2 u t As the damping was stronger for high frequencies and the r 2 damping was stronger for low frequencies a combination of the two was tested, namely u tt = c 2 u + u t r 2 u t. (18) The idea was that a combination of the two damping terms could result in a fair approximation of the damping ranging over various frequencies. The two damping terms were adjusted to optimize the damping for two di erent frequencies, 50 Hz and 1000 Hz, i.e the values of and were chosen to give correct damping for signals of 50 Hz and 1000 Hz. When this setup was tested against frequencies lying between 50 Hz and 1000 Hz, it did not result in a correct damping. This is shown in Figure 9, for the amplitude ratios for frequencies 50, 250 and 1000 Hz Hz 250 Hz 1000 Hz 250 Hz Undamped Amplitude ratio Distance (m) Figure 9: The amplitude ratio for the frequencies 50, 250 and 1000 Hz using the modeling explained above. Satisfactory damping was obtained for signals near 50 Hz and 1000 Hz, but not for signals in between. The amplitude ratio for an undamped signal is shown as a reference. 4 Conclusions and discussion It was proved hard to find a damping term able to model the atmospheric attenuation by using the damping terms investigated above. The aim was to find a way to model signals containing multiple frequencies, but no satisfactory way to do this was found. The reason for this was that the atmospheric attenuation followed a more complicated frequency dependence than expected. A closer look at the physics describing the attenuation reveals some 11
15 useful information about this. The ISO-paper [2] describes the attenuation as a sum of di erent types of absorption derived to several physical mechanisms having distinct frequency dependencies: Classical ( cl ) and rotational absorption ( rot ) caused by transport processes described by classical mechanics and molecular rotational relaxation. The classical absorption is caused by compression and expansion of the air leading to friction between particles [6]. These types of absorptions follow a quadratic frequency dependence given by cr = cl + rot = (T/T 0 ) 1/2 f 2 p a /p r, where T 0 and p r represent reference temperature and pressure. Vibrational relaxation damping caused by vibrational relaxation of oxygen( vib,o ) and nitrogen( vib,n ).The attenuation caused by these factors follow a more complicated frequency dependence given by vib,o = 2[( ) max] cf ro f 2 [1 + (f/f ro ) 2 ] 1 vib,n = 2[( ) max] cf rn f 2 [1 + (f/f rn ) 2 ] where the terms f r represent relaxation frequencies. The relaxation frequencies are properties of the nitrogen and oxygen, dependent on temperature and humidity of the air. Based on these factors the frequency dependence in the vibrational relaxation damping can change significantly (see figure 10 and figure 11). The total attenuation is thus given by = cr + vib,o + vib,n.sincethese di erent types of attenuation are derived from di erent physical mechanisms and follow varying frequency dependence, it seems logical to separate them into di erent damping terms in the time dependent wave equation as well. As can be seen in Figure 12, the frequency dependence in varies as the properties of the air changes. As the damping terms we investigated all have a fixed frequency dependence it is easy to see that they will not be suitable to model the attenuation as it will not be able to represent these air property changes. 1 12
16 Vibational relaxational damping Humidity: 10% Humidity: 50% Humidity: 90% Frequency (Hz) Figure 10: Vibrational relaxation damping vib,o caused by oxygen for three di erent atmospheres. The frequency dependence varies significantly. Vibrational relaxation damping 8 x Humidity: 10% Humidity: 50% Humidity: 90% Frequency (Hz) Figure 11: Vibrational relaxation damping vib,n caused by nitrogen for three di erent atmospheres. The frequency dependence varies significantly. Total damping Humidity: 10% Humidity: 50% Humidity: 90% Frequency (Hz) Figure 12: Total damping for three di erent atmospheres. Since the damping is caused by several factors it is di cult to model it using a single damping term. 13
17 Further research should include trying to find damping terms corresponding to the di erent types of attenuation described above. By finding representation for all three terms in the time dependent wave equation it would be easier to make the model suitable for varying atmospheric data. 14
18 References [1] M. J. Grote, Nonreflection Boundary Conditions for Time Dependent Wave Propagation, Research Report , Seminar für Angewandte Mathematik, Zürich, Switzerland, April [2] ISO b, Acoustics - Attenuation of sound during propagation outdoors. [3] Martin Almquist, Numerical wave propagation in large-scale 3-d environments, [4] Ken Mattsson, Frank Ham, Gianluca Iaccarino, Stable and accurate wave-propagation in discontinuous media, Journal of Computional Physics, 277: , [5] DELTA (Danish Electronics, Lights & Acoustics), Nord2000. Comprehensive Outdoor Sound Propagation Model. Part 2: Propagation in an Atmosphere with Refraction, Journal no. AV 1851/00, Project no. A550054, 31 December 2001, Revised 31 March [6] Gustav Myhrman, The impact on tra c voice by roundabouts, ISSN , December [7] Conny Larsson, Olof Öhlund, Variations of sound from wind turbines during di erent weather conditions. 15
Efficient wave propagation on complex domains
Center for Turbulence Research Annual Research Briefs 2006 223 Efficient wave propagation on complex domains By K. Mattsson, F. Ham AND G. Iaccarino 1. Motivation and objectives In many applications, such
More informationA Stable and Accurate Davies-like Relaxation Procedure using Multiple Penalty Terms for Lateral Boundary Conditions
A Stable and Accurate Davies-like Relaxation Procedure using Multiple Penalty Terms for Lateral Boundary Conditions Hannes Frenander Division of Computational Mathematics, Department of Mathematics, Linköping
More informationAssessment of noise prediction models for long-range sound propagation of wind turbines over water
Assessment of noise prediction models for long-range sound propagation of wind turbines over water Lukas Mylonas, Bahri Uzunoglu lukas.mylonas@gmail.com, bahri.uzunoglu@angstrom.uu.se Uppsala University
More informationNumerical Model for Prediction of Sound Propagation Emitted from Wind Turbines
Numerical Model for Prediction of Sound Propagation Emitted from Wind Turbines Jasmin Hörmeyer Institute of Structural Analysis, Leibniz Universität Hannover, Germany. Raimund Rolfes Institute of Structural
More informationA New Finite-Difference Method for General Long-Range Rotorcraft Acoustics: Initial Comparisons with Intermediate- Range Data
A New Finite-Difference Method for General Long-Range Rotorcraft Acoustics: Initial Comparisons with Intermediate- Range Data Subhashini Chitta, John Steinhoff (B.H. Goethert Professor, UTSI), Andrew Wilson
More informationStable boundary treatment for the wave equation on second-order form
Stable boundary treatment for the wave equation on second-order form Ken Mattsson Frank Ham Gianluca Iaccarino June 24, 2008 Abstract A stable and accurate boundary treatment is derived for the secondorder
More informationNUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY
NUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY Stefan Johansson Department of Information Technology Scientific Computing
More informationWell-posedness, stability and conservation for a discontinuous interface problem: an initial investigation.
Well-posedness, stability and conservation for a discontinuous interface problem: an initial investigation. Cristina La Cognata and Jan Nordström Abstract A robust interface treatment for the discontinuous
More informationISO 354 INTERNATIONAL STANDARD. Acoustics Measurement of sound absorption in a reverberation room
INTERNATIONAL STANDARD ISO 354 Second edition 2003-05-15 Acoustics Measurement of sound absorption in a reverberation room Acoustique Mesurage de l'absorption acoustique en salle réverbérante Reference
More informationEfficient outdoor sound propagation modelling in time-domain
Efficient outdoor sound propagation modelling in time-domain Timothy Van Renterghem Acoustics Group Department of Information Technology Ghent University/iMinds CEAS, X-Noise : Atm. and ground effects
More informationRECENT DEVELOPMENTS IN COUPLING TOPOGRAPHICAL AND METEOROLOGICAL EFFECTS WITH THE GREEN'S FUNCTION PARABOLIC EQUATION (GFPE): THEORY AND EXPERIMENTS
Twelfth International Congress on Sound and Vibration RECENT DEVELOPMENTS IN COUPLING TOPOGRAPHICAL AND METEOROLOGICAL EFFECTS WITH THE GREEN'S FUNCTION PARABOLIC EQUATION (GFPE): THEORY AND EXPERIMENTS
More informationDiagonal-norm upwind SBP operators
Diagonal-norm upwind SBP operators Ken Mattsson June 8, 16 Abstract High-order accurate first derivative finite difference operators are derived that naturally introduce artificial dissipation. The boundary
More informationSpectral analysis of the incompressible Navier-Stokes equations with different boundary conditions
Spectral analysis of the incompressible Navier-Stokes equations with different boundary conditions Cristina La Cognata, Jan Nordström Department of Mathematics, Computational Mathematics, Linköping University,
More informationNonlinear parabolic equation model for finite-amplitude sound propagation in an inhomogeneous medium over a non-flat, finite-impedance ground surface
Nonlinear parabolic equation model for finite-amplitude sound propagation in an inhomogeneous medium over a non-flat, finite-impedance ground surface T. Leissing a, P. A H Jean a, J. Defrance a and C.
More informationNoise in enclosed spaces. Phil Joseph
Noise in enclosed spaces Phil Joseph MODES OF A CLOSED PIPE A 1 A x = 0 x = L Consider a pipe with a rigid termination at x = 0 and x = L. The particle velocity must be zero at both ends. Acoustic resonances
More informationTHE ACOUSTIC SOURCE STRENGTH OF HIGH-ENERGY BLAST WAVES: COMBINING MEASUREMENTS AND A NON-LINEAR MODEL
Proceedings of 20 th International Congress on Acoustics, ICA 2010 23-27 August 2010, Sydney, Australia THE ACOUSTIC SOURCE STRENGTH OF HIGH-ENERGY BLAST WAVES: COMBINING MEASUREMENTS AND A NON-LINEAR
More informationUsing windscreens to improve the efficiency of noise barriers in wind: finite-difference time-domain simulations
Using windscreens to improve the efficiency of noise barriers in wind: finite-difference time-domain simulations PACS REFERENCE: 43.8.Js, 43.8.Gq, 43.5.Gf Van Renterghem Timothy; Botteldooren Dick University
More informationA simple model for estimating excess attenuation of road traffic noise
A simple model for estimating excess attenuation of road traffic noise Kohei Yamamoto and Mitsuyasu Yamashita Kobayasi Institute of Physical Research, 3-20-41, Higashimotomachi, Kokubunji, 185 Japan (Received
More informationErik Sloth Vestas. Niels Christian Møller Nielsen VESTAS Ejler Kristensen BONUS Energy Bo Søndergaard DELTA
Problems related to the use of the existing noise measurement standards when predicting noise from wind turbines and wind farms. Erik Sloth Vestas Niels Christian Møller Nielsen VESTAS Ejler Kristensen
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
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 informationIn situ measurement methods for characterising sound diffusion
Proceedings of the International Symposium on Room Acoustics, ISRA 9 August, Melbourne, Australia In situ measurement methods for characterising sound diffusion I. Schmich (), N. Brousse () () Université
More informationEstimating Community Sound Levels of Large Industrial Equipment
Estimating Community Sound Levels of Large Industrial Equipment Paper # 2001-603 Prepared by David J. Parzych Power Acoustics, Inc. 12472 Lake Underhill Rd #302, Orlando, FL 32828 ABSTRACT A simple acoustic
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationPrediction of Sound Propagation From Power Transmission Plant
Prediction of Sound Propagation From Power Transmission Plant Jingchao Sun Stockholm, 2013 Thesis for the degree of Master of Science, 30 Hp Department of Sound and Vibration The Marcus Wallenberg Laboratory
More informationStability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws
Stability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws G. Bastin http://perso.uclouvain.be/georges.bastin/ Georges Bastin Jean Michel Coron Progress in Nonlinear
More informationStable shock formation near plane-symmetry for nonlinear wave equations
Stable shock formation near plane-symmetry for nonlinear wave equations Willie WY Wong wongwwy@member.ams.org Michigan State University University of Minnesota 9 March 2016 1 Based on: arxiv:1601.01303,
More information1 Wind Turbine Acoustics. Wind turbines generate sound by both mechanical and aerodynamic
Wind Turbine Acoustics 1 1 Wind Turbine Acoustics Wind turbines generate sound by both mechanical and aerodynamic sources. Sound remains an important criterion used in the siting of wind farms. Sound emission
More informationON SOUND POWER MEASUREMENT OF THE ENGINE IN ANECHOIC ROOM WITH IMPERFECTIONS
ON SOUND POWER MEASUREMENT OF THE ENGINE IN ANECHOIC ROOM WITH IMPERFECTIONS Mehdi Mehrgou 1, Ola Jönsson 2, and Leping Feng 3 1 AVL List Gmbh., Hans-list-platz 1,8020, Graz, Austria 2 Scania CV AB, Södertälje,
More informationA Provable Stable and Accurate Davies-like Relaxation Procedure Using Multiple Penalty Terms for Lateral Boundaries in Weather Prediction
Department of Mathematics A Provable Stable and Accurate Davies-like Relaxation Procedure Using Multiple Penalty Terms for Lateral Boundaries in Weather Prediction Hannes Frenander and Jan Nordström LiTH-MAT-R--24/9--SE
More informationStandard ECMA-108 3rd Edition - December Measurement of High-frequency Noise emitted by Information Technology and Telecommunications Equipment
Standard ECMA-108 3rd Edition - December 1996 Standardizing Information and Communication Systems Measurement of High-frequency Noise emitted by Information Technology and Telecommunications Equipment
More informationAdapting a propeller noise model for aircraft at cruising altitudes
Adapting a propeller noise model for aircraft at cruising altitudes David M. BLUNT 1 ; Adrian JONES 2 ; David MEWETT 3 1 Aerospace Division, Defence Science and Technology Organisation, Australia 2 Maritime
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 007 Numerical simulation of low level sonic boom propagation through random inhomogeneous sound speed fields PACS: 43.8.Mw Ollivier, Sébastien
More informationVibrations and waves: revision. Martin Dove Queen Mary University of London
Vibrations and waves: revision Martin Dove Queen Mary University of London Form of the examination Part A = 50%, 10 short questions, no options Part B = 50%, Answer questions from a choice of 4 Total exam
More informationInsertion Loss Analysis of the Acoustic Panels with Composite Construction
ANALELE UNIVERSITĂłII EFTIMIE MURGU REŞIłA ANUL XX, NR., 13, ISSN 1453-7397 Vasile Ovidiu Insertion Loss Analysis of the Acoustic Panels with Composite Construction In order to reduce noise pollution,
More informationAn integrated approach to Acoustics - sound propagation and the phenomena affecting it
An integrated approach to Acoustics - sound propagation and the phenomena affecting it IOA LONDON BRANCH EVENING MEETING, WSP PARSONS BRINCKERHOFF PRESENTATION BY PANOS ECONOMOU OF PEMARD, 19 TH APRIL
More informationStable and high-order accurate finite difference schemes on singular grids
Center for Turbulence Research Annual Research Briefs 006 197 Stable and high-order accurate finite difference schemes on singular grids By M. Svärd AND E. van der Weide 1. Motivation and objectives The
More informationPrediction of environmental sound pressure levels due to large wind turbines
Wind Farm Noise: Paper ICA2016-168 Prediction of environmental sound pressure levels due to large wind turbines Matteo Deambrosi Papini (a), Matías Hernández Castellani (b), Alice Elizabeth González (c),
More informationComparison of Noise Test Codes when Applied to air Compressors
VDI-Berichte Nr. 1932, 2006 A 6 83 Comparison of Noise Test Codes when Applied to air Compressors Michael J. Lucas, INCE Bd. Cert., Ingersoll-Rand Company, Davidson/NC INTRODUCTION Beginning January 2004,
More informationOpen quantum systems
Chapter 4 Open quantum systems 4. Quantum operations Let s go back for a second to the basic postulates of quantum mechanics. Recall that when we first establish the theory, we begin by postulating that
More informationSUMMATION-BY-PARTS IN TIME: THE SECOND DERIVATIVE
Department of Mathematics SUMMATION-BY-PARTS IN TIME: THE SECOND DERIVATIVE Jan Nordström and Tomas Lundquist LiTH-MAT-R--2014/11--SE Department of Mathematics Linköping University S-581 83 Linköping SUMMATION-BY-PARTS
More informationCalculation methods for the physical properties of air used in the calibration of microphones
Downloaded from orbit.dtu.dk on: Jun 15, 018 Calculation methods for the physical properties of air used in the calibration of microphones Rasmussen, Knud Publication date: 1997 Document Version Publisher's
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 informationDubai Municipality LRT Noise and Vibration Impact Stage 2
Dubai Municipality LRT Noise and Vibration Impact Stage 2 Environmental Impact 11 th of April, 2004 Environmental Impact Report According to the PS007 Agreement, Environmental Impact Report Based on Noise
More informationSound Propagation in the Nocturnal Boundary Layer. Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert
Sound Propagation in the Nocturnal Boundary Layer Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert The Propagation of Sound Outdoors (over flat ground) The atmosphere is a gas under the influence
More informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationROOM RESONANCES USING WAVE BASED GEOMET- RICAL ACOUSTICS (WBGA)
ROOM RESONANCES USING WAVE BASED GEOMET- RICAL ACOUSTICS (WBGA) Panos Economou, Panagiotis Charalampous P.E. Mediterranean Acoustics Research & Development Ltd, Cyprus email: panos@pemard.com Geometrical
More informationProblem Set 9: Solutions
University of Alabama Department of Physics and Astronomy PH 53 / LeClair Fall 1 Problem Set 9: Solutions 1. Energetics of diatomic systems I. An approximate expression for the potential energy of two
More informationTHE ACOUSTIC IMPEDANCE MEASUREMNET SYSTEM USING TWO MICROPHONES
P-7 THE ACOUSTIC IMPEDANCE MEASUREMNET SYSTEM USING TWO MICROPHONES RYU, YUNSEON BRUEL & KJAER SOUND & VIBRATION MEASUREMENT A/S SKODSBORGVEJ 307 NAERUM 2850 DENMARK TEL : +45 77 41 23 87 FAX : +45 77
More informationCalculation Of Fields On The Beam Axis Of The Delta Undulator
LCLS-TN-3- Calculation Of Fields On The Beam Axis Of The Delta Undulator Zachary Wolf SLAC December 8, 203 Abstract It is not possible to do Hall probe measurements on the beam axis of the Delta undulator
More informationAMATH 351 Mar 15, 2013 FINAL REVIEW. Instructor: Jiri Najemnik
AMATH 351 Mar 15, 013 FINAL REVIEW Instructor: Jiri Najemni ABOUT GRADES Scores I have so far will be posted on the website today sorted by the student number HW4 & Exam will be added early next wee Let
More informationNUMERICAL SOLUTION OF DOPANT DIFFUSION EQUATIONS. M.J. Bainss, C.P. Please and P.K. Sweby. Mathematics Department University of Reading, Berkshire.
271 NUMERICAL SOLUTION OF DOPANT DIFFUSION EQUATIONS M.J. Bainss, C.P. Please and P.K. Sweby. Mathematics Department University of Reading, Berkshire. 272 SUMMARY A description is given of the application
More information1 4 r q. (9.172) G(r, q) =
FOURIER ANALYSIS: LECTURE 5 9 Green s functions 9. Response to an impulse We have spent some time so far in applying Fourier methods to solution of di erential equations such as the damped oscillator.
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationNumerical Methods I: Polynomial Interpolation
1/31 Numerical Methods I: Polynomial Interpolation Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 16, 2017 lassical polynomial interpolation Given f i := f(t i ), i =0,...,n, we would
More informationsingle uniform density, but has a step change in density at x = 0, with the string essentially y(x, t) =A sin(!t k 1 x), (5.1)
Chapter 5 Waves II 5.1 Reflection & Transmission of waves et us now consider what happens to a wave travelling along a string which no longer has a single uniform density, but has a step change in density
More informationInfluence of Atmospheric Conditions on Sound Propagation - Mathematical Modeling
Óbuda University e-bulletin Vol. 5, No. 1, 2015 Influence of Atmospheric Conditions on Sound Propagation - Mathematical Modeling Pavol Liptai, Miroslav Badida, Katarína Lukáčová Department of Process and
More informationContents. Page 1 of 12
Contents Calculation method for road traffic noise Description of the calculation method NMPB General comment on the similarities and differences of NMPB with the END and on the possible use as Interim
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 informationDELTA Test Report. DANAK TEST Reg. no Measurement of Sound Absorption Coefficient for Kvadrat Soft Cells Wall Panel Type Time
We help ideas meet the real world DELTA Test Report DANAK TEST Reg. no. 100 Measurement of Sound Absorption Coefficient for Kvadrat Soft Cells Wall Panel Type Time Client: Kvadrat Soft Cells A/S Page 1
More informationGlasses and Jamming: lessons from hard disks in a narrow channel
Glasses and Jamming: lessons from hard disks in a narrow channel Mike Moore School of Physics and Astronomy, University of Manchester, UK June 2014 Co-Author: Mike Godfrey, University of Manchester Mike
More informationinter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE
Copyright SFA - InterNoise 2000 1 inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering 27-30 August 2000, Nice, FRANCE I-INCE Classification: 0.0 BACKSCATTERING
More informationPhysics 115 Lecture 20. Reflection and Reverberation March 9, 2018
Physics 115 Lecture 20 Reflection and Reverberation March 9, 2018 1 Attenuation by absorption The wave intensity decreases due to internal friction among the molecules of the medium Attenuation is described
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More informationChalmers Publication Library
Chalmers Publication Library Thick barrier noise-reduction in the presence of atmospheric turbulence: Measurements and numerical modelling This document has been downloaded from Chalmers Publication Library
More informationOpen quantum systems
Chapter 4 Open quantum systems 4.1 Quantum operations Let s go bac for a second to the basic postulates of quantum mechanics. Recall that when we first establish the theory, we begin by postulating that
More informationDELTA Test Report. Measurement of Sound Absorption for AqFlex ON, one Sample as a discrete Object. Performed for Flex Acoustics
DELTA Test Report TEST Reg. no. 100 Measurement of Sound Absorption for AqFlex ON, one Sample as a discrete Object Performed for Flex Acoustics Project no.: T202611 Page 1 of 7 21 January 2013 DELTA Venlighedsvej
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
FINITE DIFFERENCE MODELING OF INFRASOUND PROPAGATION TO LOCAL AND REGIONAL DISTANCES Catherine de Groot-Hedlin Scripps Institution of Oceanography, University of California at San Diego Sponsored by Army
More informationOscillations and Waves
Oscillations and Waves Oscillation: Wave: Examples of oscillations: 1. mass on spring (eg. bungee jumping) 2. pendulum (eg. swing) 3. object bobbing in water (eg. buoy, boat) 4. vibrating cantilever (eg.
More informationMath 413/513 Chapter 6 (from Friedberg, Insel, & Spence)
Math 413/513 Chapter 6 (from Friedberg, Insel, & Spence) David Glickenstein December 7, 2015 1 Inner product spaces In this chapter, we will only consider the elds R and C. De nition 1 Let V be a vector
More informationCalculation Of Fields On The Beam Axis Of The Delta Undulator
LCLS-TN-3- Calculation Of Fields On The Beam Axis Of The Delta Undulator Zachary Wolf SLAC April, 20 Abstract It is not possible to do Hall probe measurements on the beam axis of the Delta undulator since
More informationPART VIII: ABSORPTIVE SILENCER DESIGN
PART VIII: ABSORPTIVE SILENCER DESIGN Elden F. Ray June 10, 2013 TABLE OF CONTENTS Introduction 2 Silencer Performance 4 Flow Resistance and Resistivity 7 Flow Velocity 7 Baffle Attenuation Example 7 Silencer
More informationApplication of the Lattice-Boltzmann method in flow acoustics
Application of the Lattice-Boltzmann method in flow acoustics Andreas Wilde Fraunhofer Institut für Integrierte Schaltungen, Außenstelle EAS Zeunerstr. 38, 01069 Dresden Abstract The Lattice-Boltzmann
More informationDoppler echocardiography & Magnetic Resonance Imaging. Doppler echocardiography. History: - Langevin developed sonar.
1 Doppler echocardiography & Magnetic Resonance Imaging History: - Langevin developed sonar. - 1940s development of pulse-echo. - 1950s development of mode A and B. - 1957 development of continuous wave
More informationFDTD analysis on the sound insulation performance of wall system with narrow gaps
FDTD analysis on the sound insulation performance of wall system with narrow gaps Takumi Asakura a Shinichi Sakamoto b Institute of Industrial Science, The University of Tokyo. Komaba 4-6-, Meguro-ku,
More informationDiagonalization by a unitary similarity transformation
Physics 116A Winter 2011 Diagonalization by a unitary similarity transformation In these notes, we will always assume that the vector space V is a complex n-dimensional space 1 Introduction A semi-simple
More informationSection 7.4: Inverse Laplace Transform
Section 74: Inverse Laplace Transform A natural question to ask about any function is whether it has an inverse function We now ask this question about the Laplace transform: given a function F (s), will
More informationMeasurement of sound absorption coefficient for Fraster felt SpaceCover
DELTA Test Report TEST Reg. no. 100 Measurement of sound absorption coefficient for Fraster felt SpaceCover Performed for Fraster ApS Project no.: I100645 Page 1 of 9 27 November 2014 DELTA Venlighedsvej
More informationWell-posedness, Stability and Conservation for a Discontinuous Interface Problem
Department of Mathematics Well-posedness, Stability and Conservation for a Discontinuous Interface Problem Cristina La Cognata and Jan Nordström LiTH-MAT-R--214/16--SE Department of Mathematics Linköping
More informationAPPENDIX B. Noise Primer
APPENDIX B Noise Primer NOISE PRIMER TABLE OF CONTENTS 1. INTRODUCTION...1 2. BASIC SOUND PRINCIPLES...1 2.1. SOUND AS A WAVE...1 2.2. SOUND PRESSURE LEVEL...2 2.2.1. Decibel Math...4 2.2.2. Descriptive
More informationP2.2 REDUCING THE IMPACT OF NOISE ABATEMENT PRACTICES ON AIRPORT CAPACITY BY FORECASTING SITUATIONAL DEPENDENT AIRCRAFT NOISE PROPAGATION
P2.2 REDUCING THE IMPACT OF NOISE ABATEMENT PRACTICES ON AIRPORT CAPACITY BY FORECASTING SITUATIONAL DEPENDENT AIRCRAFT NOISE PROPAGATION R. Sharman* and T. Keller Research Applications Program National
More informationFinish up undetermined coefficients (on the midterm, on WW3) Physics applications - mass springs (not on midterm, on WW4)
Today Midterm 1 - Jan 31 (one week away!) Finish up undetermined coefficients (on the midterm, on WW3) Physics applications - mass springs (not on midterm, on WW4) Undamped, over/under/critically damped
More informationGMM-based inference in the AR(1) panel data model for parameter values where local identi cation fails
GMM-based inference in the AR() panel data model for parameter values where local identi cation fails Edith Madsen entre for Applied Microeconometrics (AM) Department of Economics, University of openhagen,
More informationIncluding atmospheric propagation effects in aircraft take-off noise modeling
Including atmospheric propagation effects in aircraft take-off noise modeling Michael Arntzen 1,2 ; Martijn Hordijk 2 ; Dick G. Simons 2 1 National Aerospace Laboratory (NLR), Netherlands 2 TU Delft, Aerospace
More informationBeam Propagation Method Solution to the Seminar Tasks
Beam Propagation Method Solution to the Seminar Tasks Matthias Zilk The task was to implement a 1D beam propagation method (BPM) that solves the equation z v(xz) = i 2 [ 2k x 2 + (x) k 2 ik2 v(x, z) =
More informationChristmas Calculated Colouring - C1
Christmas Calculated Colouring - C Tom Bennison December 20, 205 Introduction Each question identifies a region or regions on the picture Work out the answer and use the key to work out which colour to
More informationTheory and Practice of Rotor Dynamics Prof. Dr. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati
Theory and Practice of Rotor Dynamics Prof. Dr. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 2 Simpul Rotors Lecture - 2 Jeffcott Rotor Model In the
More informationStandard Practices for Air Speed Calibration Testing
Standard Practices for Air Speed Calibration Testing Rachael V. Coquilla Bryza Wind Lab, Fairfield, California Air speed calibration is a test process where the output from a wind measuring instrument
More informationUnforced Mechanical Vibrations
Unforced Mechanical Vibrations Today we begin to consider applications of second order ordinary differential equations. 1. Spring-Mass Systems 2. Unforced Systems: Damped Motion 1 Spring-Mass Systems We
More information1 Some general theory for 2nd order linear nonhomogeneous
Math 175 Honors ODE I Spring, 013 Notes 5 1 Some general theory for nd order linear nonhomogeneous equations 1.1 General form of the solution Suppose that p; q; and g are continuous on an interval I; and
More informationOn predicting wind turbine noise and amplitude modulation using Amiet s theory
On predicting wind turbine noise and amplitude modulation using Amiet s theory Samuel Sinayoko 1 Jeremy Hurault 2 1 University of Southampton, ISVR, UK 2 Vestas, Isle of Wight, UK Wind Turbine Noise 2015,
More informationBoundary Conditions for a Divergence Free Velocity-Pressure Formulation of the Incompressible Navier-Stokes Equations
Boundary Conditions for a Divergence Free Velocity-Pressure Formulation of the Incompressible Navier-Stokes Equations Jan Nordström, Ken Mattsson and Charles Swanson May 5, 6 Abstract New sets of boundary
More informationLecture 12. Functional form
Lecture 12. Functional form Multiple linear regression model β1 + β2 2 + L+ β K K + u Interpretation of regression coefficient k Change in if k is changed by 1 unit and the other variables are held constant.
More informationSound Attenuation by Hearing Aid Earmold Tubing
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Sound Attenuation by Hearing Aid Earmold Tubing Mads J. Herring Jensen Widex A/S Ny Vestergaardsvej 19, DK-3500 Vaerloese, Denmark, mjj@widex.dk
More informationPRODUCT DATA. Sound Intensity Calibrator Type 3541-A. Uses and Features
PRODUCT DATA Sound Intensity Calibrator Type 3541-A Type 3541-A enables calibration of sound intensity measuring instruments by using a coupler designed especially for sound intensity calibrations. Intensity-probe
More informationIDR(s) A family of simple and fast algorithms for solving large nonsymmetric systems of linear equations
IDR(s) A family of simple and fast algorithms for solving large nonsymmetric systems of linear equations Harrachov 2007 Peter Sonneveld en Martin van Gijzen August 24, 2007 1 Delft University of Technology
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 informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Architectural Acoustics Session 1pAAa: Advanced Analysis of Room Acoustics:
More informationFan Noise Control by Enclosure Modification
Fan Noise Control by Enclosure Modification Moohyung Lee a, J. Stuart Bolton b, Taewook Yoo c, Hiroto Ido d, Kenichi Seki e a,b,c Ray W. Herrick Laboratories, Purdue University 14 South Intramural Drive,
More information