Time-dependent propagation of high energy laser beams through the atmosphere

The computation of time-dependent three-space-dimensional laser beam propagation is described. The methods are applicable to the propagation of high energy laser beams through the atmosphere in the presence of a horizontal wind and turbulence for most situations of interest. Possible cases are propagation of cw beams through stagnation zones, multi-pulse propagation, including the self-consistent treatment of pulse self-blooming, and propagation involving transonic slewing. The solution of the Maxwell wave equation in Fresnel approximation is obtained by means of a discrete Fourier transform method, which, surprisingly, gives excellent results for diffraction problems. The latter provide a stringent test for the accuracy of any solution method. Considerable use is also made of discrete Fourier transform methods in solving the hydrodynamic equations. The treatment of turbulence is based on the generation of random phase screens at each calculation step along the propagation path. In a time-dependent calculation the random phase screens can be either made to move with the wind at a given propagation position or generated anew for each successive time. This work was done under contract to the Advanced Research Projects Agency of the Department of Defense, the Army Missile Command, Huntsville, Alabama, and the U.S. Energy Research and Development Administration.     

Studies of propagation characteristics of elliptical Gaussian beams through circular apertured optical systems

An approximate analytical propagation formula of the elliptical Gaussian beams through a paraxial ABCD optical system with a circular aperture is obtained on the basis of the expansion of a hard-edged aperture into a finite sum of complex Gaussian functions. And the parameter characteristics of the truncated beams including the beam propagation factor and the kurtosis parameter are also studied in detail by using the second-order-moments method.     

Three-dimensional FDTD method for optical pulse propagation analysis in microstructured optical fibers

We propose a three-dimensional (3D) finite-difference time-domain (FDTD) method to analyze the pulse propagation characteristics in microstructured optical fibers (MOFs). The computation domain size is greatly reduced by adopting the technique of moving problem space. The propagating pulse is virtually held in the buffer cell of the problem space as simulation continues. This method is capable to investigate the temporal evolution of the propagating pulse. Spectral information can be obtained by Fourier analysis. As an example, the influence of the kerr nonlinearity on the optical pulse propagation in a Lorentz dispersive MOF is demonstrated. The model is also used to simulate the nonlinear interactions between the pump spectral broadening and third harmonic generations in a highly nonlinear fused silica nanowire with good agreement with the generalized nonlinear envelop equation (GNEE) model.     

强非局域介质中1+1维高斯型损耗空间双孤子的相互作用

研究了1+1维高斯型双光束在含小损耗的强非局域非线性介质中的传输特性。通过对该介质中光束传输遵循的非局域非线性薛定谔方程进行近似简化,得到了含小损耗强非局域非线性介质中1+1维高斯型双光束传输模型。在此基础上运用解析的方法研究了双光束传输的演化规律,得到了准双孤子解。经过进一步分析发现,在传输过程中两光束中心的轨迹为艾里函数;两光束会准周期性地碰撞、分离;随着传输距离的增大,两光束中心之间的最大距离会越来越大。另一方面,当损耗逐渐增大时,两光束的碰撞空间周期将变短,同时两光束中心之间的最大距离也越来越大。     

Elastic waves excited by a plane source on the surface of a multilayered medium

Elastic waves excited by a plane piezoelectric source with an arbitrary shape on the surface of a multilayered medium have been studied for the first time in this paper. On the basis of Abzo-zena [Geophys. J. R. Astron. Soc. 58, 91-105 (1979)] and Menke [Geophys. J. R. Astron. Soc. 59, 315-323 (1979)], the propagator matrix for the elastic wave field in multilayered medium is extended from two- to three-dimensional (3D) space. 3D elastic wave field is investigated and the displacement-stress response for the boundary conditions is obtained. The propagation of elastic wave in multilayered media is analyzed in 3D space in the frequency domain. The P-SV and SH modes corresponding to the poles are studied. The excitation and propagation of the modes are analyzed further. It is found that the propagation velocities of the P-SV and SH modes do not depend on the propagation azimuth theta in the plane parallel to the free surface of the multilayered medium while the displacement amplitudes are strongly dependent on the azimuth theta. The directional distribution functions of the modes are independent of the medium parameters and the modes and dependent on the shape and excitation fashion of the source. Finally, as an example, the displacement fields of the P-SV and SH modes excited by a rectangle source are analyzed. The displacement representation and numerical results of the directivity distribution functions for the P-SV and SH modes are obtained.     

High precision analysis of finlines on semiconductor substrate

In this study, the effect of the metallization thickness in finlines on semiconcuctor substrate is researched. The propagation parameters are computed to measure the inluence of the metallization. The theory and numerical results are presented to the propagation constant and characteristic impedance of the bilateral and unilateral finlines. The full wave analysis of the transverse transmission line — TTL method is used to determine the electromagnetic fields of the structure in Fourier transform domain — FTD. Applying the suitable boundary conditions, the moment method and expanding the fields in a set of suitable bases functions, a homogeneous matrix system is obtained and the propagation constant is computed. The characteristic impedance is obtained using the relation of the voltage in the slot and the transmitted power by the structure.Computational programs are developed to obtain numerical results to the propagation parameters composed by the propagation constant and characteristic impedance.This work received financial support by CAPES and CNPq.     

