Modeling global solar convection requires capturing both the highly stratified turbulent envelope and its coupling with the stable radiative interior. Traditional global simulations rely on anelastic approximations, which filter out acoustic waves and fail near the surface where flows become transonic. In this work, we use the fully compressible finite-volume AMR code Dyablo Whole Sun to perform global solar simulations. We overcome spherical coordinate singularities using non-cartesian mapped grids and ensure hydrostatic stability across large density stratification (N_=7) using the - well-balanced scheme. We successfully model the coupling of the convective and radiative zones in 2D and 3D in hydrodynamic setups. This fully compressible framework allows us to report the first joint observation of acoustic and gravity modes in a global solar simulation, showing excellent agreement with theoretical spectra obtained via the GYRE stellar oscillation code.
Grégoire Doebele (Sat,) studied this question.