Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
106 I. Sazonov et al.<br />
Scattering Width, dB<br />
(a)<br />
Viewing Angle, degrees<br />
(b)<br />
Fig. 8. Scattering of a plane TE wave by a circular PEC cylinder of diameter<br />
15λ showing (a) a comparison between the exact <strong>and</strong> computed scattering width<br />
distributions, (b) computed contour distribution of the total magnetic field in the<br />
complete computational domain.<br />
6.3 Scattering by a Square PEC Cylinder<br />
The next example involves scattering of a plane single frequency electromagnetic<br />
wave by a perfectly conducting cylinder of square cross section. The<br />
sides of the square are of length λ. The objective is to use this example to<br />
illustrate the accuracy of the FETD <strong>and</strong> co-volume schemes in the presence<br />
of singularities. This simple geometry means that the computational domain<br />
may be discretised using a structured mesh of square elements <strong>and</strong>, in this<br />
case, the co-volume scheme of the equations (11) <strong>and</strong> (12) reduces to the classical<br />
Yee scheme. The distribution of the scattering width obtained using the<br />
Yee scheme on a fine Cartesian grid, with 512 elements per wavelength, is<br />
taken as the benchmark solution. An unstructured mesh, termed mesh a, is<br />
generated with mesh spacing λ/16. The solution is advanced on this mesh for 8<br />
cycles using both the co-volume <strong>and</strong> the FETD schemes. In this case, the error<br />
E SW is determined as the maximum difference, in absolute value, between the<br />
computed <strong>and</strong> the benchmark scattering width distributions. Table 3 shows<br />
the values of spc, time <strong>and</strong> E SW for this grid. It is apparent that the error<br />
in the FETD scheme is an order of magnitude greater than the error in the<br />
co-volume method. This is believed to be due to the singularity in the geometry,<br />
where higher mesh resolution will be required in a scheme such as FETD.<br />
Two further unstructured meshes are generated, by reducing the spacing by<br />
a factor of 2 (termed mesh b) <strong>and</strong> 4 (termed mesh c), in the vicinity of the<br />
corners. Details of the three meshes in the region of one of the corners are<br />
shown in Figure 9. Figure 10 shows the variation in the computed error with<br />
the near corner resolution that is employed. It can be seen that the error in<br />
the FETD results decreases as the mesh is refined. It is also clear that the