虚阶Bessel函数的计算

给出了计算虚阶Bessel函数及其导数的完整算法。数值结果表明,计算精度可达10位以上有效数字。此外,所给出的逆推法可推广到复阶数的情形。     

非球形粒子的散射特性分析

本文首先利用T矩阵方法计算了复折射率吸收指数和折射指数变化时的椭球粒子和Chebyshev粒子在不同等效尺度参数下的光散射特性,并与等效的球形粒子的光散射结果进行了比较;然后分析以上两种类型非球形粒子散射特性之间的关系.结果表明:椭球粒子和Chebyshev粒子的散射特性与等效球形粒子的散射特性存在着差别,粒子的形状越偏离球形,这种差别就越大;复折射率折射指数的变化对非球形散射效率因子的影响要比吸收指数的影响更大一些;当等效尺度参数相同时,椭球粒子与等效球形粒子的散射效率因子的差别要远远大于Chebyshev粒子与等效球形粒子散射效率因子之间的差别.     

基于递归T矩阵的离散随机散射体散射特性研究

本文根据电磁场矢量球波函数多极点展开原理及矢量叠加定理提出了递归T矩阵算法的矢量形式,并且基于矢量递归T矩阵算法建立了多散射球模拟离散随机散射体散射的三维电磁散射模型.通过计算不同尺寸、随机分布散射球的散射以及分析散射球间的高阶散射效应,结果表明:矢量递归T矩阵算法具有很高的计算精度,算法中包含多散射体间的高阶散射效应,因此能够精确计算多散射体总的散射效应.本文所建模型可应用于土壤湿度探测工程中评估地表下掩埋离散随机散射体散射对雷达回波信号产生的影响.     

雨滴对电磁波的散射特性研究

由于重力和向上气流的作用,高速摄影显示较大的雨滴在降落过程中表现为一底部扁平、顶部圆滑的非球形对称体,已有的模型中BC模型(由Beard和Chuang提出)被认为能准确描述雨滴的形状。本文利用沿降落方向旋转对称的扁椭球体和均匀球体模拟降落中的雨滴,扁椭球体的短轴(旋转轴)与BC模型的相等,长轴通过与BC模型截面积等效法确定;利用T矩阵方法研究了水平取向的旋转扁椭球体和BC雨滴模型的散射特性,利用Mie理论计算了均匀球的散射特性,比较了半径分别为2 mm和3 mm大小的三种雨滴模型对94 GHz电磁波斜入射时的单次散射相矩阵,以及对频率为30 GHz~100 GHz电磁波的单次散射反照率,最后讨论了半径为0.5 mm~3.5 mm雨滴对35 GHz和94 GHz的后向散射函数。计算结果表明,利用均匀球模拟雨滴将带来较大误差,而利用截面积等效的旋转扁椭球体模型可以较好地模拟雨滴的散射特性。     

取向比对椭球气溶胶粒子散射特性的影响

利用T矩阵和离散坐标法研究了取向比对椭球粒子散射特性的影响, 计算了小尺度范围内椭球粒子的散射特征参量, 包括消光效率因子、不对称因子、单次散射反照率、散射相矩阵及双向反射函数(BRDF). 结果表明, 椭球粒子的散射特性与取向比密切相关, 粒子取向比会影响散射参量的振荡频率和振幅, 与球形粒子散射参量的相对差异也呈周期振荡趋势. 研究还发现, 某些特殊粒子尺寸的散射参量与粒子取向比基本无关. 在多次散射条件下, 分析不同取向比粒子群的BRDF随反射角和光学厚度的变化特性. 结果显示: 不同取向比粒子群的BRDF随反射角的变化趋势基本一致, 球形粒子群比非球形粒子群的BRDF曲线波动振幅更大; 球形-非球形粒子的BRDF相对差异随光学厚度和取向比的增大而减小, 随入射角的增大而增大.     

基于Fourier-Bessel级数的Bessel型超声场二次谐波近场特性研究

