In the supersymmetric quantum mechanics formalism, the shape invariance
condition provides a sufficient constraint to make a quantum mechanical problem
solvable; i.e., we can determine its eigenvalues and eigenfunctions
algebraically. Since shape invariance relates superpotentials and their
derivatives at two different values of the parameter a, it is a non-local
condition in the coordinate-parameter (x,a) space. We transform the shape
invariance condition for additive shape invariant superpotentials into two
local partial differential equations. One of these equations is equivalent to
the one-dimensional Euler equation expressing momentum conservation for
inviscid fluid flow. The second equation provides the constraint that helps us
determine unique solutions. We solve these equations to generate the set of all
known ℏ-independent shape invariant superpotentials and show that there
are no others. We then develop an algorithm for generating additive shape
invariant superpotentials including those that depend on ℏ explicitly,
and derive a new ℏ-dependent superpotential by expanding a Scarf
superpotential.Comment: 1 figure, 4 tables, 18 page