Archive for September, 2008

Practical Electricity

September 27, 2008

Light Bulbs in Series

Incandescent light bulbs are designed to output light at the rated power at the mains voltage. Light is produced by a tungsten filament which is heated to a very high temperature, about 2000°C. With the black body assumption, we have, using Wien’s displacement law, λT = constant, where λ is the wavelength in the emission that has the highest intensity. Considering that the Sun, with a temperature of about 6000°C has its peak at around 500 nm, the peak for a light bulb lies in the infra red region. This means that most of the power is radiated as heat.

Light Bulbs in Series

Light Bulbs in Series

If two bulbs are connected in series, and the system connected to the mains, voltage available to each is halved. The power input to the bulbs is given by Ohm’s law, V²/R. However, since the power is lower, the temperature is lower too, and hence the resistance is also lower. This means that the power output is not as small as it would be if the resistance variation with temperature was low. The net power would be reduced, but what about the brightness? It will be much lower—as the temperature drops, the black body curve shifts to the right on the intensity wavelength plot, and a very small fraction is now radiated in the visible region. Thus, a 100 W bulb connected with other bulbs such that the power consumed is 25 W will be much dimmer than a 25 W bulb at the rated voltage. The room temperature resistivity of tungsten (300K) is about 5 \times 10^{-8} \Omega m. This means that the value at the operating temperature of 2300 K will be around 50 \times 10^{-8} \Omega m. (This also explains the “inrush current” when a bulb is first switched on.) The working temperature can be obtained from an iterative procedure. As a first approximation, we can assume the resistance to be unchanged, and then calculate the total power. For a black body, the total emitted power is proportional to the fourth power of temperature (Stefan’s law, \sigma T^4). Since the voltage has halved, the temperature has fallen by a factor of 1.414. The new temperature gives a resistivity value of around 30. The equation is easily solved using a computer (see the source code below). As can be seen, the final temperature is around 1830 K or 1600°C. The max intensity has shifted to a longer value. Even though the original value was also in the infrared region, this has much smaller amount in the visible. Once again, assuming the black body radiation formula, we have

{e_\lambda}_b = \frac{2\pi h c^2}{\lambda^5     \left[     \exp     \left(       hc/k_BT\lambda     \right) - 1     \right]}

#include
#include

int main(int argc, char** argv){
float q, prevq, prevr, r;

prevq = 2300;
prevr = 5.5 + (prevq – 300) * 0.025;
q = prevq/sqrt(2);

while (1) {
r = 5.5 + (q – 300) * 0.025;
q = q * sqrt(prevr/r);
printf(“Resistance = %5.0f\n Temperature = %5.0f\n”, r, q);
if (fabs (q – prevq) < 5 ) break; prevq = q; prevr = r; } return 0; } [/sourcecode] The amount of light available from the bulb should be adjusted with the brightness curve for the human eye, which is most sensitive to yellow. Thus, the power available in the red region of the spectrum is not as useful as the one in the yellow region. [caption id="attachment_197" align="aligncenter" width="320" caption="The total power output from bulbs at different voltages"]The total power output from bulbs at different voltages[/caption]

The above image shows that the power output (given by the area under the curve) is lower for the second case. However, one might expect that would have significant light since it is consuming 25 W of power.

Visible output from the two bulbs shown above

Visible output from the two bulbs shown above

The above image shows the output in the visible region—300 to 700 nm—and we can see that the line from the other bulb is not visible at all. This means that it will not have any visible light that you can see. What you effectively have is an infra red lamp!.

The illumination of a light source is provided in terms of candle power. Once again, this is not a good measure, unlike something like mass or length which can be measured to a very high degree of accuracy. Just like loudness of a sound, the optical property of brightness is a psychological phenomenon. And just like the psycho-acoustical phenomenon will lead you believe some sounds in certain frequencies (that commonly used in human speech) to be louder than others, the brightness is also mediated by the brain, so it varies from individual to individual.

Electric Lines of Force

September 14, 2008

Electric lines of force represent a powerful visualization tool, first popularized by Michael Faraday. Faraday, who claimed not to understand the advanced mathematical methods of the time, introduced it to a skeptical scientific community, and its true power would not be clear for decades.

Lines of Force for a Positive Charge

Lines of Force for a Positive Charge

The above picture shows the electric lines for a positive charge. According to convention, the lines move out from a positive charge. In the presence of a negative charge, they end at that charge, while in this case, they extend to infinity.

Field Lines for a Negative Charge

Field Lines for a Negative Charge

The field lines for a negative charge are as shown. The lines point inward.

Advanced Engineering Mathematics 5e, Kreyszig, 1983, pp. 46: Experiments show that the field lines for charges placed at (-1,0) and (1,0) are circles with equation x² + (y – c)² = 1 + c², while the equipotential lines are (x + C)² + y² = C² – 1. Note that the question is framed as an experimental result, so that it need not be the way the field of a dipole actually behaves.

Note that the potential of a dipole of moment p in polar co-ordinates is given by V = \frac{1}{4\pi\epsilon_0} \frac{p\cos\theta}{r^2}. This represents and ideal dipole of zero length and infinite charge, but finite value of dipole moment, i.e., it is a “pure” dipole.

A Field Line of  a Dipole

A Field Line of a Dipole

The above picture shows one field line of a dipole consisting of two equal and opposite charges separated by a distance of 10 units.

#include
#include

struct point {
float x;
float y;
};

typedef struct point Point;

Point plrcrsn(float, float);
Point ptadd(Point, Point);
float efangle(Point, Point[]);

int main(int argc, char** argv){
Point loc[2];
Point crn, nxt;
float r0;
float ang0;
int i = 200;

loc[0].x = 0;
loc[0].y = 0;
loc[1].x = 10;
loc[1].y = 0;

r0 = 0.1;
ang0 = M_PI / 2;

crn = ptadd(loc[0],plrcrsn(r0, ang0));

while (i– > 0) {
nxt = ptadd(crn, plrcrsn(r0, efangle(crn, loc)));
printf(“%3.2f %3.2f\n”, crn.x, crn.y);
crn = nxt;
}

return 0;
}

Point plrcrsn (float len, float ang) {
Point pt;
pt.x = len * cos(ang);
pt.y = len * sin(ang);
return pt;
}

Point ptadd (Point pt1, Point pt2) {
Point pt;

pt.x = pt1.x + pt2.x;
pt.y = pt1.y + pt2.y;
return pt;
}

float efangle(Point cur, Point loc[]) {
float efx, efy;
float x, y, z;
int i;
efx = efy = 0;
int sgn = -1;
for (i = 0; i < 2; i++) { sgn *= -1; x = cur.x - loc[i].x; y = cur.y - loc[i].y; z = pow( x * x + y * y, 1.5); efx += sgn * x/z; efy += sgn * y/z; } return atan (efy/efx); } [/sourcecode] It is easy to see that the above code can be easily modified to plot the equipotentials in the plane, as they are perpendicular to the field lines. Instead of moving along the field, we move perpendicular to the field at all times. However, as noted above, the code for the line of force is error prone, as it is a differential equation, attempting to estimate the derivative. This inherently has a larger error than finding the integral. The equipotentials are, of course, easily plotted by finding V at many points and joining the points of same potential. This can be made as accurate as you want. A better method would be to use a Runge-Kutta method (of any order) to calculate the lines of force.