共轴双涡旋声束形成离轴多涡旋的定位

涡旋声束具有螺旋的相位波前,中轴线上形成声强为零的相位奇点,其所携带的轨道角动量在粒子操控领域有着广阔的应用前景。传统声涡旋只在传播轴线上形成一个拓扑荷可控的涡旋波束,这限制了声涡旋的应用灵活性。基于环形点声源阵列和相位编码技术,利用奇偶声源分别产生共轴双涡旋声束的声场叠加,在传播截面上形成了具有中心涡旋和子涡旋的离轴多涡旋声场;研究了双涡旋拓扑参数对离轴涡旋的个数、位置及拓扑荷的影响,基于声涡旋的径向声压和相位分布,确定了离轴涡旋的离轴半径,并结合声源位置推导子涡旋中心方位角的计算公式,实现离轴涡旋的精确定位。本研究突破了沿轴分布的涡旋声场只能形成单点涡旋势阱的操控局限,为利用离轴多涡旋实现多点粒子捕获提供了理论依据,促进涡旋声场在精确粒子操控和传输方面的高效应用。      

Influence of the abdominal wall on the nonlinear propagation of focused therapeutic ultrasound

This article theoretically studies the influence of inhomogeneous abdominal walls on focused therapeutic ultrasound based on the phase screen model. An inhomogeneous tissue is considered as a combination of a homogeneous medium and a phase aberration screen. Variations of acoustic parameters such as peak positive pressure, peak negative pressure, and acoustic intensity are discussed with respect to the phase screen statistics of human abdominal walls. Results indicate that the abdominal wall can result in energy loss of the sound in the focal plane. For a typical human abdominal wall with correlation length of 7.9~mm and variance of 0.36, the peak acoustic intensity radiated from a 1~MHz transmitter with a radius of 30~mm can be reduced by about 14% at the focal plane.     

Influence of Change in Inner Layer Thickness of Composite Circular Tube on Second-Harmonic Generation by Primary Circumferential Ultrasonic Guided Wave Propagation

The influence of change in inner layer thickness of a composite circular tube is investigated on second-harmonic generation(SHG) by primary circumferential ultrasonic guided wave(CUGW) propagation. Within a secondorder perturbation approximation, the nonlinear effect of primary CUGW propagation is treated as a second-order perturbation to its linear response. It is found that change in inner layer thickness of the composite circular tube will influence the efficiency of SHG by primary CUGW propagation in several aspects. In particular, with change in inner layer thickness, the phase velocity matching condition that is originally satisfied for the primary and double-frequency CUGW mode pair selected may no longer be satisfied. This will remarkably influence the efficiency of SHG by primary CUGW propagation. Theoretical analyses and numerical results show that the effect of SHG by primary CUGW propagation is very sensitive to change in inner layer thickness, and it can be used to accurately monitor a minor change in inner layer thickness of the composite circular tube.     

On the physical origin of conical bubble structure under an ultrasonic horn

The cavitation field generated by an ultrasonic horn at low frequency and high power is known to self-organize into a conical bubble structure. The physical mechanism at the origin of this bubble structure is investigated using numerical simulations and acoustic pressure measurements. The thin bubbly layer lying at horn surface is shown to act as a nonlinear thickness resonator that amplifies acoustic pressure and distorts acoustic waveform. This mechanism explains the self-stabilization of the conical bubble structure as well as the generation of shock wave and the focusing at very short distance.     

Pulling force of acoustic-vortex beams on centered elastic spheres based on the annular transducer model

To solve the difficulty of generating an ideal Bessel beam, an simplified annular transducer model is proposed to study the axial acoustic radiation force(ARF) and the corresponding negative ARF(pulling force) exerted on centered elastic spheres for acoustic-vortex(AV) beams of arbitrary orders. Based on the theory of acoustic scattering, the axial distributions of the velocity potential and the ARF for AV beams of different orders generated by the annular transducers with different physical sizes are simulated. It is proved that the pulling force can be generated by AV beams of arbitrary orders with multiple axial regions. The pulling force is more likely to exert on the sphere with a smaller k_0a(product of the wave number and the radius) for the AV beam with a bigger topological charge due to the strengthened off-axis acoustic scattering. The pulling force decreases with the increase of the axial distance for the sphere with a bigger k_0a.More pulling force areas with wider axial regions can be formed by AV beams using a bigger-sized annular transducer.The theoretical results demonstrate the feasibility of generating the pulling force along the axes of AV beams using the experimentally applicable circular array of planar transducers, and suggest application potentials for multi-position stable object manipulations in biomedical engineering.     

