Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate
Electromagnetic Scattering 107 Table 3. Simulation of scattering of a plane TE wave by a square PEC cylinder of side length λ. Mesh FETD Co-volume Speed up resolution spc time, s E SW spc time, s E SW ratio a 61 18. 2.64 45 0.4 0.21 45 b 90 27. 1.66 88 0.8 0.25 34 c 182 58. 0.38 164 1.3 0.14 44 (a) (b) (c) Fig. 9. Details of the meshes employed for the simulation of scattering of a plane TE wave by a square PEC cylinder of side length λ showing (a) mesh a, (b) mesh b, (c) mesh c. E SW,dB Near-Corner Resolution Fig. 10. Simulation of scattering of a plane TE wave by a square PEC cylinder of side length λ showing the variation in the computed error with the near corner mesh resolution. error in the FETD results on mesh c is similar to the error in the co-volume results obtained on mesh a. The constant error in the co-volume results confirm the belief that no special modification of the scheme is required in the vicinity of geometrical singularities. Table 3 also displays information about the calculations performed on meshes b and c. For this example, the co-volume scheme is faster than FETD by a factor that ranges between 34 and 45. This level of variation in the speed-up factor is probably due to the difficulty in determining exactly the small times required for the co-volume solution.
108 I. Sazonov et al. 6.4 Scattering by a PEC NACA0012 Aerofoil The next example involves the simulation of scattering of a plane single frequency wave, directed along the x-axis, by a perfectly conducting NACA0012 aerofoil of length λ. The aim of this example is to analyse the performance of the numerical schemes when the geometry exhibits high curvature. A benchmark solution is computed using an unstructured mesh with spacing λ/120. The unstructured mesh is generated, outside the aerofoil, in the region −λ ≤ x, y ≤ λ. The scattering width distributions computed on this mesh with the co-volume scheme and the FETD scheme proved to be identical. An unstructured mesh was generated to meet the spacing requirement of λ/15. Another unstructured mesh, providing better representation of the leading edge curvature, is generated by locally reducing the mesh spacing in the vicinity of the leading edge of the airfoil by a factor of 2. A view of both these meshes is shown in Figure 11. The computed scattering width distributions are compared with the benchmark distribution in Figure 12. It can be observed that the co-volume results are better on the uniform mesh and that the accuracy of the FETD results improve with the local refinement in the leading edge region. For this example, Table 4 shows the values of spc, time and E SW . The co-volume method is approximately 30 times faster than FETD for this example. 6.5 Scattering by a PEC Cavity The final example considers the simulation of scattering of a plane single frequency wave by a U-shaped PEC cavity. The thickness of the cavity walls is equal to 0.4λ, the internal cavity width is 2λ and the internal cavity length is 8λ. In the simulation, the wave is incident upon the open end of the cavity and propagates in a direction which lies at an angle θ =30 ◦ to the main axis of the cavity. An unstructured mesh is employed, with typical spacing λ/15, in the region that lies within a distance of λ from the scatterer, as (a) (b) Fig. 11. Details of the unstructured meshes employed for the simulation of scattering of a plane TE wave by a PEC NACA0012 aerofoil of length λ showing (a) the uniform mesh, (b) the locally refined mesh.
- 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 and 110: 102 I. Sazonov et al. H z 1 exact F
- Page 111 and 112: 104 I. Sazonov et al. (a) (b) Fig.
- Page 113: 106 I. Sazonov et al. Scattering Wi
- 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
- Page 161 and 162: Domain Decomposition Approach for C
- Page 163 and 164: Domain Decomposition Approach for C
Electromagnetic Scattering 107<br />
Table 3. Simulation of scattering of a plane TE wave by a square PEC cylinder of<br />
side length λ.<br />
Mesh FETD Co-volume Speed up<br />
resolution spc time, s E SW spc time, s E SW ratio<br />
a 61 18. 2.64 45 0.4 0.21 45<br />
b 90 27. 1.66 88 0.8 0.25 34<br />
c 182 58. 0.38 164 1.3 0.14 44<br />
(a) (b) (c)<br />
Fig. 9. Details of the meshes employed for the simulation of scattering of a plane<br />
TE wave by a square PEC cylinder of side length λ showing (a) mesh a, (b) mesh<br />
b, (c) mesh c.<br />
E SW,dB<br />
Near-Corner Resolution<br />
Fig. 10. Simulation of scattering of a plane TE wave by a square PEC cylinder<br />
of side length λ showing the variation in the computed error with the near corner<br />
mesh resolution.<br />
error in the FETD results on mesh c is similar to the error in the co-volume<br />
results obtained on mesh a. The constant error in the co-volume results confirm<br />
the belief that no special modification of the scheme is required in the<br />
vicinity of geometrical singularities. Table 3 also displays information about<br />
the calculations performed on meshes b <strong>and</strong> c. For this example, the co-volume<br />
scheme is faster than FETD by a factor that ranges between 34 <strong>and</strong> 45. This<br />
level of variation in the speed-up factor is probably due to the difficulty in<br />
determining exactly the small times required for the co-volume solution.