Research Article  Open Access
WenKai Wang, GuangPing Zhu, JingWei Yin, Hui Sun, "Reflection Characteristics of a NearInterface Cavity in Ice at Supercritical Incidence", Mathematical Problems in Engineering, vol. 2019, Article ID 2143487, 11 pages, 2019. https://doi.org/10.1155/2019/2143487
Reflection Characteristics of a NearInterface Cavity in Ice at Supercritical Incidence
Abstract
The longrange propagation modes in an acoustic channel under ice are basically caused by supercritical incidence. The energy distribution and transmission loss in the acoustic channel under ice are changed by a scatter in ice. The influence of a slender cylindrical cavity near and parallel to the icewater interface on the sound propagation is analyzed using FourierBessel series and SommerfeldWatson transformation. The research found that the acoustic field presents a beam in the mirror reflection direction at supercritical incidence, and the beamwidth is proportional to secant of incident angle; meanwhile, the reflected coefficient is proportional to cosine of incident angle. The reflection coefficient increases with relative depth and Helmholtz number if the incident angle is a constant.
1. Introduction
The Arctic region has important strategic significance and commercial value. The sound propagation characteristics in an acoustic channel under ice depend on the parameters of the ice canopy, which has a noticeable influence on the acoustic navigation and communication under ice [1, 2]. The ice canopy in current researches was considered as an ideal elastic layer on the macro scale, or its reflection coefficient was obtained through experiments to calculate the transmission loss, but the influence of microstructure in ice on the acoustic channel under ice has not been considered yet. The microstructure changes not only the attenuation of sound wave in the original propagation direction but also the spatial distribution of acoustic energy under ice. Therefore, it is of great significance to study the influence of microstructures in ice on the sound propagation characteristics in an acoustic channel under ice.
Acoustic experimental researches in the Arctic region can be traced back to the last century [3]. The ice area, climate, and earthquakes were monitored through the experimental data of underwater acoustic propagation [4, 5]. The ice canopy has an influence on the transmission loss in an acoustic channel under ice, and the transmission loss at different depths might exceed 10dB [6], but the roughness of the ice canopy was not considered. The BELLHOP model was introduced to calculate the transmission loss in shallow water under ice canopy and the performance of a vertical array was analyzed [7], but the inhomogeneity of the ice canopy was not contemplated. The phenomenon of overlapping sound lines [8] occurs if the sound signals propagate over a long distance under ice, and the communication can be realized by applying adaptive coherent communication technology. The transmission loss of shallow source is very large because of the scattering of the icewater mixed layer [9] and the signal can hardly be propagated for a long distance. The longrange propagation needs certain conditions [10], such as the source should be located at a specific depth; however, the constraint effect of sound velocity gradient in icewater mixed layer was not taken into account. The ice layer was considered as a homogeneous medium without microstructure in above researches, and the authors only focused on the macroscopic equivalent parameters of the ice layer, other than the influence of microstructure in ice on the macroscopic characteristics.
The scattering by microstructures in ice is an elastic scattering problem. The scattering of elastic waves is accompanied by the conversion of particle polarization modes and energy redistribution of longitudinal waves and transverse waves. Therefore, the physical characteristics and mathematical methods of the elastic scattering problem are more complicated than those in fluids [11]. Many methods were introduced to solve the problem [12], among which the FourierBessel series [13] approximating the plane halfspace surface as a circular surface of large radius is an effective method. However, it is generally used to solve halfspace problems [14, 15], such as the scattering of seismic waves near the earth's surface[16], and it has poor convergence for the scattering problem in two contiguous halfspaces. The SommerfeldWatson transformation [17, 18] is a more efficient solution to the scattering problem of a regular cylindrical obstacle, but it cannot be directly applied to the scattering problem of a plane interface [19]. In engineering applications, the calculation error can be controlled within an acceptable range by selecting reasonable parameters [20]. The combination of FourierBessel series and SommerfeldWatson transformation gives full play to the advantages of two methods, which can not only realize the fast convergence of the underwater acoustic field but also retain the main characteristics of the spatial distribution of acoustic energy. The influence of a slender cylindrical cavity near and parallel to the icewater interface on the acoustic propagation under ice is analyzed using the combination of FourierBessel series and SommerfeldWatson transformation.
2. Solutions
The ice and water media can be considered as in two contiguous halfspaces if the ice thickness and water depth are much larger than the wave train length of a short pulse source; the problem is simplified into a twodimensional model as shown in Figure 1.
A plane wave in water is incident on the icewater interface and generate reflected wave and transmitted waves and . The acoustic field in ice and water will be distorted due to the cavity as shown in Figure 1. Assume that scattering of the cavity is and , scattering of the icewater interface in ice is and , and scattering of the icewater interface in water is .
A convex approximation of the icewater interface is introduced because it is very complicated to deal with the horizontal interface using cylindrical functions. The radius of the cavity is , and the radius of the icewater interface is . A minor error will be introduced to the calculation by the geometric deviation if . Two cylindrical coordinate systems and are established, respectively, with the centers of the cavity and icewater interface. The displacement potential in ice is set as follows:
The coefficients in the scalar or vector potential functions are represented, respectively, by or , and , , , and are the known coefficients of transmitted waves [14, 15]; in subscripts represents or coordinate system, in superscripts is the integral orders of cylindrical function, or in superscripts corresponds to the parity of trigonometric function, is the longitudinal wave number, and is the transverse wave number.
Acoustic field in fluid is usually expressed by sound pressure. Therefore, the expression of the scattering field in water is set as follows:where represents the wave number in water. The SommerfeldWatson transformation is appropriate to reform the acoustic field into a new complexorder series solution to achieve fast convergence. The transformation is carried out on (2):
As shown in Figure 2, in the integrand is a singular function when is an integer, while is a singular function when is a complex, where is a singular point in complex plane and is an integral path that bypasses the secant including all the complex singular points .
It can be proved that equals the sum of residues at all the singular complex points in complex plane. The complex causes exponential decay of the acoustic field propagating along the circumferential direction. Therefore, SommerfeldWatson transformation can realize fast convergence for the problem in this paper. The displacement continuity equation shown in appendix (A.3) is substituted into (3) and a steepest descent equation is obtained:
A graphic method is used to solve the transcendental equation in (4), as shown in Figure 3.
is the observation direction and is the reflection of direct wave on the icewater interface. Then SommerfeldWatson transformation is carried out on the stress continuous equation shown in appendix (A.1) and (A.2) and another steepest descent equation is obtained:
In the same way, a graphic method is used to solve the transcendental equation of formula (5), as shown in Figure 4.
As can be seen that and coexist in Figure 4; that is, the scattering of the cavity in ice is still expressed by the classical integer series with the order , and the scattering of the icewater interface is jointly represented by the steady phase integral at the point and the complex series order . The acoustic field under ice is obtained by solving the above simultaneous equations.
3. Numerical Analysis
The medium and geometric parameters are shown in Table 1. The incident angle is defined as the included angle between the incident direction of and the normal vector of icewater interface, is the scattering angle, and is the mirror reflection direction.

