Laboratory study of the dynamics of a layer of turbulence excited by standing surface waves

The evolution of a turbulent layer excited by standing surface parametric waves is studied experimentally. The time dependence of the depth of the turbulized region is obtained. It is found that the thikness of the layer is described by a nonlinear dependence as distinguished from the previously obtained theoretical results. The dynamics of turbulent-layer embedding is described on the basis of a semiempirical theory and the assumption that the outer scale of turbulence is constant. Satisfactory agreement with experimental data is obtained. In the context of the proposed model, a relationship is found between the quantities characterizing the turbulence and surface waves. The dependence of the steepness of surface waves on the acceleration of the bottom is obtained experimentally. The data obtained in this study can be used for estimation of turbulence generation during a submarine earthquake.     

The wind-field structure in a stably stratified atmospheric boundary layer over a rough surface

Both wind turning with height and ageostrophic flow in a stably stratified atmospheric boundary layer are analyzed using a three-parameter turbulence model. For a quasi-steady state of the boundary layer, the cross-isobaric flow is determined only by turbulent stress at the surface in the direction of geostrophic wind. The “operative” prediction models, in which the first-order turbulence closure schemes are used, tend to overestimate the boundary-layer depth and underestimate the angle between the surface and geostrophic winds when compared to “research” models (schemes of high-level turbulence closure). The true value of the angle between the surface and geostrophic winds is significant for the presentation of a large-scale flow. A nocturnal low-level jet is a mesoscale phenomenon reflected in data obtained from measurements in a stably stratified atmospheric boundary layer. It is found that such jets are of great importance in transporting humidity, momentum, and air pollution. In this study, the difference between jet flows over a homogeneous underlying surface and over a spatially localized large-scale aerodynamic roughness is shown.     

Turbulent mixing in stably stratified flows of the environment: The current state of the problem (Review)

Specific features of the turbulent transfer of the momentum and heat in stably stratified geophysical flows, as well as possibilities for including them into RANS turbulence models, are analyzed. The momentum (but not heat) transfer by internal gravity waves under conditions of strong stability is, for example, one such feature. Laboratory data and measurements in the atmosphere fix a clear dropping trend of the inverse turbulent Prandtl number with an increasing gradient Richardson number, which must be reproduced by turbulence models. Ignoring this feature can cause a false diffusion of heat under conditions of strong stability and lead, in particular, to noticeable errors in calculations of the temperature in the atmospheric boundary layer. Therefore, models of turbulent transfer must include the effect of the action of buoyancy and internal gravity waves on turbulent flows of the momentum. Such a strategy of modeling the stratified turbulence is presented in the review by a concrete RANS model and original results obtained during the modeling of stratified flows in the environment. Semiempirical turbulence models used for calculations of complex turbulent flows in deep stratified bodies of water are also analyzed. This part of the review is based on the data of investigations within the framework of the large international scientific Comparative Analysis and Rationalization of Second-Moment Turbulence Models (CARTUM) project and other publications of leading specialists. The most economical and effective approach associated with modified two-parameter turbulence models is a real alternative to classical variants of these models. A class of test problems and laboratory and full-scale experiments used by the participants of the CARTUM project for the approbation of numerical models are considered.     

Simultaneous acoustic observations of turbulence and zooplankton in the ocean

Models and laboratory experiments show that zooplankton may locate food more easily in turbulent waters, but whether plankton seek or avoid turbulence in the ocean is an open question. It is difficult to measure turbulence and plankton simultaneously and with the necessary spatial resolution using traditional methods (nets and airfoil shear sensors). Acoustics is commonly used to survey zooplankton abundance and recent studies have shown that stratified turbulence can also be a significant source of sound scatter. This may seem like more of a complication than a boon for those aiming to use acoustics to observe plankton in turbulence. We present acoustic data, however, that show that zooplankton and turbulence can be observed simultaneously with a single 307 kHz sounder. The different natures of the two targets (discrete targets versus a volume effect) allow them to be distinguished. The key is sampling the same targets at multiple ranges. The volume scattering strength of a discrete target will increase as the target nears the sounder, because the volume sampled decreases. Turbulence, as a volume scattering effect, has little range dependence to its volume scattering strength.     

