五阶克尔非线性光纤放大器中光脉冲的传输特性

以变系数高阶非线性薛定谔方程为理论模型,对光脉冲在具有五阶非线性克尔效应光纤放大器中的传输特性进行了研究.在加入振幅微扰、相位扰动、噪声微扰的情况下,数值模拟了亮孤波、灰孤波和黑孤波在光纤系统中的传输稳定性,并且讨论了孤子间的相互作用.数值模拟结果表明:在有限的微扰下,三种光脉冲具有良好的传输稳定性.     

基于洛伦兹函数微扰的磁光光纤布喇格光栅线性透射谱和非线性双稳特性研究

基于非线性耦合模理论,利用逆向递推龙格-库塔法,数值研究了在磁光耦合系数中引入洛伦兹函数微扰后对磁光光纤布喇格光栅线性透射谱和非线性双稳特性的影响.研究结果表明:给磁光耦合系数引入洛伦兹微扰能在线性透射谱阻带中打开线宽极窄的一个透射窗口,微扰宽度和微扰中心位置可以影响透射窗口的位置、宽度以及峰值大小;当微扰宽度和微扰中...     

含非线性微扰项的二阶动力学系统的一阶近似守恒量的一种新求法

从一阶近似守恒量的性质出发,把受微扰系统视为未受微扰系统与微扰项的迭加,提出一种分三步求得一阶近似守恒量的新方法:先选择合适的方法求得未受微扰系统的守恒量I0,再考虑微扰项对守恒量I0的影响,最后利用一阶近似守恒量的性质求得一阶近似守恒量.用该方法研究了一实际的受非线性微扰作用的两自由度动力学系统,得到4个稳定的一阶近似守恒量.用坐标变换法和微扰法得到系统一阶近似解的表达式,并讨论4种特殊情况下的一阶近似解.     

基于洛伦兹函数微扰的磁光光纤布喇格光栅线性透射谱和非线性双稳特性研究

基于非线性耦合模理论,利用逆向递推龙格-库塔法,数值研究了在磁光耦合系数中引入洛伦兹函数微扰后对磁光光纤布喇格光栅线性透射谱和非线性双稳特性的影响.研究结果表明:给磁光耦合系数引入洛伦兹微扰能在线性透射谱阻带中打开线宽极窄的一个透射窗口,微扰宽度和微扰中心位置可以影响透射窗口的位置、宽度以及峰值大小;当微扰宽度和微扰中心位置发生变化时,光栅的双稳特性表现出明显的差异,合理地选择微扰参量可以实现对其双稳特性的优化.     

简并微扰法求解拉比模型

运用简并微扰的方法得到拉比模型在二阶微扰下的能量和—阶微扰下的波函数.所得结果在当前的实验参数范围内与数值模拟符合非常好.这一方法对求解拉比模型的系统能级具有指导意义,有益于拉比模型的进一步应用.     

连续波微扰下飞秒光孤子的传输特性研究

给出了描述连续波扰动下飞秒光脉冲在光纤系统中传输的微扰高阶非线性薛定谔方程,通过矩法和微扰理论分析了飞秒亮孤子与连续波相互作用的特性,并利用龙格-库塔积分和分步傅里叶方法进行了数值模拟.结果表明,连续波微扰对光纤孤子通讯系统是十分有害的,在实际通讯过程中应当尽量避免连续波的渗入.     

一类滞后相对转动动力学方程的分岔特性及其解析近似解

建立了一类含Davidenkov滞后环的非线性相对转动动力学方程.分别分析了该非线性相对转动自治方程和微外扰下非自治方程的分岔特性,并采用KBM法求解了滞后环指数n=2时该非线性相对转动方程在周期激励下的解析近似解.通过数值仿真,得到了几种分岔结构及外扰下全局分岔图,同时将数值解与本文KBM法求解结果进行比较,证明本文求解结果有较高的精度,为研究这一类滞后相对转动系统提供了理论参考依据. 关键词： 相对转动 滞后环 分岔 KBM法     

用参数微扰法求解非线性谐振子问题

用参数微扰法求解了一维非线性谐振子及二维非线性耦合谐振子的一级近似本征能量.计算表明,参数微扰法的基态和激发态结果比相应的无参数微扰法的一级近似甚至多级近似结果精度有显著提高.     

微扰的耦合非线性薛定谔方程的近似求解

将直接微扰方法应用于可积的含修正项的非线性薛定谔方程，通过近似解与精确解的比较确定了直接微扰方法的可靠性.继而，将该方法应用于微扰的耦合非线性薛定谔方程，并获得了该微扰方程的可靠的近似解. 关键词： 直接微扰方法 微扰 耦合非线性薛定谔方程 近似解     

内开槽螺纹回旋行波管线性理论研究

