DOCSLIB.ORG

Simulation of the Finite Difference Time Domain in Two Dimension

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Simulation of the Finite Difference Time Domain in Two Dimension World Academy of Science, Engineering and Technology International Journal of Computer and Information Engineering Vol:6, No:3, 2012 Simulation of the Finite Difference Time Domain in Two Dimension Akram G., Jasmy Y. ∂H ∂(EE + ) Abstract—The finite-difference time-domain (FDTD) method is µy + σ * = zx zy xH y (2) one of the most widely used computational methods in ∂t ∂ x electromagnetic. This paper describes the design of two-dimensional (2D) FDTD simulation software for transverse magnetic (TM) ∂ ∂E H y polarization using Berenger's split-field perfectly matched layer εzx + σ E = (3) (PML) formulation. The software is developed using Matlab ∂tx zx ∂ x programming language. Numerical examples validate the software. ∂ Keywords—Finite difference time domain (FDTD) method, Ezy ∂H ε+ σ E = − x (4) perfectly matched layer (PML), split-filed formulation, transverse ∂ty zy ∂ y magnetic (TM) polarization. I. INTRODUCTION In order to solve these equations numerically, they are discretized according to finite differencing techniques. The HE finite-difference time-domain (FDTD) method is a central difference approximation is used in this case. For flexible and powerful solution tool for solving Maxwell’s T example, in equations (1) and (3), central differencing results equations [1-5]. Further, the efficient and accurate solution of in the following equations electromagnetic wave interaction problems in unbounded regions is one of the greatest challenges of the FDTD method. + −σ* () + δ µ Hn 1 () i+1 2, j = ey i1 2, j t Hn () i +1 2, j For such problems, an absorbing boundary condition (ABC) x x (5) −σ* () + δ µ must be introduced at the outer grid boundary to simulate the (1− e y i1 2, j t ) n+1 2()+ + + n + 1 2 () + + Ezx i1 2, j 1 2 E zy i1 2, j 1 2 extension of the grid to infinity. − σ* ()i+1 2, j δ x −En+1 2() i +1 2, j − 1 2 − E n + 1 2 () i + 1 2, j − 1 2 Recently, among a few important studies on ABC was J.P. y zx zy Berenger’s Perfectly Matched Layer (PML) method [6]. In the + −σ* () + + δ ε − n 1 2 ()+ + =x i1 2, j 1 2 t n 1 2 () + + PML, he created an artificial medium with magnetic Ezx i1 2, j 1 2 e Ezx i1 2, j 1 2 conductivity (σ*) that will make the boundary condition to −σ* () + + δ ε (6) ()1− e x i1 2, j 1 2 t work as wave absorber region. This paper describes the design −Hn() i, j + 1 2 − H n () i + 1, j + 1 2 σ* ()+ + δ y y of a two-dimensional (2D) FDTD simulation software for x i1 2, j 1 2 x transverse magnetic (TM) polarized incident wave using Berenger's split-field perfectly matched layer (PML) Based on the above formulation, a 2-D FDTD code has formulation. The software is developed using Matlab been written in Matlab programming language to simulate the programming language. fields of a plane wave source in lossy media. The formulation and the FDTD code are validated by checking the numerical II. TWO-DIMENSIONAL TM POLARIZATION results for homogeneous media against the analytical solution. The Maxwell’s curl equations as modified by Berenger are The code for 2D FDTD is shown in Appendix. expressed in their time-dependent form as, In order to test and validate the FDTD code, an electromagnetic wave source interacting with dielectric International Science Index, Computer and Information Engineering Vol:6, No:3, 2012 waset.org/Publication/5827 ε σ ∂ + cylinder ( r = 4.0, = 0.12) of radius 6 cm is simulated in this ∂H (EEzx zy ) 2 µx + σ * H = − (1) section, shown in Fig. 1. A 1 mW/cm plane wave is used as ∂ty x ∂ y source, and the frequency is setup to 2.5 GHz. The reason for using dielectric cylinder is there exist an analytical solution to this problem [7]. A +y-directed plane wave source : G. Akram is with the Faculty of Health Science & Biomedical Eng., , ezy(:,jebc+2) = ezy(:,jebc+2) + 61.4*sin(omega*n*dt) Universiti Teknologi Malaysia (UTM), 81310 Skudai, Johor, Malaysia (phone: 607-55-35785; fax: 607-55-66272; e-mail: akram@ fke.utm.my). is used to generate the incident wave at j = jebc + 2 and n Y. Jasmy is with the Faculty of Health Science & Biomedical Eng., = 1. Where jebc is thickness of the PML region, n is time step Universiti Teknologi Malaysia (UTM), 81310 Skudai, Johor, Malaysia. International Scholarly and Scientific Research & Innovation 6(3) 2012 360 scholar.waset.org/1307-6892/5827 World Academy of Science, Engineering and Technology International Journal of Computer and Information Engineering Vol:6, No:3, 2012 in FDTD code and dt is 5.0 psec. Fig. 2 shows a comparison of FDTD calculated results with analytic solution along the center axis o the cylinder. The simulation was made with a code attached in Appendix. Fig. 1 Schematic representation of the FDTD computational domain Fig. 3 Test and benchmark computational domains 50 Then, the local error of the computed field in Ω , due to analytical solution T FDTD + PML ( 16 ) reflections from the PML, is obtained by subtracting the field 45 FDTD + PML ( 8 ) at any point within ΩT at a given time step from the field at the corresponding point in ΩB. Here, Ez is used to define the error 40 n n n |E| e= E − E (8) locali, j z,, Ti,, j z B i j 35 where Ez,T and Ez,B are, respectively, the FDTD computed E- 30 fields within the test and benchmark domains. Further, the global error was defined as 25 5 10 15 20 25 30 35 40 45 n n n 2 Grid Position, I e= E − E (9) global∑∑ z,, Ti,, j z B i j Fig. 2 Waveform comparison of the FDTD results (with different i j PML thicknesses) and the analytical solution for a plane wave source In Fig. 4 and 5, different numbers of PML are used outside III. PROCEDURE FOR PAPER SUBMISSION the FDTD simulation region to test the efficiency of the PML In this section, various numerical simulations are performed ABC for the lossy media. The accuracy of the PML ABC is to illustrate the effectiveness of PML code [8]. Fig. 3 shows compared with that of standard second-order Mur ABC [9]. the two FDTD computational domains used in the The numbers of PML used are 16, 8 and 4, respectively. As is experiments: a test domain, ΩT, and the much larger evident from the result, 16 PML are good enough to absorb benchmark domain ΩB. By calculating the difference between most of the incident waves at the boundaries. the FDTD solutions in the two domains at each grid-point at each time-step, a measure of the spurious reflection caused by IV. CONCLUSION the test ABC is obtained. Initially, a 2-D test domain 100×50 The FDTD method is widely used because it is simple to cells long was examined. A square cell was used, with δ implement numerically. It provides a flexible means for = 0.015 m, and the time step was chosen based on the stability directly solving Maxwell’s time-dependent curl equations by criterion δt= δ / 2 c , where c is the speed of light in vacuum. using finite differences to discretize them. This paper In all the cases examined here, the excitation was a pulse International Science Index, Computer and Information Engineering Vol:6, No:3, 2012 waset.org/Publication/5827 describes the design of two-dimensional (2D) FDTD exhibiting a very smooth transition to zero, as used in [8], and simulation software for transverse magnetic (TM) polarized defined as follows, incident wave using Berenger's split-field perfectly matched 1 10− 15cos(π n / 20) + layer (PML) formulation. The software is developed using n n ≤ 40 = ()()π− π (7) Matlab programming language. The numerical results agree Ez, T 32 6cos 2n / 20 cos 3 n / 20 50,25 very well with the analytical results for homogeneous media. 0 n > 40 International Scholarly and Scientific Research & Innovation 6(3) 2012 361 scholar.waset.org/1307-6892/5827 World Academy of Science, Engineering and Technology International Journal of Computer and Information Engineering Vol:6, No:3, 2012 0 10 iefbc = ie + 2*iebc; jefbc = je + 2*jebc; ibfbc = iefbc + 1; -10 10 jbfbc = jefbc + 1; %********************************************** % Material parameters -20 10 %********************************************** media = 2; Mur ABC eps = [1.0 4.0]; -30 PML ABC (16 cell) 10 sig = [0.0 0.12]; Global Error Global PML ABC (8 cell) PML ABC (4 cell) mur = [1.0 1.0]; sim = [0.0 0.0]; -40 10 %********************************************** % Field arrays %********************************************** -50 10 ez = zeros(ibfbc,jbfbc); %fields in main grid 0 100 200 300 400 500 ezx = zeros(ibfbc,jbfbc); Time Step ezy = zeros(ibfbc,jbfbc); Fig. 4 Global error within a 100×50 cell two-dimensional TM test hx = zeros(iefbc,jefbc); hy = zeros(iefbc,jefbc); grid -5 %********************************************** x 10 2 % Updating coefficients %********************************************** 1.5 for i = 1:media eaf = dt*sig(i)/(2.0*epsz*eps(i)); ca(i) = (1.0-eaf)/(1.0+eaf); 1 cb(i) = dt/epsz/eps(i)/dx/(1.0+eaf); haf = dt*sim(i)/(2.0*muz*mur(i)); 0.5 da(i) = (1.0-haf)/(1.0+haf); db(i) = dt/muz/mur(i)/dx/(1.0+haf); 0 end %********************************************** -0.5 PML ABC ( 16 cell ) % Geometry specification (main grid) PML ABC ( 8 cell ) %********************************************** Normalized Reflection Along J =1 J Along Reflection Normalized -1 % Initialize entire main grid to free space caezx(1:ibfbc,1:jbfbc) = ca(1); -1.5 cbezx(1:ibfbc,1:jbfbc) = cb(1); 0 20 40 60 80 100 caezy(1:ibfbc,1:jbfbc) = ca(1); Grid Position, I cbezy(1:ibfbc,1:jbfbc) = cb(1); Fig.
Load more

