1. Mathematical and Physical Foundations of Electromagnetic Fields

Introduction and Context

Electromagnetic field theory provides the conceptual and mathematical basis for electrical engineering. Power systems, electrical machines, antennas, microwave circuits, wireless communication, sensors, electromagnetic compatibility, optical devices, and computational electromagnetics all rely on the same physical language: fields defined in space and time and related through differential and integral laws.

At postgraduate level, electromagnetic theory should not be approached as a collection of formulas. Its strength lies in the connection between physical interpretation, coordinate geometry, vector calculus, material models, and experimentally measurable quantities. A correct solution therefore requires more than manipulating equations. It requires selecting an appropriate coordinate system, identifying the field quantities involved, checking dimensions, applying the correct integral theorem, and interpreting boundary and material conditions.

This chapter develops those foundations. It begins with coordinate descriptions and vector operations, proceeds to integral theorems and generalized functions, and then introduces complex notation and phasors for sinusoidal steady state. It concludes by relating field quantities to constitutive models, dimensions, applications, and professional engineering practice.

Learning Objectives

After studying this chapter, the reader should be able to:

1. Select Cartesian, cylindrical, or spherical coordinates for an electromagnetic problem.

2. Interpret scalar and vector fields and evaluate gradient, divergence, curl, and Laplacian operations.

3. Apply the line, surface, and volume integral forms of electromagnetic relations.

4. Use the divergence theorem and Stokes' theorem correctly, including their physical meaning.

5. Explain the role of Dirac delta functions and idealized sources.

6. Manipulate complex quantities and phasors while maintaining a consistent time convention.

7. Distinguish electric field, magnetic field, flux density, and source quantities.

8. Use dimensional analysis to test electromagnetic equations.

9. Select and assess constitutive models for different materials.

10. Recognize practical, numerical, ethical, and professional implications of field modeling.

Space, Coordinates, and Field Descriptions

An electromagnetic field assigns a physical quantity to every point in a region of space, possibly at every instant of time. A scalar field, such as electric potential V(r,t) or temperature, has magnitude only. A vector field, such as electric field intensity E(r,t) or magnetic flux density B(r,t), has magnitude and direction.

The position vector in Cartesian coordinates is

r = x a_x + y a_y + z a_z,

where a_x, a_y, and a_z are mutually perpendicular unit vectors. Cartesian coordinates are natural for rectangular conductors, planar interfaces, and uniform structures. The differential displacement is

d l = a_x dx + a_y dy + a_z dz,

and the volume element is dV = dx dy dz.

Cylindrical coordinates (rho, phi, z) are appropriate when the geometry has rotational symmetry about an axis. The position vector is

r = rho a_rho + z a_z,

although the azimuthal coordinate is represented through the direction of a_rho. The line and volume elements are

d l = a_rho d rho + a_phi rho d phi + a_z dz,

and dV = rho d rho d phi dz.

Spherical coordinates (r, theta, phi) are useful for isolated sources, spheres, radiation, and central-force problems. Here theta is measured from the positive z-axis and phi is the azimuthal angle. The differential displacement is

d l = a_r dr + a_theta r d theta + a_phi r sin(theta) d phi,

and dV = r squared sin(theta) dr d theta d phi.

The unit vectors in cylindrical and spherical systems vary with position. This is an important distinction from Cartesian coordinates. For example, a_rho changes direction as phi changes. Consequently, derivatives of vector fields in curvilinear coordinates must include both changes in components and changes in unit vectors. Omitting this fact is a frequent source of error.

Vector Calculus as the Language of Fields

The dot product of two vectors gives a scalar:

A dot B = |A||B| cos(alpha).

It measures the component of one vector along another and appears in flux calculations. The cross product gives a vector perpendicular to both:

A cross B = |A||B| sin(alpha) a_n.

It appears in torque, force, and electromagnetic power-flow expressions.

Gradient

For a scalar field V, the gradient is

∇V = a_x partial V/partial x + a_y partial V/partial y + a_z partial V/partial z.

The gradient points in the direction of maximum increase of V. Its magnitude is the maximum spatial rate of change. In electrostatics, E = -∇V; therefore, the electric field points toward decreasing electric potential.

Divergence

For a vector field A, divergence is defined in Cartesian coordinates by

∇ dot A = partial A_x/partial x + partial A_y/partial y + partial A_z/partial z.

