r p = r o r cos( φ ) cos( α )

Similar documents
Example 2.1. Draw the points with polar coordinates: (i) (3, π) (ii) (2, π/4) (iii) (6, 2π/4) We illustrate all on the following graph:

COMPUTATIONAL ANALYSIS OF ACOUSTIC HORNS USING AN MFS MODEL

1 Conduction Heat Transfer

FastBEM Acoustics. Verification Manual , Advanced CAE Research, LLC (ACR) Cincinnati, Ohio, USA All Rights Reserved

THE RADIATION EFFICIENCY OF FINITE SIZE FLAT PANELS

14.1. Multiple Integration. Iterated Integrals and Area in the Plane. Iterated Integrals. Iterated Integrals. MAC2313 Calculus III - Chapter 14

P773 Acoustics June, 2008

1.060 Engineering Mechanics II Spring Problem Set 3

42. Change of Variables: The Jacobian

SEAFLOOR MAPPING MODELLING UNDERWATER PROPAGATION RAY ACOUSTICS

Sound radiation from a loudspeaker, from a spherical pole cap, and from a piston in an infinite baffle1

Multi Acoustic Prediction Program (MAPP tm ) Recent Results Perrin S. Meyer and John D. Meyer

Mathematical modeling of the radiated acoustic field of ultrasonic transducers

Chapter 15. Mechanical Waves

Open Baffle w/ 2 Drivers - Acoustic Response 12/30/14. Copyright 2014 by Martin J. King. All Rights Reserved.

Copyright 2014 by Martin J. King. All Rights Reserved. Configuration : Extended Range Driver w/ Bass Driver Mounted on an Open Baffle

Sound radiation and transmission. Professor Phil Joseph. Departamento de Engenharia Mecânica

Introduction to Acoustics Exercises

Open Baffle w/ 2 Drivers - Acoustic Response 12/27/12. Copyright 2012 by Martin J. King. All Rights Reserved.

Sound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur

Physics 208: Electricity and Magnetism Final Exam, Secs May 2003 IMPORTANT. Read these directions carefully:

MP204 Electricity and Magnetism

Section 7.0 : Design of a Back Loaded Exponential Horn By Martin J. King, 08/17/04 Copyright 2004 by Martin J. King. All Rights Reserved.

CBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

Chapter 5 Oscillatory Motion

LECTURE 18: Horn Antennas (Rectangular horn antennas. Circular apertures.)

Lecture 17 - The Secrets we have Swept Under the Rug

AP Calculus AB Chapter 2 Test Review #1

Lecture notes 5: Diffraction

Muffler Transmission Loss Simple Expansion Chamber

A R T A - A P P L I C A T I O N N O T E

THE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))

JUST THE MATHS UNIT NUMBER INTEGRATION APPLICATIONS 13 (Second moments of a volume (A)) A.J.Hobson

Lecture 10 - Moment of Inertia

DIFFRACTION CORRECTION OF FREQUENCY RESPONSE FOR LOUDSPEAKER IN RECTANGULAR BAFFLE A. DOBRUCKI

Selected electrostatic problems (Lecture 14)

Sinusoids. Amplitude and Magnitude. Phase and Period. CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation

MULTIVARIABLE INTEGRATION

LECTURE 18: Horn Antennas (Rectangular horn antennas. Circular apertures.) Equation Section 18

Impedance/Reactance Problems

Chapter 7: Applications of Integration

OPAC102. The Acoustic Wave Equation

CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation

Design of Mufflers and Silencers. D. W. Herrin, Ph.D., P.E. University of Kentucky Department of Mechanical Engineering

EFFECTS OF PERMEABILITY ON SOUND ABSORPTION AND SOUND INSULATION PERFORMANCE OF ACOUSTIC CEILING PANELS

Contents. Objectives Circular Motion Velocity and Acceleration Examples Accelerating Frames Polar Coordinates Recap. Contents

Solutions to Math 41 Second Exam November 5, 2013

The Rotating Mercury Reflector

PHY 6347 Spring 2018 Homework #10, Due Friday, April 6

First Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields

CHAPTER 20 Magnetism

Calculation of Reflection and Transmission Coefficients in scuff-transmission

