Abstract

We have developed an efficient iterative algorithm for electromagnetic scattering of arbitrary but relatively smooth dielectric objects. The algorithm iteratively adapts the equivalent surface currents until the electromagnetic fields inside and outside the dielectric objects match the boundary conditions. Theoretical convergence is analyzed for two examples that solve scattering of plane waves incident upon air/dielectric slabs of semi-infinite and finite thicknesses. We applied the iterative algorithm for simulation of sinusoidally-perturbed dielectric slab on one side and the method converged for such unsmooth surfaces. We next simulated the shift in radiation pattern of a 6-inch dielectric lens for different offsets of the feed antenna on the focal plane. The result is compared to that of the Geometrical Optics (GO).

© 2011 Optical Society of America

1. Introduction

Electromagnetic scattering of dielectric objects of arbitrary shape and size is difficult to simulate with reasonable speed and accuracy. There exist many computational methods such as Finite Difference Time Domain (FDTD) [1], Finite Element (FEM) [2], and Method of Moments (MoM) [3] which can provide accurate results, but the computation times become prohibitive for electrically large objects. For example, FDTD and FEM require generation of dense computational mesh of scattering object volume, with problem freedom or unknowns scaling as M3(we have assumed the number of mesh nodes Mx = My = Mz = M to be equal for all 3 dimensions). MOM is better in terms of the number of unknowns since it deals with boundary surface of the object instead of the volume of the object. However, MOM requires the inverse of the impedance matrix, which takes a computational complexity of N3 for N unknowns. The required memory scales as N2. Although Multi-level Fast Multipole Method (MLFMM) can reduce the computational complexity to N log N [4], but its computational procedure is too complicated and still not so efficient for some ultra-large-scale object simulation. Other boundary element methods can only simulate some particular geometries efficiently [58]. Finally, methods such as Physical Optics (PO) and Geometric Optics (GO) are generally faster but less accurate.

In this paper, we describe a new iterative algorithm based on finding the equivalent surface currents (both magnetic and electrical) on the dielectric surface. Fast and accurate convergence of the iterative algorithm is theoretically validated for plane wave incidence upon dielectric slabs of semi-infinite and finite thicknesses. We then applied the iterative algorithm on two examples: 1) sinusoidally-perturbed dielectric slab on one side and 2) a 6-inch dielectric lens with different offsets that we used for a passive millimeter wave radiometer [9], both at the frequency of 150 GHz. Also, the radiation patterns for different offsets of the lens with respect to a horn antenna were obtained and compared with the GO results.

2. Equivalent Surface Current Model for Scattering Problem

Fig. 1 shows the simulation based on equivalent surface currents on an arbitrary dielctric surface. The equivalent surface currents are given by,

Js=2n^×H;Ms=2n^×E
where (E, H) are the total fields expressed by the sum of the incident (Ei, Hi) and scattered fields (Esca, Hsca),
H=Hi+Hsca;E=Ei+Esca

 

Fig. 1 The equivalent current approach to the scattering problem.

Download Full Size | PPT Slide | PDF

Now the potentials A and F are given by,

A=SdsJsexpjk|rrs|4π|rrs|
and
F=ɛSdsMsexpjk|rrs|4π|rrs|
from which we can obtain the scattered electric and magnetic fields [10],
Esca=EA+EFHsca=HA+HF
with the following definitions,
EF(r)=×F(r)ɛHF(r)=1jωμ×EF(r)=××F(r)jωμɛ
and
HA(r)=×A(r)μEA(r)=1jωɛ×HA(r)=××A(r)jωμɛ

The calculation of the electromagnetic fields in Eq. (6) and Eq. (7) can be done efficiently through the Taylor-FFT algorithm developed by the author in [11]. The problem now reduces to finding a way to obtain the equivalent surface currents (Ms, Js), which is the main subject of this paper.

3. The Iterative Algorithm

The iterative algorithm aims at obtaining the equivalent surface currents Js and Ms for the scattering problem. Fig. 1 shows that the total field of the scattering problem can be viewed as the sum of the incident and the scattered fields. One way to obtain the equivalent surface currents Js and Ms is [12]: the scattered fields Esca and Hsca outside the object are obtained when these surface currents radiate in a homogeneous medium ɛ0; while the negative scattered fields inside the object Eo and Ho are obtained when the these surface currents −Js and −Ms radiate in a homogeneous medium ɛr. In this way, Js and Ms can be obtained as [12]

Js=n^×H=n^×[Hi+Hsca]=n^×HoMs=n^×E=n^×[Ei+Esca]=n^×Eo

The algorithm begins with initial guess of (Js, Ms) through approximate incident fields Ei and Hi; then scattered fields (Esca,k, Hsca,k) of the kth iteration are calculated as follows,