A wave-resolving model for nearshore suspended sediment transport

This paper presents a wave-resolving sediment transport model, which is capable of simulating sediment suspension in the field-scale surf zone. The surf zone hydrodynamics is modeled by the non-hydrostatic model NHWAVE (Ma et al., 2012). The turbulent flow and suspended sediment are simulated in a coupled manner. Three effects of suspended sediment on turbulent flow field are considered: (1) baroclinic forcing effect; (2) turbulence damping effect and (3) bottom boundary layer effect. Through the validation with the laboratory measurements of suspended sediment under nonbreaking skewed waves and surfzone breaking waves, we demonstrate that the model can reasonably predict wave-averaged sediment profiles. The model is then utilized to simulate a rip current field experiment (RCEX) and nearshore suspended sediment transport. The offshore sediment transport by rip currents is captured by the model. The effects of suspended sediment on self-suspension are also investigated. The turbulence damping and bottom boundary layer effects are significant on sediment suspension. The suspended sediment creates a stably stratified water column, damping fluid turbulence and reducing turbulent diffusivity. The suspension of sediment also produces a stably stratified bottom boundary layer. Thus, the drag coefficient and bottom shear stress are reduced, causing less sediment pickup from the bottom. The cross-shore suspended sediment flux is analyzed as well. The mean Eulerian suspended sediment flux is shoreward outside the surf zone, while it is seaward in the surf zone.     

Experimental study of turbulent-wake expansion from a surface ship

We studied the time dynamics of a turbulent region excited by a moving surface vessel in a field experiment. The time dependences of the geometrical sizes of the turbulent region are obtained, and it is shown that the time dependence of the width of the turbulent wake at the initial stage is close to the power-law function with an exponent of 0.4 for different experiments (the depth of the wake remaining practically constant). In the semiempirical turbulence theory, we suggested a qualitative model describing the process of the initial expansion of a ship wake as a diffusion of a one-dimensional layer of turbulized liquid due to a pulsed source.     

Observational and numerical modeling methods for quantifying coastal ocean turbulence and mixing

In this review paper, state-of-the-art observational and numerical modeling methods for small scale turbulence and mixing with applications to coastal oceans are presented in one context. Unresolved dynamics and remaining problems of field observations and numerical simulations are reviewed on the basis of the approach that modern process-oriented studies should be based on both observations and models. First of all, the basic dynamics of surface and bottom boundary layers as well as intermediate stratified regimes including the interaction of turbulence and internal waves are briefly discussed. Then, an overview is given on just established or recently emerging mechanical, acoustic and optical observational techniques. Microstructure shear probes although developed already in the 1970s have only recently become reliable commercial products. Specifically under surface waves turbulence measurements are difficult due to the necessary decomposition of waves and turbulence. The methods to apply Acoustic Doppler Current Profilers (ADCPs) for estimations of Reynolds stresses, turbulence kinetic energy and dissipation rates are under further development. Finally, applications of well-established turbulence resolving particle image velocimetry (PIV) to the dynamics of the bottom boundary layer are presented. As counterpart to the field methods the state-of-the-art in numerical modeling in coastal seas is presented. This includes the application of the Large Eddy Simulation (LES) method to shallow water Langmuir Circulation (LC) and to stratified flow over a topographic obstacle. Furthermore, statistical turbulence closure methods as well as empirical turbulence parameterizations and their applicability to coastal ocean turbulence and mixing are discussed. Specific problems related to the combined wave-current bottom boundary layer are discussed. Finally, two coastal modeling sensitivity studies are presented as applications, a two-dimensional study of upwelling and downwelling and a three-dimensional study for a marginal sea scenario (Baltic Sea). It is concluded that the discussed methods need further refinements specifically to account for the complex dynamics associated with the presence of surface and internal waves.     

