Activities
In this small group activity, students integrate over non-uniform charge densities in cylindrical and spherical coordinates to calculate total charge.
Students work in small groups to use the superposition principle \[V(\vec{r}) = \frac{1}{4\pi\epsilon_0}\sum_i \frac{q_i}{\vert\vec{r}-\vec{r}_i\vert}\] to find the electrostatic potential \(V\) everywhere in space due to a pair of charges (either identical charges or a dipole). Different groups are assigned different arrangements of charges and different regions of space to consider: either on the axis of the charges or in the plane equidistant from the two charges, for either small or large values of the relevant geometric variable. Each group is asked to find a power series expansion for the electrostatic potential, valid in their group's assigned region of space. The whole class wrap-up discussion then compares and contrasts the results and discuss the symmetries of the two cases.
Students work in small groups to use the superposition principle \[V(\vec{r}) =\frac{1}{4\pi\epsilon_0}\int\frac{\rho(\vec{r}^{\,\prime})}{\vert \vec{r}-\vec{r}^{\,\prime}\vert} \, d\tau^{\prime}\] to find an integral expression for the electrostatic potential, \(V(\vec{r})\), everywhere in space, due to a ring of charge.
In an optional extension, students find a series expansion for \(V(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.
Students are prompted to consider the scalar superposition of the electric potential due to multiple point charges. First a single point charge is discussed, then four positive charges, then an electric quadrupole. Students draw the equipotential curves in the plane of the charges, while also considering the 3D nature of equipotentials.
Students use known algebraic expressions for length elements \(d\ell\) to determine all simple scalar area \(dA\) and volume elements \(d\tau\) in cylindrical and spherical coordinates.
This activity is identical to Vector Surface and Volume Elements except uses a scalar approach to find surface, and volume elements.
A pretzel is to be dipped in chocolate. The pretzel is in the shape of a quarter circle, consisting of a straight segment from the origin to the point (2,0), a circular arc from there to (0,2), followed by a straight segment back to the origin; all distances are in centimeters. The (linear) density of chocolate on the pretzel is given by \(\lambda = 3(x^ 2 + y^2 )\) in grams per centimeter. Find the total amount of chocolate on the pretzel.Main ideas
- Calculating (scalar) line integrals.
- Use what you know!
Prerequisites
- Familiarity with \(d\boldsymbol{\vec{r}}\).
- Familiarity with “Use what you know” strategy.
Warmup
It is not necessary to explicitly introduce scalar line integrals, before this lab; figuring out that the (scalar) line element must be \(|d\boldsymbol{\vec{r}}|\) can be made part of the activity (if time permits).
Props
- whiteboards and pens
- “linear” chocolate covered candy (e.g. Pocky)
Wrapup
Emphasize that students must express each integrand in terms of a single variable prior to integration.
Emphasize that each integral must be positive!
Discuss several different ways of doing this problem (see below).
Details
In the Classroom
- Make sure the shape of the pretzel is clear! It might be worth drawing it on the board.
- Some students will work geometrically, determining \(ds\) on each piece by inspection. This is fine, but encourage such students to try using \(d\vec{r}\) afterwards.
- Polar coordinates are natural for all three parts of this problem, not just the circular arc.
- Many students will think that the integral “down” the \(y\)-axis should be negative. They will argue that \(ds=dy\), but the limits are from \(2\) to \(0\). The resolution is that \(ds = |dy\,\boldsymbol{\hat x}|=|dy|=-dy\) when integrating in this direction.
- Unlike work or circulation, the amount of chocolate does not depend on which way one integrates, so there is in fact no need to integrate “down” the \(y\)-axis at all.
- Some students may argue that \(d\boldsymbol{\vec{r}}=\boldsymbol{\hat T}\,ds\Longrightarrow ds=d\boldsymbol{\vec{r}}\cdot\boldsymbol{\hat T}\), and use this to get the signs right. This is fine if it comes up, but the unit tangent vector \(\boldsymbol{\hat T}\) is not a fundamental part of our approach.
- There is of course a symmetry argument which says that the two “legs” along the axes must have the same amount of chocolate --- although some students will put a minus sign into this argument!
Subsidiary ideas
- \(ds=|d\boldsymbol{\vec{r}}|\)
Consider the electric field \begin{equation} \vec E(r,\theta,\phi) = \begin{cases} 0&\textrm{for } r<a\\ \frac{1}{4\pi\epsilon_0} \,\frac{Q}{b^3-a^3}\, \left( r-\frac{a^3}{r^2}\right)\, \hat r & \textrm{for } a<r<b\\ 0 & \textrm{for } r>b \\ \end{cases} \end{equation}
- (4pts) Use step and/or delta functions to write this electric field as a single expression valid everywhere in space.
- (4pts) Find a formula for the charge density that creates this electric field.
- (2pts) Interpret your formula for the charge density, i.e. explain briefly in words where the charge is.
Students integrate numerically to find the electric field due to a cone of surface charge, and then visualize the result. This integral can be done in either spherical or cylindrical coordinates, giving students a chance to reason about which coordinate system would be more convenient.
Students work in small groups to use Coulomb's Law \[\vec{E}(\vec{r}) =\frac{1}{4\pi\epsilon_0}\int\frac{\rho(\vec{r}^{\,\prime})\left(\vec{r}-\vec{r}^{\,\prime}\right)}{\vert \vec{r}-\vec{r}^{\,\prime}\vert^3} \, d\tau^{\prime}\] to find an integral expression for the electric field, \(\vec{E}(\vec{r})\), everywhere in space, due to a ring of charge.
In an optional extension, students find a series expansion for \(\vec{E}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.
Students work in small groups to use the Biot-Savart law \[\vec{B}(\vec{r}) =\frac{\mu_0}{4\pi}\int\frac{\vec{J}(\vec{r}^{\,\prime})\times \left(\vec{r}-\vec{r}^{\,\prime}\right)}{\vert \vec{r}-\vec{r}^{\,\prime}\vert^3} \, d\tau^{\prime}\] to find an integral expression for the magnetic field, \(\vec{B}(\vec{r})\), due to a spinning ring of charge.
In an optional extension, students find a series expansion for \(\vec{B}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.
Consider the magnetic field \[ \vec{B}(s,\phi,z)= \begin{cases} 0&0\le s<a\\ \alpha \frac{1}{s}(s^4-a^4)\, \hat{\phi}&a<s<b\\ 0&s>b \end{cases} \]
- (2pts) Use step and/or delta functions to write this magnetic field as a single expression valid everywhere in space.
- (4pts) Find a formula for the current density that creates this magnetic field.
- (2pts) Interpret your formula for the current density, i.e. explain briefly in words where the current is.
Consider a column of atoms each of mass \(M\) at temperature \(T\) in a uniform gravitational field \(g\). Find the thermal average potential energy per atom. The thermal average kinetic energy is independent of height. Find the total heat capacity per atom. The total heat capacity is the sum of contributions from the kinetic energy and from the potential energy. Take the zero of the gravitational energy at the bottom \(h=0\) of the column. Integrate from \(h=0\) to \(h=\infty\). You may assume the gas is ideal.
(4pts) Consider a thin charged rod of length \(L\) standing along the \(z\)-axis with the bottom end on the \(xy\)-plane. The charge density \(\lambda\) is constant. Find the electric field at the point \((0,0,2L)\).
(4pts) Find the electric field around an infinite, uniformly charged, straight wire, starting from the following expression for the electrostatic potential: \begin{equation*} V(\vec r)=\frac{2\lambda}{4\pi\epsilon_0}\, \ln\left( \frac{ s_0}{s} \right) \end{equation*}
Consider the finite line with a uniform charge density from class.
- Write an integral expression for the electric field at any point in space due to the finite line. In addition to your usual physics sense-making, you must include a clearly labeled figure and discuss what happens to the direction of the unit vectors as you integrate.Consider the finite line with a uniform charge density from class.
- Perform the integral to find the \(z\)-component of the electric field. In addition to your usual physics sense-making, you must compare your result to the gradient of the electric potential we found in class. (If you want to challenge yourself, do the \(s\)-component as well!)
The gravitational field due to a spherical shell of matter (or equivalently, the electric field due to a spherical shell of charge) is given by: \begin{equation} \vec g = \begin{cases} 0&\textrm{for } r<a\\ -G \,\frac{M}{b^3-a^3}\, \left( r-\frac{a^3}{r^2}\right)\, \hat r & \textrm{for } a<r<b\\ -G\,\frac{M}{r^2}\, \hat r & \textrm{for } r>b \\ \end{cases} \end{equation}
This problem explores the consequences of the divergence theorem for this shell.
- Using the given description of the gravitational field, find the divergence of the gravitational field everywhere in space. You will need to divide this question up into three parts: \(r<a\), \(a<r<b\), and \(r>b\).
- Briefly discuss the physical meaning of the divergence in this particular example.
- For this gravitational field, verify the divergence theorem on a sphere, concentric with the shell, with radius \(Q\), where \(a<Q<b\). ("Verify" the divergence theorem means calculate the integrals from both sides of the divergence theorem and show that they give the same answer.)
- Briefly discuss how this example would change if you were discussing the electric field of a uniformly charged spherical shell.
Problem
Nuclei of a particular isotope species contained in a crystal have spin \(I=1\), and thus, \(m = \{+1,0,-1\}\). The interaction between the nuclear quadrupole moment and the gradient of the crystalline electric field produces a situation where the nucleus has the same energy, \(E=\varepsilon\), in the state \(m=+1\) and the state \(m=-1\), compared with an energy \(E=0\) in the state \(m=0\), i.e. each nucleus can be in one of 3 states, two of which have energy \(E=\varepsilon\) and one has energy \(E=0\).
Find the Helmholtz free energy \(F = U-TS\) for a crystal containing \(N\) nuclei which do not interact with each other.
Find an expression for the entropy as a function of temperature for this system. (Hint: use results of part a.)
- Indicate what your results predict for the entropy at the extremes of very high temperature and very low temperature.
Students write python programs to compute the potential due to a square of surface charge, and then to visualize the result. This activity can be used to introduce students to the process of integrating numerically.
- Students need to understand that the surface represents the electric potential in the center of a parallel plate capacitor. Try doing the activity Electric Potential of Two Charged Plates before this activity.
- Students should know that
- objects with like charge repel and opposite charge attract,
- object tend to move toward lower energy configurations
- The potential energy of a charged particle is related to its charge: \(U=qV\)
- The force on a charged particle is related to its charge: \(\vec{F}=q\vec{E}\)
Students work in small groups to use the superposition principle \[\vec{A}(\vec{r}) =\frac{\mu_0}{4\pi}\int\frac{\vec{J}(\vec{r}^{\,\prime})}{\vert \vec{r}-\vec{r}^{\,\prime}\vert}\, d\tau^{\prime}\] to find an integral expression for the magnetic vector potential, \(\vec{A}(\vec{r})\), due to a spinning ring of charge.
In an optional extension, students find a series expansion for \(\vec{A}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.
The electrostatic potential due to a point charge at the origin is given by: \begin{equation*} V=\frac{1}{4\pi\epsilon_0} \frac{q}{r} \end{equation*}
- (2pts) Find the electric field due to a point charge at the origin as a gradient in rectangular coordinates.
- (2pts) Find the electric field due to a point charge at the origin as a gradient in spherical coordinates.
- (2pts) Find the electric field due to a point charge at the origin as a gradient in cylindrical coordinates.
Students, pretending they are point charges, move around the room so as to make an imaginary magnetic field meter register a constant magnetic field, introducing the concept of steady current. Students act out linear \(\vec{I}\), surface \(\vec{K}\), and volume \(\vec{J}\) current densities. The instructor demonstrates what it means to measure these quantities by counting how many students pass through a gate.
Students will estimate the work done by a given electric field. They will connect the work done to the height of a plastic surface graph of the electric potential.
Students use a plastic surface representing the potential due to a charged sphere to explore the electrostatic potential, equipotential lines, and the relationship between potential and electric field.
These lecture notes from the ninth week of https://paradigms.oregonstate.edu/courses/ph441 cover phase transformations, the Clausius-Clapeyron relation, mean field theory and more. They include a number of small group activities.
Students consider the relation (1) between the angular momentum and magnetic moment for a current loop and (2) the force on a magnetic moment in an inhomogeneous magnetic field. Students make a (classical) prediction of the outcome of a Stern-Gerlach experiment.
In this small group activity, students determine various aspects of local points on an elliptic hill which is a function of two variables. The gradient is emphasized as a local quantity which points in the direction of greatest change at a point in the scalar field.
Begin by prompting the students to brainstorm different ways to represent a three dimensional scalar field on a 2-D surface (like their paper or a whiteboard). The students use a pre-made Sage code or a Mathematica worksheet to visualize the electrostatic potential of several distributions of charges. The computer algebra systems demonstrate several different ways of plotting the potential.
In this small group activity, students determine various aspects of local points on an elliptic hill which is a function of two variables. The gradient is emphasized as a local quantity which points in the direction of greatest change at a point in the scalar field.
Students use prepared Sage code to predict the gradient from contour graphs of 2D scalar fields.
Students set up and compute a scalar surface integral.
Students are asked to review:in preparation for an in-class quiz.
- Addition of matrices
- Multiplication of a matrix by a scalar
- Matrix multiplication
- Finding the determinant of a matrix
Students hold rulers and meter sticks to represent a vector field. The instructor holds a hula hoop to represent a small area element. Students are asked to describe the flux of the vector field through the area element.
Students explore path integrals using a vector field map and thinking about integration as chop-multiply-add.
Student discuss how many paths can be found on a map of the vector fields \(\vec{F}\) for which the integral \(\int \vec{F}\cdot d\vec{r}\) is positive, negative, or zero. \(\vec{F}\) is conservative. They do a similar activity for the vector field \(\vec{G}\) which is not conservative.
Students explore the effects of putting a point charge at various places inside, outside, and on the surface of a cubical Gaussian surface. The Mathematica worksheet or Sage activity shows the electric field due to the charge, then plots the the flux integrand on the top surface of the box, calculates the flux through the top of the box, and the value of the flux through the whole cube.
Students calculate the flux from the vector field \(\vec{F} = C\, z\, \hat{z}\) through a right cone of height \(H\) and radius \(R\) .
Students predict from graphs of simple 2-d vector fields whether the divergence is positive, negative, or zero in various regions of the domain using the geometric definition of the divergence of a vector field at a point as flux per unit volume (here: area) through an infinitesimal box surrounding that point. Optionally, students can use computer algebra to verify their predictions.
Students predict from graphs of simple 2-d vector fields whether the curl is positive, negative, or zero in various regions of the domain using the definition of the curl of a vector field at a point as the maximum circulation per unit area through an infinitesimal box surrounding that point. Optionally, students can use computer algebra to verify their predictions.