将Fourier-Bessel级数引入KZK方程的求解，用于计算黏滞媒质中零阶Bessel型超声场的二次谐波声场，得到其级数形式的解析解，并由此得出二次谐波声场在近场分布的一个新结论.设声源表面声压分布为J₀(α₀r)，则二次谐波声压在近场的径向分布服从J²₀(α₀r)函数规律.这一结论合理解释了相关的实验结果，表明二次谐波声场在近场和远场有不同的径向分布，从而解决了非线性Bessel型超声场二次谐波的近场分布问题.研究还发现二次谐波声场具有类似基波声场的有限衍射特性.给出了一个数值计算和仿真实例.     

太赫兹波在雾中的多重散射特性

太赫兹(THz)波作为微波和毫米波的延伸,它所提供的通信带宽和容量远大于毫米波。在随机介质中传播时,不但会发生时域和空域的形变,介质中的粒子还会对入射波发生散射,这些都会使得脉冲信号发生衰减。根据Mie理论与随机离散分布粒子的波传播与散射理论,计算了THz波信号入射下雾滴粒子的消光系数。结合雾滴粒子谱分布,得到了雾媒质的平均体系散射特性,采用蒙特卡罗法得到了平流雾对THz信号的多重散射特性,计算了THz波段信号对平流雾的透过率与反射率,分析了THz波段信号的前向、后向散射特性随散射角的分布。结果表明,低能见度大气环境中,雾对THz波产生的吸收和衰减不容忽视。相关研究结果对THz在大气传输、通信等方面的应用具有重要意义。     

包容性非球形气溶胶粒子的散射特性及其对激光传输的影响研究

针对局部地域大气气溶胶的生成、演化的特点,在分析典型对称结构的包容性非球形气溶胶粒子的形态学特征基础上,利用T矩阵法进行了椭球形粒子含非球形粒子的包容性粒子光散射特性计算和分析.以1 064nm波长激光传输为例分别计算了两种典型分布状况的包容性非球形粒子的散射相函数等散射特性,并针对稠密气溶胶大气路径上激光传输的多散射现象,利用蒙特卡罗方法进行了典型大气状况下包容性粒子多次散射对激光信号大气传输中脉冲展宽的应用分析,在传输3km和5km距离时多次散射效应造成激光脉冲峰值半宽分别为3μs和4μs左右.指出气溶胶的基本散射特性与粒子的形态学特征、尺度因子、化学成分及颗粒密度等有密切的关系,其散射特性对激光信号的大气传输具有明显的影响.     

电子碰撞原子弹性散射的扭曲波解

使用Ｓｌａｔｅｒ型轨道参数描述束缚电子，连续电子与靶之间的相互作用，采用静势近似，充分利用Ｒ－矩阵程序包中的现有资源，首次在“二区”近似图象下，建立了我们自己的计算扭曲库仑波函数的计算方法，并编制了标准程序。     

新型压缩感知计算模型分析三维电大目标电磁散射特性

为提高基于压缩感知技术的矩量法在三维电大目标双站电磁散射问题中的计算效率和稳定性,提出新的稀疏、测量和重构方法,构建一种新型压缩感知计算模型.不同于基于欠定方程的传统的压缩感知计算模型,新型计算模型首先采用按行均匀抽取阻抗矩阵的方法构造测量矩阵以获得稳定的计算结果;然后,基于Foldy-Lax方程生成多阶特征基函数并作为稀疏基对感应电流进行稀疏转换;再依据少数低阶特征基函数足以近似表征感应电流的先验条件,将恢复算法简化为最小二乘法;最后,将矩阵方程转换为一个超定系统并采用最小二乘法解出电流系数.与传统的计算模型相比,新型计算模型不仅可以获得更加稳定的精确解,还可以显著提高电大目标双站散射问题的求解效率和计算精度.数值仿真结果证明了新方法的可行性和高效性.     