Attraction between pancake vortices in the crossing lattices of layered superconductors

The intervortex interaction is investigated in very anisotropic layered superconductors in tilted magnetic field. In such a case, the crossing lattice of Abrikosov vortices (AVs) and Josephson vortices (JVs) appears. The interaction between pancake vortices, forming the AVs and JVs, produces the deformation of the AV line. It is demonstrated that, as a result of this deformation, a long range attraction between AVs is induced. This phenomenon is responsible for the dense vortex chain formation. The vortex structure in the weak perpendicular magnetic field is the vortex chain phase, where only a small part of JVs is occupied by AVs.     

板厚对无障薄板声辐射特性影响的分析

分析无障薄板的声辐射特性时通常忽略板厚对格林函数的影响而采用双层势计算.本文考虑板厚引起的声辐射阻抗,采用混合势计算结构表面声压与振速,并分析板厚对声辐射参数的影响。根据交界相容性条件,采用边界积分方程分别表示平板上下表面的声压和振速,并合并同类项.进一步将结构的动力方程代入混合势形式的振速方程中,离散声压差值和板的位移为振动模态叠加的形式,获得二重积分形式的声辐射阻抗,从而求解振动模态系数,确定声辐射特性参数.以水下简支矩形板为例计算对比了声辐射参数,并讨论了其对板厚的敏感性。结果表明:板厚引起的声辐射阻抗对声辐射参数的大小影响较小,但随着频率的增加致使共振频率发生较大偏移;在相同阶数的共振频率范围内,板厚度越大,采用双层势计算的误差越大。      

Generation and Radiation of Waves at Combined Frequencies in a Warm Inhomogeneous Plasma Layer

We investigate generation and radiation of waves at combined frequencies from an arbitrary inhomogeneous, isotropic plasma layer, when electromagnetic waves are obliquely incident on it to interact with a surface wave at the plasma boundary. This interaction has been described for both S- and P-polarized waves. We consider a warm plasma layer with thickness very small as compared to the wavelengths of oscillations. It is shown that generated waves are strongly amplified, compared to cold plasma, when phase velocities of generated waves approaches the electron thermal velocity. Waves are not emitted when P-polarized waves are incident perpendicular onto the plasma layer.     

Acoustic radiation from the submerged circular cylindrical shell treated with active constrained layer damping

Based on the transfer matrix method of exploring the circular cylindrical shell treated with active constrained layer damping(i.e., ACLD), combined with the analytical solution of the Helmholtz equation for a point source, a multi-point multipole virtual source simulation method is for the first time proposed for solving the acoustic radiation problem of a submerged ACLD shell. This approach, wherein some virtual point sources are assumed to be evenly distributed on the axial line of the cylindrical shell, and the sound pressure could be written in the form of the sum of the wave functions series with the undetermined coefficients, is demonstrated to be accurate to achieve the radiation acoustic pressure of the pulsating and oscillating spheres respectively. Meanwhile, this approach is proved to be accurate to obtain the radiation acoustic pressure for a stiffened cylindrical shell. Then, the chosen number of the virtual distributed point sources and truncated number of the wave functions series are discussed to achieve the approximate radiation acoustic pressure of an ACLD cylindrical shell. Applying this method, different radiation acoustic pressures of a submerged ACLD cylindrical shell with different boundary conditions, different thickness values of viscoelastic and piezoelectric layer, different feedback gains for the piezoelectric layer and coverage of ACLD are discussed in detail. Results show that a thicker thickness and larger velocity gain for the piezoelectric layer and larger coverage of the ACLD layer can obtain a better damping effect for the whole structure in general. Whereas, laying a thicker viscoelastic layer is not always a better treatment to achieve a better acoustic characteristic.     

Liquid flows and vibration characteristics of straight-tube cylindrical shells

The influence of a moving fluid confined by a solid circular cylindrical shell on the propagation of acoustic waves generated by sources located on the circular cylindrical shell is examined. An expression for the acoustic pressure in a moving fluid is derived including azimuthal asymmetry effects in the general case, where the fluid velocity points along the cylindrical shell axis and can be written as an infinite power series expansion in the radial co-ordinate. Secondly, continuity of pressure and normal velocity at the liquid-shell interface is imposed to (a) derive a set of coupled differential equations governing the possible vibrational modes of the shell and (b) determine dispersion relations, i.e., mode propagation constants β as a function of frequency as well as changes in β values accomodated by flow. In the remaining part of the paper, phase speed changes with flow and transit-time differentials of circular cylindrical shell vibrations are discussed with special emphasis to flow measurement properties.     

