Ultrasound propagation in shallow pond

Jupiter notebook specific imports

import os
os.chdir('../../../')
import warnings
warnings.filterwarnings('ignore')

PyWaveProp imports

from uwa.source import GaussSource
from uwa.environment import UnderwaterEnvironment, Bathymetry, sound_speed_mps
from uwa.sspade import UWASSpadeComputationalParams, uwa_ss_pade
from uwa.vis import AcousticPressureFieldVisualiser2d
from scipy.interpolate import interp1d

Preparing environment

max_range_m = 100
env = UnderwaterEnvironment(
    sound_speed_profile_m_s=lambda x, z: 1500,
    bottom_profile=Bathymetry(
        ranges_m=[0, max_range_m],
        depths_m=[4, 1]),
    bottom_sound_speed_m_s = 1700,
    bottom_density_g_cm = 1.5,
    bottom_attenuation_dm_lambda = 0.1
)

Horizontally oriented antenna pattern

Preparing transmitting antenna

src = GaussSource(
    freq_hz=31.4E3,
    depth_m=2,
    beam_width_deg=20,
    elevation_angle_deg=0,
    multiplier=0.1
)

Calculating the acoustics pressure field

params = UWASSpadeComputationalParams(
    max_range_m=max_range_m,
    max_depth_m=6,
    dx_m=0.1, # output grid steps affects only on the resulting field, NOT the computational grid
    dz_m=0.1,
    precision=0.1
)

field1 = uwa_ss_pade(
    src=src,
    env=env,
    params=params
)

Visualising the results

Two dimensional distribution of the field amplitude

vis1 = AcousticPressureFieldVisualiser2d(field=field1, env=env, label="Horizontal")

vis1.plot2d(min_val=-60, max_val=0, grid=True, show_terrain=True).show()

60 deg. elevation angle

src.elevation_angle_deg = -60

field2 = uwa_ss_pade(
    src=src,
    env=env,
    params=params
)

vis2 = AcousticPressureFieldVisualiser2d(field=field2, env=env, label="Vertical")
vis2.plot2d(min_val=-60, max_val=0, grid=True, show_terrain=True).show()

vis1.plot_hor(0.5, [-70, -10], vis2).show()

vis1.plot_hor(2, [-70, -10], vis2).show()