Esca,k=Esca,k1+n^×1Y[Hi,k+Hsca,k1Ho,k1]Hsca,k=Hsca,k1n^×1Z[Ei,k+Esca,k1Eo,k1]
where Y¯=12[1η0+1ηr], Z¯=12[η0+ηr], with η0=μ0ɛ0 and ηr=μrɛr being the wave impedances for vacuum and dielectric materials respectively; so we have the equivalent surface currents (Js,k, Ms,k) of the kth iteration according to Eq. (8),
Js,k=Js,k1+1Z[Ei+Esca,k1Eo,k1]Ms,k=Ms,k1+1Y[Hi+Hsca,k1Ho,k1]

The iteration continues until it converges for a given criterion. Fig. 2 shows the detailed procedures for the iterative algorithm, which are summarized below,

  1. Make the initial guess of the total field (E0, H0): although there is no specific requirement of the initial guess, it is preferable to use reasonable value so that the iteration converges fast. The simplest guess is to begin with the incident field (Ei, Hi);
  2. Update the equivalent surface currents (Ms,k, Js,k) of the kth iteration according to Eq. (10).
  3. Update the scattered fields on both sides (inside and outside) of the dielectric object.
  4. Correct the total field at kth iteration (Ek, Hk) according to Eq. (2).
  5. If the correction is small compared to some criterion, the iterative algorithm converges and go to step 6) below; otherwise, repeat step 2) to step 4) until the algorithm converges.
  6. Calculate the far field pattern.

 

Fig. 2 The iterative algorithm implementation procedures (also see text for details).

Download Full Size | PPT Slide | PDF

4. Convergence Analysis of the Iterative Algorithm

Before we show the numerical results for the lens simulation, let us first give a general analysis of the convergence of the iterative algorithm and apply it to two basic examples: i.e., air/dielectric interface and dielectric slab with finite thickness. Both examples assume plane wave incidence.

4.1. General Analysis

Following Eq. (10), for plane wave incidence of arbitrary number of dielectric slabs, the iterative algorithm follows the correction procedure shown in Fig. 3, i.e., the kth iterative equivalent surface currents Ms,k, Js,k can be expressed using the previous ones,

[Ms,kJs,k]=t[Miη0ηrJi]+C=[Ms,k1Js,k1]
where
C==[α=β=θ=γ=];[MiJi]=[Ei×n^n^×Hi]

 

Fig. 3 The equivalent currents approach to the scattering problem.

Download Full Size | PPT Slide | PDF

From Eq. (11) we can obtain the convergence of the iterative algorithm by letting k → ∞,

[MsJs]=tk=0[α=β=θ=γ=]k[Miη0ηrJi]+[α=β=θ=γ=][Ms,1Js,1]

Eq. (13) can be solved by expressing the matrix in its eigenvalue-eigenvector form,

[α=β=θ=γ=]=Q=[Λ=100Λ=2]Q=1[α=β=θ=γ=]k=Q=[Λ=1k00Λ=2k]Q=1[α=β=θ=γ=]=0=

Eq. (13) converges if all the diagonal elements in matrices Λ̿1,2 are smaller than unit 1, and the result is,

[MsJs]=tQ=[11Λ=10011Λ=2]Q=1[Miη0ηrJi]

4.2. Semi-infinite dielectric slab

For a plane wave (assuming Ei = x̂Ei,x; Hi = ŷHi,x) incident upon an air-dielectric (ɛr) interface, let us denote the total field of the (k − 1)th iteration as (Ek−1, Hk−1) = (x̂Ek−1, ŷHk−1). The equivalent surface currents are (we choose the surface normal to be = −),

Ms,k1=(z^)×Ek1=y^Ek1;Js,k1=(z^)×Hk1=x^Hk1
Then the scattered field is given by,
Esca,k=x^[12Ek1η02Hk1];Hsca,k=y^[12η0Ek1+12Hk1]
Eo,k=x^[12Ek1+ηr2Hk1];Ho,k=y^[12ηrEk1+12Hk1]

The updated equivalent surface currents are given by Eq. (11) and we find that α̿ = β̿ = θ̿ = γ̿ = 0, so

[MsJs]=t[Miη0ηrJi]
or
[EH]=t[Eiη0ηrHi]
as expected. Interestingly, for infinite dielectric slab, only one iteration is required to obtain the result.

4.3. Finite dielectric slab

Now let us extend the problem to the finite dielectric slab. In this case, we have surface currents on both left and right planes, with surface normals of = −+ = −(− and + denote the left and right sides respectively). Similarly, let us denote the total field of the (k − 1)th iteration as follows (Ek1±,Hk1±)=(x^Ek1±,y^Hk1±). The equivalent surface currents are,

