The finite-volume time-domain algorithm using least square method in solving Maxwell's equations

ShiYan,LiangChang-Hong

IF: 4.645

2007-01-01

Journal of Computational Physics

Abstract:A finite-volume time-domain algorithm using least square method with a well-posed perfectly matched layer (PML) has been developed for the time-domain solution of Maxwellu0027s equations. This algorith...

What problem does this paper attempt to address?

The Finite-Volume Time-Domain Algorithm Using Least Square Method in Solving Maxwell’s Equations

Yan Shi,Chang-Hong Liang

DOI: https://doi.org/10.1016/j.jcp.2007.05.033

IF: 4.645

2007-01-01

Journal of Computational Physics

Abstract:A finite-volume time-domain algorithm using least square method with a well-posed perfectly matched layer (PML) has been developed for the time-domain solution of Maxwell’s equations. This algorithm uses the unstructured grids to obtain good computational efficiency and geometric flexibility. A novelty cell-wise data reconstruction scheme based on least square method is derived to achieve second-order spatial accuracy. A well-posed PML is applied to truncate computational domain by absorbing outgoing electromagnetic waves. The explicit Runge–Kutta scheme is employed to solve the semi-discrete Maxwell’s equations. Several numerical results are presented to illustrate the efficiency and accuracy of the algorithm.
Finite-Difference Time-Domain Simulation of the Maxwell–Schrödinger System

Christopher Ryu,Aiyin Liu,Weng C. Chew,Wei Sha,Lijun Jiang

DOI: https://doi.org/10.1109/aps.2016.7695823

2016-01-01

IEEE Journal on Multiscale and Multiphysics Computational Techniques

Abstract:A thorough study on the finite-difference time-domain (FDTD) simulation of the Maxwell-Schrödinger system in the semi-classical regime is given. For the Maxwell part which is treated classically, an alternative approach directly using the vector and scalar potentials (A and Φ) is taken. This approach is stable in the long wavelength regime and removes the need to extract the potentials at every time step. Perfectly matched layer (PML) is developed for this system. For the Schrödinger and quantum mechanical part, minimal coupling is applied to couple the charges to the electromagnetic potentials. Simulation results agree with the theory of quantum coherent states.
The Finite-Volume Time-Domain Electromagnetic Solver

Yan Shi,Chang-Hong Liang

DOI: https://doi.org/10.1109/icmmt.2007.381393

2007-01-01

Abstract:The paper proposes a finite-volume time-domain (FVTD) algorithm using least square (LS) method with a well-posed perfectly matched layer (PML). This approach uses the unstructured grids to attain computational efficiency and geometric flexibility. A novelty cell-wise data reconstruction scheme based on least square method is derived to achieve second-order spatial accuracy. Several numerical results are presented to illustrate the efficiency and accuracy of the algorithm.
A Novel High‐order Time‐domain Scheme for Three‐dimensional Maxwell's Equations

Zhi-Xiang Huang,Wei E. I. Sha,Xian-Liang Wu,Ming-Sheng Chen

DOI: https://doi.org/10.1002/mop.21563

IF: 1.311

2006-01-01

Microwave and Optical Technology Letters

Abstract:A novel high‐order time‐domain scheme with a four‐stage optimized symplectic integrator propagator is presented for 3D electromagnetic scattering problems. The scheme is nondissipative and does not require more storage than the classical finite‐difference time‐domain (FDTD) method. The numerical results show the scheme has better stability and more efficiency than the classical FDTD method. © 2006 Wiley Periodicals, Inc. Microwave Opt Technol Lett 48: 1123–1125, 2006; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/mop.21563
The High-Order Symplectic Finite-Difference Time-Domain Scheme

Wei E. I. Sha,Xian-Liang Wu,Zhi-Xiang Huang,Ming-Sheng Chen

DOI: https://doi.org/10.5772/9395

2010-01-01

