共查询到10条相似文献，搜索用时 109 毫秒
1.
In numerical models of ocean circulation, it is widespread practice to split the fast and slow motions into barotropic and baroclinic subsystems, respectively. In the case of the baroclinic equations, the dependent variables can either be (1) slowly-varying baroclinic quantities, obtained from splitting the original flow variables into barotropic and baroclinic components, or (2) the original unsplit variables, which can vary on both the fast and slow time scales. In the second case, the variables in each layer are adjusted after each (long) baroclinic time step to ensure compatibility with the results produced from the barotropic equations. The second approach can be applied to the layer thickness equation to ensure exact conservation of mass within each layer. In the case of the momentum equations, the second approach amounts to replacing unresolved fast portions of Coriolis and pressure forcing with time averages of well-resolved forcing from the barotropic system. In this study, both approaches for the momentum equations are evaluated, in several test problems, by comparing to analytical solutions or to solutions computed with an unsplit code that uses short time steps. The two methods give very similar results in some simple problems for which analytical solutions are known. However, in some eddying double-gyre simulations, the formulation with unsplit variables requires a significant reduction in the baroclinic time step in order to avoid numerical difficulties that include grid noise and inaccurate representation of the flow field. In contrast, the formulation with split variables does not display such difficulties, and in those same examples it can be used with zero explicit horizontal viscosity. All of these computations employ a two-level time-stepping method that was previously developed by the author. 相似文献
2.
Three Dimensional Baroclinic Numerical Model for Simulating Fresh and Salt Water Mixing in the Yangtze Estuary 总被引：1，自引：5，他引：1
For simulating fresh and salt water mixing in estuaries, a three dimensional nonlinear baroclinic numerical model is developed, in which the gradients of horizontal pressure contain the gradient of barotropic pressure arising from the gradi-ent of tidal level and the gradient of baroclinic pressure due to the gradient of salinity. The Eulerian-Lagrangian method is employed to descretize both the momentum equations of tidal motion and the equation of salt water diffusion so as to im-prove the computational stability and accuracy. The methods to provide the boundary conditions and the initial conditions are proposed, and the criterion for computational stability of the salinity fields is presented. The present model is used for modeling fresh and salt water mixing in the Yangtze Estuary. Computations show that the salinity distribution has the characteristics of partial mixing pattern, and that the present model is suitable for simulalion of fresh and salt waler mixing in ihe Yanglze Esluary. 相似文献
3.
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. 相似文献
4.
一个两时间层分裂显格式海洋环流模式(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. 相似文献
5.
Baroclinic terms have been implemented in a three-dimensional fully hydrodynamic model developed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] modifying its momentum equations to account for density gradients and utilizing the scalar (salinity, temperature, etc.) conservation equation (SCE) and a state equation for the calculation of density. In the solution of advection-diffusion terms of the governing Navier-Stokes equations (NSE) and SCE, a symmetric splitting method was applied to ensure the long-term stability of simulations. Correction terms proposed by Ruddic et al. (1995) were applied to SCE to ensure the conservation of the scalar quantity. In the presence of baroclinic terms, the zero gradient pressure in the vertical direction in the vicinity of surface and bottom boundaries assumed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] created spurious currents. This problem was solved by assuming a hydrostatic pressure variation at those boundaries. The ability of extended model was validated by comparing its results with an experimental test case. The simulation of hydrodynamic and salt intrusion at Anzali Port located at the southern coasts of Caspian Sea in Iran was carried out by the model with both barotropic and baroclinic modes. The simulated results with baroclinic mode show a better agreement with measured data as compared to the results of barotropic mode that clearly demonstrate the significance of baroclinic terms in the simulation of cyclic intrusion of salt wedge into the Port Basin. 相似文献
6.
A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia. 相似文献
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.
《Ocean Modelling》2010,31(4):310-322
A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia. 相似文献
9.
A non-hydrostatic algorithm for the Regional Oceanic Modeling System (ROMS) is proposed. It is based on a decomposition technique for hydrostatic and non-hydrostatic pressure. The algorithm has a pressure-correction scheme with split-explicit time-stepping for baroclinic and barotropic vertical modes with a free surface. The algorithm implementation requires solving a Poisson equation for a non-hydrostatic pressure that has a non-symmetric matrix in discrete form. The efficiency of a different class of solvers and preconditioners were tested. The algorithm is successfully implemented with several examples where non-hydrostatic effects are important. These include standing external gravity waves; strongly nonlinear internal wave generation and transformation; stratified shear instability and its associated mixing; and nonlinear internal tidal generation over a ridge. The corresponding changes in the pre-processing and post-processing infrastructure in the existing hydrostatic ROMS code were performed to implement parallel elliptic solvers and a new set of dynamical equations. 相似文献
10.
A finite-difference system of motion equations for a baroclinic ocean with an irregular bottom is developed on the basis of the box method and the method of indefinite coefficients. The numerical model is shown to possess the property of conserving momentum, mass, and energy. In the case of barotropic motion, the finite-difference scheme retains potential enstrophy. A numerical experiment on the adaptation of hydrophysical fields in the Equatorial Atlantic in summer was carried out within the framework of the model developed. We managed to obtain a series of features of the equatorial circulation when the coefficients of turbulent exchange and diffusion were small enough.Translated by Mikhail M. Trufanov. 相似文献