Baroclinic tides result from the interaction of barotropic tides with topography in stratified oceans. They play an important role in driving deep ocean mixing. In this research, investigations of the dynamics of baroclinic tides and internal solitary waves (ISWs) in the northern South China Sea (SCS) are conducted, mainly by means of the Massachusetts Institute of Technology general circulation model (MITgcm). Firstly, simulations of internal wave generation at the Luzon Strait (LS) are carried out. By conducting three-dimensional (3D), high-resolution experiments, it was found that the generated wave field features a multi-modal structure: large, pronounced ISWs of first mode (amplitude ~120 m) and second mode (amplitude ~120 m) were reproduced. The two north-south aligned ridges in the LS contribute together to the generation of the second mode ISWs, whereas the easternmost ridge of the two is responsible for the first mode ISWs. It was found that multiple generation mechanisms of internal waves could occur in this region, and overall it belongs to a mixed lee wave regime. A specific type of short internal waves arose during the 3D simulation. These ride on a second mode ISW with similar phase speed, trailing a first mode ISW. The short waves possess wavelengths of ~1.5 km and amplitudes of ~20 m, and only show up in the upper layer up to a depth of ~500 m. Scrutiny of the generation process showed that these short waves appear in two distinct regions and are produced due to two mechanisms, namely, the disintegration of an inclined baroclinic bore near the LS, and the overtaking of a second mode ISW in the deep water by a faster first mode ISW. Robust evidence has been sought from satellite imagery and by solving the theoretical Taylor-Goldstein Equation to verify their existence. The effects of superposition of multiple tidal harmonics (diurnal and semidiurnal) on the resultant ISW generation were investigated. It was first found that, by analyzing historical observational data, the occurrence of ISWs in the far-field always follow strong semidiurnal barotropic tidal peaks in the LS, regardless of whether it is the maximum for the diurnal or total tidal strength. However, modelling results of MITgcm and a linear internal tide generation model demonstrate that the diurnal tidal harmonics modulate the arrival time and amplitude of the propagating ISWs. Specifically, it leads to the emergence of the so-called A and B type ISWs and an alternation and transition between the two. Secondly, the shoaling process of ISWs in the northern SCS slope-shelf area is investigated. A series of two-dimensional (2D) experiments are set up to study the shoaling of a large-amplitude second mode concave ISW over a linear slope that resembles the SCS slope. Modelling results show that a strong transformation of the wave profile starts to take place when the wave is approaching the shelf break. A convex type wave is born at the trailing edge of the incident wave and gradually disintegrates into a group of ISWs due to the steepening of the rear wave profile. The frontal face of the wave gets flatter when travelling on the slope, but forms a steep structure right above the shelf break. However, this steep structure shows no tendency to evolve into an ISW: instead, it gets increasingly flat again while evolving on the shelf. The trailing convex wave packet travels faster and merges with the frontal concave wave. Finally, a wave packet with rank-ordered convex ISWs moves forward steadily on the shelf. Energy transfer to the ambient modes is evident, as both first mode and higher modes are clearly seen during and after the shoaling process. First mode ISW evolution is studied too by performing 3D, high-resolution experiments over the wide northern SCS slope and shelf area. It was found that the wave profiles change drastically near the shelf break and the Dongsha Atoll. In agreement with satellite imagery, the wavefront of the leading ISW becomes more spatially oblique with respect to its original orientation as it progresses westward due to the inclination of the slope in the topography. Wave disintegration is prominent in the shallow water zone, and wave polarity reverses near the turning point (at the 130 m isobath), which is consistent with the predictions of weakly nonlinear theory. A series of 2D experiments were set up to inspect the effects of rotation on the shoaling ISW. The results indicate that under the rotation, upon reaching the continental shelf, one shoaling ISW could disintegrate into one ISW packet and one secondary solibore that contains a number of rank-ordered waves with much shorter wavelength than an ISW. The secondary solibore is very pronounced in the northern portion of the northern SCS slope and shelf, but could hardly be discerned in the southern portion, which is consistent with the outcome of 3D simulations.

Due to advances in technology and sophistication of many efficient algorithms, accurate numerical results can be achieved by using highly efficient computational software for research in wave action on coastal and harbor structures. These advances have benefitted the research in the physical phenomena of internal wave generation, propagation and breaking, which are some of the important topics in oceanography. In this study, the Flow-3D CFD (Computational Fluid Dynamics) software is used to simulate internal solitary wave motion in a density stratified fluid, in which the upper and lower layers are fresh and brine water, respectively. An internal solitary wave (ISW) is produced numerically by gravitational collapse mechanism in a numerical wave flume of 0.7 x 0.5 x 12.0 m (height x width x length ). The ISW in depression is then allow to propagate and across four different bottom obstacles (long uniform slope, trapezoidal section with short platform and isosceles triangle), in order to explore its waveform evolution and flow field distribution. This study also describes the setting and operation of the Flow-3D software, its application to the internal wave experiment, as well as verification of the numerically simulated results using previous laboratory experimental data. In this study, the lifting speed for the sluice gate was vital for not only the amplitude of an ISW, but also the speed of wave propagation in the flume. The result showed that the faster the gate opening, the faster propagation speeds and larger amplitude for the ISW so generated. Conversely, a slower gate opening led to weak wave speed and small amplitude to an ISW. Upon analyzing the results, we have found that the relationship between the speed of the gate opening and the wave propagation speed are not linear. Moreover, preliminary analysis and discussion are given for the ISW propagation over an obstacle (uniform long slope, trapezoidal section with short platform, and isosceles triangle), particularly on waveform evolution, vortex motions and flow field variations. It is believed that we can gain a better and thorough understanding of the internal wave characteristics, compared to physical laboratory experiments, if the numerical tool is applied with very fine grids and detailed analysis on the numerical outputs.