Abstract:The book chapter will aim at introducing the background knowledge, basic theories, supporting techniques, numerical results, and future research for the high-order symplectic finite-difference time-domain scheme. The theories of symplectic geometry and Hamiltonian are reviewed in Section 2 followed by the symplectiness of Maxwell's equations presented in Section 3. Next, the numerical stability and dispersion analyses are given in Section 4. Then, in Section 5, we will make a tour of the supporting techniques but do not discuss them in detail. These techniques involve source excitation, perfectly matched layer, near-to-far-field transformation, inhomogeneous boundary treatments, and parameter extractions. The numerical results on propagation, scattering, and guided-wave problems are shown in Section 6. The high-order symplectic finite-difference time-domain scheme demonstrates the powerful advantages and potentials for the time-domain solution of Maxwell's equations, especially for electrically-large objects and for long-term simulation. Finally, the conclusion and future research are summarized in Section 7.
A Novel Matrix-Free Finite Element Method for Time-Harmonic Maxwell’s Equations

Yi-Yao Wang,Qiwei Zhan,Haoqiang Feng,Yin-Da Wang,Bozhao Sun,Wen-Yan Yin

DOI: https://doi.org/10.1109/tap.2024.3360700

IF: 5.7

2024-01-01

IEEE Transactions on Antennas and Propagation

Abstract:A frequency-domain finite element method (FEM) with the low-storage matrix-free feature is proposed for efficient analysis of electromagnetic fields. As opposed to conventional frequency-domain FEM which requires the matrix assembly before solving, the proposed matrix-free algorithm avoids the assembly and storage of the global matrix. In this approach, the global sparse matrix-vector (SpMV) multiplication is decomposed into element-wise matrix-vector (MV) multiplications. Additionally, the sum factorization technique is applied on tensorial basis functions to reduce the complexities of local MV multiplications. The numerical results demonstrate the superiority of our algorithm over traditional FEM in terms of both memory and time consumption. Besides, the improvement is more profound when higher-order basis functions are considered. A speedup of more than 10 against matrix-assembled solvers is observed. Given that memory transfer is bounded more than computation resources in modern supercomputer architectures, our algorithm is more friendly to high-performance computing platforms.
Decomposition Methods for Time‐domain Maxwell's Equations

Zhi-Xiang Huang,Wei E. I. Sha,Xian-Liang Wu,Ming-Sheng Chen

DOI: https://doi.org/10.1002/fld.1569

IF: 1.847

2008-01-01

International Journal for Numerical Methods in Fluids

Abstract:Decomposition methods based on split operators are proposed for numerical integration of the time‐domain Maxwell's equations for the first time. The methods are obtained by splitting the Hamiltonian function of Maxwell's equations into two analytically computable exponential sub‐propagators in the time direction based on different order decomposition methods, and then the equations are evaluated in the spatial direction by the staggered fourth‐order finite‐difference approximations. The stability and numerical dispersion analysis for different order decomposition methods are also presented. The theoretical predictions are confirmed by our numerical results. Copyright © 2007 John Wiley & Sons, Ltd.
Solution of Time-Domain Maxwell Equation with PML by Using Modified Laguerre Polynomials

LianLin Li,Fang Li

DOI: https://doi.org/10.1007/s11432-008-0044-y

2008-01-01

Abstract: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.
Solution of time-harmonic Maxwell's equations by a domain decomposition method based on PML transmission conditions

Sahar Borzooei,Victorita Dolean,Pierre-Henri Tournier,Claire Migliaccio

DOI: https://doi.org/10.48550/arXiv.2211.04912

2022-11-09

Numerical Analysis

Abstract:Numerical discretization of the large-scale Maxwell's equations leads to an ill-conditioned linear system that is challenging to solve. The key requirement for successive solutions of this linear system is to choose an efficient solver. In this work we use Perfectly Matched Layers (PML) to increase this efficiency. PML have been widely used to truncate numerical simulations of wave equations due to improving the accuracy of the solution instead of using absorbing boundary conditions (ABCs). Here, we will develop an efficient solver by providing an alternative use of PML as transmission conditions at the interfaces between subdomains in our domain decomposition method. We solve Maxwell's equations and assess the convergence rate of our solutions compared to the situation where absorbing boundary conditions are chosen as transmission conditions.
An Approach to Solving Maxwell's Equations in Time Domain

Hongli Yang,Xianyang Zeng,Xinyuan Wu

DOI: https://doi.org/10.1016/j.jmaa.2022.126678

IF: 1.417

2023-01-01

Journal of Mathematical Analysis and Applications

