LONG RANGE PROPAGATION OF FINITE AMPLITUDE ACOUSTIC WAVES IN AN OCEAN WAVEGUIDE Kaelig Castor, Peter Gerstoft, Philippe Roux, W.A. Kuperman,, Marine Physical Laboratory, Scripps Institution of Oceanography B. Edward McDonald,, Naval Research Laboratory, Acoustics Division
0.5 dB/λ, 1800 kg/m3, 1550 m/s
Linear
0.6
200
-20
250
-25
300 0
-30
0.5
50
1 1.5 2 Reduced time (s)
2.5
Nonlinear
0
dB 0
Non-linear
-5
100
-10
150
-15
200
-20
250
-25
300 0
-30
0.5
1 1.5 2 Reduced time (s)
2.5
1500 1495 1490
1
0.2 0
1485
1470
10
20
30 40 50 Frequency (Hz)
60
0.2
10
20
30 40 50 Frequency (Hz)
60
0.6
0.4
0.2
0
70
10
Deep water, r=0km
70
20
30 40 50 Frequency (Hz)
60
r=100km, source near bottom interface (zs=4.5km)
70
1
1
2
Non-linear 4
0.6
1465 1500 0.4 1510
1
1520 1530 phase speed (m/s)
1540
0.2
10
20
30 40 50 Frequency (Hz)
60
0.8
Non-linear Deep water
0.6
0.4
0.2
Nonlinear coupling sea bottom
0.8
0.6
Parametric low frequency
0.4
0.2
70
0
10
20
30 40 50 Frequency (Hz)
60
0
70
10
20
30 40 50 Frequency (Hz)
60
70
Frequency-mode coupling
Time series, deep water ocean waveguide Linear, range=100km
WK31 Dec 99
McDonald & Kuperman, J. Acoust. Soc. Am. 81, 1406-1417, 1987.
∫
+∞
∇ ⊥2 p dx +
δ 2
∂ 2x p + 0.02 α ( dB / λ )
c0 2
∫
+∞ 0
p ( x ') dx' x − x'
Deep water time series (10 Hz) 50 km Nonlinear Linear
Bottom attenuation
Diffraction : Crank-Nicholson
Linear, range=3000km
Integration region (moving spatial frame) : x
min
< x < x
- more uniform modal distribution - self-refraction
−5 −10
−0.5
0
20
40
60
80
100
120
140
160
180
200
Time (s)
The nonlinear NPE code and a time domain finite difference code (Cabrillo) and are used to generate movies. NPE shallow water simulation NPE deep water simulation
MOVIES! Depth (km)
Time series, shallow to Deep water
x 10
2000 km
-8
reduced time (s)
2
-8
Shallow to deep water (10 Hz) Linear, range=100km
-4
1
2
0.5
-100
0 -200
-0.5 -1
200 m 4 km
Nonlinear, range=100km
-300 0
480 Range (m)
960
In Depth FDTD modeling: Using Cabrillo Linear
•Staggered Fourier pseudo spectral method •Fourier spectral methods requires less grid points than classical FD (l/2 vs l/10) •Can model both acoustic, elastic and poroelastic (Biot) media. •FD method can better model variations in sound speed (including bathymetry) than classical ocean acoustic propagation codes. Any grid point can have different properties! •The wrap around is destroyed by tapering of the grid (last 30). Shallow to deep water Shallow water to T-phase station.
50 km Nonlinear
500 km Linear, range=3000km
Nonlinear, range=3000km Depth (km)
Nonlinear steepening Multivalued waveform no physical sense Shock dissipation Shock wave formation
0
0
- harmonic generation (2f, 3f, 4f,...)
- shock dissipation
5
−20
max
- parametric interaction (f1 ± f2)
10
30
−15
Broader frequency spectrum Structure changes during propagation
15
40
0
500 km
Nonlinear effects Additional frequencies
20
50
0.5
1000 km
Important for an ocean waveguide: - Refraction included. - Porous medium attenuation in the bottom
25
10
moving frame time incremental step ∆ t = ∆ x/c0
c0
30
70
20
Nonlinear, range=3000km
Thermoviscous dissipation : Finite difference scheme Refraction + Nonlinear steepening : Second order upwind flux corrected transport scheme [B. E. McDonald, J. Comp. Phys. 56, 448-460, (1984)] NPE moving frame
WK30 Dec 99
80
60
Frequency (Hz)
For both shallow and deep water the NPE is propagating the field the first 20 km where nonlinearities are strong. An adiabatic normal mode code is used for propagating the field to longer ranges. 4 km
Nonlinear Progressive wave Equation (NPE) x
Increased sea-bottom coupling (nonlinearities can be candidate for T-wave formation)
But, we also work with real data!
Non-linear, range=100km
Long range propagation
⎡ ⎤ c β ∂ t p = −∂ x ⎢c1 p + p2 ⎥ − 0 2 ρ0c0 ⎦ 2 ⎣
Nonlinear
Linear
Deep water, r=100km
Nonlinear, r=100km
0.8
0
Non-linear Shallow water
0.4
0
2
3
1480 1475
Shallow water fondamental frequency first harmonic
1 0.4
0.6
Depth-averaged spectrum
0.8
Depth-averaged spectrum
-15
Depth-averaged spectrum
150
Depth (m )
-10
Linear, r=100km
1
-5
100
Depth-averaged spectrum
0
Linear
1
0.8
Depth (m)
Nonlinearities can give effects similar to long-range propagation in random media: → modal distribution + arrival time structure → increased sea-bottom coupling (NL candidate for T-phase generation)
dB
Linear
0
1
0.8
Nonlinear
Linear
Shallow water, r=100km
2nd
1st mode
1000 km
0 2nd
1st mode
2000 km -1 reduced time (s) 1 -1 1 Only few modes excited: little time spread, energy close to sound speed axis. The Nonlinear propagation excite higher order modes
Depth (km)
Frequency-Mode Coupling
- frequency distribution → parametric and harmonic generation - modal distribution → larger number of modes excited → tendency toward modal equipartition → changes in the arrival time structure
Parametric low frequency might not appear in a shallow water waveguide
Depth-averaged spectrum
Range=100km
r=100km, source on sofar axis (zs=800m)
Normalized depth-averaged spectrum Shallow water, r=0km
Depth (m )
What is the main difference between shallow and deep water ? In shallow water, => lower number of modes (bottom interaction), faster time-separation of modes. => lower geometrical spreading, higher amplitudes, stronger nonlinear effects How can we identify at long ranges a nonlinear acoustic propagation path ? Redistribution of the energy during the propagation
Source spectrum
Source 50Hz narrowband, depth 100m Source over density Rm=3.5 10-3.
50
Conclusion
Deep water ocean waveguide: influence of parametric mode conversion on sea-bottom coupling
Source time series 1500 m/s
group speed (m /s)
Abstract The Nonlinear Progressive Wave Equation (NPE) [McDonald and Kuperman, 1987] computer code was coupled with a linear normal mode code in order to study propagation from a high intensity source in either shallow or deep water. Simulations using the coupled NPE/linear code are used to study both harmonic (high frequency) and parametric (low frequency) generation and propagation in shallow or deep water with long-range propagation paths. Included in the modeling are both shock dissipation and linear attenuation in the bottom.
Shallow water Pekeris waveguide 200 m
Depth-averaged spectrum
Background: The nonlinear progressive wave equation (NPE) [McDonald and Kuperman, JASA, 1987] was developed to obtain accurate and affordable simulations of shock propagation in the deep ocean out to convergence zone ranges.
Sponsored by Defense Threat Reduction Agency Contract No. DTRA01-00-C-0084
0.8 0
Range (km)
1.5
•Cair=0 m/s, Cwater=1500 m/s, Cbot=2500 m/s •Source at 70 m depth •Field tapered outside red box
0
Range (km)
1.5