Theory of complex-coordinate transformation acoustics for non-Hermitian metamaterials

Transformation acoustics(TA)has emerged as a powerful tool for designing several intriguing conceptual devices,which can manipulate acoustic waves in a flexible manner,yet their applications are limited in Hermitian materials.In this work,we propose the theory of complex-coordinate transformation acoustics(CCTA)and verify the effectiveness in realizing acoustic non-Hermitian metamaterials.Especially,we apply this theory for the first time to the design of acoustic parity-time(PT)and antisymmetric parity-time(APT)metamaterials and demonstrate two distinctive examples.First,we use this method to obtain the exceptional points(EPs)of the PT/APT system and observe the spontaneous phase transition of the scattering matrix in the transformation parameter space.Second,by selecting the Jacobian matrix's constitutive parameters,the PT/APT-symmetric system can also be configured to approach the zero and pole of the scattering matrix,behaving as an acoustic coherent perfect absorber and equivalent laser.We envision our proposed CCTAbased paradigm to open the way for exploring the non-Hermitian physics and finding application in the design of acoustic functional devices such as absorbers and amplifiers whose material parameters are hard to realize by using the conventional transformation method.

High-order Bragg forward scattering and frequency shift of low-frequency underwater acoustic field by moving rough sea surface

Acoustic scattering modulation caused by an undulating sea surface on the space-time dimension seriously affects underwater detection and target recognition.Herein,underwater acoustic scattering modulation from a moving rough sea surface is studied based on integral equation and parabolic equation.And with the principles of grating and constructive interference,the mechanism of this acoustic scattering modulation is explained.The periodicity of the interference of moving rough sea surface will lead to the interference of the scattering field at a series of discrete angles,which will form comb-like and frequency-shift characteristics on the intensity and the frequency spectrum of the acoustic scattering field,respectively,which is a high-order Bragg scattering phenomenon.Unlike the conventional Doppler effect,the frequency shifts of the Bragg scattering phenomenon are multiples of the undulating sea surface frequency and are independent of the incident sound wave frequency.Therefore,even if a low-frequency underwater acoustic field is incident,it will produce obvious frequency shifts.Moreover,under the action of ideal sinusoidal waves,swells,fully grown wind waves,unsteady wind waves,or mixed waves,different moving rough sea surfaces create different acoustic scattering processes and possess different frequency shift characteristics.For the swell wave,which tends to be a single harmonic wave,the moving rough sea surface produces more obvious high-order scattering and frequency shifts.The same phenomena are observed on the sea surface under fully grown wind waves,however,the frequency shift slightly offsets the multiple peak frequencies of the wind wave spectrum.Comparing with the swell and fully-grown wind waves,the acoustic scattering and frequency shift are not obvious for the sea surface under unsteady wind waves.

Scheme of negative acoustic radiation force based on a multiple-layered spherical structure

Acoustic radiation force(ARF), as an important particle manipulation method, has been extensively studied in recent years. With the introduction of the concept of “acoustic tweezers”, negative acoustic radiation has become a research hotspot. In this paper, a scheme of realizing negative ARF based on the multiple-layered spherical structure design is proposed. The specific structure and design idea are presented. Detailed theoretical calculation analysis is carried out.Numerical simulations have been performed to verify the correctness of this prediction. The conjecture that the suppression of backscattering can achieve negative ARF is verified concretely, which greatly expands the application prospect and design ideas of the ARF. This work has laid a theoretical foundation for realizing precise control of the structure.

Pipeline thickness estimation using the dispersion of higher-order SH guided waves