Acoustic scattering of a Bessel vortex beam by a rigid fixed spheroid

Partial-wave series representation of the acoustic scattering field of high-order Bessel vortex beams by rigid oblate and prolate spheroids using the modal matching method is developed. The method, which is applicable to slightly elongated objects at low-to-moderate frequencies, requires solving a system of linear equations which depends on the partial-wave index n$n$ and the order of the Bessel vortex beam m$m$ using truncated partial-wave series expansions (PWSEs), and satisfying the Neumann boundary condition for a rigid immovable surface in the least-squares sense. This original semi-analytical approach developed for Bessel vortex beams is demonstrated for finite oblate and prolate spheroids, where the mathematical functions describing the spheroidal geometry are written in a form involving single angular (polar) integrals that are numerically computed. The transverse (θ=π/2$\theta =\pi /2$) and 3D scattering directivity patterns are evaluated in the far-field for both prolate and oblate spheroids, with particular emphasis on the aspect ratio (i.e., the ratio of the major axis over the minor axis of the spheroid) not exceeding 3:1, the half-cone angle β$\beta$ and order m$m$ of the Bessel vortex beam, as well as the dimensionless size parameter kr₀$k{r}_0$. Periodic oscillations in the magnitude plots of the far-field scattering form function are observed, which result from the interference of the reflected waves with the circumferential (Franz’) waves circumnavigating the surface of the spheroid in the surrounding fluid. Moreover, the 3D directivity patterns illustrate the far-field scattering from the spheroid, that vanishes in the forward (θ=0$\theta =0$) and backward (θ=π$\theta =\pi$) directions. Particular applications in underwater acoustics and scattering, acoustic levitation and the detection of submerged elongated objects using Bessel vortex waves to name a few, would benefit from the results of the present investigation.     

Instantaneous axial force of a high-order Bessel vortex beam of acoustic waves incident upon a rigid movable sphere

The present investigation examines the instantaneous force resulting from the interaction of an acoustical high-order Bessel vortex beam (HOBVB) with a rigid sphere. The rigid sphere case is important in fluid dynamics applications because it perfectly simulates the interaction of instantaneous sound waves in a reduced gravity environment with a levitated spherical liquid soft drop in air. Here, a closed-form solution for the instantaneous force involving the total pressure field as well as the Bessel beam parameters is obtained for the case of progressive, stationary and quasi-stationary waves. Instantaneous force examples for progressive waves are computed for both a fixed and a movable rigid sphere. The results show how the instantaneous force per unit cross-sectional surface and unit pressure varies versus the dimensionless frequency ka (k is the wave number in the fluid medium and a is the sphere’s radius), the half-cone angle β and the order m of the HOBVB. It is demonstrated here that the instantaneous force is determined only for (m, n) = (0, 1) (where n is the partial-wave number), and vanishes for m > 0 because of symmetry. In addition, the instantaneous force and normalized amplitude velocity results are computed and compared with those of a rigid immovable (fixed) sphere. It is shown that they differ significantly for ka values below 5. The proposed analysis may be of interest in the analysis of instantaneous forces on spherical particles for particle manipulation, filtering, trapping and drug delivery. The presented solutions may also serve as a method for comparison to other solutions obtained by strictly numerical or asymptotic approaches.     

Quasi Bessel beam by Billet's N-split lens

This study utilizes the focal property of a classical Billet's split lens to create more focal points by splitting the lens. This approach distributes the focal points circularly on the focal plane. This study explores the characteristics of beam propagation and analytically derives the asymptotic characteristics of beam propagation based on the stationary phase approximation and the moment-free Filon-type method. Results show that the unique Billet's N-split lens can generate a quasi Bessel beam if the number of splitting N is large enough, e.g., N ≧ 24. This study also explores the diffraction efficiency of corresponding quasi Bessel beam and the influence of aperture size. The potential advantage of proposed split lens approach is that, unlike the classical means of annular aperture, this simple lens approach allows a much larger throughput in creating the Bessel beam and hence the Bessel beam could have more optical energy.     

