截断各向异性等离子体的修正各向异性完全匹配层吸收边界

基于磁化等离子体本构方程,提出了一种截断各向异性色散介质的修正的各向异性完全匹配层(M-UPML)吸收边界算法。通过等效相对介电常数,将UPML推广到截断各向异性等离子体介质的情形,并推导了其时域有限差分方法(FDTD)迭代式。用该方法计算了半空间磁化与非磁化等离子体的反射系数,计算结果与解析解相一致,表明该吸收边界具有良好的吸收效果。     

截断各向异性等离子体的修正各向异性完全匹配层吸收边界

 基于磁化等离子体本构方程,提出了一种截断各向异性色散介质的修正的各向异性完全匹配层(M-UPML)吸收边界算法。通过等效相对介电常数,将UPML推广到截断各向异性等离子体介质的情形,并推导了其时域有限差分方法(FDTD)迭代式。用该方法计算了半空间磁化与非磁化等离子体的反射系数,计算结果与解析解相一致,表明该吸收边界具有良好的吸收效果。     

一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法

基于伸展坐标完全匹配层方程和辅助微分方程 方法, 给出了一种在时域有限差分(FDTD) 计算中适用于常见色散介质的通用边界条件算法. 该算法适用于任意阶的FDTD空间差分, 并且由于所采用的D-H方程独立于计算区域, 所以可以被用于截断任意电介质. 数值试验结果表明, 与卷积完全匹配层 算法相比较, 所提出的吸收边界条件算法不仅通用性强、 计算复杂度低、 计算时间短, 并且吸收效果有明显的提高.     

一种色散介质FDTD通用吸收边界

由单轴各向异性介质所满足的场方程出发,并根据相位匹配原理以得到介质中的色散关系和无反射条件.结合频域到时域的转换关系（即用/t代替jω）和移位算子时域有限差分（SO-FDTD）方法提出了一种适用于各向同性常见色散介质模型,包括德拜模型、洛伦兹模型、德鲁德模型等的通用UPML吸收边界.数值结果表明了该算法的通用性和高效性. 关键词： 时域有限差分 吸收边界 移位算子 色散介质     

Accurate absorbing boundary conditions for anisotropic elastic media. Part 2: Untilted non-elliptic anisotropy

With the ultimate goal of devising effective absorbing boundary conditions (ABCs) for general elastic media, we investigate the accuracy aspects of local ABCs designed for untilted non-elliptic anisotropy in the frequency domain (time-harmonic analysis). While simple space–time transformations are available to treat the wavemodes with opposing phase and group velocities present in elliptic anisotropic media, no such transformations are known to exist for the case of non-elliptic anisotropy. In this paper, we use the concept of layer groupings along with an unconventional stretching of the finite element mesh to guarantee the accuracy of local ABCs designed to treat all propagating wavemodes, even those with opposing phase and group velocities. The local ABC used here is the perfectly matched discrete layer (PMDL) which is a simple variant of perfectly matched layers (PMLs) that is also equivalent to rational approximation-based local ABCs (rational ABCs); it inherits the straightforward approximation properties of rational ABCs along with the versatility of PML. The approximation properties of PMDL quantified through its reflection matrix allow us to (a) show that it is impossible to design an accurate PMDL with wavenumber-independent parameters, (b) theoretically demonstrate the ability of wavenumber-dependent parameters to ensure accuracy, and finally (c) design a practical though unconventional stretching of the finite element PMDL mesh that facilitates the implementation of wavenumber-dependent parameters. The validity of this work is demonstrated through a series of numerical experiments.     

有耗介质空间完全匹配层吸收边界条件及其应用

针对Gedney提出的完全匹配层(PML)无法用于有耗各向同性计算域的截断及其对倏逝波的衰减不理想等问题,提出了一种扩展方法。扩展的PML的主要思想是在各向异性的PML中引入与有耗介质空间相一致的复介电常数和复磁导率,使之可以与有耗介质计算域相匹配。通过给PML的张量介电常数、张量磁导率增加衰减因子以加速倏逝波的衰减。构造了PML吸收效果验证模型,数值结果证明了扩展的PML在处理有耗介质计算域截断问题中的有效性。利用该吸收边界条件,采用时域有限差分法计算了电磁脉冲作用下地面铺设电缆的电磁脉冲响应,计算结果和试验结果取得了较好的一致。     

二维单步交替方向隐式时域有限差分法吸收边界性能分析

给出了两种适用于二维单步交替方向隐式时域有限差分(2-D Leapfrog ADI-FDTD)方法的吸收边界：Mur边界和卷积完全匹配层(CPML)边界。单步交替方向隐式时域有限差分(Leapfrog ADI-FDTD)方法是一种无条件稳定的全隐式差分算法,由于二维空间Leapfrog ADI-FDTD的迭代同时存在显式和隐式方程,故而不同电磁分量的边界条件也存在差异。从原理出发,推导了适用于2-D Leapfrog ADI-FDTD方法的CPML边界条件,并与一阶Mur边界进行比较,利用自由空间的反射误差来表征两种边界的吸收性能,简要总结了两种吸收边界的优缺点。     

各向异性磁等离子体的辅助方程FDTD算法

