OS-V: 1110 Vibrating Sphere: Exterior Acoustic Analysis using Infinite Elements (IE) and Adaptive Perfectly Matched Layer (APML) Methods

Acoustic modeling in finite and semi-infinite domains is essential in the prediction of quantities such as external and radiated noise in vibro-acoustic problems.

APML is a popular way of modeling these domains. If sound pressure at microphone locations is to be calculated because of sound propagating through sections of the fluid domain and through panels, this method shows the fidelity of various vibrating sound sources, such as speakers, as it allows prediction of radiated noise.

Model Files

Before you begin, copy the file(s) used in this problem to your working directory.

Benchmark Model

Figure 1.

The sphere is of 1 m radius.

All the nodes of the model are constrained to six degrees of freedom (123456), along with an enforced velocity of a 1.0 m/s amplitude on SPCD via RLOAD1 in DOF 3.

The loading frequencies at which the responses are calculated are specified using the FREQ1 entry starting from 50 Hz to 160 Hz in increments of 10 Hz.

For APML, the entire vibrating structure is enclosed with an acoustic cavity mesh and further adding a layer of PML elements on this enclosed acoustic cavity mesh. A minimum of one layer of fluid elements will be defined on the surface of the structural domain of interest. Then, the APML elements, CACPML3 and CACPML4 will only be defined on the topmost surface of the fluid elements (Figure 2).

Similarly for IE, the vibrating structure is enclosed with an acoustic cavity mesh and a layer of Infinite Elements (CACINF3 and CACINF4) on the enclosed acoustic cavity mesh is added.

Figure 2.

Units: m, s, Pa, kg/m3

Material

Sphere shell is aluminum which is specified using MAT1 Bulk Data Entry. Fluid material properties (bulk modulus, speed of sound, fluid density) are specified for the fluid cavity elements on the MAT10 Bulk Data Entry. For this model, the fluid is assumed to be air.

Sound pressure can be measured on the receiver grid points (microphone locations). These microphones are located at coordinates z = -4 m and z = 10 m as fluid grids 9003 and 9000, respectively in the frequency range 50 to 160 Hz.

Analytical Calculation

The spatially dependent amplitudes of the field quantities are identified as:

Φ= v ^ C a 3 e ika cosθ 2 ka 2 2ika d dr e ikr r MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuOPdyKaey ypa0ZaaSaaaeaaceWG2bGbaKaadaWgaaWcbaGaam4qaaqabaGccaWG HbWaaWbaaSqabeaacaaIZaaaaOGaamyzamaaCaaaleqabaGaeyOeI0 IaamyAaiaadUgacaWGHbaaaOGaci4yaiaac+gacaGGZbGaeqiUdeha baGaaGOmaiabgkHiTmaabmaabaGaam4AaiaadggaaiaawIcacaGLPa aadaahaaWcbeqaaiaaikdaaaGccqGHsislcaaIYaGaamyAaiaadUga caWGHbaaamaalaaabaGaamizaaqaaiaadsgacaWGYbaaamaabmaaba WaaSaaaeaacaWGLbWaaWbaaSqabeaacaWGPbGaam4Aaiaadkhaaaaa keaacaWGYbaaaaGaayjkaiaawMcaaaaa@5989@

Sound pressure is defined as:

p=iωρ*Φ MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCaiabg2 da9iaadMgacqaHjpWDcqaHbpGCcaGGQaGaeuOPdyeaaa@3E92@

Where,
ρ
Density of medium
c MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaaaa@36DB@
Speed of sound in the medium
k MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaaaa@36DB@
Wave number (=circular loading frequency/speed of sound)
a MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyaaaa@36DA@
Radius of the sphere
v ^ C MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmODayaaja WaaSbaaSqaaiaadoeaaeqaaaaa@37F3@
Prescribed oscillatory velocity
r MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaaaa@36DB@
Radial distance of the microphone location
θ
Polar angle of microphone location (angle measured about the +ve z axis)

Results

Sound Pressure versus Frequency for two microphone locations is plotted.

Adaptive Perfectly Matched Layer (APML), Infinite Elements (IE) and Analytical results are compared.

The goal of this project is to understand how OptiStruct methods, APML and IE correlate with the analytical method and to check the accuracy of the OptiStruct results with the Analytical results.
Table 1. Average Percentage Difference of Sound Pressure between APML and IE with Analytical
APML IE
Z = 10 8.55% 24.46%
Z = -4 7.79% 11.15%
From the plots below, APML and IE correlate with the analytical method.
  • Microphone Pressure versus Frequency at z = 10 m
    Figure 3. Real
    Figure 4. Imaginary
  • Microphone Pressure versus Frequency at z = 10 m
    Figure 5. Real
    Figure 6. Imaginary
1 Pierce, Allan D. Acoustics: An Introduction to Its Physical Principles and Applications. Springer, 2019