Classification of highly turbid Jiaojiang Estuary

ThisstudywassupportedbytheNationalNaturalScienceFoundationofChinaundercontractNo.49276274,theZhejiangProvinceNaturalScienceFoundationundercontractNo.490013,theChina-Australiabilateralscienceandtechnologyprogram,theAustralianInstituteofMarineScience,theModellingLaboratoryoftheMarineScienceintheSecondInstituteofOceanographyoftheStateOceanicAdministration.INTRODUCTIONTheJiaojiangEstuaryis1ocatedintheeasterncoastofChina,2OokmfromthesouthoftheChangjiangRiver(YangtzeRiver),linkedin…     

Downwind sail aerodynamics: A CFD investigation with high grid resolution

A computational fluid dynamics (CFD) code was applied to an America's Cup Class Yacht to investigate sailing performance in a downwind configuration. Apparent wind angles at 45°, 105° and 120° are reported, sailed with mainsail and asymmetrical spinnakers. Numerical results are in good agreement with wind tunnel data. A large mesh investigation was performed, ranging from 60,000 elements up to 37 million elements, which shows a converging trend to the experimental values with differences smaller than 3% in both lift and drag. The most commonly used turbulence models in sail applications were tested and the results are presented here in two meshes with 1 million elements and 6.5 millions, respectively. All turbulence models over-estimate forces more than solving the Navier–Stokes system without any additional equations, hence turbulence models do not increase solution accuracy according to these results.     

On the application of a turbulence closure modified model to the description of the density jump evolution in a stably stratified medium

The self-similar turbulent density jump evolution has been studied in the scope of a turbulence closure modernized theory which takes into account the anisotropy and mutual transformation of the turbulent fluctuation kinetic and potential energy for a stably stratified fluid. The numerical calculation, performed using the equations for the average density and kinetic and potential energies of turbulent fluctuations, indicates that the vertical profiles of the buoyancy frequency, turbulence scale, and kinetic and potential energies drastically change when the turbulence anisotropy is strong. The vertical profiles of the corresponding energy and spatial discontinuity parameters, calculated at a weaker anisotropy, indicate that similar drastic changes are absent and a qualitative agreement exists with the known analytical solution, which describes the density jump evolution in a freshwater basin and was obtained previously [5, 8] in the scope of a turbulence local-similarity hypothesis applied in combination with the budget equation for the turbulent fluctuation kinetic energy.     

Three-parameter model of turbulence for the atmospheric boundary layer over an urbanized surface

A modified three-parameter model of turbulence for a thermally stratified atmospheric boundary layer (ABL) is presented. The model is based on tensor-invariant parametrizations for the pressure-strain and pressure-temperature correlations that are more complete than the parametrizations used in the Mellor-Yamada model of level 3.0. The turbulent momentum and heat fluxes are calculated with explicit algebraic models obtained with the aid of symbol algebra from the transport equations for momentum and heat fluxes in the approximation of weakly equilibrium turbulence. The turbulent transport of heat and momentum fluxes is assumed to be negligibly small in this approximation. The three-parameter E ? ε ? 2> model of thermally stratified turbulence is employed to obtain closed-form algebraic expressions for the fluxes. A computational test of a 24-h ABL evolution is implemented for an idealized two-dimensional region. Comparison of the computed results with the available observational data and other numerical models shows that the proposed model is able to reproduce both the most important structural features of the turbulence in an urban canopy layer near the urbanized ABL surface and the effect of urban roughness on a global structure of the fields of wind and temperature over a city. The results of the computational test for the new model indicate that the motion of air in the urban canopy layer is strongly influenced by mechanical factors (buildings) and thermal stratification.     

A numerical study of swash flows generated by bores

