Group delay generally varies with frequency in many systems that are of interest to us. For example, low-pass and high-pass filters all have nonlinear group delay near their cut-off frequencies. So group delay is far from being constant in many systems, although there are of course regions where it will be approximately constant too.
If we implement a Butterworth low-pass or high-pass filter in the digital domain, it will have the same group delay characteristics as the analog version of the filter. Of course, due to DSP latency, there will also be a time delay added to the signal because of the time it takes to perform the computations and then output the results through the DAC.
There is a well-known class of digital filters that have a constant group delay, owing to their linear-phase characteristics. These are known as Parks–McClellan linear-phase FIR digital filters, and procedures for their design were published back in 1973. The relevant paper,
A computer program for designing optimum FIR linear phase digital filters, describes their characteristics quite nicely.
That was an interesting plot, which I hadn't seen before. Thanks for sharing. It was interesting to see that the delta phase in the passband was within ±4 degrees. That's not a lot of variation. At 4kHz, a 4-degree phase difference corresponds to a spatial shift of 0.85 mm.