首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
The Finite Element Ocean circulation Model (FEOM) is applied to study the sensitivity of density driven overflows to the vertical discretization and bottom topography representation using the dynamics of overflow mixing and entrainment (DOME) setup. FEOM allows for hybrid grids combining σ, z + σ, full cell, partly shaved cell and fully shaved cell grids within the same numerical kernel thus isolating as far as possible effects of mesh geometry from those of model numerics. The sensitivity of diapycnal mixing, entrainment, plume thickness and plume meridional distribution to vertical discretization and partly to the subgrid process parameterization is explored. It is shown that simulations on pure σ grids or the combination of z + σ resolve the overflow processes best in terms of downslope plume propagation, plume thickness and dilution, and also have the least resolution dependence. Grids using z-levels generate excessive spurious mixing when resolution is insufficient. Applying partial cells improves the plume representation, but still requires higher horizontal and vertical resolution to converge to the σ grid results. It is demonstrated that increasing lateral viscosity causes the plume thickness to reduce whereas increasing lateral diffusivity has opposite effect. When keeping the Prandtl number constant, the increase in diffusivity and viscosity leads to an increase in mixing and plume thickness on z-level grids and also on σ-grids when lateral dissipation is oriented along geopotential surfaces. Using the along σ- diffusion helped to obtain correct plume thickness and entrainment on σ grids. Increasing the vertical mixing coefficients leads to an increase in diapycnal mixing and in downslope penetration as well.  相似文献   

2.
Jan O. Backhaus   《Ocean Modelling》2008,22(3-4):114-127
This is the first part of a publication that describes the generation of adaptive grids (this part), and simulations with vector-ocean-model (VOM) in unstructured grids resulting from the adaptation (part II). A static vertical adaptive grid in z-coordinates allows improving the approximation of topography and vertical resolution at slopes. Adaptive grids use elements from a set of grid sizes by multiplying a basic smallest cell size with powers of two, as in cell division. Grids with locally isotropic vertical resolution at surface, seabed, and slopes can be generated whereby resolution decreases towards the ocean interior. The adaptation to topography yields unstructured grids that are organised in a one-dimensional vector by column-wise storage of cells, discarding land cells. The vector storage suggested the model’s name. Grids are generated by an iterative procedure that relies on rules, i.e. criteria and directives to control the grid structure in favour of a good representation of physics and smooth numerical operations. The directives govern vertical resolution at sea surface and seabed, and at slopes. For the latter vertical resolution is extended in the horizontal. In the ocean interior horizontal distances between changes in grid size can be controlled for the sake of smooth numerics. The use of a z-grid that avoids transformation errors, the depth-independence of vertical resolution, and the lateral extension of vertical resolution at slopes towards the ocean interior are the most significant differences of adaptive grids in comparison to vertical coordinate transformations. Unstructured grids do not rely on a smoothing of topography and can be used within any of the horizontal Arakawa-grids. For the same topography directives allow creating various grids as demonstrated for a shelf-ocean topography. The number of cells per column in two unstructured grids generated for the North Atlantic may locally well exceed typical layer numbers in conventional model matrices. But the domain average is similar to layer numbers of today’s ocean models. Thus, with the same investment of cells per domain a higher resolution in slope regions can be achieved by unstructured grids as compared to conventional z-grids.  相似文献   

