Multiscale Modeling of Thin-Wire Coupling Problems Using Hybridization of Finite Element and Dipole Moment Methods and GPU Acceleration

Mittra, Raj
Kuzuoğlu, Mustafa
In this article, a hybrid numerical method, called finite element method (FEM) + dipole moment (DM), is presented for efficient solution of multiscale electromagnetic radiation and scattering problems that involve structures with fine features, such as thin-wire antennas or objects. In this method, the FEM is hybridized with the DM approach to help ease certain computational burdens, such as mesh refinement, ill-conditioning, memory overload, and long computation times, when solving multiscale problems with adequate numerical precision. Although the FEM+DM method is referred to as a single algorithm above, it is actually comprised of several different hybridization and/or implementation approaches. Two hybridization approaches are described, each tailored to the type of problem, e.g., radiation or scattering. Furthermore, both iterative and noniterative (self-consistent) implementations of the FEM+DM method are discussed. One of the important characteristics of the proposed method is that it can easily be parallelized to accelerate the speed of computation, thanks to the use of DMs for handling fine features as well as for modeling mutual interactions. The hybrid method with different hybridization/implementation approaches has been implemented in an in-house code by using graphics processing unit programming in MATLAB. Several numerical results are generated by the proposed method and are compared to those obtained from a commercial electromagnetic solver.


Numerical Solution of Multi-scale Electromagnetic Boundary Value Problems by Utilizing Transformation-Based Metamaterials
Ozgun, Ozlem; Kuzuoğlu, Mustafa (2011-06-23)
We present numerical solution techniques for efficiently handling multi-scale electromagnetic boundary value problems having fine geometrical details or features, by utilizing spatial coordinate transformations. The principle idea is to modify the computational domain of the finite methods (such as the finite element or finite difference methods) by suitably placing anisotropic metamaterial structures whose material parameters are obtained by coordinate transformations, and hence, to devise easier and effic...
Questioning Degree of Accuracy Offered by the Spectral Element Method in Computational Electromagnetics
Mahariq, I.; KURT, HAMZA; Kuzuoğlu, Mustafa (2015-07-01)
In this paper, a comparison amongst the spectral element method (SEM), the finite difference method (FDM), and the first-order finite element method (FEM) is presented. For the sake of consistency, the comparison is carried out on one-dimensional and two-dimensional boundary value problems based on the same measure of error in order to emphasize on the high accuracy gained by the SEM. Then, the deterioration in the accuracy of the SEM due to the elemental deformation is demonstrated. Following this, we try ...
Recursive Two-Way Parabolic Equation Approach for Modeling Terrain Effects in Tropospheric Propagation
Ozgun, Ozlem (2009-09-01)
The Fourier split-step method is a one-way marching-type algorithm to efficiently solve the parabolic equation for modeling electromagnetic propagation in troposphere. The main drawback of this method is that it characterizes only forward-propagating waves, and neglects backward-propagating waves, which become important especially in the presence of irregular surfaces. Although ground reflecting boundaries are inherently incorporated into the split-step algorithm, irregular surfaces (such as sharp edges) in...
Numerical Modeling of Electromagnetic Scattering from Periodic Structures by Transformation Electromagnetics
ÖZGÜN, ÖZLEM; Kuzuoğlu, Mustafa (2016-09-22)
The transformation electromagnetics is applied to the modeling of electromagnetic scattering from periodic structures in conjunction with the finite element method with periodic boundary conditions. In a unit cell of periodic structure, a uniform mesh is used over a flat surface and the arbitrary periodic surface is modeled by a coordinate transformation. The major advantage of this approach is that arbitrary geometries can be handled by using a single and simple mesh. Therefore, repeated computations (such...
Monte Carlo simulations of Helmholtz scattering from randomly positioned array of scatterers by utilizing coordinate transformations in finite element method
ÖZGÜN, ÖZLEM; Kuzuoğlu, Mustafa (2015-07-01)
Electromagnetic scattering from randomly distributed array of scatterers is numerically analyzed by Monte Carlo simulations by utilizing coordinate transformations in the context of finite element method solution of Helmholtz equation. The major goal in proposed approaches is to place transformation media into computational domain by employing the form invariance property of Maxwell's equations under coordinate transformations, and hence avoiding repeated mesh generation process in multiple realizations of ...
Citation Formats
Ö. ÖZGÜN, R. Mittra, and M. Kuzuoğlu, “Multiscale Modeling of Thin-Wire Coupling Problems Using Hybridization of Finite Element and Dipole Moment Methods and GPU Acceleration,” IEEE JOURNAL ON MULTISCALE AND MULTIPHYSICS COMPUTATIONAL TECHNIQUES, pp. 155–166, 2020, Accessed: 00, 2020. [Online]. Available: