Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate
Electromagnetic Scattering 103 Table 1. Propagation of a plane harmonic TE wave in a waveguide. δ ≈ λ/15 Scheme spc time, s C A Yee 21 7 0.99613 1.00 Co-volume 46 29 0.99850 1.00 FETD 44 3151 1.0015 0.723 δ ≈ λ/30 Scheme spc time, s C A Yee 43 61 0.99896 1.00 Co-volume 106 263 0.99964 1.00 FETD 89 23040 1.0008 0.96 h , he Fig. 4. Propagation of a plane harmonic TE wave in a waveguide showing variation of the computed phase velocity with ∆t/〈h〉 (∆t/〈h e〉 for FETD). Solid symbols and solid line: δ ≈ λ/15; open symbols and dotted line: δ ≈ λ/30. Here 〈h〉 is the averaged Voronoï edge length, 〈h e〉 is the averaged minimal triangle height. 6.2 Scattering by a Circular PEC Cylinder The second example is the simulation of scattering of a plane single frequency TE wave by a perfectly conducting circular cylinder of diameter λ. The objective is to use this example to illustrate the order of accuracy that can be achieved with the co-volume solution technique and the FETD technique on unstructured meshes. The problem is solved on a series of unstructured meshes, with mesh spacings ranging from λ/8 toλ/128. The minimum distance from the rectangular PML to the cylinder is λ. When the spacing is
104 I. Sazonov et al. (a) (b) Fig. 5. Scattering of a plane TE wave by a circular PEC cylinder of diameter λ showing (a) an unstructured mesh with δ ≈ λ/16, (b) the corresponding computed total magnetic field. Scattering Width, dB Viewing Angle, degrees Fig. 6. Scattering of a plane TE wave by a circular PEC cylinder of diameter λ showing a comparison between the computed and analytical scattering width distributions. λ/16, the mesh employed, excluding the PML region, and the corresponding distribution of the computed total magnetic field is shown in Figure 5. The computed scattering width distributions are compared to the exact distribution in Figure 6. For each simulation undertaken, the error, E SW , in the solution is determined as the maximum difference, in absolute value, between the computed and analytical scattering width distributions. The variation of this computed error, with the number of elements per wavelength, λ/δ, forboth the FETD and co-volume schemes, is shown in Figure 7. It can be observed that a convergence rate of around O(δ 2 ) is obtained with both methods on these unstructured meshes, indicating that second-order accuracy is achieved. It is likely that the error in the FETD results is slightly less because the approach adopted for the evaluation of the scattering width integral requires an interpolation, in the co-volume scheme, to obtain all the field components at
- Page 59 and 60: 48 E.J. Dean and R. Glowinski S:T=
- Page 61 and 62: 50 E.J. Dean and R. Glowinski minim
- Page 63 and 64: 52 E.J. Dean and R. Glowinski 6 On
- Page 65 and 66: 54 E.J. Dean and R. Glowinski Fig.
- Page 67 and 68: 56 E.J. Dean and R. Glowinski and
- Page 69 and 70: 58 E.J. Dean and R. Glowinski 7 Num
- Page 71 and 72: 60 E.J. Dean and R. Glowinski Fig.
- Page 73 and 74: 62 E.J. Dean and R. Glowinski Assum
- Page 75 and 76: Higher Order Time Stepping for Seco
- Page 77 and 78: u n+1 h Optimal Higher Order Time D
- Page 79 and 80: Optimal Higher Order Time Discretiz
- Page 81 and 82: Optimal Higher Order Time Discretiz
- Page 83 and 84: Optimal Higher Order Time Discretiz
- Page 85 and 86: Optimal Higher Order Time Discretiz
- Page 87 and 88: Optimal Higher Order Time Discretiz
- Page 89 and 90: Optimal Higher Order Time Discretiz
- Page 91 and 92: Optimal Higher Order Time Discretiz
- Page 93 and 94: Optimal Higher Order Time Discretiz
- Page 95 and 96: Optimal Higher Order Time Discretiz
- Page 97 and 98: Optimal Higher Order Time Discretiz
- Page 99 and 100: Optimal Higher Order Time Discretiz
- Page 101 and 102: Optimal Higher Order Time Discretiz
- Page 103 and 104: 96 I. Sazonov et al. To provide a p
- Page 105 and 106: 98 I. Sazonov et al. In the first s
- Page 107 and 108: 100 I. Sazonov et al. Fig. 1. An ex
- Page 109: 102 I. Sazonov et al. H z 1 exact F
- Page 113 and 114: 106 I. Sazonov et al. Scattering Wi
- Page 115 and 116: 108 I. Sazonov et al. 6.4 Scatterin
- Page 117 and 118: 110 I. Sazonov et al. (a) (b) Fig.
- Page 119 and 120: 112 I. Sazonov et al. [MHP96] K. Mo
- Page 121 and 122: 114 R. Sanders and A.M. Tesdall I R
- Page 123 and 124: 116 R. Sanders and A.M. Tesdall imp
- Page 125 and 126: 118 R. Sanders and A.M. Tesdall alo
- Page 127 and 128: 120 R. Sanders and A.M. Tesdall (a)
- Page 129 and 130: 122 R. Sanders and A.M. Tesdall D C
- Page 131 and 132: 124 R. Sanders and A.M. Tesdall 8.6
- Page 133 and 134: 126 R. Sanders and A.M. Tesdall 0.3
- Page 135 and 136: 128 R. Sanders and A.M. Tesdall [TR
- Page 137 and 138: 132 S. Lapin et al. Ω R γ Ω 2 Γ
- Page 139 and 140: 134 S. Lapin et al. ∫ ∂ 2 ∫
- Page 141 and 142: 136 S. Lapin et al. Ω R γ Fig. 3.
- Page 143 and 144: 138 S. Lapin et al. 4 Energy Inequa
- Page 145 and 146: 140 S. Lapin et al. 5 Numerical Exp
- Page 147 and 148: 142 S. Lapin et al. Fig. 6. Contour
- Page 149 and 150: 144 S. Lapin et al. Fig. 9. Obstacl
- Page 151 and 152: Domain Decomposition and Electronic
- Page 153 and 154: Domain Decomposition Approach for C
- Page 155 and 156: Domain Decomposition Approach for C
- Page 157 and 158: Domain Decomposition Approach for C
- Page 159 and 160: Domain Decomposition Approach for C
104 I. Sazonov et al.<br />
(a)<br />
(b)<br />
Fig. 5. Scattering of a plane TE wave by a circular PEC cylinder of diameter λ<br />
showing (a) an unstructured mesh with δ ≈ λ/16, (b) the corresponding computed<br />
total magnetic field.<br />
Scattering Width,<br />
dB<br />
Viewing Angle, degrees<br />
Fig. 6. Scattering of a plane TE wave by a circular PEC cylinder of diameter<br />
λ showing a comparison between the computed <strong>and</strong> analytical scattering width<br />
distributions.<br />
λ/16, the mesh employed, excluding the PML region, <strong>and</strong> the corresponding<br />
distribution of the computed total magnetic field is shown in Figure 5. The<br />
computed scattering width distributions are compared to the exact distribution<br />
in Figure 6. For each simulation undertaken, the error, E SW , in the solution<br />
is determined as the maximum difference, in absolute value, between the<br />
computed <strong>and</strong> analytical scattering width distributions. The variation of this<br />
computed error, with the number of elements per wavelength, λ/δ, forboth<br />
the FETD <strong>and</strong> co-volume schemes, is shown in Figure 7. It can be observed<br />
that a convergence rate of around O(δ 2 ) is obtained with both methods on<br />
these unstructured meshes, indicating that second-order accuracy is achieved.<br />
It is likely that the error in the FETD results is slightly less because the approach<br />
adopted for the evaluation of the scattering width integral requires an<br />
interpolation, in the co-volume scheme, to obtain all the field components at