轴棱锥对无衍射光束的线聚焦特性

基于汉克尔波理论和衍射积分理论详细分析了轴棱锥对无衍射光束的线聚焦特性, 提出了一种产生周期性局域空心光束的新方法, 即无衍射贝塞尔光束经过轴棱锥聚焦后产生具有塔尔博特效应的局域空心光束. 数值模拟了无衍射贝塞尔光束照射轴棱锥后, 沿传输距离变化的光强分布及一个周期内光强的演变和局域空心光束的形成过程. 设计实验系统, 由He-Ne激光经过一套光学系统后透过轴棱锥, 产生近似无衍射贝塞尔光束, 再由第二个轴棱锥对产生的无衍射贝塞尔光束进行线聚焦, 在第二个轴棱锥后面由显微镜观测到周期性局域空心光束, 并用CCD照相机拍摄了两个周期内的光斑图, 实验结果和理论分析相符合. 研究结果可用于多层面微粒的操控, 对周期性局域空心光束在光学微操控领域的应用具有重要的指导意义. 关键词： 贝塞尔光束 轴棱锥 线聚焦 局域空心光束     

非稳腔主动式直接获取纳秒近似无衍射贝塞尔绿光

通过对基于轴棱锥的贝塞尔谐振腔和贝塞尔-高斯谐振腔的研究,设计了一台腔内倍频Nd:YAG纳秒近似无衍射贝塞尔绿光激光器. 非稳激光器谐振腔由轴棱锥和凸面镜组成. 实验采用单灯脉冲抽运激光增益介质Nd:YAG晶体,倍频晶体选用KTP. 当抽运电压为350 V时,由主动式直接产生纳秒近似零阶贝塞尔绿光,其脉冲宽度为55.1 ns,波长为532 nm,线宽为0.8 nm,近似无衍射零阶贝塞尔绿光的中心光斑直径为192 μm. 利用Fresnel-Kirchhoff 衍射积分和Fox-Li迭代法,通过数值计算得出 关键词： 主动式 纳秒贝塞尔光束 腔内倍频 轴棱锥     

The generation of Bessel beams at millimetre-wave frequencies by use of an axicon

An axicon is used to generate a Bessel beam at 90 GHz in the millimetre-wave region of the spectrum. The Bessel beam has a central intensity maximum of approximately 4 mm in diameter that is maintained over a propagation distance greater than 60 mm.     

Transverse (lateral) instantaneous force of an acoustical first-order Bessel vortex beam centered on a rigid sphere

In a recent report [F.G. Mitri, Z.E.A. Fellah, Ultrasonics 51 (2011) 719-724], it has been found that the instantaneous axial force (i.e. acting along the axis of wave propagation) of a Bessel acoustic beam centered on a sphere is only determined for the fundamental order (i.e. m = 0) but vanishes when the beam is of vortex type (i.e. m > 0, where m is the order (or helicity) of the beam). It has also been recognized that for circularly symmetric beams (such as Bessel beams of integer order), the transverse (lateral) instantaneous force should vanish as required by symmetry. Nevertheless, in this commentary, the present analysis unexpectedly reveals the existence of a transverse instantaneous force on a rigid sphere centered on the axis of a Bessel vortex beam of unit magnitude order (i.e. |m| = 1) not reported in [F.G. Mitri, Z.E.A. Fellah, Ultrasonics 51 (2011) 719-724]. The presence of the transverse instantaneous force components of a first-order Bessel vortex beam results from mathematical anti-symmetry in the surface integrals, but vanishes for the fundamental (m = 0) and higher-order Bessel (vortex) beams (i.e. |m| > 1). Here, closed-form solutions for the instantaneous force components are obtained and examples for the transverse components for progressive waves are computed for a fixed and a movable rigid sphere. The results show that only the dipole (n = 1) mode in the scattering contributes to the instantaneous force components, as well as how the transverse instantaneous force per unit cross-sectional surface varies versus the dimensionless frequency ka (k is the wave number in the fluid medium and a is the sphere’s radius), and the half-cone angle β of the beam. Moreover, the velocity of the movable sphere is evaluated based on the concept of mechanical impedance. The proposed analysis may be of interest in the analysis of transverse instantaneous forces on spherical particles for particle manipulation and rotation in drug delivery and other biomedical or industrial applications.     