Ms,k1±=y^Ek1;Js,k1±=x^Hk1
Then the scattered field is given by,
Esca,k±=x^[Ek1±2±η02Hk1±(Ek12±η02Hk1)expjkd]Hsca,k±=y^[±Ek1±2η0+Hk1±2+(Ek12η0Hk12)expjkd]
Eo,k±=x^[Ek1±2ηr2Hk1±+(Ek12±ηr2Hk1)expjkrd]Ho,k±=y^[Ek1±2ηr+Hk1±2+(±Ek12ηr+Hk12)expjkrd]
where k and kr are the wave vectors of the air and the dielectric.

From Eq. (10), we have,

[Ms,k+Js,k+Ms,kJs,k]=t[Mi+η0ηrJi+Miη0ηrJi]+C=[Ms,k1+Js,k1+Ms,k1Js,k1]
where,
C==[00𝒫η0+𝒫rηr2Y=𝒫+𝒫r2Y=00𝒫+𝒫r2Z=η0𝒫+ηr𝒫r2Z=𝒫η0+𝒫rηr2Y=𝒫+𝒫r2Y=00𝒫+𝒫r2Z=η0𝒫+ηr𝒫r2Z=00]

The eigenvalues of matrix C̿ are,

[Λ=100Λ=2]=[r𝒫r0000r𝒫r0000r𝒫0000r𝒫]
where r=ηrηηr+η is the reflection coefficient; 𝒫 = expjkd and 𝒫r = expjkrd. Since |r𝒫r| = |r𝒫| ≤ 1, the iterative algorithm converges according to Eq. (15).

5. Numerical Results

We applied the iterative algorithm on two examples: 1) Gaussian beam propagation through a dielectric slab with sinusoidal shape on one side and 2) dielectric lens simulation with different offsets of feed horn on the focal plane.

5.1. Dielectric slab with sinusoidal shape on one side

Fig. 4 shows the schematic of the simulation: a Gaussian beam source is 1” away from the dielectric slab,

Ex(x,y)=expx2+y2w2
with w = 5λ = 10 mm in the simulation. The sinusoidal shape is given by
z(x,y)=10λ+λcos(k40x)cos(k40y)

 

Fig. 4 The application of iterative algorithm on dielectric slab with sinusoidal shape on one side.

Download Full Size | PPT Slide | PDF

Fig. 5 shows the convergence of tangential fields along x-direction on the sinusoidal surface of the dielectric slab after 7 iterations. As can be seen from the Fig. 5, the accuracy is up to −120 dB, good enough for most applications.

 

Fig. 5 E// and H// of outside (air) and inside (dielectric) on the sinusoidal surface of the dielectric slab along x-direction after 7 iterations. Both amplitude and real part are shown.

Download Full Size | PPT Slide | PDF

5.2. Dielectric lens simulation

Lens simulation is important to obtain accurate radiation pattern [1315]. We apply the iterative algorithm on the simulation of a 6-inch dielectric lens at the frequency of 150 GHz for feed antenna at different offsets on the focal plane (see Fig. 6). The antenna pattern is approximated by a Gaussian beam as,

Ex(x,y)=exp(xxoff)2+y2w2
where xoff is the antenna offset.

 

Fig. 6 The application of iterative algorithm on horn-lens radiation pattern simulation.

Download Full Size | PPT Slide | PDF

The simulated lens is a thin spherical lens given by,

z(x,y)=R2x2y2(Rh)
where R = 3.4″ is the radius of the lens and h = 1.8″ is the thickness of the lens. The focal length F is given by the lens maker’s formula,
FRn1=2R=6.8
with n = 1.5 for our simulation.

Fig. 7 shows the tangential component E// of the electric field and the tangential component H// of the magnetic field inside and outside the dielectric surface after the algorithm runs for 7 iterations. The corresponding radiation patterns are shown in Fig. 8, together with the GO results for comparison. The GO result is obtained by using Snell’s law and Fresnel’s law on the plane surface of the lens to obtain the field inside the lens; the obtained field then is propagated to the convex surface of the lens, where again the Snell’s law and Fresnel’s law are applied to obtain the transmitted field outside the lens.

 

Fig. 7 E// and H// of outside (air) and inside (dielectric) on the convex surface of the dielectric lens along x-direction after 7 iterations. Only amplitude is shown.

Download Full Size | PPT Slide | PDF

 

Fig. 8 Radiation patterns for different antenna offsets along x-direction: from 0λ to 4λ, 1λ increment: lines are results from iterative algorithm and dots are GO results.

Download Full Size | PPT Slide | PDF

6. Discussion

