Wave Propagation Research Articles

A novel hybrid SBM-MFS methodology for acoustic wave propagation problems

In this paper, a novel hybrid meshless approach that combines the singular boundary method (SBM) and the method of fundamental solutions (MFS) to deal with two-dimensional (2D) exterior acoustic wave propagation problems is proposed and studied. The methodology is particularly devised to solve problems with complex boundary geometries containing geometric singularities such as corners and sharp edges. It employs the SBM to model intricate segments of these geometries and the MFS for the smooth ones. The proposed hybrid SBM-MFS method is studied in a 2D context in the framework of three benchmark examples involving acoustic radiation problems of circular-, square- and L-shaped objects in a full-space acoustic medium. In addition, the applicability of the proposed hybrid SBM-MFS methodology to predict the acoustic performance of a T-shaped thin barrier is also investigated. These examples are specifically designed to assess the feasibility, validity and accuracy of the hybrid SBM-MFS approach in comparison with the available analytical solutions and alternative numerical strategies such as the MFS, the SBM and the boundary element method (BEM). Numerical simulations demonstrate that the proposed method can match the level of accuracy of the MFS while keeps the robustness of the SBM when dealing with complex geometries, overcoming one the most important drawbacks of the traditional MFS. Moreover, the proposed hybrid SBM-MFS naturally avoids the non-uniqueness problem arising at the fictitious eigenfrequencies associated with the corresponding interior problems, a feature that neither the SBM nor the BEM inherently possesses.

Anti-scattering propagation in multiple-bend valley phononic crystals

Valley topological phononic crystals (PCs) have attracted wide attention due to the topological properties of their edge states. In general, valley interface states can exist in the interfaces that are constructed by opposite valley topological phases. Here we study the anti-scattering propagation properties of edge states in a single valley PC. We present that the edge states can exist in different boundary terminations with different band dispersions. The boundary transport behaviors of acoustic waves along the two designed PCs are demonstrated numerically. The results show that the chiral edge states are immune against additional scatterers that preserve the valley pseudospins, but the backscattering can happen when intervalley scattering is included. Nevertheless, the anti-scattering propagation in complex multiple-bend structures can be realized by the smooth transition between the edge states and the valley interface states. Similar to the designed frequency-selective device, more prospective applications can be anticipated in the manipulation of acoustic wave propagation.

Initial Experimental Results of Producing Multicharged Ions Efficiently by Lower Hybrid Resonance Heating with Exciting Helicon Waves on Electron Cyclotron Resonance Ion Source

Based on experimental results, we considered the accessibility conditions for wave propagation in the magnetized plasma in an electron cyclotron resonance (ECR) ion source (ECRIS), and proposed introducing electromagnetic waves with frequencies lower than those of ECR. We generate the lower hybrid resonance (LHR) in the ECRIS by the electric field of the extra-ordinary mode (X-mode) on the helicon wave which additionally heat the electrons, and then it was conducted to improve the efficiency of multicharged ion generation successfully in over density condition. The electron energy distribution function obtained by the probe measurement shows an increase in the high-energy region in the corresponding LHR region. In this paper, for the first time we describe the initial experimental results on enhanced production of multicharged ions by the LHR of the low-frequency RF electromagnetic waves into the ECRIS.

Wave Reflection in the Solar Atmosphere

We present evidence supporting wave reflection in the lower solar chromosphere based on helioseismic analysis of multi-height Doppler data from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager and the Magneto-Optical filters at Two Heights II instrument. This evidence is derived through a wave propagation model that incorporates both upward- and downward-traveling (reflected) waves. Moreover, we find that the height of the reflecting region varies with magnetic field strengths in a way that suggests a connection with the plasma β ∼ 1 region. We measure an effective reflection coefficient of 13% in a magnetically quiet region of the Sun.

Flexural-gravity wave interaction with undulating bottom topography in the presence of uniform current: An asymptotic approach

