View source for EM.Ferma Tutorial Lesson 7: Analyzing A Microstrip Transmission Line

{{projectinfo|Tutorial| Analyzing A Micrsotrip Transmission Line|Ferma L7 Fig title.png|In this project, you will analyze a microstrip transmission line using EM.Ferma's quasi-static simulation engine.|

*Fixed-Potential PEC Object
*Dielectric Object
*2D Solution Plane
*Quasi-Static Analysis
*Characteristic Impedance
*Effective Permittivity
*Variables
*Parametric Sweep
|All versions|{{download|http://www.emagtech.com/downloads/ProjectRepo/EMFerma_Lesson7.zip EMFerma_Lesson7}} }}

== What You Will Learn ==

In this tutorial you will learn how to define and characterize a 2D transmission line structure like a microstrip line using quasi-static analysis. [[EM.Ferma]] provides three transmission line wizards:

#Microstrip Wizard
#CPW Wizard
#Coaxial Wizard

You will also learn how to set up and perform a parametric sweep of your design variables.

[[Image:Back_icon.png|30px]] '''[[EM.Ferma  | Back to EM.Ferma Manual]]'''

[[Image:Back_icon.png|30px]] '''[[EM.Cube#EM.Ferma_Documentation | Back to EM.Ferma Tutorial Gateway]]'''

[[Image:Download2x.png|30px]] '''[http://www.emagtech.com/downloads/ProjectRepo/EMFerma_Lesson7.zip Download projects related to this tutorial lesson]'''

== Getting Started ==

Open the [[EM.Cube]] application and switch to [[EM.Ferma]]. Start a new project with the following parameters: 

<div class="noprint" style="float:left;margin-right:10px">
{| border="1" class="wikitable"
|+ Starting Parameters
|-
! Name
| EMFerma_Lesson7
|-
! Length Units
| Millimeters
|-
! Frequency Units
| N/A
|-
! Center Frequency
| N/A
|-
! Bandwidth
| N/A
|}
</div>

<div class="noprint" style="clear:both"></div>

== A Note on Quasi-Static Solution of Transmission Line Problems ==

Most popular two-conductor transmission line types like microstrip, stripline, coplanar waveguide (CPW), <i>etc</i>., feature a single dominant TEM or TEM-like propagating mode. Therefore, they can be reasonably modeled using quasi-static methods. [[EM.Ferma]] takes the cross section of an infinitely long 2D transmission line structure and encloses it by a PEC domain box. Shielded transmission lines are naturally well-represented using this approach. To model open-boundary transmission lines like a microstrip line, you must place the domain boundary walls far enough from the signal line.  

[[EM.Ferma]] can perform 2D electrostatic and quasi-static simulations. For this purpose, you need to define a "2D Solution Plane". The cross section of your physical structure on this 2D solution plane is taken as the geometry to be solved. In other words, [[EM.Ferma]] assumes that this cross sectional geometry is invariably extended to the infinity from both sides of the solution plane.

To compute the characteristic impedance Z<sub>0</sub> and effective permittivity &epsilon;<sub>eff</sub> of the transmission line, [[EM.Ferma]] first analyzes the structure "As Is" and finds the electric field everywhere in the computational domain. Gauss' law is used to compute the total charge on the signal line, from which the capacitance C is calculated. Then, all the dielectric parts are replaced with air (&epsilon;<sub>r</sub> = 1). The electrostatic problem is solve one more time and the "air-filled" capacitance C<sub>a</sub> is calculated. The characteristic impedance and effective permittivity of the transmission line are then calculated from the following relations:

<math> Z_0 = \frac{1}{c_0 \sqrt{CC_a} } </math>  

<math> \epsilon_{eff} = \frac{C}{C_a} </math>  

where c<sub>0</sub> = 3&times;10<sup>8</sup> m/s is the speed of light in the free space.

==Creating the Microstrip Line Geometry ==

A microstrip line consists of a metallic strip printed on top of a dielectric substrate with a conductor backing on its other side. Click on the <b>Microstrip Wizard</b> [[Image:MicrostripWizardIconx.png]] button of the Wizard Toolbar or select the menu item '''Tools &rarr; Transmission Line Wizards &rarr; Microstrip Line'''.

<table>
<tr>
<td> [[Image:Ferma L7 Fig1.png|thumb|left|720px|EM.Ferma's wizard toolbar.]] </td>
</tr>
</table>

The geometry of a microstrip line appears at the center of the project workspace. 

<table>
<tr>
<td> [[Image:Ferma L7 Fig3.png|thumb|left|720px|The microstrip geometry created by the wizard in the project workspace.]] </td>
</tr>
</table>

The geometry created by the wizard consists of three objects belonging to three different material groups:

{| class="wikitable"
|-
! scope="col"| Object Name
! scope="col"| Object Type
! scope="col"| Material Group Name
|scope="col"| Material  Group Type
|-
| ANCHOR
| Rectangle Strip
| STRIP
| Fixed-potential PEC 
|-
| Ground
| Rectangle Strip
| GROUND
| Fixed-potential PEC 
|-
| Substrate
| Box
| SUBSTRATE 
| Dielectric material
|}

Open the property dialogs of the three material groups and examine their parameters. The fixed voltage of the "STRIP" PEC group is 1V, while the fixed voltage of the "GROUND" PEC group is 0V. The geometry created by the wizard is fully parameterized. Open the variables dialog and review the list of the variables used for the definition of the microstrip line. Change the value of "strip_wid" to 4.62893mm. You will notice that the "sub_size" will change to 46.2893mm.

<table>
<tr>
<td> [[Image:Ferma L7 Fig4.png|thumb|left|550px|The variables dialog showing the parameters of the microstrip line.]] </td>
</tr>
</table>

Open the property dialog of the rectangle strip object called "ANCHOR", which represents the metallic strip on the top side of the dielectric substrate. From the variables dialog shown above, note that the width of the transmission line has been defined as a function of another variable called "z0", which represents the characteristics impedance. The current value of z0 is 50&Omega;. In other words, the Y-dimension of the rectangle strip "ANCHOR" is controlled by the value you set for the variable "z0". Also note that since [[EM.Ferma]] uses the finite difference method for quasi-static analysis of 2D transmission lines, it requires a substrate of finite size. For this reason, the wizard has set the widths of the substrate and ground equal to 10 times the width of the microstrip.     

<table>
<tr>
<td> [[Image:Ferma L7 Fig4A.png|thumb|left|480px|The property dialog of the rectangle strip object called "ANCHOR".]] </td>
</tr>
</table>

==Examining the Computational Domain & the 2D Solution Plane ==

So far in all the previous tutorial lessons, you dealt with 3D domain boxes and performed 3D electrostatic or 3D magneostatic simulations of 3D physical structures. [[EM.Ferma]] also lets you perform 2D electrostatic analysis of 2D structures. By a 2D structure we mean a geometry that is invariant along one of the three principal axes and has an infinite extent along that direction. The cross section of a 2D structure is identical at any point along the so-called invariance axis. Therefore, you can build a 3D transmission line geometry just like other structures you built earlier and then specify a "2D Solution Plane" that would capture the cross section of your transmission line geometry and would pass its mesh to the simulation engine. In [[EM.Ferma]], 2D solution planes are defined based on existing field sensor observables. Since in [[EM.Ferma]], a field sensor plane extends across the entire computational domain, it can also serve to capture the cross section of your transmission line geometry.

Open the domain settings dialog. You will see that the wizard has set the domain offest values in the &plusmn;X and &plusmn;Y directions equal to zero. The -Z-offset is zero, too. 

<table>
<tr>
<td> [[Image:Ferma L7 Fig5.png|thumb|left|480px|EM.Ferma's domain settings dialog.]] </td>
</tr>
</table>

Remember that all the six walls of the domain box are held at zero potential by default (due to the Dirichlet boundary condition). This means that all six walls are assumed to be fix-potential PEC surface with a fixed zero voltage. You may pick a nonzero value for the -Z-offset parameter, but this will create a redundant air region between the object "Ground" and the bottom domain wall, both of which have zero voltages. The table below lists all the domain offset values.

{| class="wikitable" 
|-
! Offset !! Value
|-
| -X || 0
|-
| +X || 0
|-
| -Y || 0
|-
| +Y || 0
|-
| -Z || 0
|-
| +Z || 20mm
|-
|}

Also, note that the wizard defined a vertical X-directed field sensor observable called "FS1", which is centered at (0, 0, 0). 

== Running a Quasi-Static Simulation of the Microstrip Line ==

Open the mesh settings dialog. You will find that the wizard has set the mesh cell size to a very fine value of &Delta;x = &Delta;y = &Delta;z = 0.2mm. Since the width of your microstrip is about 4.6mm, this will place 23 cells across the width of the transmission line. To get accurate results, you typically need to set a very high resolution mesh for 2D transmission line structures. 

<table>
<tr>
<td> [[Image:Ferma L7 Fig7A.png|thumb|left|480px|High resolution fixed-cell mesh settings.]] </td>
</tr>
</table>

Run a quasi-static "Analysis" of your 2D transmission line structure. At the end of the simulation, the output message window reports the computed values of the characteristic impedance and effective permittivity of the transmission line:

Z0: 48.605429 Ohms 

Epsilon_Effective: 1.946599

<table>
<tr>
<td> [[Image:Ferma_L7_FigValues.png|thumb|left|600px|The computed values of the characteristic impedance and effective permittivity of the transmission line in the output message window.]] </td>
</tr>
</table>

Visualize the electric field and potential distribution on the Sensor_1 plane. From the vector plot, you can clearly see the field confinement inside the dielectric substrate underneath the micrsotrip as well as the fringing fields at the edges of the metal strip.   

<table>
<tr>
<td>
[[Image:Ferma L7 Fig8.png|thumb|left|640px|The intensity plot of the total electric field distribution in the YZ plane.]]
</td>
</tr>
<tr>
<td>
[[Image:Ferma L7 Fig8A.png|thumb|left|640px|The vector plot of the total electric field distribution in the YZ plane.]]
</td>
</tr>
<tr>
<td>
[[Image:Ferma L7 Fig9.png|thumb|left|640px|The intensity plot of the electric potential distribution in the YZ plane.]]
</td>
</tr>
</table>

== Verifying Your Simulation Results ==

To verify the computed results for the characteristic impedance Z<sub>0</sub> and &epsilon;<sub>eff</sub>, let's take a look at the analytical formulas available for the characteristic impedance and effective permittivity of a microstrip line:

<math> \epsilon_{eff} = \frac{\epsilon_r + 1}{2} +  \frac{\epsilon_r - 1}{2}  \left( 1+ 12\frac{h}{w} \right) ^{-\frac{1}{2}} </math>  

<math> Z_0 = \frac{120\pi}{\sqrt{\epsilon_{eff}}}  \left[ \frac{w}{h} + 1.393 + 0.677 ln \left(  \frac{w}{h} + 1.444 \right) \right] ^{-1} </math>  

where w and h are the microstrip width and substrate thickness, respectively, and &epsilon;<sub>r</sub> is the relative permittivity of the dielectric. The above formulas are valid for the case of w/h &ge; 1. For your microstrip line defined in this project, you have w = 4.7mm, h = 1.5mm and &epsilon;<sub>r</sub> = 2.2. For these values, the above formulas give: Z<sub>0</sub> = 49.54&Omega; and &epsilon;<sub>eff</sub> = 1.87, which are not far from the simulated results.

If you have installed [[RF.Spice A/D]] on your computer, you can also produce similar results in its Device Manager. In that case, you can access the [[RF.Spice A/D]] Device Manager directly from the '''Tool Menu''' of [[EM.Ferma]].  

<table>
<tr>
<td>
[[Image:Ferma L7 Fig10.png|thumb|left|500px|The microstrip line calculator of RF.Spice A/D.]]
</td>
</tr>
</table>

== Running a Parametric Sweep of the Microstrip Width ==

You saw earlier that the variable "sub_size" was defined as a function of the variable "strip_wid". This means that "sub_size" is a dependent variable, while "strip_wid" is an independent variable. You can easily change the definition of any variable and turn it into an independent or dependent variable. Open the variable dialog and change the values of two variables "h" and "strip_wid" according to the table below:

{| class="wikitable"
|-
! scope="col"| Variable Name
! scope="col"| Original Definition
! scope="col"| New Definition
|-
| h
| 0.0015*to_meters
| 1.5 
|-
| strip_wid
| 4.62893
| 4.6 
|}

<table>
<tr>
<td> [[Image:Ferma L7 Fig11.png|thumb|left|550px|The variables dialog showing the modified definition of some variables.]] </td>
</tr>
</table>

In this part of this tutorial lesson, you will vary the microstrip width and see how it affects the characteristic impedance and effective permitivity of your transmission line. In a sweep simulation, one or more parameters are varied, and the simulation engine is run for each parameter set. Open the Simulation Run dialog and choose the '''Parametric Sweep''' option from the '''Simulation Mode''' drop-down list. Click on the {{key|Settings}} button next to this drop-down list to open the Parametric Sweep Settings dialog. 

<table>
<tr>
<td> [[Image:Ferma L7 Fig12.png|thumb|left|480px|Selecting parametric sweep as the simulation mode in EM.Ferma's run dialog.]] </td>
</tr>
</table>

The sweep variables list is initially empty. 

<table>
<tr>
<td> [[Image:Ferma L7 Fig13.png|thumb|left|720px|The parametric sweep settings dialog.]] </td>
</tr>
</table>

On the left side of this dialog, you see a list of all the available independent variables. Select "strip_wid" from the left table and use the right arrow {{key|-->}} button to move it to the right table. Another dialog titled "Sweep Variable" opens up. You have to set the start, stop and step values of your sweep variable. By default, the sweep variable is of uniform type. Enter 3, 5, and 0.2 for the start, stop and step values, respectively. This will create a value list of {3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0}. Close the sweep variable dialog and then close the sweep settings dialog to return to the simulation run dialog.   

<table>
<tr>
<td> [[Image:Ferma L7 Fig14.png|thumb|left|560px|The sweep variable settings dialog.]] </td>
</tr>
</table>

<table>
<tr>
<td> [[Image:Ferma L7 Fig15.png|thumb|left|720px|The parametric sweep settings dialog showing "strip_wid" as the sweep variable.]] </td>
</tr>
</table>

Run the sweep simulation. It may take a while as a total of eleven individual 2D electrostatic simulations must be completed. At the end of the parametric sweep, open the data manager and plot the data files "SP1_Z0_Sweep.DAT" and "SP1_EpsEff_Sweep.DAT". You should see graphs like the figures below.   

<table>
<tr>
<td> [[Image:Ferma L7 Fig17.png|thumb|left|480px|The graph of variation of characteristic impedance as a function of microstrip width.]] </td>
</tr>
<tr>
<td> [[Image:Ferma L7 Fig18.png|thumb|left|480px|The graph of variation of effective permittivity as a function of microstrip width.]] </td>
</tr>
</table>

<p>&nbsp;</p>
[[Image:Top_icon.png|30px]] '''[[#What_You_Will_Learn | Back to the Top of the Page]]'''

[[Image:Back_icon.png|30px]] '''[[EM.Cube#EM.Ferma_Documentation | Back to EM.Ferma Tutorial Gateway]]'''