Cavitation-induced shock wave behaviour in different liquids

This paper follows our earlier work where a strong high frequency pressure peak has been observed as a consequence of the formation of shock waves due to the collapse of cavitation bubbles in water, excited by an ultrasonic source at 24 kHz. We study here the effects of liquid physical properties on the shock wave characteristics by replacing water as the medium successively with ethanol, glycerol and finally a 1:1 ethanol–water solution. The pressure frequency spectra obtained in our experiments (from more than 1.5 million cavitation collapsing events) show that the expected prominent shockwave pressure peak was barely detected for ethanol and glycerol, particularly at low input powers, but was consistently observed for the 1:1 ethanol–water solution as well as in water, with a slight shift in peak frequency for the solution. We also report two distinct features of shock waves in raising the frequency peak at MHz (inherent) and contributing to the raising of sub-harmonics (periodic). Empirically constructed acoustic pressure maps revealed significantly higher overall pressure amplitudes for the ethanol–water solution than for other liquids. Furthermore, a qualitative analysis revealed that mist-like patterns are developed in ethanol–water solution leading to higher pressures.     

Weakly nonlinear focused ultrasound in viscoelastic media containing multiple bubbles

To facilitate practical medical applications such as cancer treatment utilizing focused ultrasound and bubbles, a mathematical model that can describe the soft viscoelasticity of human body, the nonlinear propagation of focused ultrasound, and the nonlinear oscillations of multiple bubbles is theoretically derived and numerically solved. The Zener viscoelastic model and Keller–Miksis bubble equation, which have been used for analyses of single or few bubbles in viscoelastic liquid, are used to model the liquid containing multiple bubbles. From the theoretical analysis based on the perturbation expansion with the multiple-scales method, the Khokhlov–Zabolotskaya–Kuznetsov (KZK) equation, which has been used as a mathematical model of weakly nonlinear propagation in single phase liquid, is extended to viscoelastic liquid containing multiple bubbles. The results show that liquid elasticity decreases the magnitudes of the nonlinearity, dissipation, and dispersion of ultrasound and increases the phase velocity of the ultrasound and linear natural frequency of the bubble oscillation. From the numerical calculation of resultant KZK equation, the spatial distribution of the liquid pressure fluctuation for the focused ultrasound is obtained for cases in which the liquid is water or liver tissue. In addition, frequency analysis is carried out using the fast Fourier transform, and the generation of higher harmonic components is compared for water and liver tissue. The elasticity suppresses the generation of higher harmonic components and promotes the remnant of the fundamental frequency components. This indicates that the elasticity of liquid suppresses shock wave formation in practical applications.     

Method of measuring one-dimensional photonic crystal period-structure-film thickness based on Bloch surface wave enhanced Goos-Hänchen shift

This paper puts forward a novel method of measuring the thin period-structure-film thickness based on the Bloch surface wave (BSW) enhanced Goos-Hänchen (GH) shift in one-dimensional photonic crystal (1DPC). The BSW phenomenon appearing in 1DPC enhances the GH shift generated in the attenuated total internal reflection structure. The GH shift is closely related to the thickness of the film which is composed of layer-structure of 1DPC. The GH shifts under multiple different incident light conditions will be obtained by varying the wavelength and angle of the measured light, and the thickness distribution of the entire structure of 1DPC is calculated by the particle swarm optimization (PSO) algorithm. The relationship between the structure of a 1DPC film composed of TiO₂ and SiO₂ layers and the GH shift, is investigated. Under the specific photonic crystal structure and incident conditions, a giant GH shift, 5.1×10³ times the wavelength of incidence, can be obtained theoretically. Simulation and calculation results show that the thickness of termination layer and periodic structure bilayer of 1DPC film with 0.1-nm resolution can be obtained by measuring the GH shifts. The exact structure of a 1DPC film is innovatively measured by the BSW-enhanced GH shift.     

Spin waves and phase transition on a magnetically frustrated square lattice with long-range interactions