We compared the efficiency of our algorithm in terms of the memory costs and computational time to different computational electromagnetics methods, with result shown in Table 1. Without loss of generality, we chose methods of FDTD, FEM and MOM using MLFMM (MOM-MLFMM) for comparison. In Table 1, we have assumed the dimension sizes are equal in all 3 dimension, Lx = Ly = Lz = L, which also means the numbers of mesh cells are the same Mx = My = Mz = M. We also assumed that the boundary surface of the object is proportional to the square of dimension, i.e., SL2 (e.g., a simple cube). It is straight forward that the memory for FDTD and FEM is at the order of 𝒪(M3) [1, 2]. The computation time for FDTD scales as 𝒪(M3 × M) = 𝒪(M4) since it takes ∝ M time step to propagate a pulse through the scattering object of length L [1]. The computation time for FEM using an iterative solver can be as low as 𝒪(M3) [2]. It is noted that the FDTD computation time is M times as much as that of the FEM. This is because FDTD can give 𝒪(M) frequency points after the time domain has been transformed to the frequency domain using FFT, while FEM only gives one frequency point. The memory for MOM-MLFMM can achieve the order of 𝒪(M2 log M2) and the computation time using MLFMM scales as 𝒪(M2 log M2) [4]. At last, the memory for our algorithm is 𝒪(M2). The computation time is 𝒪(M2 log M2) when fast field propagation method like FFT-based method [11] is used. Now let’s apply the above analysis to our two examples. In example 1, Lx = Ly = 127λ, Lz = 11λ. In FDTD and FEM, assuming the discretization size ΔL = 1/10λ, we have, Mx = My = 1270, Mz = 110. The memory for both FDTD and FEM is 𝒪(MxMyMz) = 𝒪(1.8 × 108). The computation time for FDTD and FEM are 𝒪(MxMyMz × max (Mx, My, Mz)) = 𝒪(2.3 × 1011) and 𝒪(MxMyMz) = 𝒪(1.8 × 108) respectively. The memory and computation time for MOM-MLFMM are the same, i.e., 𝒪(MxMy log (MxMy)) = 𝒪(2.3 × 107). The memory and computation time for our algorithm are 𝒪(MxMy) = 𝒪(1.6 × 106) and 𝒪(MxMy log(MxMy)) = 𝒪(2.3 × 107) respectively. In example 2, we used the following discretization: Mx = My = 760, Mz = 230 for ΔL = 1/10λ. The result is shown in Table 1. Obviously, our algorithm is more efficient in memory requirement than or at least as efficient as other methods used here for comparison. In terms of computation time, our algorithm is at least as efficient as MOM-MLFMM (only a few iterations are needed: 7 in both of our examples) and is much simpler to be implemented.

Tables Icon

Table 1. Memory and Computation Time (C. T.) for Different Methods

In updating the surface currents (Ms, Js) in Eq. (10), we have implicitly used the geometric optics approximation, which means our algorithm is more efficient for relatively smooth surface, i.e., the smoother the surface, the fewer the iterations it takes to converge. The rule of thumb of this is that good convergence could be obtained if each of the principal radius of curvature of the object surface is greater than a few wavelengths [16].

7. Conclusion

We have developed an efficient iterative algorithm to simulate the electromagnetic scattering by dielectric objects of different shapes, including a dielectric lens. The new algorithm iteratively adapts the equivalent surface currents to match the boundary conditions inside and outside the dielectric objects. The iterative algorithm converges after just a few iterations. Theoretical convergence analysis has been done for plane wave incidence upon dielectric slab. The algorithm has been applied to two examples, i.e., dielectric slab with perturbation on one side and a 6-inch dielectric lens with different offsets of a feed antenna on the focal plane. A comparison has been made among our algorithm and other commonly used computational electromagnetics methods, which shows that our algorithm is more efficient or at least as efficient as the state-of-the-art method like MOM-MLFMM and our algorithm is straightforward and simpler to be implemented.

References and links

1. A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd. ed. (Artech House, 2005).

2. J. L. Volakis, L. C. Kempel, and A. Chatterjee, Finite Element Method Electromagnetics (IEEE Computer Society Press, 1998). [CrossRef]  

3. R. F. Harrington, Field Computation by Moment Methods (Wiley IEEE Press, 1993). [CrossRef]  

4. W. C. Chew, J. M. Jin, E. Michielssen, and J. Song, Fast and Efficient Algorithm in Computational Electromagnetics (Artech House Publisher, 2001).

5. C. M. Kelso, P. D. Flammer, J. A. DeSanto, and R. T. Collins, “Integral equations applied to wave propagation in two dimensions: modeling the tip of a near-field scanning optical microscope,” J. Opt. Soc. Am. A 18(8), 1993–2001 (2001). [CrossRef]  

6. Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009). [CrossRef]  

7. M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antennas Propag. 51(6), 1142–1149 (2003). [CrossRef]  

8. X. An and Z. Q. Lu, “An efficient finite element-boundary integral method solving electromagnetic scattering problems,” Microwave Opt. Technol. Lett. 51(9), 2065–2071 (2009). [CrossRef]  

9. N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010). [CrossRef]  

10. C. A. Balanis, Advanced Engineering Electromagnetics, (John Wiley & Sons, 1989).

11. S. Liao and R. J. Vernon, “A fast algorithm for computation of electromagnetic wave propagation in half-space,” IEEE Trans. Antennas Propag. 57(7), 2068–2075 (2009). [CrossRef]  

12. S. B. Sorensen and K. Pontoppidan, Lens analysis methods for quasioptical systems, in The 2nd European Conference on Antennas and Propagation (EuCAP 2007), Edinburgh, UK, 11–16 Nov. 2007.

13. J. P. Thakur, W.-G. Kim, and Y.-H. Kim, “Large aperture low aberration aspheric dielectric lens antenna for W-band quasi-optics,” PIER 103, 57–65 (2010). [CrossRef]  

14. Z. X. Wang and W. B. Dou, “Full-wave analysis of monopulse dielectric lens antennas at W-band,” Int. J. Infrared Millim. Waves 31, 151–161 (2010).

15. A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antennas Propag.56, 485–492 (2008). [CrossRef]  

16. S. Liao and R. J. Vernon, “On the image approximation for electromagnetic wave propagation and PEC scattering in cylindrical harmonics,” Prog. Electromagn. Res. 66, 65–88 (2006). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. A. Taflove, and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd. ed. (Artech House, 2005).
  2. J. L. Volakis, L. C. Kempel, and A. Chatterjee, Finite ElementMethod Electromagnetics (IEEE Computer Society Press, 1998).
    [CrossRef]
  3. R. F. Harrington, Field Computation by Moment Methods (Wiley IEEE Press, 1993).
    [CrossRef]
  4. W. C. Chew, J. M. Jin, E. Michielssen, and J. Song, Fast and Efficient Algorithm in Computational Electromagnetics (Artech House Publisher, 2001).
  5. C. M. Kelso, P. D. Flammer, J. A. DeSanto, and R. T. Collins, “Integral equations applied to wave propagation in two dimensions: modeling the tip of a near-field scanning optical microscope,” J. Opt. Soc. Am. A 18(8), 1993–2001 (2001).
    [CrossRef]
  6. Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
    [CrossRef]
  7. M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
    [CrossRef]
  8. X. An, and Z. Q. Lu, “An efficient finite element-boundary integral method solving electromagnetic scattering problems,” Microw. Opt. Technol. Lett. 51(9), 2065–2071 (2009).
    [CrossRef]
  9. N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
    [CrossRef]
  10. C. A. Balanis, Advanced Engineering Electromagnetics, (John Wiley & Sons, 1989).
  11. S. Liao, and R. J. Vernon, “A fast algorithm for computation of electromagnetic wave propagation in half-space,” IEEE Trans. Antenn. Propag. 57(7), 2068–2075 (2009).
    [CrossRef]
  12. S. B. Sorensen, and K. Pontoppidan, Lens analysis methods for quasioptical systems, in The 2nd European Conference on Antennas and Propagation (EuCAP 2007), Edinburgh, UK, 11–16 Nov. 2007.
  13. J. P. Thakur, W.-G. Kim, and Y.-H. Kim, “Large aperture low aberration aspheric dielectric lens antenna for W-band quasi-optics,” PIER 103, 57–65 (2010).
    [CrossRef]
  14. Z. X. Wang, and W. B. Dou, “Full-wave analysis of monopulse dielectric lens antennas atW-band,” Int. J. Infrared Millim. Waves 31, 151–161 (2010).
  15. A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antenn. Propag. 56, 485–492 (2008).
    [CrossRef]
  16. S. Liao, and R. J. Vernon, “On the image approximation for electromagnetic wave propagation and PEC scattering in cylindrical harmonics,” Prog. Electromagn. Res. 66, 65–88 (2006).
    [CrossRef]

2010

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

J. P. Thakur, W.-G. Kim, and Y.-H. Kim, “Large aperture low aberration aspheric dielectric lens antenna for W-band quasi-optics,” PIER 103, 57–65 (2010).
[CrossRef]

Z. X. Wang, and W. B. Dou, “Full-wave analysis of monopulse dielectric lens antennas atW-band,” Int. J. Infrared Millim. Waves 31, 151–161 (2010).

2009

