Rapid intensification of tropical cyclones is a challenge for forecasters. In 2017, Hurricane Maria intensified to a Category 5 storm within 24 hours and devastated Puerto Rico. The official forecast and all computer models were unable to predict it. Hurricane Dorian had been predicted as a tropical storm; unexpectedly, it intensified into a Category 5 storm and destroyed the Bahamas. Soloviev et al. (2017) suggested that under the assumption of constant enthalpy exchange coefficient, rapid cyclone intensification and decay can be related to the drag coefficient dependence on wind speed including an “aerodynamic drag well” around 60 m/s. This concept is in general terms consistent with Emanuel’s (1988) theory of maximum potential intensity of a tropical cyclone and its extension by Lee et al. (2019). The influence of sea spray is still a significant uncertainty. In order to study the effect of spray on dynamics of tropical cyclones, we have implemented a Volume of Fluid to Discrete-Phase Model (VOF to DPM). This model re-meshes the areas with increased gradients or curvature, which are suspicious for the interface instability. The generated water particles that satisfy the condition of asphericity are converted into Lagrangian particles. The size distribution of spray measured in air-sea interaction facilities is used for the model verification. Due to dynamic remeshing, VOF to DPM resolves spray particle radius from ten micrometers to a few millimeters, which correspond to spume. Results of the numerical simulation show a dramatic increase of spume generation under major tropical cyclones. Though sub-micrometer and micrometer scale spray particles are not resolved in this simulation, they are likely less significant in the momentum exchange at the air-sea interface than spume. These results are expected to contribute to the parameterization and proper treatment of spray in forecasting models, including cases of rapid intensification and rapid decline of tropical cyclones.