基于耦合波理论，运用阻抗微扰法得到了内开槽螺纹波导的色散方程及耦合条件，并从相对论电子运动方程入手得到了自洽非线性方程组，通过一阶线性近似导出了注-波互作用线性关系及色散方程.运用数值计算及3维粒子模拟软件冷腔分析对该波导冷腔色散特性及注-波互作用色散特性进行了研究. 关键词： 内开槽螺纹波导 阻抗微扰 线性理论 自洽非线性方程组     

基于空域分解的非线性有限声束快速计算

研究非线性有限声束的一种快速数值计算算法.理论研究表明非线性有限声束存在多种耦合:谐波之间的全耦合,各场点沿轴向的递推耦合以及沿径向的局部耦合,因此可以通过声场的径向空间区域分割提高计算效率,采用多线程实现并行计算.对非线性高斯聚焦波的计算结果表明,当声场计算规模与分割线程数合理匹配时,算法能够显著提高计算速度同时保证计算精度,计算结果与理论分析相符.     

Analytical and numerical prediction of harmonic sound power in the inlet of aero-engines with emphasis on transonic rotation speeds

Tone noise radiated through the inlet of a turbofan is mainly due to rotor-stator interactions at subsonic regimes (approach flight), and to the shock waves attached to each blade at supersonic helical tip speeds (takeoff). The axial compressor of a helicopter turboshaft engine is transonic as well and can be studied like turbofans at takeoff. The objective of the paper is to predict the sound power at the inlet radiating into the free field, with a focus on transonic conditions because sound levels are much higher. Direct numerical computation of tone acoustic power is based on a RANS (Reynolds averaged Navier–Stokes) solver followed by an integration of acoustic intensity over specified inlet cross-sections, derived from Cantrell and Hart equations (valid in irrotational flows). In transonic regimes, sound power decreases along the intake because of nonlinear propagation, which must be discriminated from numerical dissipation. This is one of the reasons why an analytical approach is also suggested. It is based on three steps: (i) appraisal of the initial pressure jump of the shock waves; (ii) 2D nonlinear propagation model of Morfey and Fisher; (iii) calculation of the sound power of the 3D ducted acoustic field. In this model, all the blades are assumed to be identical such that only the blade passing frequency and its harmonics are predicted (like in the present numerical simulations). However, transfer from blade passing frequency to multiple pure tones can be evaluated in a fourth step through a statistical analysis of irregularities between blades. Interest of the analytical method is to provide a good estimate of nonlinear acoustic propagation in the upstream duct while being easy and fast to compute. The various methods are applied to two turbofan models, respectively in approach (subsonic) and takeoff (transonic) conditions, and to a Turbomeca turboshaft engine (transonic case). The analytical method in transonic appears to be quite reliable by comparison with the numerical solution and with available experimental data.     

非线性声波方程二次谐波高阶近似解及实验验证

非线性声波方程的传统二次谐波二阶近似解不够精确,对测量材料的非线性系数造成较大误差.利用摄动法展开非线性声波方程得到一系列非齐次偏微分方程,根据低阶解的性质拟定高阶特解的形式,通过符号计算求出高次谐波特解,对所有二次谐波成分求和最终获得二次谐波的高阶近似解.在水中开展非线性声学实验验证高阶近似解,结果表明:二次谐波相对...     

The effect of turbulence damping on acoustic wave propagation in tubes

The attenuation of sound due to the interaction between a low Mach number turbulent boundary layer and acoustic waves can be significant at low frequencies or in narrow tubes. In a recent publication by the present authors the acoustics of charge air coolers for passenger cars has been identified as an interesting application where turbulence attenuation can be of importance. Favourable low-frequency damping has been observed that could be used for control of the in-duct sound that is created by the engine gas exchange process. Analytical frequency-dependent models for the eddy viscosity that controls the momentum and thermal boundary layers are available but are restricted to thin acoustic boundary layers. For cases with cross-sections of a few millimetres a model based on thin acoustic boundary layers will not be applicable in the frequency range of interest.In the present paper a frequency-dependent axis-symmetric numerical model for interaction between turbulence and acoustic waves is proposed. A finite element scheme is used to formulate the time harmonic linearized convective equations for conservation of mass, momentum and energy into one coupled system of equations. The turbulence is introduced with a linear model for the eddy viscosity that is added to the shear viscosity. The proposed model is validated by comparison with experimental data from the literature.     

Extended characteristics method in nonlinear geometric acoustics

An asymptotic “extended characteristics method” is developed for solving nonlinear Riemann-type wave equations as applied to calculating the ray pattern of intense spatially modulated waves in weakly inhomogeneous media. The method makes it possible to avoid the singularity related to the foci of the initial wavefront, calculate the displacement in foci caused by the inhomogeneity of the medium, and thus calculate the ray pattern and intensity of the acoustic field. The beauty of the method is an exact nonlinear transfer equation for the field along the ray and the construction of its general solution for an arbitrary form of inhomogeneity. It is shown that the method is applicable to calculating the spatial structure of intense focused waves and wave beams outside the focal region in a nonlinear geometric acoustics approximation.     

不确定声场分析的二阶区间摄动有限元法

