共查询到20条相似文献，搜索用时 31 毫秒
1.
针对杭州湾独特的喇叭型强潮河口湾的特点,基于Blumberg等(1996)的ECOMSED模式,引入动边界技术,建立杭州湾三维动边界的潮流模型.模型以正交曲线坐标下三维非线性水动力方程为基本方程,应用Mellor和Yamada的2.5阶湍流闭合模型计算紊动黏滞系数,嵌入Grant和Madsen的底边界层模型考虑波浪对底部应力的作用,采用干湿网格法模拟潮流漫滩过程;综合考虑径流,风应力,密度流和M2,S2,K1,O1四个主要分潮和M4,S4,MS4三个浅水分潮的作用,从而提高杭州湾潮流模拟的精度.通过验潮站调和常数和多次海流连续观测资料的验证,表明该文建立的模型可以更好的用于杭州湾流场的预报模拟. 相似文献
2.
Constant flows, as well as oscillatory tidal flow, play an important role in the long-term dispersion of water in the Seto Inland Sea. Two kinds of numerical model (1-line and 2-line models) of the Seto Inland Sea have been developed to determine the role of density-induced currents, one type of the constant flow, in water dispersion in the Inland Sea. The seasonal variations of temperature, salinity and density fields are simulated and the density-induced current field is predicted at the same time. It is found that the most appropriate value of the longitudinal eddy diffusion coefficient,K
x, is 5×106–7×106 cm2sec–1. The value of the overall mean dispersion coefficient is of the order of 107cm2sec–1 (Hayami and Unoki, 1970). Consequently, it is suggested that 50–70% of the total dispersion in the Seto Inland Sea can be attributed to currents other than density-induced currents,i.e., tidal currents, tide-induced currents and wind-driven currents.In winter, both density and velocity fields, calculated using the 1-line model, satisfy the conditions for the existence of a coastal front in Kii Channel and in the eastern Iyo-nada. 相似文献
3.
The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water. 相似文献
4.
《Estuarine, Coastal and Shelf Science》2003,57(5-6):873-885
In the present work we explore the impact of assimilating local tide-gauge and altimetric data on the quality of predicting the major Adriatic tides (M2 and K1). To that end we compute optimal tidal open boundary conditions for a 3D high-resolution finite-element model by using an incremental assimilation formalism. The essence of the method is the use of two dynamical models where the solution in the complex 3D high-resolution model is sought via assimilation of prediction errors into the simpler 2D model with explicit inverse. In the central numerical experiment, harmonic constants from 12 tide gauges are assimilated and the results are analysed at 31 locations, hence 19 independent ones. The data assimilation contributes to the reduction of maximum amplitude error from 5.6 to 0.5 cm for M2 and from 3.9 to 0.1 cm for K1. The assimilation procedure is repeated by assimilating suitably processed Topex/Poseidon altimeter data, again validating the outcome at 31 tide gauge locations. The result was very similar to the gauge-data assimilation outcome. The model output is also validated with the current data, not used in the assimilation. At two locations and at three depths the model was able to reproduce the major and the minor semi-axes of tidal ellipses, as well as their orientations very well. 相似文献
5.
Hydrographic data show that the meridional deep current at 47°N is weak and southward in northeastern North Pacific; the strong
northward current expected for an upwelling in a flat-bottom ocean is absent. This may imply that the eastward-rising bottom
slope in the Northeast Pacific Basin contributes to the overturning circulation. After analysis of observational data, we
examine the bottom-slope effect using models in which deep water enters the lower deep layer, upwells to the upper deep layer,
and exits laterally. The analytical model is based on geostrophic hydrostatic balance, Sverdrup relation, and vertical advection–diffusion
balance of density, and incorporates a small bottom slope and an eastward-increasing upwelling. Due to the sloping bottom,
current in the lower deep layer intensifies bottomward, and the intensification is weaker for larger vertical eddy diffusivity
(K
V), weaker stratification, and smaller eastward increase in upwelling. Varying the value of K
V changes the vertical structure and direction of the current; the current is more barotropic and flows further eastward as
K
V increases. The eastward current is reproduced with the numerical model that incorporates the realistic bottom-slope gradient
and includes boundary currents. The interior current flows eastward primarily, runs up the bottom slope, and produces an upwelling.
The eastward current has a realistic volume transport that is similar to the net inflow, unlike the large northward current
for a flat bottom. The upwelling water in the upper deep layer flows southward and then westward in the southern region, although
it may partly upwell further into the intermediate layer. 相似文献
6.
Three-dimensional structure of tidal current in the East China Sea and the Yellow Sea 总被引：9，自引：0，他引：9
A three-dimensional tidal current model is developed and applied to the East China Sea (ECS), the Yellow Sea and the Bohai
Sea. The model well reproduces the major four tides, namely M2, S2, K1 and O1 tides, and their currents. The horizontal distributions of the major four tidal currents are the same as those calculated
by the horizontal two-dimensional models. With its high resolutions in the horizontal (12.5 km) and the vertical (20 layers),
the model is used to investigate the vertical distributions of tidal current. Four vertical eddy viscosity models are used
in the numerical experiments. As the tidal current becomes strong, its vertical shear becomes large and its vertical profile
becomes sensitive to the vertical eddy viscosity. As a conclusion, the HU (a) model (Davieset al., 1997), which relates the vertical eddy viscosity to the water depth and depth mean velocity, gives the closest results to
the observed data. The reproduction of the amphidromic point of M2 tide in Liaodong Bay is discussed and it is concluded that it depends on the bottom friction stress. The model reproduces
a unique vertical profile of tidal current in the Yellow Sea, which is also found in the observed data. The reason for the
reproduction of such a unique profile in the model is investigated. 相似文献
7.
AbstractIn this paper, two analytical models used worldwide to assess salinity variation in alluvial estuaries are applied to the Ashtamudi estuary, a Ramsar site, southwest coast of India; and Bouregreg estuary, in northwest part of Morocco. The estuaries’ bathymetry is described by an exponential function. Both models are quite similar and use a predictive equation for the dispersion in the estuary mouth (D0). The major difference between the two models is the use of the constant value of K?=?0.5 for the Van der Burgh coefficient (K) and the introduction of the correction factor ζ, which is a function of damping (δ) and shape (γ). The performance of these two models was evaluated by comparing their results with field measurements; this revealed that both analytical models apply well to both the estimation of salinity distribution and the prediction of salt intrusion in the Ashtamudi and Bouregreg estuaries (Ashtamudi: RMSE = 0.60–1.22 ppt; Bouregreg: RMSE = 0.92–2.71 ppt). One model agrees more with the field measurements of salinity distribution along the estuaries axis; the second underestimate and overestimate some values of salinity distribution along the estuaries. Possibly, the constant value of K?=?0.5 for the Van der Burgh coefficient (K) has applicability limits for the estuaries under tidal conditions. The specifying of the parameterization may be a field of research. 相似文献
8.
M. Zhou 《Estuarine, Coastal and Shelf Science》1998,46(6):811-825
Baroclinic circulation in highly stratified and partially stratified estuaries is characterised by a two-layer flow: a bottom salt- water inflow and a surface brackish-water outflow. Tidal period variation of the thicknesses of a two-layer flow is observed to be associated with mixing, bottom stress and hydraulic characteristics of superposed tidal and gravity currents. Here, both analytical two-layer hydraulic equations with weak friction and a numerical model including a turbulence closure were utilised to understand the mechanism of the layer tendency within a two-layer flow under different barotropic flow conditions. It has been found that in the weak bottom friction case, a gravity current has two critical solutions at the layer thickness equal to 0·5Hand 0·292H. The layer thickness towards a particular critical solution is dependent on the sign of the bottom stress, i.e. when the bottom stress is opposite (favor) to the bottom gravity current, its layer thickness converges to 0·5H(0·292H). In the case of strong bottom stress and mixing opposing the gravity current, the solutions of the gravity current layer thickness at 0·5Hand 0·292Hwill not be valid. Both mixing and vorticity produced by bottom stress erode the halocline, and produce a high velocity core in the mid-depth, which leads to the thickness of a bottom gravity current greater than 0·5H. These internal hydraulic tendency and mixing processes, varying with time-dependent barotropic tidal current forcing, determine the tidal period variation of the gravity current structure. 相似文献
9.
-In order to avoid prescribing open boundary condition on the upstream side of the Hangzhou Bay, in numerical simulation of the tides and residual currents of the Bay, a 1-D model for the Qiantang River is connected to the 2-D model for the Hangzhou Bay. The harmonic constants of diurnal constituent [ (K1 O1)/2],semidiurnal constituent (M2) and shallow water constituent (M4) are obtained. The results produced by the combined model are in better agreement with the observed ones than those produced solely by the original 2-D model. The combined model gives much more reliable results for tide-induced residual water level and current. 相似文献
10.
This study has been undertaken to quantify the tank wall effects on resistance estimation of ship models. Given the finite width of a tank, the flow around a ship model has been numerically modelled and the pressure and pressure related drag have been estimated.Since the model runs at speeds essentially in the laminar and transient speed range, an inviscid model has been chosen for obtaining the pressure drag component in the numerical studies. Grid dependency study has been done to optimize the mesh in the control volume for the numerical studies. An unstructured grid consisting of hexahedral cells has been used in the volume of fluid (VOF) model. The model chosen is a medium speed, ocean going barge and the residuary resistance has been obtained for different tank width conditions. The tank width has been defined using a non-dimensionalized parameter W/B (tank width W, model width B) ratio. The study shows that the residuary resistance obtained at W/B=5.0 is free from tank wall influence for the chosen model. The findings of the study have been compared by testing two geosim models under the same tank width conditions. The residuary resistance values have been compared with numerical results. The combined numerical experimental approach provides interesting results of consistency for comparison. The tank wall influences suggested by the numerical study are well quantified in the experimental study and give useful guideline for limiting wall influences. 相似文献
11.
In hydrodynamic models of marine and estuarine currents the use of boundary-fitted curvilinear grids not only makes the model grids fit to the coastline and bathymetry well, but also makes the kinetic boundary conditions simple and more accurate. Because of these advantages, a three-dimensional tidal model with boundary-fitted curvilinear grids has been developed to simulate both tide and current in estuarine and shelf water. The basic idea is to use a set of coupled σ-stretched and elliptic transformations to map the physical space into a corresponding transformed space such that all boundaries are coincident with co-ordinate lines and the transformed grids are rectangular. The hydrodynamic equations in the transformed space are solved in a rectangular mesh of the transformed grid system. The application of the model in the simulation of the M2,S2 , K1and O1tidal waves in the Bohai Sea, China, shows that the numerical results are in good agreement with the observations. The method of boundary-fitted curvilinear grids is effective for improving simulation accuracy of current in the estuarine and shallow seas, especially at the coastal regions where the current was usually impractical by a uniform Cartesian grid system with the shoreline and bathymetry represented by numerous stair-steps. 相似文献
12.
13.
Tang Yuxiang 《海洋学报(英文版)》1989,8(2):179-190
-A two-dimensional.nonlinear numerical model is used to study the residual current generated by tides in the East China Sea (ECS)and the South Huanghai Sea (SHS). At first, the principal semidiurnal lunar tide (M2)and the tidal current are derived in these areas. The results obtained with the model are strongly supported by the observational results available. Then, the tide-induced residual flow is determined by using the currents generated by the tidal input. The main features of the residual current in ECS and SHS are presented by analyzing the calculated results. Some of the problems are discussed such as the cause of generating residual current and the contribution of the residual current to the observed current. 相似文献
14.
15.
Wave–current flow is a phenomenon that is present in many practical engineering situations. Over the past several decades, this type of flow has been increasingly investigated under controlled laboratory conditions. This paper presents a numerical study of wave–current flow in the ocean basin of the LabOceano (COPPE/UFRJ). A homogeneous multiphase model based on the RANS equations and the k–ɛ turbulence model implemented in ANSYS-CFX code were used. A cross section of the ocean basin was represented. A regular wave with a height of 0.08 m and a period of 1.80 s (i.e., a wave steepness of H/L = 0.016), propagating on favourable currents, was simulated. The behaviour of the free surface elevation over time and the streamlines along the basin for wave and wave–current flows were presented. The numerical results were compared to the non-viscous theory given by the Rayleigh equation applied to the problem of wave–current interaction. Good agreement was found between the wave length estimated by the numerical results and the analytical solutions, with a deviation of less than 2%. 相似文献
16.
John B. Christoffersen 《Ocean Engineering》1985,12(5):387-423
Two simple two-layer eddy viscosity models, which facilitate analytical solutions, are presented in order to describe the velocity field and associated shear stress in a combined current wave motion. The models, which have the same eddy viscosity in the current boundary layer, but different eddy viscosities in the wave boundary layer, cover together the whole rough turbulent regime. Straightforward definitions are made for the wave friction factor and the current friction factor for the combined motion, which are in accordance with the results for pure waves and pure currents. In this way one avoids the fictitious reference velocities and elliptic integrals which e.g. Grant and Madsen (1978, 1979) experienced. The two friction factors turn out to be functions of four dimensionless parameters. A detailed calculation procedure is presented. Comparison with laboratory experiments yields promising results. A new relation connecting dissipation and bed shear stress is also developed. 相似文献
17.
The flow around a cylinder, surface mounted or partly buried has been examined at Reynolds number of 1.3×104 and 2.6×104 by using the flow visualization method in the case of steady current. The lengths of the upstream and downstream separation regions were determined and the characteristics of the flow area were discussed for the cylinder having different burial depth to the diameter ratio (G/D). For comparsion, some numerical results obtained from the application of the FLUENT, computational fluid dynamics (CFD) software package have also been presented. Both experimental and numerical results indicated that the lengths of the separation regions near the upstream and downstream of the cylinder decreased with the increasing burial ratio (G/D=0, 0.20, 0.50). The results have significant implications both for the illustration of the potential for using CFD software and highlighting the need for data of two dimensional flow around the buried cylinder. 相似文献
18.
Sensitivity analysis of a coastal marine ecosystem 总被引：1，自引：0，他引：1
An ecosystem with four compartments, i.e. PO4-P, phytoplankton, zooplankton and detritus, in a coastal region was considered. Sensitivity analysis of a phosphorus flow model in well-mixed water without currents concludes that (1) the maximal photosynthetic rate,V
m
, in Michaelis-Menten relation plays an important role in the distribution of biomass among compartments but the half saturation constant,K
s
, is not so important, (2) the natural death rate of phytoplankton is important for the ecosystem, (3) the natural death rate of zooplankton is also effective on the ecosystem.A numerical experiment was also performed on the ecosystem with four compartments are also studied using a dynamical barotropic model of tidal currents of Mikawa Bay (Japan). Diffusion coefficient in diffusion equation plays the role of a linear smoothing parameter in the horizontal distribution of compartment. On the other hand, perturbations of biological parameters cause nonlinear variations in the horizontal distribution of compartment. 相似文献
19.
Wilton Z. Arruda Doron Nof James J. O&#x;Brien 《Deep Sea Research Part I: Oceanographic Research Papers》2004,51(12):2073-2090
A nonlinear theory for the generation of the Ulleung Warm Eddy (UWE) is proposed. Using the nonlinear reduced gravity (shallow water) equations, it is shown analytically that the eddy is established in order to balance the northward momentum flux (i.e., the flow force) exerted by the separating western boundary current (WBC). In this scenario, the presence of β produces a southward (eddy) force balancing the northward momentum flux imparted by the separating East Korean Warm Current (EKWC).It is found that, for a high Rossby number EKWC (i.e., highly nonlinear current), the eddy radius is roughly 2Rd/ε1/6 (here ε≡βRd/f0, where Rd is the Rossby radius), implying that the UWE has a scale larger than that of most eddies (Rd). This solution suggests that, in contrast to the familiar idea attributing the formation of eddies to instabilities (i.e., the breakdown of a known steady solution), the UWE is an integral part of the steady stable solution. The solution also suggests that a weak WBC does not produce an eddy (due to the absence of nonlinearity).A reduced gravity numerical model is used to further analyze the relationship between β, nonlinearity and the eddy formation. First, we show that a high Rossby number WBC which is forced to separate from the wall on an f plane does not produce an eddy near the separation. To balance the northward momentum force imparted by the nonlinear boundary current, the f plane system moves constantly offshore, producing a southward Coriolis force. We then show that, as β is introduced to the problem, an anticyclonic eddy is formed. The numerical balance of forces shows that, as suggested by the analytical reasoning, the southward force produced by the eddy balances the northward flow force imparted by the boundary current. We also found that the observed eddy scale in the Japan/East Sea agrees with the analytical estimate for a nonlinear current. 相似文献
20.
Numerical study of current fields near the Changjiang Estuary and impact of Quick-EnKF assimilation 总被引：1，自引：0，他引：1
A 30-d current numerical simulation is running for the Yangshan Port, the Changjiang Estuary, the Hangzhou Bay and their adjacent seas using a finite volume coastal ocean model (FVCOM), with Changjiang River runoff and wind effect being considered. At the open boundary, this model is driven by the water level obtained from prediction including eight main partial tides. After the harmonic analysis, the cotidal chart and the iso-amplitude line as well as the current ellipse distribution map are displayed to illustrate the propagation property of a tidal wave. Horizontal velocity of both the U and V components coincides with the actual measurement, which shows that the model result is credible to describe the hydrodynamic pattern in this sea area. On this basis, real-time current data from high-frequency radar is assimilated with the implementation of quick ensemble Kalman filter, which takes the variation tendency of the state vector to compute the analysis field, instead of integrating the field for N (the number of ensemble) times as it used to in the standard EnKF, aiming at raising the efficiency of computation, reducing the error of prediction and at the same time, improving the forecast effect. 相似文献