Browsing by Subject "Surface integral equations"

Now showing 1 - 19 of 19

Open Access
Accurate solutions of scattering problems involving low-contrast dielectric objects with surface integral equations
(Institution of Engineering and Technology, 2007) Ergül, Özgür; Gürel, Levent
We present the stabilization of the surface integral equations for accurate solutions of scattering problems involving low-contrast dielectric objects. Unlike volume formulations, conventional surface formulations fail to provide accurate results for the scattered fields when the contrast of the object is small. Therefore, surface formulations are required to be stabilized by extracting the nonradiating parts of the equivalent currents. In addition to previous strategies for the stabilization, we introduce a novel procedure called field-based stabilization (FBS) based on using fictitious incident fields and rearranging the right-hand-side of the equations. The results show that the formulations using FBS provide accurate results even for scattering problems involving extremely low-contrast objects, while the extra cost due to the stabilization procedure is negligible.
Open Access
Analysis of double-negative materials with surface integral equations and the multilevel fast multipole algorithm
(IEEE, 2011) Ergül O.; Gürel, Levent
We present a fast and accurate analysis of double-negative materials (DNMs) with surface integral equations and the multilevel fast multipole algorithm (MLFMA). DNMs are commonly used as simplified models of metamaterials at resonance frequencies and are suitable to be formulated with surface integral equations. However, realistic metamaterials and their models are usually very large with respect to wavelength and their accurate solutions require fast algorithms, such as MLFMA. We consider iterative solutions of DNMs with MLFMA and we investigate the accuracy and efficiency of solutions when DNMs are formulated with two recently developed formulations, namely, the combined tangential formulation (CTF) and the electric and magnetic current combined-field integral equation (JMCFIE). Numerical results on canonical objects are consistent with previous results in the literature on ordinary objects. © 2011 IEEE.
Open Access
Approximate Schur preconditioners for efficient solutions of dielectric problems formulated with surface integral equations
(IEEE, 2009-07) Malas, Tahir; Gürel, Levent
We propose direct and iterative versions of approximate Schur preconditioners to increase robustness and efficiency of iterative solutions of dielectric problems formulated with surface integral equations. The performance of these preconditioners depends on the availability of fast and approximate solutions to reduced matrix systems. We show that sparse-approximate-inverse techniques provide a suitable mechanism for this purpose. The proposed preconditioners are demonstrated to significantly improve convergence rates of dielectric problems formulated with two different surface integral equations. ©2009 IEEE.
Open Access
Comparison of integral-equation formulations for the fast and accurate solution of scattering problems involving dielectric objects with the multilevel fast multipole algorithm
(Institute of Electrical and Electronics Engineers, 2009) Ergül, Özgür; Gürel, Levent
We consider fast and accurate solutions of scattering problems involving increasingly large dielectric objects formulated by surface integral equations. We compare various formulations when the objects are discretized with Rao-Wilton-Glisson functions, and the resulting matrix equations are solved iteratively by employing the multilevel fast multipole algorithm (MLFMA). For large problems, we show that a combined-field formulation, namely, the electric and magnetic current combined-field integral equation (JMCFIE), requires fewer iterations than other formulations within the context of MLFMA. In addition to its efficiency, JMCFIE is also more accurate than the normal formulations and becomes preferable, especially when the problems cannot be solved easily with the tangential formulations.
Open Access
Contamination of the accuracy of the combined-field integral equation with the discretization error of the magnetic-field integral equation
(Institute of Electrical and Electronics Engineers, 2009) Gürel, Levent; Ergül, Özgür
We investigate the accuracy of the combined-field integral equation (CFIE) discretized with the Rao-Wilton-Glisson (RWG) basis functions for the solution of scattering and radiation problems involving three-dimensional conducting objects. Such a low-order discretization with the RWG functions renders the two components of CFIE, i.e., the electric-field integral equation (EFIE) and the magnetic-field integral equation (MFIE), incompatible, mainly because of the excessive discretization error of MFIE. Solutions obtained with CFIE are contaminated with the MFIE inaccuracy, and CFIE is also incompatible with EFIE and MFIE. We show that, in an iterative solution, the minimization of the residual error for CFIE involves a breakpoint, where a further reduction of the residual error does not improve the solution in terms of compatibility with EFIE, which provides a more accurate reference solution. This breakpoint corresponds to the last useful iteration, where the accuracy of CFIE is saturated and a further reduction of the residual error is practically unnecessary.
Open Access
Discretization error due to the identity operator in surface integral equations
(ELSEVIER, 2009-05-03) Ergül, Özgür; Gürel, Levent
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.
Open Access
Fast and accurate analysis of large-scale composite structures with the parallel multilevel fast multipole algorithm
(Optical Society of America, 2013) Ergül, Özgür; Gürel, Levent
Accurate electromagnetic modeling of complicated optical structures poses several challenges. Optical metamaterial and plasmonic structures are composed of multiple coexisting dielectric and/or conducting parts. Such composite structures may possess diverse values of conductivities and dielectric constants, including negative permittivity and permeability. Further challenges are the large sizes of the structures with respect to wavelength and the complexities of the geometries. In order to overcome these challenges and to achieve rigorous and efficient electromagnetic modeling of three-dimensional optical composite structures, we have developed a parallel implementation of the multilevel fast multipole algorithm (MLFMA). Precise formulation of composite structures is achieved with the so-called "electric and magnetic current combined-field integral equation." Surface integral equations are carefully discretized with piecewise linear basis functions, and the ensuing dense matrix equations are solved iteratively with parallel MLFMA. The hierarchical strategy is used for the efficient parallelization of MLFMA on distributed-memory architectures. In this paper, fast and accurate solutions of large-scale canonical and complicated real-life problems, such as optical metamaterials, discretized with tens of millions of unknowns are presented in order to demonstrate the capabilities of the proposed electromagnetic solver.
Open Access
Fast and accurate analysis of optical metamaterials using surface integral equations and the parallel multilevel fast multipole algorithm
(IEEE, 2013) Ergül, Özgür; Gürel, Levent
We present fast and accurate simulations of optical metamaterials using surface integral equations and the multilevel fast multipole algorithm (MLFMA). Problems are formulated with the electric and magnetic current combined-field integral equation and solved iteratively with MLFMA, which is parallelized using the hierarchical strategy on distributed-memory architectures. Realistic metamaterials involving dielectric, perfectly conducting, and plasmonic regions of finite extents are solved rigorously with the developed implementation without any periodicity assumptions.
Open Access
Fast and accurate solutions of extremely large integral-equation problems discretised with tens of millions of unknowns
(The Institution of Engineering and Technology, 2007) Gürel, Levent; Ergül, Özgür
The solution of extremely large scattering problems that are formulated by integral equations and discretised with tens of millions of unknowns is reported. Accurate and efficient solutions are performed by employing a parallel implementation of the multilevel fast multipole algorithm. The effectiveness of the implementation is demonstrated on a sphere problem containing more than 33 million unknowns, which is the largest integral-equation problem ever solved to our knowledge.
Open Access
Fast and accurate solutions of scattering problems involving dielectric objects with moderate and low contrasts
(IEEE, 2007-08) Ergül, Özgür; Gürel, Levent
We consider the solution of electromagnetic scattering problems involving relatively large dielectric objects with moderate and low contrasts. Three-dimensional objects are discretized with Rao-Wilton-Glisson functions and the scattering problems are formulated with surface integral equations. The resulting dense matrix equations are solved iteratively by employing the multilevel fast multipole algorithm. We compare the accuracy and efficiency of the results obtained by employing various integral equations for the formulation of the problem. If the problem size is large, we show that a combined formulation, namely, electric-magnetic current combined-field integral equation, provides faster iterative convergence compared to other formulations, when it is accelerated with an efficient block preconditioner. For low-contrast problems, we introduce various stabilization procedures in order to avoid the numerical breakdown encountered in the conventional surface formulations. © 2007 IEEE.
Open Access
Hierarchical parallelization of the multilevel fast multipole algorithm (MLFMA)
(IEEE, 2013) Gürel, Levent; Ergül, Özgür
Due to its O(N log N) complexity, the multilevel fast multipole algorithm (MLFMA) is one of the most prized algorithms of computational electromagnetics and certain other disciplines. Various implementations of this algorithm have been used for rigorous solutions of large-scale scattering, radiation, and miscellaneous other electromagnetics problems involving 3-D objects with arbitrary geometries. Parallelization of MLFMA is crucial for solving real-life problems discretized with hundreds of millions of unknowns. This paper presents the hierarchical partitioning strategy, which provides a very efficient parallelization of MLFMA on distributed-memory architectures. We discuss the advantages of the hierarchical strategy over previous approaches and demonstrate the improved efficiency on scattering problems discretized with millions of unknowns. © 1963-2012 IEEE.
Open Access
A hierarchical partitioning strategy for an efficient parallelization of the multilevel fast multipole algorithm
(IEEE, 2009) Ergül, Özgür; Gürel, Levent
We present a novel hierarchical partitioning strategy for the efficient parallelization of the multilevel fast multipole algorithm (MLFMA) on distributed-memory architectures to solve large-scale problems in electromagnetics. Unlike previous parallelization techniques, the tree structure of MLFMA is distributed among processors by partitioning both clusters and samples of fields at each level. Due to the improved load-balancing, the hierarchical strategy offers a higher parallelization efficiency than previous approaches, especially when the number of processors is large. We demonstrate the improved efficiency on scattering problems discretized with millions of unknowns. In addition, we present the effectiveness of our algorithm by solving very large scattering problems involving a conducting sphere of radius 210 wavelengths and a complicated real-life target with a maximum dimension of 880 wavelengths. Both of the objects are discretized with more than 200 million unknowns.
Open Access
Iterative solutions of hybrid integral equations for coexisting open and closed surfaces
(IEEE, 2009) Ergül, Özgür; Gürel, Levent
We consider electromagnetics problems involving composite geometries with coexisting open and closed conductors. Hybrid integral equations are presented to improve the efficiency of the solutions, compared to the conventional electric-field integral equation. We investigate the convergence characteristics of iterative solutions of large composite problems with the multilevel fast multipole algorithm. Following a thorough study of how the convergence characteristics depends on the problem geometry, formulation, and iterative solvers, we provide concrete guidelines for efficient solutions.
Open Access
Novel electromagnetic surface integral equations for highly accurate computations of dielectric bodies with arbitrarily low contrasts
(Journal of Computational Physics, 2008) Ergül O.; Gürel, Levent
We present a novel stabilization procedure for accurate surface formulations of electromagnetic scattering problems involving three-dimensional dielectric objects with arbitrarily low contrasts. Conventional surface integral equations provide inaccurate results for the scattered fields when the contrast of the object is low, i.e., when the electromagnetic material parameters of the scatterer and the host medium are close to each other. We propose a stabilization procedure involving the extraction of nonradiating currents and rearrangement of the right-hand side of the equations using fictitious incident fields. Then, only the radiating currents are solved to calculate the scattered fields accurately. This technique can easily be applied to the existing implementations of conventional formulations, it requires negligible extra computational cost, and it is also appropriate for the solution of large problems with the multilevel fast multipole algorithm. We show that the stabilization leads to robust formulations that are valid even for the solutions of extremely low-contrast objects. © 2008 Elsevier Inc. All rights reserved.
Open Access
Rigorous analysis of double-negative materials with the multilevel fast multipole algorithm
(Applied Computational Electromagnetics Society, Inc., 2012) Ergül, Özgür; Gürel, Levent
We present rigorous analysis of double-negative materials (DNMs) with surface integral equations and the multilevel fast multipole algorithm (MLFMA). Accuracy and efficiency of numerical solutions are investigated when DNMs are formulated with two recently developed formulations, i.e., the combined tangential formulation (CTF) and the electric and magnetic current combined-field integral equation (JMCFIE). Simulation results on canonical objects are consistent with previous results in the literature on ordinary objects. MLFMA is also parallelized to solve extremely large electromagnetics problems involving DNMs.
Open Access
Scattering analysis of silver nanoparticles for solar cell applications using integral equations
(IEEE, 2018) Uysal, İsmail Enes; Ülkü, H. A.; Bağcı, H.; Gülseren, Oğuz
Plasmonic nanoparticles (NPs) can be used to improve the efficiency of solar cells. Analysis of electromagnetic scattering from NPs is often carried out under the assumptions that they reside in air and have 'ideal' shapes (sphere, cube, etc.) However, in a realistic setup, nanoparticles are fabricated on a substrate and their shape and size cannot be controlled precisely. In this work, a surface integral equation solver is used to accurately characterize the scattering from a realistic system, where silver hemispheres of varying sizes are fabricated on an indium tin-oxide substrate. Results obtained by the solver are compared to the experimental results obtained for a similar system.
Open Access
Solutions of large integral-equation problems with preconditioned MLFMA
(IEEE, 2007) Ergül, Özgür; Malas, Tahir; Ünal, Alper; Gürel, Levent
We report the solution of the largest integral-equation problems in computational electromagnetics. We consider matrix equations obtained from the discretization of the integral-equation formulations that are solved iteratively by employing parallel multilevel fast multipole algorithm (MLFMA). With the efficient parallelization of MLFMA, scattering and radiation problems with millions of unknowns are easily solved on relatively inexpensive computational platforms. For the iterative solutions of the matrix equations, we are able to obtain accelerated convergence even for ill-conditioned matrix equations using advanced preconditioning schemes, such as nested preconditioned based on an approximate MLFMA. By orchestrating these diverse activities, we have been able to solve a closed geometry formulated with the CFIE containing 33 millions of unknowns and an open geometry formulated with the EFIE containing 12 millions of unknowns, which are the largest problems of their classes, to the best of our knowledge.
Open Access
Solutions of large-scale electromagnetics problems using an iterative inner-outer scheme with ordinary and approximate multilevel fast multipole algorithms
(2010) Ergül, A.; Malas, T.; Gürel, Levent
We present an iterative inner-outer scheme for the efficient solution of large-scale electromagnetics problems involving perfectly-conducting objects formulated with surface integral equations. Problems are solved by employing the multilevel fast multipole algorithm (MLFMA) on parallel computer systems. In order to construct a robust preconditioner, we develop an approximate MLFMA (AMLFMA) by systematically increasing the efficiency of the ordinary MLFMA. Using a flexible outer solver, iterative MLFMA solutions are accelerated via an inner iterative solver, employing AMLFMA and serving as a preconditioner to the outer solver. The resulting implementation is tested on various electromagnetics problems involving both open and closed conductors. We show that the processing time decreases significantly using the proposed method, compared to the solutions obtained with conventional preconditioners in the literature.
Open Access
Stabilization of integral-equation formulations for the accurate solution of scattering problems involving low-contrast dielectric objects
(Institute of Electrical and Electronics Engineers, 2008) Ergül, Özgür; Gürel, Levent
The solution of scattering problems involving low-contrast dielectric objects with three-dimensional arbitrary shapes is considered. Using the traditional forms of the surface integral equations, scattered fields cannot be calculated accurately if the contrast of the object is low. Therefore, we consider the stabilization of the formulations by extracting the nonradiating parts of the equivalent currents. We also investigate various types of stable formulations and show that accuracy can be improved systematically by eliminating the identity terms from the integral-equation kernels. Traditional and stable formulations are compared, not only for small scatterers but also for relatively large problems solved by employing the multilevel fast multipole algorithm. Stable and accurate solutions of dielectric contrasts as low as 104 are demonstrated on problems involving more than 250000 unknowns.