Modeling the propagation of nonlinear three-dimensional acoustic beams in inhomogeneous media

A three-dimensional model of the forward propagation of nonlinear sound beams in inhomogeneous media, a generalized Khokhlov-Zabolotskaya-Kuznetsov equation, is described. The Texas time-domain code (which accounts for paraxial diffraction, nonlinearity, thermoviscous absorption, and absorption and dispersion associated with multiple relaxation processes) was extended to solve for the propagation of nonlinear beams for the case where all medium properties vary in space. The code was validated with measurements of the nonlinear acoustic field generated by a phased array transducer operating at 2.5 MHz in water. A nonuniform layer of gel was employed to create an inhomogeneous medium. There was good agreement between the code and measurements in capturing the shift in the pressure distribution of both the fundamental and second harmonic due to the gel layer. The results indicate that the numerical tool described here is appropriate for propagation of nonlinear sound beams through weakly inhomogeneous media.     

牙齿组织光热动态特性仿真与试验研究

调制激光作用牙齿组织发生散射形成光子密度波, 而由于光热效应产生热波, 基于一维介质辐射传输漫射近似方程与一维热传导方程建立了调制激光作用牙齿组织半透明混合介质的一维热波数学模型. 利用该模型仿真分析了牙齿龋损特性参数(牙釉质龋损层光吸收系数、散射系数、热扩散系数及龋损深度)对光热辐射动态响应特性的影响与规律. 利用红外探测器(HgCdTe, 2–12 μm)记录808 nm半导体激光激发牙齿组织产生的热波信号, 由锁相放大器计算热波信号的幅值与相位. 通过频率扫描试验获得了牙齿组织的光热动态响应, 利用多参数最佳统计拟合方法得到了牙齿组织特性参数. 结果表明光热辐射测量对牙齿组织不均匀性和龋损特性均具有较高敏感性与特异性.     

环形稀疏声源阵列声涡旋的优化与操控

提出和优化了环形稀疏声源阵列涡旋声场的形成方法,进行了涡旋声压,相位和振动速度的理论推导和特性分析。通过改变声源数、声源频率和观测距离等参数,对声涡旋场的模拟结果证明,声涡旋半径主要由信号波长决定,声源频率越高,声涡旋半径越小,能量越集中;涡旋声压随着声源数的增加而线性提高,随着传播距离的增加而减小;声涡旋环形相位分布的线性度随着声源数、声源波长和传播距离的增加而提高。建立了3/4/6/8声源实验系统,径向声压和相位分布以及环形相位分布的测量结果和模拟结果具有较好的一致性,同时悬吊泡沫盘的旋转证明了声涡旋角动量的传输。本研究中声涡旋场的特性分析和参数优化为粒子操控及其在生物医学中的应用提供了依据。      

Effect of time profile of laser pulse on the excitation of sound signals in a liquid with a series of laser pulses

We have studied the sound generation with high repetition rate pulsed laser. We have solved the inhomogeneous wave equation for acoustic pressure in a liquid generated by a laser, using Green’s function formalism and convolution technique. To obtain the maximum pressure of the sound waves, we found the conditions on repetition rate and on period of laser pulse of various shapes. Our analysis shows that the sound generated in a liquid with a series of laser pulses is highly affected by the time profile of the pulses besides other parameters, namely laser beam diameter, laser beam optical wavelength, repetition rate and period of laser pulse. This effect is pronounced particularly in frequency domain. We found that the noise of higher harmonics in the generated sound can be greatly removed with the proper choice of the time profile of the laser pulses. It is found that the pressure is generated around the fundamental frequency for the half-sine and rectangular pulses, with the proper choice of repetition rate and period of pulse. The application of the present analysis for underwater communication is pointed out.     

圆孔非线性声阻抗三维时域仿真与特性分析

为研究有限振幅声波作用下圆孔的非线性声学特性,提出了基于三维时域计算流体动力学(CFD)仿真的圆孔非线性声阻抗提取方法,通过求解层流方程来模拟声信号在圆孔及上下游的传播,以及采用横向周期性边界条件来考虑高穿孔率时圆孔之间相互作用的影响。研究了不同幅值声波作用下孔径、厚度和穿孔率对声阻抗的影响规律,通过对质点振速幅值、频率和板厚等组成的无量纲参量进行非线性回归分析,得到了圆孔非线性声阻抗的拟合公式,并将其转换为可考虑多频声波影响的时域模型。最后结合声阻抗时域模型和有限差分方法计算了直通穿孔管消声器在小振幅和有限振幅声波作用下的传递损失,通过与实验测量结果的比较,验证了拟合公式的准确性和实用性。     

