Numerical simulation of parametric sound generation and its application to length-limited sound beam

In this study we propose a simulation model for predicting the nonlinear sound propagation of ultrasound beams over a distance of a few hundred wavelengths, and we estimate the beam profile of a parametric array. Using the finite-difference time-domain method based on the Yee algorithm with operator splitting, axisymmetric nonlinear propagation was simulated on the basis of equations for a compressible viscous fluid. The simulation of harmonic generation agreed with the solutions of the Khokhlov–Zabolotskaya–Kuznetsov equation around the sound axis except near the sound source. As an application of the model, we estimated the profiles of length-limited parametric sound beams, which are generated by a pair of parametric sound sources with controlled amplitudes and phases. The simulation indicated a sound beam with a narrow truncated array length and a width of about one-quarter to half that of regular a parametric beam. This result confirms that the control of sound source conditions changes the shape of a parametric beam and can be used to form a torch like low-frequency sound beam.     

Model equation for strongly focused finite-amplitude sound beams

A model equation that describes the propagation of sound beams in a fluid is developed using the oblate spheroidal coordinate system. This spheroidal beam equation (SBE) is a parabolic equation and has a specific application to a theoretical prediction on focused, high-frequency beams from a circular aperture. The aperture angle does not have to be small. The theoretical background is basically along the same analytical lines as the composite method (CM) reported previously [B. Ystad and J. Berntsen, Acustica 82, 698-706 (1996)]. Numerical examples are displayed for the amplitudes of sound pressure along and across the beam axis when sinusoidal waves are radiated from the source with uniform amplitude distribution. The primitive approach to linear field analysis is readily extended to the case where harmonic generation in finite-amplitude sound beams becomes significant due to the inherent nonlinearity of the medium. The theory provides the propagation and beam pattern profiles that differ from the CM solution for each harmonic component.     

非局域非线性介质中光束传输的拉盖尔-高斯变分解

光束在非局域非线性介质中的传输过程由非局域非线性薛定谔方程描述.1+2D非局域非线性薛定谔方程可以转化为圆柱坐标系下的变分问题.通过展开介质响应函数并合理假设试探解求解变分方程,得到光束在强非局域非线性介质中的拉盖尔-高斯解.满足一定条件时,拉盖尔-高斯光束将形成光孤子或退化为高斯光束. 关键词： 非局域非线性介质 强非局域性 变分法 拉盖尔-高斯光束     

Nonlinear evolution equations for surface plasmons for nano-focusing at a Kerr/metallic interface and tapered waveguide

Maxwell's equations for a metallic and nonlinear Kerr interface waveguide at the nanoscale can be approximated to a (1+1) D Nonlinear Schrodinger type model equation (NLSE) with appropriate assumptions and approximations. Theoretically, without losses or perturbations spatial plasmon solitons profiles are easily produced. However, with losses, the amplitude or beam profile is no longer stationary and adiabatic parameters have to be considered to understand propagation. For this model, adiabatic parameters are calculated considering losses resulting in linear differential coupled integral equations with constant definite integral coefficients not dependent on the transverse and longitudinal coordinates. Furthermore, by considering another configuration, a waveguide that is an M–NL–M (metal–nonlinear Kerr–metal) that tapers, the tapering can balance the loss experienced at a non-tapered metal/nonlinear Kerr interface causing attenuation of the beam profile, so these spatial plasmon solitons can be produced. In this paper taking into consideration the (1+1)D NLSE model for a tapered waveguide, we derive a one soliton solution based on He's Semi-Inverse Variational Principle (HPV).     

Sound propagation in two-dimensional waveguide with circular wavefront

Sound propagation through a waveguide is generally modeled by the Webster horn equation which assumes a planar pressure wavefront. However, most of the sources are non-planar in nature. In this work, a 1-D model is derived for sound propagation through a 2-D waveguide with circular wavefront. The model is derived from the 2-D Helmholtz equation using the weighted residual method. The model assumes a uniform pressure across the angular coordinate at a given radial distance. A 2-D finite element model is used to validate the results for different waveguide geometries and it shows good agreement.     

电子束空间极限电流的非线性理论研究

从电子束一维稳态传输的电势非线性Poisson方程出发,推导了圆柱波导内实心束和环形束空间极限电流与电子入射电势的依赖关系,给出了数值求解方法和解的特征,分析比较了数值计算与现有解析公式及粒子模拟的结果.考虑电子的横向运动,对数值方法进行了二维修正,计算结果与粒子模拟结果具有很好的一致性.说明利用数值方法计算电子电势非线性微分方程能够得到更精确的电子束在圆柱波导内的空间极限电流;另外,对其他形状的波导,尤其是难以得到解析式的情况,根据实际几何结构设置边界条件,数值方法可以方便地给出束流传输特性,对设计新型结构的高功率微波器件提供理论指导. 关键词： 相对论电子束 圆柱波导 空间极限电流 束流传输     

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.     

