首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 234 毫秒
1.
The influence of high vertical velocity gradients in the Black Sea Rim Current on the intensity of the vertical turbulent exchange is demonstrated on the basis of numerical modeling based on CTD data. The vertical turbulent exchange is confirmed by the anomalous distribution of the hydrochemical parameters in the redox layer. A system of equations for the kinetic energy of the turbulence and dissipation rate (k-? model) is used for the calculation of the coefficient of the vertical turbulent viscosity (diffusivity).  相似文献   

2.
The results of simulating the hydrophysical fields of the Black Sea with a resolution of 1.64 × 1.64 km for January–September 2006 with the use of real atmospheric forcing are analyzed. Both vertical turbulent momentum exchange and vertical turbulent heat and salt diffusions are parameterized using the Mellor-Yamada level 2.5 scheme. The results of this numerical experiment are compared with similar data obtained with a horizontal resolution of 5 km. The features of the meso- and submesoscale dynamics of waters in individual sea regions are given. Possible physical mechanisms of forming meso- and submesoscale vortices are studied on the basis of energy analysis. It is shown that, in the absence of significant wind forcing, the main contribution to kinetic energy is made by the buoyancy force and wind-field inhomogeneities result in significant variations in both total vertical viscosity and total vertical diffusion.  相似文献   

3.
The frequency spectrum of surface elevations in the presence of wind waves is well known. On this basis, one can estimate the frequency spectrum of vertical velocities in sea-surface waves. Owing to liquid incompressibility, the spectrum of horizontal velocities should have the same frequency dependence. The use of the dispersion equation for waves on the surface of a heavy liquid allows one to obtain to the spatial spectrum of velocities. Therefore, one can estimate the spatial structure function of the velocity field. For short waves and large depths, the structure function increases as r 1/2, where r is the distance between the points of observations. For long waves and shallow depths h, this increase is proportional to r. The coefficient of turbulent mixing K(r) of pollution spots of size r on the sea surface is now estimated as the product of the spot size and the rms difference of velocities. As a result, depending on r and h, the exponent in the r n dependence of K(r) may vary between 1.25 and 1.5. This outcome provides an explanation for a scatter in the values of the exponent n, a phenomenon that has been observed by many experimentalists.  相似文献   

4.
Estimates of vertical turbulent diffusion coefficient (K t ) in the Black Sea pycnohalocline have been obtained from data of simultaneous observations of seawater temperature, salinity, density, and horizontal current velocity, obtained in the northeastern part of the Black Sea during 2013–2014 with a moored Aqualog profiler. A Munk and Andersson (1948) type parameterization, adapted for the Black Sea environment, is proposed for calculating K t . Strong short-period (several days) variability of turbulent exchange is revealed, induced by vertical shear variations of the current velocity.  相似文献   

5.
Two numerical studies (Endoh, 1977;Harashima et al., 1978) have been proposed on a front formed by a coupling effect of cooling of the sea surface and inflow of the fresh water in a vertical two-dimensional plane without the rotation of the earth. It is, however, not easy to interpret their numerical results. A simple interpretation will be proposed by an analytical study in this paper.It is found that local convection due to the density inversion, which is expressed by the convective adjustment of the vertical diffusion coefficient in the actual numerical calculations, plays an important role on the front formation.The characteristics of the front is also clarified in the case of steady state. Namely, simple functional dependences are obtained of the position and the width of the front, the horizontal and the vertical velocities and the distribution of the buoyancy and the salinity in the neighborhood of the front on the horizontal coordinate, the cooling rate, the eddy coefficients of diffusion and viscosity, the water depth and the vertically averaged horizontal fluxes of buoyancy and salinity.  相似文献   

6.
Based on the fluid motion equations, the physical meaning of eddy viscosity coefficient and the rationality of the Boussinesq hypothesis are discussed in this paper. The effect of the coefficient on numerical stability is analyzed briefly. A semi-enclosed rectangular sea area, with an orthogonal spur dike, is applied in a 2-D numerical model to study the effect of horizontal eddy viscosity coefficient (AH). The computed result shows that AH has little influence on the tidal level and averaged flow velocity, but has obvious influence on the intensity and the range of return flow around near the spur dike. Correspondingly, a wind-driven current pool and an annular current are applied in a 3-D numerical model respectively to study the effect of vertical eddy viscosity coefficient (Av). The computed result shows that the absolute value of Av is inversely proportional to that of horizontal velocity, and the vertical gradient value of Av determines the vertical distribution of horizontal velocity. The distrib  相似文献   