针对一阶区间摄动有限元法在声场参数不确定程度增大时误差过大的缺陷,在二阶Taylor展开的基础上推导了声学二阶区间摄动有限元法,并将其应用于区间不确定声场的声压响应分析。该方法先对声学区间有限元方程的声压响应向量进行二阶Taylor展开,获取声压响应的二阶近似响应向量;再根据二次函数极值定理获得声压响应向量的上下界。二维管道声场与轿车声腔模型的数值分析算例表明,与一阶区间摄动有限元法相比,二阶区间摄动有限元法有效提高了计算精度。因此二阶区间摄动有限元适合不确定度更大的区间不确定声场声压响应分析,具有良好的工程应用前景。      

Computer models in room acoustics: The ray tracing method and the auralization algorithms

Computer algorithms are described for constructing virtual acoustic models of various rooms that should satisfy some specific sound quality criteria. The algorithms are based on the ray tracing method, which, in the general case, allows calculation of the amplitude of an acoustic ray that survived multiple reflections from arbitrary curved surfaces. As a result, calculations of room acoustics are reduced to tracing the trajectories of all the acoustic rays in the course of their propagation with multiple reflections from reflecting surfaces to the point of their complete decay. For this approach to be used, the following physical properties of a room should be known: the geometry of the reflecting surfaces, the absorption and diffusion coefficients on each of these surfaces, and the decay law for rays propagating in air. The proposed models allow for the solution of the important problem of architectural acoustics called the auralization problem, i.e., to predict how any given audio segment will sound in any given hall on the basis of computer simulation alone, without any full-scale testing in specific halls.     

谐振管内非线性驻波的有限体积数值算法

为了扩展谐振管内非线性驻波在工程中的应用, 以及克服现有数值计算方法仅局限于求解直圆柱形和指数形谐振管内非线性驻波的问题. 根据变截面的非稳态可压缩热黏性流体Navier-Stokes方程和空间守恒方程, 并基于求解压力速度耦合方程的半隐式算法和交错网格技术, 构建一种能够计算任意形状轴对称谐振管受活塞驱动时内部非线性驻波的有限体积算法. 分别对圆柱形、指数形和圆锥形谐振管内的非线性驻波进行仿真计算. 通过与现有试验结果以及数值仿真结果的对比, 验证了该方法的正确性.并获得除驻波声压之外的另外一些新的物理结果, 包括速度、密度、温度的瞬时变化.在直圆柱形谐振管内产生冲击声压波, 速度波形中出现钉状结构.而在指数形和圆锥形谐振管内产生高声压幅值的驻波, 没有出现冲击波, 速度波形中均未发现钉状结构. 计算结果表明谐振管内非线性驻波的物理属性与谐振管形状之间有密切关系.     

Nonlinear pulsed ultrasound beams radiated by rectangular focused diagnostic transducers

A numerical model for simulating nonlinear pulsed beams radiated by rectangular focused transducers, which are typical of diagnostic ultrasound systems, is presented. The model is based on a KZK-type nonlinear evolution equation generalized to an arbitrary frequency-dependent absorption. The method of fractional steps with an operator-splitting procedure is employed in the combined frequency-time domain algorithm. The diffraction is described using the implicit backward finite-difference scheme and the alternate direction implicit method. An analytic solution in the time domain is employed for the nonlinearity operator. The absorption and dispersion of the sound speed are also described using an analytic solution but in the frequency domain. Numerical solutions are obtained for the nonlinear acoustic field in a homogeneous tissue-like medium obeying a linear frequency law of absorption and in a thermoviscous fluid with a quadratic frequency law of absorption. The model is applied to study the spatial distributions of the fundamental and second harmonics for a typical diagnostic ultrasound source. The nonlinear distortion of pulses and their spectra due to the propagation in tissues are presented. A better understanding of nonlinear propagation in tissue may lead to improvements in nonlinear imaging and in specific tissue harmonic imaging. Published in Russian in Akusticheskiĭ Zhurnal, 2006, Vol. 52, No. 4, pp. 560–570. This article was translated by the authors.     

Lie Symmetry and Nonlinear Instability in Computation of KdV Solitons

The soliton calculation method put forward by Zabusky and Kruskal has played an important role in the development of soliton theory, however numerous numerical results show that even though the parameters satisfy the linear stability condition, nonlinear instability will also occur. We notice an exception in the numerical calculation of soliton, gain the linear stability condition of the second order Leap-frog scheme constructed by Zabusky and Kruskal, and then draw the perturbed equation with the finite difference method. Also, we solve the symmetry group of the KdV equation with the knowledge of the invariance of Lie symmetry group and then discuss whether the perturbed equation and the conservation law keep the corresponding symmetry. The conservation law of KdV equation satisfies the scaling transformation, while the perturbed equation does not satisfy the Galilean invariance condition and the scaling invariance condition. It is demonstrated that the numerical simulation destroy some physical characteristics of the original KdV equation. The nonlinear instability in the calculation of solitons is related to the breaking of symmetry.