Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

Wave generation by the falling rock in the two-dimensional wave tank

Wave generation by the falling rock in the two-dimensional wave tank is experimentally and numerically studied, where the numerical model utilizes the boundary element method to solve the fully nonlinear potential flow theory. The wave profiles at different times are measured in the laboratory, which are also used to test the numerical model. Comparisons show that the experimental and numerical results are in good agreement, and the numerical model can be used to simulate the wave generation due to the submarine rock falling. Further numerical tests on the influences of the rock size, density, initial position and the falling angle on the wave elevation of the generated waves are performed, respectively. The results show that the size and density of the rock have strong effects on the maximum elevation of the generated wave, while the effects of the initial position and the falling angle of the rock are also significant. When the size or the density of the rock increases, the maximum elevation of the generated wave increases. The same effect on the generated wave would be produced if the initial position of the rock becomes closer to the surface, or the falling angle between the falling route and the vertical direction turns larger. In addition, the present numerical tests reveal that the submarine rock falling provides a new generation method for the breaking wave in the wave tank.

Determination of fractional energy loss of waves in nearshore waters using an improved high-order Boussinesq-type model

Fractional energy losses of waves due to wave breaking when passing over a submerged bar are studied systematically using a modified numerical code that is based on the high-order Boussinesq-type equations. The model is first tested by the additional experimental data, and the model's capability of simulating the wave transformation over both gentle slope and steep slope is demonstrated. Then, the model's breaking index is replaced and tested. The new breaking index, which is optimized from the several breaking indices, is not sensitive to the spatial grid length and includes the bottom slopes. Numerical tests show that the modified model with the new breaking index is more stable and efficient for the shallow-water wave breaking. Finally, the modified model is used to study the fractional energy losses for the regular waves propagating and breaking over a submerged bar. Our results have revealed that how the nonlinearity and the dispersion of the incident waves as well as the dimensionless bar height (normalized by water depth) dominate the fractional energy losses. It is also found that the bar slope (limited to gentle slopes that less than 1:10) and the dimensionless bar length (normalized by incident wave length) have negligible effects on the fractional energy losses.

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42,134-149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313-351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313-351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971-5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift. (C) 2008 Elsevier Ltd. All rights reserved.

Influence of mean water depth and a subsurface sandbar on the onset and strength of wave breaking

Wave breaking in the open ocean and coastal zones remains an intriguing yet incompletely understood process, with a strong observed association with wave groups. Recent numerical study of the evolution of fully nonlinear, two-dimensional deep water wave groups identified a robust threshold of a diagnostic growth-rate parameter that separated nonlinear wave groups that evolved to breaking from those that evolved with recurrence. This paper investigates whether these deep water wave-breaking results apply more generally, particularly in finite-water-depth conditions. For unforced nonlinear wave groups in intermediate water depths over a flat bottom, it was found that the upper bound of the diagnostic growth-rate threshold parameter established for deep water wave groups is also applicable in intermediate water depths, given by k(0) h greater than or equal to 2, where k(0) is the mean carrier wavenumber and h is the mean depth. For breaking onset over an idealized circular arc sandbar located on an otherwise flat, intermediate-depth (k(0) h greater than or equal to 2) environment, the deep water breaking diagnostic growth rate was found to be applicable provided that the height of the sandbar is less than one-quarter of the ambient mean water depth. Thus, for this range of intermediate-depth conditions, these two classes of bottom topography modify only marginally the diagnostic growth rate found for deep water waves. However, when intermediate-depth wave groups ( k(0) h greater than or equal to 2) shoal over a sandbar whose height exceeds one-half of the ambient water depth, the waves can steepen significantly without breaking. In such cases, the breaking threshold level and the maximum of the diagnostic growth rate increase systematically with the height of the sandbar. Also, the dimensions and position of the sandbar influenced the evolution and breaking threshold of wave groups. For sufficiently high sandbars, the effects of bottom topography can induce additional nonlinearity into the wave field geometry and associated dynamics that modifies the otherwise robust deep water breaking-threshold results.

Applications of improved 2-D numerical wave tanks in wave breaking

Song and Banner (2002, henceforth referred to as SB02) used a numerical wave tank (developed by Drimer and Agnon, and further refined by Segre, henceforth referred to as DAS) to study the wave breaking in the deep water, and proposed a dimensionless breaking threshold that based on the behaviour of the wave energy modulation and focusing during the evolution of the wave group. In this paper, two modified DAS models are used to further test the SB02's results, the first one (referred to MDAS1) corrected many integral calculation errors appeared in the DAS code, and the second one (referred to MDAS2) replaced the linear boundary element approximation of DAS into the cubic element on the free surface. Researches show that the results of MDAS1 are the same with those of DAS for the simulations of deep water wave breaking, but, the different values of the wavemaker amplitude, the breaking time and the maximum local average energy growth rate delta(max) for the marginal breaking cases are founded by MDAS2 and MDAS1. However, MDAS2 still satisfies the SB02' s breaking threshold. Furthermore, MDAS1 is utilized to study the marginal breaking case in the intermediate water depth when wave passes over a submerged slope, where the slope is given by 1 : 500, 1 : 300, 1 : 150 or 1 : 100. It is found that the maximum local energy density U increases significantly if the slope becomes steeper, and the delta(max) decreases weakly and increases intensively for the marginal recurrence case and marginal breaking case respectively. SB02's breaking threshold is still valid for the wave passing over a submerged slope gentler than 1 : 100 in the intermediate water depth.

