Introduction

The following tutorial presents an electrostatic application. This example looks at a spherical capacitor formed of a solid conductor sphere, marked with 1 in the figure, and a hollow spherical conductor shell, marked with 3 in the figure, where the region between the conductors is a dielectric material, marked with 2 in the figure. The aim is to reproduce an electric potential distribution using the finite element method and compare the result to an analytical solution of the capacitance. The following sketch shows the geometry of the capacitor in 3D, where the different colors represent the material of each region.

Typically a 3D electrostatic equation is used to create a partial differential equation model for a 3D capacitor. Since the geometry and the boundary conditions, in this case, are rotationally symmetric about any axis, including the -axis, a 2D axisymmetric electrostatic equation can be used to model the capacitor. An axisymmetric model has the advantage that the computational cost in both time and memory is much less than in the case of solving a full 3D model.

Electrostatic equation

Electrostatics can be modeled with a Poisson equation. The 2D axisymmetric Poisson equation is given as:

(-ϵ·r)+(-ϵ)=ρ

where is the electric potential in , is the absolute permittivity in of the material and in is the volume charge density.

The axisymmetric Poisson equation uses a truncated cylindrical coordinate system in 2D with independent variables instead of the cylindrical coordinates . The cylindrical coordinate variable disappears because the system is rotationally symmetric about the -axis. For further reference see DiffusionPDETerm.

The symbols and corresponding units used throughout this tutorial are summarized in the Nomenclature section.

Model

The model consists of a multilayered ball. The inner most layer, the inner conductor, extends to a radius of and has a charge of . The next layer is a dielectric layer and extends from to r₂. This is the layer of interest. The out layer, a second conductor extends from r₂ to r₃. This outer layer is connected to ground and has a potential of .

Both conductors have different but constant electrostatic conditions. This implies that the potential on each surface is also constant, though different. The goal of this simulation is to compute the distribution of the electric potential in in the dielectric layer.

Domain

The domain of the dielectric is an axially symmetric annulus. The external boundaries represent the interfaces between the dielectric and the conductors.

First, the dimensions of the domain are defined, and then an element mesh is created with ToElementMesh.

will be used for a three dimensional visualization.

Parameters setup

Next, the material parameters and the 2D axisymmetric electrostatic equation are set up. These include the relative permittivity , a dimensionless number, and the charge .

The ElectrostaticPDEComponent function can produce the axisymmetric form of the electrostatic equation. To do so the parameter "RegionSymmetry" is set to "Axisymmetric".

Boundary conditions

The boundary conditions considered here include both DirichletCondition and NeumannValue. The ElectricPotentialCondition specifies a ground boundary potential of at the outer boundary of the annulus at . The use of point element markers will be used here to specify where the ElectricPotentialCondition is applied.

The outer boundary has a point element marker of 3 assigned to it. Because the marker assignment algorithm cannot know what markers are supposed to be assigned at corners a unique marker is assigned there and thus allows for a choice whether that corner node should be part of the Dirichlet boundary condition or not.

In the top left corner is a corner point element marker of 8, and in the bottom left corner is a point element marker of 7, where the electric potential condition is also to be applied.

Point element markers are used for DirichletCondition. For NeumannValue boundary element markers are to be used. The ElectricFluxDensityValue specifies a surface charge density in at the boundary of the annulus at , which represents the charge that has the inner conductor.

On a boundary that falls on the symmetry axis, like the -axis in this case, a Neumann 0 value needs to be specified because it is a Symmetry Boundary Condition. As Neumann 0 boundary conditions are the natural default boundary conditions, these can be left out.

Model evaluation

Compute the numeric solution of the PDE with NDSolveValue.

Visualization

To visualize the full 3D solution from the axisymmetric model, one approach is to apply the interpolating function to create a 3D representation of the data. Specifically, a revolution plot of the data can be generated. Using RegionPlot3D, a plot can be created that shows the three-dimensional region where pred is the equation of a sphere with radius , and a color function is specified to map from the 3D data points to the 2D solution function.

Verification

The definition of a capacitance is given as

C =

where is the charge that has the inner conductor sphere and is the potential difference between the conductors.

For a spherical capacitor the capacitance is

C =

where is the relative permittivity of the dielectric layer which is delimited by the radius and .

See [Hayt & Buck, 2012] for further reference. Alternative methods to compute the capacitance can be found in the Capacitance section of the electrostatics monograph.

Another way to verify the 2D axisymmetric model is to compare it with a full 3D model.

The graph shows that the 2D axisymmetric model matches the full 3D model.

Nomenclature

References

1.  Hayt, W. H., & Buck, J. A. (2019). Engineering electromagnetics. McGraw-Hill Education.