Recommended publications
  • Cvel-08-011.3
    TECHNICAL REPORT: CVEL-08-011.3 Survey of Current Computational Electromagnetics Techniques and Software T. Hubing, C. Su, H. Zeng and H. Ke Clemson University June 6, 2009 This report is an update of the 2008 report, CVEL-08-011.2. This work was supported by the ARMY RESEARCH OFFICE (W. Devereux Palmer) under the auspices of the U.S. Army Research Office Scientific Services Program administered by Battelle (Delivery Order 0271, Contract No. W911NF-07-D-0001). Table of Contents I. Introduction ............................................................................................................................. 3 A. Purpose of This Report ................................................................................................................. 3 B. Categorizing CEM Modeling Tools .............................................................................................. 4 1. Time vs. Frequency Domain ........................................................................................ 4 2. IE vs. PDE .................................................................................................................... 4 3. 2D vs. 3D ...................................................................................................................... 5 4. Static, Quasi-static, Full-wave or Asymptotic .............................................................. 5 5. Circuit-based Field Solvers .......................................................................................... 6 6. Linear vs. Higher-order Methods
    [Show full text]
  • Design of an Ultrahigh Birefringence Photonic Crystal Fiber with Large Nonlinearity Using All Circular Air Holes for a Fiber-Optic Transmission System
    hv photonics Article Design of an Ultrahigh Birefringence Photonic Crystal Fiber with Large Nonlinearity Using All Circular Air Holes for a Fiber-Optic Transmission System Shovasis Kumar Biswas 1,* ID , S. M. Rakibul Islam 1, Md. Rubayet Islam 1, Mohammad Mahmudul Alam Mia 2, Summit Sayem 1 and Feroz Ahmed 1 ID 1 Department of Electrical and Electronic Engineering, Independent University Bangladesh, Dhaka 1229, Bangladesh; [email protected] (S.M.R.I.); [email protected] (M.R.I.); [email protected] (S.S.); [email protected] (F.A.) 2 Department of Electronics & Communication Engineering, Sylhet International University, Sylhet 3100, Bangladesh; [email protected] * Correspondence: [email protected]; Tel.: +880-1717-827-665 Received: 11 August 2018; Accepted: 29 August 2018; Published: 5 September 2018 Abstract: This paper proposes a hexagonal photonic crystal fiber (H-PCF) structure with all circular air holes in order to simultaneously achieve ultrahigh birefringence and high nonlinearity. The H-PCF design consists of an asymmetric core region, where one air hole is a reduced diameter and the air hole in its opposite vertex is omitted. The light-guiding properties of the proposed H-PCF structure were studied using the full-vector finite element method (FEM) with a circular perfectly matched layer (PML). The simulation results showed that the proposed H-PCF exhibits an ultrahigh birefringence of 3.87 × 10−2, a negative dispersion coefficient of −753.2 ps/(nm km), and a nonlinear coefficient of 96.51 W−1 km−1 at an excitation wavelength of 1550 nm. The major advantage of our H-PCF design is that it provides these desirable modal properties without using any non-circular air holes in the core and cladding region, thus making the fiber fabrication process much easier.
    [Show full text]
  • Recent Advances in Perfectly Matched Layers in Finite Element Applications
    Turk J Elec Engin, VOL.16, NO.1 2008, c TUB¨ ITAK˙ Recent Advances in Perfectly Matched Layers in Finite Element Applications Ozlem¨ OZG¨ UN¨ and Mustafa KUZUOGLU˘ Department of Electrical and Electronics Engineering, Middle East Technical University 06531, Ankara-TURKEY Abstract We present a comparative evaluation of two novel and practical perfectly matched layer (PML) implementations to the problem of mesh truncation in the finite element method (FEM): locally-conformal PML,and multi-center PML techniques. The most distinguished feature of these methods is the simplicity and flexibility to design conformal PMLs over challenging geometries,especially those with curvature discontinuities,in a straightforward way without using artificial absorbers. These methods are based on specially- and locally-defined complex coordinate transformations inside the PML region. They can easily be implemented in a conventional FEM by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. After overviewing the theoretical bases of these methods,we present some numerical results in the context of two- and three-dimensional electromagnetic radiation/scattering problems. Key Words: Finite element method (FEM); perfectly matched layer (PML); locally-conformal PML; multi-center PML; complex coordinate stretching; electromagnetic scattering. 1. Introduction Perfectly matched layers (PMLs) have been popular in the finite element method (FEM) and the finite difference time domain method (FDTD) for solving the domain truncation problem in electromagnetic radiation and scattering problems.The PML approach is based on the truncation of computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence.The main advantage of the PML is its close proximity and conformity to the surface of the radiator or scatterer, implying that the white space can be minimized.
    [Show full text]
  • On the Analysis of Perfectly Matched Layers for a Class of Dispersive
    On the analysis of perfectly matched layers for a class of dispersive media and application to negative index metamaterials Eliane Bécache, Patrick Joly, Valentin Vinoles To cite this version: Eliane Bécache, Patrick Joly, Valentin Vinoles. On the analysis of perfectly matched layers for a class of dispersive media and application to negative index metamaterials. 2016. hal-01327315v3 HAL Id: hal-01327315 https://hal.archives-ouvertes.fr/hal-01327315v3 Submitted on 12 Jul 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution - NonCommercial - NoDerivatives| 4.0 International License On the analysis of perfectly matched layers for a class of dispersive media and application to negative index metamaterials Éliane Bécache1,a, Patrick Joly1,b and Valentin Vinoles ∗2,c 1 Laboratoire Poems (UMR 7231 CNRS/Inria/ENSTA ParisTech), ENSTA ParisTech, 828, Boulevard des Maréchaux, 91762 Palaiseau, France 2École Polytechnique Fédérale de Lausanne, SB MATHAA CAMA, Station 8, CH-1015 Lausanne, Switzerland [email protected] [email protected] [email protected] Abstract This work deals with Perfectly Matched Layers (PMLs) in the context of dispersive media, and in particular for Negative Index Metamaterials (NIMs).
    [Show full text]
  • Thesis the Conformal Perfectly Matched Layer
    THESIS THE CONFORMAL PERFECTLY MATCHED LAYER FOR ELECTRICALLY LARGE CURVILINEAR HIGHER ORDER FINITE ELEMENT METHODS IN ELECTROMAGNETICS Submitted by Aaron P. Smull Department of Electrical and Computer Engineering In partial fulfillment of the requirements For the Degree of Master of Science Colorado State University Fort Collins, Colorado Summer 2017 Master’s Committee: Advisor: Branislav Notaros Ali Pezeshki Donald Estep Copyright by Aaron Patrick Smull 2017 All Rights Reserved ABSTRACT THE CONFORMAL PERFECTLY MATCHED LAYER FOR ELECTRICALLY LARGE CURVILINEAR HIGHER ORDER FINITE ELEMENT METHODS IN ELECTROMAGNETICS The implementation of open-region boundary conditions in computational electromagnetics for higher order finite element methods presents a well known set of challenges. One such boundary condition is known as the perfectly matched layer. In this thesis, the generation of perfectly matched layers for arbitrary convex geometric hexahedral meshes is discussed, using a method that can be implemented without differential operator based absorbing boundary conditions or coupling to boundary integral equations. A method for automated perfectly matched layer element generation is presented, with geometries based on surface projections from a convex mesh. Material parameters are generated via concepts from transformation electromagnetics, from complex-coordinate transformation based conformal PML’s in existing literature. A material parameter correction algorithm is also presented, based on a modified gradient descent optimization algorithm Numerical results are presented with comparison to analytical results and commercial software, with studies on the effects of discretization error of the effectiveness of the perfectly matched layer. Good agreement is found between simulated and analytical results, and between simulated results and commercial software. ii ACKNOWLEDGEMENTS I am grateful most of all to my family, whose support and love through the years has brought me continued fortune in my education.
    [Show full text]
  • Time-Domain Modeling of Elastic and Acoustic Wave Propagation in Unbounded Media, with Application to Metamaterials
    Time-domain modeling of elastic and acoustic wave propagation in unbounded media, with application to metamaterials by Hisham Assi A thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Graduate Department of Electrical and Computer Engineering collaborative program with the Institute of Biomaterial and Biomedical Engineering University of Toronto c Copyright 2016 by Hisham Assi Abstract Time-domain modeling of elastic and acoustic wave propagation in unbounded media, with application to metamaterials Hisham Assi Doctor of Philosophy Graduate Department of Electrical and Computer Engineering collaborative program with the Institute of Biomaterial and Biomedical Engineering University of Toronto 2016 Perfectly matched layers (PML) are a well-developed method for simulating wave propa- gation in unbounded media enabling the use of a reduced computational domain without having to worry about spurious boundary reflections. Many PML studies have been reported for both acoustic waves in fluids and elastic waves in solids. Nevertheless, further studies are needed for improvements in the fields of formulation, stability, and inhomogeneity of PMLs. This thesis introduces new second-order time-domain PML formulations for modeling mechanical wave propagation in unbounded solid, fluid, and coupled fluid-solid media. It also addresses certain stability issues, and demonstrates application of these formulations. Using a complex coordinate stretching approach a PML for the time-domain anisotropic elastic wave equation in two dimensional space is compactly formulated with two second- order equations along with only four auxiliary equations. This makes it smaller than existing formulations, thereby simplifying the problem and reducing the computational burden. A simple method is proposed for improving the stability of the discrete PML problem for a wide range of otherwise unstable anisotropic elastic media.
    [Show full text]
  • Perfectly Matched Layers for T,Φ Formulation
    Proceedings - 49 - 12th International IGTE Symposium Perfectly Matched Layers for T, Φ Formulation Gernot Matzenauer, Oszkar´ B´ıro,´ Karl Hollaus, Kurt Preis and Werner Renhart Institute for Fundamentals and Theory in Electrical Engineering, Kopernikusgasse 24/3, A-8010 Graz, Austria E-mail: [email protected] Abstract: Perfectly Matched Layers (PMLs) are used for II. PROOF OF THE PERFECT MATCHING PROPERTY reflectionless truncation of the problem boundaries in FEM applications. The basic concept behind PMLs is to create an Based on [2] it is shown below that the PML interface artificial material with a complex and diagonal anisotropic also works with any lossless material in front of the PML. permittivity and permeability. For the A,V formulation Similarly to (2), the permittivity and permeability tensors PMLs are well known. In this paper the method of perfectly of the PML are assumed to be: matched layers is extended to truncate any lossless medium and the method is implemented for the T, Φ formulation. [] = 0r[Λ], Keywords: Perfectly Matched Layer, Finite Element Method, T, Φ formulation, Perfectly Matched Anisotropic [µ] = µ0µr[Λ]. (4) Absorber, Absorbing Boundary Condition The question is how to choose the parameters a, b, c, µr and r, so that no reflection occurs on an interface between I. INTRODUCTION ∗ the PML and a medium with permeability µr and permit- ∗ In the present paper, an artificial anisotropic lossy tivity r. The wave equations for the electric and the mag- material has been applied to a 3D edge finite-element T, Φ netic field in the PML are: formulation [1] to form perfectly matched layers.
    [Show full text]
  • Benchmarking Five Numerical Simulation Techniques for Computing Resonance Wavelengths and Quality Factors in Photonic Crystal Membrane Line Defect Cavities
    Downloaded from orbit.dtu.dk on: Oct 11, 2021 Benchmarking five numerical simulation techniques for computing resonance wavelengths and quality factors in photonic crystal membrane line defect cavities de Lasson, Jakob Rosenkrantz; Frandsen, Lars Hagedorn; Gutsche, Philipp; Burger, Sven; Kim, Oleksiy S.; Breinbjerg, Olav; Ivinskaya, Aliaksandra; Wang, Fengwen; Sigmund, Ole; Häyrynen, Teppo Total number of authors: 13 Published in: Optics Express Link to article, DOI: 10.1364/OE.26.011366 Publication date: 2018 Document Version Publisher's PDF, also known as Version of record Link back to DTU Orbit Citation (APA): de Lasson, J. R., Frandsen, L. H., Gutsche, P., Burger, S., Kim, O. S., Breinbjerg, O., Ivinskaya, A., Wang, F., Sigmund, O., Häyrynen, T., Lavrinenko, A. V., Mørk, J., & Gregersen, N. (2018). Benchmarking five numerical simulation techniques for computing resonance wavelengths and quality factors in photonic crystal membrane line defect cavities. Optics Express, 26(9), 11366-11392. https://doi.org/10.1364/OE.26.011366 General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.
    [Show full text]
  • Stable Perfectly Matched Layers for a Cold Plasma in a Strong Background Magnetic Field Eliane Bécache, Patrick Joly, Maryna Kachanovska
    Stable Perfectly Matched Layers for a Cold Plasma in a Strong Background Magnetic Field Eliane Bécache, Patrick Joly, Maryna Kachanovska To cite this version: Eliane Bécache, Patrick Joly, Maryna Kachanovska. Stable Perfectly Matched Layers for a Cold Plasma in a Strong Background Magnetic Field . Journal of Computational Physics, Elsevier, 2017, 341, pp.76-101. 10.1016/j.jcp.2017.03.051. hal-01397581 HAL Id: hal-01397581 https://hal.archives-ouvertes.fr/hal-01397581 Submitted on 16 Nov 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Stable Perfectly Matched Layers for a Cold Plasma in a Strong Background Magnetic FieldI Eliane´ B´ecache, Patrick Joly, Maryna Kachanovska1 POEMS (ENSTA-INRIA-CNRS-Universit´eParis Saclay), ENSTA ParisTech, 828, Boulevard des Mar´echaux,91762 Palaiseau Cedex Abstract This work addresses the question of the construction of stable perfectly matched layers (PML) for a cold plasma in the infinitely large background magnetic field. We demonstrate that the traditional, B´erenger'sperfectly matched layers are unstable when applied to this model, due to the presence of the backward propagating waves. To overcome this instability, we use a combination of two techniques presented in the article.
    [Show full text]
  • (PML) Method for the Numerical Modeling and Simulation of Elastic Wave Propagation in Unbounded Domains Florent Pled, Christophe Desceliers
    Review and Recent Developments on the Perfectly Matched Layer (PML) Method for the Numerical Modeling and Simulation of Elastic Wave Propagation in Unbounded Domains Florent Pled, Christophe Desceliers To cite this version: Florent Pled, Christophe Desceliers. Review and Recent Developments on the Perfectly Matched Layer (PML) Method for the Numerical Modeling and Simulation of Elastic Wave Propagation in Unbounded Domains. Archives of Computational Methods in Engineering, Springer Verlag, 2021, 10.1007/s11831-021-09581-y. hal-03196974 HAL Id: hal-03196974 https://hal-upec-upem.archives-ouvertes.fr/hal-03196974 Submitted on 19 Apr 2021 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Archives of Computational Methods in Engineering manuscript No. (will be inserted by the editor) Review and recent developments on the perfectly matched layer (PML) method for the numerical modeling and simulation of elastic wave propagation in unbounded domains Florent Pled · Christophe Desceliers Received: 18 November 2020 / Accepted: 24 March 2021 Abstract This review article revisits and outlines the perfectly matched layer (PML) method and its various formulations developed over the past 25 years for the numerical modeling and simulation of wave propagation in unbounded me- dia.
    [Show full text]
  • A Reflectionless Discrete Perfectly Matched Layer
    A Reflectionless Discrete Perfectly Matched Layer Albert Chern Institute of Mathematics, Technical University of Berlin, 10623 Berlin, Germany Abstract Perfectly Matched Layer (PML) is a widely adopted non-reflecting boundary treatment for wave simulations. Re- ducing numerical reflections from a discretized PML has been a long lasting challenge. This paper presents a new discrete PML for the multi-dimensional scalar wave equation which produces no numerical reflection at all. The reflectionless discrete PML is discovered through a straightforward derivation using Discrete Complex Analysis. The resulting PML takes an easily-implementable finite difference form with compact stencil. In practice, the discrete waves are damped exponentially in the PML, and the error due to domain truncation is maintained at machine zero by a moderately thick PML. The numerical stability of the proposed PML is also demonstrated. Keywords: Perfectly matched layers, absorbing boundary, discrete complex analysis, scalar wave equation 1. Introduction Perfectly Matched Layer (PML) [1–3] is a state-of-the-art numerical technique for absorbing waves at the boundary of a computation domain, highly demanded by simulators of wave propagations in free space. More precisely, PML is an artificial layer attached to the boundary where the wave equation is modified into a set of so-called PML equations. Requiring simulations only local in space and time, the layer effectively damps the waves while creating no reflection at the transitional interface (Fig.1). Although the reflectionless property is well established in the continuous theory, to the best knowledge all existing discretizations of PML equations still produce numerical reflections at the interface due to discretization error (see e.g.
    [Show full text]
  • Improving the Performance of Perfectly Matched Layers by Means of Hp-Adaptivity C
    Improving the Performance of Perfectly Matched Layers by Means of hp-Adaptivity C. Michler, L. Demkowicz, J. Kurtz, D. Pardo Institute for Computational Engineering and Sciences (ICES), The University of Texas at Austin, Austin, Texas 78712 Received 10 February 2006; accepted 31 January 2007 Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.20252 We improve the performance of the Perfectly Matched Layer (PML) by using an automatic hp-adaptive discretization. By means of hp-adaptivity, we obtain a sequence of discrete solutions that converges expo- nentially to the continuum solution. Asymptotically, we thus recover the property of the PML of having a zero reflection coefficient for all angles of incidence and all frequencies on the continuum level. This allows us to minimize the reflections from the discrete PML to an arbitrary level of accuracy while retaining optimal computational efficiency. Since our hp-adaptive scheme is automatic, no interaction with the user is required. This renders tedious parameter tuning of the PML obsolete. We demonstrate the improvement of the PML performance by hp-adaptivity through numerical results for acoustic, elastodynamic, and electromagnetic wave-propagation problems in the frequency domain and in different systems of coordinates. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 832–858, 2007 Keywords: Perfectly Matched Layer; hp-adaptivity; exterior boundary-value problem; acoustic scattering; elastodynamic wave propagation; electromagnetic scattering I. INTRODUCTION Most wave-propagation problems arising in acoustics, elastodynamics, and electromagnetics are posed on a spatially unbounded domain. Since domain-based discretization methods such as finite elements can handle only bounded domains, a truncation of the unbounded physical domain to a bounded computational domain is necessary.
    [Show full text]