FDTD方法吸收边界条件的研究及应用

用时域有限差分法(FDTD)求解电磁散射问题中,吸收边界条件的设置起着关键性作用.通过时间和空间上的递推算法对时域有限差分法中的两种吸收边界条件:Mur吸收边界条件和完全匹配层(PML)的吸收效果进行了比较和分析.同时,引入参数对PML的差分方程进行了优化,避免了将电磁场分裂为两个分量进行计算,进而降低了计算内存开销.实验结果证明PML具有更优越的吸收性能.最后,在FDTD算法中应用PML吸收层对一圆柱形导体的雷达散射截面积(RCS)进行数值仿真,验证了FDTD算法在计算雷达散射截面积(RCS)上的有效性.     

Unconditionally stable higher order perfectly matched layer applied to terminate anisotropic magnetized plasma

An effective and efficient absorbing boundary condition applied to terminate anisotropic magnetized plasma is proposed in this article. As its mathematical derivations are based on the high order perfectly matched layer (HO‐PML) and Crank‐Nicolson finite‐difference time‐domain method, it has the advantages of reducing late‐time reflections, attenuating evanescent waves, absorbing low‐frequency propagation waves and overcoming the Courant‐Friedrichs‐Levy condition. To validate the efficiency and effectiveness of the proposed HO‐PML, a ridge waveguide model composed of vacuum and the anisotropic magnetized plasma, a dielectric waveguide model composed of the dielectric and the anisotropic magnetized plasma are employed. The numerical results show that the proposal can not only maintain the stability of the algorithm with the increment of the time step but also further enhance the absorbing performance.     

Development of PML absorbing boundary conditions for computational aeroacoustics: A progress review

Recent advances in the development of perfectly matched layer (PML) as absorbing boundary conditions for computational aeroacoustics are reviewed. The PML methodology is presented as a complex change of variables. In this context, the importance of a proper space-time transformation in the PML technique for Euler equations is emphasized. A unified approach for the derivation of PML equations is offered that involves three essential steps. The three-step approach is illustrated in details for the PML of linear and non-linear Euler equations. Numerical examples are also given that include non-reflecting boundary conditions for a ducted channel flow and mixing layer roll-up vortices.     

基于正交轨迹的共形PML数值仿真方法

从材料的观点看,完全匹配层是各向异性的介质.场在完全匹配层遵从麦克斯韦方程.通过设计磁导率和介电常数张量参数来实现任意形状的共形完全匹配层.为匹配层动态稳定,可应用二维正交轨迹网格和FDTD方法实现了数值共形完全匹配层仿真,数值仿真结果表明基于正交轨迹网格的共形完全匹配层为一个共形FDTD或FEM方法提供了一个有效地吸收边界条件.     

蝶形天线的计算机仿真设计

蝶形天线是在脉冲型探地雷达中广泛采用的一种宽带天线。本文采用时域有限差分算法(FDTD)结合PML吸收边界条件分析了蝶形天线在高斯脉冲激励下的时域特性,通过傅立叶变换,计算出天线的方向图和在不同频率下的输入阻抗,结果表明FDTD算法用于分析蝶形天线是有效的。     

Perfectly matched layers for time-harmonic elastodynamics of unbounded domains: theory and finite-element implementation

One approach to the numerical solution of a wave equation on an unbounded domain uses a bounded domain surrounded by an absorbing boundary or layer that absorbs waves propagating outwards from the bounded domain. A perfectly matched layer (PML) is an unphysical absorbing layer model for linear wave equations that absorbs, almost perfectly, outgoing waves of all non-tangential angles-of-incidence and of all non-zero frequencies. This paper develops the PML concept for time-harmonic elastodynamics in Cartesian coordinates, utilising insights obtained with electromagnetics PMLs, and presents a novel displacement-based, symmetric finite-element implementation of the PML for time-harmonic plane-strain or three-dimensional motion. The PML concept is illustrated through the example of a one-dimensional rod on elastic foundation and through the anti-plane motion of a two-dimensional continuum. The concept is explored in detail through analytical and numerical results from a PML model of the semi-infinite rod on elastic foundation, and through numerical results for the anti-plane motion of a semi-infinite layer on a rigid base. Numerical results are presented for the classical soil–structure interaction problems of a rigid strip-footing on a (i) half-plane, (ii) layer on a half-plane, and (iii) layer on a rigid base. The analytical and numerical results obtained for these canonical problems demonstrate the high accuracy achievable by PML models even with small bounded domains.     