S. Liao, and R. J. Vernon, “A fast algorithm for computation of electromagnetic wave propagation in half-space,” IEEE Trans. Antenn. Propag. 57(7), 2068–2075 (2009).
[CrossRef]

Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
[CrossRef]

X. An, and Z. Q. Lu, “An efficient finite element-boundary integral method solving electromagnetic scattering problems,” Microw. Opt. Technol. Lett. 51(9), 2065–2071 (2009).
[CrossRef]

2008

A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antenn. Propag. 56, 485–492 (2008).
[CrossRef]

2006

S. Liao, and R. J. Vernon, “On the image approximation for electromagnetic wave propagation and PEC scattering in cylindrical harmonics,” Prog. Electromagn. Res. 66, 65–88 (2006).
[CrossRef]

2003

M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
[CrossRef]

2001

An, X.

X. An, and Z. Q. Lu, “An efficient finite element-boundary integral method solving electromagnetic scattering problems,” Microw. Opt. Technol. Lett. 51(9), 2065–2071 (2009).
[CrossRef]

Boriskin, A. V.

A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antenn. Propag. 56, 485–492 (2008).
[CrossRef]

Chan, C. H.

M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
[CrossRef]

Chien, H.-T.

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

Collins, R. T.

DeSanto, J. A.

Dou, W. B.

Z. X. Wang, and W. B. Dou, “Full-wave analysis of monopulse dielectric lens antennas atW-band,” Int. J. Infrared Millim. Waves 31, 151–161 (2010).

Elmer, T. W.

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

Flammer, P. D.

Godi, G.

A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antenn. Propag. 56, 485–492 (2008).
[CrossRef]

Gopalsami, N.

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

Heifetz, A.

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

Kelso, C. M.

Kim, W.-G.

J. P. Thakur, W.-G. Kim, and Y.-H. Kim, “Large aperture low aberration aspheric dielectric lens antenna for W-band quasi-optics,” PIER 103, 57–65 (2010).
[CrossRef]

Kim, Y.-H.

J. P. Thakur, W.-G. Kim, and Y.-H. Kim, “Large aperture low aberration aspheric dielectric lens antenna for W-band quasi-optics,” PIER 103, 57–65 (2010).
[CrossRef]

Koehl, E. R.

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

Lee, J. H.

Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
[CrossRef]

Li, S. Q.

M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
[CrossRef]

Liao, S.

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

S. Liao, and R. J. Vernon, “A fast algorithm for computation of electromagnetic wave propagation in half-space,” IEEE Trans. Antenn. Propag. 57(7), 2068–2075 (2009).
[CrossRef]

S. Liao, and R. J. Vernon, “On the image approximation for electromagnetic wave propagation and PEC scattering in cylindrical harmonics,” Prog. Electromagn. Res. 66, 65–88 (2006).
[CrossRef]

Lin, Y.

Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
[CrossRef]

Liu, J.

Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
[CrossRef]

Liu, Q. H.

Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
[CrossRef]

Lu, Z. Q.

X. An, and Z. Q. Lu, “An efficient finite element-boundary integral method solving electromagnetic scattering problems,” Microw. Opt. Technol. Lett. 51(9), 2065–2071 (2009).
[CrossRef]

Nosich, A. I.

A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antenn. Propag. 56, 485–492 (2008).
[CrossRef]

Raptis, A. C.

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

Sauleau, R.

A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antenn. Propag. 56, 485–492 (2008).
[CrossRef]

Simsek, E.

Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
[CrossRef]

Thakur, J. P.

J. P. Thakur, W.-G. Kim, and Y.-H. Kim, “Large aperture low aberration aspheric dielectric lens antenna for W-band quasi-optics,” PIER 103, 57–65 (2010).
[CrossRef]

Tsang, L.

M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
[CrossRef]

Vernon, R. J.

S. Liao, and R. J. Vernon, “A fast algorithm for computation of electromagnetic wave propagation in half-space,” IEEE Trans. Antenn. Propag. 57(7), 2068–2075 (2009).
[CrossRef]

S. Liao, and R. J. Vernon, “On the image approximation for electromagnetic wave propagation and PEC scattering in cylindrical harmonics,” Prog. Electromagn. Res. 66, 65–88 (2006).
[CrossRef]

Wang, Z. X.

Z. X. Wang, and W. B. Dou, “Full-wave analysis of monopulse dielectric lens antennas atW-band,” Int. J. Infrared Millim. Waves 31, 151–161 (2010).

Xia, M. Y.

M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
[CrossRef]

Zhang, B.

M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
[CrossRef]

IEEE Microw. Wirel. Compon. Lett.