The breaking of internal solitary waves (ISWs) of depression shoaling upon a uniformly sloping boundary in a smoothed two-layer density field was investigated using high-resolution two-dimensional simulations. The simulations were performed for a wide range of boundary slopes S∈[0.01,0.3] and wave slopes. Over steep slopes (S≥0.1), three distinct breaking processes were observed; surging, plunging and collapsing breakers which are associated with reflection, convective instability and boundary layer separation, respectively. Over mild slopes (S≤0.05), nonlinearity varies gradually and the wave fissions into a train of waves of elevation after it passes through the turning point where solitary waves reverse polarity. The dynamics of each breaker type were investigated and the predominance of a particular mechanism was associated with a relative developmental timescale. The breaker type was characterized in wave slope S_w versus S space and the reflection coefficient (R), modeled as a function of the internal Iribarren number, was in agreement with other studies. The same 2D model was applied to investigate boundary layer separation-driven global instability, which is shown to play an important role in breaking of shoaling ISWs. The simulations were conducted with waves propagating over a flat bottom and shoaling over relatively mild (S=0.05) and steep (S=0.1) slopes. Combining the results over flat and sloping boundaries, a unified criterion for vortex shedding is proposed, which depends on the momentum thickness Reynolds number and the non-dimensionalized ISW-induced pressure gradient at the point of separation. The criterion is generalized to a form that may be readily computed from field data and compared to published laboratory and field observations. During vortex shedding, the bed shear stress, vertical velocity and near-bed Reynolds stress were elevated, implying potential for sediment re-suspension. Laboratory experiments were also performed to study three-dimensionality (3D) of global instability. Our results agree with previous laboratory experiments, using the proposed pressure gradient parameter and Reynolds number. The 3D effects prevent the vortices from ascending as high as they do in 2D simulations. The instabilities were not able to re-suspend sediments with 20 µm median diameters, which suggests applying lighter sediments, as finer sediments will be cohesive and dynamically different. / Thesis (Ph.D, Civil Engineering) -- Queen's University, 2011-12-23 15:03:29.76

Les ondes internes sont observées partout dans le monde. Elles ont un rôle important dans la mise en place de la chaine trophique, et elles peuvent avoir un impact sur les activités humaines. Ainsi, dans le Golfe de Guinée, des ondes internes solitaires (OIS) ont provoqué des incidents sur des plateformes pétrolières. L’objectif de la thèse est d’identifier leurs zones de génération, leurs directions de propagation, ainsi que les conditions d’environnement susceptibles de les déclencher. La thèse a comporté deux aspects : d’une part des mesures de courant étaient disponibles et d’autre part la mise en place de modèles numériques. Ces mesures ont montré le passage d’OIS. En outre, les zones de génération potentielles sont identifiées en haut du talus continental. Leur répartition mensuelle a montré une variabilité saisonnière, due à la Mousson Ouest Africaine qui modifie les conditions de stratification au cours de l’année.La modélisation a été réalisée avec HYCOM. Une maquette hydrostatique régionale du Golfe de Guinée a été mise en place. Elle a d’abord été validée en configuration barotrope (océan homogène). Une configuration bicouche a ensuite été testée, et la théorie linéaire des ondes internes a été vérifiée sur ses résultats. Des termes non hydrostatiques ont ensuite été ajoutés dans les équations horizontales du courant de HYCOM, pour modéliser des OIS. L’impact des paramètres physiques, en particulier la force de Coriolis et le forçage barotrope, ainsi que l’effet d’un courant moyen, a été étudié sur des configurations bicouches radiales. Enfin, un canyon a été ajouté au milieu du talus pour évaluer qualitativement les effets d’une bathymétrie 2D sur les OIS. / Internal waves are observed in different locations. They have a key role in the set up of the trophic chain, and they can impact human activities.For instance, in the Gulf of Guinea, internal solitary waves (ISWs) have caused hazards on offshore oil platforms. The aim of the thesis is to identify the generation zones and the directions of propagation of these ISWs, and the environmental conditions (tides, stratification) that can generate them. One the one hand measurements of currents were available, and on the other hand numerical models were used.Measurements of currents were analyzed. Packets of ISWs were identified. Moreover, the potential generation zones of the ISWs, particularly the top of the shelf break, were inferred from their direction of propagation. The monthly repartition of the packets of ISWs revealed a seasonal variability, due to the West African Monsoon, which modified the conditions of the stratification throughout the year.Concerning the modelling, the numerical model HYCOM was used. One regional hydrostatic model of the Gulf of Guinea was built. This regional model was validated for a barotropic configuration (homogeneous ocean) first. Then, a two-layered configuration was set up, and the linear theory of the internal waves was confirmed on the results.Non hydrostatic terms were added to the horizontal velocity equations in HYCOM, to model ISWs. The influence of different physical parameters, such as the Coriolis force, the barotropic forcing, and a mean current was studied on two-layered radial configurations. Finally, a canyon in the middle of the shelf break was set up to assess qualitatively the effect of a 2D bathymetry on the ISWs.