Thickness measurement plays an important role in the monitoring of pipeline corrosion damage. However, the requirement for prior knowledge of the shear wave velocity in the pipeline material for popular ultrasonic thickness measurement limits its widespread application. This paper proposes a method that utilizes cylindrical shear horizontal(SH) guided waves to estimate pipeline thickness without prior knowledge of shear wave velocity. The inversion formulas are derived from the dispersion of higher-order modes with the high-frequency approximation. The waveform of the example problems is simulated using the real-axis integral method. The data points on the dispersion curves are processed in the frequency domain using the wave-number method. These extracted data are then substituted into the derived formulas. The results verify that employing higher-order SH guided waves for the evaluation of thickness and shear wave velocity yields less than1% error. This method can be applied to both metallic and non-metallic pipelines, thus opening new possibilities for health monitoring of pipeline structures.

Energy Walk-off in Nonlinear Acoustics

It is found that acoustic internal conical refraction occurs for both fundamental and second harmonic of transverse wave propagating dong threefold axis of crystals.Energy walk-off which usually does not exist in acoustic collinear harmonic generation appears as well.

REPRESENTATION THEOREM OF H^(P) FUNCTIONS AND THEIR APPLICATION TO ACOUSTICS

The correct analytic continuation formulae of H^(P) functions on the upper half-plane are given and used to extrapolate the transfer function of a spherical piezoelectric receiver.The results obtained by extrapolation are in good agreement with the theoretical ones.

Response characteristics of drill-string guided wave in downhole acoustic telemetry

Modeling of a drill-string acoustic channel has been an important topic in downhole telemetry for a long time.The propagation of drill-string guided waves in the borehole contains excitation,attenuation,and mode conversion issues that have not been considered by existing modeling methods.In this article,we formulate a hybrid modeling method to investigate the response characteristics of a fundamental-mode drill-string wave in various borehole environments.This hybrid method provides channel functions,including transmitting and receiving deployments,periodicity of the structure,and formation property changes.The essential physics of the drill-string wave propagation is captured with a one-dimensional model.The analytical solutions of the wavefield in multilayered cylindrical structures are introduced into a propagation matrix to express drill-string-wave interactions with the borehole environments.The effectiveness of the proposed method is confirmed through comparison with the finite-difference method.In addition,by designing numerical models,we investigate the conversion effect of the drill-string wave at the tool joint.We demonstrate that the conversion intensity of the drill-string wave is positively correlated not only with the cross-sectional area of the tool joint but also with the wave impedance of the outer formation.Hard formation outside the borehole reduces the energy leakage while intensifying the conversion of drill-string waves to Stoneley waves,and the opposite is true for the drill string in an infinite fluid.The converted Stoneley waves interfere with the drill-string waves,resulting in variations of bandgap distribution,which challenges the reliability of the data transmission.

Sound absorbing properties of spiral metasurfaces inspired by micro-perforated plates

As a kind of classical low-frequency sound-absorbing material,the microperforated plate(MPP)has been widely used.Here,we inspired by the sound absorption mechanism of the MPP,a spiral metasurface(SM)is designed and the analytical solution of acoustic impedance and sound absorption coefficient are obtained.The relationship between the sound absorption properties of the MPP and the SM with their own structures is systematically studied,and the analytical solutions are used to optimise the structure.It is concluded that the MPP and the SM of the same thickness achieve effective absorption in the frequency range between 390-900 Hz and 1920-4266 Hz,with a total thickness less than 1/6 of the wavelength.Meanwhile,the numerical calculation shows that the MPP and SM can match well with the background medium in the effective rang.Our study provides new insights into the design methods of sound-absorbing materials and is potentially suitable for many acoustic engineering applications.

Basis functions for shallow-water temperature profiles based on the internal-wave eigenmodes

The shallow-water temperature profile is typically parameterized using a few empirical orthogonal function(EOF)coefficients.However,when the experimental area is poorly known or highly variable,the adaptability of the EOFs will be significantly reduced.In this study,a new set of basis functions,generated by combining the internal-wave eigenmodes with the average temperature gradient,is developed for characterizing the temperature perturbations.Temperature profiles recorded by a thermistor chain in the South China Sea in 2015 are processed and analyzed.Compared to the EOFs,the new set of basis functions has higher reconstruction accuracy and adaptability;it is also more stable in ocean regions that have internal waves.

