The interaction of localised solitary waves with large-scale, time-varying dispersive mean flows subject to non-convex flux is studied in the framework of the modified Korteweg–de Vries (mKdV) equation, a canonical model for internal gravity wave propagation and potential vorticity fronts in stratified fluids. The effect of large amplitude, dynamically evolving mean flows on the propagation of localised waves – essentially ‘soliton steering’ by the mean flow – is considered. A recent theoretical and experimental study of this new type of dynamic soliton–mean flow interaction for convex flux has revealed two scenarios where the soliton either transmits through the varying mean flow or remains trapped inside it. In this paper, it is demonstrated that the presence of a non-convex cubic hydrodynamic flux introduces significant modifications to the scenarios for transmission and trapping. A reduced set of Whitham modulation equations is used to formulate a general mathematical framework for soliton–mean flow interaction with non-convex flux. Solitary wave trapping is stated in terms of crossing modulation characteristics. Non-convexity and positive dispersion – common for stratified fluids – imply the existence of localised, sharp transition fronts (kinks). Kinks play dual roles as a mean flow and a wave, imparting polarity reversal to solitons and dispersive mean flows, respectively. Numerical simulations of the mKdV equation agree with modulation theory predictions. The mathematical framework developed is general, not restricted to completely integrable equations like mKdV, enabling application beyond the mKdV setting to other fluid dynamic contexts subject to non-convex flux such as strongly nonlinear internal wave propagation that is prevalent in the ocean.
more »
« less
The sandpaper theory of flow–topography interaction for multilayer shallow-water systems
Seafloor roughness profoundly influences the pattern and dynamics of large-scale oceanic flows. However, these kilometre-scale topographic patterns are unresolved by global numerical Earth system models and will remain subgrid for the foreseeable future. To properly represent the effects of small-scale bathymetry in analytical and coarse-resolution numerical models, we develop the stratified ‘sandpaper’ theory of flow–topography interaction. This model, which is based on the multilayer shallow-water framework, extends its barotropic antecedent to stratified flows. The proposed theory is successfully tested on the configuration representing the interaction of a zonal current with a corrugated cross-flow ridge.
more »
« less
- Award ID(s):
- 2241625
- PAR ID:
- 10574263
- Publisher / Repository:
- CUP
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 988
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Recent studies reveal the dramatic impact of seafloor roughness on the dynamics and stability of broad oceanic flows. These findings motivate the development of parameterizations that concisely represent the effects of small-scale bathymetric patterns in theoretical and coarse-resolution numerical circulation models. The previously reported quasi-geostrophic ‘sandpaper’ theory of flow–topography interactiona prioriassumes gentle topographic slopes and weak flows with low Rossby numbers. Since such conditions are often violated in the ocean, we now proceed to formulate a more general model based on shallow-water equations. The new version of the sandpaper model is validated by comparing roughness-resolving and parametric simulations of the flow over a corrugated seamount.more » « less
-
The interaction between deformable surfaces and oscillatory driving is known to produce complex secondary time-averaged flows due to inertial and elastic nonlinearities. Here, we revisit the problem of oscillatory flow in a cylindrical tube with a deformable wall, and analyse it under a long-wave theory for small deformations, but for arbitrary Womersley numbers. We find that the oscillatory pressure does not vary linearly along the length of a deformable channel, but instead decays exponentially with spatial oscillations. We show that this decay occurs over an elasto-visco-inertial length scale that depends on the material properties of the fluid and the elastic walls, the geometry of the system, and the frequency of the oscillatory flow, but is independent of the amplitude of deformation. Inertial and geometric nonlinearities associated with the elastic deformation of the channel drive a time-averaged secondary flow. We quantify the flow using numerical solutions of the perturbation theory, and gain insight by developing analytic approximations. The theory identifies a complex non-monotonic dependence of the time-averaged flux on the elastic compliance and inertia, including a reversal of the flow. Finally, we show that our analytic theory is in excellent quantitative agreement with the three-dimensional direct numerical simulations of Pandeet al.(Phys. Rev. Fluids, vol. 8, no. 12, 2023, 124102).more » « less
-
Simulations of strongly stratified turbulence often exhibit coherent large-scale structures called vertically sheared horizontal flows (VSHFs). VSHFs emerge in both two-dimensional (2D) and three-dimensional (3D) stratified turbulence with similar vertical structure. The mechanism responsible for VSHF formation is not fully understood. In this work, the formation and equilibration of VSHFs in a 2D Boussinesq model of stratified turbulence is studied using statistical state dynamics (SSD). In SSD, equations of motion are expressed directly in the statistical variables of the turbulent state. Restriction to 2D turbulence facilitates application of an analytically and computationally attractive implementation of SSD referred to as S3T, in which the SSD is expressed by coupling the equation for the horizontal mean structure with the equation for the ensemble mean perturbation covariance. This second-order SSD produces accurate statistics, through second order, when compared with fully nonlinear simulations. In particular, S3T captures the spontaneous emergence of the VSHF and associated density layers seen in simulations of turbulence maintained by homogeneous large-scale stochastic excitation. An advantage of the S3T system is that the VSHF formation mechanism, which is wave–mean flow interaction between the emergent VSHF and the stochastically excited large-scale gravity waves, is analytically understood in the S3T system. Comparison with fully nonlinear simulations verifies that S3T solutions accurately predict the scale selection, dependence on stochastic excitation strength, and nonlinear equilibrium structure of the VSHF. These results constitute a theory for VSHF formation applicable to interpreting simulations and observations of geophysical examples of turbulent jets such as the ocean’s equatorial deep jets.more » « less
-
Abstract Flow over a surface can be stratified by imposing a fixed mean vertical temperature (density) gradient profile throughout or via cooling at the surface. These distinct mechanisms can act simultaneously to establish a stable stratification in a flow. Here, we perform a series of direct numerical simulations of open-channel flows to study adaptation of a neutrally stratified turbulent flow under the combined or independent action of the aforementioned mechanisms. We force the fully developed flow with a constant mass flow rate. This flow forcing technique enables us to keep the bulk Reynolds number constant throughout our investigation and avoid complications arising from the acceleration of the bulk flow if a constant pressure gradient approach were to be adopted to force the flow instead. When both stratification mechanisms are active, the dimensionless stratification perturbation number emerges as an external flow control parameter, in addition to the Reynolds, Froude, and Prandtl numbers. We demonstrate that significant deviations from the Monin–Obukhov similarity formulation are possible when both types of stratification mechanisms are active within an otherwise weakly stable flow, even when the flux Richardson number is well below 0.2. An extended version of the similarity theory due to Zilitinkevich and Calanca shows promise in predicting the dimensionless shear for cases where both types of stratification mechanisms are active, but the extended theory is less accurate for gradients of scalar. The degree of deviation from neutral dimensionless shear as a function of the vertical coordinate emerges as a qualitative measure of the strength of stable stratification for all the cases investigated in this study.more » « less
An official website of the United States government