Propagation characteristics of laser beams with amplitude modulations and phase fluctuations in apertured fractional Fourier transforming system

The apertured fractional Fourier transforming system is introduced and applied to treat the propagation of Gaussian beams with amplitude modulations and phase fluctuations. Based on the treatment that a rectangular function can be expanded into an approximate sum of complex Gaussian functions with finite numbers, the analytical expressions for the mutual intensity distribution of laser beams with amplitude modulations and phase fluctuations passing through the apertured fractional Fourier transforming system are obtained. Some numerical simulations are illustrated for their propagation properties.     

On dynamic plane elasticity problems of 2D quasicrystals

In this Letter, the dynamic plane elasticity problems of 2D quasicrystals is considered. By use of the Fourier transform and matrix transformations the system is reduced to uncoupled ordinary differential equations. Fourier images of Green's functions for dynamic plane elasticity problems of 2D dodecagonal, pentagonal and decagonal quasicrystals are obtained explicitly by the suggested method.     

一种捕捉多向刚性爆轰波的激波定位随机投影方法

提出一种激波定位随机投影方法,利用当地随机投影,引入激波指示器并提出恰当的判定当地投影区域的规则,使爆轰波多向推进的问题得以解决。数值实验表明,激波定位随机投影方法不仅适用于爆轰波的多向推进,还能高效求解刚性反应方程,易于推广到三维情况。     

光束通过硬边光阑光学系统变换的改进算法

 使用傅里叶级数与有限个复高斯函数之和的乘积作为硬边光阑的窗口函数的近似式，提出了含硬边光阑近轴ABCD光学系统对光束变换的改进算法。以高斯光束的变换为例，推导出了一个表示为初等函数之和的近似解析传输公式。给出了数值计算结果，并与直接积分Collins公式和使用Wen方法的结果作了比较。研究表明，改进算法能兼顾计算精度和机时，并能用于远场和近场。而使用Wen方法时，在近场（小于0.12倍菲涅尔距离）有明显误差。     

拉曼增益对孤子传输特性的影响

利用考虑拉曼增益效应的非线性薛定谔方程, 在忽略光纤损耗的情况下, 采用基于MATLAB的分步傅里叶数值算法, 得出线性算符和非线性算符具体的表达式, 分步作用于光孤子脉冲传输方程, 仿真模拟了光孤子在光纤中传输时的演变. 与不考虑拉曼增益的光孤子在光纤中传输相对比, 探析了拉曼增益对孤子传输特性的影响.拉曼增益会破坏孤子的传输周期, 导致孤子在光纤中传输时快速衰减, 并且影响程度和输入孤子的脉冲峰值功率大小有关, 拉曼增益对基态孤子和高阶孤子的影响也不相同. 关键词： 拉曼增益 孤子 对称分步傅里叶法 非线性薛定谔方程     

Solutions of the full- and quasi-vector wave equations based on H-field for optical waveguides by using mapped Galerkin method

The mapped Galerkin method in solving the full-vector and quasi-vector wave equations in terms of transverse magnetic fields (H-formulation) for optical waveguides with step-index profiles is described. By transforming the whole x-y space onto a unit square and using two-dimensional Fourier series expansion, the modal distributions and propagation constants for optical waveguides are obtained in the absence of boundary truncation. Results for step-index circular fiber, buried rectangular waveguide, and optical rib waveguide are presented. Solutions are good agreed with exact solutions and numerical results by using vector nonlinear iterative method, Fourier operator transform method, and vector beam propagation method.     

Laser amplification with a twist: traveling-wave propagation and gain functions from throughout the cochlea

Except at the handful of sites explored by the inverse method, the characteristics-indeed, the very existence-of traveling-wave amplification in the mammalian cochlea remain largely unknown. Uncertainties are especially pronounced in the apex, where mechanical and electrical measurements lack the independent controls necessary for assessing damage to the preparation. At a functional level, the form and amplification of cochlear traveling waves are described by quantities known as propagation and gain functions. A method for deriving propagation and gain functions from basilar-membrane mechanical transfer functions is presented and validated by response reconstruction. Empirical propagation and gain functions from locations throughout the cochlea are obtained in mechanically undamaged preparations by applying the method to published estimates of near-threshold basilar membrane responses derived from Wiener-kernel (chinchilla) and zwuis analysis (cat) of auditory-nerve responses to broadband stimuli. The properties of these functions, and their variation along the length of the cochlea, are described. In both species, and at all locations examined, the gain functions reveal a region of positive power gain basal to the wave peak. The results establish the existence of traveling-wave amplification throughout the cochlea, including the apex. The derived propagation and gain functions resemble those characteristic of an active optical medium but rotated by 90 degrees in the complex plane. Rotation of the propagation and gain functions enables the mammalian cochlea to operate as a wideband, hydromechanical laser analyzer.     