3.
基于2000年8月在南海调查航次得到的水文资料,首次采用广义随底坐标形式的改进POM模式对南海夏季环流进行了数值研究.用正交曲线性水平网格覆盖观测区域,在垂向上对近表海面层次采用近似z坐标,而近底层则为随底坐标.在计算海区实际地形及假设的水平均匀而垂直层化的密度分布下,实施的两个数值计算试验表明,本模式采用的垂直坐标方案比传统的σ坐标方案优越,随底坐标模式因压力梯度项在起伏地形下产生的系统计算误差将变得十分的微小.在南海2000年夏季环流的实际计算中,首先对观测资料进行了60d的诊断计算,然后在诊断已得到的动力场结果基础上,又进行了10d左右的预报运行得到半诊断结果.从计算结果来看,它依赖于参数Cvis与Cdif的选择,特别是参数Cvis,文中取值为Cvis=Cdif=008.比较诊断与半诊断两个计算过程的结果,它们在定性上较为一致,在定量上有些差别.这是因为半诊断计算的方法对密度场作适当的动力调整,使其与地形、风场等更加匹配.在大尺度环流结构不受影响的情况下,尽可能地消除了小尺度噪声,可使计算得到的流场更为清晰.2000年8月南海计算区域环流的最大特点是多涡结构,其中有些反气旋暖涡和气旋式冷涡相间分布.在越南东南海域自表层至1000m水层稳定存在着一个显著的反气旋暖涡,其中心位置在11°51'N,112°07'E(诊断计算),水平尺度约为300km.此暖涡以东存在一个气旋式冷涡,这两个冷、暖涡是研究海区夏季环流的重要环流特征之一.在计算区域东北部夏季环流以反气旋环流系统为主;在计算区域东南部夏季环流以气旋系统为主;南海夏季环流分布,明显出现西部强化特征.  相似文献   

4.
MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Numerical simulations are carried out using three of the vertical coordinate types available in MPAS-Ocean, including z-star with partial bottom cells, z-star with full cells, and sigma coordinates. The results are first benchmarked against other models, including the MITgcm’s z-coordinate model and HIM’s isopycnal coordinate model, which are used to set the base case used for this work. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. Lastly, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates.  相似文献   

5.
A 1/8° global version of the Navy Coastal Ocean Model (NCOM) is described with details of its formulation, implementation, and configuration of the vertical coordinate. NCOM is a baroclinic, hydrostatic, Boussinesq, free-surface ocean model that allows its vertical coordinate to consist of σ coordinates for the upper layers and z-levels below a user-specified depth. This flexibility allows implementation of a hybrid σ–z coordinate system that is expected to mitigate some of the weaknesses that can be associated with either pure coordinate option. For the global NCOM application, the σ–z coordinate is used to allow terrain-following σ coordinates in the upper ocean, providing better resolution and topographic fidelity in shelf regions where flow is most sensitive to its representation. Including z coordinates for deeper regions efficiently maintains high near-surface vertical resolution in the open ocean. Investigation into the impact of the selected coordinate system focuses on differences between atmospherically-forced free-running (no assimilation) global solutions using σ–z and pure z coordinates. Comparisons with independent temperature observations indicate that global NCOM using the σ–z coordinate has improved skill relative to its z coordinate implementation. Among other metrics, we show that in comparison with time series of surface temperature from National Oceanic Data Center (NODC) buoys, mostly located in coastal regions, root mean squared differences (RMSD) improved for 63% and correlation improved for 71% of the stations when σ–z coordinates were used instead of pure z. For the exclusively open-ocean Tropical Atmosphere-Ocean (TAO) buoys, differences between the simulations were small, with the σ–z showing smaller RMSD for 45% of the stations and higher correlation for 65% of the stations. Additional comparisons using temperature profile observations further confirm a tendency for improved performance using the hybrid σ–z coordinates.  相似文献   

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.
Recent observations over the Sigsbee Escarpment in the Gulf of Mexico have revealed extremely energetic deep currents (near 1 m s−1), which are trapped along the escarpment. Both scientific interest and engineering needs demand dynamical understanding of these extreme events, and can benefit from a numerical model designed to complement observational and theoretical investigations in this region of complicated topography. The primary objective of this study is to develop a modeling methodology capable of simulating these physical processes and apply the model to the Sigsbee Escarpment region. The very steep slope of the Sigsbee Escarpment (0.05–0.1) limits the application of ocean models with traditional terrain-following (sigma) vertical coordinates, which may represent the very complicated topography in the region adequately, can result in large truncation errors during calculation of the horizontal pressure gradient. A new vertical coordinate system, termed a vanishing quasi-sigma coordinate, is implemented in the Navy Coastal Ocean Model for application to the Sigsbee Escarpment region. Vertical coordinate surfaces for this grid have noticeably gentler slopes than a traditional sigma grid, while still following the terrain near the ocean bottom. The new vertical grid is tested with a suite of numerical experiments and compared to a classical sigma-layer model. The numerical error is substantially reduced in the model with the new vertical grid. A one-year, realistic, numerical simulation is performed to simulate strong, deep currents over the Escarpment using a very-high-resolution nested modeling approach. The model results are analyzed to demonstrate that the deep-ocean currents in the simulation replicate the prominent dynamical features of the observed intense currents in the region.  相似文献   