Q. H. Liu, Y. Lin, J. Liu, J. H. Lee, and E. Simsek, “A 3-D spectral integral method (SIM) for surface integral equations,” IEEE Microw. Wirel. Compon. Lett. 19(2), 62–64 (2009).
[CrossRef]

IEEE Trans. Antenn. Propag.

M. Y. Xia, C. H. Chan, S. Q. Li, B. Zhang, and L. Tsang, “An efficient algorithm for electromagnetic scattering from rough surfaces using a single integral equation and multilevel sparse-matrix canonical-grid method,” IEEE Trans. Antenn. Propag. 51(6), 1142–1149 (2003).
[CrossRef]

S. Liao, and R. J. Vernon, “A fast algorithm for computation of electromagnetic wave propagation in half-space,” IEEE Trans. Antenn. Propag. 57(7), 2068–2075 (2009).
[CrossRef]

A. V. Boriskin, G. Godi, R. Sauleau, and A. I. Nosich, “Small hemielliptic dielectric lens antenna analysis in 2-D: boundary integral equations versus geometrical and physical optics,” IEEE Trans. Antenn. Propag. 56, 485–492 (2008).
[CrossRef]

Int. J. Infrared Millim. Waves

Z. X. Wang, and W. B. Dou, “Full-wave analysis of monopulse dielectric lens antennas atW-band,” Int. J. Infrared Millim. Waves 31, 151–161 (2010).

J. Opt. Soc. Am. A

Microw. Opt. Technol. Lett.

X. An, and Z. Q. Lu, “An efficient finite element-boundary integral method solving electromagnetic scattering problems,” Microw. Opt. Technol. Lett. 51(9), 2065–2071 (2009).
[CrossRef]

PIER

J. P. Thakur, W.-G. Kim, and Y.-H. Kim, “Large aperture low aberration aspheric dielectric lens antenna for W-band quasi-optics,” PIER 103, 57–65 (2010).
[CrossRef]

Proc. SPIE

N. Gopalsami, S. Liao, E. R. Koehl, T. W. Elmer, A. Heifetz, H.-T. Chien, and A. C. Raptis, “Passive millimeter wave imaging and spectroscopy system for terrestrial remote sensing,” Proc. SPIE 7670, 767003 (2010).
[CrossRef]

Prog. Electromagn. Res.

S. Liao, and R. J. Vernon, “On the image approximation for electromagnetic wave propagation and PEC scattering in cylindrical harmonics,” Prog. Electromagn. Res. 66, 65–88 (2006).
[CrossRef]

Other

A. Taflove, and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd. ed. (Artech House, 2005).

J. L. Volakis, L. C. Kempel, and A. Chatterjee, Finite ElementMethod Electromagnetics (IEEE Computer Society Press, 1998).
[CrossRef]

R. F. Harrington, Field Computation by Moment Methods (Wiley IEEE Press, 1993).
[CrossRef]

W. C. Chew, J. M. Jin, E. Michielssen, and J. Song, Fast and Efficient Algorithm in Computational Electromagnetics (Artech House Publisher, 2001).

C. A. Balanis, Advanced Engineering Electromagnetics, (John Wiley & Sons, 1989).

S. B. Sorensen, and K. Pontoppidan, Lens analysis methods for quasioptical systems, in The 2nd European Conference on Antennas and Propagation (EuCAP 2007), Edinburgh, UK, 11–16 Nov. 2007.

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1
Fig. 1

The equivalent current approach to the scattering problem.

Fig. 2
Fig. 2

The iterative algorithm implementation procedures (also see text for details).

Fig. 3
Fig. 3

The equivalent currents approach to the scattering problem.

Fig. 4
Fig. 4

The application of iterative algorithm on dielectric slab with sinusoidal shape on one side.

Fig. 5
Fig. 5

E// and H// of outside (air) and inside (dielectric) on the sinusoidal surface of the dielectric slab along x-direction after 7 iterations. Both amplitude and real part are shown.

Fig. 6
Fig. 6

The application of iterative algorithm on horn-lens radiation pattern simulation.

Fig. 7
Fig. 7

E// and H// of outside (air) and inside (dielectric) on the convex surface of the dielectric lens along x-direction after 7 iterations. Only amplitude is shown.

Fig. 8
Fig. 8

Radiation patterns for different antenna offsets along x-direction: from 0λ to 4λ, 1λ increment: lines are results from iterative algorithm and dots are GO results.

Tables (1)

Tables Icon

Table 1 Memory and Computation Time (C. T.) for Different Methods

Equations (31)

Equations on this page are rendered with MathJax. Learn more.

