Para algunas observaciones sobre cómo es la historia general, es posible que desee consultar el artículo del Princeton Companion sobre ecuaciones diferenciales parciales. Personalmente, no tengo mucho que decir al respecto, así que solo mencionaré dos ejemplos.
La ecuación de Korteweg-de Vries es de interés para las personas que trabajan en sistemas integrables (hay otros ejemplos en ese enlace). Entiendo que tiene algunas conexiones interesantes con la geometría algebraica.
Otro ejemplo (realmente una familia de ejemplos) es el siguiente. Sea $G$ un grupo de Lie que actúa suavemente sobre una variedad suave $M$. Esto induce un mapa del álgebra de Lie $\mathfrak{g}$ al álgebra de Lie de campos vectoriales suaves en $M$, lo que induce una acción del álgebra envolvente universal $U(\mathfrak{g})$ de $\mathfrak{g}$ en $C^{\infty}(M)$ mediante operadores diferenciales. Elija $D \in U(\mathfrak{g})$. Entonces puedes escribir la ecuación diferencial $$Df = 0$$
donde $f \in C^{\infty}(M)$. Esta ecuación diferencial debería tener un orden igual al grado de $D$ (con respecto a la filtración usual en $U(\mathfrak{g})$). El espacio de soluciones a esta ecuación diferencial es actuado por el subgrupo de $G$ que conmuta con $D$. En particular, si $D$ yace en el centro, el espacio de soluciones es actuado por todo $G$.
Para un ejemplo muy simple, deja que $\mathbb{R}$ actúe sobre sí mismo por traslación. Luego $\mathfrak{g}$ está generado por la diferenciación ordinaria $D$, el álgebra envolvente universal $U(\mathfrak{g})$ son polinomios en $D$, y las ecuaciones diferenciales correspondientes que puedes escribir son todas EDO lineales homogéneas con coeficientes constantes, todas las cuales admiten una acción de $G$ (nuevamente por traslación). Escribir esta acción explícitamente da generalizaciones de las fórmulas de adición de ángulos (que obtienes en el caso especial $D^2 + 1$).
Para un ejemplo no abeliano, deja que $\text{SO}(3)$ actúe en $\mathbb{R}^3$ por rotación. Luego la imagen de $\mathfrak{g} = \mathfrak{so}(3)$ en campos vectoriales en $\mathbb{R}^3$ está generada por
$$L_x = y \frac{\partial}{\partial z} - z \frac{\partial}{\partial y}$$ $$L_y = z \frac{\partial}{\partial x} - x \frac{\partial}{\partial z}$$ $$L_z = x \frac{\partial}{\partial y} - y \frac{\partial}{\partial x}.$$
Estos campos vectoriales generan rotación alrededor de los ejes $x, y, z$ respectivamente. El operador de Casimir $$L^2 = L_x^2 + L_y^2 + L_z^2$$
resulta generar el centro de $U(\mathfrak{so}(3))$, por lo que el espacio de soluciones de la ecuación diferencial $p(L^2) f = 0$ para cualquier polinomio $p$ es actuado por todo $\text{SO}(3)$ y tiene orden $2 \deg p$.