I am going to be a TA for a course on partial differential equations in the coming semester. I recall in my undergraduate years taking a course on PDE, following Evans' book. One of the key components of the course was the theory of distributions. Unfortunately, the subject was treated very abstractly, and few examples were given, let alone treated with great detail.
I'm sure that I'm not alone in my passion for a bank of examples. As a consequence, I'd be interested to know of your favourite examples that can appear in a course of this nature.
I should mention that I'm not proficient in PDEs, so I can only give my perspective as a student who has taken an introductory (non rigorous) course in PDEs (but having much more real/functional analysis background).
The very first example of distributions is of course given by locally Lebesgue integrable functions: for any open $U\subset\Bbb{R}^n$, we have the canonical injection $L^1_{\text{loc}}(U)\to \mathcal{D}'(U)$ given by integration against test functions. These of course include all continuous functions, and much more.
Next, we have the dirac point masses, and also analogues for submanifolds. Given a point $p\in\Bbb{R}^n$, we can define $\delta_p$ by setting $\langle\delta_p,\phi\rangle:=\phi(p)$ (for $\phi\in C^{\infty}(\Bbb{R}^n)$). This is of course a very basic distribution to introduce because it appears very naturally when talking about fundamental solutions of PDEs. Along these lines, if $M\subset\Bbb{R}^n$ is any embedded submanifold, we can define $\delta_M$ by setting $\langle \delta_M,\phi \rangle:=\int_M\phi(x)\,dV(x)$, where $dV$ is the volume element of $M$ (induced by the pullback of Riemannian metric of $\Bbb{R}^n$ to $M$). A very simple example is of course to let $M=S_r(p)$ be a sphere in $\Bbb{R}^n$ of radius $r$, center $p$. This example might be relevant if you're discussing the fundamental solution of the wave equation in $\Bbb{R}^3$.
Elaborating slightly on (1), some important $L^1_{\text{loc}}$ functions are of course the fundamental solutions for the Laplace equation: \begin{align} L(x)&= \begin{cases} \frac{-1}{(n-2)A_n}\frac{1}{|x|^{n-2}}&\text{if $n\geq 3$}\\\\ \frac{1}{2\pi}\log|x|&\text{if $n=2$} \end{cases} \end{align} where $A_n$ is the surface area of the unit sphere in $\Bbb{R}^n$. Also, the fundamental solution of the heat equation: \begin{align} H(x,t)&=\begin{cases} \frac{1}{(4\pi t)^{n/2}}e^{-\frac{|x|^2}{4t}}&\quad \text{if $t>0$}\\ 0&\quad \text{if $t\leq 0$} \end{cases} \end{align} Clearly, the Laplace and Heat equations are very important, so it makes sense to talk about these functions (which are smooth at most places), and it of course be good to show that indeed their distributional derivatives satisfy $\Delta L=\delta$ and $(\partial_t-\Delta)H=\delta$ (if this isn't already done in lectures). I found it very instructive when my TA went through the calculation for showing $H$ really is the fundamental solution. These calculations made it sink in how derivatives of distributions work.
Of course one of the most important operations on (tempered) distributions is the Fourier transform, which is very relevant when dealing with the basic three PDEs (Laplace, Heat, and Wave), so you may want to calculate some Fourier transforms (especially of functions which are not in $L^1(\Bbb{R}^n)$) rigorously.
As a side note, the heat/Weierstrass kernel $H$ above can also be used to give a proof of Weierstrass' approximation theorem (not sure how relevant this may be for your course, but I find this interesting).