Divergence measures the net outward flow per unit volume in the limit of a shrinking volume. Positive divergence indicates a local source; negative divergence indicates a local sink. In cylindrical and spherical coordinates, scale factors must be included. For example, in cylindrical coordinates,

∇ dot A = (1/rho) partial(rho A_rho)/partial rho + (1/rho) partial A_phi/partial phi + partial A_z/partial z.

Curl

Curl measures local circulation. In Cartesian coordinates,

∇ cross A = determinant | a_x a_y a_z; partial/partial x partial/partial y partial/partial z; A_x A_y A_z |.

A field with zero curl in a simply connected region is irrotational and can often be represented as the gradient of a scalar potential. The magnetic field, however, may have nonzero curl where current or changing electric flux is present.

Laplacian

For a scalar field,

∇ squared V = ∇ dot ∇V.

In a charge-free homogeneous electrostatic region, the potential satisfies Laplace's equation, ∇²V = 0. In a region containing volume charge density rho_v, it satisfies Poisson's equation, ∇²V = -rho_v/epsilon. These equations are central to analytical and numerical field solutions.

Integral Theorems and Their Physical Meaning

Differential expressions describe behavior at a point, whereas integral expressions describe accumulated behavior over a path, surface, or volume. The two forms are connected by theorems of vector calculus.

The divergence theorem states

∫_V (∇ dot A) dV = ∮_S A dot dS,

where S is the closed surface bounding volume V. It converts volume divergence into net outward flux. Gauss's law is an electromagnetic application: the total electric flux through a closed surface equals the enclosed charge.

Stokes' theorem states

∫_S (∇ cross A) dot dS = ∮_C A dot dl,

where C is the closed contour bounding S. It converts the curl through a surface into circulation around its boundary. Faraday's law and the Ampere-Maxwell law can be written in this form.

The line integral of a vector field along a path is ∫ A dot dl. In electrostatics, the circulation of E around every closed path is zero, so the electrostatic field is conservative. The surface integral ∫ A dot dS measures flux through a surface. Orientation matters: reversing the surface normal changes the sign of the flux.

A common analytical mistake is applying the divergence theorem to an open surface or Stokes' theorem without identifying a compatible boundary curve. The geometry of the integration domain is part of the physical problem, not merely a mathematical detail.

Sources, Distributions, and Idealization

Real charge and current distributions may be continuous, localized, or concentrated in very small regions. Volume charge density is denoted by rho_v and has units C/m³; surface charge density by rho_s, with units C/m²; and line charge density by rho_l, with units C/m.

The Dirac delta distribution represents an idealized concentrated source. In one dimension,

∫ delta(x - x_0) f(x) dx = f(x_0).

A point charge q at r_0 can be represented by

rho_v(r) = q delta(r - r_0),

where the three-dimensional delta function has units m⁻³. The delta function is not an ordinary function with a finite peak; it is a distribution defined through its integral action. It enables compact source descriptions and appears in Green-function methods, impulse responses, scattering theory, and computational formulations.

Idealization is useful but must be interpreted carefully. A point charge is a mathematical model whose validity depends on the observation scale. At dimensions comparable to the physical charge distribution, the model may no longer be adequate. Similar cautions apply to perfect conductors, infinitesimally thin current sheets, and perfectly lossless materials.

Electromagnetic Field Quantities

The electric field intensity E is force per unit positive test charge and has units V/m or N/C. Electric flux density D is related to electric charge and has units C/m². In vacuum, D = epsilon_0 E.

The magnetic flux density B has units tesla (T), equivalent to Wb/m². Magnetic field intensity H has units A/m. In vacuum, B = mu_0 H. The distinction between B and H is essential in magnetic materials because their relationship includes magnetization.

Current density J has units A/m². The total current through a surface is I = ∫ J dot dS. Charge conservation is expressed by the continuity equation,

∇ dot J = -partial rho_v/partial t.

This equation states that charge cannot disappear locally; a decrease in charge within a volume must be accompanied by outward current.

Maxwell's equations in differential form are

∇ dot D = rho_v,

∇ dot B = 0,

∇ cross E = -partial B/partial t,

∇ cross H = J + partial D/partial t.

The second equation expresses the absence of isolated magnetic monopoles in classical electromagnetism. The displacement-current term, partial D/partial t, permits the Ampere-Maxwell law to remain compatible with charge conservation and is indispensable in capacitive and radiating systems.

Complex Notation and Phasors

Sinusoidal steady-state analysis becomes simpler when real sinusoids are represented by complex exponentials. With the convention

x(t) = Re{X exp(j omega t)},

