求解波动方程的准粒子离散修正辛算法

针对波动方程求解,在Hamilton体系下建立了对空间离散的准粒子体系,该准粒子体系实现简单,物理意义明确;在时间离散方面,构造了一种适合高效声波模拟的修正辛格式,该格式是在常规的二阶Partitioned Runge-Kutta(PRK)基础之上构造而成,其具有三阶时间精度,从理论上分析了修正辛格式的数值稳定性和频散性能.数值结果表明,本文提出的方法在计算时间,计算精度和计算存储量等各方面性能都有相应改善。      

气动声学中声传播问题的辛算法研究

引入辛算法对气动声学中的声传播问题进行了数值研究。采用Hamilton系统描述理想气体的声波方程,时间离散采用辛可分Runge-Kutta方法,空间离散采用近似解析方法,构造声波方程的保辛格式。将辛算法和有限差分算法分别在数值频散和计算效率等方面进行了对比分析,研究结果表明:辛算法能够有效地抑制数值频散,在计算效率方面具有明显的优越性。声传播特性模拟结果表明辛算法能够准确地模拟点源声辐射、声波干涉、反射及衍射现象。     

孤立波方程的保结构算法

讨论了孤立波方程的保结构差分算法，以一些经典的孤立波方程为例，如KdV，sine-Gordon，K-P方程，给出了它们的辛和多辛结构，说明辛和多辛算法的可适用性．提出局部守恒算法和广义保结构算法的概念，它们是保结构算法的概念自然推广．还给出一种能系统构造局部守恒格式的复合方法．数值例子说明，保结构数值能很好模拟各种孤立波的演化。     

大气动力学方程的Hamilton算法

将大气动力学方程组写成正则算子方程的形式,通过引入泊松括号,并利用原方程组的无穷个不变量,深入研究其Hamilton性质,在忽略摩擦和外源强迫的情况下,证明了大气动力体系是一个Hamilton系统,从而构造出求解它的辛算法,并用数值试验检验了该算法.     

Birkhoff意义下Hénon-Heiles方程的离散变分计算

在Birkhoff意义下研究了非线性不可积Hamilton系统——Hénon-Heiles方程的离散变分计算方法,并和辛算法及Runge-Kutta方法相比较,说明在Birkhoff意义下采用离散变分算法研究非线性不可积系统的动力学行为是合理和可行的. 关键词： Hénon-Heiles方程 离散变分方法 自治Birkhoff方程     

相对论速调管中间腔与调制电子束间的非线性互作用

在强流相对论电子束驱动的相对论速调管放大器中, 由于强流和高场强的影响, 尤其是中间腔具有高Q值, 微波腔与电子束之间的非线性作用很明显, 严重影响器件性能. 根据麦克斯韦方程组以及电子在微波场作用下运动方程给出了中间腔的束-波互作用自洽方程. 从这些方程出发, 研究了调制深度和调制频率对间隙电压幅度和相位的影响. 对比常规速调管的等效电路模型, 自洽公式给出的间隙电压幅值同粒子模拟结果更接近, 尤其是较高调制深度的情况. 同时器件带宽随调制深度的增加而变宽, 这也同粒子模拟结果一致. 由此设计了一个S波段高增益相对论放大器, 在LTD (长脉冲螺旋线)加速器上开展了相应的实验工作, 实验上获得了1.1 GW的输出功率, 器件增益49 dB. 关键词： 相对论速调管 非线性互作用 自洽方程     

二维非定常Sine-Gordon方程辛算法及其孤子数值模拟

在矩形域[-a,a]×[-a,a]内对微分算子L=(ə²)/(əx²)+(ə²)/(əy²)用5点差分格式将二维非定常Sine Gordon方程离散化为一个2×799²阶非线性Hamilton系统.对该系统使用Euler中心格式,得到一个非线性方程组.对此方程组建立迭代解法并给出了这个迭代方法的收敛条件和收敛速度.Sine Gordon方程单孤子和双孤子的数值模拟试验显示该辛算法是有效的.     

辛和非辛算法求解薛定谔方程误差分析

