Partial Differential Equation Toolbox | ![]() ![]() |
AC Power Electromagnetics
AC power electromagnetics problems are found when studying motors, transformers and conductors carrying alternating currents.
Let us start by considering a homogeneous dielectric, with coefficient of dielectricity and magnetic permeability µ, with no charges at any point. The fields must satisfy a special set of the general Maxwell's equations; see [2]:
In the absence of current, we can eliminate H from the first set and E from the second set and see that both fields satisfy wave equations with wave speed
We move on to studying a charge-free homogeneous dielectric, with coefficient of dielectrics , magnetic permeability µ, and conductivity
. The current density then is
and the waves are damped by the Ohmic resistance,
The case of time harmonic fields is treated by using the complex form, replacing
The plane case of this PDE Toolbox mode has Ec = (0,0,Ec) and J = (0,0,Jejt), and the magnetic field
The scalar equation for Ec becomes
This is the equation used by the PDE Toolbox in the AC Power Electromagnetics mode. It is a complex Helmholtz's equation, describing the propagation of plane electromagnetic waves in imperfect dielectrics and good conductors ( »
). A complex permittivity
ccan be defined as
c =
-j
/
. The conditions at material interfaces with abrupt changes of
and µ are the natural ones for the variational formulation and need no special attention.
The PDE parameters that have to be entered into the PDE Specification dialog box are the angular frequency , the magnetic permeability µ, the conductivity
, and the coefficient of dielectricity
.
The boundary conditions associated with this mode are a Dirichlet boundary condition, specifying the value of the electric field Ec on the boundary, and a Neumann condition, specifying the normal derivative of Ec. This is equivalent to specifying the tangential component of the magnetic field H:
Interesting properties that can be computed from the solution -- the electric field E -- are the current density J = E and the magnetic flux density
The electric field E, the current density J, the magnetic field H and the magnetic flux density B are available for plots. Additionally, the resistive heating rate
is also available. The magnetic field and the magnetic flux density can be plotted as vector fields using arrows.
Example
The example demonstrates the skin effect when AC current is carried by a wire with circular cross section. The conductivity of copper is 57 · 106, and the permeability is 1, i.e., µ = 410-7. At the line frequency (50 Hz) the
2
-term is negligible.
Due to the induction, the current density in the interior of the conductor is smaller than at the outer surface where it is set to JS = 1, a Dirichlet condition for the electric field, Ec = 1/. For this case an analytical solution is available,
R is the radius of the wire, r is the distance from the centerline, and J0(x) is the first Bessel function of zeroth order.
Using the Graphical User Interface
Start the pdetool
GUI and select the application mode AC Power Electromagnetics. Draw a circle with radius 0.1 to represent a cross section of the conductor, and proceed to the boundary mode to define the boundary condition. Use the Select All option to select all boundaries and enter 1/57E6
into the r edit field in the Boundary Condition dialog box to define the Dirichlet boundary condition (E = J/).
Open the PDE Specification dialog box and enter the PDE parameters. The angular frequency = 2 ·
· 50.
Initialize the mesh and solve the equation. Due to the skin effect, the current density at the surface of the conductor is much higher than in the conductor's interior. This is clearly visualized by plotting the current density J as a 3-D plot. To improve the accuracy of the solution close to the surface, you need to refine the mesh. Open the Solve Parameters dialog box and check the Adaptive mode option. Also, set the maximum numbers of triangles to Inf
, the maximum numbers of refinements to 1, and use the triangle selection method that picks the worst triangles. Recompute the solution several times. Each time the adaptive solver refines the area with the largest errors. The number of triangles is printed on the command line. The mesh below is the result of successive adaptions and contains 1548 triangles.
The solution of the AC power electromagnetics equation is complex. The plots show the real part of the solution (a warning message is issued), but the solution vector, which can be exported to the main workspace, is the full complex solution. Also, you can plot various properties of the complex solution by using the user entry. imag(u)
and abs(u)
are two examples of valid user entries.
The skin effect is an AC phenomenon. Decreasing the frequency of the alternating current results in a decrease of the skin effect. Approaching DC conditions, the current density is close to uniform (experiment using different angular frequencies).
The Current Density in an AC Wire
![]() | Magnetostatics | Conductive Media DC | ![]() |