Jupiter's magnetosheath is a natural yet complex laboratory for analyzing compressible plasma turbulence. Recent observations by the Juno mission provide a promising opportunity for the first time to reckon the energy cascade rate in the magnetohydrodynamic scales in the vicinity of Jupiter's space. In the present work, a two-dimensional model is constructed for a whistler wave that is nonlinearly coupled with a wave magnetic field via ion density perturbation. The dynamics of whistler wave propagating in the direction of the magnetic field are derived within the limit of the two-fluid modeling approach. The magnetic field localization along with magnetic field spectra and spectral slope variations are estimated to realize the turbulence generation and energy cascade from large to small scales in the Jovian magnetosheath region. The simulated magnetic field spectrum in the wave number (in the unit of ion inertial length ρi) consists of turbulence in the inertial range with a spectral slope of -1.4 and a spectral knee at kρi = 1. Subsequently, the spectral slope increases to -2.6 and the spectrum becomes steeper. The simulated magnetic field spectrum in the wave number is further translated into the frequency domain using the whistler wave dispersion relation and by considering the Taylor frozen-in condition. The analytically estimated magnetic field spectrum slopes, i.e., -1.8 and -4.2 at low and high frequencies are further compared with recent Juno mission observations. The comparison further affirms the existence of Kolmogorov scaling, a spectral knee, and steepening in the spectrum at high frequencies. Furthermore, it is found that the two-fluid model can reasonably simulate the turbulence effects in Jovian magnetosheath in terms of magnetic field spectral distribution in wave number and frequency domains.