水声工程中基于均匀圆阵列的涡旋声波性能分析*

声波作为信息和能量的载体,一直以来在水下通信中被广泛采用,但尚未解决带宽窄、速率低的问题。在光学领域和电磁波领域,轨道角动量都表征了螺旋相位结构的自然属性;通过引入轨道角动量到声学领域中,水声通信系统的传输能力以及频谱效率都得到扩展。基于换能器圆阵列产生涡旋声波进行分析和检测,研究涡旋声波波束的阵列产生方法,给出涡旋声波波束在水下传播的特性。在主轴方向,采用均匀圆阵列产生不同拓扑模式的涡旋声波波束,确定轨道角动量拓扑模式与换能器阵列之间的对应关系;为生成不同拓扑模式下的涡旋声波,研究阵列单元数目、阵列半径、传输频率等对涡旋声波的影响。通过研究发现模式数越高,涡旋声波主瓣波束角越大,主瓣峰值越小。阵列半径越大,主瓣波束角越大,而主瓣峰值则随着阵列半径的增大而减小;频率越高,主瓣波束角越小,主瓣峰值变化不大;阵列单元数对主瓣波束角无影响,但与主瓣峰值成正比关系,阵列单元数越多,主瓣峰值越大。     

Analytical and numerical calculations of optimum design frequency for focused ultrasound therapy and acoustic radiation force

Focused ultrasound therapy relies on acoustic power absorption by tissue. The stronger the absorption the higher the temperature increase is. However, strong acoustic absorption also means faster attenuation and limited penetration depth. Hence, there is a trade-off between heat generation efficacy and penetration depth. In this paper, we formulated the acoustic power absorption as a function of frequency and attenuation coefficient, and defined two figures of merit to measure the power absorption: spatial peak of the acoustic power absorption density, and the acoustic power absorbed within the focal area. Then, we derived “rule of thumb” expressions for the optimum frequencies that maximized these figures of merit given the target depth and homogeneous tissue type. We also formulated a method to calculate the optimum frequency for inhomogeneous tissue given the tissue composition for situations where the tissue structure can be assumed to be made of parallel layers of homogeneous tissue. We checked the validity of the rules using linear acoustic field simulations. For a one-dimensional array of 4 cm acoustic aperture, and for a two-dimensional array of 4 × 4 cm² acoustic aperture, we found that the power absorbed within the focal area is maximized at 0.86 MHz, and 0.79 MHz, respectively, when the target depth is 4 cm in muscle tissue. The rules on the other hand predicted the optimum frequencies for acoustic power absorption as 0.9 MHz and 0.86 MHz, respectively for the 1D and 2D array case, which are within 6% and 9% of the field simulation results. Because radiation force generated by an acoustic wave in a lossy propagation medium is approximately proportional to the acoustic power absorption, these rules can be used to maximize acoustic radiation force generated in tissue as well.     

Boundary layer effects on sound in a circular duct

Boundary layer effects on an acoustic field in a unidirectional flow with transverse shear are studied. The acoustic pressure variation in the direction normal to that of the flow is governed in the boundary layer by a second order differential equation. The problem in the boundary layer is reduced from a two point boundary value problem to a one point boundary value problem by transforming the governing equation into the Riccati equation. The Riccati equation is easily integrated with standard numerical procedures. The integration process yields the effective admittance of the wall-boundary layer combination. The acoustic field in the uniform flow is then determined for this effective admittance. Further complications imposed by the boundary layer are thus eliminated. The simplicity of the technique allows calculation of the propagation and decay constants in a circular duct over a wide range of parameters and duct modes.     

超短强激光与固体薄膜靶作用产生keV相干X射线数值模拟研究

用一维粒子模拟研究了超短强激光脉冲与两个固体薄膜靶作用产生X射线的一个方案, 特别研究了该方案中产生相对论电子层的源靶的厚度和密度分布对产生X射线的能谱、能量转换效率的影响. 数值模拟发现当产生高能电子层的源靶的厚度d与产生的X射线的波长λ/4γ_x²相当或者更小(λ 是入射激光波长)时, 才能产生准单色的X射线光谱, 否则产生的光谱有极大展宽, 且最高频率下降很快. 另外, 当薄膜靶前面存在不均匀预等离子体时, X射线光谱会明显变差. 关键词： 相对论电子层 Doppler 频移 相干X射线