Piecewise smooth functions are functions that are smooth (differentiable) on each piece of their domain, but may have a finite number of discontinuities at the boundaries between pieces. They are useful in mathematical modeling where different conditions apply in different regions, allowing for flexibility in capturing complex phenomena.