where X is a complex phasor, differentiation becomes multiplication by j omega:

dx/dt corresponds to j omega X.

For example, if E(t) = E_0 cos(omega t + phi), its phasor is E = E_0 exp(j phi), provided peak rather than rms magnitude is being used. Phasors eliminate repeated trigonometric manipulation and convert linear time-invariant differential equations into algebraic equations.

The time convention must be declared. Under exp(j omega t), a derivative corresponds to +j omega; under exp(-j omega t), it corresponds to -j omega. Both conventions are valid, but mixing them produces incorrect signs in propagation constants, impedance, and Maxwell equations.

For two phasors A and B, the time-average of their instantaneous product in power calculations is commonly expressed as

P_avg = (1/2) Re{A cross B*},

for peak-valued field phasors. The vector A cross B* is associated with the complex Poynting vector. If rms phasors are used, the factor one-half is omitted. This convention must be stated in engineering reports and software documentation.

Constitutive Models and Material Response

Maxwell's equations alone do not specify how a particular material responds. Constitutive relations provide that information. The simplest linear, isotropic, homogeneous model is

D = epsilon E,

B = mu H,

J = sigma E,

where epsilon is permittivity, mu is permeability, and sigma is conductivity.

In a linear isotropic material, epsilon, mu, and sigma are scalar constants. In anisotropic media they are tensors, so D need not be parallel to E. In dispersive media, the parameters depend on frequency; in the time domain, the response may involve convolution and memory. In nonlinear materials, epsilon or mu depends on field magnitude, and superposition no longer applies.

The relative permittivity is epsilon_r = epsilon/epsilon_0 and relative permeability is mu_r = mu/mu_0. The loss tangent, often used for dielectric materials, compares dissipative and stored electric effects. For a simple dielectric, tan(delta) is related to conductivity and frequency through an effective complex permittivity.

A conducting medium is often characterized by the ratio sigma/(omega epsilon). If this ratio is much greater than one, conduction current dominates displacement current; if it is much less than one, displacement current dominates. This comparison is frequency-dependent and should not be treated as an intrinsic classification independent of operating conditions.

Boundary conditions follow from Maxwell's equations. At an interface, the normal component of D changes according to the surface charge density, while the normal component of B is continuous. The tangential component of E is continuous in the absence of an impulsive magnetic flux, and the tangential component of H changes according to surface current density. These conditions are frequently more useful than solving the full field equations in layered structures.

Dimensional Analysis and Verification

Dimensional analysis is an efficient method for detecting errors. The SI dimensions include:

• E: V/m or N/C
• D: C/m²
• H: A/m
• B: T or Wb/m²
• J: A/m²
• epsilon: F/m
• mu: H/m
• sigma: S/m

For example, epsilon E has dimensions (C/Vm)(V/m) = C/m², confirming that D = epsilon E is dimensionally consistent. In the wave equation, terms such as ∇²E and mu epsilon partial²E/partial t² must have identical dimensions. Dimensional consistency does not prove that an equation is physically correct, but dimensional inconsistency proves that it is wrong.

A sound verification procedure includes checking units, limiting cases, symmetry, boundary conditions, sign conventions, and conservation laws. For instance, a solution for an isolated point charge should show radial symmetry and decrease with distance according to the inverse-square law. A proposed magnetic field with nonzero net flux through a closed surface violates ∇ dot B = 0.

Applied Scenario: Field Modeling of a Shielded Cable

Background

Consider a coaxial cable with an inner conductor of radius a and an outer conducting shield of inner radius b. The cable carries a sinusoidal signal and is used to connect sensitive measurement equipment.

Situation and Problem

The design team must estimate electric field, magnetic field, characteristic impedance, and shielding behavior while deciding whether a lossless or lossy model is appropriate.

Analysis

The cylindrical symmetry suggests cylindrical coordinates. In the dielectric region a < rho < b, the electrostatic potential depends only on rho. Applying Laplace's equation gives a logarithmic potential distribution, and the electric field is radial with magnitude proportional to 1/rho. The voltage between conductors is obtained from the radial line integral of E.

For a current I on the inner conductor and return current on the shield, Ampere's law gives an azimuthal magnetic field in the dielectric proportional to I/rho. The electric and magnetic fields are orthogonal, so their cross product points along the cable axis and represents transmitted power.