在哈密尔顿体系下,提出气体声波传播的一种新的谐振子模型,并引入群论确定气体声波传播过程中的分子振动模式、能级简并.新模型将气动声学声传播问题与分子振动关联起来。由于发展高效的薛定谔方程的数值计算方法,有利于联系分子的性质来解释声的传播.本文从此出发,用二阶有限差分格式和生成函数法构造的二阶辛格式分别计算一维定态谐振子势场和含时谐振子势场的薛定谔方程,分析了数值解的误差以及传播能量误差.结果表明辛算法具有明显的优势.     

SRLW方程的多辛Fourier谱格式及其守恒律

通过引进正则动量，将对称正则长波方程（简称SRLW方程）转化成多辛形式的方程组，它具有多辛守恒律；介绍了空间方向满足周期边界条件的函数的Fourier谱方法；对SRLW方程的多辛方程组在空间方向利用Fourer谱方法，时间方向上应用Euler中点格式离散，得到其多辛Fourier拟谱格式；证明此格式的一些离散守恒律．用此格式模拟了SRLW方程的单个孤立波，还模拟了多个孤立波的追赶、碰撞和分离过程．     

弹性波通过一维复合材料系统的透射性质

提出了不同结构的一维弹性波复合材料系统模型,包括一维周期结构声子晶体、标准Fibonacci准周期结构声子晶体、广义Fibonacci准周期结构声子晶体以及完全无序结构的复合材料系统. 采用模式匹配理论法,数值计算了弹性波通过一维复合材料系统的透射系数. 计算结果表明,利用特殊的准周期结构声子晶体可获得比周期结构声子晶体更宽的带隙范围,准周期结构排列的复合材料系统相当于在周期结构中引入了缺陷体一样,带隙内出现了丰富的局域模式. 对弹性波/声波在复合材料系统中局域态性质的研究有助于弹性波/声波滤波器、导波器 关键词： 弹性波复合材料 局域化     

Nonlinear interaction of plane ultrasonic waves with an interface between rough surfaces in contact

A theoretical investigation of the nonlinear interaction between an acoustic plane wave and an interface formed by two rough, nonconforming surfaces in partial contact is presented. The macroscopic elastic properties of such a nonlinear interface are derived from micromechanical models accounting for the elastic interaction that is characteristic of spherical bodies in contact. These results are used to formulate set of boundary conditions for the acoustic field, which are to be enforced at the imperfect interface. The scattering problem is solved for plane wave incidence by using a simple perturbation approach and the harmonic balance method. Sample results are presented for arbitrary wave polarization and angle of incidence. The relative magnitude of the nonlinear signals and their potential use toward the nondestructive evaluation of imperfect interfaces are assessed. In particular, attention is drawn to the enhanced nonlinear response of an interface insonified by a shear vertical wave in the neighborhood of the longitudinal critical angle. The motivation for this investigation is provided by the need to develop nondestructive methods to detect and localize small, partially closed cracks in metals with coarse microstructures.     

Bubble interaction dynamics in Lagrangian and Hamiltonian mechanics

Two models of interacting bubble dynamics are presented, a coupled system of second-order differential equations based on Lagrangian mechanics, and a first-order system based on Hamiltonian mechanics. Both account for pulsation and translation of an arbitrary number of spherical bubbles. For large numbers of interacting bubbles, numerical solution of the Hamiltonian equations provides greater stability. The presence of external acoustic sources is taken into account explicitly in the derivation of both sets of equations. In addition to the acoustic pressure and its gradient, it is found that the particle velocity associated with external sources appears in the dynamical equations.     

Slow two-level systems in point contacts

A great variety of experiments, like heat release measurements, acoustic measurements, and transport measurements on mesoscopic samples, have proved that two-level systems (TLSs) play a crucial role in the low-temperature thermal and electric properties of disordered systems. This paper is aimed at reviewing the role of slow TLSs in point contacts. First the theory of point contacts is summarized, concentrating on the discussion of different point-contact models, and on the different regimes of electron flow in the contact, mainly focusing on the ballistic and diffusive limit. The Boltzmann equation is solved in both regimes, and the position dependence of the electrical potential is determined. Then the scattering processes in point contacts are investigated, particularly concentrating on the scattering on slow TLSs. If the electron-assisted transitions between the two states are negligible the electron–two-level system interaction can be treated with a simplified Hamiltonian. The scattering on such slow TLSs causes non-linearity in the current–voltage characteristics of the point contact, which can be determined using Fermi's Golden Rule. These calculations are presented showing both the contribution of elastic and inelastic scattering, and including the dependence on the position of the TLS, and on the effect of high-frequency irradiation. These results are used to discuss the differences between these slow TLSs and the fast centres which may be described by the two-channel Kondo model. The available experimental results are analysed, distinguishing between the effects due to the different types of TLSs.     

