Numerical study of a round buoyant jet under the effect of JONSWAP random waves

This paper presents a numerical study on the hydrodynamic behaviours of a round buoyant jet under the effect of JONSWAP random waves. A three-dimensional large eddy simulation (LES) model is developed to simulate the buoyant jet in a stagnant ambient and JONSWAP random waves. By comparison of velocity and concentration fields, it is found that the buoyant jet exhibits faster decay of centerline velocity, wider lateral spreading and larger initial dilution under the wave effect, indicating that wave dynamics improves the jet entrainment and mixing in the near field, and subsequently mitigate the jet impacts in the far field. The effect of buoyancy force on the jet behaviours in the random waves is also numerically investigated. The results show that the wave effect on the jet entrainment and mixing is considerably weakened under the existence of buoyancy force, resulting in a slower decay rate of centerline velocity and a narrower jet width for the jet with initial buoyancy.     

A note on the horizontal momentum exchange in combined waves and current

The horizontal exchange of momentum due to the organized motion in combined waves and current has been analyzed. The combination of the vertical orbital wave motion and the mean current gives a periodic variation in the horizontal velocity in addition to the wave orbital motion. This periodic variation, combined with the wave orbital motion, gives a significant contribution to the momentum exchange. Two examples are considered, the interaction of a pure wave motion and a current normal to the direction of wave propagation, and a wave driven longshore current with an undertow velocity profile. It is demonstrated that the new contribution changes the resulting momentum exchange considerably.     

Numerical study on three-dimensional waves produced by a bottom jet

The eruption of an underwater volcano can initiate the disturbances of the sea surface and subsequently generate a group of outward-propagating tsunamis. The theme of this study is to introduce a three-dimensional (3D) fully nonlinear wave model for the simulation of wave motions induced by a bottom jet. A boundary-fitted coordinate system is utilized to conveniently adjust grids according to the transient moving free surface. The governing Laplace equation of the velocity potential is solved by an implicit finite-difference scheme while a mixed explicit/implicit iteration procedure is applied to solve the free-surface boundary conditions. In addition, a set of generalized Boussinesq equations are solved for comparison with the fully nonlinear model. Good agreements in comparisons with the existing numerical and analytical solutions are achieved for cases investigated. Waves induced by three types of bottom jets: namely (1) sudden eruption, (2) initial transient, and (3) periodic transient are discussed in this paper. For the case of sudden erupted jet, a system of 3D outgoing waves as the cylindrical wave pattern are presented and discussed. For the initial transient types, it shows the transition in the incipient stage has a great influence on the initial rising of the water surface and the induced leading waves. Furthermore, an interesting up-down phenomenon in the center of disturbed free surface due to the type of periodic jet is revealed.     

Note on conservation equations for nonlinear surface waves

Euler's equations of motion in conjunction with the dynamic boundary condition are manipulated to obtain exact (and approximate) alternative momentum equations for nonlinear irrotational surface waves. The Airy and Boussinesq equations are re-derived as demonstrative examples. A fully nonlinear version of the improved Boussinesq equations is presented as a new application of the proposed equations. Further use of the equations in developing depth-integrated wave models, which are not necessarily restricted to finite depths, is also pointed out.     

A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

Numerical Modeling of a Round Jet Discharged into Random Waves

Coastal disposal of waste water can be idealized as the problem of a jet under random waves. Understanding of this phenomenon is important for engineering design and environmental impact assessment. The present study aims to simulate such phenomenon by using a 3D numerical model based on the solution of the spatially filtered and σ-transformed Navier–Stokes equations with dynamic sub-grid scale model of turbulence. The numerical solution procedures are split into three steps: advection, diffusion and pressure propagation, and a Lagrange–Euler method is used to track the free surface. Cases of vertical jet in stagnant water, pure random waves and vertical jet in random waves are simulated with the same grid system for comparative study. Different methods of generating jet inflow turbulence have been tested and the method of jet azimuthal modes is found to be the optimum. The numerical results reproduce the distinct characteristics of jet in waves, including faster decay of centerline velocity, wider lateral spreading and the occurrence of wave tractive mechanism.     