The material model determines accuracy. For a low-loss dielectric, a model with real epsilon and negligible sigma may be sufficient for preliminary impedance estimation. At high frequency, dielectric loss, conductor surface resistance, dispersion, connector discontinuities, and shield imperfections must be included. A finite-element or finite-difference solver can use the analytical cylindrical solution as a benchmark.

Approach and Outcome

The team first derives the ideal field using symmetry and integral laws, checks units and boundary conditions, and then introduces losses through complex material parameters. Measurements of attenuation and impedance are compared with the model. Disagreement is investigated through connector geometry and frequency-dependent material data rather than being hidden by arbitrary parameter adjustment.

Learning Points

The scenario illustrates that coordinate selection reduces mathematical complexity, constitutive data controls predictive quality, analytical solutions are valuable verification tools, and a field model should be refined only when the intended decision requires the additional complexity.

Professional Practice, Risks, and Emerging Tools

Electromagnetic analysis is used in insulation coordination, high-voltage design, electromagnetic compatibility, wireless systems, biomedical equipment, electric vehicles, and semiconductor packaging. Engineers increasingly combine symbolic derivation, laboratory measurement, and numerical simulation. Finite-element, finite-difference time-domain, method-of-moments, and multiphysics platforms implement the same field concepts introduced here.

Numerical models introduce their own risks: inadequate mesh resolution, incorrect boundary conditions, poorly characterized materials, numerical dispersion, and mistaken interpretation of complex amplitudes. Software output is not automatically evidence of physical correctness. Mesh-convergence studies, analytical benchmarks, energy checks, and experimental validation are professional best practices.

Machine-learning methods are beginning to assist surrogate modeling, inverse design, parameter extraction, and rapid optimization. They do not remove the need for Maxwell-consistent formulations or reliable training data. An apparently accurate prediction outside the range of validated data may be unsafe.

Engineers must also consider safety and governance. Field exposure, insulation failure, interference with medical or communication equipment, and unintended radiation have consequences beyond numerical error. Designs should comply with applicable technical standards and institutional safety procedures. Material data, measurement uncertainty, assumptions, and software settings should be documented so that results are reproducible and auditable.

Common mistakes include confusing E with D or B with H, using a surface integral where a closed-surface integral is required, ignoring coordinate scale factors, mixing peak and rms phasors, changing time conventions midway, treating a distribution as an ordinary function, and applying a linear isotropic constitutive law to a strongly anisotropic or nonlinear material.

Best-Practice Workflow

A reliable field analysis generally follows this sequence:

1. Define the physical region, sources, materials, frequency range, and required outputs.

2. Exploit geometric symmetry and select coordinates accordingly.

3. Identify whether a scalar, vector, differential, integral, static, transient, or phasor formulation is most appropriate.

4. State source distributions, constitutive relations, boundary conditions, and time conventions.

5. Derive or compute the fields while preserving units.

6. Check conservation laws, limiting behavior, symmetry, and interface conditions.

7. Compare with analytical solutions, measurements, or independent numerical methods.

8. Report assumptions, uncertainties, validity range, and safety implications.

Chapter Summary

Electromagnetic field theory combines physical laws with mathematical descriptions of space, sources, materials, and time variation. Cartesian, cylindrical, and spherical coordinates enable geometry-appropriate analysis. Gradient, divergence, curl, and Laplacian operations describe local field behavior, while line, surface, and volume integrals describe accumulated effects. The divergence theorem and Stokes' theorem connect these descriptions and provide the basis for integral forms of Maxwell's equations.

Charge and current may be modeled as volume, surface, line, or idealized point distributions. Complex notation and phasors simplify sinusoidal analysis, provided the time and amplitude conventions remain consistent. Constitutive relations distinguish material response from universal field laws and may be linear, nonlinear, isotropic, anisotropic, homogeneous, inhomogeneous, or dispersive. Dimensional analysis, boundary conditions, conservation laws, analytical benchmarks, and experimental validation are essential safeguards. These foundations support later study of electrostatics, magnetostatics, wave propagation, transmission lines, antennas, and computational electromagnetics.

Key Terms

Coordinate system: A mathematical scheme for locating points and expressing vector components.

Scalar field: A position- and possibly time-dependent quantity with magnitude only.

Vector field: A quantity having magnitude and direction at each point.

Gradient: The vector spatial derivative of a scalar field.

Divergence: The local net outward flux density of a vector field.

Curl: The local circulation density of a vector field.

Laplacian: The divergence of the gradient, or a second-order spatial differential operator.

Flux: The surface integral of a field through a specified surface.