The application of Gaussian beam. Formalism to optical propagation in nonlinear media

The nonlinear propagation problem of optical beams in media with an intensity dependent indices of refraction is solved in the case of a fundamental gaussian beam. This solution is made possible by assuming everywhere along the propagation path a quadratic index profile. The quadratic constant is obtained self-consistently from the second term in a Taylor expansion of the local gaussian intensity. The solution constitutes an extension of the well known ABCD formalism to nonlinear propagation.     

新型反馈式轴向同轴虚阴极振荡器

提出一种利用低磁场导引环形电子束在同轴波导内轴向激励虚阴极振荡的新型虚阴极振荡器.利用反馈式结构,使虚阴极在由反馈环、阳极箔和同轴波导构成的高品质因子准谐振腔内形成.作用腔内环形电子束激励同轴波导TM₀₁模式,在相对低品质因子的同轴提取区转化为TEM主模输出.用25维KARAT粒子模拟软件研究得到束波功率转换效率12%,输出微波平均27 GW,中心频率38 GHz. 关键词： 环形电子束 轴向导引磁场 同轴波导 反馈环     

新型反馈式轴向同轴虚阴极振荡器

从电子束一维稳态传输的电势非线性Poisson方程出发,推导了圆柱波导内实心束和环形束空间极限电流与电子入射电势的依赖关系,给出了数值求解方法和解的特征,分析比较了数值计算与现有解析公式及粒子模拟的结果.考虑电子的横向运动,对数值方法进行了二维修正,计算结果与粒子模拟结果具有很好的一致性.说明利用数值方法计算电子电势非线性微分方程能够得到更精确的电子束在圆柱波导内的空间极限电流;另外,对其他形状的波导,尤其是难以得到解析式的情况,根据实际几何结构设置边界条件,数值方法可以方便地给出束流传输特性,对设计新型结构的高功率微波器件提供理论指导.     

Harmonic propagation of finite amplitude sound beams: experimental determination of the nonlinearity parameter B/A

We propose a method to determine the nonlinearity parameter B/A of a liquid from the spatial evolution of harmonic components. We describe an analytical model, in the parabolic and quasi-linear approximations, that predicts the continuous finite amplitude sound beam propagation radiated by a plane piston source. This model takes into consideration attenuation, diffraction and nonlinear effects. The fundamental and second harmonic ultrasonic fields are expressed as the superposition of Gaussian beams. Axial propagation curves are then compared with those obtained by direct numerical solution of the transformed beam equation using the finite difference method, and with experimental results. Accurate measurements of pressure levels for the nonlinearly generated harmonics in water are performed along and across the propagation axis for different pressure values delivered at the piston surface. Experimental results, for water and ethanol, are in agreement with those of our model, which allows us to obtain the expected value of the nonlinearity parameter B/A.     

The contribution of a lateral wave in simulating low-frequency sound fields in an irregular waveguide with a liquid bottom

A further development of a previously proposed approach to calculating the sound field in an arbitrarily irregular ocean is presented. The approach is based on solving the first-order causal mode equations, which are equivalent to the boundary-value problem for acoustic wave equations in terms of the cross-section method. For the mode functions depending on the horizontal coordinate, additional terms are introduced in the cross-section equations to allow for the multilayer structure of the medium. A numerical solution to the causal equations is sought using the fundamental matrix equation. For the modes of the discrete spectrum and two fixed low frequencies, calculations are performed for an irregular two-layer waveguide model with fluid sediments, which is close to the actual conditions of low-frequency sound propagation in the coastal zone of the oceanic shelf. The calculated propagation loss curves are used as an example for comparison with results that can be obtained for the given waveguide model with the use of adiabatic and one-way propagation approximations.     

Simulation of low-frequency sound propagation in an irregular shallow-water waveguide with a fluid bottom

A further development of a previously proposed approach to calculating the sound field in an arbitrarily irregular ocean is presented. The approach is based on solving the first-order causal mode equations, which are equivalent to the boundary-value problem for acoustic wave equations in terms of the cross-section method. For the mode functions depending on the horizontal coordinate, additional terms are introduced in the cross-section equations to allow for the multilayer structure of the medium. A numerical solution to the causal equations is sought using the fundamental matrix equation. For the modes of the discrete spectrum and two fixed low frequencies, calculations are performed for an irregular two-layer waveguide model with fluid sediments, which is close to the actual conditions of low-frequency sound propagation in the coastal zone of the oceanic shelf. The calculated propagation loss curves are used as an example for comparison with results that can be obtained for the given waveguide model with the use of adiabatic and one-way propagation approximations.     

