In this activity, students apply the Stefan-Boltzmann equation and the principle of energy balance in steady state to find the steady state temperature of a black object in near-Earth orbit.
This small group activity using surfaces combines practice with the multivariable chain rule while emphasizing numerical representations of derivatives.
Students work in small groups to measure partial derivatives in both rectangular and polar coordinates, then verify their results using the chain rule.
The whole class wrap-up discussion emphasizes the relationship between a directional derivative in the \(r\)-direction and derivatives in \(x\)- and \(y\)-directions using the chain rule.
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.
Consider a system consisting of four point charges arranged on the corners of a square in 3D Cartesian space of coordinates \((x,y,z)\).
Write a python function that returns the potential at any point in space caused by four equal point charges forming a square. Make the sides of the square parallel to the \(x\) and \(y\) axes and on the \(z=0\) plane.
To do this you will need the expression for a the potential due to a single point charge \(V= \frac{k_Cq}{r}\) where \(r\) is the distance from the point charge. You will also need to use the fact that the total potential is the sum of the potentials due to each individual point charge.
It is important that we ask students first to create a function for the potential, and only then try to visualize the potential. This allows students to reason about the computation for a single point in space (defined in their choice of coordinate systems).
Once you have written the above function, use it to plot the electrostatic potential versus position along the three cartesian axes.
Since the students have already written a function for their potential, they can create a plot by creating an array for \(x\) (or \(y\), or \(z\)), and then passing that array to their function, along with scalars for the other two coordinates. Many students will discover this simply by modifying an example script they find on the web, replacing \(\sin(x)\) or similar with their function. It is well worth showing this easier approach to students who attempt who attempt to write a loop in order to compute the potential at each point in space.
We ask students to explicitly plot the potential along axes because students seldom spontaneously think to create a 1D plot such as this.
Label your axes.
Work out the first non-zero term in a power series approximation for the potential at large \(x\), small \(x\), etc. Plot these approximations along with your computed potential, and verify that they agree in the range that you expect. Useful 1\(^{st}\) order Taylor expansions are:
\begin{eqnarray}
\sqrt{1+\epsilon} &\sim& 1+\frac{\epsilon}{2} \\
\frac{1}{1+\epsilon} &=& 1-\epsilon
\end{eqnarray}
where \(\epsilon\) is a small quantity.
This may need to be omitted on the first Tuesday of class, since students probably will not yet have seen power series approximations. It may work in this case to at least talk about what is expected at large distance, since "it looks like a point charge" is reasoning students do make.
Students struggle with the \(x\) approximations (assuming the square is in the xy plane). Each pair will probably need to have a little lecture on grouping terms according to the power of \(x\), and keeping only those terms for which they have every instance.
Extra fun
Create one or more different visualizations of the electrostatic potential. For example a 2D representation in the \(z=0\) plane.
More extra fun
Create a plot of the potential along a straight line that is not one of the axes. Hint: start from a line on the \(z=0\) plane, then try a random straight line. You can use your browser for help.
Even more extra fun
Move the charges around (e.g., off the \(z=0\) plane) and see what happens to your graphs
Dipole fun
Repeat the above (especially the limiting cases!) for four point charges in which half are positive and half negative, with the positive charges neighbors.
Common visualizations for 2D slices of space include contour plots, color plots, and "3D plots". Another option (less easy) would be to visualize an equipotential surface in 3 dimensions. It is worth reminding students to consider other planes than those at \(x=0\), \(y=0\), and \(z=0\).
Quadrupole fun
Repeat the above (especially the limiting cases!) for four point charges in which half are positive and half negative, with the positive charges diagonal from one another. It will help in this case to place the charges on the axes (rotating the square by 45 degrees), since otherwise the potential on each axis will be zero.
Students solve numerically for the potential due to a spherical shell of charge. Although this potential is straightforward to compute using Gauss's Law, it serves as a nice example for numerically integrating in spherical coordinates because the correct answer is easy to recognize.
Students implement a finite-difference approximation for the kinetic energy operator as a matrix, and then use numpy to solve for eigenvalues and eigenstates, which they visualize.
Students compute probabilities and averages given a probability density in one dimension. This activity serves as a soft introduction to the particle in a box, introducing all the concepts that are needed.
Students draw the 3D graphs of equations using three variables. They make choices for drawing a stack of curves in parallel planes and a curve in a perpendicular plane (e.g. substituting in values for \(x\), \(y\), or \(z\). )
Students find matrix elements of the position operator \(\hat x\) in a sinusoidal basis. This allows them to express this operator as a matrix, which they can then numerically diagonalize and visualize the eigenfunctions.
Students compute inner products to expand a wave function in a sinusoidal basis set. This activity introduces the inner product for wave functions, and the idea of approximating a wave function using a finite set of basis functions.
This small group activity is designed to help students visual the process of chopping, adding, and multiplying in single integrals.
Students work in small groups to determine the volume of a cylinder in as many ways as possible.
The whole class wrap-up discussion emphasizes the equivalence of different ways of chopping the cylinder.
This small group activity has students reasoning about how the Planck distribution shifts when the temperature is doubled. This leads to a qualitative argument for the Stefan-Boltzmann law.