Characteristic analysis of scattering field in two-layer media by Green's function

The problem of three-dimensional(3D) acoustic scattering in a complex medium has aroused considerable interest of researchers for many years. An ultrasonic scattered field calculating technique is proposed to study the scattering echo from strongly scattered materials in a two-layer medium in this work. Firstly, with the high frequency stationary phase method,the Green's function of two-layer fluid media is derived. And then based on the idea of integral equation discretization,the Green's function method is extended to two-layer fluid media to derive the scattering field expression of defects in a complex medium. With this method, the scattering field of 3D defect in a two-layer medium is calculated and the characteristics of received echoes are studied. The results show that this method is able to solve the scattering P wave field of 3D defect with arbitrary shape at any scattering intensity in two-layer media. Considering the circumstance of waterimmersion ultrasonic non-destructive test(NDT), the scattering sound field characteristics of different types of defects are analyzed by simulation, which will help to optimize the detection scheme and corresponding imaging method in practice so as to improve the detection quality.

Unidirectional Negative Refraction at an Exceptional Point of Acoustic P T-Symmetric Systems

We demonstrate a method to realize unidirectional negative refraction in an acoustic parity-time(P T)-symmetric system, which is composed of a pair of metasurfaces sandwiching an air gap. The pair of metasurfaces possesses loss and gain modulations. The unidirectional negative refraction, which is strictly limited to the case of incident wave imposing on the loss end of the metasurface, is demonstrated at the exception point(EP) in this P T-symmetric system, while the incidence from the other side leads to strong reflection. Based on rigorous calculations, we explicitly show the underlying mechanism of this model to achieve unidirectional wave scatterings around the EP in the parametric space. In addition, the perfect imaging of a point source in the three-dimensional space, as a signature of negative refraction, is simulated to provide a verification of our work. We envision that this work may sharpen the understanding of P T-symmetric structures and inspire more acoustic functional devices.

Distribution of chromophytic phytoplankton in the Western Subarctic Gyre of Pacific Ocean revealed by morphological observation and rbc L gene sequences

Western Subarctic Gyre(WSG),which possesses distinctive differences in oceanographic and biogeochemical processes,is situated in the northwest subarctic Pacific.The WSG is characterized by high nutrient and low chlorophyll.We carried out a field investigation in this area in summer 2020 and performed microscopic observation,cytometric counting,and RuBisCO large subunit(rbc L)gene analysis to understand the community structure and spatial distribution of chromophytic phytoplankton better.Microscopic method revealed that total phytoplankton(>10μm,including Bacillariophyta,Dinoflagellata,Ochrophyta,and Chlorophyta)abundances ranged(0.6×10^(3))-(167.4×10^(3))cells/L with an increasing trend from south to north.Dinoflagellates and Pennatae diatoms dominated the phytoplankton assemblages in the southern and northern stations,respectively.Major chromophytic phytoplankton groups derived from rbc L genes included Haptophyta,Ochrophyta,Bacillariophyta,as well as rarely occurring groups,such as Xanthophyta,Cyanobacteria,Dinoflagellata,Rhodophyta,and Cryptophyta.At the phylum level,Haptophyta was the most abundant phylum,accounting for approximately 30.80%of the total obtained operational taxonomic units in all samples.Ochrophyta and Bacillariophyta were the second and third most abundant phylum,and their relative abundance was 20.26% and 19.60%,respectively.Further,redundancy analysis showed that high proportion of diatoms(e.g.,microscopic and rbc L methods)was positively correlated with nutrients(e.g.,dissolved inorganic nitrogen(DIN),dissolved inorganic phosphorous,and dissolved silicate(DSi))and negatively correlated with temperature and salinity.The proportion of Ochrophyta,Rhodophyta,and Cyanobateria identified by rbc L genes was positively correlated with salinity and temperature and showed negative correlation to nutrients.This work is the first molecular study of phytoplankton accomplished in the WSG,and our results show some discrepancies between morphological observation and rbc L gene sequences,which highlight the necessity of combining the microscopic and molecular methods to reveal the diversity of phytoplankton in marine environment.

