Global existence for the two-dimensional Kuramoto–Sivashinsky equation with a shear flow

We consider the Kuramoto–Sivashinsky equation (KSE) on the two-dimensional torus in the presence of advection by a given background shear flow. Under the assumption that the shear has a finite number of critical points and there are linearly growing modes only in the direction of the shear, we prove global existence of solutions with data in $$L^2$$ L 2 , using a bootstrap argument. The initial data can be taken arbitrarily large.

Suppressed Thermocline Mixing in the Center of Anticyclonic Eddy in the North South China Sea

Direct microstructure observations and fine-scale measurements of an anticyclonic mesoscale eddy were conducted in the northern South China Sea in July 2020. An important finding was that suppressed turbulent mixing in the thermocline existed at the center of the eddy, with an averaged diapycnal diffusivity at least threefold smaller than the peripheral diffusivity. Despite the strong background shear and significant wave–mean flow interactions, the results indicated that the lack of internal wave energy in the corresponding neap tide period during measurement of the eddy’s center was the main reason for the suppressed turbulent mixing in the thermocline. The applicability of the fine-scale parameterization method in the presence of significant wave–mean flow interactions in a mesoscale eddy was evaluated. Overprediction via fine-scale parameterization occurred in the center of the eddy, where the internal waves were inactive; however, the parameterization results were consistent with microstructure observations along the eddy’s periphery, where active internal waves existed. This indicates that the strong background shear and wave–mean flow interactions affected by the mesoscale eddy were not the main contributing factors that affected the applicability of fine-scale parameterization in the northern South China Sea. Instead, our results showed that the activity of internal waves is the most important consideration.

Balanced ellipsoidal vortex equilibria in a background shear flow at finite Rossby number

We consider a uniform ellipsoid of potential vorticity (PV), where we exploit analytical solutions derived for a balanced model at the second order in the Rossby number, the next order to quasi-geostrophic (QG) theory, the so-called QG+1 model. We consider this vortex in the presence of an external background shear flow, acting as a proxy for the effect of external vortices. For the QG model the system depends on four parameters, the height-to-width aspect ratio of the vortex, $h/r$ , as well as three parameters characterising the background flow, the strain rate, $\gamma$ , the ratio of the background rotation rate to the strain, $\beta$ , and the angle from which the flow is applied, $\theta$ . However, the QG+1 model also depends on the PV, as well as the Prandtl ratio, $f/N$ ( $f$ and $N$ are the Coriolis and buoyancy frequencies, respectively). For QG and QG+1 we determine equilibria for different values of the background flow parameters for increasing values of the imposed strain rate up to the critical strain rate, $\gamma _c$ , beyond which equilibria do not exist. We also compute the linear stability of this vortex to second-order modes, determining the marginal strain $\gamma _m$ at which ellipsoidal instability erupts. The results show that for QG+1 the most resilient cyclonic ellipsoids are slightly prolate, while anticyclonic ellipsoids tend to be more oblate. The highest values of $\gamma _m$ occur as $\beta \to 1$ . For large values of $f/N$ , changes in the marginal strain rates occur, stabilising anticyclonic ellipsoids while destabilising cyclonic ellipsoids.

Phase Structure of Internal Gravity Waves in the Ocean with Shear Flows

Purpose. The description of the internal gravity waves dynamics in the ocean with background fields of shear currents is a very difficult problem even in the linear approximation. The mathematical problem describing wave dynamics is reduced to the analysis of a system of partial differential equations; and while taking into account the vertical and horizontal inhomogeneity, this system of equations does not allow separation of the variables. Application of various approximations makes it possible to construct analytical solutions for the model distributions of buoyancy frequency and background shear ocean currents. The work is aimed at studying dynamics of internal gravity waves in the ocean with the arbitrary and model distributions of density and background shear currents. Methods and Results. The paper represents the numerical and analytical solutions describing the main phase characteristics of the internal gravity wave fields in the stratified ocean of finite depth, both for arbitrary and model distributions of the buoyancy frequency and the background shear currents. The currents are considered to be stationary and horizontally homogeneous on the assumption that the scale of the currents' horizontal and temporal variability is much larger than the characteristic lengths and periods of internal gravity waves. Having been used, the Fourier method permitted to obtain integral representations of the solutions under the Miles – Howard stability condition is fulfilled. To solve the vertical spectral problem, proposed is the algorithm for calculating the main dispersion dependences that determine the phase characteristics of the generated wave fields. The calculations for one real distribution of buoyancy frequency and shear flow profile are represented. Transformation of the dispersion surfaces and phase structures of the internal gravitational waves’ fields is studied depending on the generation parameters. To solve the problem analytically, constant distribution of the buoyancy frequency and linear dependences of the background shear current on depth were used. For the model distribution of the buoyancy and shear flow frequencies, the explicit analytical expressions describing the solutions of the vertical spectral problem were derived. The numerical and asymptotic solutions for the characteristic oceanic parameters were compared. Conclusions. The obtained results show that the asymptotic constructions using the model dependences of the buoyancy frequency and the background shear velocities’ distribution, describe the numerical solutions of the vertical spectral problem to a good degree of accuracy. The model representations, having been applied for hydrological parameters, make it possible to describe qualitatively correctly the main characteristics of internal gravity waves in the ocean with the arbitrary background shear currents.

Phase Structure of Internal Gravity Waves in the Ocean with Shear Flows