The dynamic processes of bore propagation over a uniform slope are studied numerically using a 2-D Reynolds Averaged Navier–Stokes (RANS) solver, coupled to a non-linear k − ε turbulence closure and a volume of fluid (VOF) method. The dam-break mechanism is used to generate bores in a constant depth region. Present numerical results for the ensemble-averaged flow field are compared with existing experimental data as well as theoretical and numerical results based on non-linear shallow water (NSW) equations. Reasonable agreement between the present numerical solutions and experimental data is observed. Using the numerical results, small-scale bore behaviors and flow features, such as the bore collapse process near the still-water shoreline, the ‘mini-collapse’ during the runup phase and the ‘back-wash bore’ in the down-rush phase, are described. In the case of a strong bore, the evolution of the averaged turbulence kinetic energy (TKE) over the swash zone consists of two phases: in the region near the still-water shoreline, the production and the dissipation of TKE are roughly in balance; in the region farther landwards of the still-water shoreline, the TKE decay rate is very close to that of homogeneous grid turbulence. On the other hand, in the case of a weak bore, the bore collapse generated turbulence is confined near the bottom boundary layer and the TKE decays at a much slower rate.     

Large Eddy Simulation for Plunge Breaker and Sediment Suspension

Breaking waves are a powerful agent for generating turbulence that plays an important role in many fluid dynamical processes, particularly in the mixing of materials. Breaking waves can dislodge sediment and throw it into suspension, which will then be carried by wave-induced steady current and tidal flow. In order to investigate sediment suspension by breaking waves, a numerical model based on large-eddy-simulation (LES) is developed. This numerical model can be used to simulate wave breaking and sediment suspension. The model consists of a free-surface model using the surface marker method combined with a two-dimensional model that solves the flow equations. The turbulence and the turbulent diffusion are described by a large-eddy-simulation (LES) method where the large turbulence features are simulated by solving the flow equations, and a subgrid model represents the small-scale turbulence that is not resolved by the flow model. A dynamic eddy viscosity subgrid scale stress model has been used for the     

Mixing processes relevant to phytoplankton dynamics in lakes

Active turbulence in lakes is confined to the surface mixed layer, to boundary layers on the lake sides and bottom, and to turbulent patches in the interior. The density stratification present in most lakes fundamentally alters the pathways connecting external mechanical energy inputs, for example by wind, with its ultimate fate as dissipation to heat; the density stratification supports internal waves and intrusions that distribute the input energy throughout the lake. Intrusions may be viewed as internal waves with zero horizontal wavenumber and are formed each time localised mixing occurs in a stratified fluid. Intrusions are also formed in the epilimnion by differential heating or cooling and by differential deepening. The fraction of lake volume below the diurnal mixed layer that is subject to active turbulence is very small, probably of the order of 1% or less in small to medium‐sized lakes. By contrast, in the surface mixed layer, turbulence is less intermittent and maintains phytoplankton in suspension and controls their exposure to the underwater solar flux. Nutrient transport to individual cells depends not only on the cell Reynolds number but also on the Peclet number, which, if large, implies enhanced mass transfer above purely diffusive values.     

End wall effects in experiments on water intrusion along the interface between two homogeneous layers

End wall effect puts an inherent limitation on tank experiments, especially when the problems in a stratified fluid are dealt with. During experiments involving horizontal intrusion along the interface between two homogeneous layers, a curious phenomenon was found, i.e., the tip of the intruding water wedge continues to extend for a short time after supply is stopped, but then it begins to retreat in the cases of relatively high Reynolds numbers. The cause of this retreat of the wedge was investigated and was shown to be attributable to the initial disturbance generated near the mouth of the feeder at the start of water supply which propagates along the interface layer and reflects at the end of the tank as a bulge of the interface layer. The retreat of the intruding wedge would not occur in a sufficiently long tank, and so the cause of the retreat can be considered as one kind of end wall effect in a tank.     

Evolution and breaking of a propagating internal wave in stratified ocean