THE ACOUSTIC POWER RADIATED BY A CIRCULAR MEMBRANE EXCITED FOR VIBRATION BOTH BY MEANS OF THE EDGE AND BY EXTERNAL SURFACE LOAD

Design of Partial Enclosures. D. W. Herrin, Ph.D., P.E. University of Kentucky Department of Mechanical Engineering

ENGI 4430 Parametric Vector Functions Page dt dt dt

Chapter 17 : Fourier Series Page 1 of 12

The choice of origin, axes, and length is completely arbitrary.

University of Illinois at Chicago Department of Physics. Electricity & Magnetism Qualifying Examination

Integration to Compute Volumes, Work. Goals: Volumes by Slicing Volumes by Cylindrical Shells Work

The Scattering of Light by Small Particles. Advanced Laboratory, Physics 407 University of Wisconsin Madison, Wisconsin 53706

WAVE ACOUSTICS. F = k x (6.1) where k is the constant of proportionality or the spring constant.

Chapter 5. (b) One route to solving the above equation is to use the relationship between the observer divergence and the observer time derivative.

PHY 114 Summer Midterm 2 Solutions

cancel each other out. Thus, we only need to consider magnetic field produced by wire carrying current 2.

Lab Reports Due on Monday, 11/24/2014

Chapter 10 Sound in Ducts

Since x + we get x² + 2x = 4, or simplifying it, x² = 4. Therefore, x² + = 4 2 = 2. Ans. (C)

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx

Electro Magnetic Field Dr. Harishankar Ramachandran Department of Electrical Engineering Indian Institute of Technology Madras

Review of Fundamental Equations Supplementary notes on Section 1.2 and 1.3

θ + mgl θ = 0 or θ + ω 2 θ = 0 (2) ω 2 = I θ = mgl sinθ (1) + Ml 2 I = I CM mgl Kater s Pendulum The Compound Pendulum

Kinematics. Vector solutions. Vectors

Unit #13 - Integration to Find Areas and Volumes, Volumes of Revolution

Math Exam IV - Fall 2011

Problems (F/M): Part 2 - Solutions (16 pages; 29/4/17)

Three Loose Ends: Edge Focusing; Chromaticity; Beam Rigidity.

Continuum Limit and Fourier Series

Problem Set 7: Solutions

Light as a Transverse Wave.

Acoustic radiation by means of an acoustic dynamic stiffness matrix in spherical coordinates

6.17 The Lossy Voice-Coil Inductance

Computing the Acoustic Field of a Radiating Cavity by the Boundary Element - Rayleigh Integral Method (BERIM)

Moment of a force (scalar, vector ) Cross product Principle of Moments Couples Force and Couple Systems Simple Distributed Loading

Q ( q(m, t 0 ) n) S t.

Chapter 27 Magnetism 1/20/ Magnets and Magnetic Fields Magnets and Magnetic Fields Magnets and Magnetic Fields

Transmission Loss Assessment for a Muffler by Boundary Element Method Approach

IIT JEE PAPER The question paper consists of 3 parts (Chemistry, Mathematics and Physics). Each part has 4 sections.

EXERCISES. Let us see an example of graphical methods for solving static analysis associated with numerical evaluation of the results. 2 1 E 3 P=?

Moments of Inertia (7 pages; 23/3/18)

Part IB GEOMETRY (Lent 2016): Example Sheet 1

ERRATA AND ADDITIONS FOR "ENGINEERING NOISE CONTROL" 4th Edn. First printing April 23, 2018

Wave Phenomena Physics 15c. Lecture 17 EM Waves in Matter

General Physics I. Lecture 8: Rotation of a Rigid Object About a Fixed Axis. Prof. WAN, Xin ( 万歆 )

Some negative results on the use of Helmholtz integral equations for rough-surface scattering

Wave Phenomena Physics 15c

Wave Phenomena Physics 15c. Lecture 8 LC Transmission Line Wave Reflection