The sound pressure level of underwater scattering field vs. and for and is shown in Figure 5.
As can be seen from Figure 5, the scattering field presents a fast fluctuating interference pattern due to the joint action of the cavity and icewater interface if and . The acoustic energy in the direction of mirror reflection is much higher than that in other scattering angles, and the reflected wave presents a beam at supercritical incidence. The reflection coefficient of the icewater interface with a cavity in this paper is marked by , and the reflection coefficient without a cavity is marked by .
The longrange propagation modes in an acoustic channel under ice are basically caused by supercritical incidence, and the scattering of a nearinterface cavity in ice changes the spatial distribution of the acoustic energy under ice and thereby affects the transmission loss. Therefore, this paper focuses on the acoustic characteristics of reflection coefficient at supercritical incidence.
3.1. Angular Characteristics of R
The amplitude of inhomogeneous plane wave in ice attenuates with depth at supercritical incidence, and the acoustic energy returns to water after propagating for a distance along the horizontal direction. Therefore, the shallow surface on the cavity plays a major role in the scattering field, and a variable crosssection waveguide is formed by a section of the shallow surface on the cavity and the icewater interface, as shown in Figure 6.
During the propagation of the inhomogeneous plane wave in the variable crosssection waveguide for a certain distance, the horizontal wave number changes when the inhomogeneous plane wave returns to water, which causes a deflection of the reflection direction.
It is very difficult to find an analytical solution of the acoustic field in the variable crosssection waveguide, so some reasonable approximations are introduced: (1) The longitudinal component of the inhomogeneous plane wave attenuates much more with depth, so only the contribution of transverse component in the inhomogeneous plane wave is considered, and the mode conversion of elastic wave on the waveguide interface is ignored; (2) the upper boundary of the variable crosssection waveguide is approximated to a horizontal step shape, assuming that the acoustic field in the waveguide can be expressed by the following form with introducing a minor error:
The incident angle at supercritical incidence is very large on icewater interface and the beamwidth is small. It can be approximately considered as that , so the beamwidth is equal to
The beamwidth depends on the incident angle and the horizontal wave number variation , which is not significant compared with if the geometry parameters and incident frequency are constants. The relation between beamwidth and incident angle is shown in Figure 7.
As can be seen from Figure 7, the relationship between the beamwidth and the incident angle predicted by (7) is consistent with the trend of numerical results. There are some differences between the theoretical prediction and the numerical results because the incident angle and the scattering angle are discrete values.
Assume that the effective horizontal area of the variable crosssection waveguide to deflect the reflection is and the sound pressure in the beam is approximately sinusoidal and its peak appears in the mirror reflection direction. It can be known from the conservation of energy that
where is sound intensity and , in farfield. The wave impedance of the reflected wave equals the characteristic impedance of the medium if the wave propagates to a sufficient distance in the mirror reflection direction, and the expression of reflection coefficient can be written as follows:
It is found that the effective horizontal area and horizontal wave velocity variation are not significant compared with , and the reflection coefficient is approximately proportional to the cosine of incident angle, as shown in Figure 8.
As can be seen from Figure 8, the reflection coefficients and are obviously different; especially is less than 1 at supercritical incidence, and it decreases with the increase of the incident angle. The predicted trend of at supercritical incidence according to (9) is in good agreement with the actual result, which is consistent with the phenomenon that the beamwidth increases with the incidence angle at supercritical angle, as shown in Figure 5.
In summary, the effective horizontal area and horizontal wave velocity variation is not significant compared with the incident angle if the geometry parameters and incident frequency are constants, and the beamwidth is proportional to secant of incident angle, while the reflection coefficient is proportional to cosine of incident angle.
3.2. Frequency Characteristics of R
In practical application, the equipment and operations required to measure the beamwidth are very complicated, while it is easier to measure the frequency characteristic of reflection coefficient at a fixed incident angle.
The radical formula in (9) is defined as , which characterizes the influence of geometry parameters and incident frequency in the reflection coefficient and its relation with the relative depth and Helmholtz number is shown in Figure 9.
(a)
(b)
(c)
(d)
(e)
As can be seen from Figure 9, and incident angle are hardly irrelevant, and increases with Helmholtz number. The variation of horizontal wave number is approximated to the first term of its Taylor series, deriving , where is the spatial derivative of horizontal wave number to horizontal location and is the horizontal migration length of inhomogeneous plane wave in the variable crosssection waveguide. The maximum of depends on , and an implicit function in appendix (A.8) is determined. Its approximate explicit function has the expression because of . It could be proved that , , and dimensions are measured in terms of the wavelength in acoustics, so the expression of reflection coefficient is expressed in the dimensionless Helmholtz number and relative depth:
where
The numerical results and theoretical predictions of are shown in Figure 10. The numerical results are marked with different linetypes, and the solid lines with corresponding colors are the theoretical predictions.
As can be seen from Figure 10, the theoretical predictions are in very good agreement with the numerical results. increases with the relative depth ; that is, the deflection effect of the variable crosssection waveguide and the beamwidth of reflected wave become smaller if the ice cavity is away from the icewater interface, and the reflection coefficient correspondingly becomes greater. This feature can also be seen through the formula (10) and (11), which means the scattering of the ice cavity becomes weaker with the increasing relative depth and incident frequency.
4. Conclusions
In this paper, the FourierBessel series and SommerfeldWatson transformation method are used to analyze the influence of a slender cylindrical cavity near and parallel to an icewater interface on the sound propagation in an acoustic channel under ice. Through research, it is found that, due to the joint action of the cavity and the icewater interface, the underwater acoustic field presents a fast fluctuating interference pattern, and what is more, the acoustic field in the mirror reflection direction has the following characteristics at supercritical incidence:
(a) The beamwidth is proportional to secant of incident angle.
(b) The reflection coefficient is proportional to cosine of incident angle.
(c) The reflection coefficient increases with Helmholtz number and relative depth if incident angle is a constant.
The scattering of the inhomogeneous plane wave in ice by a cavity changes the spatial distribution of underwater acoustic energy, making the transmission loss during longdistance communication much larger than that of the icewater interface without a cavity. The research provides a theoretical basis for an indepth understanding of the acoustic channel under ice, as well as supporting the detection and identification of underwater targets and communication technology under ice.
Appendix
(1) Boundary Continuous Equation (Matrix Form). Normal stress continuity equation:
Tangential stress continuity equation:
Normal displacement continuous equation:wherein is the density of water and is the second parameter of the Lamé constant of ice. The matrix element expression is as follows:whereinand
wherein
(2) SommerfeldWatson Transformation for the Scattering Field under Ice:
wherein
(3) SommerfeldWatson Transformation for Boundary Continuity Equations:wherein
(4) The Relation between the BeamWidth and the Horizontal Wave Number :
wherein is the polar angle in the coordinate system and . The equation corresponding to beamwidth is , obtaining a quintic equation:
wherein
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research is finally supported by the National Key Research and Development Plan (2018YFC1405900), the Key Project of National Natural Science Foundation (61631008), and the Special Fund for Basic Scientific Research Service Fees of Central Universities (heu10500170006).
References
 J. Gobat and C. Lee, “Acoustic navigation and communication for highlatitude ocean research,” AGU Fall Meeting Abstracts, 2006. View at: Google Scholar
 A. K. Morozov, L. Freitag, and S. Singh, “Underwater acoustic technologies for longrange navigation and communications in the Arctic,” in Proceedings of the Undersea Acoustic Measurement Conference, 2011. View at: Google Scholar
 G. Jin, J. F. Lynch, and R. Pawlowicz, “Effects of sea ice cover on acoustic ray travel times, with applications to the Greenland Sea tomography experiment,” Journal of the Acoustical Society of America, vol. 94, no. 94, pp. 1044–1057, 1993. View at: Google Scholar
 P. N. Mikhalevsky, A. B. Baggeroer, A. Gavrilov, and M. Slavinsky, “Experiment tests use of acoustics to monitor temperature and ice in Arctic Ocean,” Eos, Transactions, American Geophysical Union, vol. 76, no. 27, pp. 265–269, 1995. View at: Publisher Site  Google Scholar
 P. N. Mikhalevsky, A. N. Gavrilov, and A. B. Baggeroer, “Transarctic acoustic propagation experiment and climate monitoring in the Arctic,” IEEE Journal of Oceanic Engineering, vol. 24, no. 2, pp. 183–201, 1999. View at: Publisher Site  Google Scholar
 F. Desharnais, G. J. Heard, G. R. Ebbeson, D. J. Thomson, D. J. Schillinger, and R. Seyan, “Acoustic propagation in an arctic shallowwater environment,” in Proceedings of the Oceans 2003  Marine Technology Society, Institute of Electrical and Electronics Engineers Conference, MTS/IEEE, pp. 10–17, USA, September 2003. View at: Google Scholar
 H. Schmidt and T. Schneider, “Acoustic communication and navigation in the new Arctic — A model case for environmental adaptation,” in Proceedings of the 2016 IEEE Third Underwater Communications and Networking Conference (UComms), pp. 1–4, Lerici, Italy, August 2016. View at: Publisher Site  Google Scholar
 L. Freitag, P. Koski, A. Morozov, S. Singh, and J. Partan, “Acoustic communications and navigation under Arctic ice,” in Proceedings of the OCEANS 2012 MTS/IEEE Hampton Roads Conference: Harnessing the Power of the Ocean, IEEE, Hampton Roads, Va, USA, October 2012. View at: Google Scholar
 L. Freitag, K. Ball, J. Partan, P. Koski, and S. Singh, “Long range acoustic communications and navigation in the Arctic,” in Proceedings of the MTS/IEEE Washington, OCEANS 2015, IEEE, Washington, DC, USA, October 2015. View at: Google Scholar
 L. Freitag, P. Koski, S. Singh, T. Maksym, and H. Singh, “Acoustic Communications under Shallow ShoreFast Arctic Ice,” in Proceedings of the OCEANS 2017  Anchorage, pp. 1–5, IEEE, Anchorage, Alaska, USA, September 2017. View at: Google Scholar
 L. Godinho, P. AmadoMendes, and A. Pereira, “A Hybrid analyticalnumerical model based on the method of fundamental solutions for the analysis of sound scattering by buried shell structures,” Mathematical Problems in Engineering, vol. 2011, Article ID 710623, 22 pages, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 J. L. Fleming and J. Moser, “A fourier approach to model electromagnetic fields scattered by a buried rectangular cavity,” Mathematical Problems in Engineering, vol. 2009, Article ID 937327, 22 pages, 2009. View at: Google Scholar
 Q.F. He, Q. Zhang, Y. Luo, and L. Sun, “Sinusoidal frequency modulation fourierbessel series for multicomponent sfm signal estimation and separation,” Mathematical Problems in Engineering, 2017. View at: Google Scholar
 C. H. Lin, V. W. Lee, M. I. Todorovska, and M. D. Trifunac, “Zerostress cylindrical wave functions around a circular underground tunnel in a flat elastic halfspace: incident P waves,” Soil Dynamics and Earthquake Engineering, vol. 30, no. 10, pp. 879–894, 2010. View at: Google Scholar
 V. W. Lee and W.Y. Liu, “Twodimensional scattering and diffraction of P and SVwaves around a semicircular canyon in an elastic halfspace: an analytic solution via a stressfree wave function,” Soil Dynamics and Earthquake Engineering, vol. 63, pp. 110–119, 2014. View at: Publisher Site  Google Scholar
 C. A. Davis, V. W. Lee, and J. P. Bardet, “Transverse response of underground cavities and pipes to incident SV waves,” Earthquake Engineering & Structural Dynamics, vol. 30, no. 3, pp. 383–410, 2010. View at: Google Scholar
 J. Pumplin, “Application of sommerfeldwatson transformation to an electrostatics problem,” American Journal of Physics, vol. 37, no. 7, pp. 737–739, 1969. View at: Publisher Site  Google Scholar
 M. L. Rumerman, “Increased accuracy in the application of the SommerfeldWatson transformation to scattering of acoustic waves from spherical shells,” The Journal of the Acoustical Society of America, vol. 94, no. 5, pp. 2502–2509, 1993. View at: Publisher Site  Google Scholar  MathSciNet
 M. L. Rumerman, “Increased accuracy in the application of the SommerfeldWatson transformation to scattering of acoustic waves from spherical shells,” The Journal of the Acoustical Society of America, vol. 94, no. 5, pp. 2739–2750, 1993. View at: Publisher Site  Google Scholar  MathSciNet
 W. Fan, G. Y. Zheng, and J. Fan, “Analysis of scattering wave from a waterfilled cylindrical shell using SommerfeldWatson transformation,” Technical Acoustics, vol. 28, no. 8, pp. 142–144, 2009. View at: Google Scholar
Copyright
Copyright © 2019 WenKai Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.