Difference between revisions of "EM.Ferma Tutorial Lesson 2: Analyzing A Parallel Plate Capacitor"
(17 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
{{projectinfo|Tutorial| Analyzing A Parallel Plate Capacitor|Ferma L2 Fig title.png|In this project, you will analyze parallel plate capacitor structures made up of metal and dielectric parts.| | {{projectinfo|Tutorial| Analyzing A Parallel Plate Capacitor|Ferma L2 Fig title.png|In this project, you will analyze parallel plate capacitor structures made up of metal and dielectric parts.| | ||
− | *[[CubeCAD]] | + | *[[Building_Geometrical_Constructions_in_CubeCAD | CubeCAD]] |
*PEC Objects | *PEC Objects | ||
*Dielectric Objects | *Dielectric Objects | ||
Line 14: | Line 14: | ||
In this tutorial you will construct structures involving metal and dielectric objects in [[EM.Ferma]], visualize their electric field and potential and calculate their capacitance. Specifically, you will build and analyze air-filled and dielectric-filled parallel plate capacitors. | In this tutorial you will construct structures involving metal and dielectric objects in [[EM.Ferma]], visualize their electric field and potential and calculate their capacitance. Specifically, you will build and analyze air-filled and dielectric-filled parallel plate capacitors. | ||
+ | |||
+ | [[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_Lesson2.zip Download projects related to this tutorial lesson]''' | ||
== Getting Started == | == Getting Started == | ||
Line 80: | Line 86: | ||
Activate each PEC group by right-clicking on their name in the navigation tree and selecting '''Activate''' from the contextual menu. Then draw two rectangle objects with the location and dimensions given in the table below. | Activate each PEC group by right-clicking on their name in the navigation tree and selecting '''Activate''' from the contextual menu. Then draw two rectangle objects with the location and dimensions given in the table below. | ||
+ | |||
+ | <table> | ||
+ | <tr> | ||
+ | <td> | ||
+ | [[Image:Ferma_L2_Activate.png|thumb|left|480px|Activating PEC group in the navigation tree.]] | ||
+ | </td> | ||
+ | </tr> | ||
+ | </table> | ||
{| class="wikitable" | {| class="wikitable" | ||
Line 106: | Line 120: | ||
|} | |} | ||
− | To draw a rectangle, click the '''Rectangle Strip''' button the Object Toolbar or select the menu item <b>Object → Surface → Rectangle Strip</b>. | + | To draw a rectangle, click the '''Rectangle Strip''' [[Image:RectangleStripIconx.png]] button of the Object Toolbar or select the menu item <b>Object → Surface → Rectangle Strip</b>. |
<table> | <table> | ||
Line 138: | Line 152: | ||
== Computational Domain, Boundary Conditions & Observables == | == Computational Domain, Boundary Conditions & Observables == | ||
− | You saw in the previous tutorial lesson that [[EM.Ferma]] automatically places a domain box around your objects. By default, this is a PEC box consisting of six perfect metal plates held at zero potential. In other words, the Dirichlet boundary conditions are applied to your boundary-value problem. The default domain box is placed at a certain offset distance from the largest bounding box enclosing all of your geometric objects in the project workspace. [[EM.Ferma]]'s default domain offset is 20 project units. You can change this offset value individually along the six axial directions. Open the Domain Settings Dialog by clicking the <b>Domain | + | You saw in the previous tutorial lesson that [[EM.Ferma]] automatically places a domain box around your objects. By default, this is a PEC box consisting of six perfect metal plates held at zero potential. In other words, the Dirichlet boundary conditions are applied to your boundary-value problem. The default domain box is placed at a certain offset distance from the largest bounding box enclosing all of your geometric objects in the project workspace. [[EM.Ferma]]'s default domain offset is 20 project units. You can change this offset value individually along the six axial directions. Open the Domain Settings Dialog by clicking the <b>3D Static Domain</b> [[Image:fdtd_domainsettings2.png]] button of Simulate Toolbar or by selecting the menu item <b>Simulate → Computational Domain → Domain Settings…</b> or using the keyboard shortcut {{key|Ctrl+A}}). |
In this project, the bottom capacitor plate is held at zero potential, the same as the six domain boundary walls. On the other hand, the top capacitor plate has a potential of 1V. If the domain walls are placed too close to the top plate, strong field lines will establish between them. You need to place the domain walls far enough from the top plate to make them "invisible". Set the six offset values according to the table below: | In this project, the bottom capacitor plate is held at zero potential, the same as the six domain boundary walls. On the other hand, the top capacitor plate has a potential of 1V. If the domain walls are placed too close to the top plate, strong field lines will establish between them. You need to place the domain walls far enough from the top plate to make them "invisible". Set the six offset values according to the table below: | ||
Line 232: | Line 246: | ||
<tr> | <tr> | ||
<td> | <td> | ||
− | [[Image:Ferma L2 Fig9.png|thumb|left|640px|The electric field distribution on the vertical | + | [[Image:Ferma L2 Fig9.png|thumb|left|640px|The electric field distribution on the vertical Sensor_2 (YZ) plane.]] |
</td> | </td> | ||
</tr> | </tr> | ||
<tr> | <tr> | ||
<td> | <td> | ||
− | [[Image:Ferma L2 Fig10.png|thumb|left|640px|The electric potential distribution on the vertical | + | [[Image:Ferma L2 Fig10.png|thumb|left|640px|The electric potential distribution on the vertical Sensor_2 (YZ) plane.]] |
</td> | </td> | ||
</tr> | </tr> | ||
</table> | </table> | ||
− | Also, open the Data Manager and plot the data file | + | Also, open the Data Manager and plot the data file. As you should have expected, the results show a quite uniform electric field at the center of the parallel plate capacitor. |
== Verifying the Simulation Results == | == Verifying the Simulation Results == | ||
Line 258: | Line 272: | ||
<math> E_z = -\nabla \Phi = -\frac{d\Phi}{dz} = -500 V/m </math> | <math> E_z = -\nabla \Phi = -\frac{d\Phi}{dz} = -500 V/m </math> | ||
− | Open the Data Manager and plot the data files “Sensor_1_X_ETotal.DAT” and “Sensor_1_X_EPotential.DAT” | + | Open the Data Manager and plot the data files “Sensor_1_X_ETotal.DAT” and “Sensor_1_X_EPotential.DAT” as shown below. You can clearly see a wide region of constant potential Ψ = 0.5V and of constant electric field strength |E| = 500V/m under the two metal plates. |
<table> | <table> | ||
Line 316: | Line 330: | ||
<table> | <table> | ||
+ | <tr> | ||
<td> | <td> | ||
[[Image:Ferma L2 Fig20.png|thumb|left|720px|The voltage line and flux box defined for the parallel plate capacitor structure.]] | [[Image:Ferma L2 Fig20.png|thumb|left|720px|The voltage line and flux box defined for the parallel plate capacitor structure.]] | ||
Line 322: | Line 337: | ||
</table> | </table> | ||
− | + | Freeze the physical objects to see other objects hidden behind, beneath or inside it. In the figure below, the two plates have been shown in freeze state. | |
<table> | <table> | ||
Line 337: | Line 352: | ||
<tr> | <tr> | ||
<td> | <td> | ||
− | [[Image: | + | [[Image:Ferma_L2_Fig22size.png|thumb|left|480px|Increasing the mesh resolution in the mesh settings dialog.]] |
</td> | </td> | ||
</tr> | </tr> | ||
</table> | </table> | ||
− | Run a new electrostatic analysis of the structure with the new field integral observable. At the end of the simulation, open the data manager and select the data file "FI_1_Capacitance.DAT". You can use the {{key|View}} button of the data manager to view the contents of this file in the data manager's built-in spreadsheet. If you do so, you will see a zero value. This is because the value of the capacitance is extremely small. While in the data manager dialog, select "FI_1_Capacitance.DAT" and click the {{key|Notepad}} button. This opens the Windows Notepad, where you will see a value of "1.49196e-011" or 14.92pF. | + | Run a new electrostatic analysis of the structure with the new field integral observable. At the end of the simulation, open the data manager and select the data file "FI_1_Capacitance.DAT". You can use the {{key|View}} button of the data manager to view the contents of this file in the data manager's built-in spreadsheet. If you do so, you will see a zero value. This is because the value of the capacitance is extremely small. While in the data manager dialog, select "FI_1_Capacitance.DAT" and click the {{key|Notepad}} or {{key|Open}} button. This opens the Windows Notepad, where you will see a value of "1.49196e-011" or 14.92pF. |
Using the well-known parallel-plate capacitance formula, you will get: | Using the well-known parallel-plate capacitance formula, you will get: | ||
Line 357: | Line 372: | ||
<tr> | <tr> | ||
<td> | <td> | ||
− | [[Image:Ferma L2 Fig13.png|thumb|left| | + | [[Image:Ferma L2 Fig13.png|thumb|left|600px|The new dielectric dialog.]] |
</td> | </td> | ||
</tr> | </tr> | ||
Line 375: | Line 390: | ||
<tr> | <tr> | ||
<td> | <td> | ||
− | [[Image:Ferma L2 Fig15.png|thumb|left| | + | [[Image:Ferma L2 Fig15.png|thumb|left|600px|The property dialog of the dielectric material group with Mica selected.]] |
</td> | </td> | ||
</tr> | </tr> | ||
Line 388: | Line 403: | ||
</table> | </table> | ||
− | With the Box tool selected, click on a blank space in the project workspace and drag the mouse to draw the rectangular base of your box object. A property dialog pops up at the lower right corner of your screen. As you drag the mouse, you will see that the X-dimension and Y-dimension of your new object continuously change. When the base reaches the desired size or something close to that, click the mouse. Next, you have to give the right height to your box. Drag the mouse upward until you reach the desired height (Z-dimension). Then, click once more (i.e. drop the mouse). At this time, the drawing of your box is complete. You can always enter values for the dimensions or LCS coordinates at any time. While the new Dielectric_1 group is active, draw a box object according to the table below: | + | With the Box tool selected, click on a blank space in the project workspace and drag the mouse to draw the rectangular base of your box object. A property dialog pops up at the lower right corner of your screen. As you drag the mouse, you will see that the X-dimension and Y-dimension of your new object continuously change. When the base reaches the desired size or something close to that, click the mouse. Next, you have to give the right height to your box. Drag the mouse upward until you reach the desired height (Z-dimension). Then, click once more (<i>i.e.</i> drop the mouse). At this time, the drawing of your box is complete. You can always enter values for the dimensions or LCS coordinates at any time. While the new Dielectric_1 group is active, draw a box object according to the table below: |
{| class="wikitable" | {| class="wikitable" | ||
Line 435: | Line 450: | ||
<p> </p> | <p> </p> | ||
− | [[Image:Back_icon.png| | + | [[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]]''' |
Latest revision as of 19:26, 22 July 2019
Contents
- 1 What You Will Learn
- 2 Getting Started
- 3 Creating the Metal Plate Objects
- 4 Computational Domain, Boundary Conditions & Observables
- 5 Running an Electrostatic Analysis of Your Capacitor
- 6 Verifying the Simulation Results
- 7 Defining a Field Integral Observable for Computing the Capacitance
- 8 Adding a Dielectric Filling to Your Capacitor
What You Will Learn
In this tutorial you will construct structures involving metal and dielectric objects in EM.Ferma, visualize their electric field and potential and calculate their capacitance. Specifically, you will build and analyze air-filled and dielectric-filled parallel plate capacitors.
Back to EM.Ferma Tutorial Gateway
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:
Creating the Metal Plate Objects
in EM.Ferma, perfect electric conductors (PEC) represent equipotential surfaces, bodies or objects. A PEC group has a color and a "Fixed Potential" parameter expressed in Volts. All the objects belonging to the same PEC group have the same color and the same voltage. Therefore, a PEC object with a nonzero voltage can be used as a source in your electrostatic problem. In that case, EM.Ferma solve the Laplace equation in a charge-free medium as opposed to the Poisson equation which you solved earlier in Tutorial Lesson 1.
Since the two metal plates of a capacitor have different potentials, they must belong to two different PEC groups. Create two PEC groups in the navigation tree according to the table below. To do so, right-click on Fixed-Potential PEC objects and select Insert New PEC… from the contextual menu. Set the specified voltage value in the PEC dialog.
Material Group | Material Type | Color | Fixed Potential |
---|---|---|---|
PEC_1 | PEC | Rosy red (default) | 0V |
PEC_2 | PEC | Orange | 1V |
Activate each PEC group by right-clicking on their name in the navigation tree and selecting Activate from the contextual menu. Then draw two rectangle objects with the location and dimensions given in the table below.
Part | Object Type | Material Type | Dimensions | Coordinates | Rotation Angles |
---|---|---|---|---|---|
Rect_Strip_1 | Rectangle Strip | PEC_1 | 50mm × 50mm | (0, 0, -1mm) | (0°, 0°, 0°) |
Rect_Strip_2 | Rectangle Strip | PEC_2 | 50mm × 50mm | (0, 0, 1mm) | (0°, 0°, 0°) |
To draw a rectangle, click the Rectangle Strip button of the Object Toolbar or select the menu item Object → Surface → Rectangle Strip.
With the rectangle strip tool selected, click on a blank space in the project workspace and drag the mouse to draw the planar rectangle object. A property dialog pops up at the lower right corner of your screen. As you drag the mouse, you will see that the X-dimension and Y-dimension of your new object continuously change. You can enter any values for the side dimensions or the X, Y, Z coordinates of the center of local coordinate system (LCS) at any time.
Your physical structure should look like this:
Computational Domain, Boundary Conditions & Observables
You saw in the previous tutorial lesson that EM.Ferma automatically places a domain box around your objects. By default, this is a PEC box consisting of six perfect metal plates held at zero potential. In other words, the Dirichlet boundary conditions are applied to your boundary-value problem. The default domain box is placed at a certain offset distance from the largest bounding box enclosing all of your geometric objects in the project workspace. EM.Ferma's default domain offset is 20 project units. You can change this offset value individually along the six axial directions. Open the Domain Settings Dialog by clicking the 3D Static Domain button of Simulate Toolbar or by selecting the menu item Simulate → Computational Domain → Domain Settings… or using the keyboard shortcut Ctrl+A).
In this project, the bottom capacitor plate is held at zero potential, the same as the six domain boundary walls. On the other hand, the top capacitor plate has a potential of 1V. If the domain walls are placed too close to the top plate, strong field lines will establish between them. You need to place the domain walls far enough from the top plate to make them "invisible". Set the six offset values according to the table below:
Location | Offset | Type |
---|---|---|
-X | 20mm | PEC |
+X | 20mm | PEC |
-Y | 20mm | PEC |
+Y | 20mm | PEC |
-Z | 10mm | PEC |
+Z | 30mm | PEC |
Next, right-click on the Boundary Conditions item under the "Computational Domain" section of the navigation tree and open its settings dialog from the contextual menu. As you can see from the figure below, all the six walls of the domain box have a zero voltage.
Finally, define two field sensor observables: the first Z-directed and called "Sensor_1", and the second X-directed and called "Sensor_2". Leave the default center coordinates of (0, 0, 0) for both sensor planes.
Running an Electrostatic Analysis of Your Capacitor
At this time, your project is ready for simulation. Open the mesh settings dialog and change the Fixed Mesh Cell Size to Δx = Δy = Δz = 0.5mm, i.e. a uniform mesh along all the three principal directions. The mesh resolution was chosen so that four cells are placed between the two metal plates across the gap region.
Run an electrostatic "Analysis" of your structure and visualize the electric field and electric potential results on both sensor planes.
Also, open the Data Manager and plot the data file. As you should have expected, the results show a quite uniform electric field at the center of the parallel plate capacitor.
Verifying the Simulation Results
Here we review a simple, approximate, analytical solution of the capacitor problem. Given the large size of the metal plates and their X- and Y-symmetries, one can assume that the electric potential is only a function of the z coordinate at the center. The Laplace equation can then be written as:
[math]\nabla^2 \Phi = \frac{d^2\Phi}{dz^2} = 0 [/math]
subject to the boundary conditions: Φ(z=-0.001m) = 0V and Φ(z = 0.001m) = 1V. The solution to this equation is:
[math]\Phi (z) = Az+B = 500z +0.5 [/math]
The electric field is then found to be:
[math] E_z = -\nabla \Phi = -\frac{d\Phi}{dz} = -500 V/m [/math]
Open the Data Manager and plot the data files “Sensor_1_X_ETotal.DAT” and “Sensor_1_X_EPotential.DAT” as shown below. You can clearly see a wide region of constant potential Ψ = 0.5V and of constant electric field strength |E| = 500V/m under the two metal plates.
Defining a Field Integral Observable for Computing the Capacitance
EM.Ferma lets you define a number of post-processed observables called "Field Integrals", which are calculated by integrating field quantities over paths, loops, open or closed surfaces, etc. In this tutorial lesson, you will define a "Capacitance" observable. The capacitance is defined as:
[math] C = \frac{\Phi_E}{V} = \frac{\int\int_{S_o} \epsilon \mathbf{E(r)} . \mathbf{ds} } {\int_{\Gamma} \mathbf{E(r)} . \mathbf{dl} } [/math]
where Γ represents an open curve (path), So represents a closed surface like a box, E(r) is the electric field, ΦE is the electric flux passing through So, and V is the voltage along the path Γ. In EM.Ferma, So is defined as a box characterized by two opposite corner points, and Γ is defined a line segment characterized by two endpoints.
To define the capacitance observable, right-click on the Static Field Integrals item in the "Observables" section of the navigation tree and select Insert New Field Integral… from the contextual menu. At the top of the field integral dialog, select Capacitance from the drop-down list labeled Integral Type. Also, choose turquoise as the color of the observable.
This dialog has two more sections, one for defining the domain of the electric flux surface integral, and the other for defining the path of the voltage line integral. Define the two integral domains as follows:
Integral Domain | Domain Type | Point 1 Coordinates | Point 2 Coordinates |
---|---|---|---|
Electric Flux Box | Closed Surface | (-30mm, -30mm, 0) | (-30mm, -30mm, 5mm) |
Voltage Line | Open Curve Path | (0, 0, -1mm) | (0, 0, 1mm) |
For this part, delete the two field sensor observables "Sensor_1" and "Sensor_2". Right-click on the name of each one in the navigation tree and select Delete from the contextual menu. Or right-click on the topmost "Near-Field Sensors" item in the navigation tree and select Delete All from the contextual menu. Note that the flux box has been defined to enclose the upper metal plate only. The voltage path has been defined as a vertical line joining the centers of the two metal plates.
Freeze the physical objects to see other objects hidden behind, beneath or inside it. In the figure below, the two plates have been shown in freeze state.
Before running a new simulation, open the mesh settings dialog once again and change the mesh cell size along the Z-axis to 0.2mm. This will generate a total of 10 cells across the the gap region of the capacitor and will provide a much better resolution.
Run a new electrostatic analysis of the structure with the new field integral observable. At the end of the simulation, open the data manager and select the data file "FI_1_Capacitance.DAT". You can use the View button of the data manager to view the contents of this file in the data manager's built-in spreadsheet. If you do so, you will see a zero value. This is because the value of the capacitance is extremely small. While in the data manager dialog, select "FI_1_Capacitance.DAT" and click the Notepad or Open button. This opens the Windows Notepad, where you will see a value of "1.49196e-011" or 14.92pF.
Using the well-known parallel-plate capacitance formula, you will get:
[math] C = \epsilon_0 \frac{A}{d} = 8.8542\times 10^{-12} \left( \frac{0.05^2}{0.002} \right) = 1.1067\times 10^{-11} F = 11.06pF [/math]
In the above formula, A is the surface area of the capacitor plates and d is the spacing between them. Note that the discrepancy is due to the assumption by the formula that the electric field is constant, uniform and strictly confined between the two parallel plates of the capacitor. However, from the field distribution plots you got earlier, it is evident that this assumption might not be very accurate.
Adding a Dielectric Filling to Your Capacitor
Next, you will add a cubic dielectric filler block that will be sandwiched between the two metal plates of your capacitor structure. To do so, you need to create a dielectric material group. Right-click on the Dielectric/Magnetic Materials item in the navigation tree and select Insert New Dielectric... from the contextual menu. This opens up the New Dielectric dialog. You can choose a name for your material and assign material properties like permittivity.
Or you can use EM.Ferma's material database. Click the Material button to open the Materials List. For this project, choose "Mica" with a permittivity of εr = 5.4. In the Materials Dialog, you can simply type in the letter "M" and it will take you to the materials whose name starts with that letter.
To draw a box object, click the Box button of the Object Toolbar or select the menu item Object → Solid → Box.
With the Box tool selected, click on a blank space in the project workspace and drag the mouse to draw the rectangular base of your box object. A property dialog pops up at the lower right corner of your screen. As you drag the mouse, you will see that the X-dimension and Y-dimension of your new object continuously change. When the base reaches the desired size or something close to that, click the mouse. Next, you have to give the right height to your box. Drag the mouse upward until you reach the desired height (Z-dimension). Then, click once more (i.e. drop the mouse). At this time, the drawing of your box is complete. You can always enter values for the dimensions or LCS coordinates at any time. While the new Dielectric_1 group is active, draw a box object according to the table below:
Part | Object Type | Material Type | Dimensions | Coordinates | Rotation Angles |
---|---|---|---|---|---|
Box_1 | Box | Dielectric_1 | 50mm × 50mm × 2mm | (0, 0, -1mm) | (0°, 0°, 0°) |
Your dielectric-filled parallel plate capacitor should now look like this:
With the same domain box, mesh settings and observables from the previous part, run a new electrostatic analysis of your structure. The electric field and potential distributions will be very similar to the case of air-filled capacitor. That is why you didn't initiate new field sensor observables for this part.
Open the data manager and view the contents of the data file "FI_1_Capacitance.DAT" in Windows Notepad. This time you will see that the value of the capacitance has increased to "6.79847e-011" or 67.98pF. The parallel-plate capacitance formula with a dielectric spacer would yield:
[math] C = \epsilon_r \epsilon_0 \frac{A}{d} = 5.4 \times 8.854\times 10^{-12} \left( \frac{0.05^2}{0.002} \right) = 5.976\times 10^{-11} F = 59.76pF [/math]
The better agreement in this case (12% deviation) is clearly due to the fact that the dielectric filler helps increased confinement of the electric field in the space between the parallel plates.