Surface manifestations of internal waves investigated by a subsurface buoyant Jet: 1. The mechanism of internal-wave generation

In a large test reservoir with artificial temperature stratification at the Institute of Applied Physics, Russian Academy of Sciences, we have performed a major laboratory simulation of the nonstationary dynamics of buoyant turbulent jets generated by wastewater flows from underwater collector diffusers. The interaction of buoyant jets with the pycnocline leads to an active generation of internal waves. An analysis of the dependence of wave amplitude on the control parameter proportional to the rate of liquid flow from the collector diffuser has indicated that this dependence is adequately described by a function that is characteristic for the presence in the Hopf bifurcation system, which occurs for a soft actuation mode of self-oscillations of the globally instable mode. To check the conditions for the actuation of the globally instable mode, we have performed an auxiliary experiment in a small reservoir with a salt stratification formulated similar to the experiment in the big reservoir. Using the particle image velocimetry (PIV) method, we have measured the velocity field in the buoyant jet and constructed the profiles of transverse velocity in several sections. When the jet approaches the pycnocline, a counterflow is generated at the edges. A stability analysis for the resulting profiles of flow velocities performed by the method of normal modes has revealed that, for the jet portions with counterflow, the condition of absolute instability by the Briggs criterion for axisymmetric jet oscillations is satisfied, which testifies to the fact that the globally instable mode is actuated. The estimates for oscillation frequencies of the globally instable mode are well consistent quantitatively with the measured spectrum of jet oscillations.     

Initial Dilution of A Vertical Round Non-Buoyant Jet in Wavy Cross-Flow Environment

The phenomenon of wastewater discharged into coastal waters can be simplified as a turbulent jet under the effect of waves and currents. Previous studies have been carried out to investigate the jet behaviors under the current only or the wave only environment. To obtain better understanding of the jet behaviors in a realistic situation, a series of physical experiments on the initial dilution of a vertical round jet in the wavy cross-flow environment are conducted. The diluted processes of the jet are recorded by a high-resolution camcorder and the concentration fields of the jet are measured with a peristaltic suction pumping system. When the jet is discharged into the wavy cross-flow environment, a distinctive phenomenon, namely "effluent clouds", is observed. According to the quantitative measurements, the jet width in the wavy cross-flow environment increases more significantly than that does in the cross-flow only environment, indicating that the waves impose a positive effect on the enhancement of jet initial dilution. In order to generalize the experimental findings, a comprehensive velocity scale ua and a characteristic length scale l are introduced. Through dimensional analysis, it is found that the dimensionless centerline concentration trajectories cy/l is in proportion to 1/3 power of the dimensionless downstream distance x/l, and the dimensionless centerline dilution 2c aS Q/(u l) is proportional to the square of the dimensionless centerline trajectory cy/l. Several empirical equations are then derived by using the Froude number of cross-flow Frc as a reference coefficient. This paper provides a better understanding and new estimations of the jet initial dilution under the combined effect of waves and cross-flow current.     

Modeling wave�Ccurrent bottom boundary layers beneath shoaling and breaking waves

