A pseudospectral solver for the coupled YBJ-QG system in 2D
The coupled YBJ-QG system describes the evolution of a NIW field in the presence of a QG eddy field [1] as well as the back-reaction of the waves on the evolution of the QG field [2] [3]. The 2D version of this system, where we make a vertical plane wave assumption for the NIW field, was formulated in [4]. The model equations are:
where
We use dedalus to solve these equations on a periodic domain.
The test file contains code to run the two examples from [4].
rocha_anim.mp4
(cf. [4] Figure 1.)
[1] Young, W. R., Ben Jelloul, M. (1997). Propagation of near-inertial oscillations through a geostrophic flow Journal of Marine Research, 55(4), 735-766.
[2] Xie, J-H., and Vanneste, J. (2015). A generalised-Lagrangian-mean model of the interactions between near-inertial waves and mean flow Journal of Fluid Mechanics, 774, 143-169.
[3] Wagner, G. L., and W. R. Young. (2016). A three-component model for the coupled evolution of near-inertial waves, quasi-geostrophic flow and the near-inertial second harmonic Journal of Fluid Mechanics, 802, 806-837.
[4] Rocha, C., Wagner, G., & Young, W. (2018). Stimulated generation: Extraction of energy from balanced flow by near-inertial waves Journal of Fluid Mechanics, 847, 417-451.