7.
The characteristics of wave and turbulence velocities created by a broad-banded irregular wave train breaking on a 1:35 slope were studied in a laboratory wave flume. Water particle velocities were measured simultaneously with wave elevations at three cross-shore locations inside the surf zone. The measured data were separated into low-frequency and high-frequency time series using a Fourier filter. The measured velocities were further separated into organized wave-induced velocities and turbulent velocity fluctuations by ensemble averaging. The broad-banded irregular waves created a wide surf zone that was dominated by spilling type breakers. A wave-by-wave analysis was carried out to obtain the probability distributions of individual wave heights, wave periods, peak wave velocities, and wave-averaged turbulent kinetic energies and Reynolds stresses. The results showed that there was a consistent increase in the kurtosis of the vertical velocity distribution from the surface to the bottom. The abnormally large downward velocities were produced by plunging breakers that occurred from time to time. It was found that the mean of the highest one-third wave-averaged turbulent kinetic energy values in the irregular waves was about the same as the time-averaged turbulent kinetic energy in a regular wave with similar deep-water wave height to wavelength ratio. It was also found that the correlation coefficient of the Reynolds stress varied strongly with turbulence intensity. Good correlation between u′ and w′ was obtained when the turbulence intensity was high; the correlation coefficient was about 0.3–0.5. The Reynolds stress correlation coefficient decreased over a wave cycle, and with distance from the water surface. Under the irregular breaking waves, turbulent kinetic energy was transported downward and landward by turbulent velocity fluctuations and wave velocities, and upward and seaward by the undertow. The undertow in the irregular waves was similar in vertical structure but lower in magnitude than in regular waves, and the horizontal velocity profiles under the low-frequency waves were approximately uniform.  相似文献   

8.
The summertime phytoplankton bloom near the Kerguelen Plateau is in marked contrast to the low-chlorophyll conditions typical of the Southern Ocean and is thought to arise from natural iron fertilisation. The mechanisms of iron supply to the euphotic zone in this region are poorly understood, and numerical studies of iron transport have until now omitted fine-scale (sub-mesoscale) dynamics which have been shown to significantly increase vertical transport in other parts of the ocean.We present the first sub-mesoscale-resolving study of the flow and vertical transport in this region. The modelled transport and flow structure agree well with observations. We find that an increase in horizontal resolution from mesoscale-resolving (1/20°) to 1/80° resolves sub-mesoscale filamentary frontal structures in which vertical velocities are dramatically higher and are consistent with available observations. Lagrangian tracking shows that water is advected to the surface from much greater depth in the sub-mesoscale-resolving experiment, and that vertical exchange is far more rapid and frequent. This study of sub-mesoscale vertical velocities sets the foundation for subsequent investigation of iron transport in this environment.  相似文献   

9.
The numerical analysis of the stationary field of current velocity on the upper boundary of the bottom boundary layer in the Barents Sea is performed on the basis of a simplified model taking into account the fields of wind velocity and density of water for the principal periods of the seasonal cycle and the bottom topography. The analysis is based on the climatic BarKode database and the data on the wind velocity over the Barents Sea for the last 50 yr. The numerical results demonstrate that the field of bottom currents is fairly nonuniform and the current velocities vary from several fractions of 1 cm/sec to 5 cm/sec in the zones with noticeable slopes of the bottom. The estimates of the thickness of the bottom boundary layer are obtained for the constant coefficient of bottom friction C f = 0.04. In the major part of the water area of the Barents Sea, the thickness of the bottom boundary layer is close to 1 m. In the regions with significant slopes of the bottom, it increases to 2–2.5 m and, in the two zones of intensification of the bottom currents, becomes as large as 5 m. The maximum estimate of the coefficient of turbulent viscosity is close to 5 cm2/sec. The mean value of the coefficient of vertical density diffusion K S is equal to 2.34 cm2/sec and its standard deviation is equal to 1.52 cm2/sec. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 4, pp. 31–49, September–October, 2007.  相似文献   

10.
The k-ε turbulence model which relates the eddy viscosity to turbulent kinetic energy, k, and to the rate of its dissipation, ε, and determines the distribution of these two quantities from modeled transport equations, is employed in calculating the vertical structure of wave-induced turbulent flows in two-dimensional estuaries. The empirical constants in this model are given the standard values cited in the literature and used successfully for calculating a large variety of steady flows; they are not tuned to the oscillating flows. The free surface elevation and the eddy-viscosity distribution are also calculated as a part of the solution. First, calculations performed for oscillatory laminar flows are compared with analytic solutions to ensure the proper performance of the numerical scheme. The turbulent flow in a laboratory-estuary model with one end closed, and in the Humber Estuary is simulated with the numerical model, and the results are compared with the corresponding flume and field measurements. The influence of frequency and roughness is thereby investigated. The results indicate that the time and space variation of the velocity field and the phase lag between the surface slope and the horizontal velocity can be predicted satisfactorily in wave-induced turbulent flows.  相似文献   

