An axisymmetric, nonhydrostatic version of the ZetaC model, with the GFDL radiation package on a 10 degrees domain, has been used to produce a large number of simulations covering a vast range of environmental conditions (surface temperature, tropopause temperature, latitude, surface heat flux and other parameters were varied). In the moist simulations we use the Lin microphysics. The horizontal resolution of the model is 1 km or 500m thus convection is fully resolved. Simulations are run for about 25 days, producing a steady state.
Our analysis of the numerical simulations suggests that the entropy distribution in the boundary layer cannot be determined solely from the surface fluxes. In particular, turbulent mixing across entropy and angular momentum surfaces is associated with a strengthening of the entropy gradient near the radius of maximum wind. In our simulations, this thermodynamic effect did result in a 5 to 10 percent increase in the maximum wind speed. Numerical simulations show that the maximum tangential winds are stronger than theory predicts, which means that the MPI theory is not really an upper bound for intensity. The superintensity can be explained by the presence of ageostrophic winds within the eyewall inflow region.