8.
A three dimensional hydrostatic finite volume ocean model has been developed to solve the integral dynamical equations. Since the basic (integral) equations are solved for finite volumes rather than grid points, the flux conservation is easily enforced, even on arbitrary meshes. Both upwind and high-order combined compact schemes can be incorporated into the model to increase computational stability and accuracy. This model uses a highly distorted grid system near the boundary. The lateral boundaries of each finite volume are perpendicular to x and y axes and the two vertical boundaries are not purely horizontal. Four types of finite volumes are designed to follow the terrain with four (Type-A), three (Type-B), two (Type-C), and one (Type-D) vertices in the lower surface. Such a terrain-following grid discretization has superior features to z- and σ-coordinate systems. The accuracy of this model was tested.  相似文献   

9.
A static adaptive grid approximates the topography and defines the vertical resolution in Vector-Ocean-Model (VOM). The adaptation to topography creates unstructured grids, which are organised in a one-dimensional vector by column-wise storage of only wet cells. The model’s name reflects this data structure. The intention of VOM is better resolving flow and stratification near topographic boundaries in Z-coordinates. This is the second part of a publication that describes the generation of adaptive grids (part I), and simulations with VOM in unstructured grids (this part). Adaptive grids generated for a synthetic topography in part I include shelf, continental slope, and ocean. Three of those grids are here utilised in upwelling simulations. Under the same forcing increased vertical resolution at seabed and slopes yields a significant increase in flow energy as compared to coarser grids. Results allow explaining the surface intensification of a continental slope jet by vertical displacements of water masses in the seabed Ekman layer. Results in unstructured grids are almost identical to reference simulations in equidistant grids where the respective smallest grid size of unstructured grids was used. Negative effects of grids on predicted flow and stratification are absent also over particularly rough topography, as demonstrated by using vertical velocity as most sensitive indicator. In a further simulation an overflow governed by the advection of water mass properties is presented to demonstrate the conservation properties of the model. After 5 months of simulation the predicted domain average temperature deviated by 10−8 from the initial temperature field. Compared to equidistant grids the advection/diffusion scheme looses about one order of magnitude in accuracy when used in an unstructured grid. The results of VOM, being defined in Z-coordinates, are void of coordinate transformation errors. In an arbitrary topography unforced zero-flow remains quiescent in a stratification that only varies in the vertical. VOM due to its depth-independent vertical resolution appears particularly suitable for simulations of ocean-shelf exchange.  相似文献   

10.
Observations of turbulence, stratification, and mean current were made using a microstructure profiler and an acoustic Doppler current profiler (ADCP) during four cruises at a central location in the Ariake Sea, under weakly and strongly stratified conditions. Continuous measurements of the dissipation rate of turbulent kinetic energy (TKE), ε, were made. These revealed that frictional bed turbulence with quarterdiurnal variation in the bottom boundary layer (BBL) was one of the most energetic sources of vertical mixing in the sea. Thickness of the BBL was strongly confined by the stable stratification. We investigate a relationship between the BBL height h and the Ozmidov scale. We present a systematic argument that describes the vertical structure and characteristic scales of velocity and turbulence inside the frictional BBL, where the stratification persisted. Considerable deviation of observed vertical shear from the law of the wall indicated a modification of turbulent scales by the stratification. Shear stress calculated from the velocity data using vertical integration of the equation of motion was found to decrease approximately linearly with height. The TKE production rate P, estimated using the shear stress, was highly correlated with the dissipation rate. The buoyancy contribution to TKE balance in the BBL was quantified in terms of the flux Richardson number R f as R f?=?0.12.  相似文献   