带矢量ABC底面的共形完全匹配层

共形完全匹配层是一种有耗各向异性媒质组成的凸且光滑的壳体,其底面一般是PEC面或PMC面,但是PEC面或PMC面会对原散射场产生反射;为了减少底面反射,将CPML原有的PEC（或PMC）底面改为矢量ABC吸收边界,并给出了带矢量ABC底面的CPML泛函公式。通过数值算例证明,这种带矢量ABC底面的CPML边界不仅减少了底面反射,而且吸收效果好,计算精度高。     

Hybrid Spectral Difference Methods for Elliptic Equations on Exterior Domains with the Discrete Radial Absorbing Boundary Condition

The hybrid spectral difference methods (HSD) for the Laplace and Helmholtz equations in exterior domains are proposed. We consider the fictitious domain method with the absorbing boundary conditions (ABCs). The HSD method is a finite difference version of the hybridized Galerkin method, and it consists of two types of finite difference approximations; the cell finite difference and the interface finite difference. The fictitious domain is composed of two subregions; the Cartesian grid region and the boundary layer region in which the radial grid is imposed. The boundary layer region with the radial grid makes it easy to implement the discrete radial ABC. The discrete radial ABC is a discrete version of the Bayliss–Gunzburger–Turkel ABC without pertaining any radial derivatives. Numerical experiments confirming efficiency of our numerical scheme are provided.     

单极天线的计算机仿真设计

aa输入阻抗,在PML边界条件采用指数差分的形式,减少了误差,与文献比较,获得比较好的结果。     

A coordinate transformation approach for efficient repeated solution of Helmholtz equation pertaining to obstacle scattering by shape deformations

A computational model is developed for efficient solutions of electromagnetic scattering from obstacles having random surface deformations or irregularities (such as roughness or randomly-positioned bump on the surface), by combining the Monte Carlo method with the principles of transformation electromagnetics in the context of finite element method. In conventional implementation of the Monte Carlo technique in such problems, a set of random rough surfaces is defined from a given probability distribution; a mesh is generated anew for each surface realization; and the problem is solved for each surface. Hence, this repeated mesh generation process places a heavy burden on CPU time. In the proposed approach, a single mesh is created assuming smooth surface, and a transformation medium is designed on the smooth surface of the object. Constitutive parameters of the medium are obtained by the coordinate transformation technique combined with the form-invariance property of Maxwell’s equations. At each surface realization, only the material parameters are modified according to the geometry of the deformed surface, thereby avoiding repeated mesh generation process. In this way, a simple, single and uniform mesh is employed; and CPU time is reduced to a great extent. The technique is demonstrated via various finite element simulations for the solution of two-dimensional, Helmholtz-type and transverse magnetic scattering problems.     

Design of absorbing material distribution for sound barrier using topology optimization

A topology optimization approach based on the boundary element method (BEM) and the optimality criteria (OC) method is proposed for the optimal design of sound absorbing material distribution within sound barrier structures. The acoustical effect of the absorbing material is simplified as the acoustical impedance boundary condition. Based on the solid isotropic material with penalization (SIMP) method, a topology optimization model is established by selecting the densities of absorbing material elements as design variables, volumes of absorbing material as constraints, and the minimization of sound pressure at reference surface as design objective. A smoothed Heaviside-like function is proposed to help the SIMP method to obtain a clear 0–1 distribution. The BEM is applied for acoustic analysis and the sensitivities with respect to design variables are obtained by the direct differentiation method. The Burton–Miller formulation is used to overcome the fictitious eigen-frequency problem for exterior boundary-value problems. A relaxed form of OC is used for solving the optimization problem to find the optimal absorbing material distribution. Numerical tests are provided to illustrate the application of the optimization procedure for 2D sound barriers. Results show that the optimal distribution of the sound absorbing material is strongly frequency dependent, and performing an optimization in a frequency band is generally needed.     

