共查询到20条相似文献，搜索用时 937 毫秒
1.
An unstructured-grid, finite-volume, nonhydrostatic, parallel coastal ocean simulator 总被引：4，自引：4，他引：4
A finite-volume formulation is presented that solves the three-dimensional, nonhydrostatic Navier–Stokes equations with the Boussinesq approximation on an unstructured, staggered, z-level grid, with the goal of simulating nonhydrostatic processes in the coastal ocean with grid resolutions of tens of meters. In particular, the code has been developed to simulate the nonlinear, nonhydrostatic internal wave field in the littoral ocean. The method is based on the formulation developed by Casulli, in that the free-surface and 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 with the second-order Adams–Bashforth method, while the pressure-correction method is employed for the nonhydrostatic pressure in order to achieve overall second-order temporal accuracy. Advection of momentum is accomplished with an Eulerian discretization which conserves momentum in cells that do not contain the free surface, and scalar advection is discretized in a way that ensures consistency with continuity, thereby ensuring local and global mass conservation using a velocity field that conserves volume on a local and global basis. The nonhydrostatic pressure field is solved efficiently using a block-Jacobi preconditioner, and while stability is limited by the internal gravity wave speed and vertical advection of momentum, applications requiring relatively small time steps due to accuracy or stability constraints are run efficiently on parallel computers, since the present formulation is written entirely with the message-passing interface (MPI). The ParMETIS libraries are employed in order to achieve a load-balanced parallel partitioning that minimizes interprocessor communication, and the grid is reordered to optimize per-processor performance by limiting cache misses while accessing arrays in memory. Test cases demonstrate the ability of the code to efficiently and accurately compute the nonhydrostatic lock exchange and internal waves in idealized as well as real domains, and we evaluate the parallel efficiency of the code using up to 32 processors. 相似文献
2.
强潮河口盐水入侵对饮用水源地危害极大。基于平面二维水动力盐度模型, 对典型强潮河口—钱塘江的水动力及盐水入侵过程进行了数值模拟研究。结果表明枯水径流时盐度变化与潮位过程曲线类似, 潮差对盐度大小影响显著, 径流量的增加将逐渐减小其相似程度。当流量增加到一定程度后, 继续增加的一定径流量所产生的抑咸效果减弱, 水资源有效利用率降低, 此时允许水源地盐度超标并改从蓄淡避咸水库取水可有效节约水资源。盐度平面分布显示, 盐水入侵在强潮河口弯道处受涨潮流主流线影响明显, 靠近主流线一岸的盐度大于对岸, 单从盐水入侵角度考虑, 强潮河口弯道段的取水口应设置在远离涨潮流主流线一岸。钱塘江河口盐度数值模拟对于研究减轻盐水入侵对水源地危害的措施具有指导意义。 相似文献
3.
A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary. 相似文献
4.
南海环流的三维数值模拟 总被引：1，自引：0，他引：1
本文采用Backhaus的三维海流模式,运用半隐式及C-网格方法求解基本方程,对南海各季的平均海流流场进行了数值模拟。将模拟结果与已有的研究结果进行比较,其主要流系基本上是相符的,夏季表层基本上为一反气旋型环流;冬季则转变成气旋型;在冬季,从50m层起开始显露出“南海暖流”的存在。这些表明南海海流的某些主要特征基本上已经被此模型成功地再现出来,同时此模型又给出了垂直方向各层的海流情况,在目前尚缺乏深层实测海流资料的情况下,上述深层海流的模拟结果有一定参考意义。 相似文献
5.
强潮河口三维无结构网格盐度数学模型 总被引：3，自引：0，他引：3
采用平面无结构网格和垂向直接分层,建立了强潮河口盐度三维斜压力数学模型。无结构网格使模型能够适应河口复杂的边界,并可以根据需要进行局部加密。Casulli的半隐式法离散模型水位梯度和垂项紊动扩散项,克服了水位梯度和垂项紊动扩散对模型稳定性的影响,Semi-Lagrangian法用于模型对流项的离散,使模型具有"无条件"稳定性。采用水池风生流和盐度异重流算例检验了模型,并应用模型模拟了强潮河口(瓯江口)的潮流盐度运动,计算和实测结果进行比较,表明该模型是准确、可靠的。 相似文献
6.
《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. 相似文献
7.
ECOM模式在丁字湾的应用 总被引：3，自引：3，他引：3
应用河口、陆架和海洋模式 (ECOM模型 ) ,引入干湿网格法模拟潮滩涨落的改进 ,并建立丁字湾及近岸海域的三维变动边界潮流模型。该模型考虑了湾口拦门沙、湾内水道和人工围海等地形特点。计算结果与实测值比较符合良好 ,较好地刻画出丁字湾 M2 分潮潮流场的时空分布特点。 相似文献
8.
四叉树网格下的椭圆型缓坡方程数值模型研究 总被引：1，自引：1，他引：0
波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。 相似文献
9.
斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式 总被引：1，自引：0，他引：1
依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍 相似文献
10.
A three-dimensional finite-difference hydrodynamic model has been developed using σ-coordinate for the vertical dimension. An explicit scheme for temporal integration and a staggered grid for spatial discretization have been adopted. The model has been tested against analytical or literature cases for wind and tide induced circulation. Results are in good agreement both with analytical solutions under idealised conditions and with results from the model of Shankar et al. (1996). 相似文献
11.
The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, a 3-D numerical model on an unstructured triangular grid has been developed. The σ coordinate transformation, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simu... 相似文献
12.
13.
本文基于2001年4月东中国海区域实测水文资料,应用三维有限元模式FEOM(Finite Element Ocean Model)对东中国海三维环流系统进行了数值计算分析。模式水平网格系统采用单节点线性三角形网格,垂向使用z坐标,观测温、盐度场通过客观分析法插值得到初始条件,分别进行了诊断计算和强诊断计算,计算结果表明:(1)改进逆方法可以很好地反演研究区域流函数和流量分布,为数值模拟提供优质可靠的开边界条件。(2)有限元模式在网格自由度方面和对研究区域的完整覆盖方面优势明显,高分辨率的垂向z坐标也可以较好地拟合海底地形,从而得到较高分辨率的三维数值模拟结果。(3)通过诊断计算,模拟再现了东中国海春季环流的多涡结构,分辨出了台湾暖流、黄海沿岸流、黄海暖流等流系。(4)比较诊断与强诊断两个计算结果,它们在定性上较为一致,在定量上有些差别,总体分布强诊断计算结果更为合理。 相似文献
14.
A wetting and drying algorithm is considered and implemented in a three-dimensional sigma–z coordinate model of ocean thermo- and hydrodynamics. The algorithm is tested in two idealized experiments simulating the run-up of a tidal wave on the coast and in a realistic experiment simulating the evolution of the Caspian Sea coastline in the 20th century. 相似文献
15.
Diapycnal mixing is important in oceanic circulation. An inverse method in which a semi-explicit scheme is applied to discretize the one-dimensional temperature diffusion equation is established to estimate the vertical temperature diffusion coefficient based on the observed temperature profiles. The sensitivity of the inverse model in the idealized and actual conditions is tested in detail. It can be found that this inverse model has high feasibility under multiple situations ensuring the stability of the inverse model, and can be considered as an efficient way to estimate the temperature diffusion coefficient in the weak current regions of the ocean. Here, the hydrographic profiles from Argo floats are used to estimate the temporal and spatial distribution of the vertical mixing in the north central Pacific based on this inverse method. It is further found that the vertical mixing in the upper ocean displays a distinct seasonal variation with the amplitude decreasing with depth, and the vertical mixing over rough topography is stronger than that over smooth topography It is suggested that the high-resolution profiles from Argo floats and a more reasonable design of the inverse scheme will serve to understand mixing processes. 相似文献
16.
Dmitry S. Dukhovskoy Steven L. Morey Paul J. Martin James J. O&#x;Brien Cortis Cooper 《Ocean Modelling》2009,28(4):250-265
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. 相似文献
17.
18.
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. 相似文献
19.
《Ocean Modelling》2004,6(2):125-150
We describe a three-dimensional (3D) finite-element ocean model designed for investigating the large-scale ocean circulation on time scales from years to decades. The model solves the primitive equations in the dynamical part and the advection–diffusion equations for temperature and salinity in the thermodynamical part. The time-stepping is implicit. The 3D mesh is composed of tetrahedra and has a variable resolution. It is based on an unstructured 2D surface mesh and is stratified in the vertical direction. The model uses linear functions for horizontal velocity and tracers on tetrahedra, and for surface elevation on surface triangles. The vertical velocity field is elementwise constant. An important ingredient of the model is the Galerkin least-squares stabilization used to minimize effects of unresolved boundary layers and make the matrices to be inverted in time-stepping better conditioned. The model performance was tested in a 16-year simulation of the North Atlantic using a mesh covering the area between 7° and 80° N and providing variable horizontal resolution from 0.3° to 1.5°. 相似文献
20.
We document the accuracy and convergence of solutions for a z-coordinate primitive-equation model of internal tide generation and propagation. The model, which is based on MOM3 numerics, is linearized around a state of rest to facilitate comparison with analytic estimates of baroclinic generation at finite-amplitude topography in a channel forced by barotropic tidal flow at its boundaries. Unlike the analytical model, the numerical model includes mixing of both buoyancy and momentum, and several definitions of “baroclinic conversion” are possible. These are clarified by writing out the energetics of the linearized equations in terms of barotropic kinetic energy, baroclinic kinetic energy, and available potential energy. The tidal conversion computed from the model, defined as the rate of conversion of barotropic kinetic energy into available potential energy, agrees well with analytical predictions. A comparison of different treatments of bottom topography (full-cells, partial-cells, and ghost-cells) indicates that the partial-cell treatment is the most accurate in this application. Convergence studies of flow over a smooth supercritical ridge show that the dissipation along tidal characteristics is, apparently, an integrable singularity. When the ocean bottom is not smooth, the accuracy and convergence of the model depend on the power spectrum of the topography. A numerical experiment suggests that the power spectrum of the resolved topography must roll off faster than k−2 to obtain convergent results from a linear numerical model of this type. 相似文献