This study presents a high-fidelity direct numerical simulation (DNS) framework tailored for investigating turbulent flows through complex porous structures. It employs a compressible Navier–Stokes solver based on the spectral difference (SD) method, with immersed boundary conditions (IBCs) implemented via the Brinkman penalisation technique and integrated using a Strang splitting approach. A pressure gradient scaling (PGS) strategy is incorporated to improve computational efficiency. To provide realistic inflow conditions, synthetic turbulence is injected at the inlet using a random Fourier modes method. The methodology is validated in several stages. First, the IBC approach is tested against results from a body-fitted mesh, showing strong agreement in the mean velocity field. Next, the effectiveness of the PGS technique is demonstrated by comparing scaled and unscaled simulations, both of which yield consistent velocity fields and spectral content. Finally, the full DNS-SD framework is benchmarked against finite volume method results from the literature, successfully reproducing key turbulence characteristics, including two-point correlations. The validated solver is ultimately applied to simulate turbulent flow through a complex porous geometry. The results illustrate the robustness of the approach and highlight its potential for advancing the understanding of turbulence in porous materials.