The evolution and breaking of a propagating internal wave are directly numerically simulated using a pseudo-spectral method. The mechanism of PSI （ parametric subharmonic instability） involved in the evolution is testified clearly. It dominates gradually in nonlinear resonant interactions. As a consequence, the energy cascades to a second plant wave packet which has lower frequencies and higher wavenumbers than that of the primary wave. With the growth of this wave packet, wave breaking occurs and causes strongly nonlinear regime, i.e. stratified turbulence. The strong mixing and intermittent of the turbulence can be learned from the evolution of the total energy and kurtosis of vorticity vs. time. Some statistic properties of the stratified turbulence are also analyzed, including the spectra of KE （kinetic energy） and PE （potential energy）. The results show that the PE spectra display a wavenumber range scaling as 0. 2 N^4ky^-3 （N is the Brunt - Vaisala frequency, k, is the vertical wavenumber）, which is called buoyancy sub-range. However, the KE spectra cannot satisfy the negative cubic law of vertical wavenumber, which have a much larger downtrend than that of the PE spectra, for the potential energy is transferred more efficiently toward small scales than the kinetic energy. The Cox number of diapycnal diffusivity is also calculated, and it shows a good consistency with the observations and deductions in the ocean interior, during the stage of the stratified turbulence maintaining a fairly active level.     

珠江河口底边界层憩流动力特征

通过在珠江黄茅海河口进行的25 h定点连续观测,分析了潮流底边界层在憩流时刻的动力特征。结果表明,憩流期非恒定性增强,改变了平均流和湍流结构。主要表现:1)平均流结构出现三种不同类型,分别是对数分布(摩擦力主导)、垂向均匀分布(惯性主导)和过渡状态(摩擦力和惯性力平衡);2)湍流各向异性增强,湍谱惯性子区带宽减小。各湍流特征参数在转流时刻达到最小值,表明惯性作用抑制了湍流的发育。     

Numerical Experiment on Two-Dimensional Line Thermal

The time evolution of a two-dimensional line thermal-a turbulent flow produced by an initial element with signifi-cant buoyancy released in a large water body, is numerically studied with the two-equation k - e model for turbulence closure. The numerical results show that the thermal is characterized by a vortex pair flow and a kidney shaped concentra-tion structure with double peak maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time around 10. There are two regions in the flow field of a line thermal: a mixing region where the concentration of tracer fluid is high and the flow is turbulent and rotational with a pair of vortex eyes, and an ambient region where the concentration is zero and the flow is potential and well-described by a model of doublet with strength very close to those given by early experimental and analytical studies. The added virtual mass coeffi-cient of the thermal motion is found to be approximat     

On the turbulent prandtl number in a stably stratified atmospheric boundary layer

The results obtained from both atmospheric and laboratory measurements and from LES data show that, in the stably stratified flows of the atmospheric boundary layer, turbulent mixing occurs at gradient Richardson numbers Ri _g that significantly exceed one: the inverse turbulent Prandtl number Pr _t ⁻¹ decreases with an increase in the thermal flow stability. The decreasing trend of the inverse turbulent Ptandtl number is reproduced in a stably stratified atmospheric boundary layer in agreement with measurement data with the aid of an improved three-parameter turbulence model. In this model, a modified model that takes into account the effect of stratification in the expression for the time scale of the scalar field is used for the pressure-scalar correlation.     

Nonlinear modeling of liquid sloshing in a moving rectangular tank

A nonlinear liquid sloshing inside a partially filled rectangular tank has been investigated. The fluid is assumed to be homogeneous, isotropic, viscous, Newtonian and exhibit only limited compressibility. The tank is forced to move harmonically along a vertical curve with rolling motion to simulate the actual tank excitation. The volume of fluid technique is used to track the free surface. The model solves the complete Navier–Stokes equations in primitive variables by use of the finite difference approximations. At each time step, a donor–acceptor method is used to transport the volume of fluid function and hence the locations of the free surface. In order to assess the accuracy of the method used, computations are verified through convergence tests and compared with the theoretical solutions and experimental results.