A preliminary study on the acoustic properties of seafloor sediment in the southern U-boundary of the South China Sea

The acoustic properties of seafloor sediment are essential parameters in the exploration of marine resources,ocean scientific research and ocean engineering.Seafloor sediment samples were collected at the southern U-boundary of the South China Sea(SCS),and the acoustic and physical properties were measured in the laboratory.The correlation between physical and sound speed ratio(SSR)was discussed,and SSR-physical property empirical regressions in the Sunda Shelf were established for the first time.Compared with the northern continental shelf of SCS,the Sunda Shelf are mainly silty and sand sediment,and the SSR ranges from 0.9949 to 1.0944,which has higher SSR than the northern continental shelf,implies that the Sunda Shelf is a high SSR area.Since the same kind of sediment has different physical properties,the single physical parameter of sediment cannot fully represent the acoustic properties of sediment,therefore,the multiple parameter prediction model should develop in the future to improve the prediction precision.

Robust Space-Time Adaptive Track-Before-Detect Algorithm Based on Persymmetry and Symmetric Spectrum

Underwater monopulse space-time adaptive track-before-detect method,which combines space-time adaptive detector(STAD)and the track-before-detect algorithm based on dynamic programming(DP-TBD),denoted as STAD-DP-TBD,can effectively detect low-speed weak targets.However,due to the complexity and variability of the underwater environment,it is difficult to obtain sufficient secondary data,resulting in a serious decline in the detection and tracking performance,and leading to poor robustness of the algorithm.In this paper,based on the adaptive matched filter(AMF)test and the RAO test,underwater monopulse AMF-DP-TBD algorithm and RAO-DP-TBD algorithm which incorporate persymmetry and symmetric spectrum,denoted as PSAMF-DP-TBD and PS-RAO-DP-TBD,are proposed and compared with the AMF-DP-TBD algorithm and RAO-DP-TBD algorithm based on persymmetry array,denoted as P-AMF-DP-TBD and P-RAO-DP-TBD.The simulation results show that the four methods can work normally with sufficient secondary data and slightly insufficient secondary data,but when the secondary data is severely insufficient,the P-AMF-DP-TBD and P-RAO-DP-TBD algorithms has failed while the PSAMF-DP-TBD and PS-RAO-DP-TBD algorithms still have good detection and tracking capabilities.

基于吉布斯采样的稀疏水声信道估计方法

The estimation of sparse underwater acoustic(UWA)channels can be regarded as an inference problem involving hidden variables within the Bayesian framework.While the classical sparse Bayesian learning(SBL),derived through the expectation maximization(EM)algorithm,has been widely employed for UWA channel estimation,it still differs from the real posterior expectation of channels.In this paper,we propose an approach that combines variational inference(VI)and Markov chain Monte Carlo(MCMC)methods to provide a more accurate posterior estimation.Specifically,the SBL is first re-derived with VI,allowing us to replace the posterior distribution of the hidden variables with a variational distribution.Then,we determine the full conditional probability distribution for each variable in the variational distribution and then iteratively perform random Gibbs sampling in MCMC to converge the Markov chain.The results of simulation and experiment indicate that our estimation method achieves lower mean square error and bit error rate compared to the classic SBL approach.Additionally,it demonstrates an acceptable convergence speed.

Matched-field inversion of sound speed profile in shallow water using a parallel genetic algorithm