Abstract:In this paper, we propose a new direct approach to solving Maxwell's equations in time domain. Completely different from the classic one, this approach is derived and analysed based on a strategy which is independent of the well-known Fourier (or Laplace) transformation and the corresponding inverse transformation. The analysis first decouples Maxwell's equations on the basis of a potential theory and then uses an operator-variation-of-constants formula (the Duhamel Principle) to achieve the expression of exact solutions of Maxwell's equations. We further illustrate the simplicity and multi-purpose applicability of the new approach in linear, homogeneous, isotropic and lossless (bounded or unbounded and source-free or source) medium (in rectangular, cylindrical or spherical coordinates). Moreover, we make a study on impulse electromagnetics. With this new approach, this paper bridges the gap between the time domain electromagnetics and the frequency domain electromagnetics.
An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Qiwei Zhan,Yiyao Wang,Yuan Fang,Qiang Ren,Shiyou Yang,Wen-Yan Yin,Qing Huo Liu

DOI: https://doi.org/10.1109/tap.2021.3111639

IF: 5.7

2022-01-01

IEEE Transactions on Antennas and Propagation

Abstract:This article presents a stabilized nodal discontinuous Galerkin pseudospectral time-domain (DG-PSTD) algorithm for fully anisotropic electromagnetic waves. This solver permits arbitrary high-order basis functions and adaptive hexahedral elements, thus very efficient for large-scale wave propagation in complex media. Maxwell’s equations are reformulated in a unified hyperbolic form, where a localized anisotropic Riemann solver is derived to serve as a numerical flux to exchange information across adjacent elements in the DG-PSTD scheme. This local analysis method also helps impose the time-domain anisotropic plane wave incidence in the total/scattering field framework. Numerical validations and applications demonstrate the efficiency, accuracy, and capability of this new high-order solver for 3-D large-scale generally anisotropic electromagnetic media.
Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

Z. X. Huang,X. L. Wu,W. E. I. Sha,B. Wu

DOI: https://doi.org/10.1155/2012/956431

2012-01-01

International Journal of Antennas and Propagation

Abstract:Optimized operator splitting methods for numerical integration of the time domain Maxwell’s equations in computational electromagnetics (CEM) are proposed for the first time. The methods are based on splitting the time domain evolution operator of Maxwell’s equations into suboperators, and corresponding time coefficients are obtained by reducing the norm of truncation terms to a minimum. The general high-order staggered finite difference is introduced for discretizing the three-dimensional curl operator in the spatial domain. The detail of the schemes and explicit iterated formulas are also included. Furthermore, new high-order Padé approximations are adopted to improve the efficiency of the proposed methods. Theoretical proof of the stability is also included. Numerical results are presented to demonstrate the effectiveness and efficiency of the schemes. It is found that the optimized schemes with coarse discretized grid and large Courant-Friedrichs-Lewy (CFL) number can obtain satisfactory numerical results, which in turn proves to be a promising method, with advantages of high accuracy, low computational resources and facility of large domain and long-time simulation. In addition, due to the generality, our optimized schemes can be extended to other science and engineering areas directly.
An enhanced finite difference time domain method for two dimensional Maxwell's equations

Timothy Meagher,Bin Jiang,Peng Jiang

DOI: https://doi.org/10.1002/num.22467

2023-07-13

Abstract:An efficient finite-difference time-domain (FDTD) algorithm is built to solve the transverse electric 2D Maxwell's equations with inhomogeneous dielectric media where the electric fields are discontinuous across the dielectric interface. The new algorithm is derived based upon the integral version of the Maxwell's equations as well as the relationship between the electric fields across the interface. It is an improvement over the contour-path effective-permittivity algorithm by including some extra terms in the formulas. The scheme is validated in solving the scattering of a dielectric cylinder with exact solution from Mie theory and is then compared with the above contour-path method, the usual staircase and the volume-average method. The numerical results demonstrate that the new algorithm has achieved significant improvement in accuracy over the other methods. Furthermore, the algorithm has a simple structure and can be merged into any existing FDTD software package very easily.

Computational Physics,Numerical Analysis
A high-order finite volume method for Maxwell's equations in heterogeneous and time-varying media

Damian P. San Roman Alerigi,David I. Ketcheson,Boon S. Ooi

2023-07-22

Abstract:We develop a finite volume method for Maxwell's equations in materials whose electromagnetic properties vary in space and time. We investigate both conservative and non-conservative numerical formulations. High-order methods accurately resolve fine structures that develop due to the varying material properties. Numerical examples demonstrate the effectiveness of the proposed method in handling temporal variation and its efficiency relative to traditional 2nd-order FDTD.

