24.02.2014 Views

Partial Differential Equations - Modelling and ... - ResearchGate

Partial Differential Equations - Modelling and ... - ResearchGate

Partial Differential Equations - Modelling and ... - ResearchGate

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!