通过对辅助方程的差分迭代，给出了一种新的各向异性磁化等离子体的辅助方程时域有限差分(FDTD)算法.并与电流密度卷积(JEC)算法进行了比较.通过计算磁化等离子体平板对平行于磁场传播的电磁波的反射和透射系数，验证了该算法的精度与JEC算法基本相同，而计算效率提高数倍. 关键词： 色散介质 时域有限差分算法 各向异性 磁化等离子体     

Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

声学中用于时域有限差分法的一种高效吸收边界

本文提出一种用于声场时域有限差分法(FDTD)的高效吸收边界条件。通过伪逆矩阵对截断边界附近网格线之间的转移矩阵作最小二乘估计,所得到的吸收边界条件性能明显优于多种常用方法。新方法的数值计算效率高,算法稳定。对点源辐射,尖劈衍射和圆柱体散射问题的FDTD数值计算结果显示了这种方法的优异性能。     

An algorithm for solving the optical problem for stratified anisotropic media

An algorithm for solving the Maxwell equations for propagation of light through anisotropic stratified media is considered. The algorithm uses the Berreman matrices of order 4 × 4. In contrast to the numerical methods suggested by Berreman, the new method is exact. The Sylvester theorem for calculating functions of a matrix and the Laguerre method for determining eigenvalues provide the basis for an algorithm with an efficiency comparable to that of the algorithms based on analytic solutions, which exist only in the case of uniaxial media. The method suggested in this paper allows for the analysis of complex optical systems where the effects of biaxiality, magnetic anisotropy, and optical activity play an important role.     

Phonon absorbing boundary conditions for molecular dynamics

With the goal of minimizing the domain size for molecular dynamics (MD) simulations, we develop a new class of absorbing boundary conditions (ABCs) that mimic the phonon absorption properties of an unbounded exterior. The proposed MD-ABCs are extensions of perfectly matched discrete layers (PMDLs), originally developed as an absorbing boundary condition for continuous wave propagation problems. Called MD-PMDL, this extension carefully targets the absorption of phonons, the high frequency waves, whose propagation properties are completely different from continuous waves. This paper presents the derivation of MD-PMDL for general lattice systems, followed by explicit application to one-dimensional and two-dimensional square lattice systems. The accuracy of MD-PMDL for phonon absorption is proven by analyzing reflection coefficients, and demonstrated through numerical experiments. Unlike existing MD-ABCs, MD-PMDL is local in both space and time and thus more efficient. Based on their favorable properties, it is concluded that MD-PMDL could provide a more effective alternative to existing MD-ABCs.     

Flux limiting embedded boundary technique for electromagnetic FDTD

A general approach for incorporating embedded boundaries into an electromagnetic finite difference time domain (FDTD) code is presented. This algorithm is shown to satisfy Gauss’s law and enforces no magnetic monopoles while maintaining a globally second-order result (first-order at physical boundaries), with no added time-step restriction. Theoretically predicted superior results are shown with an 11% time-step reduction from the Courant stability limit. This is achieved through a physics-based flux limiting scheme near physical boundaries. Stability, local truncation error and energy conservation analysis are also provided.     

The kinetic boundary layer for the Fokker-Planck equation with absorbing boundary

The Fokker-Planck equation for the distribution of position and velocity of a Brownian particle is a particularly simple linear transport equation. Its normal solutions and an apparently complete set of stationary boundary layer solutions can be determined explicitly. By a numerical algorithm we select linear combinations of them that approximately fulfill the boundary condition for a completely absorbing plane wall, and that approach a linearly increasing position space density far from the wall. Various aspects of these approximate solutions are discussed. In particular we find that the extrapolated asymptotic density reaches zero at a distance x_M beyond the wall. We find x_M=1.46 in units of the velocity persistence length of the Brownian particle. This study was motivated by certain problems in the theory of diffusion-controlled reactions, and the results might be used to test approximate theories employed in that field.     

Absorbing boundary conditions for scalar waves in anisotropic media. Part 2: Time-dependent modeling

With the ultimate goal of devising effective absorbing boundary conditions (ABCs) for general anisotropic media, we investigate the well-posedness and accuracy aspects of local ABCs designed for the transient modeling of the scalar anisotropic wave equation. The ABC analyzed in this paper is the perfectly matched discrete layers (PMDL), a simple variant of perfectly matched layers (PML) that is also equivalent to rational approximation based ABCs. Specifically, we derive the necessary and sufficient condition for the well-posedness of the initial boundary value problem (IBVP) obtained by coupling an interior and a PMDL ABC. The derivation of the reflection coefficient presented in a companion paper (S. Savadatti, M.N. Guddati, J. Comput. Phys., 2010, doi:10.1016/j.jcp.2010.05.018) has shown that PMDL can correctly identify and accurately absorb outgoing waves with opposing signs of group and phase velocities provided the PMDL layer lengths satisfy a certain bound. Utilizing the well-posedness theory developed by Kreiss for general hyperbolic IBVPs, and the well-posedness conditions for ABCs derived by Trefethen and Halpern for isotropic acoustics, we show that this bound on layer lengths also ensures well-posedness. The time discretized form of PMDL is also shown to be theoretically stable and some instability related to finite precision arithmetic is discussed.