Computational Physics,Numerical Analysis,Optics
Pml-Fdtd Method In Prolate Spheroidal Coordinates

Maoyu Zhang,Jianguo Wang

2010-01-01

Abstract:This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) formulation of Maxwell's equation and perfectly matched layers (PML) in 3-D prolate spheroidal coordinates. The perfectly matched layer formulation relies on the complex stretching coordinates approach. Numerical simulation result of the perfectly matched layer finite-difference time-domain (PML-FDTD) method shows very good agreement with reference solution computed on an extremely large space.
Implicit High-Order Unconditionally Stable Complex Envelope Algorithm for Solving the Time-Dependent Maxwell's Equations

Shuqi Chen,Weiping Zang,Axel Schuelzgen,Jinjie Liu,Lin Han,Yong Zeng,Jianguo Tian,Feng Song,Jerome V. Moloney,Nasser Peyghambarian

DOI: https://doi.org/10.1364/ol.33.002755

2008-01-01

Abstract:Based on the Padé approximation and multistep method, we propose an implicit high-order unconditionally stable complex envelope algorithm to solve the time-dependent Maxwell's equations. Unconditional numerical stability can be achieved simultaneously with a high-order accuracy in time. As we adopt the complex envelope Maxwell's equations, numerical dispersion and dissipation are very small even at comparatively large time steps. To verify the capability of our algorithm, we compare the results of the proposed method with the exact solutions.
Computation of electromagnetic scattering and radiation using a time-domain finite-volume discretization procedure

Alireza H. Mohammadian,Vijaya Shankar,William F. Hall

DOI: https://doi.org/10.1016/0010-4655(91)90199-u

IF: 4.717

1991-11-01

Computer Physics Communications

Abstract:Maxwell's equations in the time domain are first cast into a conservation form. They are then solved using a finite-volume discretization procedure proven very successful in solving some hyperbolic partial differential equations in computational fluid dynamics such as Euler or Navier-Stokes equations. The Lax-Wendroff explicit scheme is used to solve the discrete Maxwell's equations and as a result, second-order accuracy is achieved in both time and space. Multizoning is used to facilitate treatment of objects with internal and/or external sophistication. Body-fitted coordinates are used to map each zone independently from the rest of the object, allowing for nonmatching grid lines at the interface. Body-fitted coordinates facilitate accurate implementation of various boundary conditions on the object. The density of the computational mesh can vary and may be specified according to the local wavelength in the material. Also in the exterior region and away from the body, the mesh density is gradually decreased to reduce the total number of the grid points and the unknown field vectors. The formulation accounts for any variations in the material properties as a function of space, time, or frequency. In addition to perfect conductors, impedence surfaces and resistive sheets are handled. At the outer boundary of the computational domain a first-order, nonreflecting outer boundary condition is used. The Lax-Wendroff scheme, while providing second-order accuracy in the time and space, does not require, for its stability, an interlaced mesh for the electric and magnetic fields or the material properties. The method can be used for both time-harmonic and pulse excitations. At the end of the time-marching, the computed field in the time domain near the object and tangential to a conveniently chosen contour is stored. This information is converted into the frequency domain using fast Fourier transforms. If far field quantities such as radar cross-section or antenna radiation pattern are desired, they may be computed from the frequency-domain information using a standard free-space Green's function.

physics, mathematical,computer science, interdisciplinary applications
AN UNCONDITIONALLY STABLE PRECISE INTEGRATION TIME DOMAIN METHOD FOR THE NUMERICAL SOLUTION OF MAXWELL'S EQUATIONS IN CIRCULAR CYLINDRICAL COORDINATES

Xiaotian Zhao,Xikui Ma,Yanzhen Zhao

DOI: https://doi.org/10.2528/pier06120505

2007-01-01

