The University of St Andrews

Research@StAndrews:FullText >
University of St Andrews Research >
University of St Andrews Research >
University of St Andrews Research >

Please use this identifier to cite or link to this item:
This item has been viewed 21 times in the last year. View Statistics

Files in This Item:

File Description SizeFormat
final.pdf739.12 kBAdobe PDFView/Open
Title: Numerical simulation of shear-induced instabilities in internal solitary waves
Authors: Carr, Magda
King, Stuart Edward
Dritschel, David Gerard
Keywords: Internal waves
Solitary waves
Stratified flows
QA Mathematics
Issue Date: 25-Sep-2011
Citation: Carr , M , King , S E & Dritschel , D G 2011 , ' Numerical simulation of shear-induced instabilities in internal solitary waves ' Journal of Fluid Mechanics , vol 683 , pp. 263-288 . , 10.1017/jfm.2011.261
Abstract: A numerical method that employs a combination of contour advection and pseudo-spectral techniques is used to simulate shear-induced instabilities in an internal solitary wave (ISW). A three-layer configuration for the background stratification, in which a linearly stratified intermediate layer is sandwiched between two homogeneous ones, is considered throughout. The flow is assumed to satisfy the inviscid, incompressible, Oberbeck–Boussinesq equations in two dimensions. Simulations are initialized by fully nonlinear, steady-state, ISWs. The results of the simulations show that the instability takes place in the pycnocline and manifests itself as Kelvin–Helmholtz billows. The billows form near the trough of the wave, subsequently grow and disturb the tail. Both the critical Richardson number (Ric) and the critical amplitude required for instability are found to be functions of the ratio of the undisturbed layer thicknesses. It is shown, therefore, that the constant, critical bound for instability in ISWs given in Barad & Fringer (J. Fluid Mech., vol. 644, 2010, pp. 61–95), namely Ric = 0.1 ± 0.01 , is not a sufficient condition for instability. It is also shown that the critical value of Lx/λ required for instability, where Lx is the length of the region in a wave in which Ri < 1/4 and λ is the half-width of the wave, is sensitive to the ratio of the layer thicknesses. Similarly, a linear stability analysis reveals that δiTw (where δi is the growth rate of the instability averaged over Tw, the period in which parcels of fluid are subjected to Ri < 1/4) is very sensitive to the transition between the undisturbed pycnocline and the homogeneous layers, and the amplitude of the wave. Therefore, the alternative tests for instability presented in Fructus et al. (J. Fluid Mech., vol. 620, 2009, pp. 1–29) and Barad & Fringer (J. Fluid Mech., vol. 644, 2010, pp. 61–95), respectively, namely Lx/λ ≥ 0.86 and δiTw > 5 , are shown to be valid only for a limited parameter range.
Version: Postprint
Description: This work was supported by the UK Engineering and Physical Sciences Research Council [grant number EP/F030622/1]
Status: Peer reviewed
ISSN: 0022-1120
Type: Journal article
Rights: This is the author's version of this article. The published version (c)Cambridge University Press is available from
Appears in Collections:University of St Andrews Research
Applied Mathematics Research
Scottish Oceans Institute Research

This item is protected by original copyright

Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.


DSpace Software Copyright © 2002-2012  Duraspace - Feedback
For help contact: | Copyright for this page belongs to St Andrews University Library | Terms and Conditions (Cookies)