共查询到20条相似文献，搜索用时 15 毫秒
1.
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. 相似文献
2.
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. 相似文献
3.
岛屿岛礁海域海浪能谱模型研究进展 总被引：2，自引：0，他引：2
波浪能谱模型在岛屿岛礁海域的波浪预报研究和海洋工程中应用广泛,但存在模式计算格点无法充分体现岛屿岛礁的复杂地形特征和很难刻画波浪受到岛屿岛礁影响发生变形物理过程等两个关键问题。多重网格嵌套方案、岛屿次网格地形效应计算方案以及非结构网格、无网格、动态自适应四叉树网格等技术在体现岛屿岛礁复杂地形方面取得了较好的效果;将相位解析模型与波浪能谱模型优势互补是提高能谱模型对岛屿近岸波浪变形物理过程计算能力的一个有效方法。开展球坐标系下波作用密度谱方程的自适应四叉树网格求解方法研究,借鉴相位解析模型最新成果完善能谱模式的绕射、反射、底摩擦等物理过程,是提高岛屿岛礁海域海浪精细预报技术水平的前沿性、探索性研究方向。 相似文献
4.
5.
强潮河口三维无结构网格盐度数学模型 总被引：3，自引：0，他引：3
采用平面无结构网格和垂向直接分层,建立了强潮河口盐度三维斜压力数学模型。无结构网格使模型能够适应河口复杂的边界,并可以根据需要进行局部加密。Casulli的半隐式法离散模型水位梯度和垂项紊动扩散项,克服了水位梯度和垂项紊动扩散对模型稳定性的影响,Semi-Lagrangian法用于模型对流项的离散,使模型具有"无条件"稳定性。采用水池风生流和盐度异重流算例检验了模型,并应用模型模拟了强潮河口(瓯江口)的潮流盐度运动,计算和实测结果进行比较,表明该模型是准确、可靠的。 相似文献
6.
《Ocean Modelling》2011,40(3-4):248-261
Accurate representation of geostrophic and hydrostatic balance is an essential requirement for numerical modelling of geophysical flows. Potentially, unstructured mesh numerical methods offer significant benefits over conventional structured meshes, including the ability to conform to arbitrary bounding topography in a natural manner and the ability to apply dynamic mesh adaptivity. However, there is a need to develop robust schemes with accurate representation of physical balance on arbitrary unstructured meshes. We discuss the origin of physical balance errors in a finite element discretisation of the Navier–Stokes equations using the fractional timestep pressure projection method. By considering the Helmholtz decomposition of forcing terms in the momentum equation, it is shown that the components of the buoyancy and Coriolis accelerations that project onto the non-divergent velocity tendency are the small residuals between two terms of comparable magnitude. Hence there is a potential for significant injection of imbalance by a numerical method that does not compute these residuals accurately. This observation is used to motivate a balanced pressure decomposition method whereby an additional “balanced pressure” field, associated with buoyancy and Coriolis accelerations, is solved for at increased accuracy and used to precondition the solution for the dynamical pressure. The utility of this approach is quantified in a fully non-linear system in exact geostrophic balance. The approach is further tested via quantitative comparison of unstructured mesh simulations of the thermally driven rotating annulus against laboratory data. Using a piecewise linear discretisation for velocity and pressure (a stabilised P1P1 discretisation), it is demonstrated that the balanced pressure decomposition method is required for a physically realistic representation of the system. 相似文献
7.
An unstructured-grid procedure for SWAN is presented. It is a vertex-based, fully implicit, finite difference method which can accommodate unstructured meshes with a high variability in geographic resolution suitable for representing complicated bottom topography in shallow areas and irregular shoreline. The numerical solution is found by means of a point-to-point multi-directional Gauss–Seidel iteration method requiring a number of sweeps through the grid. The approach is stable for any time step while permitting local mesh refinements in areas of interest. A number of applications are shown to verify the correctness and numerical accuracy of the unstructured version of SWAN. 相似文献
8.
《Coastal Engineering》2004,51(7):557-579
In this paper, a Reynolds Averaged Navier–Stokes (RANS) model was developed to simulate the vortex generation and dissipation caused by progressive waves passing over impermeable submerged double breakwaters. The dynamics of the turbulence are described by introducing a k–ɛ model with Boussinesq closure. The Height Function (HF) is implemented to define the free-surface configuration. The governing equations are discretized by means of a finite volume method based on a staggered grid system with variable width and height. The feasibility of the numerical model was verified through a series of comparisons of numerical results with the existing analytical solutions and the experimental data. The good agreements demonstrate the satisfactory performance of the developed numerical model. The flow separation mechanism both near the upstream and the downstream edges of the obstacles demonstrates the physical and expected nature of development of the flow. The present model provides an accurate and efficient tool for the simulation of flow field and wave transformation near coastal structures without breaking. 相似文献
9.
A BFG model for calculation of tidal current and diffusion of pollutants in nearshore areas 总被引：4，自引：0，他引：4
ＡＢＦＧｍｏｄｅｌｆｏｒｃａｌｃｕｌａｔｉｏｎｏｆｔｉｄａｌｃｕｒｒｅｎｔａｎｄｄｉｆｆｕｓｉｏｎｏｆｐｏｌｌｕｔａｎｔｓｉｎｎｅａｒｓｈｏｒｅａｒｅａｓ￥ＳｈｉＦｅｎｇｙａｎ；ａｎｄＺｈｅｎｇ；Ｌｉａｎｙｕａｎ（ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙ... 相似文献
10.
This paper presents the development and application of two-dimensional and three-dimensional oil trajectory and fate models for coastal waters. In the two-dimensional model, the oil slick is divided into a number of small grids and the properties of each grid due to spreading, advection, turbulent diffusion, evaporation and dissolution are studied. This model can predict the movement of the oil slick on the water surface. In order to simulate the distribution of oil particles in the water column, a three-dimensional oil fate model is developed based on the mass transport equation and the concentration distribution of oil particles can be solved. A comparison of numerical results with the observed data shows good conformity. 相似文献
11.
Large eddy simulations of the flow around a circular cylinder at high Reynolds numbers are reported. Five Reynolds numbers were chosen, such that the drag crisis was captured. A total of 18 cases were computed to investigate the effect of gridding strategy, turbulence modelling, numerical schemes and domain width on the results. It was found that unstructured grids provide better resolution of key flow features, when a ‘reasonable’ grid size is to be maintained.When using coarse grids for large eddy simulation, the effect of turbulence models and numerical schemes becomes more pronounced. The dynamic mixed Smagorinsky model was found to be superior to the Smagorinsky model, since the model coefficient is allowed to dynamically adjust based on the local flow and grid size. A blended upwind-central convection scheme was also found to provide the best accuracy, since a fully central scheme exhibits artificial wiggles, due to dispersion errors, which pollute the solution.Mean drag, fluctuating lift Strouhal number and base pressure are compared to experiments and empirical estimates for Reynolds numbers ranging from 6.31 × 104 to 5.06 × 105. In terms of the drag coefficient, the drag crisis is well captured by the present simulations, although the other integral quantities (rms lift and Strouhal number) show larger discrepancies. For the lowest Reynolds number, the drag is seen to be more sensitive to the domain width than the spanwise grid spacing, while at the higher Reynolds numbers the grid resolution plays a more important role, due to the larger extent of the turbulent boundary layer. 相似文献
12.
《Ocean Modelling》2008,20(2):170-182
Most ocean models in current use are built upon structured meshes. It follows that most existing tools for extracting diagnostic quantities (volume and surface integrals, for example) from ocean model output are constructed using techniques and software tools which assume structured meshes. The greater complexity inherent in unstructured meshes (especially fully unstructured grids which are unstructured in the vertical as well as the horizontal direction) has left some oceanographers, accustomed to traditional methods, unclear on how to calculate diagnostics on these meshes. In this paper we show that tools for extracting diagnostic data from the new generation of unstructured ocean models can be constructed with relative ease using open source software. Higher level languages such as Python, in conjunction with packages such as NumPy, SciPy, VTK and MayaVi, provide many of the high-level primitives needed to perform 3D visualisation and evaluate diagnostic quantities, e.g. density fluxes. We demonstrate this in the particular case of calculating flux of vector fields through isosurfaces, using flow data obtained from the unstructured mesh finite element ocean code ICOM, however this tool can be applied to model output from any unstructured grid ocean code. 相似文献
13.
SELFE: A semi-implicit Eulerian–Lagrangian finite-element model for cross-scale ocean circulation 总被引：1，自引：0，他引：1
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. 相似文献
14.
15.
预估修正法对河口海岸海洋模式稳定性的提高 总被引：4，自引：0，他引：4
河口海岸海洋模式时间差分为欧拉前差格式，在涡动粘滞系数较小的情况下，模式存在着弱不稳定性。为提高模式的稳定性，采用预估修正法对模式中科氏力项作半隐式处理。首先在理论上证明了模式稳定性的提高，随后设计一个高分辨率网格，考虑实际岸线和水深、边界通量、密度梯度力和风应力，把模式应用于冬季渤、黄、东海环流的数值模拟，数值试验也证明了模式稳定性的提高。改进后的河口海岸海洋模式较成功地模拟出了冬季渤、黄、东海环流（黑潮、台湾暖流、对马暖流、黄海暖流和沿岸流），与实际观测和研究结果较为一致。 相似文献
16.
Numerical modeling of nonlinear resonance of semi-enclosed water bodies: Description and experimental validation 总被引：1，自引：0，他引：1
Inigo J. Losada Jose M. Gonzalez-Ondina Gabriel Diaz-HernandezErnesto M. Gonzalez 《Coastal Engineering》2008
The purpose of this work is to a present a numerical model to solve a set of modified Boussinesq equations to analyse nonlinear resonance of semi-enclosed water bodies. The equations are solved on a finite element unstructured grid in order to achieve an optimal mesh resolution with the local geometry. The model is able to simulate long time lapses and realistic forcing in real bathymetries with a reasonable computational cost. To validate the numerical results, a set of experiments was carried out in a physical model of two adjacent elongated basins. Comparisons between numerical and experimental results for different geometries and nonlinear conditions show that the model is able to simulate with an excellent agreement the transient nonlinear resonant process. 相似文献
17.
18.
When a steep bottom slope exists, it is well known that conventional methods for calculating horizontal diffusion in sigma-coordinate coastal ocean models causes spurious transport (e.g. salinity, temperature, and sediments) and currents. In this study, a second-order accurate finite-difference algorithm and program have been developed to reduce the spurious numerical diffusion errors. In the proposed algorithm, the finite differencing is performed in the x-z coordinate system to approximate the horizontal gradient. Each variable in the finite differential formation is calculated in the sigma-coordinate grid cells using a second-order Lagrangian interpolation polynomial. In conjunction with a stepwise bottom boundary condition, numerical experiments show that the proposed finite-difference scheme considerably reduces numerical errors compared to conventional approaches when dealing with horizontal diffusion over steep topography, which often occurs in coastal oceans and navigation channels. 相似文献
19.
on the evolution equation for water waves,a mathematical model for wave propagation in large mild-slope areas is derived.The model is solved by the finite difference method with the staggered grid system.The computational results are in good agreement with experimental data and show that the model can obtain better results with relatively coarser grids.The model can be used to simulate water wave propagation in large coastal areas and can be efficiently solved without much programming effort. 相似文献
20.
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. 相似文献