Numerical Simulation of Waves Generated by Seafloor Movements

Waves generated by vertical seafloor movements are simulated by use of a fully nonlinear two-dimensional numerical wave tank. In the source region, the seafloor lifts to a designated height by a generation function. The numerical tests show that file linear theory is only valid for estimating the wave behaviors induced by the seafloor movements with a small amplitude, and the fully nonlinear numerical model should be adopted in the simulation of the wave generation by the large amplitude seafloor movements. Without the background surface waves, many numerical tests on the stable maximum elevations eta(max)(0) are carried out by both the linear theory and the fully nonlinear model. The results of two models are compared and analyzed. For the fully nonlinear model, the influences of the amplitudes and the horizontal lengths on eta(max)(0) are stronger than that of the characteristic duration times. Furthermore, results reveal that there are significant differences between the linear theory and the fully nonlinear model. When the influences of the background surface waves are considered, the corresponding numerical analyses reveal that with the fully nonlinear model the eta(max)(0) near-linearly varies with the wave amplitudes of the surface waves, and the eta(max)(0) has significant dependences on the wave lengths and the wave phases of the surface waves. In addition, the differences between the linear theory and the fully nonlinear model are still obvious, aid these differences are significantly affected by The wave parameters of the background surface waves, such as the wave amplitude, the wave length and the wave phase.

Micro-shock waves generated inside a fluid jet impinging on plane wall

An Nd:glass laser pulse (18 ns, 1.38 J) is focused in a tiny area of about 100-mum diam under ambient conditions to produce micro-shock waves. The laser is focused above a planar surface with a typical standoff distance of about 4 mm, The laser energy is focused inside a supersonic circular jet of carbon dioxide gas produced by a nozzle with internal diameter of 2.9 mm and external diameter of 8 mm, Nominal value of the Mach number of the jet is around 2 with the corresponding pressure ratio of 7.5 (stagnation pressure/static pressure at the exit of the nozzle), The interaction process of the micro-shock wave generated inside the supersonic jet with the plane wall is investigated using double-pulse holographic interferometry. A strong surface vortex field with subsequent generation of a side jet propagating outward along the plane wail is observed. The interaction of the micro-shock wave with the cellular structure of the supersonic jet does not seem to influence the near surface features of the flowfield. The development of the coherent structures near the nozzle exit due to the upstream propagation of pressure waves seems to be affected by the outward propagating micro-shock wave. Mach reflection is observed when the micro-shock wave interacts with the plane wall at a standoff distance of 4 mm, The Mach stem is slightly deflected, indicating strong boundary-layer and viscous effects near the wall. The interaction process is also simulated numerically using an axisymmetric transient laminar Navier-Stokes solver. Qualitative agreement between experimental and numerical results is good.

Dynamic Characteristics of Spherically Converging Detonation Waves

The spherically converging detonation wave was numerically investigated by solving the one-dimensional multi-component Euler equations in spherical coordinates with a dispersion-controlled dissipative scheme. Finite rate and detailed chemical reaction models were used and numerical solutions were obtained for both a spherical by converging detonation in a stoichiometric hydrogen-oxygen mixture and a spherically focusing shock in air. The results showed that the post-shock pressure approximately arises to the same amplitude in vicinity of the focal point for the two cases, but the post-shock temperature level mainly depends on chemical reactions and molecular dissociations of a gas mixture. While the chemical reaction heat plays an important role in the early stage of detonation wave propagation, gas dissociations dramatically affect the post-shock flow states near the focal point. The maximum pressure and temperature, non-dimensionalized by their initial value, are approximately scaled to the propagation radius over the initial detonation diameter. The post-shock pressure is proportional to the initial pressure of the detonable mixture, and the post-shock temperature is also increased with the initial pressure, but in a much lower rate than that of the post-shock pressure.

An Analytical Solution for Three-Dimensional Diffraction of Plane P-Waves by a Hemispherical Alluvial Valley with Saturated Soil Deposits

An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot's dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot's dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.

The Interaction Between Shock Waves and Foam in a Shock Tube

An experimental study and a numerical simulation were conducted to investigate the mechanical and thermodynamic processes involved in the interaction between shock waves and low density foam. The experiment was done in a stainless shock tube (80mm in inner diameter, 10mm in wall thickness and 5360mm in length). The velocities of the incident and reflected compression waves in the foam were measured by using piezo-ceramic pressure sensors. The end-wall peak pressure behind the reflected wave in the foam was measured by using a crystal piezoelectric sensor. It is suggested that the high end-wall pressure may be caused by a rapid contact between the foam and the end-wall surface. Both open-cell and closed-cell foams with different length and density were tested. Through comparing the numerical and experimental end-wall pressure, the permeability coefficients a and 0 are quantitatively determined.

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

Gas exchange by bubbles in waves

A simple and feasible model feet the calculation of the gas transfer by bubble clouds is proposed in this article. N-2, O-2, and CO2 transferred by bubble clouds are obtained. At wind speed of 10 m/s, the calculated supersaturation of dissolved oxygen is 1.93-3.89% in agreement with the measurement.