11.
A complication of finite-volume triangular C-grid methods is the numerical emergence of horizontal divergence errors that lead to grid-scale oscillations in vertical velocity. Nonlinear feedback via advection of momentum can lead to numerical instability in velocity modes via positive feedback with spurious vertical velocities induced by horizontal divergence truncation error. Existing strategies to mitigate divergence errors such as direct divergence averaging and increased diffusion do not completely mitigate horizontal vertical velocity oscillations. We present a novel elliptic filtering approach to mitigate this spurious error and more accurately represent vertical velocities via improved calculation of horizontal divergences. These results are applied to laminar curved channel flows, demonstrating the applicability of the method to reproduce secondary flow features.  相似文献   

12.
We propose an improvement of the algorithm of joint assimilation of the data on climatic temperature, salinity, and altimetric sea level in a model of circulation. Unlike the previous works, the variances of the forecast errors of temperature and salinity and the cross-covariance functions of of the forecast errors of salinity-level and temperature-level depend on the dynamics of waters. It is shown that the structure of the fields of cross-covariance functions in the upper mixed layer is formed by the vertical turbulent diffusion of the variances of forecast errors of temperature and salinity. At greater depths, these statistical characteristics are mainly determined by the vertical advection. We compared the results of calculations with and without taking into account the dynamics of the statistical characteristics. The analysis of the influence of the dynamics of these characteristics makes it possible to reconstruct the mutually adapted climatic fields of temperature, salinity, and horizontal and vertical current velocities in the Black Sea with the assimilation of data in the numerical model in each time step. Translated from Morskoi Gidrofizicheskii Zhurnal, No. 4, pp. 18–31, July–August, 2008.  相似文献   

13.
The long-term evolution of the Black Sea dynamics (1980–2020) is reconstructed by numerical simulation. The model of the Black Sea circulation has 4.8 km horizontal spatial resolution and 40 levels in z-coordinates. The mixing processes in the upper layer are parameterized by Mellor-Yamada turbulent model. For the sea surface boundary conditions, atmospheric forcing functions were used, provided for the Black Sea region by the Euro mediterranean Center on Climate Change (CMCC) from the COSMO-CLM regional climate model. These data have a spatial resolution of 14 km and a daily temporal resolution. To evaluate the quality of the hydrodynamic fields derived from the simulation, they were compared with in-situ hydrological measurements and similar results from physical reanalysis of the Black Sea.  相似文献   

14.
A fundamental mechanism of generation of the tidal residual flow, the steady or quasi-steady flow induced in the tidal current system, is studied by numerical methods. The model basin is a very simple one, a rectangular basin of 5m×10m of constant depth and with a cape of 4 m length jutting out at a right angle from the center of the longer side wall. This basin has the same topography as that studied byYanagi (1976) by means of the hydraulic model experiments.The steady, circular, horizontal current is found to be induced through the following processes. Horizontal friction at the coast makes the vorticity of vertical component in the oscillating flow. Self-interaction of this flow causes the vorticity transfer to the steady flow in frequency domain. This vorticity transfer is confined in the narrow coastal boundary layer. The steady flow advects the transferred vorticity and makes itself develop fully wide over the bay. In other words, there are two kinds of cascade-up, one with regard to time scale and the other with regard to horizontal space scale.When the tidal range, the tidal period and the horizontal eddy viscosity change under the condition that the model geometry is fixed, the nondimensional parameter which controlls the steady flow is found to be the Reynolds number of the oscillating flow.  相似文献   