双焦透镜对Bessel光束传输的影响

研究非轴对称象散光学元件双焦透镜对无衍射Bemel光束的聚焦特性.基于空间域中的广义惠更斯.菲涅耳衍射理论衍射积分,用4×4阶矩阵来描述非轴对称光学系统,导出Bessel光通过双焦透镜更为普遍和广义的光强分布表达式,并对其光强进行数值模拟.分析双焦透镜在x,y面上的焦距取不同值时对Bessel光的影响.导出的广义光强分布表达式包括了普通透镜、双焦透镜和柱透镜的所有情况.研究发现,当象散越大,图形畸变程度越大,导致得到的光学bottle beam的光束质量下降,bottle壁变形,囚禁微观粒子的能力也随之下降.     

Evanescent Bessel beam generation through filtering highly focused cylindrical vector beams with a defect mode one-dimensional photonic crystal

We propose a simple setup for generating evanescent Bessel beams using the defect mode of one-dimensional (1D) photonic crystal. The angular selectivity provided by the defect mode mimics the role of an axicon for Bessel beam generation. When an azimuthally polarized beam is strongly focused onto a 1D defect mode photonic crystal interface, an evanescent Bessel beam of the first-order is produced, while an evanescent Bessel beam of the zeroth-order will be created under a radially polarized beam illumination. Switching between a donut shape and a solid focal distribution can be easily realized by controlling the polarization of the illumination. Such a versatile evanescent Bessel beam generation may find potential applications in optical trapping.     

Interaction of a high-order Bessel beam with a submerged spherical ultrasound contrast agent shell - Scattering theory

Background and objective

Acoustic scattering properties of ultrasound contrast agents are useful in extending existing or developing new techniques for biomedical imaging applications. A useful first step in this direction is to investigate the acoustic scattering of a new class of acoustic beams, known as helicoidal high-order Bessel beams, to improve the understanding of their scattering characteristics by an ultrasound contrast agent, which at present is very limited.

Method

The transverse acoustic scattering of a commercially available albuminoidal ultrasound contrast agent shell filled with air or a denser gas such as perfluoropropane and placed in a helicoidal Bessel beam of any order is examined numerically. The shell is assumed to possess an outer radius a = 3.5 microns and a thickness of ∼105 nm. Moduli of the total and resonance transverse acoustic scattering form functions are numerically evaluated in the bandwidth 0 < ka? 3, which corresponds to a frequency bandwidth of 0-205 MHz that covers a wide range of applications for imaging with contrast agents. Particular attention is paid to the shell’s material, the content of its interior hollow region and the fluid surrounding its exterior. The contrast agent shell is assumed to be immersed in an ideal compressible fluid so the viscous corrections are not considered. Analytical equations are derived and numerical calculations of the total and resonance form functions are performed with particular emphasis on the effect of varying the half-cone angle, the order of the helicoidal Bessel beam as well as the fluid that fills the interior hollow space.

Results and conclusion

It is shown that shell wave resonance modes can be excited on an encapsulated micro-bubble. The forward and backscattering vanish for a helicoidal high-order Bessel beam. Additionally, the fluid filling the inner core affects the shell’s response significantly. Moreover, there is no monopole contribution to the axial scattering of a helicoidal Bessel beam of order m ? 1 so that the dynamics of contrast agents would be significantly altered. The main finding of the present theory is the suppression or enhancement for a particular resonance that may be used to advantage in imaging with ultrasound contrast agents for clinical applications.