首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Fine grids with small spacing in boundary-fitted coordinates are normally used to treat thecomputation of fluid dynamics for estuarine areas and tidal flats.However,the adoption of Cartesian com-ponents of velocity vectors in this kind of non-orthogonal coordinates will definitely result in a diffculty insolving implicitly the transformed momentum equations,and also complicate the wet-dry point judgementused for flood areas.To solve this problem,equations in terms of generalized contravariant velocityvectors in curvilinear coordinates are derived in the present study,by virtue of which,anAlternative-Direction-Implicit numerical scheme in non-orthogonal grids would then be easily obtained,and wet-dry point judgement would as well be largely simplified.A comparison is made between the explic-it scheme and implicit scheme,showing that the present model is accurate and numerically stable for com-putations of fluid dynamics for estuarine areas and tidal flats.  相似文献   

2.
ABFGmodelforcalculationoftidalcurrentanddiffusionofpollutantsinnearshoreareas¥ShiFengyan;andZheng;Lianyuan(StateKeyLaboratory...  相似文献   

3.
Unstructured-grid models grounded on semi-implicit, finite-volume, Eulerian–Lagrangian algorithms, such as UnTRIM and ELCIRC, have enjoyed considerable success recently in simulating 3D estuarine and coastal circulation. However, opportunities for improving the accuracy of this type of models were identified during extensive simulations of a tightly coupled estuary–plume–shelf system in the Columbia River system. Efforts to improve numerical accuracy resulted in SELFE, a new finite-element model for cross-scale ocean modeling. SELFE retains key benefits, including computational efficiency of existing semi-implicit Eulerian–Lagrangian finite-volume models, but relaxes restrictions on grids, uses higher-order shape functions for elevation, and enables superior flexibility in representing the bathymetry. Better representation of the bathymetry is enabled by a novel, “localized” vertical grid that resembles unstructured grids. At a particular horizontal location, SELFE uses either S coordinates or SZ coordinates, but the equations are consistently solved in Z space. SELFE also performs well relative to volume conservation and spurious oscillations, two problems that plague some finite-element models. This paper introduces SELFE as an open-source code available for community use and enhancement. The main focus here is on describing the formulation of the model and on showing results for a range of progressively demanding benchmark tests. While leaving details to separate publications, we also briefly illustrate the superior performance of SELFE over ELCIRC in a field application to the Columbia River estuary and plume.  相似文献   

4.
A three-dimensional semi-implicit finite volume numerical model has been developed and applied to study tidal circulation and salinity stratification in the region of Oujiang River Estuary, China. The model employs horizontally unstructured grids and boundary-fitted coordinate system in the vertical direction. Governing equations consisting of continuity, momentum, and transport equations are all solved in the integral form of the equations, which provides a better representation of the conservative laws for mass, momentum, and transport in the coastal region with complex geometry and bottom bathymetry. The model performance was firstly quantified with skill assessment statistics on the choice of different parameters and validated with observed tidal elevation, current velocity, direction and salinity data over a spring–neap tidal cycle collected in 2006. Numerical results show that the model with wetting–drying capability successfully simulated the tidal currents and salinity fields with a reasonable accuracy and indicate that the Oujiang River Estuary is a macrotidal estuary with strong tidal mixing. In addition, the model results also show that the Oujiang River Estuary is a well-mixed estuary during spring tide. Then, the numerical simulations were performed to compare the hydrodynamic process and salinity distribution before and after a river training, which was conducted by blocking the south branch of the Oujiang River mouth. The results reveal that with the only north access to the sea, the influence of the blocking project on the flood discharge capacity is limited and the incremental velocity is beneficial to the navigation channel maintenance, although it will cause some scour to the embankment. Furthermore, the redistribution of tidal prism passing in or out the north branch makes a little severe salinity intrusion during high tide or low tide. However, the salinity intrusion is still within acceptable range, although it can cause some adverse effect on water intaking of production and life. The variations of salinity levels in Yueqing Bay situated at the north of the river mouth are not obvious, so the blocking project will not bring damage to local aquiculture. However, significant changes of salinity happen inside or outside of the south branch, so enough attention need to be paid to the changes of environment caused by the salinity variation after the blocking project. Overall, by weighing advantages and disadvantages of the blocking project, it is feasible and the model can be considered as a tool for managing and studying estuarine circulation.  相似文献   

5.
海南岛莺歌海近岸的潮汐不对称与潮致余流研究   总被引:3,自引:1,他引:2  
潮汐不对称与潮致余流在河口海岸区的物质输运中扮演着重要角色。已往的研究表明,在驻波占主导的河口海湾中,涨落潮的历时不对称与流速不对称有较为良好的对应关系。而潮致余流主要由地形与潮波的非线性作用所致。本研究以海南岛莺歌海附近为代表,结合实际观测与数值模型,研究复杂地形的开阔近岸区的潮汐不对称与潮致余流。结果表明,在莺歌海近岸区,涨落潮历时不对称皆表现为涨潮历时短于落潮历时,而流速不对称则出现复杂的空间变化。对流速不对称的机制分解表明,研究区的流速不对称主要由K1、O1与M2的相互作用,以及潮余流与各潮汐分潮的相互作用所控制。其中前者产生涨潮流速大于落潮流速的涨潮优势,而后者则与余流的方向相对应,出现多个涨潮优势与落潮优势的区域。总体而言,研究区的流速不对称由余流与各潮汐分潮的相互作用所决定。这表明,采用涨落潮历时的不对称来确定潮汐不对称的方法在开阔近海区可能并不适用。对潮致余流的研究表明,研究区的欧拉余流远大于斯托克斯余流。欧拉余流表现为多个顺时针与逆时针的涡流。涡流分布与地形具有较好的对应关系,潮流沙脊区多发育顺时针涡流,而深槽区则以发育逆时针涡流为主。摩擦力在涡流的发育中起着重要作用。  相似文献   

6.
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.  相似文献   

7.
In order to determine robustness of the models based on the viscous shallow-water equations in the boundary-fitted coordinates the solution characteristics of the boundary-value problem are compared for different forms of equations in contravariant, covariant and Cartesian presentations of the velocity. Distinctions between these forms associated with presentation of advection and eddy viscosity in their approximation on a curvilinear grid are established. A boundary-value problem for the tidal dynamics computation in the Strait of Messina (Mediterranean Sea) is formulated with use of different forms of equations in the boundary-fitted coordinates. An estimation and comparison of the modelling results are performed.  相似文献   

8.
Current, sea level and bed-load transport are investigated in the Lower Piscataqua River section of the Great Bay Estuary, New Hampshire, USA—a well-mixed and geometrically complex system with low freshwater input, having main channel tidal currents ranging between 0.5 and 2 m s−1. Current and sea level forced by the M2M4M6 tides at the estuarine mouth are simulated by a vertically averaged, non-linear, time-stepping finite element model. The hydrodynamic model uses a fixed boundary computational domain and accounts for flooding–drying of tidal flats by making use of a groundwater component. Inertia terms are neglected in comparison with pressure gradient and bottom friction terms, which is consistent with the observed principal dynamic balance for this section of the system. The accuracy of hydrodynamic predictions in the study area is demonstrated by comparison with four tidal elevation stations and two cross-section averaged current measurements. Simulated current is then used to model bed-load transport in the vicinity of a rapidly growing shoal located in the main channel of the lower system. Consisting of coarse sand and gravel, the shoal must be dredged every five to eight years. Two approaches are taken—an Eulerian parametric method in which nodal bed-load flux vectors are averaged over the tidal cycle and a Lagrangian particle tracking approach in which a finite number of sediment particles are released and tracked. Both methods yield pathways and accumulations in agreement with the observed shoal formation and the long-term rate of sediment accumulation in the shoal area.  相似文献   

9.
The nonhydrostatic boundary problem for an arbitrary three-dimensional domain with a seamount is considered. The problem is integrated into curvilinear boundary-fitted coordinates on a nonuniform grid. In order to identify nonhydrostatic effects the grid is condensed on the slopes of the seamount preserving a coarse resolution in the rest of the domain, where the problem is solved in the hydrostatic approximation. Calculation results for the nonhydrostatic tidal dynamics and hydrology of the Strait of Messina in the area of a seamount are presented.  相似文献   

10.
The surface M 2 tide in the Canadian Arctic Archipelago (CAA) is reproduced on the basis of the QUODDY-4 three-dimensional finite-element hydrodynamic model. Particular emphasis has been placed on comparing model estimates for the amplitudes and phases of tidal elevations and the parameters of ellipses (major semiaxis and eccentricity) of the barotropic tidal current velocity with observational data. We present their spatial distributions and the distributions of averaged (over a tidal cycle) values of the density, horizontal transfer, and dissipation rate of barotropic tidal energy. It is found that the CAA is a much less effective dissipator of barotropic tidal energy than the World Ocean.  相似文献   

11.
北部湾茅尾海是中国南方重要的经济开发区,兴建有滨海新城和多个港口码头,近年来围垦开发严重。本研究通过建立二维水动力数值模型系统,分析重要潮汐动力参数,对比研究1985年与2020年间钦州湾围垦和港口建设对当地水动力环境的影响。结果表明:经过围填海和港口工程之后,茅尾海的潮差变化较小,略微增加了0.05 m左右;全日分潮K1、O1及半日分潮M2、S2是影响钦州湾潮汐动力较大的驱动力,围垦后在茅尾海内海地区都略微增加了0.02~0.03 m,其中K1、O1是影响茅尾海的关键潮汐动力参数,敏感性测试分析表明三墩公路建设、钦州港海岸围垦和核电厂导堤建设对茅尾海潮差增加贡献率大致占60%、20%和10%;同时,围垦对束窄钦州湾航道具有一定的优化效应,围垦后钦州湾外湾三条水道峰值通量都明显增加,形成航道束水攻沙效果,对通航和维护主航道稳定性具有一定优势。因此,仅从潮汐动力参数变化角度分析,目前的围垦和港口工程迎合了当地河势特征,对潮汐动力场扰动较小,具有优化局部水动力场环境和提升通航安全性作用。  相似文献   

12.
A three-dimensional multi-level hydrodynamic model has recently been developed and applied to tidal motion in Singapore’s coastal waters. This paper describes a series of numerical experiments to evaluate the sensitivity of the tidal currents and elevations to model parameters. The results show that the predicted tidal elevations are insensitive to three model parameters: horizontal eddy viscosity coefficient (Smagorinsky constant, ch), bottom friction coefficient (cb) and internal friction coefficient (cv), whereas the effects of these parameters are quite different for tidal current velocities. The velocities are slightly reduced with an increase in ch and cb. The bottom friction effects on velocity profiles increase with water depth. The effect of cv might be significant for the tidal velocities at all levels. The velocities at upper layers of the water column decrease with the increase in cv, whereas the velocities at the bottom layer show the reverse trend. The effects of three model parameters on the magnitude and phase of the simulated currents are in the order (from strong to weak) of cv, cb and ch.  相似文献   

13.
—In this paper,integrating the Yangtze Estuary with the Hangzhou Bay,a 2-D velocity fieldmodel is established.In the model,fine self-adaptive grids are employed to adapt to the complicated coast-al shape.The hydrodynamic equations satisfied by two contravariant components of velocity vector andsurface elevation in non-orthogonal curvilinear coordinates are used.In each momentum equation thecoefficients before the two partial derivatives of surface elevation with respect to variables of alternative di-rection coordinates have different orders of magnitude,i.e.,the derivative with the larger coefficient mayplay a more important role than that with the smaller one.With this advantage,the ADI scheme can beeasily employed.The hydrodynamic factors include tidal current,river runoff and wind-induced current.In terms of tidal current,seven main constituents in the area are considered in the open boundaries.Theverifications of surface elevation process and current velocity process in the spring tide and in the n  相似文献   

14.
伶仃洋河口横向环流   总被引:6,自引:4,他引:2  
易侃  龚文平 《海洋学报》2015,37(3):1-14
横向环流存在于大多数河口中,在河口的动量平衡和物质输运过程中起着重要的作用,其形成的机制众多。本文利用三维水动力模型EFDC(Environmental Fluid Dynamics Code),结合1989年和2002年的实测数据对伶仃洋河口进行研究,得到了伶仃洋河口在洪季和枯季时期不同的水动力和盐度分布结构。在文中选取了具有代表性的两个断面,针对洪枯季、大小潮以及涨落潮等不同时期的横向流结构和盐度剖面分布特点进行比较和分析。结果显示,受局地的径、潮流条件和地形约束影响,横向环流结构并不能充分形成,在不同的区域表现出截然不同的结构特点。伶仃洋河口横向环流主要是在径、潮流与地形三者之间的相互作用下形成的,不同的区域的形成机制存在差异。综合分析后,本文认为伶仃洋河口的横向环流在该区域的动力过程中相对较为次要。  相似文献   

15.
An accurate prediction of ocean tides in southeast Alaska is developed using a regional, barotropic ocean model with a finite difference scheme. The model skill is verified by the observational tidal harmonics in southeast Alaska including Glacier Bay. The result is particularly improved in Glacier Bay compared to the previous model described by Foreman et al. (2000). The model bathymetry dominates the model skill. We re-estimate tidal energy dissipation in the Alaska Panhandle and suggest a value for tidal energy dissipation of 3.4 GW associated with the M2 constituent which is 1.5 times the estimation of Foreman et al. (2000). A large portion of the M2 energy budget entering through Chatham Strait is dissipated in the vicinity of Glacier Bay. Moreover, it is shown that the developed model has the potential to correct the ocean tide loading effect in geodetic data more efficiently than the model of Foreman et al. (2000), especially around Glacier Bay.  相似文献   

16.
Two models, a spectral refraction model (Longuet-Higgins) and a parabolic equation method (PEM) refraction-diffraction model (Kirby), are used to simulate the propagation of surface gravity waves across the Southern California Bight. The Bight contains numerous offshore islands and shoals and is significantly larger (≈ 300 km by 300 km) than regions typically studied with these models. The effects of complex bathymetry on the transformation of incident wave directional spectra, S0(f,θ0), which are very narrow in both frequency and direction are difficult to model accurately. As S0(f,θ0) becomes broader in both dimensions, agreement between the models improves and the spectra predicted at coastal sites become less sensitive to errors in the bathymetry grid, to tidal changes in the mean water depth, and to uncertainty in S0(f,θ0) itself. The smoothing associated with even relatively narrow (0.01 Hz-5° bandwidth) S0(f,θ0) is usually sufficient to bring the model predictions of shallow water energy into at least qualitative agreement. However, neither model is accurate at highly sheltered sites. The importance of diffraction degrades the predictions of the refraction model, and a positive bias [O (10%) of the deep ocean energy] in the refraction-diffraction model estimates, believed to stem from numerical “noise” (Kirby), may be comparable to the low wave energy. The best agreement between the predicted spectra generally occurs at moderately exposed locations in deeper waters within the Bight, away from shallow water diffractive effects and in the far-field of the islands. In these cases, the differences between the models are small, comparable to the errors caused by tidal fluctuations in water depth as waves propagate across the Bight. The accuracy of predicted energies at these sites is likely to be limited by the uncertainty in specifying S0(f,θ0).  相似文献   

17.
泰国湾及邻近海域潮汐潮流的数值模拟   总被引:2,自引:0,他引:2  
本文基于FVCOM(Finite-Volume Coastal Ocean Model)模式,模拟了泰国湾及其周边海域K1、O1、M2和S2四个主要分潮。采用47个验潮站实测调和常数与模拟结果进行比较,所得4个分潮的均方差分别为4.06cm、3.76cm、8.22cm和4.71cm,符合良好。根据计算结果分析了泰国湾及其周边海域的潮汐、潮流的分布特征和潮波的传播特征。数值试验表明,现有的数字水深资料(ETOPO1,ETOPO5,DBDB-V)的准确度不足以合理地模拟泰国湾潮波。  相似文献   

18.
Numericalsimulationofthree-dimensionaltidalcurrentintheBobalSea¥DenZhenxing;YangLianwu;andJoseOzer(ReceivedMay1,1993;accepted...  相似文献   

19.
许朋柱  毛锐 《海洋与湖沼》1993,24(1):109-113
为了能够精确地拟合天然水域不规则的水陆边界,并能保持有限差分数值计算在矩形网格中进行,本文利用Thompson的数值网格生成技术,在长江口南支七丫口至横沙河段设计了一个椭圆型边界拟合坐标系,并在这个坐标系下建立了平面二维潮流数学模型。通过数值求解此数学模型,实现了计算区域内现状潮流的一个数值模拟,模拟结果表明,模式的设计是成功的。  相似文献   

20.
The tidal regime of Shark Bay, Western Australia   总被引:1,自引:0,他引:1  
A non-linear hydrodynamic model is used to describe the tidal dynamics of Shark Bay, Western Australia. The model is forced by tidal elevations generated by M2, S2, K1 and O1 constituent data at the open boundaries. The absence of suitable boundary data required a ‘calibration’ of the boundary condition against the known constituent data from within the model domain. The model provides a good match to the available field data, and allows the surface-level and current response to be resolved over the entire domain. Due to a near quarter-wave resonance of the semi-diurnal tide along the eastern Hopeless Reach, which increases the semi-diurnal tide by a factor of 2, the tidal characteristics on each of the Reaches are different: on the eastern Hopeless Reach the tides are mainly semi-diurnal while on the western Freycinet Reach the tides are mainly diurnal. The tidal range is also higher along Hopeless Reach. Tidal harmonics, generated by non-linearity, are important in the shallow regions. The tidal wave is shown to propagate as a progressive wave into the Bay. Substantial phase-lag, attenuation and dissipation occur over the Faure Sill, a major shallow region of the eastern reach of the Bay. Non-linear generation of the M4 and MS4 tides is also significant in this region. Depth-averaged residual currents are presented, which show a tidally generated circulation that is enhanced in regions of complex topography. Estimates of tidal dissipation indicate that although the total dissipation is small on a global scale, the areal average is comparable with the Gulf of Carpentaria and approximately one-quarter of the value estimated for the Patagonian Shelf.  相似文献   

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

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