Using an asymptotic method, this article deals with flexural-gravity wave scattering with undulating bottom topography, including the effect of uniform currents. The interest in this problem lies in developing second-order solutions using the Fourier transform, which minimises the error gap between first and second-order solutions. The present method allows the physical processes involved in the sea-bed topography, uniform current, plate-covered surface, and wave interaction to be studied. Specifically, we observe Bragg resonance between the flexural-gravity waves and the bottom ripples, which are associated with the reflection of incident wave energy. We examine the effects of wave current and emphasise how crucial the asymptotic expansion method is to the emergence of the current response. We demonstrate that bottom topography dominates the effects of Bragg resonance for depth Froude numbers valued at 0.8 or less. Further, most reflected wave components have their frequencies shifted by the current, and wave action conservation causes reflected wave energy to be enhanced for following currents. Using the Joint North Sea Wave Observation Project spectrum and the discrete Fourier transform, the theory derived in the frequency domain is shown in the time domain to analyse wave propagation through the whole system.

Band gap study of periodic piezoelectric micro-composite laminated plates by finite element method and its application in feedback control

Phononic crystals are materials or structures with periodic variations of elastic constants and density. Elastic waves in phononic crystals cannot propagate within a certain frequency range, called the band gap, due to the interference of their internal structure. This article investigates how fiber orientation affects the band gap properties of periodic micro-composite laminated plates. We examined the impact of five common composite laminate configurations on the system’s band gap and discovered that the antisymmetric configuration had superior high-frequency band gap performance. Based on this finding, we further explored how the antisymmetric angle influenced the system’s band gap. Since the study was conducted at a micro scale, micro structural effects were considered. Therefore, we used the modified coupled stress elastic dynamics and a four-node quadrilateral non-conforming element to discretize the micro composite laminated plate model, which had the nodal compatibility of high-order elasticity theory. Moreover, we employed feedback control to the proposed structure to dynamically adjust the band gap width according to different practical needs. This research supports the design of periodic micro-composite laminated plates and the piezoelectric feedback control system, which can control the vibration and wave propagation behaviors of micro structures using coupled piezoelectric sensors and actuators. The piezoelectric feedback control employed two control methods: direct proportional control and acceleration control. This article compared the single and multiple control methods of these two approaches to determine the optimal control strategy.

Effect of bathymetry data on tsunami wave ray tracing in the western Banten sea

The lack of tsunami wave warning during the Anak Krakatau collapse in 2018 resulted in devastating damage to coastal areas. The tsunami wave arrived in the coastal areas approximately half an hour after the collapse. As a tsunami wave can travel abruptly, an early warning system should warn faster. However, a warning system based on a sophisticated hydrodynamic model in real-time would take time to conduct the numerical tsunami wave simulation. Hence, this study proposes a fast and reliable estimation of tsunami wave propagation through a classic ray tracing analysis. We use two ray-tracing methods to investigate the tsunami wave propagation from Anak Krakatau to the western Banten Sea. The first method follows Snell’s law, considering dispersive waves, and the second assumes non-dispersive waves based on the ray tracing equations, considering the Earth’s sphericity. Both methods are quantitatively evaluated by comparing the travel time measured at Anyer and Marina Jambu. This study finds that non-dispersive wave tracing performs a shorter computational time and slightly better prediction of tsunami wave propagation than dispersive-based wave tracing, with a relative absolute difference of the travel time of 17.9–26.7% in Anyer and 3.6–5.2% in Marina Jambu. This study also confirms the importance of bathymetry validity in wave ray tracing. Two regional bathymetry datasets with a mean difference of less than 5 m result in different wave ray tracing, in which one dataset does not produce the wave ray path toward the Panaitan Island. Based on bathymetric surveys in Anyer and Marina Jambu, the national bathymetry data (BATNAS) has shown its superiority to being used as a bathymetry in the ray-tracing process, with correlations of 81% and 93% in both areas compared to other available bathymetry datasets. We summarize that reliable bathymetry data and the non-dispersive ray tracing method can be used as an initial estimation of tsunami wave propagation efficiently.

Analysis of shock wave propagation in two-layered blood flow model via Lie symmetry

In this paper, we study the two-layered blood flow model within the framework of Lie symmetry analysis to investigate a range of closed-form solutions. Our investigation involves optimal classification, revealing the existence of six optimal subalgebras within the model. Further, we streamline the model by reducing it into a system of ordinary differential equations (ODEs), thereby deriving several exact solutions. Our analysis establishes the presence of the nonlinear self-adjointness property in the governing equation, leading to the development of various conservation laws. Furthermore, leveraging one of the exact solutions we obtained, we delve into the behaviour of characteristic shocks, C1 waves, and their interactions, providing a comprehensive understanding of their dynamics.

Research on the upper limit of constant strain rates and dynamic fracture behavior for hybrid woven composite

