"

Physics II

Class Notes, 4/04/99


In an electrostatic field, such as that of the two charges shown below, a test charge can be moved in a direction perpendicular to the electric field without any work being done by or on the charge.

Along an equipotential curve, a test charge can be moved with 0 work involved.

Video Clip #01

Video Clip #02

Video Clip #03

We have written a QBASIC program to plot field lines and equipotential lines.   The program is instructive, and will be outlined here:

We begin by setting up the coordinates of the screen window and asking the user for the charges and the coordinates of the charges, then plotting two circles representing the charges:

WINDOW (-13, -10)-(13, 10)

PRINT "Screen dimensions are -13 to +13 horizontally and -10 to 10 vertically."

INPUT "charge 1: ", q1
INPUT "coordinates of charge 1: ", x1, y1
INPUT "charge 2: ", q2
INPUT "coordinates of charge 2: ", x2, y2


CIRCLE (x1, y1), .05
CIRCLE (x2, y2), .05

We proceed to plot the equipotentials and the fieldlines. 

chois = 2                     

**chois tells the plotting routine whether to plot fieldlines or equipotentials (chois=2 plots equipotentials, chois = 1 plots field lines) **

FOR xx = -13 TO 13    ** xx is the x coordinate on the screen; y is just below 0 **


y = -.00001
x = xx
PSET (x, y)  ** (x,y) is the starting point for the equipotential**
CALL fieldLine ** routine will plot an equipotential**


NEXT

**we set chois=1 so the plotting routine will know to plot fieldlines**
chois = 1
q = q1

** from a starting small angle to `pi get starting points for 24 field lines around q1, then plot the field lines below the y axis from each starting point**
FOR theta = 3.14 / 48 TO 3.15 STEP 3.14 / 24

** the point (x,y) will lie at distance 1 from the charge in the direction -theta from the charge **
x = x1 + COS(theta)
y = y1 - SIN(theta)  ** - sin(theta) puts us below the y axis **
PSET (x, y)
CALL fieldLine  ** this routine will plot the field line or equipotential from the starting point **
NEXT

** do likewise for q2 **
q = q2
FOR theta = 3.14 / 24 TO 3.15 STEP 3.14 / 24 * ABS(q1 / q2)
x = x2 + COS(theta)
y = y2 - SIN(theta)
PSET (x, y)
CALL fieldLine
NEXT

Video Clip #04

The plotting routine is as follows:

SUB fieldLine
SHARED x, y, x1, y1, x2, y2, q1, q2, chois, q  ** these variables are shared with the main routine **

DO UNTIL xq = 999

** the following if-then-else routine cuts off the plotting of the fieldline; you don't need to understand it unless you want to learn to program**

IF choise = 1 THEN
    IF ABS(x) > 13 THEN xq = 999
    IF ABS(y) > 10 THEN xq = 999
ELSE
    IF ABS(x) > 13 THEN xq = 999
    IF ABS(y) > 10 THEN xq = 999
IF y > 0 THEN xq = 999
END IF

** We find the field strength and direction for the field at (x,y) due to q1 **
r1x = x1 - x
r1y = y1 - y  ** x and y distances from (x,y) to charge **
r1 = SQR(r1x ^ 2 + r1y ^ 2) **Phythagorean Theorem for distance to charge **

sine1 = r1y / r1
cosine1 = r1x / r1  ** sine and cosine of angle from charge to field point (x,y) **
fieldStrength1 = q1 / r1 ^ 2  **Coulomb's Law **

** In the same way we find the field strength and direction for the field at (x,y) due to q1 **
r2x = x2 - x
r2y = y2 - y

r2 = SQR(r2x ^ 2 + r2y ^ 2)
sine2 = r2y / r2
cosine2 = r2x / r2
fieldStrength2 = q2 / r2 ^ 2

**We find the magnitude and direction of the total field at (x,y)**

fieldX = fieldStrength1 * cosine1 + fieldStrength2 * cosine2
fieldY = fieldStrength1 * sine1 + fieldStrength2 * sine2  ** x and y components of total field **
fieldStrength = SQR(fieldX ^ 2 + fieldY ^ 2) ** Pythagorean Theorem for vector magnitude**
sine = fieldY / fieldStrength
cosine = fieldX / fieldStrength  ** sine and cosine of angle made by field at (x,y) **

SELECT CASE chois
    CASE 1:

        ** for field line, move (x,y) .1 unit in direction away from charges **
        x = x - .1 * cosine * q / ABS(q)
        y = y - .1 * sine * q / ABS(q) 
    CASE ELSE:

        ** for equipotential, move (x,y) 1 unit in direction perpendicular to field line          (note negative reciprocal slope) **
        x = x + .1 * sine * q / ABS(q)
        y = y - .1 * cosine * q / ABS(q)
END SELECT

LINE -(x, y), 0  ** plot from last point to new point **

You can run the program, which is named FIELDS2, by clicking on its name.

Video Clip #05

Video Clip #06

The figure below depicts the electric field lines emanating from two equal and unlike charges.

The 'flattened circles' are at every point perpendicular to the field lines (note that the aspect ratio in the figure below might result in a slight distortions of these angles).

The figure below depicts two equal and unlike charges.

The figure below shows these two unequal charges as seen from a greater distance.

In the figure below we see equal like charges.

At large distances the field due to these two charges will become more and more nearly radial, approaching that of a doubled charge at the midpoint between the two charges.

In the figure below we see unequal like charges.

The figure below depicts a capacitor consisting of two parallel plates, each with charge density +-`sigma.

The voltage across the plates is the work per unit charge required to move a test charge from one plate to the other.

We can find the charge to voltage ratio, or the capacitance, on a capacitor consisting of parallel plates of area A, with uniform separation d.

Video Clip #07

We thus see that the capacitance of a capacitor is directly proportional to the area of its plates and inversely proportional to the separation between the plates.

If we make a capacitor out of two pieces of aluminum foil, with a piece of paper sandwiched between them to keep the charges apart, we can estimate the capacitance of the system.

= 9 * 10-12 C2 / (N m2) * 600 m

= 5.4 * 10-9 C2 / J

= 5.4 * 10-9 C / (J/C)

= 5.4 * 10-9 C / V

= 5.4 * 10-9 Farad.

"