共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper examines the simplification strategy of retaining only the nonhydrostatic effect of local acceleration in a three-dimensional fully nonhydrostatic model regarding the submesoscale wave phenomenon in the ocean.Elaborate scale analysis of the vertical component of the Reynold-averaged Navier–Stokes(RANS) equation was performed, confirming the rationalization of this simplification. Then, the simplification was implemented in a RANS equation-based nonhydrostatic model NHWAVE(nonhydrostat... 相似文献
2.
Coastal upwelling meanders and filaments are common features off eastern ocean boundaries. Their growth is reinvestigated herein using a nonhydrostatic three-dimensional model and a reduced-gravity model, with the objective of assessing contributions from two mechanisms that emerge in the nonhydrostatic regime. The first mechanism is caused by the vertical projection of the Coriolis force in the momentum equation. It is found that the vertical Coriolis force often acts as a restoring force against numerical damping off eastern ocean boundaries and thus enhances the growth of meanders and filaments. The second mechanism arises from unstable ocean stratification when the cold upwelled water intrudes seaward over the warm layer. The unstable stratification, albeit transient, further enhances the growth of meanders and filaments. It is concluded that although nonhydrostatic effects do not change our understanding of how meanders and filaments grow, the realism can be enhanced using a nonhydrostatic model insofar as meanders and filaments off eastern ocean boundaries are concerned. 相似文献
3.
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. 相似文献
4.
5.
A new method of assimilating sea surface height (SSH) data into ocean models is introduced and tested. Many features observable by satellite altimetry are approximated by the first baroclinic mode over much of the ocean, especially in the lower (but non-equatorial) and mid latitude regions. Based on this dynamical trait, a reduced-dynamics adjoint technique is developed and implemented with a three-dimensional model using vertical normal mode decomposition. To reduce the complexity of the variational data assimilation problem, the adjoint equations are based on a one-active-layer reduced-gravity model, which approximates the first baroclinic mode, as opposed to the full three-dimensional model equations. The reduced dimensionality of the adjoint model leads to lower computational cost than a traditional variational data assimilation algorithm. The technique is applicable to regions of the ocean where the SSH variability is dominated by the first baroclinic mode. The adjustment of the first baroclinic mode model fields dynamically transfers the SSH information to the deep ocean layers. The technique is developed in a modular fashion that can be readily implemented with many three-dimensional ocean models. For this study, the method is tested with the Navy Coastal Ocean Model (NCOM) configured to simulate the Gulf of Mexico. 相似文献
6.
We explore the efficacy of “super parameterization” (SP) in ocean modeling in which local 2-d non-hydrostatic plume-resolving fine-grained (FG) models are embedded at each vertical column of a coarse-grained (CG) hydrostatic model. A general multi-scale algorithm is described in which tendencies from the FG models are projected onto the CG model which in turn constrains the average state of the FG models, coupling the two models together. The approach is tested in the context of models of open ocean deep convection and compared with a pure hydrostatic, coarse resolution model using convective adjustment (HYD) and a full 3-d non-hydrostatic plume-resolving simulation (NH). The SP model is found to be greatly superior to HYD at much less computational cost than the fully non-hydrostatic calculation. 相似文献
7.
An algorithm is proposed for solving three-dimensional ocean hydrodynamics equations without hydrostatic approximation and traditional simplification of Coriolis acceleration. It is based on multicomponent splitting of the modified model with artificial compressibility. The original system of equations is split into two subsystems describing the transport of three velocity components and adjustment of the density and velocity fields. At the adjustment stage, the horizontal velocity components are represented as a sum of the depth means and deviations; the two corresponding subsystems are derived. For barotropic dynamics, the compressibility effect is represented as the boundary condition at the free surface, while for the baroclinic subsystem, it is introduced as ε-regularization of the continuity equation. Then, the baroclinic equations are split into two subsystems describing the hydrostatic and nonhydrostatic dynamics. The nonhydrostatic dynamics is computed at a separate splitting stage. The algorithm is included into the Institute of Numerical Mathematics of the Russian Academy of Sciences model based on “primitive” equations and verified by solving the hydrodynamics problem for the Sea of Marmara. 相似文献
8.
We present a simplified method for solving the local equilibrium carbonate chemistry in numerical ocean biogeochemistry models. Compared to the methods typically used, the scheme is fast, efficient and compact. The accuracy of the solution is dictated by the number of species retained in the expression for alkalinity and there is almost no computational penalty for retaining minor contributions. We demonstrate that this scheme accurately reproduces the results of the commonly used method in the context of a three-dimensional global ocean carbon cycle model. Using this model we also show that neglecting the regional variations in surface dissolved inorganic phosphorus and silicic acid concentrations can lead to significant systematic bias in regional estimates of air–sea carbon fluxes using such models. 相似文献
9.
10.
Through a suite of three-dimensional, high-resolution numerical modeling experiments, we examine the role of nonhydrostatic effects on O(1 km) submesoscale processes at ocean fronts, with particular focus on the vertical velocity field. Several differences between nonhydrostatic and hydrostatic models are pointed out using a framework that enables precise comparison, but it is difficult to identify categorical differences between the model solutions at the grid resolutions afforded. The instantaneous vertical velocity structure is sensitive to the model choice and, even more so, to grid resolution, but the average vertical flux is similar in both hydrostatic and nonhydrostatic cases.When a frontal region with horizontal density gradients is perturbed by wind, a profusion of submesoscale, O(1 km), secondary circulation features develops in the upper 50 m. Narrow, elongated cells of intense up- and down-welling are found to occur close to the surface, overlying broader regions of weaker up- and down-welling associated with the mesoscale meanders of the baroclinically unstable front. The submesoscale down-welling is considerably stronger than up-welling and is concentrated in 1–2 km width filaments within which velocities can attain magnitudes as high as 200 m day−1. The submesoscale features are found to be robust at horizontal grid resolutions varying between 1 and 0.25 km and exist even in the hydrostatic model. Submesoscale circulation is difficult to observe or resolve in coarser resolution circulation models, but is likely to play a significant role in the exchange of energy and properties between the surface ocean and thermocline. Possible mechanisms for the generation of these features are investigated in a follow-on paper. 相似文献
11.
A novel high resolution model without open boundary conditions applied to the China Seas: first investigation on tides 总被引：4，自引：1，他引：3
We developed a Global Ocean Circulation and Tide Model (GOCTM) with coarse grids in the open deep ocean degrading ‘smoothly’ into the highly resolved China Seas (CS) of refined grids to study the tides and circulation there.GOCTM is based on the framework of the Finite Volume approach for better mass conservation through improved transports across the discrete individual control volume.It also takes a full advantage of the geometric flexibility of unstructured mesh using a realistic global topography including the Arctic Ocean.The CS are given a special focus by refining the unstructured grids,but they are embedded into global domain naturally.Furthermore,GOCTM not only successfully avoids the treatment of the open boundaries,but also optimizes the trade-off between computational cost and model accuracy.Meanwhile,GOCTM is driven by the astronomical tide-generating potential and the secondary tide-generating potential directly,together with the wind stress and heat flux.GOCTM succeeds in reproducing the global eight principal tidal harmonic constants.Particularly,the simulated tidal results in the CS are improved compared to some other regional models with the discrepancy of 3.9 cm for M 2 tide.This idea of GOCTM can also be referred for other regional ocean study. 相似文献
12.
《Ocean Modelling》2002,4(2):121-135
Numerical studies of surface ocean fronts forced by inhomogeneous buoyancy loss show nonhydrostatic convective plumes coexisting with baroclinic eddies. The character of the vertical overturning depends sensitively on the treatment of the vertical momentum equation in the model. It is less well known how the frontal evolution over scales of O(10 km) is affected by these dynamics. Here, we compare highly resolved numerical experiments using nonhydrostatic and hydrostatic models and the convective-adjustment parametrization. The impact of nonhydrostatic processes on average cross-frontal transfer is weak compared to the effect of the O(1 km) scale baroclinic motions. For water-mass distribution and formation rate nonhydrostatic dynamics have similar influence to the baroclinic eddies although adequate resolution of the gradients in forcing fluxes is more important. The overall implication is that including nonhydrostatic surface frontal dynamics in ocean general circulation models will have only a minor effect on scales of O(1 km) and greater. 相似文献
13.
数值模式与统计模型相耦合的近岸海浪预报方法 总被引：2，自引：2，他引：0
针对数值模式和统计模型预报近岸海浪存在的局限性,构建了数值模式和统计模型相耦合的近岸海浪预报框架,在模式计算格点和近岸预报目标点之间定义一个海浪能量密度谱传递系数,通过经验正交函数分解和卡尔曼滤波方法建立传递系数的统计预报模型并与数值模式进行耦合。经过对近岸波浪观测站1a的预报试验表明:该方法能够提高近岸海浪有效波高预报精度,有效波高的均方根误差降低了约0.16m,平均相对误差降低约9%。进一步试验和分析发现,该方法的预报有效时间小于24h,将海浪能量密度谱经过分解后得到的基本模态反映了近岸波侯的主要特征,海浪能量密度谱传递系数的变化体现了波侯的季节变化特点。 相似文献
14.
基于推板造波理论和摇板造波理论,在Open FOAM平台上采用重叠网格技术建立黏性数值波浪水槽,并使用一种结合SIMPLE算法和PISO算法的PIMPLE算法对数值模型进行求解。利用开发的数值模型通过数值收敛性测试和网格独立性测试分别重点研究了时间步长、库朗数和网格尺寸对数值精度和计算效率的影响。并对比研究了此数值模型分别嵌入层流模型和湍流模型的计算精度和计算效率。实现的规则波和二阶有限振幅波与理论结果和试验结果吻合,验证了此黏性数值波浪水槽的造波和主动消波功能。基于二维数值波浪水槽,进一步研究了三维数值造波,数值计算结果与理论结果吻合良好。研究结果不仅验证了重叠网格在二维和三维两相流体域中求解运动物体与流场交互的可靠性和正确性,而且为使用此黏性数值波浪水槽解决更复杂的海洋工程问题提供了依据。 相似文献
15.
一个两时间层分裂显格式海洋环流模式(MASNUM)及其检验 总被引：1，自引：0，他引：1
HAN Lei 《海洋学报(英文版)》2014,33(11):11-35
A two-time-level, three-dimensional numerical ocean circulation model（named MASNUM） was established with a two-level, single-step Eulerian forward-backward time-differencing scheme. A mathematical model of large-scale oceanic motions was based on the terrain-following coordinated, Boussinesq, Reynolds-averaged primitive equations of ocean dynamics. A simple but very practical Eulerian forward-backward method was adopted to replace the most preferred leapfrog scheme as the time-differencing method for both barotropic and baroclinic modes. The forward-backward method is of second-order of accuracy, computationally efficient by requiring only one function evaluation per time step, and free of the computational mode inherent in the three-level schemes. This method is superior to the leapfrog scheme in that the maximum time step of stability is twice as large as that of the leapfrog scheme in staggered meshes thus the computational efficiency could be doubled. A spatial smoothing method was introduced to control the nonlinear instability in the numerical integration. An ideal numerical experiment simulating the propagation of the equatorial Rossby soliton was performed to test the amplitude and phase error of this new model. The performance of this circulation model was further verified with a regional（northwest Pacific） and a quasi-global（global ocean simulation with the Arctic Ocean excluded） simulation experiments. These two numerical experiments show fairly good agreement with the observations. The maximum time step of stability in these two experiments were also investigated and compared between this model and that model which adopts the leapfrog scheme. 相似文献
16.
17.
The aim of the current paper is to investigate hydrodynamic characteristics of the artificial upwelling induced by ocean currents. Experiments were performed in a flume at different density difference heads, horizontal current velocities and upwelling pipe diameters. A three-dimensional computational fluid dynamics (CFD) model was employed on wider range of parameters for further analysis. The performance of the numerical model has been confirmed by the experimental findings. The present results show that the volume flow rate of current-induced artificial upwelling is influenced by geometrical parameters and inclination angle of the pipe, the horizontal current velocity and vertical distribution of water density. In ideal two-layer density stratified water, the critical current velocity to generate upwelling linearly increases with the increase of the density difference, and the maximum rising height for upwelling is inversely proportional to the density difference. Feasibility analysis was taken by using current and density profiles of the East China Sea near Dongji Islands, which provides an useful reference for engineering practice. 相似文献
18.
Shenn-Yu Chao Ping-Tung Shaw Ming-Kuang Hsu Ying-Jang Yang 《Journal of Oceanography》2006,62(6):811-823
We have investigated the reflection and diffraction of first-mode and second-mode solitary waves by an island, using a three-dimensional
nonhydrostatic numerical model. The model domain consists of a circular island 15 km in diameter in an ocean 300 m deep. We
use prescribed density anomalies in an initially motionless ocean to produce highly energetic internal solitary waves; their
subsequent propagation is subject to island perturbations with and without the effect of earth’s rotation. In addition to
reflected waves, two wave branches pass around the island and reconnect behind it. Island perturbations to the first-mode
and second-mode waves are qualitatively similar, but the latter is more profound because of the longer contact time and, in
the presence of earth’s rotation, the scale compatibility between Rossby radius of the second baroclinic mode and the island
diameter. Without earth’s rotation, reflected and diffracted waves are symmetrical relative to the longitudinal axis passing
through the island center. With earth’s rotation, the current following the wave front veers to the right due to Coriolis
deflection. For a westward propagating incoming wave, the deflection favors northward wave propagation in the region between
the crossover point and the island, shifting the wave reconnection point behind the island northward. It also displaces the
most visible part of the reflected waves to the southeast. In the presence of earth’s rotation, a second-mode incoming wave
produces island-trapped internal Kelvin waves, which are visible after the passage of the wave front. 相似文献
19.
The application of very large floating structure (VLFS) to the utilization of ocean space and exploitation of ocean resources has become one of the issues of great interest in international ocean engineering field. Owing to the advantage of simplicity in structure and low cost of construction and maintenance, box-type VLFS can be used in the calm water area near the coast as the structure configuration of floating airport. In this paper, a 3D linear hydroelastic theory is used to study the dynamic response of box-type VLFS in sinusoidal regular waves. A beam model and a 3D FEM model are respectively employed to describe the dynamic characteristics of the box-type structure in vacuum. A hydrodynamic model (3D potential theory of flexible body) is applied to investigate the effect of different dry models on the hydroelastic response of box-type structure. Based on the calculation of hydroelastic response in regular waves, the rigid body motion displacement, flexible deflection, and the short term and long 相似文献
20.
A non-hydrostatic ocean model using an effective Poisson solver is developed. The Poisson solver is a combination of the multigrid method, the Krylov-subspace method, and the sparse approximate inverse. Its numerical cost only linearly increases with total number of computational cells, and it also has high parallel computing efficiency. The numerical cost of the non-hydrostatic model described in the present paper remains only twice of that of a hydrostatic model, even with non-smooth topography and with a huge number of computational grid cells on massively parallelized computer systems. Therefore, it has a potential to expand the applicability of non-hydrostatic ocean models. We also present the preliminary result of the high-resolution non-hydrostatic experiment on Ice Shelf Water overflow in the southern Weddell Sea, which shows good agreement with observations in terms of the pathway of dense water and velocity field. 相似文献