On the generation and evolution of numerically simulated large-amplitude internal gravity wave packets

Numerical simulations of internal gravity wave (IGW) dynamics typically rely on wave velocity and density fields which are either generated through forcing terms in the governing equations or are explicitly introduced as initial conditions. Both approaches are based on the associated solution to the inviscid linear internal wave equations and, thus, assume weak-amplitude, space-filling waves. Using spectral multidomain-based numerical simulations of the two-dimensional Navier-Stokes equations and focusing on the forcing-driven approach, this study examines the generation and subsequent evolution of large-amplitude IGW packets which are strongly localized in the vertical in a linearly stratified fluid. When the vertical envelope of the forcing terms varies relatively rapid when compared to the vertical wavelength, the associated large vertical gradients in the Reynolds stress field drive a nonpropagating negative horizontal mean flow component in the source region. The highly nonlinear interaction of this mean current with the propagating IGW packet leads to amplification of the wave, a significant distortion of its rear flank, and a substantial decay of its amplitude. Scaling arguments show that the mean flow is enhanced with a stronger degree of localization of the forcing, larger degree of hydrostaticity, and increasing wave packet steepness. Horizontal localization results in a pronounced reduction in mean flow strength mainly on account of the reduced vertical gradient of the wave Reynolds stress. Finally, two techniques are proposed toward the efficient containment of the mean flow at minimal computational cost. The findings of this study are of particular value in overcoming challenges in the design of robust computational process studies of IGW packet (or continuously forced wave train) interactions with a sloping boundary, critical layer, or caustic, where large wave amplitudes are required for any instabilities to develop. In addition, the detailed description of the dynamics of the wave generation region may offer a first probe into the underlying physics of disruptive near-source nonlinearities observed in laboratory experiments of persistently forced IGW beams. Finally, the question arises as to whether a highly vertically compact IGW packet, which has propagated far from its source but still maintains its original structure and amplitude, can indeed occur in nature.