The boundary layer characteristics beneath waves transforming on a natural beach are affected by both waves and wave-induced currents, and their predictability is more difficult and challenging than for those observed over a seabed of uniform depth. In this research, a first-order boundary layer model is developed to investigate the characteristics of bottom boundary layers in a wave–current coexisting environment beneath shoaling and breaking waves. The main difference between the present modeling approach and previous methods is in the mathematical formulation for the mean horizontal pressure gradient term in the governing equations for the cross-shore wave-induced currents. This term is obtained from the wave-averaged momentum equation, and its magnitude depends on the balance between the wave excess momentum flux gradient and the hydrostatic pressure gradient due to spatial variations in the wave field of propagating waves and mean water level fluctuations. A turbulence closure scheme is used with a modified low Reynolds number k-ε model. The model was validated with two published experimental datasets for normally incident shoaling and breaking waves over a sloping seabed. For shoaling waves, model results agree well with data for the instantaneous velocity profiles, oscillatory wave amplitudes, and mean velocity profiles. For breaking waves, a good agreement is obtained between model and data for the vertical distribution of mean shear stress. In particular, the model reproduced the local onshore mean flow near the bottom beneath shoaling waves, and the vertically decreasing pattern of mean shear stress beneath breaking waves. These successful demonstrations for wave–current bottom boundary layers are attributed to a novel formulation of the mean pressure gradient incorporated in the present model. The proposed new formulation plays an important role in modeling the boundary layer characteristics beneath shoaling and breaking waves, and ensuring that the present model is applicable to nearshore sediment transport and morphology evolution.     

Local balance in the air-sea boundary processes

A new growth equation for wind waves of simple spectrum is presented upon three basic concepts. The period and the wave height of significant waves in dimensionless forms, which are considered to correspond to the peak frequency and the energy level, respectively, are used as representative quantities of wind waves. One of the three basic concepts is the concept of local balance, and the other two concern the acquisition of wave energy and the dissipation of wave energy, respectively. It is shown from some actual data that the equation, together with two universal constants concerning the acquisition and the dissipation of wave energy (B=6.2×10^?2 andK=2.16×10^?5, respectively), is applied universally to wide ranges of wind waves from those in a wind-wave tunnel to fully developed sea in the open ocean. “The three-second power law for wind waves of simple spectrum”, and a few relations as the lemmas, are derived, such that the mean surface transport due to the orbital motion of wind waves is always proportional to the friction velocity in wind, and that the steepness is inversely proportional to the root of the wave age. It is also derived that the portion of wind stress which directly enters the wind waves decreases exponentially with increasing wave age and is 7.5 % of the total wind stress for very young waves. Also, equations are presented as to the increase of momentum of drift current, and as to the supply of turbulent energy by wind waves into the upper ocean.     

Simulation of directional waves in a numerical basin by a desingularized integral equation approach

A numerical solution is developed to investigate the generation and propagation of small-amplitude water waves in a semi-infinite rectangular wave basin. The three-dimensional wave field is produced by the prescribed “snake-like” motion of an array of segmented wave generators located along the wall at one end of the tank. The solution technique is based on the boundary element approach and uses an appropriate three-dimensional Green function which explicitly satisfies the tank-wall boundary conditions. The Green function and its derivatives which appear in the integral equation formulation can be shown to be slowly convergent when the source and field points are in close proximity. Therefore, when computing the velocity potentials on the wave generators, the source points are chosen outside the fluid domain, thereby ensuring the rapid convergence of these functions and rendering the integral equations non-singular. Numerical results are shown which illustrate the influence of the various wavemaker and basin parameters on the generated wave field. Finally, the complete wave field produced by the diffraction of oblique waves by a vertical circular cylinder in a basin is presented.     

Flow under standing waves: Part 1. Shear stress distribution,energy flux and steady streaming

The conditions for energy flux, momentum flux and the resulting streaming velocity are analysed for standing waves formed in front of a fully reflecting wall. The exchange of energy between the outer wave motion and the near bed oscillatory boundary layer is considered, determining the horizontal energy flux inside and outside the boundary layer. The momentum balance, the mean shear stress and the resulting time averaged streaming velocities are determined. For a laminar bed boundary layer the analysis of the wave drift gives results similar to the original work of Longuet–Higgins from 1953. The work is extended to turbulent bed boundary layers by application of a numerical model. The similarities and differences between laminar and turbulent flow conditions are discussed, and quantitative results for the magnitude of the mean shear stress and drift velocity are presented. Full two-dimensional simulations of standing waves have also been made by application of a general purpose Navier–Stokes solver. The results agree well with those obtained by the boundary layer analysis. Wave reflection from a plane sloping wall is also investigated by using the same numerical model and by physical laboratory experiments. The phase shift of the reflected wave train is compared with theoretical and empirical models.     