Divergence theorem: A theorem relating volume divergence to flux through a closed surface.

Stokes' theorem: A theorem relating surface curl to circulation around a boundary curve.

Dirac delta distribution: An idealized distribution defined by its integral sampling property.

Electric field intensity: Force per unit positive charge, measured in V/m.

Electric flux density: Electric flux per unit area, measured in C/m².

Magnetic flux density: Magnetic flux per unit area, measured in tesla.

Constitutive relation: A material-dependent relation between field quantities.

Phasor: A complex amplitude representing a sinusoidal quantity.

Permittivity: A material parameter relating D to E.

Permeability: A material parameter relating B to H.

Continuity equation: The local mathematical expression of charge conservation.

Loss tangent: A measure comparing dissipative and reactive dielectric behavior.

Review Questions

Short Answer Questions

1. Why are cylindrical coordinates advantageous for an ideal coaxial cable?

2. What physical interpretation is assigned to the divergence of a vector field?

3. State the relationship between a sinusoidal derivative and its phasor under the exp(j omega t) convention.

4. Distinguish electric field intensity E from electric flux density D.

5. Why must the dimensions of both sides of an electromagnetic equation agree?

Long Answer Questions

1. Explain gradient, divergence, curl, and Laplacian, including one electromagnetic application of each.

2. Derive the physical significance of the divergence theorem and Stokes' theorem in field theory.

3. Discuss volume, surface, line, and point-source representations of charge.

4. Explain linear, isotropic, anisotropic, nonlinear, and dispersive constitutive models.

5. Present a systematic workflow for validating a numerical electromagnetic simulation.

Analytical / Application Questions

1. A point charge is located at the origin. Which coordinate system is most appropriate, and what symmetry should the electric field exhibit? Explain how this symmetry can be used to select a Gaussian surface.

2. A dielectric has conductivity sigma and operates at angular frequency omega. Explain how the ratio sigma/(omega epsilon) influences whether conduction or displacement current dominates.

3. A simulation reports a nonzero net magnetic flux through a closed surface. Identify the governing law that is violated and propose at least two possible modeling or numerical causes.

22. Multiple Choice Questions with Answers

Question 1

Which coordinate system is most appropriate for the ideal field analysis of a coaxial cable?

A. Cartesian

B. Cylindrical

C. Spherical

D. Bipolar only

Correct Answer: B. Cylindrical

Question 2

The divergence of a vector field most directly represents:

A. Local circulation

B. Maximum scalar increase

C. Net outward flux per unit volume

D. The field's time average

Correct Answer: C. Net outward flux per unit volume

Question 3

Which theorem converts the surface integral of curl over a surface into a line integral around its boundary?

A. Divergence theorem

B. Stokes' theorem

C. Green's theorem only

D. Parseval's theorem

Correct Answer: B. Stokes' theorem

Question 4

Under the phasor convention exp(j omega t), differentiation with respect to time corresponds to multiplication by:

A. 1/(j omega)

B. -j omega

C. j omega

D. omega squared

Correct Answer: C. j omega

Question 5

A point charge represented as a volume charge density requires a three-dimensional delta distribution with units of:

A. m

B. m^-1

C. m^-2

D. m^-3

Correct Answer: D. m^-3

Question 6

Which equation expresses the absence of magnetic monopoles in classical electromagnetism?

A. ∇ dot D = rho_v

B. ∇ dot B = 0

C. ∇ cross E = -partial B/partial t

D. ∇ cross H = J + partial D/partial t

Correct Answer: B. ∇ dot B = 0

Question 7

If sigma/(omega epsilon) is much greater than one, which current generally dominates?

A. Displacement current

B. Conduction current

C. Radiation current only

D. Magnetization current only

Correct Answer: B. Conduction current

Question 8

Which relation is dimensionally consistent for a linear isotropic material?

A. D = mu E

B. B = epsilon H

C. D = epsilon E

D. J = mu E

Correct Answer: C. D = epsilon E

Question 9

For peak-valued phasors, the time-average power-flow density is commonly calculated using:

A. Re{E cross H*}/2

B. Im{E dot H}/2

C. E cross H without conjugation only

D. |E| + |H|

Correct Answer: A. Re{E cross H*}/2

Question 10

A proposed magnetic field has nonzero net flux through a closed surface. The first fundamental check should involve:

A. ∇ dot B = 0

B. ∇ dot E = 0 in all materials

C. J = sigma H

D. D = mu E

Correct Answer: A. ∇ dot B = 0