We investigate the effects of long-range interactions on the spin wave spectra and the competition between magnetic phases on a frustrated square lattice with large spin S. Applying the spin wave theory and assisted with symmetry analysis, we obtain analytical expressions for spin wave spectra of competing Neel and (π, 0) stripe states of systems containing any-order long-range interactions. In the specific case of long-range interactions with power-law decay, we find surprisingly that the staggered long-range interaction suppresses quantum fluctuation and enlarges the ordered moment, especially in the Neel state, and thus extends its phase boundary to the stripe state. Our findings illustrate the rich possibilities of the roles of long-range interactions, and advocate future investigations in other magnetic systems with different structures of interactions.     

Wave reflection and refraction in triclinic crystalline media

Summary In this paper, the reflection and refraction of a plane wave at an interface between two half-spaces composed of triclinic crystalline material is considered. It is shown that due to incidence of plane wave three types of waves, namely quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated governed by the propagation condition involving the acoustic tensor. A simple procedure has been presented for the calculation of all the three phase velocities of the quasi waves. It has been established that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation. Relations are established between directions of motion and propagation, respectively. The expressions for reflection and refraction coefficients of qP, qSV and qSH waves are obtained. Numerical results of reflection and refraction coefficients are presented for different types of anisotropic media and for different types of incident waves. Graphical representations have been made for incident qP waves, and for incident qSV and qSH waves numerical data are presented in tables.The work was completed while the author was visiting the University of Kaiserslautern, Department of Geomathematics as Visiting Professor. The Author is grateful to Professor Dr. W. Freeden for providing DAAD fellowship and all the facilities for conducting research, as well as to Dr. V.Michel for various discussions about the research work and also for all kinds of help during his stay at Kaiserslautern, Germany. This award is very gratefully acknowledged.     

NORMAL EXPANSION THEORY FOR PENETRATION OF PROJECTILE AGAINST CONCRETE TARGET

Based on the equations which describe the dynamic behavior of material under high-velocity and high-pressure shock, corresponding equations at shock front whose surface is general space curve surface were established. For concrete material, a normal expansion theory was proposed by which some deceleration about time history of the projectile can be analytically given. This normal expansion theory is not only suitable for spherical and cylindrical-nose projectile, but also suitable for other general nose projectile, for example conical nose or ogive-nose. And it is not only suitable for perpendicular shock but also for oblique shock.     

机动方舱核爆炸冲击波动力响应和安定性分析

本文以核爆炸当量、高度、距离为参数，采用等超压圆给出了普通方舱和防核加固方舱结构在核冲击波下的安全域。研究了限制舱体倾斜的压力条件。本文还以冲击波压力的时间函数为激励，求解了方舱的动力响应。从试验结果推论和计算模拟两方面，分析了方舱承受冲击波的能力。     

弹性半空间二维出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性半空间自由波场时域计算的一维化有限元方法。首先利用Snell定律确定平面波沿水平方向的传播规律,在用有限元法对弹性半空间进行离散化时,竖向单元尺寸根据波动有限元模拟精度要求确定,而水平向有限元网格尺寸根据水平向波的传播规律和采用的离散时间步长确定,使得有限元离散模型中任意节点的运动可以用水平向相邻节点的运动表示,从而将二维有限元节点运动方程组化为一维的形式。求解此一维方程组,可得到弹性半空间中一列节点的运动,再根据行波的传播规律,可确定全空间自由波场。理论分析和数值算例表明,该方法具有较高的精度和良好的稳定性。     

内含运动物体的非定常流场计算

本文改造了高阶精度的MUSCL(MonotonicUpstream-centered Scheme for Conservation Laws)差分格式,对含有膛口装置和运动弹丸的火炮膛口流场进行了数值模拟,得到了由空气冲击波和射流近场结构所组成的全流场流动图象。     

用应力波方法,测定粘弹性材料的动态力学性能

本文用粘弹性材料中应力波的弥散及衰减特性的实验数据来计算材料的动态力学性能参数。用塑料杆的纵波实验,测定了材料的拉压复柔量,并探讨了用横波实验测定材料剪切复模量的方法及有关问题。由测定结果在理论分析中应用,讨论了结果的可靠性。