Femtet Help/Manual
 

Home / Examples / Acoustic Analysis [Mach] / Example 18: Scattering from Rigid Sphere

Example 18: Scattering from Rigid Sphere

 

General

 

  • Apply plane waves.

  • The reflected wave from the rigid sphere is compared with the theoretical wave.
     

  • Obtain this session's project file. (Right-click and choose 'Save link as')

Analysis Space

Item

Settings

Analysis Space

3D

Model Unit

m

 

 

 

Analysis Conditions

Item

Settings

Solver

Acoustic Analysis

Analysis Type

Harmonic Analysis

Options

Incident Wave

 

 

Tab

Item

Settings

Harmonic Analysis

Frequency

Single Frequency

200 [Hz]

Frequency Sweep

Single Frequency

Options

Select

[Enable each port's individual weight setting for the field display].

 

 

Incident Wave (plane wave)

Propagation Direction

(0, 0, 1)

Pressure (Sound Pressure)

1 [Pa]

Select [Sound Pressure [Pa]].

Model

A hole is provided in the body of air. The hole does not need any setting. It is treated as a rigid wall, which has been set for the outer boundary condition.

Body Attribute and Material Property Setting

Body Number/Type

Body Attribute Name

Material Name

Body 4/Solid body

Body_Attribute_001

000_Air *

* Available from the material DB

 

Boundary Conditions

Boundary Condition Name/Topology

Tab

Boundary Condition Type

Settings

Open/Face

Acoustic

Open Boundary

 

Outer Boundary Condition

Acoustic

Rigid Wall

 

 

 

Results

How to separate incident and reflected waves, shown at the top of this page, is explained below.

Operations follow as shown in Example 17. Fig. 1 shows the sound pressure distribution [Pa] at a phase of 0 deg.

 

 

 

 

 

 

Fig. 1 Result Field Displays (sound pressure at phase of 0 deg)

 

 

 

Fig. 2 Result Field Graph (sound pressure at phase of 0 deg)

 

Using the contour diagram in Fig. 1, Fig. 2 is plotted with Θ as the horizontal axis.

The coordinate system is displayed on the right side of Fig. 2. The graph with Θ as its horizontal axis is created under the condition that the distance from the center of the sphere is 3 m.

(To draw the graph, select [On Arc] for the location for evaluation in the Graph Setting dialog box.) The Graph Setting dialog box is explainedhere.

 

 

Next, we compare the results with the theoretical values. The synthetic wave, resulting from the incident and reflected waves, is expressed by the theoretical formula below. [1]

 

 

r,Θ: Observation Point

k: Wavenumber

a: Radius of rigid sphere

P0: Amplitude of incident wave

jn(x): Spherical Bessel function

hn(x): Hankel function of the first kind

Pn(x): Legendre polynominal

 

 

 

 

Fig. 3 Comparison of absolute sound pressure results with theoretical values

 

Fig. 3 is created using Excel.