Unitary transformation method for solving generalized Jaynes-Cummings models

Two fully quantized generalized Jaynes-Cummings models for the interaction of a two-level atom with radiation field are treated, one involving intensity dependent coupling and the other involving multiphoton interaction between the field and the atom. The unitary transformation method presented here not only solves the time dependent problem but also allows a determination of the eigensolutions of the interacting Hamiltonian at the same time.     

极性晶体中与形变势相互作用的表面极化子

有不少的极性晶体,电子与表面光学(SO)声子耦合强,但与表面声学(SA)声子耦合弱.研究电子与SO声子耦合强,与SA声子耦合弱的极性晶体中与形变势相互作用的表面极化子的性质.采用改进了的线性组合算符和微扰法导出了极性晶体中与形变势相互作用的表面极化子的有效哈密顿量.在计及电子在反冲效应中发射和吸收不同波矢的声子之间的相互作用时,讨论对表面极化子的有效哈密顿量、有效质量和有效相互作用势的影响.     

驻波声场中单分散细颗粒的相互作用特性

利用外加声场促进悬浮在气相中的细颗粒发生相互作用,进而引起颗粒的碰撞和凝并,使得颗粒平均粒径增大、数目浓度降低,是控制细颗粒排放的重要技术途径.为探究驻波声场中单分散细颗粒的相互作用,建立包含曳力、重力、声尾流效应的颗粒相互作用模型,采用四阶经典龙格-库塔算法和二阶隐式亚当斯插值算法对模型进行求解.将数值模拟得到的颗粒声波夹带速度和相互作用过程与相应的解析解和实验结果进行对比,验证模型的准确性.进而研究颗粒初始条件和直径对相互作用特性的影响.结果表明,初始时刻颗粒中心连线越接近声波波动方向、颗粒位置越接近波腹点,颗粒间的声尾流效应就越强,颗粒发生碰撞所需要的时间就越短.研究还发现,颗粒直径对颗粒相互作用的影响取决于初始时刻颗粒中心连线偏离声波波动方向的程度.当偏离较小时,颗粒直径越大,颗粒发生碰撞所需要的时间越短;当偏离很大时,直径较小的颗粒能够发生碰撞,而直径较大的颗粒则无法发生碰撞.     

Influence of acoustic phonons on dephasing of free excitons in quantum wells

A theory of the dephasing rate of quasi-2D free excitons due to acoustic phonon interaction at low exciton densities is presented. Both deformation potential and piezoelectric couplings are considered for the exciton–phonon interaction in quantum wells. Using the derived interaction Hamiltonian obtained recently by us, exciton linewidth and dephasing rate are calculated as a function of the exciton density, exciton temperature, exciton momentum and lattice temperature.     

Unbiased sampling of globular lattice proteins in three dimensions

We present a Monte Carlo method that allows efficient and unbiased sampling of Hamiltonian walks on a cubic lattice. Such walks are self-avoiding and visit each lattice site exactly once. They are often used as simple models of globular proteins, upon adding suitable local interactions. Our algorithm can easily be equipped with such interactions, but we study here mainly the flexible homopolymer case where each conformation is generated with uniform probability. We argue that the algorithm is ergodic and has dynamical exponent z=0. We then use it to study polymers of size up to 64(3)=262 144 monomers. Results are presented for the effective interaction between end points, and the interaction with the boundaries of the system.     

Interaction of a Surface Acoustic Wave with a Two-dimensional Electron Gas

When a surface acoustic wave （SAW） propagates on the surface of a GaAs semiconductor, coupling between electrons in the two-dimensional electron gas beneath the interface and the elastic host crystal through piezoelectric interaction will attenuate the SAW. The coupling coemcient is ~alculated for the SAW propagating along an arbitrary direction. It is found that the coupling strength is strongly dependent on the propagating direction. When the SAW propagates along the [011] direction, the coupling becomes quite weak.