Understanding the statistics of ocean geostrophic turbulence is of utmost importance in understanding its interactions with the global ocean circulation and the climate system as a whole. Here, a study of eddy-mixing entropy in a forced-dissipative barotropic ocean model is presented. Entropy is a concept of fundamental importance in statistical physics and information theory; motivated by equilibrium statistical mechanics theories of ideal geophysical fluids, we consider the effect of forcing and dissipation on eddy-mixing entropy, both analytically and numerically. By diagnosing the time evolution of eddy-mixing entropy it is shown that the entropy provides a descriptive tool for understanding three stages of the turbulence life cycle: growth of instability, formation of large scale structures and steady state fluctuations. Further, by determining the relationship between the time evolution of entropy and the maximum entropy principle, evidence is found for the action of this principle in a forced-dissipative flow. The maximum entropy potential vorticity statistics are calculated for the flow and are compared with numerical simulations. Deficiencies of the maximum entropy statistics are discussed in the context of the mean-field approximation for energy. This study highlights the importance entropy and statistical mechanics in the study of geostrophic turbulence.
A two-phase Smoothed Particle Hydrodynamics (SPH) model has been developed on the basis of GPUSPH, which is an open-source implementation of the weakly compressible SPH method on graphics processing units, to investigate oil dispersion under breaking waves. By assuming that the multiple phases are immiscible, the two-phase model solves the same set of governing equations for both phases. Density in each phase is preserved by renormalization, and the harmonic mean of viscosities is used in the transition zone. Interfacial surface tension effect between the oil and the water is considered by a numerical surface tension model. The model is first used to simulate a single oil drop rising through still water. The numerical model predicts realistic shape change of the oil drop during the rising process caused by the buoyancy force. Next it is applied to reproduce a laboratory experiment on oil dispersion under breaking waves conducted at Johns Hopkins University. Several high-speed cameras were used to record the interaction between breaking waves and the oil. Comparison with the laboratory measurements shows that GPUSPH is able to reproduce well the pre- & post-breaking wave in the laboratory, however, oil dispersion predicted by GPUSPH only match part of the laboratory observation. Several factors (e.g., 3D & chaotic nature of breaking waves, numerical setup) cause the discrepancy.
We formulate a general criterion for the exact preservation of the "lake at rest" solution in general mesh-based and meshless numerical schemes for the strong form of the shallow-water equations with bottom topography. The main idea is a careful mimetic design for the spatial derivative operators in the momentum flux equation that is paired with a compatible averaging rule for the water column height arising in the bottom topography source term. We prove consistency of the mimetic difference operators analytically and demonstrate the well-balanced property numerically using finite difference and RBF-FD schemes in the one- and two-dimensional cases.