Abstract:The extension of an unconditionally stable precise integration time domain method for the numerical solutions of Maxwell's equations to circular cylindrical coordinate system is presented in this paper.In contrast with the conventional cylindrical finite-difference time-domain method, not only can it remove the Courant stability condition constraint, but also make the numerical dispersion independent of the time-step size.Moreover, the first-order absorbing boundary condition can be introduced into the proposed method successfully, whereas the alternating-direction-implicit finitedifference time-domain method may become instable for open region radiation problems terminated with absorbing boundary conditions.Theoretical proof of the unconditional stability is mentioned and the numerical results are presented to demonstrate the effectiveness of the proposed method in solving electromagnetic-field problem.
Application of the RBF-based Meshless Method to Solve 2-D Time Domain Maxwell&#x2019;s Equations

Sheng-Jian Lai,Bing‐Zhong Wang,Yong Duan

DOI: https://doi.org/10.1109/icmmt.2008.4540505

2008-01-01

Abstract:Radial basis functions (RBF), as a meshless technique, is widely applied to solve partial differential equations. In this paper, a meshless RBF method is applied to time domain Maxwell's equations and calculates a two-dimensional (2-D) cavity case. The main idea is that the fields in the space domain are expanded into a series of radial basis functions and are treated with a meshless method procedure, and the time derivatives are still tackled with the customary difference scheme. The 2-D cavity numerical experiment has been used to validate the propose technique.
A 3-D Precise Integration Time-Domain Method Without the Restraints of the Courant-Friedrich-levy Stability Condition for the Numerical Solution of Maxwell's Equations

Xikui Ma,Xintai Zhao,Yanzhen Zhao

DOI: https://doi.org/10.1109/tmtt.2006.877427

IF: 4.3

2006-01-01

IEEE Transactions on Microwave Theory and Techniques

Abstract:In this paper, a new three-dimensional time-domain method for solving vector Maxwell's equations, called the precise-integration time-domain (PITD) algorithm, is proposed in order to eliminate the Courant-Friedrich-Levy (CFL) condition restraint. The new algorithm is based on the precise-integration technique. It is shown that this method is quite stable even when the CFL condition is not satisfied. Although the memory requirement of the PITD method is much larger than that of the finite-difference time-domain (FDTD) method, this new algorithm is very appealing since the time step used in the simulation is no longer restricted by stability. As a result, computation speed can be improved. Therefore, if the minimum cell size in the computational domain is required to be much smaller than the wavelength, this new algorithm is more efficient than the FDTD scheme. Theoretical proof of the unconditional stability is shown and numerical results are presented to demonstrate the effectiveness and efficiency of the method. It is found that the accuracy of the PITD is independent of the time-step size.

The finite-volume time-domain algorithm using least square method in solving Maxwell's equations

The Finite-Volume Time-Domain Algorithm Using Least Square Method in Solving Maxwell’s Equations

Finite-Difference Time-Domain Simulation of the Maxwell–Schrödinger System

The Finite-Volume Time-Domain Electromagnetic Solver

A Novel High‐order Time‐domain Scheme for Three‐dimensional Maxwell's Equations

The High-Order Symplectic Finite-Difference Time-Domain Scheme

A Novel Matrix-Free Finite Element Method for Time-Harmonic Maxwell’s Equations

Decomposition Methods for Time‐domain Maxwell's Equations

Solution of Time-Domain Maxwell Equation with PML by Using Modified Laguerre Polynomials

Solution of time-harmonic Maxwell's equations by a domain decomposition method based on PML transmission conditions

An Approach to Solving Maxwell's Equations in Time Domain

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

An enhanced finite difference time domain method for two dimensional Maxwell's equations

A high-order finite volume method for Maxwell's equations in heterogeneous and time-varying media

Pml-Fdtd Method In Prolate Spheroidal Coordinates

Implicit High-Order Unconditionally Stable Complex Envelope Algorithm for Solving the Time-Dependent Maxwell's Equations

Computation of electromagnetic scattering and radiation using a time-domain finite-volume discretization procedure

AN UNCONDITIONALLY STABLE PRECISE INTEGRATION TIME DOMAIN METHOD FOR THE NUMERICAL SOLUTION OF MAXWELL'S EQUATIONS IN CIRCULAR CYLINDRICAL COORDINATES

Application of the RBF-based Meshless Method to Solve 2-D Time Domain Maxwell&amp;#x2019;s Equations

A 3-D Precise Integration Time-Domain Method Without the Restraints of the Courant-Friedrich-levy Stability Condition for the Numerical Solution of Maxwell's Equations

Application of the RBF-based Meshless Method to Solve 2-D Time Domain Maxwell’s Equations