11.
This chapter describes two separate but complementary research echosounder systems originally developed by Krupp Atlas Elektronik GmbH for the new German oceanographic research vessel Meteor.The Hydrosweep is a Hydrographic wide-swathe sweep survey echosounder for both shallow and deep water applications providing accurate bathymetric surveys and terrain-following navigation capabilities.The Parasound system is a hull-mounted dual channel parametric narrow-beam deep sea survey and sub-bottom profiling echosounder enabling particularly high vertical and horizontal resolution of seabed features.  相似文献   

12.
A simple model of lower trophic level ecosystem has been created to analyze possible environmental control of primary production in eight sub-areas of the Seto Inland Sea. The primary production rates observed by Hashimotoet al. (1997a) in these sub-areas are well reproduced by the model, including horizontal processes such as horizontal transport of nutrients and vertical processes such as vertical mixing, light intensity and sinking of particulate matter. Without taking account of horizontal processes the model also successfully reproduces the observed primary production rates in some areas, but if fails to reproduce those in the others. This shows that the relative importance of the horizontal transport on the primary production differs are by area. Two time scales,T z andT H, are introduced to explain this difference.T z is a vertical cycling time of material, which is defined as the time during which the stock of the material in the water column is utilized for primary production;T H is the horizontal transit time of the material. The relative importance of the horizontal process is well explained by theT H/T z ratio; that is, the horizontal transport process is important in the areas where this ratio is small. Further the possible mechanisms of nutrient supply for the primary production in each sub-area are investigated using this model.  相似文献   

13.
A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C -unstructured non-orthogonal grid in the horizontal direction and z -level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohai Sea. The results are in good agreement both with the analytical solutions and measurement results.  相似文献   

14.
Plant-flow interactions on the surface of tidal wetlands result in flow characteristics that are profoundly different from non-vegetated flows. Reductions in mean flow velocity and turbulence, especially the vertical components, limit vertical mixing and may impact a wide range of processes including geochemical exchanges at the sediment water interface, larval recruitment and dispersion, and sediment deposition and retention. The goal of this paper is to quantify horizontal and vertical components of velocity, turbulence intensity and total turbulent kinetic energy in Spartina alterniflora canopies in southeastern North Carolina and to relate flow characteristics to particulate transport on the marsh surface. Another aim of this paper is to assess the extent to which the distribution of standing biomass affects mean flow and turbulence by comparing S. alterniflora data to other canopy types and through a series of canopy manipulations which altered canopy height and stem densities.The results of this study indicate that flow velocity, turbulence intensity, and total turbulent kinetic energy (TKE) are significantly reduced within the vegetated canopy and that this reduction is inversely related to the amount of biomass present in the water column. Within the canopy, approximately 50% of the initial mean velocity and TKE is reduced within 5 m of the canopy edge. Within the canopy, mean velocity and TKEhoriz usually exceeded vertical velocity or TKEvert and the vertical components of flow were attenuated more strongly than the horizontal. These results suggest that within the vegetation, turbulence contributes more to lateral advection than to vertical mixing. As a result, total suspended solid concentrations were shown to decrease logarithmically with distance from the canopy edge and to decrease at a faster rate in more densely vegetated regions of the canopy (i.e. lower TKEvert) as compared to areas of sparser vegetation (i.e. higher TKEvert).  相似文献   