Local absorbent boundary condition for non-linear hyperbolic problems with unknown Riemann invariants

Generally, in problems where the Riemann invariants (RI) are known (e.g. the flow in a shallow rectangular channel, the isentropic gas flow equations), the imposition of non-reflective boundary conditions is straightforward. In problems where Riemann invariants are unknown (e.g. the flow in non-rectangular channels, the stratified 2D shallow water flows) it is possible to impose that kind of conditions analyzing the projection of the Jacobians of advective flux functions onto normal directions of fictitious surfaces or boundaries. In this paper a general methodology for developing absorbing boundary conditions for non-linear hyperbolic advective–diffusive equations with unknown Riemann invariants is presented. The advantage of the method is that it is very easy to implement in a finite element code and is based on computing the advective flux functions (and their Jacobian projections), and then, imposing non-linear constraints via Lagrange multipliers. The application of the dynamic absorbing boundary conditions to typical wave propagation problems with unknown Riemann invariants, like non-linear Saint-Venant system of conservation laws for non-rectangular and non-prismatic 1D channels and stratified 1D/2D shallow water equations, is presented. Also, the new absorbent/dynamic condition can handle automatically the change of Jacobians structure when the flow regime changes from subcritical to supercritical and viceversa, or when recirculating zones are present in regions near fictitious walls.     

Perfectly Matched Layers for Time-Harmonic Second Order Elliptic Problems

The main goal of this work is to give a review of the Perfectly Matched Layer (PML) technique for time-harmonic problems. Precisely, we focus our attention on problems stated in unbounded domains, which involve second order elliptic equations writing in divergence form and, in particular, on the Helmholtz equation at low frequency regime. Firstly, the PML technique is introduced by means of a simple porous model in one dimension. It is emphasized that an adequate choice of the so called complex absorbing function in the PML yields to accurate numerical results. Then, in the two-dimensional case, the PML governing equation is described for second order partial differential equations by using a smooth complex change of variables. Its mathematical analysis and some particular examples are also included. Numerical drawbacks and optimal choice of the PML absorbing function are studied in detail. In fact, theoretical and numerical analysis show the advantages of using non-integrable absorbing functions. Finally, we present some relevant real life numerical simulations where the PML technique is widely and successfully used although they are not covered by the standard theoretical framework.     

Long-time behavior of PML absorbing boundaries for layered periodic structures

In this work we consider a special case of the Perfectly Matched Layer (PML) divergence which is observed by the simulation of the planar periodic structures such as photonic crystal slabs or antenna arrays. This divergence is caused by an excitation of long-living artefact evanescent waves in these structures by an incident external pulse. We study the application of the known remedies to this problem: increasing the distance between the structure and PML, employing the κ parameter, employing non-PML absorbers. We also suggest a new simple and effective solution, where the usual PML is backed by an additional absorbing layer.     

Computing the Bidirectional Scattering of a Microstructure Using Scalar Diffraction Theory and Path Tracing

Most models for bidirectional surface scattering by arbitrary explicitly defined microgeometry are either based on geometric optics and include multiple scattering but no diffraction effects or based on wave optics and include diffraction but no multiple scattering effects. The few exceptions to this tendency are based on rigorous solution of Maxwell's equations and are computationally intractable for surface microgeometries that are tens or hundreds of microns wide. We set up a measurement equation for combining results from single scattering scalar diffraction theory with multiple scattering geometric optics using Monte Carlo integration. Since we consider an arbitrary surface microgeometry, our method enables us to compute expected bidirectional scattering of the metasurfaces with increasingly smaller details seen more and more often in production. In addition, we can take a measured microstructure as input and, for example, compute the difference in bidirectional scattering between a desired surface and a produced surface. In effect, our model can account for both diffraction colors due to wavelength-sized features in the microgeometry and brightening due to multiple scattering. We include scalar diffraction for refraction, and we verify that our model is reasonable by comparing with the rigorous solution for a microsurface with half ellipsoids.     