(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves

BEHAVIOR OF THE REGULARIZED SAMPLING INVERSE SCATTERING METHOD AT INTERNAL RESONANCE FREQUENCIES

10550 PRACTICE FINAL EXAM SOLUTIONS. x 2 4. x 2 x 2 5x +6 = lim x +2. x 2 x 3 = 4 1 = 4.

Transcription:

Section 4. : Sound Radiation Pattern from the Mouth of a Horn In the previous section, the acoustic impedance at the mouth of a horn was calculated. Distributed simple sources were used to model the mouth geometry moving as a piston in an infinite baffle. The output of all the simple sources was integrated to calculate the total acoustic impedance of the horn s mouth. One additional property of the mouth of the horn needs to be derived before the complete acoustic response of a horn can be simulated. At any position in the listening environment, the distances from any two individual simple sources being used to represent the mouth will be slightly different. As the frequency increases, these small differences in path length start to represent larger phase differences leading to cancellation and reinforcement of the summed SPL at the listening position. Expressions for the SPL, as functions of distance and angular location with respect to the horn s axis, will be derived for circular and rectangular mouth geometries. Circular Horn Mouth : To derive the SPL radiation pattern of a circular horn mouth, the piston in an infinite baffle analogy is assumed again. Figure 4.1 shows the geometry definition used to derive the far field SPL at an arbitrary position in space for a vibrating circular piston in an infinite baffle. The first step in calculating the radiation pattern is to express r p in terms of. At distances in the far field, the relationship between these vectors can be stated as follows r p r cos( φ cos( α where φ and α are the angles that r and form with the X axis. Since the response of the circular piston is axisymmetric, only the solution in the X-Z plane is required to describe the complete 3D radiation pattern. For the simple source at the tip of r, the pressure generated at a point in the far field X-Z plane is shown below. dp ( r p I ( I ( ω t k r p 2 ρ fu rms dr r d φ e r p Substituting the expression for r p and formulating an integral to sum all of the point sources on the piston s surface yields the following expression. p ( I ( Ik ( I ω t 2 ρ fu rms e e a 2 π e ( Ikr cos( φ cos( α r dφ dr Page 1 of 9

Figure 4.1 : Circular Piston Geometry Definition Z ro rp zo r Y (xo,yo X Evaluating the integral is fairly straight forward requiring trigonometry and careful substitutions. After evaluating the double integral, and making all of the appropriate substitutions, the following well know expression results. p ( r, θ, t I 2 ρ f π a 2 u rms 2 J 1 ( ka sin( θ ka sin( θ ( Ik ( I ω t e e The angle θ is measured between and the Z axis. Page 2 of 9

Rectangular Horn Mouth : Section 4. : Sound Radiation Pattern from the Mouth of a Horn To derive the SPL radiation pattern of a rectangular horn mouth, the piston in an infinite baffle analogy is assumed again. Figure 4.2 shows the geometry definition used to derive the far field SPL at an arbitrary position in space for a vibrating rectangular piston in an infinite baffle. The first step in calculating the radiation pattern is to express r p in terms of. At distances in the far field, the relationship between these vectors can be stated as follows r p ( x cos( α + y cos( β where α and β are the angles that forms with the X axis and Y axis respectively. Since the response of the rectangular piston is symmetric in the X-Z and X-Y planes only the patterns in these two planes are required to describe the complete 3D radiation pattern. For the simple source at the tip of r, the pressure generated at a point in the far field is shown below. dp ( r p I ( I ( ω t k r p 2 ρ fu rms dx dy e r p Since 1/ r p ~ 1/ in the far field. dp ( I 2 ρ fu rms e ( I ω t ( Ik ( Ikx cos( α ( Iky cos( β e e e dx dy Making two more substitutions k x k cos( α and k x k sin( θ x k y k cos( β k y k sin( θ y the differential pressure can be simplified to the following expression. The angles θ x and θ y are measured from the Z axis. Page 3 of 9

Figure 4.2 : Rectangular Piston Geometry Definition Z ro rp zo Y r (xo,yo X dp ( I 2 ρ fu rms e ( I ω t ( Ik ( Ik x x ( Ik y y e e e dx dy Summing the pressure produced by all of the simple sources, making up the rectangular surface, leads to the following integral. p (, θ x, θ y I 2 ρ fu rms e ( I ω t ( Ik e a a ( Ik x x e b dx b ( Ik y y e dy Page 4 of 9

This integral is really very straight forward to evaluate and the following result is obtained. p (, θ x, θ y I 2 ρ fu rms e ( I ω t ( Ik sin( k x a e 2 a 2 b k x a sin( k y b k y b The pressure distribution forms two different patterns in the X-Z and Y-Z planes as determined by the aspect ratio between the rectangle s dimensions. Summary of Sample Radiation Patterns : To get a feel for the radiation pattern of the sound pressure emanated by the mouth of a horn, the sample problem from the previous section was used again. Assuming a constant mouth area equal to a 32 inch by 32 inch square, the equivalent radius a L for a circular mouth of equal area was calculated. The radius of this equivalent area circular piston is 18.54 inches. The following plots show the radiation patterns of a simple source, a circular horn mouth, a square horn mouth, and a rectangular horn mouth all calculated in the far field and at a frequency of 2 Hz. The results have all been normalized to the value produced by the simple source. Figure 4.3 presents the results for the circular mouth and the square mouth. Notice that at this relatively low frequency the two different shaped mouths produce equivalent pressure radiation patterns. Also notice that compared to the simple source the sound field is starting to become directional. Figures 4.4 and 4.5 present the results for the circular mouth and two different rectangular mouths. The aspect ratios, width divided by height, are.2 and.1 respectively. In the horizontal plane, the narrow width allows the rectangular mouth s radiation pattern to approach the simple source result while the circular mouth is becoming directional. In the vertical plane, the two rectangular mouths have become very directional. At low frequencies, below 1 Hz, these horn mouths are non directional and can be modeled as a simple source. Above 1 Hz, the shape of the horn mouth starts to generate a radiation pattern caused by the cancellation of pressure due to differing path lengths (really phase differences between the distributed simple sources used to model the mouth geometry. Some final comments on the data presented in Figures 4.3, 4.4, and 4.5 are worth restating. The on-axis response of any horn mouth geometry will be equal to the equivalent simple source response. The off-axis response can become a complicated 3D radiation pattern generated by a distribution of simple sources used to represent the mouth geometry. Using a low aspect ratio horn mouth, the radiation pattern in one plane can be made to approach the simple source result while in the perpendicular plane it becomes highly directional at a very low frequency. Page 5 of 9

In this section, and the proceeding one, there are some interesting acoustic results presented related to the horn mouth which give rise to a few options to ponder when designing a horn. The horn mouth is going to be a key variable in the design of successful and consistent horn geometries. Understanding the acoustic impedance and the sound pressure radiation pattern, for a given mouth geometry, is critical. Page 6 of 9

Figure 4.3 : Pressure Radiation Patterns From a Simple Source, a Circular Horn Mouth, and a Square Horn Mouth (Aspect Ratio width / height 1. Calculated in the Far Field at a Frequency of 2 Hz. black curve simple source red curve circular horn mouth blue curve square horn mouth (horizontal plane green curve square horn mouth (vertical plane 9 12 1 6 15.5 3 p p ( r, θ, ω ( p r r, θ x, deg, ω 18 ( p r r, deg, θ y, ω 21 33 24 3 27 θ, θ, θ x, θ y Page 7 of 9

Figure 4.4 : Pressure Radiation Patterns From a Simple Source, a Circular Horn Mouth, and a Rectangular Horn Mouth (Aspect Ratio width / height.2 Calculated in the Far Field at a Frequency of 2 Hz. black curve simple source red curve circular horn mouth blue curve square horn mouth (horizontal plane green curve square horn mouth (vertical plane 9 12 1 6 15.5 3 p p ( r, θ, ω ( p r r, θ x, deg, ω 18 ( p r r, deg, θ y, ω 21 33 24 3 27 θ, θ, θ x, θ y Page 8 of 9

Figure 4.5 : Pressure Radiation Patterns From a Simple Source, a Circular Horn Mouth, and a Rectangular Horn Mouth (Aspect Ratio width / height.1 Calculated in the Far Field at a Frequency of 2 Hz. black curve simple source red curve circular horn mouth blue curve square horn mouth (horizontal plane green curve square horn mouth (vertical plane 9 12 1 6 15.5 3 p p ( r, θ, ω ( p r r, θ x, deg, ω 18 ( p r r, deg, θ y, ω 21 33 24 3 27 θ, θ, θ x, θ y Page 9 of 9