Bottom friction and its effects on periodic long wave propagation

A new set of Boussinesq-type equations describing the free surface evolution and the corresponding depth-integrated horizontal velocity is derived with the bottom boundary layer effects included. Inside the boundary layer the eddy viscosity gradient model is employed to characterize Reynolds stresses and the eddy viscosity is further approximated as a linear function of the distance measured from the seafloor. Boundary-layer velocities are coupled with the irrotational velocity in the core region through boundary conditions. The leading order boundary layer effects on wave propagation appear in the depth-integrated continuity equation to account for the velocity deficit inside the boundary layer. This formulation is different from the conventional approach in which a bottom stress term is inserted in the momentum equation. An iterative scheme is developed to solve the new model equations for the free surface elevation, depth-integrated velocity, the bottom stress, the boundary layer thickness and the magnitude of the turbulent eddy viscosity. A numerical example for the evolution of periodic waves propagating in one-dimensional channel is discussed to illustrate the numerical procedure and physics involved. The differences between the conventional approach and the present formulation are discussed in terms of the bottom frictional stress and the free surface profiles.     

Boundary layer approach in the modeling of breaking solitary wave runup

The boundary layer is very important in the relation between wave motion and bed stress, such as sediment transport. It is a known fact that bed stress behavior is highly influenced by the boundary layer beneath the waves. Specifically, the boundary layer underneath wave runup is difficult to assess and thus, it has not yet been widely discussed, although its importance is significant. In this study, the shallow water equation (SWE) prediction of wave motion is improved by being coupled with the k–ω model, as opposed to the conventional empirical method, to approximate bed stress. Subsequently, the First Order Center Scheme and Monotonic Upstream Scheme of Conservation Laws (FORCE MUSCL), which is a finite volume shock-capturing scheme, is applied to extend the SWE range for breaking wave simulation. The proposed simultaneous coupling method (SCM) assumes the depth-averaged velocity from the SWE is equivalent to free stream velocity. In turn, free stream velocity is used to calculate a pressure gradient, which is then used by the k–ω model to approximate bed stress. Finally, this approximation is applied to the momentum equation in the SWE. Two experimental cases will be used to verify the SCM by comparing runup height, surface fluctuation, bed stress, and turbulent intensity values. The SCM shows good comparison to experimental data for all before-mentioned parameters. Further analysis shows that the wave Reynolds number increases as the wave propagates and that the turbulence behavior in the boundary layer gradually changes, such as the increase of turbulent intensity.     

有界赤道大洋波包解及其年际年代际变率

Linearized shallow water perturbation equations with approximation in an equatorial β plane are used to obtain the analytical solution of wave packet anomalies in the upper bounded equatorial ocean. The main results are as follows. The wave packet is a superposition of eastward travelling Kelvin waves and westward travelling Rossby waves with the slowest speed, and satisfies the boundary conditions of eastern and western coasts, respectively.The decay coefficient of this solution to the north and south sides of the equator is inversely proportional only to the phase velocity of Kelvin waves in the upper water. The oscillation frequency of the wave packet, which is also the natural frequency of the ocean, is proportional to its mode number and the phase velocity of Kelvin waves and is inversely proportional to the length of the equatorial ocean in the east-west direction. The flow anomalies of the wave packet of Mode 1 most of the time appear as zonal flows with the same direction. They reach the maximum at the center of the equatorial ocean and decay rapidly away from the equator, manifested as equatorially trapped waves. The flow anomalies of the wave packet of Mode 2 appear as the zonal flows with the same direction most of the time in half of the ocean, and are always 0 at the center of the entire ocean which indicates stagnation, while decaying away from the equator with the same speed as that of Mode 1. The spatial structure and oscillation period of the wave packet solution of Mode 1 and Mode 2 are consistent with the changing periods of the surface spatial field and time coefficient of the first and second modes of complex empirical orthogonal function(EOF)analysis of flow anomalies in the actual equatorial ocean. This indicates that the solution does exist in the real ocean, and that El Ni?o-Southern Oscillation(ENSO) and Indian Ocean dipole(IOD) are both related to Mode 2.After considering the Indonesian throughflow, we can obtain the length of bounded equatorial ocean by taking the sum of that of the tropical Indian Ocean and the tropical Pacific Ocean, thus this wave packet can also explain the decadal variability(about 20 a) of the equatorial Pacific and Indian Oceans.     