奇异摄动反应扩散方程数值模拟的粒子群优化算法

针对Shishkin网格方法在数值求解奇异摄动反应扩散方程时,网格过度点参数的选取具有不确定性的缺陷,提出了一种用粒子群优化(PSO)算法估计Shishkin网格参数的方法。首先基于有限差分方法,构造了以误差范数最小为目标的无约束优化问题,并用PSO算法进行了求解。该方法克服了人为选择参数的缺陷。实验结果表明：与单纯形算法相比,PSO算法在优化Shishkin网格参数时能够收敛到全局最优解;而且在最优网格参数下,奇异摄动反应扩散方程的数值结果在边界层的精度也得到了明显提高,进一步说明了所提方法的有效性和可行性。     

Discretization error due to the identity operator in surface integral equations

We consider the accuracy of surface integral equations for the solution of scattering and radiation problems in electromagnetics. In numerical solutions, second-kind integral equations involving well-tested identity operators are preferable for efficiency, because they produce diagonally-dominant matrix equations that can be solved easily with iterative methods. However, the existence of the well-tested identity operators leads to inaccurate results, especially when the equations are discretized with low-order basis functions, such as the Rao-Wilton-Glisson functions. By performing a computational experiment based on the nonradiating property of the tangential incident fields on arbitrary surfaces, we show that the discretization error of the identity operator is a major error source that contaminates the accuracy of the second-kind integral equations significantly.     

A parallel fictitious domain multigrid preconditioner for the solution of Poisson’s equation in complex geometries

A parallel multilevel preconditioner based on domain decomposition and fictitious domain methods has been presented for the solution of the Poisson equation in complicated geometries. Rectangular blocks with matching grids on interfaces on a structured rectangular mesh have been used for the decomposition of the problem domain. Sloping sides or curved boundary surfaces are approximated using stepwise surfaces formed by the grid cells. A seven-point stencil based on the central difference scheme has been used for the discretization of the Laplacian for both interior and boundary grid points, and this results in a symmetric linear algebraic system for any type of boundary condition. The preconditioned conjugate gradient method has been used for the solution of this symmetric system. The multilevel preconditioner for the CG is based on a V-cycle multigrid applied to the Poisson equation on a fictitious domain formed by the union of the rectangular blocks used for the domain decomposition. Numerical results are presented for two typical Poisson problems in complicated geometries—one related to heat conduction, and the other one arising from the LES/DNS of incompressible turbulent flow over a packed array of spheres. These results clearly show the efficiency and robustness of the proposed approach.     

Long Time Behavior of the Perfectly Matched Layer Equations in Computational Electromagnetics

We investigate the long time behavior of two unsplit PML methods for the absorption of electromagnetic waves. Computations indicate that both methods suffer from a temporal instability after the fields reach a quiescent state. The analysis reveals that the source of the instability is the undifferentiated terms of the PML equations and that it is associated with a degeneracy of the quiescent systems of equations. This highlights why the instability occurs in special cases only and suggests a remedy to stabilize the PML by removing the degeneracy. Computational results confirm the stability of the modified equations and is used to address the efficacy of the modified schemes for absorbing waves.     

Solution of time-domain Maxwell equation with PML by using modified Laguerre polynomials

In this work, a stable numerical algorithm proposed by Chung et al. for the time-domain Maxwell equations is generalized. The time-domain Maxwell equations are solved by expressing the transient behaviors in terms of the modified Laguerre polynomials, and then the original equations of the initial value and boundary value can be transformed into a series of problems independent of the time variable. In this case the method of finite difference （FD）, the finite element method （FEM）, the method of moment （MoM）, etc. or the combination of these methods can be used to solve the problems. Finally, a numerical model is provided for the scattering problem with perfect matched layer （PML） by using FD. The comparison between the results of the proposed method and FDTD is presented to verify the proposed new method.