Based on the analysis of stress wave propagation at both ends of the specimen, a formula for the upper limit of the constant strain rate is derived. The upper limit of constant strain rate is applicable to all the ratios of wave impedance, providing a broader applicability. The in-plane dynamic compressive behavior was obtained for 2D plain weave carbon/aramid hybrid composites at strain rates ranging from 195-790 s−1. The results validate the prediction of the upper limit of the constant strain rate. The strength exhibited sensitive to strain rate. A nonlinear constitutive model described the dynamic constitutive relationship effectively under different strain rates. The distribution of delamination failure significantly influenced the final failure shape of the specimens. Microscopic observations indicated that carbon fiber yarns and aramid fiber yarns exhibited different fracture modes.

A damage constitutive model for consistent description of static and fatigue behaviors of ultra-high performance concrete containing coarse aggregate

A suitable constitutive model for describing the damage behavior of ultra-high performance concrete (UHPC) under various loading histories plays a vital role in analyzing its structural performance. In this work, a consistent elastoplastic damage model is developed for UHPC containing coarse aggregate (UHPC-CA) subjected to static and fatigue loads, in which the driving and alleviating effects induced by the inclusions of coarse aggregate and steel fiber are considered. For damage evolution, based on the acting mechanism of the non-uniformity of stress wave propagation determined by the loading rate on the mechanical responses of UHPC-CA, the loading rate is skillfully integrated into the damage evolution rule to achieve the consistent description of mechanical behaviors of UHPC-CA under static and fatigue loading conditions. Regarding plasticity growth, an empirical plastic deformation model is adopted to improve the computational efficiency of structural nonlinear analysis. To verify the applicability of the proposed model, a user-defined UMAT subroutine is further developed for the subsequent numerical implementation. The comprehensive comparisons between the numerical predictions and independent experimental results at both the material level and structural level solidly demonstrate the capacity of the consistent model to capture the main features concerning the mechanical performance of UHPC-CA subjected to different loading paths.

Visualization of slit defect by scanning nonlinear airborne ultrasound source technique using compressed sensing

A high-speed, non-contact non-destructive testing method using a scanning airborne ultrasound source and a nonlinear harmonic method has been developed for plate-like structures. The testing time depends on the scanning speed and the number of measurement points. To solve this problem, we propose using compressed sensing with this method. In this paper, we investigated visualizing Lamb wave propagation using the proposed compressed sensing method. In addition, we detected a slit defect by using the Lamb wave propagation image. Our results demonstrated that compressed sensing could be used to reduce the testing time.

Finite element modeling of ultrasonic attenuation and dispersion in graded multilayer polymer composites

In this paper, finite element models of ultrasonic wave propagation in Cu–PMMA composites are established. The effects of second-phase particles, ultrasonic properties, and graded interfaces on the ultrasonic propagation behavior are investigated, and the contributions of particle-independent scattering, particle interactions, and matrix viscoelasticity to the ultrasonic attenuation are quantitatively evaluated. The results show that there is no obvious coupling between particle scattering and matrix viscoelasticity in Cu–PMMA composites, and the longitudinal wave speed has little effect with the variation of particle size, ultrasonic frequency, and graded interface, while the variation of the acoustic attenuation coefficient is related to the disorder of ultrasonic energy propagation direction. In the intra-layer uniform model, with the increase in Cu particle size and ultrasonic frequency, the scattering effect of Cu particles on the ultrasonic waves is enhanced, resulting in a significant increase in the acoustic attenuation coefficient. In the inter-layer graded multilayer model, there is a strong reverse energy propagation between the layers, causing the acoustic attenuation coefficient to increase significantly.

Elastic wave demultiplexer with frequency dependent topological valley Hall edge states

Valley Hall topological insulators (VHTIs) hold great promise for enhancing the manipulation of elastic wave propagation by their intrinsic topologically protected mechanism. Different from most of VHTIs designed by the deterministic Dirac degeneracy, a more flexible design of VHTIs is proposed by the accidental Dirac degeneracy to steer elastic wave propagation. Based on the accidental Dirac degeneracy, a kind of hexagonal phononic crystal is designed to independently control the topological phase transitions at different frequency ranges. Consequently, a two-channel topological demultiplexer is designed with the function of frequency separation for flexural waves, and its effectivity is verified by numerical simulations and experimental testing. Comparing with traditional designs of demultiplexers, the VHTIs-based demultiplexer possesses a series of advantages in robust performance, easy fabrication and low energy leakage, and sheds light on developing new generation of elastic wave devices.