15.
A coastal ocean -coordinate model of Monterey Bay (MOB) with realistic bottom topography and coastlines is developed using the Princeton Ocean Model (POM) and grid generation technique (GGT) to study the horizontal pressure gradient errors associated with the MOB steep topography. The submarine canyon in MOB features some of the steepest topography encountered anywhere in the world oceans. The MOB grids are designed using the EAGEAL View and GENIE++ grid generation systems. A grid package developed by Ly and Luong (1993) is used in this study to couple grids to the model. The MOB model is tested with both orthogonal and curvilinear nearly-orthogonal (CNO) grids. The CNO grid has horizontal resolution which varies from 300 m to 2 km, while the resolution of the orthogonal grid is uniform with x = 1.25 km and y = 1.38 km. These grids cover a domain of 180 × 160 km with the same number of grid points of 131 × 131. Vertical resolutions of 25, 35 and 45 vertical sigma levels are tested. The error in the MOB are evaluated in terms of mean kinetic energy and velocity against various grids, vertical, horizontal resolution and distributions, and bottom topography smoothing. Simulations with various grids show that GGT can be used as another tool in reducing -coordinate errors in coastal ocean modeling besides increasing resolution and smoothing bottom topography. Topographical smoothing not only reduces topographic slope, but changes realistic topography. A CNO grid with a high grid density packed along steep slopes and Monterey Submarine Canyon reduces the errors by 40% compared to a rectangular grid with the same number of grid points. The CNO grid is more efficient than the rectangular grid, since it has most of its grids over water. The simulations show that the presented MOB -coordinate model can be used with a confidence regarding horizontal pressure gradient error.  相似文献   

16.
Studies of mesoscaleT, S characteristics of deep waters in the Black Sea carried out during 1978–1988 are generalized. Transformed lower Bosphorus waters are spread in the layer from 400 m to the bottom following density isosurfaces corresponding to this water. This is the way that interlayers ofT, S inversions are formed; their vertical thickness reaches 100–200 m in the Anatolian flow. The streamwise decrease of theT, S maxima and the location of isotherms in the cross-sections of these tongues permit the estimation of the order of the values of the coefficients of the vertical and horizontal turbulent diffusion.Translated by Mikhail M. Trufanov.  相似文献   

17.
Based on the Boussinesq assumption,derived are couple equations of free surface elevationand horizontal velocities for horizontal irrotational flow,and analytical expressions of the correspondingpressure and vertical velocity.After the free surface elevation and horizontal velocity at a certain depth areobtained by numerical method,the pressure and vertical velocity distributions can be obtained by simplecalculation.The dispersion at different depths is the same at the O(ε)approximation.The waveamplitude will decrease with increasing time due to viscosity,but it will increase due to the matching ofviscosity and the bed slope.thus,flow is unstable.Numerical or analytical results show that the waveamplitude.velocity and length will increase as the current increases along the wave direction.but theamplitude will increase.and the wave velocity and length will decrease as the water depth decreases.  相似文献   

18.
In this study,characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical twodimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propa-gates over breakwater. When wave crest propagates over breakwater,the anticlockwise vortex may generate. On the contrary,when wave hollow propagates over breakwater,the clockwise vortex may generate. Meanwhile,the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy,turbulent dissi-pation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.  相似文献   

19.
Q.Y. Zhang   《Ocean Engineering》2006,33(2):137-151
A comparison of two three-dimensional numerical modeling systems for tidal elevations and velocities in the coastal waters is presented. The two modeling systems are: (1) the Princeton Ocean Model (POM) and (2) the MIKE 3 flow model. The model performance results for Singapore's coastal waters show that the predicted tidal elevations from the two hydrodynamic modeling systems are almost identical and are in very good agreement with field measurement data. The simulated tidal current velocities match well with field measurement data at the selected stations, but it seems that the POM provides the slightly better simulation, compared to the MIKE 3 flow model. The depth profiles of the velocities obtained from the two modeling systems may be greatly different at some time, due to the vertical diffusion coefficient calculated from different turbulent sub-models in the two modeling systems. The POM generally predicts larger peak tidal velocities. The maximum speed differences for the model results from the two modeling systems occur in the top and differ from time to time and from location to location, reaching up to 20%.  相似文献   

20.
A three-dimensional,first order turbulence closure,thermal diffusion model is described inthis paper.The governing equations consist of an equation of continuity,three components of momentum,conservation equations for salt,temperature and subgridscale energy,and an equation of state.In the mod-el,according to the hypothesis of Kolmogorov and Prandtl,the viscosity coefficient of turbulent flow ofhomogeneous fluid is related to the local turbulent energy,and the horizontal and vertical exchangecoefficients of mass,heat and momentum are computed with the introduction of subgridscale turbulenceenergy.The governing equations are solved by finite difference techniques.This model is applied to theJiaozhou bay to predict thermal pollution by the Huangdao power plant.An instantaneous tidal currentfield is computed,then the distribution of temperature increment is predicted,and finally the effect of windstress on thermal discharge is discussed.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号