A fast-scanning Fourier transform 2D IR interferometer

Ultrafast two-dimensional infrared spectroscopy (2D IR) allows for the characterization of vibrational couplings and chemical dynamics. The fastest method of acquiring a Fourier transform 2D IR data set involves spectrally dispersing the signal field onto an infrared array detector. However, use of this method carries disadvantages, including the high cost of IR arrays and the decrease of signal intensity due to dispersion. As an alternative, we demonstrate a readily implemented full time-domain 2D IR detection method in which data from a pulsed laser source is rapidly acquired by scanning an interferometer delay at constant velocity. The stage's position is determined with high accuracy on a shot-to-shot basis by quadrature detection of HeNe tracer interference fringes.     

等束腰宽度超短脉冲光束传输的脉冲修正方法

利用角谱分析和傅里叶变换的方法,得到一种描述几个周期的等束腰宽度脉冲光束传输的脉冲修正方法。以准单色光束传输的结果为出发点,通过对准单色光束的解进行泰勒级数展开,得到了一种相对简单的修正方法,可以精确的描述具有任意时间波形和横向光束分布的不短于一个周期的超短脉冲光束的传输行为。给出等束腰宽度超短脉冲的近似解,具体研究高斯脉冲光束的传输特性,分析几种不同的频谱对脉冲光束传输行为的影响。     

Spectral solvers for spherical elliptic problems

A new family of direct spectral solvers for the 3D Helmholtz equation in a spherical gap and inside a sphere for nonaxisymmetric problems is presented. A variational formulation (no collocation) is adopted, based on the Fourier expansion and the associated Legendre functions to represent the angular dependence over the sphere and using basis functions generated by Legendre or Jacobi polynomials to represent the radial structure of the solution. In the present method, boundary conditions on the polar axis and at the sphere center are not required and never mentioned, by construction. The spectral solution of the vector Dirichlet problem is also considered, by employing a transformation that uncouples the spherical components of the Fourier modes and that is implemented here for the first time. The condition numbers of the matrices involved in the scalar solvers are computed and the spectral convergence of all the proposed solution algorithms is verified by numerical tests.     

基于OpenMP的飞秒强激光在空气中传输并行计算

轴对称超短强激光在空气中的传输可以用2D+1维非线性薛定谔方程来描述,该方程一般可用FCN方法进行求解,即在时间上应用快速傅里叶变换方法、在空间横截面上采用Crank-Nicholson差分法。但由于计算非常耗时,传统串行的FCN方法只能计算强激光在百m量级距离上传输。基于OpenMP设计了求解2D+1维非线性薛定谔方程的并行方法及其程序。数值模拟结果证明了程序的正确性和较高的并行效率。当线程数为15时并行加速比为12。此并行方法可应用于模拟长距离的超短超强激光的传输。     

Reflection and transmission coefficients for guided waves reflected by defects in viscoelastic material plates

The paper presents a Fourier transform-based signal processing procedure for quantifying the reflection and transmission coefficients and mode conversion of guided waves diffracted by defects in plates made of viscoelastic materials. The case of the S(0) Lamb wave mode incident on a notch in a Perspex plate is considered. The procedure is applied to numerical data produced by a finite element code that simulates the propagation of attenuated guided modes and their diffraction by the notch, including mode conversion. Its validity and precision are checked by the way of the energy balance computation and by comparison with results obtained using an orthogonality relation-based processing method.     

A variational mode expansion mode solver

A variational approach for the semivectorial modal analysis of dielectric waveguides with arbitrary piecewise constant rectangular 2D cross-sections is developed. It is based on a representation of a mode profile as a superposition of modes of constituting slab waveguides times some unknown continuous coefficient functions, defined on the entire coordinate axis. The propagation constant and the lateral functions are found from a variational principle. It appears that this method with one or two modes in the expansion preserves the computational efficiency of the “standard” effective index method while providing more accurate estimates for propagation constants, and well-defined continuous approximations for mode profiles. By including a larger number of suitable trial fields, the present approach can also serve as a technique for rigorous semivectorial mode analysis.     

Propagation of elliptical Gaussian beam passing through apertured paraxial optical systems

The propagation of elliptical Gaussian beam passing through paraxial optical systems with aperture is investigated analytically by using tensor method. The approximate formula for propagation of elliptical Gaussian beam through hard apertured optical systems is derived based on the fact that the circ function can be expanded into a finite sum of complex Gaussian functions. The derived formula provides a convenient tool for treating the propagation and transformation of elliptical Gaussian beam through apertured optical systems. As an application example, the propagation properties of elliptical Gaussian beam through apertured fractional Fourier systems are discussed.