TS.sphere2: Calculate the analytical solution for a weakly scattering...

View source: R/TSsphere.R

TS.sphere2R Documentation

Calculate the analytical solution for a weakly scattering sphere

Description

Scattered pressure from an incident plane wave upon a fluid sphere. Based on: Jech, J. M., Horne, J. K., Chu, D., Demer, D. A., Francis, D. T., Gorska, N., ... & Reeder, D. B. (2015). "Comparisons among ten models of acoustic backscattering used in aquatic ecosystem research." The Journal of the Acoustical Society of America, 138(6), 3742-3764. modified from: Anderson, V. C. (1950). "Sound scattering from a fluid sphere." The Journal of the Acoustical Society of America, 22(4), 426-431.

Usage

TS.sphere2(f, r, a, c, h, g, rho)

Arguments

f

Frequency in Hz (s^-1)

r

Range in m from center of sphere

a

Radius of sphere

c

Soundspeed in surrounding fluid m/s

h

Soundspeed contrast inside/surrounding fluid

g

Density contrast inside/surrounding fluid

rho

Density of surrounding fluid

Examples

fs <- as.list(seq(10,400, by=1)*1000) #Frequencies
r <- 10 #range
a <- 0.01 # radius
c <- 1477.4 #soundspeed surrounding fluid
rho <- 1026.8 #density surrounding fluid
g <- 1028.9/rho #density contrast
h <- 1480.3/c #soundspeed contrast
TS <- sapply(fs,TS.sphere2,r=r,a=a,c=c,h=h,g=g,rho=rho)
plot(fs,TS, type="l", xlab="Frequency [Hz]",ylab="TS [dB re m2]")

AustralianAntarcticDivision/ZooScatR documentation built on Aug. 13, 2022, 1:21 a.m.