Here's a question I'm not even remotely qualified to ask: could it be that supermassive black holes cause a contraction of local spacetime in a way that gives the impression that local space is static while distant space expands? I think this mechanism would be largely indistinguishable from classic hubble expansion, at least on local scales.
That's a good question and is not exactly easy to answer. The Einstein field equations of General Relativity are highly non-linear which makes it difficult to talk about superpositions of solutions to the equations and to compare local physics (at the level of black holes) to cosmological effects (e.g. the expansion of the universe) and transfer results and insights between them. Let me explain. When we talk about co…
> [link to wikipedia's de Sitter-Schwarzschild page]
Let's call it Schwarzschild-de Sitter (SdS, for short, and for ease of literature-searching).
More generally there is the McVittie family of metrics of a massive object in a dynamical spacetime. SdS is the limiting case of McVittie where the spacetime is stationary and the central mass is compact, spherically symmetric, 0-angular-momentum, and uncharged. (One could alternatively say that McVittie is a generalized time-dependent SdS.)
> black hole solutions assume a perfect vacuum
Not quite, but this is mainly a specialist quibble, since the most widely known theoretical black holes are vacuum or electrovac spacetimes. See for example the Kerr-Vaidya black hole solutions, which have either a incoming radiation ("null dust") field or an outgoing one (or both) falling onto resp. shining out of ("roughly Hawking") a spinning black hole. There are many other nonvacuum solutions with a compact central mass (which can look more or less black-hole-like), both exact and non-exact.
The important feature of these is asymptotic behaviour, as you touch on in your second-last paragraph, since if the influence of the central mass fades with distance, and the sources are kept distant from each other, that lets us ignore (but see below) the difficulties in combining two or more exact solutions of the Einstein Field Equations into a new exact solution.
Linearized gravity is usually applicable and sufficient, and if not one can obtain corrections using post-Newtonian theory. We don't really need numrel unless mass-ratios are small and compactness is extreme. See the handy diagram at https://en.wikipedia.org/wiki/Post-Newtonian_expansion#/medi...>.
See also Ellis 2010 (Chapter 2, section 3 on inexact solutions, notably his complaint at the bottom of p. 34 to the top of p. 35) https://doi.org/10.1017/CBO9780511622724.002>, which is handily also at s c y h o b. Re his complaint see also Visser 2014 on horizons: https://arxiv.org/abs/1407.7295>.