We present GTC MEGARA high-dispersion integral field spectroscopic observations of the nova remnant QU Vul, which provide a comprehensive 3D view of this nova shell. The tomographic analysis of the H𝛼 emission reveals a complex physical structure characterized by an inhomogeneous and clumpy distribution of the material within this shell. The overall structure can be described as a prolate ellipsoid with an axial ratio of 1.4±0.2, a major axis inclination with the line of sight of 12 • ± 6 • , and polar and equatorial expansion velocities ≈560 km s −1 and 400±60 km s −1 , respectively. The comparison of the expansion velocity on the plane of the sky with the angular expansion implies a distance of 1.43±0.23 kpc. The ionized mass is found to be ≈ 2 × 10 −4 M , noting that the information on the 3D distribution of material within the nova shell has allowed us to reduce the uncertainty on its filling factor. The nova shell is still in its free expansion phase, which can be expected as the ejecta mass is much larger than the swept-up circumstellar medium mass. The 3D distribution and radial velocity of material within the nova shell provide an interpretation of the so-called "castellated" line profiles observed in early optical spectra of nova shells, which can be attributed to knots and clumps moving radially along different directions.