-
Sub-structure characteristic mode analysis of microstrip antennas using a global multi-trace formulation
Authors:
Ran Zhao,
Yuyu Lu,
Guang Shang Cheng,
Wei Zhu,
Jun Hu,
Hakan Bagci
Abstract:
A characteristic mode (CM) method that relies on a global multi-trace formulation (MTF) of surface integral equations is proposed to compute the modes and the resonance frequencies of microstrip patch antennas with finite dielectric substrates and ground planes. Compared to the coupled formulation of electric field and Poggio-Miller-Chang-Harrington-Wu-Tsai integral equations, global MTF allows fo…
▽ More
A characteristic mode (CM) method that relies on a global multi-trace formulation (MTF) of surface integral equations is proposed to compute the modes and the resonance frequencies of microstrip patch antennas with finite dielectric substrates and ground planes. Compared to the coupled formulation of electric field and Poggio-Miller-Chang-Harrington-Wu-Tsai integral equations, global MTF allows for more direct implementation of a sub-structure CM method. This is achieved by representing the coupling of the electromagnetic fields on the substrate and ground plane in the form of a numerical Green function matrix, which yields a more compact generalized eigenvalue equation. The resulting sub-structure CM method avoids the cumbersome computation of the multilayered medium Green function (unlike the CM methods that rely on mixed-potential integral equations) and the volumetric discretization of the substrate (unlike the CM methods that rely on volume-surface integral equations), and numerical results show that it is a reliable and accurate approach to predicting the modal behavior of electromagnetic fields on practical microstrip antennas.
△ Less
Submitted 23 October, 2023; v1 submitted 22 December, 2022;
originally announced December 2022.
-
A Coupled Hybridizable Discontinuous Galerkin and Boundary Integral Method for Analyzing Electromagnetic Scattering
Authors:
Ran Zhao,
Ming Dong,
Liang Chen,
Jun Hu,
Hakan Bagci
Abstract:
A coupled hybridizable discontinuous Galerkin (HDG) and boundary integral (BI) method is proposed to efficiently analyze electromagnetic scattering from inhomogeneous/composite objects. The coupling between the HDG and the BI equations is realized using the numerical flux operating on the equivalent current and the global unknown of the HDG. This approach yields sparse coupling matrices upon discr…
▽ More
A coupled hybridizable discontinuous Galerkin (HDG) and boundary integral (BI) method is proposed to efficiently analyze electromagnetic scattering from inhomogeneous/composite objects. The coupling between the HDG and the BI equations is realized using the numerical flux operating on the equivalent current and the global unknown of the HDG. This approach yields sparse coupling matrices upon discretization. Inclusion of the BI equation ensures that the only error in enforcing the radiation conditions is the discretization. However, the discretization of this equation yields a dense matrix, which prohibits the use of a direct matrix solver on the overall coupled system as often done with traditional HDG schemes. To overcome this bottleneck, a "hybrid" method is developed. This method uses an iterative scheme to solve the overall coupled system but within the matrix-vector multiplication subroutine of the iterations, the inverse of the HDG matrix is efficiently accounted for using a sparse direct matrix solver. The same subroutine also uses the multilevel fast multipole algorithm to accelerate the multiplication of the guess vector with the dense BI matrix. The numerical results demonstrate the accuracy, the efficiency, and the applicability of the proposed HDG-BI solver.
△ Less
Submitted 19 June, 2023; v1 submitted 6 October, 2022;
originally announced October 2022.
-
A Time Domain Volume Integral Equation Solver to Analyze Electromagnetic Scattering from Nonlinear Dielectric Objects
Authors:
Sadeed Bin Sayed,
Rui Chen,
Huseyin Arda Ulku,
Hakan Bagci
Abstract:
A time domain electric field volume integral equation (TD-EFVIE) solver is proposed for analyzing electromagnetic scattering from dielectric objects with Kerr nonlinearity. The nonlinear constitutive relation that relates electric flux and electric field induced in the scatterer is used as an auxiliary equation that complements TD-EFVIE. The ordinary differential equation system that arises from T…
▽ More
A time domain electric field volume integral equation (TD-EFVIE) solver is proposed for analyzing electromagnetic scattering from dielectric objects with Kerr nonlinearity. The nonlinear constitutive relation that relates electric flux and electric field induced in the scatterer is used as an auxiliary equation that complements TD-EFVIE. The ordinary differential equation system that arises from TD-EFVIE's Schaubert-Wilton-Glisson (SWG)-based discretization is integrated in time using a predictor-corrector method for the unknown expansion coefficients of the electric field. Matrix systems that arise from the SWG-based discretization of the nonlinear constitutive relation and its inverse obtained using the Pade approximant are used to carry out explicit updates of the electric field and the electric flux expansion coefficients at the predictor and the corrector stages of the time integration method. The resulting explicit marching-on-in-time (MOT) scheme does not call for any Newton-like nonlinear solver and only requires solution of sparse and well-conditioned Gram matrix systems at every step. Numerical results show that the proposed explicit MOT-based TD-EFVIE solver is more accurate than the finite-difference time-domain method that is traditionally used for analyzing transient electromagnetic scattering from nonlinear objects.
△ Less
Submitted 2 February, 2023; v1 submitted 17 July, 2022;
originally announced August 2022.
-
Calculation of Photocarrier Generation from Optical Absorption for Time-domain Simulation of Optoelectronic Devices
Authors:
Liang Chen,
Hakan Bagci
Abstract:
Photocarrier generation rate in optoelectronic materials is often calculated using the Poynting vector in the frequency domain. However, this approach is not accurate in time-domain simulations of photoconductive devices because the instantaneous Poynting vector does not distinguish between power flux densities of optical and low-frequency electromagnetic fields. The latter is generated by photocu…
▽ More
Photocarrier generation rate in optoelectronic materials is often calculated using the Poynting vector in the frequency domain. However, this approach is not accurate in time-domain simulations of photoconductive devices because the instantaneous Poynting vector does not distinguish between power flux densities of optical and low-frequency electromagnetic fields. The latter is generated by photocurrents and is not supposed to contribute to the photocarrier generation since the corresponding photon energy is smaller than the bandgap energy of the optoelectronic material. In this work, an optical absorption-based model is proposed to accurately calculate the generation rate in time-domain simulations. The proposed approach considers the material dispersion near the optical frequency corresponding to the bandgap energy of the optoelectronic material. The instantaneous optical absorption is calculated from the polarization current density associated with the dispersion model. Then, the optical absorption is used to calculate the generation rate. Numerical examples show that the proposed approach is more accurate than the Poynting vector-based method in calculating the instantaneous optical absorption. The proposed method is further validated against experimental results by modeling a photoconductive device. In the multiphysics simulation, the Poynting vector-based method overestimates the carrier generation rate and even generates divergent carrier densities when the low-frequency fields are strong, while the proposed method produces results that match with experimental measurements well.
△ Less
Submitted 7 October, 2024; v1 submitted 12 February, 2021;
originally announced February 2021.
-
A Memory-efficient Implementation of Perfectly Matched Layer with Smoothly-varying Coefficients in Discontinuous Galerkin Time-Domain Method
Authors:
Liang Chen,
Mehmet Burak Ozakin,
Shehab Ahmed,
Hakan Bagci
Abstract:
Wrapping a computation domain with a perfectly matched layer (PML) is one of the most effective methods of imitating/approximating the radiation boundary condition in Maxwell and wave equation solvers. Many PML implementations often use a smoothly-increasing attenuation coefficient to increase the absorption for a given layer thickness, and, at the same time, to reduce the numerical reflection fro…
▽ More
Wrapping a computation domain with a perfectly matched layer (PML) is one of the most effective methods of imitating/approximating the radiation boundary condition in Maxwell and wave equation solvers. Many PML implementations often use a smoothly-increasing attenuation coefficient to increase the absorption for a given layer thickness, and, at the same time, to reduce the numerical reflection from the interface between the computation domain and the PML. In discontinuous Galerkin time-domain (DGTD) methods, using a PML coefficient that varies within a mesh element requires a different mass matrix to be stored for every element and therefore significantly increases the memory footprint. In this work, this bottleneck is addressed by applying a weight-adjusted approximation to these mass matrices. The resulting DGTD scheme has the same advantages as the scheme that stores individual mass matrices, namely higher accuracy (due to reduced numerical reflection) and increased meshing flexibility (since the PML does not have to be defined layer by layer) but it requires significantly less memory.
△ Less
Submitted 18 July, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.
-
A hybridizable discontinuous Galerkin method for simulation of electrostatic problems with floating potential conductors
Authors:
Liang Chen,
Ming Dong,
Ping Li,
Hakan Bagci
Abstract:
In an electrostatic simulation, an equipotential condition with an undefined/floating potential value has to be enforced on the surface of an isolated conductor. If this conductor is charged, a nonzero charge condition is also required. While implementation of these conditions using a traditional finite element method (FEM) is not straightforward, they can be easily discretized and incorporated wi…
▽ More
In an electrostatic simulation, an equipotential condition with an undefined/floating potential value has to be enforced on the surface of an isolated conductor. If this conductor is charged, a nonzero charge condition is also required. While implementation of these conditions using a traditional finite element method (FEM) is not straightforward, they can be easily discretized and incorporated within a discontinuous Galerkin (DG) method. However, DG discretization results in a larger number of unknowns as compared to FEM. In this work, a hybridizable DG (HDG) method is proposed to alleviate this problem. Floating potential boundary conditions, possibly with different charge values, are introduced on surfaces of each isolated conductor and are weakly enforced in the global problem of HDG. The unknowns of the global HDG problem are those only associated with the nodes on the mesh skeleton and their number is much smaller than the total number of unknowns required by DG. Numerical examples show that the proposed method is as accurate as DG while it improves the computational efficiency significantly.
△ Less
Submitted 25 September, 2020; v1 submitted 29 May, 2020;
originally announced June 2020.
-
Efficient Discontinuous Galerkin Scheme for Analyzing Nanostructured Photoconductive Devices
Authors:
Liang Chen,
Kostyantyn Sirenko,
Ping Li,
Hakan Bagci
Abstract:
Incorporation of plasmonic nanostructures in the design of photoconductive devices (PCDs) has significantly improved their optical-to-terahertz conversion efficiency. However, this improvement comes at the cost of increased complexity for the design and simulation of these devices. Indeed, accurate and efficient modeling of multiphysics processes and intricate device geometries of nanostructured P…
▽ More
Incorporation of plasmonic nanostructures in the design of photoconductive devices (PCDs) has significantly improved their optical-to-terahertz conversion efficiency. However, this improvement comes at the cost of increased complexity for the design and simulation of these devices. Indeed, accurate and efficient modeling of multiphysics processes and intricate device geometries of nanostructured PCDs is challenging due to the high computational cost resulting from multiple characteristic scales in time and space. In this work, a discontinuous Galerkin (DG)-based unit-cell scheme for efficient simulation of PCDs with periodic nanostructures is proposed. The scheme considers two physical stages of the device and models them using two coupled systems: a system of Poisson and drift-diffusion equations describing the nonequilibrium steady state, and a system of Maxwell and drift-diffusion equations describing the transient stage. A "potential-drop" boundary condition is enforced on the opposing boundaries of the unit cell to mimic the effect of the bias voltage. Periodic boundary conditions are used for carrier densities and electromagnetic fields. The unit-cell model described by these coupled equations and boundary conditions is discretized using DG methods. Numerical results demonstrate that the proposed DG-based unit-cell scheme has the same accuracy in predicting the THz photocurrent as the DG framework that takes into account the whole device, while it significantly reduces the computational cost.
△ Less
Submitted 14 April, 2021; v1 submitted 29 May, 2020;
originally announced June 2020.
-
Multiphysics Simulation of Plasmonic Photoconductive Devices using Discontinuous Galerkin Methods
Authors:
Liang Chen,
Hakan Bagci
Abstract:
Plasmonic nanostructures significantly improve the performance of photoconductive devices (PCDs) in generating terahertz radiation. However, they are geometrically intricate and result in complicated electromagnetic (EM) field and carrier interactions under a bias voltage and upon excitation by an optical EM wave. These lead to new challenges in simulations of plasmonic PCDs, which cannot be addre…
▽ More
Plasmonic nanostructures significantly improve the performance of photoconductive devices (PCDs) in generating terahertz radiation. However, they are geometrically intricate and result in complicated electromagnetic (EM) field and carrier interactions under a bias voltage and upon excitation by an optical EM wave. These lead to new challenges in simulations of plasmonic PCDs, which cannot be addressed by existing numerical frameworks. In this work, a multiphysics framework making use of discontinuous Galerkin (DG) methods is developed to address these challenges. The operation of the PCD is analyzed in stationary and transient states, which are described by coupled systems of the Poisson and stationary drift-diffusion (DD) equations and the time-dependent Maxwell and DD equations, respectively. Both systems are discretized using DG schemes. The nonlinearity of the stationary system is accounted for using the Gummel iterative method while the nonlinear coupling between the time-dependent Maxwell and DD equations is tackled during time integration. The DG-based discretization and the explicit time marching help in handling space and time characteristic scales that are associated with different physical processes and differ by several orders of magnitude. The accuracy and applicability of the resulting multiphysics framework are demonstrated via simulations of conventional and plasmonic PCDs.
△ Less
Submitted 12 September, 2020; v1 submitted 8 December, 2019;
originally announced December 2019.
-
Steady-state Simulation of Semiconductor Devices using Discontinuous Galerkin Methods
Authors:
Liang Chen,
Hakan Bagci
Abstract:
Design of modern nanostructured semiconductor devices often calls for simulation tools capable of modeling arbitrarily-shaped multiscale geometries. In this work, to this end, a discontinuous Galerkin (DG) method-based framework is developed to simulate steady-state response of semiconductor devices. The proposed framework solves a system of Poisson equation (in electric potential) and drift-diffu…
▽ More
Design of modern nanostructured semiconductor devices often calls for simulation tools capable of modeling arbitrarily-shaped multiscale geometries. In this work, to this end, a discontinuous Galerkin (DG) method-based framework is developed to simulate steady-state response of semiconductor devices. The proposed framework solves a system of Poisson equation (in electric potential) and drift-diffusion equations (in charge densities), which are nonlinearly coupled via the drift current and the charge distribution. This system is decoupled and linearized using the Gummel method and the resulting equations are discretized using a local DG scheme. The proposed framework is used to simulate geometrically intricate semiconductor devices with realistic models of mobility and recombination rate. Its accuracy is demonstrated by comparing the results to those obtained by the finite volume and finite element methods implemented in a commercial software package.
△ Less
Submitted 14 October, 2019;
originally announced October 2019.
-
Extreme Scale FMM-Accelerated Boundary Integral Equation Solver for Wave Scattering
Authors:
Mustafa Abduljabbar,
Mohammed Al Farhan,
Noha Al-Harthi,
Rui Chen,
Rio Yokota,
Hakan Bagci,
David Keyes
Abstract:
Algorithmic and architecture-oriented optimizations are essential for achieving performance worthy of anticipated energy-austere exascale systems. In this paper, we present an extreme scale FMM-accelerated boundary integral equation solver for wave scattering, which uses FMM as a matrix-vector multiplication inside the GMRES iterative method. Our FMM Helmholtz kernels treat nontrivial singular and…
▽ More
Algorithmic and architecture-oriented optimizations are essential for achieving performance worthy of anticipated energy-austere exascale systems. In this paper, we present an extreme scale FMM-accelerated boundary integral equation solver for wave scattering, which uses FMM as a matrix-vector multiplication inside the GMRES iterative method. Our FMM Helmholtz kernels treat nontrivial singular and near-field integration points. We implement highly optimized kernels for both shared and distributed memory, targeting emerging Intel extreme performance HPC architectures. We extract the potential thread- and data-level parallelism of the key Helmholtz kernels of FMM. Our application code is well optimized to exploit the AVX-512 SIMD units of Intel Skylake and Knights Landing architectures. We provide different performance models for tuning the task-based tree traversal implementation of FMM, and develop optimal architecture-specific and algorithm aware partitioning, load balancing, and communication reducing mechanisms to scale up to 6,144 compute nodes of a Cray XC40 with 196,608 hardware cores. With shared memory optimizations, we achieve roughly 77% of peak single precision floating point performance of a 56-core Skylake processor, and on average 60% of peak single precision floating point performance of a 72-core KNL. These numbers represent nearly 5.4x and 10x speedup on Skylake and KNL, respectively, compared to the baseline scalar code. With distributed memory optimizations, on the other hand, we report near-optimal efficiency in the weak scalability study with respect to both the logarithmic communication complexity as well as the theoretical scaling complexity of FMM. In addition, we exhibit up to 85% efficiency in strong scaling. We compute in excess of 2 billion DoF on the full-scale of the Cray XC40 supercomputer.
△ Less
Submitted 27 March, 2018;
originally announced March 2018.