Feasibility study of the application of magnetic resonance elastography (MRE) to measure oral posture related changes in stiffness of zygomaticus major muscle

PurposeDevelopment of a technique for measuring the mechanical properties of zygomaticus major (ZM) may aid advances in clinical treatments for correcting abnormal oral posture. The objective of this work was to demonstrate the feasibility of measuring the stiffness of ZM using an MR elastography technique that incorporates a custom local driver and a phase-gradient (PG) inversion. Methods2D MRE investigations were performed for 3 healthy subjects using a vibration frequency of 90 Hz to test the prediction that the stiffness of ZM would be greater in the mouth-open compared to the mouth-closed position. MRE wave images were acquired along the long axis of ZM and processed using a 2D spatial-temporal directional filter applied in the direction of wave propagation along the long axis of the muscle. Stiffness measurements were obtained by applying the PG technique to a 1D-profile drawn in the phase image of the first harmonic of the wave images and a one-tailed paired t-test was used to compare the ZM stiffness between the two mouth postures for each participant (p < 0.05). ResultsThe mean stiffness and standard deviation (SD) of ZM across the three participants in the mouth-closed and mouth-open postures was 6.75 kPa (SD 3.36 kPa) and 15.5 kPa (SD 5.15 kPa), respectively. Changes of ZM stiffness were significantly greater in the mouth-open than the mouth-closed posture (p = 0.038). ConclusionThe feasibility of using the PG MRE technique to measure stiffness changes in a small muscle such as ZM for different mouth postures has been demonstrated. Further investigations are required in a larger cohort of participants to investigate the sensitivity and reproducibility of the technique for potential clinical application as well as in health and beauty related studies.

Experimental study on decoupled charge blasting-induced crack propagation with parabolic shaped charge

In this study, a new parabolic shaped charge structure is proposed to realize bidirectional blasting in coal mines. To examine the fracture characteristics of rock mass induced by shaped charge structure under explosion load, a dynamic caustic experimental system is adopted. Further, the blasting-induced damage distribution characteristics of rock mass caused by different decoupling coefficients are investigated. Using polymethyl methacrylate (PMMA) as the experimental material, a two-dimensional model is established. AUTODYN, an explicit dynamics simulation program, is used to examine the initial crack generation process under shaped charge blasting. The stress wave propagation rule and the relationship between blasting damage and decoupling coefficient are explored. The results show that the parabolic shaped charge structure can control the explosion energy distribution and reduce the surrounding rock damage. When the decoupling coefficient is 2, the main crack length is optimum. During the crack propagation process, the growth rate and stress intensity factor of the main crack of the specimen show an oscillatory downward trend. Two stress concentration points are generated on either side near the main crack, forming the initial crack. When the decoupling coefficient is less than 2, long cracks exist in both the shaped and non-shaped charged directions, and the energy gathering effect is poor. When the decoupling coefficient is greater than 2, the energy gathering effect is obvious, but the main crack length is shorter.

Effects of dust size distribution and non-Maxwellian electrons on shock waves in a dusty plasma

We present a study of dust acoustic shock waves in a non-Maxwellian plasma with dust charge fluctuations, which are seen to cause a dissipation term in fluid model, and consequently shocks are generated. In particular, we focus on dust acoustic waves as affected by various dust size distributions. Two distinct dust size distributions—the polynomial and the power law distributions—have been used. For analytical investigation of nonlinear wave propagation in complex plasmas, a reductive perturbation approach is used to obtain the Burgers equation. A dusty plasma system with non-Maxwellian Kappa distribution is considered and it is shown that the amplitude of a shock wave, for the dust size distribution is larger than that for the mono-sized counterpart, while the shock width manifests an opposite trend. Furthermore, the shock wave speed is also affected by the dust size distributions as well as by the nature of velocity distribution function. To benchmark our findings, we apply the proper limit on the spectral index, i.e., κ→∞, and retrieve the Maxwellian results. The current findings are crucial for comprehending respective shock distributions for a plasma system exhibiting non-thermal characteristics and having dust size distributions.

Variability and dynamics of the Kuroshio and Mindanao Current during the 2010-2011 La Niña and in late 2012

