EM.Tempo Tutorial Lesson 2: Analyzing Scattering From A Sphere
- 1 What You Will Learn
- 2 Getting Started
- 3 Drawing a Spherical Solid Object
- 4 Defining the Plane Wave Source
- 5 Defining the Project Observables
- 6 Running the FDTD Simulation
- 7 Visualizing the RCS of the PEC Sphere
- 8 Visualizing the Near Fields of the PEC Sphere
- 9 Adding a Dielectric Coating to Your PEC Sphere
- 10 Visualizing the RCS and Near Fields of the Coated Sphere
What You Will Learn
In this tutorial you will learn to draw spherical objects from the ground up. You will define material groups of different types, will set up a plane wave source to excite your physical structure and define field sensor and radar cross section (RCS) observables for your target.
For this tutorial lesson, set the following parameters in the new project dialog:
Drawing a Spherical Solid Object
|In EM.Tempo, all the new objects you draw are assumed to be perfect electric conductors (PEC) by default unless you define a material group of a different type.|
As soon as you draw the sphere, a default blue domain box is placed around it at an offset of 0.25λ0 away from the largest bounding box of your physical structure. Similar to the previous tutorial lesson, convolutional perfectly matched layer (PML) boundary conditions are assumed at all the six faces of the domain box. This is required for modeling an open-boundary problem.
Defining the Plane Wave Source
You will introduce a plane wave source for this project. Since FDTD is a finite-domain technique, it uses a finite-sized source box to excite the physical structure. To define a plane wave source, right-click on the Plane Waves item of the Sources section in the navigation tree and select Insert New Source… from the contextual menu. The Plane Wave Dialog opens up with a number of default settings.
The default plane wave box is places 0.2λ0 away from the largest bounding box of your entire physical structure. This is usually a good choice, and you will keep it for this project. Your plane wave source has an Amplitude of 1V/m and a zero Phase. You will keep the default TMz polarization. In the “Incident Angle” section of the dialog you need to enter the elevation θ and azimuth φ angles in the standard spherical coordinate system. Accept the default values: θ = 180o and φ = 0o, which represent an X-polarized normally incident plane wave. The above settings mean that the plane wave source will illuminate the spherical target with a modulated Gaussian pulse waveform centered at 10GHz with a bandwidth of 1GHz, where the envelope of the electric field reaches a maximum value of 1V/m. You will see a magenta wireframe box in the project workspace that completely encloses the sphere.
|In EM.Tempo, plane waves are characterized by their unit propagation vector. Therefore, for downward-looking and upward-looking normally incident plane waves along the ∓Z-axis, θ must be 180° and 0°, respectively.|
Defining the Project Observables
To calculate the radar cross section of your structure, you need to define an RCS observable. To do so, right click on the Radar Cross Sections item in the Observables section of the navigation tree, and select Insert New RCS… Accept most of the default values, but set both the Theta and Phi angle increments to 3°. The RCS box appears as a light blue box around your physical structure.
For this project, you will also define two field sensors at XY and YZ planes to observe the electric and magnetic field distributions. To define the field sensors, right-click on the Near Field Sensors item in the “Observables” section of the navigation tree and select Insert New Observable… In the Field Sensor Dialog, enter the point (0, 0, 0) for Coordinates and select Z from the Direction drop-down list for the XY plane. Then define Sensor_2 which is X-directed and centered at (0,0,0) for the YZ plane. Two orthogonal, purple boxes appear around your structure inside the computational domain.
Running the FDTD Simulation
Before running the FDTD simulation, check out the FDTD mesh of your physical structure. For this project, you will leave the mesh density at its default value of 20 cells/λeff. You will see that the curvature of the sphere is fairly well captured by EM.Tempo's adaptive mesh generator even at this relatively low mesh density. For this project, you will accept the default convergence criteria, that is a Power Threshold of -30dB and a Maximum Number of Time Steps equal to 10,000. Run the FDTD simulation and wait until it converges.
Visualizing the RCS of the PEC Sphere
At the end of the FDTD simulation, the Radar Cross Sections part of the navigation tree is populated with three 3D RCS plots: Theta RCS, Phi RCS and the total RCS. In order to better view the RCS plots, you have to either hide the target object or alternatively “freeze” it. Freezing an object replaces it with a wireframe outline of it, so you can still recognize its location and orientation in the project workspace. To freeze an object or a group of objects, right-click on their name and select Freeze from the contextual menu. To unfreeze a frozen object and get back to its normal display, repeat the same procedure and remove the check mark in front of the “Freeze” item in the contextual menu. Note that the maximum RCS is 0.009687m2 in the forward-scatter direction.
|Frozen objects cannot be selected using the mouse or cannot be highlighted using mouse-over.|
Next, open the data manager and plot the data files "RCS_1_RCS_Polar_YZ.ANG" and "RCS_1_RCS_Polar_ZX.ANG". Note that the polar plots of the RCS in the two principal planes are different.
Visualizing the Near Fields of the PEC Sphere
After the end of the FDTD simulation, 14 near-field plots are added under the nodes of each of your two field sensors: “Sensor_1” at XY plane and “Sensor_2” at YZ plane. The figures below show the total electric and magnetic fields at the XY and YZ planes. Note the null fields inside the PEC sphere.
Adding a Dielectric Coating to Your PEC Sphere
In this part of the tutorial lesson, you will analyze the scattering characteristics of a dielectric-coated PEC sphere. First, you need to define a new dielectric material group in the same project. Right-click on the Dielectric item under the “Physical Structure” section of the navigation tree and select Insert New Dielectric… The dielectric dialog opens up with the vacuum as the default material type. You can simply enter a name for your material and enter values for its constitutive parameters: relative permittivity εr, relative permeability μr, electric conductivity σ and magnetic conductivity σm. Alternatively, you can click the Material button of this dialog to access EM.Tempo’s Materials List.
Open the materials list and select "Polycarbonate" from the list. Click the OK button of the dialog to return to the dielectric dialog. You will see that the material properties are filled: εr = 3.2, μr = 1, σ = σm = 0. The last defined material group remains as the active material group. That means that you can now draw new objects under this material group. Draw a new sphere of radius R = 20mm center at (0, 0, 0). The new dielectric sphere in light green color is concentric with the old red PEC sphere, which is completely enclosed inside. If you hover your mouse over the larger sphere, its color becomes translucent and you will be able to see the enclosed smaller sphere.
EM.Tempo's FDTD solver is a volumetric solver. This means Maxwell's Equations are solved numerically at every single cell inside the computational domain. Although the two spheres are concentric, the cells inside the interior PEC sphere cannot have dielectric properties. EM.Tempo has a material hierarchy rule that takes care of object overlaps. This feature greatly facilitates the construction of physical structures with different material compositions.
|If two objects share overlap regions, the PEC material always supersedes the dielectric material.|
Examine the default FDTD mesh of the dielectric-coated PEC target and run an FDTD analysis using the default convergence criteria.
Visualizing the RCS and Near Fields of the Coated Sphere
Visualize the 3D RCS of your new target and plots its 2D polar graphs. Note how much the ZX-plane RCS has changed due to the dielectric coating. Also note that the maximum RCS has increase almost five-fold to 0.04803m2 in the forward-scatter direction.
Similarly, plot the electric and magnetic field distributions of the new target and compare them to the case of the PEC target without the dielectric coating.