15.
We present the results of six dye tracer experiments that measured the mixing and circulation at the shelfbreak front on the New England Shelf. The last three were conducted during the New England Shelfbreak Productivity Experiment (NESPEX) with concurrent isopycnal float deployments. The results are consistent with the Chapman and Lentz [Chapman, D.C., and Lentz, S.J. (1994). Trapping of a coastal density front by the bottom boundary layer. Journal of Physical Oceanography, 24, 1465–1479.] model prediction of the separation and upwelling along the shelfbreak front of bottom boundary layer (BBL) water forced by an Ekman buoyancy flux, but show considerable variability. Cross-shelf velocities at the detachment point are 2–3 × 10−2 m/s. But seaward, over the slope region, dye tagged water was sheared from the main patch into small filaments that upwelled along the front with cross-shelf speeds up to 0.1 m/s. Cross-shelf diffusion was of order 10 m2/s in the mixed bottom layer and 1 m2/s in the interior along the front. Within the stratified front, the mean vertical diffusivity was Kz  4 × 10−6 m2/s. The dispersion of shelfwater in the slope region is effected by turbulent flow with advective speeds exceeding the small scale diffusive mixing. The mean flux of the detached BBL water is sufficient to account for the net loss of shelf water during its transit from Cape Cod to Cape Hatteras.  相似文献   

16.
Realistic representation of sea ice in ocean models involves the use of a non-linear free-surface, a real freshwater flux and observance of requisite conservation laws. We show here that these properties can be achieved in practice through use of a rescaled vertical coordinate “z*” in z-coordinate models that allows one to follow undulations in the free-surface under sea ice loading. In particular, the adoption of “z*” avoids the difficult issue of vanishing levels under thick ice.Details of the implementation within MITgcm are provided. A high resolution global ocean sea ice simulation illustrates the robustness of the z* formulation and reveals a source of oceanic variability associated with sea ice dynamics and ice-loading effects. The use of the z* coordinate allows one to achieve perfect conservation of fresh water, heat and salt, as shown in extended integration of coupled ocean sea ice atmospheric model.  相似文献   

17.
The influence of horizontal mixing on the thermal structure of the equatorial Pacific Ocean is examined based on a sigma coordinate model.In general,the distributions of the temperature and currents si...  相似文献   

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

19.
垂直分辨率对长江口海域M2分潮模拟的影响   总被引：1，自引：0，他引：1  
基于EFDC(Environmental Fluid Dynamics Code)模式建立了长江口及其邻近海域的三维水动力学模型, 研究模型的垂直分辨率对该海域M2分潮模拟的影响。结果表明:垂直分辨率的变化对M2分潮传播方向的模拟结果影响较小, 但其可通过底摩擦和湍流耗散两个计算过程来影响潮能通量的模拟结果, 最终对长江口和杭州湾内的M2分潮振幅产生显著的影响。最底层厚度较大时, 上层自由水体的高流速特征在最底层过于明显, 进而导致计算的底摩擦应力偏高, 此时提高底层的垂直分辨率会降低底摩擦对能量的耗散。另一方面, 垂直湍流混合作用会随垂直分辨率的增加而增强, 所以垂直分辨率增加到一定程度后, 上层自由水体的高流速会经由增强的湍流混合而更多的传入底层, 使计算的底摩擦应力随垂直分辨率的提高而有重新增加的趋势, 进而又增强底摩擦对潮能的耗散。  相似文献   

20.
《Ocean Modelling》2002,4(2):173-205
A three-dimensional hydrostatic model is presented that combines a generalised vertical co-ordinate system with an efficient implicit solution technique for the free surface. The model is capable of maintaining high resolution in the surface and/or bottom boundary layers as well as dealing with steep topography. Horizontal diffusion is calculated using the Smagorinsky formulation and a k–ε turbulence model is used in the vertical. In addition the model uses higher-order advection routines. An important aspect in three-dimensional models is the choice of vertical discretisation. If one is mostly interested in problems which are governed by boundary layer flows, a terrain following or sigma co-ordinate system seems attractive. This paper focuses on the development of a generalised sigma-type grid in a three-dimensional hydrostatic model. The generalised grid offers a wide range of possibilities including grid refinement toward the bed or surface, a mixed layer transformation, and a constant layer transformation where the lowermost or uppermost grid cells can be specified to have a constant height above the bed or below the surface. A number of tests are presented which show that the model is capable of simulating both shallow nearshore, estuarine flows as well as large-scale geophysical flows. These include an extreme flooding event in the shallow North Sea and the Odden ice tongue formation in the Greenland Sea.  相似文献   

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

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