Abstract Variability of two Pacific western boundary currents (WBCs)—the Kuroshio and the Mindanao Current, during the strong 2010-2011 La Niña event is investigated using ship-based hydrographic observations and moored current meter data collected off the east coasts of the Philippines. The geostrophic currents calculated using the hydrographic data show that, during the 2010-2011 La Niña winter, the Kuroshio decreased by ∼ 10 Sv (1 Sv= 106 m3 s−1) whereas the Mindanao Current increased by ∼ 5-10 Sv, relative to the normal winter in late 2012. The interannual variability based on the hydrographic data is confirmed by moored current meter measurements and satellite altimeter geostrophic currents. A coastally trapped Kelvin wave model is used to explain the interannual variability of the two WBCs during the different ENSO phases. The good comparison of the simulated sea level anomalies around the Philippines with the altimeter data suggests that the interannual variability of the WBCs is associated with Kelvin wave propagation from the Sulawesi-Sulu seas clockwise around the Philippine archipelago. We identified that the Kelvin waves are excited by downwelling equatorial Rossby waves propagating into the Indonesian seas during the La Niña. The transport anomalies of the WBCs are comparable to the total meridional transport anomalies integrated across the interior North Pacific Ocean, suggesting the importance of the WBCs in the heat charge-discharge processes of the western Pacific warm pool during ENSO events.

Iterative addition of parallel non-local effects to full wave ICRF finite element models in axisymmetric tokamak plasmas

The current response of a hot magnetized plasma to a radio-frequency wave is non-local, turning the electromagnetic wave equation into an integro-differential equation. Non-local physics gives rise to wave physics and absorption processes not observed in local media. Furthermore, non-local physics alters wave propagation and absorption properties of the plasma. In this work, an iterative method that accounts for parallel non-local effects in 2D axisymmetric tokamak plasmas is developed, implemented, and verified. The iterative method is based on the finite element method and Fourier decomposition, with the advantage that this numerical scheme can describe non-local effects while using a high-fidelity antenna and wall representation, as well as limiting memory usage. The proposed method is implemented in the existing full wave solver FEMIC and applied to a minority heating scenario in ITER to quantify how parallel non-local physics affect wave propagation and dissipation in the ion cyclotron range of frequencies (ICRF). The effects are then compared to a reduced local plane wave model, both verifying the physics implemented in the model, as well as estimating how well a local plane wave approximation performs in scenarios with high single pass damping. Finally, the new version of FEMIC is benchmarked against the ICRF code TORIC.

On the application of Krylov subspace spectral methodologies to poroacoustic shock formation in an exponential class of inhomogeneous gases

In this communication, the propagation of poroacoustic acceleration waves in a class of inhomogeneous gases whose ambient mass density varies exponentially is considered. Employing the tools of singular surface theory, we first determine the evolution of both the jump amplitudes and the locations/velocities of their associated wave-fronts, along with a variety of related analytical results. We then turn to what have become known as Krylov subspace spectral (KSS) methods to numerically simulate the evolution of the full waveforms under consideration. These simulations are not only performed quite efficiently, since KSS allows the use of “large” CFL numbers, but also quite accurately, in the sense of capturing theoretically-predicted features of the solution profiles, since KSS customizes the computation of the components (of a solution) corresponding to the different frequencies involved. The presentation concludes with a discussion of possible acoustics-related follow-on studies.

Effect of Signal Features and Model Variables on Energy-Traced Arrival Time Picking of Acoustic Signals Used for Structural Damage Detection

To monitor damage developments in structures, various structural health monitoring methods based on different principles are used. The common aspect of elastic wave-based methods is to place appropriate sensors on the structure, to detect acoustic wave propagation and to analyze these signals the sensors transformed. The arrival time of these recorded signals to the sensors is the most significant parameter used to determine critical information such as the time and location of the damage. Therefore, the accurate calculation of the arrival time affects the accuracy of the damage detection. In this study, effects of the signal-to-noise ratio (SNR), sampling frequency, length of the signal, and length of the focal window on determining the arrival time of the signals to the sensors were investigated. For this purpose, an energy-traced arrival time picking approach (Akaike Information Criterion, AIC), which is the frequently used method in the literature, has been applied to a typical acoustic signal originated from a concrete cracking. The results of the study suggest the necessity of noise elimination, the optimum level of data logging and the ratios of focal window lengths for accurate time of arrival detection in the field monitoring of the structures using acoustic methods.