J s = 2 n ^ × H ; M s = 2 n ^ × E
H = H i + H sca ; E = E i + E sca
A = S d s J s exp j k | r r s | 4 π | r r s |
F = ɛ S d s M s exp j k | r r s | 4 π | r r s |
E sca = E A + E F H sca = H A + H F
E F ( r ) = × F ( r ) ɛ H F ( r ) = 1 j ω μ × E F ( r ) = × × F ( r ) j ω μ ɛ
H A ( r ) = × A ( r ) μ E A ( r ) = 1 j ω ɛ × H A ( r ) = × × A ( r ) j ω μ ɛ
J s = n ^ × H = n ^ × [ H i + H sca ] = n ^ × H o M s = n ^ × E = n ^ × [ E i + E sca ] = n ^ × E o
E sca , k = E sca , k 1 + n ^ × 1 Y [ H i , k + H sca , k 1 H o , k 1 ] H sca , k = H sca , k 1 n ^ × 1 Z [ E i , k + E sca , k 1 E o , k 1 ]
J s , k = J s , k 1 + 1 Z [ E i + E sca , k 1 E o , k 1 ] M s , k = M s , k 1 + 1 Y [ H i + H sca , k 1 H o , k 1 ]
[ M s , k J s , k ] = t [ M i η 0 η r J i ] + C = [ M s , k 1 J s , k 1 ]
C = = [ α = β = θ = γ = ] ; [ M i J i ] = [ E i × n ^ n ^ × H i ]
[ M s J s ] = t k = 0 [ α = β = θ = γ = ] k [ M i η 0 η r J i ] + [ α = β = θ = γ = ] [ M s , 1 J s , 1 ]
[ α = β = θ = γ = ] = Q = [ Λ = 1 0 0 Λ = 2 ] Q = 1 [ α = β = θ = γ = ] k = Q = [ Λ = 1 k 0 0 Λ = 2 k ] Q = 1 [ α = β = θ = γ = ] = 0 =
[ M s J s ] = t Q = [ 1 1 Λ = 1 0 0 1 1 Λ = 2 ] Q = 1 [ M i η 0 η r J i ]
M s , k 1 = ( z ^ ) × E k 1 = y ^ E k 1 ; J s , k 1 = ( z ^ ) × H k 1 = x ^ H k 1
E sca , k = x ^ [ 1 2 E k 1 η 0 2 H k 1 ] ; H sca , k = y ^ [ 1 2 η 0 E k 1 + 1 2 H k 1 ]
E o , k = x ^ [ 1 2 E k 1 + η r 2 H k 1 ] ; H o , k = y ^ [ 1 2 η r E k 1 + 1 2 H k 1 ]
[ M s J s ] = t [ M i η 0 η r J i ]
[ E H ] = t [ E i η 0 η r H i ]
M s , k 1 ± = y ^ E k 1 ; J s , k 1 ± = x ^ H k 1
E sca , k ± = x ^ [ E k 1 ± 2 ± η 0 2 H k 1 ± ( E k 1 2 ± η 0 2 H k 1 ) exp j k d ] H sca , k ± = y ^ [ ± E k 1 ± 2 η 0 + H k 1 ± 2 + ( E k 1 2 η 0 H k 1 2 ) exp j k d ]
E o , k ± = x ^ [ E k 1 ± 2 η r 2 H k 1 ± + ( E k 1 2 ± η r 2 H k 1 ) exp j k r d ] H o , k ± = y ^ [ E k 1 ± 2 η r + H k 1 ± 2 + ( ± E k 1 2 η r + H k 1 2 ) exp j k r d ]
[ M s , k + J s , k + M s , k J s , k ] = t [ M i + η 0 η r J i + M i η 0 η r J i ] + C = [ M s , k 1 + J s , k 1 + M s , k 1 J s , k 1 ]
C = = [ 0 0 𝒫 η 0 + 𝒫 r η r 2 Y = 𝒫 + 𝒫 r 2 Y = 0 0 𝒫 + 𝒫 r 2 Z = η 0 𝒫 + η r 𝒫 r 2 Z = 𝒫 η 0 + 𝒫 r η r 2 Y = 𝒫 + 𝒫 r 2 Y = 0 0 𝒫 + 𝒫 r 2 Z = η 0 𝒫 + η r 𝒫 r 2 Z = 0 0 ]
[ Λ = 1 0 0 Λ = 2 ] = [ r 𝒫 r 0 0 0 0 r 𝒫 r 0 0 0 0 r 𝒫 0 0 0 0 r 𝒫 ]
E x ( x , y ) = exp x 2 + y 2 w 2
z ( x , y ) = 10 λ + λ cos ( k 40 x ) cos ( k 40 y )
E x ( x , y ) = exp ( x x o f f ) 2 + y 2 w 2
z ( x , y ) = R 2 x 2 y 2 ( R h )
F R n 1 = 2 R = 6.8

Metrics