Propagation properties of a sinh-Gaussian beam in a Kerr medium

The effect of the Kerr nonlinearity on a sinh-Gaussian beam is investigated by using the nonlinear Schrödinger equation (NLS). Based on the method of moments, the evolution of the sinh-Gaussian beam width in the root-mean-square (RMS) sense is analytically described. Numerical simulations indicate that the central parts of each lobe of the sinh-Gaussian beam give rise to initially radial compression and the beam profile redistribution occurs during propagation even though the RMS beam width broadens. The partial collapse of the center parts of each lobe of the beam will appear below the threshold for a global collapse as predicted by the method of moments. The sinh-Gaussian beams eventually convert into sin-Gaussian type beams in Kerr media with low and moderate initial power.     

一类新平顶光束的束传播因子

 基于迭加原理，引入了利用双曲余弦-高斯光束的相干迭加而得的新光束，结果表明这类迭加光场的场强分布在适当参数条件下可具有平顶结构，并且分析了这类迭加所得的平顶光束的束传播因子与束参数的关系。     

Time-reversed sound beams of finite amplitude.

Numerical simulations based on the nonlinear parabolic wave equation are used to investigate time reversal of sound beams radiated by unfocused and focused sources. Emphasis is placed on nonlinear propagation distortion in the time-reversed beam, and specifically its effect on field reconstruction. Distortion of this kind, due to amplification during time reversal, has been observed in recent experiments [A. P. Brysev et al., Acoust. Phys. 44, 641-650 (1998)]. Effects of diffraction introduced by time-reversal mirrors with finite apertures are also considered. It is shown that even in the presence of shock formation, the ability of time reversal to retarget most of the energy on the source or focal region of the incident beam is quite robust.     

复宗量Laguerre-Gauss光束在强非局域非线性介质中的传输

研究得到了偏离束腰入射的复宗量Laguerre-Gauss光束在强非局域非线性介质中传输的解析表达式,并且得到了其二阶矩束宽的解析解.通过例子研究了偏离束腰入射的复宗量Laguerre-Gauss光束在强非局域非线性介质中传输性质.结果表明：非(0, m)模的复宗量Laguerre-Gauss光束的光束形状随着传输而发生改变,并以Δz=πz_c为周期做周期性演化.而(0,m)模复宗量Laguerre-Gauss光束在演化过程中则形状保持不变,仅改变光束宽度;不论功率多大,在偏离束腰入射条件下总是表现为呼吸子;只有当其为束腰入射,并且入射功率等于临界功率时才能形成孤子. 关键词： 强非局域非线性 复宗量Laguerre-Gauss光束 二阶矩束宽 空间光孤子     

Beam interaction in one-dimensional inhomogeneous Kerr nonlinear arrays

The interactions between two parallel and co-polarized beams in weakly coupled cubic focusing nonlinear waveguide arrays with transverse inhomogeneous modulation of the refractive index were investigated by means of beam propagation method. Results show that the in-phase beams attract each other and coalesce when the inhomogeneous parameter is smaller than a critical value. For the out-of-phase beams, oscillations exist at low power level, and the larger inhomogeneous parameter, the larger period and amplitude of oscillation. At a high power level, solitonlike propagation of weak coupling occurs. The inhomogeneity of waveguide arrays provide a flexible way to control the interactions of beams, and may find potential applications in all-optical systems.     

强非局域非线性介质中的涡旋光孤子

满足强非局域条件时，光束在非局域非线性介质的传输过程由Snyder-Mitchell模型描述.在旋转柱坐标系下求解了Snyder-Mitchell模型，得到涡旋光孤子的自相似旋转解析解.结果表明，涡旋光孤子解在径向是惠特克函数与幂函数的乘积，光束的光强呈环形分布，光束绕光束中心旋转. 关键词： 非局域非线性介质 强非局域性 涡旋光孤子 惠特克函数     

Coupled tapering/uptapering of soliton pairs in nonlinear media

Phenomenon of coupled tapering/uptapring of two mutually incoherent beams coaxially co-propagating in a nonlinear medium with small gain or loss has been investigated in this paper using standard parabolic equation approach (PEA) and the results are compared with the results obtained by Beam Propagation Method (BPM), i.e., by direct simulations of the underlying Nonlinear Schrödinger Equation (NLSE). The PEA results are shown to be in excellent agreement with the BPM results. It is seen that both beams of the pair induce uptapering in each other in presence of losses and tapering in presence of gain. When the medium offers gain to the first beam and losses to the other, both beams taper. When the medium offers gain/absorption to only one of the two beams, the beam undergoes self-tapering/self-uptapering and induces a taperd/uptaperd waveguide. The other beam (for which the medium is lossless) uptapers/tapers due to the taperd/uptaperd waveguide created by the first beam.