Applied Electromagnetics Using QuickField and MATLAB

Chapter 7 - Transient Magnetics: Time-dependent Maxwell’s Equations

In This Chapter

  • Time-dependent Maxwell’s Equations
  • Finite Difference Time-domain Method
  • Vector Potential Formulation
  • Time-dependent Magnetics in QuickField
  • Pulsed Eddy Current Nondestructive Testing Applications
  • Maxwell’s Equations with Monopole Source Term

In the previous chapter we considered the special case of Maxwell’s equations with time-harmonic sources. If the field sources have arbitrary time dependence, then we must work with the time-dependent curl and divergence equations

where Je = σE + Jsource and (μ, ε) = ( μr μ0 r ε0 ).

Current Sheet Above a Conducting Half-Space

The simplest example of transient current induction is that produced by a magnetic field tangential to a conducting half-space that is switched on at t = 0.

This field can be thought of as being produced by an infinite sheet with a uniform surface current density that is parallel to the conducting half-space. Neglecting the displacement current we have

inside the metal and taking the curl of this equation and using

we obtain

For an external field B(t) = B0θ(t)j, the differential equation with initial conditions becomes

where θ(t) is zero for t < 0 and 1 for t ≥ 0. Taking the Laplace transform of both the differential equation and initial conditions gives

where By(0) = 0. The solution to this equation is

Applying the boundary conditions we have A = 0 and giving

Performing the inverse Laplace transform in MATLAB

gives the field in the time-domain

where the complementary error function is given by

The eddy current is given by μ0Jx= -∂By/∂z so that using the identity

finally gives the eddy current

Equation 7.5 may be solved numerically using an explicit finite difference numerical method

where the spatial and time derivatives are evaluated using central difference and forward difference schemes, respectively. The magnetic field is then advanced in time at each grid point according to

This explicit method is stable for Δt μ0σΔz2/2. Once the field is known, the current density may be calculated according to

Implicit methods, such as the Crank-Nicholson method, may also be used for greater stability and accuracy. The following MATLAB program calculates the penetration of magnetic fields and the induced currents inside the metal:

The output of this program is shown in Figure 7.1. Axis labels and the figure perspectives were adjusted in the window show.

FIGURE 7.1 MATLAB simulation of the spatial and time-dependent current (left) and magnetic field (right) distributions inside a conducting half-space.

 

UNLIMITED FREE
ACCESS
TO THE WORLD'S BEST IDEAS

SUBMIT
Already a GlobalSpec user? Log in.

This is embarrasing...

An error occurred while processing the form. Please try again in a few minutes.

Customize Your GlobalSpec Experience

Category: Electric Brakes
Finish!
Privacy Policy

This is embarrasing...

An error occurred while processing the form. Please try again in a few minutes.