Isotropic Coordinates for Generalized Schwarzschild-like Black Holes with Anisotropic Media

Real black holes do not live in vacuum. Dark-energy-like components, topological-defect networks, anisotropic fluids, and extended field configurations can imprint measurable environmental effects on strong-gravity observables. In this regime, when Einstein's equations cannot be solved analytically, numerical relativity is the standard tool. Within numerical relativity, the 3+1 decomposition provides a practical evolution framework, making clean initial data and horizon-regular coordinates essential. We consider a broad class of static, spherically symmetric generalized Schwarzschild-like solutions with multiple non-interacting anisotropic fluid sources and derive the coordinate transformation from Schwarzschild-like to isotropic coordinates with conformally flat spatial slices. The isotropic form removes spatial-sector coordinate pathologies at the horizon, clarifies geometric quantities, and enables the construction of well-posed initial data on t=const hypersurfaces, suitable for the Hamiltonian and conformal formulations of numerical relativity and for perturbation theory. The isotropic backgrounds we develop make it straightforward to separate environmental effects from intrinsic strong-gravity signals and meet the growing interest in non-vacuum black-hole phenomenology in scattering, lensing, and waveform modeling.