A sound speed profile plays an important role in shallow water sound propagation.Concurrent with in-situ measurements,many inversion methods,such as matched-field inversion,have been put forward to invert the sound speed profile from acoustic signals.However,the time cost of matched-field inversion may be very high in replica field calculations.We studied the feasibility and robustness of an acoustic tomography scheme with matched-field processing in shallow water,and described the sound speed profile by empirical orthogonal functions.We analyzed the acoustic signals from a vertical line array in ASIAEX2001 in the East China Sea to invert sound speed profiles with estimated empirical orthogonal functions and a parallel genetic algorithm to speed up the inversion.The results show that the inverted sound speed profiles are in good agreement with conductivity-temperature-depth measurements.Moreover,a posteriori probability analysis is carried out to verify the inversion results.

Thickness Determination for a Two-Layered Composite of a Film and a Plate by Low-Frequency Ultrasound

We present an ultrasonic method for determining the thickness of a composite consisting of a soft thin film attached to a hard plate substrate, by resonance spectra in the low frequency region, The interrogating waves can be incident only to the two-layered composite from the substrate side. The reflection spectra are obtained by FFT analysis of the compressive pulsed echoes from the composite, and the thicknesses of the film and the substrate are simultaneously inversed by the simulated annealing method from the resonant frequencies knowing other acoustical parameters in prior. The sensitivity of the method to individual thickness, its convergence and stability against experimental noises are studied, Experiment with interrogating wavelength 4 times larger than the film thickness in a sample of a polymer film （0.054mm） on an aluminium plate （6.24mm） verifies the validity of the method. The average relative errors in the measurement of the thicknesses of the film and the substrate are found to be -4.1% and -0.62%, respectively.

Effects of rough surface on sound propagation in shallow water

Underwater acoustic applications depend critically on the prediction of sound propagation, which can be significantly affected by a rough surface, especially in shallow water. This paper aims to investigate how randomly fluctuating surface influences transmission loss(TL) in shallow water. The one-dimension wind-wave spectrum, Monterey–Miami parabolic equation(MMPE) model, Monte Carlo method, and parallel computing technology are combined to investigate the effects of different sea states on sound propagation. It is shown that TL distribution properties are related to the wind speed,frequency, range, and sound speed profile. In a homogenous waveguide, with wind speed increasing, the TLs are greater and more dispersive. For a negative thermocline waveguide, when the source is above the thermocline and the receiver is below that, the effects of the rough surface are the same and more significant. When the source and receiver are both below the thermocline, the TL distributions are nearly the same for different wind speeds. The mechanism of the different TL distribution properties in the thermocline environment is explained by using ray theory. In conclusion, the statistical characteristics of TL are affected by the relative roughness of the surface, the interaction strength of the sound field with the surface, and the changes of propagating angle due to refraction.

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses. In the finite element model, the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method, and the dynamical interaction between both crack faces is modeled using a contact-impact theory. In the simulations, the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis, especially, the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail. Meanwhile, the related ultrasonic infrared images are also obtained experimentally, and the theoretical simulation results are in agreement with that of the experiments. The results show that, by using the theoretical method, a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.

Geoacoustic Inversion for Bottom Parameters in the Deep-Water Area of the South China Sea

Bottom acoustic parameters play an important role in sound field prediction. Acoustic parameters in deep water are not well understood. Bottom acoustic parameters are sensitive to the transmission-loss （TL） data in the shadow zone of deep water. We propose a multiple-step fill inversion method to invert sound speed, density and attenuation in deep water. Based on a uniform liquid hMf-space bottom model, sound speed of the bottom is inverted by using the long range TL at low frequency obtained in an acoustic propagation experiment conducted in the South China Sea （SCS） in summer 2014. Meanwhile, bottom density is estimated combining with the Hamilton sediment empirical relationship. Attenuation coefficients at different frequencies are then estimated from the TL data in the shadow zones by using the known sound speed and density as a constraint condition. The nonlinear relationship between attenuation coefficient and frequency is given in the end. Tile inverted bottom parameters can be used to forecast the transmission loss in the deep water area of SCS very we//.