A Bingham-plastic model for fluid mud transport under waves and currents

Simplified equations of fluid mud motion, which is described as Bingham-Plastic model under waves and currents, are presented by order analysis. The simplified equations are non-linear ordinary differential equations which are solved by hybrid numerical-analytical technique. As the computational cost is very low, the effects of wave current parameters and fluid mud properties on the transportation velocity of the fluid mud are studied systematically. It is found that the fluid mud can move toward one direction even if the shear stress acting on the fluid mud bed is much smaller than the fluid mud yield stress under the condition of wave and current coexistence. Experiments of the fluid mud motion under current with fluctuation water surface are carried out. The fluid mud transportation velocity predicted by the presented mathematical model can roughly match that measured in experiments.     

Numerical simulation of wave-induced laminar boundary layers

The three-dimensional numerical model with σ-coordinate transformation in the vertical direction is applied to the simulation of surface water waves and wave-induced laminar boundary layers. Unlike most of the previous investigations that solved the simplified one-dimensional boundary layer equation of motion and neglected the interaction between boundary layer and outside flow, the present model solves the full Navier–Stokes equations (NSE) in the entire domain from bottom to free surface. A non-uniform mesh system is used in the vertical direction to resolve the thin boundary layer. Linear wave, Stokes wave, cnoidal wave and solitary wave are considered. The numerical results are compared to analytical solutions and available experimental data. The numerical results agree favorably to all of the experimental data. It is found that the analytical solutions are accurate for both linear wave and Stokes wave but inadequate for cnoidal wave or solitary wave. The possible reason is that the existing analytical solutions for cnoidal and solitary waves adopt the first-order approximation for free stream velocity and thus overestimate the near bottom velocity. Besides velocity, the present model also provides accurate results for wave-induced bed shear stress.     

Generation of surface and internal waves in a bounded basin by a moving baric front

We consider a plane problem of generation of surface and internal waves in a bounded rotating basin of variable depth by a front of atmospheric pressure moving over the basin. The fluid is assumed to be two-layer. The system of nonlinear long-wave equations is solved numerically by the method of finite differences for the distribution of depths corresponding to a zonal section of the Black-Sea basin. It is shown that the baric front moving over the basin generates barotropic and baroclinic oscillations of the fluid. The intensity of disturbances depends on the velocity of motion and the width of the front. There exists a velocity of motion of the front for which internal waves are generated especially efficiently. When the front leaves the basin, we observe the formation of a packet of internal waves propagating from one lateral boundary of the basin to the other boundary with reflections from the boundaries.     

The Mathematical and Parameterized Expressions for Wave Induced Excess Flow of Momentum

Wave induced excess flow of momentum(WIEFM)is the averaged flow of momentum over a wave period due to wave presence,which may also be called 3-D radiation stress.In this paper,the 3-D current equations with WIEFM are derived from the averaged Navier-Stokes equations over a wave period,in which the velocity is separated into the large-scale background velocity,the wave particle velocity and the turbulent fluctuation velocity.A concept of wave fluctuating layer(WFL)is put forward,which is the vertical column from the wave trough to wave ridge.The mathematical expressions of WIEFM in WFL and below WFL are given separately.The parameterized expressions of WIEFM are set up according to the linear wave theory.The integration of WIEFM in the vertical direction equals the traditional radiation stress(namely 2-D radiation stress)given by Longuet-Higgins and Stewart.