Purpose. The description of the internal gravity waves dynamics in the ocean with background fields of shear currents is a very difficult problem even in the linear approximation. The mathematical problem describing wave dynamics is reduced to the analysis of a system of partial differential equations; and while taking into account the vertical and horizontal inhomogeneity, this system of equations does not allow separation of the variables. Application of various approximations makes it possible to construct analytical solutions for the model distributions of buoyancy frequency and background shear ocean currents. The work is aimed at studying dynamics of internal gravity waves in the ocean with the arbitrary and model distributions of density and background shear currents. Methods and Results. The paper represents the numerical and analytical solutions describing the main phase characteristics of the internal gravity wave fields in the stratified ocean of finite depth, both for arbitrary and model distributions of the buoyancy frequency and the background shear currents. The currents are considered to be stationary and horizontally homogeneous on the assumption that the scale of the currents' horizontal and temporal variability is much larger than the characteristic lengths and periods of internal gravity waves. Having been used, the Fourier method permitted to obtain integral representations of the solutions under the Miles – Howard stability condition is fulfilled. To solve the vertical spectral problem, proposed is the algorithm for calculating the main dispersion dependences that determine the phase characteristics of the generated wave fields. The calculations for one real distribution of buoyancy frequency and shear flow profile are represented. Transformation of the dispersion surfaces and phase structures of the internal gravitational waves’ fields is studied depending on the generation parameters. To solve the problem analytically, constant distribution of the buoyancy frequency and linear dependences of the background shear current on depth were used. For the model distribution of the buoyancy and shear flow frequencies, the explicit analytical expressions describing the solutions of the vertical spectral problem were derived. The numerical and asymptotic solutions for the characteristic oceanic parameters were compared. Conclusions. The obtained results show that the asymptotic constructions using the model dependences of the buoyancy frequency and the background shear velocities’ distribution, describe the numerical solutions of the vertical spectral problem to a good degree of accuracy. The model representations, having been applied for hydrological parameters, make it possible to describe qualitatively correctly the main characteristics of internal gravity waves in the ocean with the arbitrary background shear currents.

Semi-Dilute Dumbbells: Solutions of the Fokker–Planck Equation

Spring bead models are commonly used in the constitutive equations for polymer melts. One such model based on kinetic theory—the finitely extensible nonlinear elastic dumbbell model incorporating a Peterlin closure approximation (FENE-P)—has previously been applied to study concentration-dependent anisotropy with the inclusion of a mean-field term to account for intermolecular forces in dilute polymer solutions for background profiles of weak shear and elongation. These investigations involved the solution of the Fokker–Planck equation incorporating a constitutive equation for the second moment. In this paper, we extend this analysis to include the effects of large background shear and elongation beyond the Hookean regime. Further, the constitutive equation is solved for the probability density function which permits the computation of any macroscopic variable, allowing direct comparison of the model predictions with molecular dynamics simulations. It was found that if the concentration effects at equilibrium are taken into account, the FENE-P model gives qualitatively the correct predictions, although the over-shoot in extension in comparison to the infinitely dilute case is significantly underpredicted.

Diffusive Staircases in Shear: Dynamics and Heat Transport

Arctic staircases mediate the heat transport from the warm water of Atlantic origin to the cooler waters of the Arctic mixed layer. For this reason, staircases have received much due attention from the community, and their heat transport has been well characterized for systems in the absence of external forcing. However, the ocean is a dynamic environment with large-scale currents and internal waves being omnipresent, even in regions shielded by sea-ice. Thus, we have attempted to address the effects of background shear on fully developed staircases using numerical simulations. The code, which is pseudo-spectral, evolves the governing equations for a Boussinesq fluid with temperature and salinity in a shearing coordinate system. We find that—– unlike many other double-diffusive systems—the sheared staircase requires three-dimensional simulations to properly capture the dynamics. Our simulations predict shear patterns that are consistent with observations and show that staircases in the presence of external shear should be expected to transport heat and salt at least twice as efficiently as in the corresponding non-sheared systems. These findings may lead to critical improvements in the representation of micro-scale mixing in global climate models.

Multi-Scale Localized Perturbation Method in OpenFOAM

A modified set of governing differential equations for geophysical fluid flows is derived. All of the simulation fields are decomposed into a nominal large-scale background state and a small-scale perturbation from this background, and the new system is closed by the assumption that the perturbation is one-way coupled to the background. The decomposition method, termed the multi-scale localized perturbation method (MSLPM), is then applied to the governing equations of stratified fluid flows, implemented in OpenFOAM, and exercised in order to simulate the interaction of a vertically-varying background shear flow with an axisymmetric perturbation in a turbulent ocean environment. The results demonstrate that the MSLPM can be useful in visualizing the evolution of a perturbation within a complex background while retaining the complex physics that are associated with the original governing equations. The simulation setup may also be simplified under the MSLPM framework. Further applications of the MSLPM, especially to multi-scale simulations that encompass a large range of spatial and temporal scales, may be beneficial for researchers.

Analytical Approximations of Dispersion Relations for Internal Gravity Waves Equation with Shear Flows

The problem of internal gravity waves fields in a stratified medium of finite depth is considered for model distributions of background shear currents. For the analytical solution of the problem, a constant distribution of the buoyancy frequency and various linear dependences of the background shear current on depth were used. The dispersion dependences are obtained, which are expressed in terms of the modified Bessel function of the imaginary index. Under the Miles–Howard stability condition and large Richardson numbers, the Debye asymptotics of the modified Bessel function of the imaginary index were used to construct analytical solutions. The dispersion equation is solved using the proposed analytical approximation. The properties of the dispersion equation are studied and the main analytical characteristics of the dispersion curves are investigated depending on the parameters of background shear flows.