共查询到20条相似文献，搜索用时 15 毫秒
1.
《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°. 相似文献
2.
A. A. Slepyshev 《Izvestiya Atmospheric and Oceanic Physics》2016,52(3):301-308
Free internal waves are considered in a Boussinesq approximation in the situation when horizontal eddy viscosity and diffusion in a vertically inhomogeneous flow are taken into account. The dispersion relation and wave damping factor are found in a linear approximation. The Stokes drift velocity is determined in the second order of smallness based on the wave amplitude. It has been indicated that the Stokes drift velocity, transverse with respect to the wave propagation direction, differs from zero if the flow-rate transverse component depends on the vertical coordinate. Vertical momentum fluxes differ from zero and can be comparable with or exceed the corresponding turbulent fluxes if eddy viscosity and diffusion are taken into account. 相似文献
3.
Effects of adjusting vertical resolution on the eastern equatorial Pacific cold tongue 总被引：3，自引：0，他引：3
The vertical resolution of LICOM1.0 (LASG/IAP Climate System Ocean Model, version 1.0) is adjusted by increasing the level amount within the upper 150 m while keeping the total of levels. It is found that the eastern equatorial Pacific cold tongue is sensitive to the adjustment. Compared with the simulation of the original level scheme, the adjusting yields a more realistic strucature of cold tongue extending from the coast of Peru to the equator, as well as a temperature minimum at Costa Rica coast, north of the cold tongue. In the original scheme experiment, the sharp heating by net surface heat flux at the beginning of spin-up leads to a great warming in the eastern equatorial Pacific Ocean. The weak vertical advection due to a too thick mixed layer in the coarse vertical structure also accounts for the warm bias. The fact that most significant improvements of the upper 50 m temperature appear at the region of the thinnest mixed layer indicates the necessity of fine vertical resolution for the eastern equatorial Pacific Ocean. However, the westward extension of equatorial cold tongue, a defect in the original scheme, gets even more serious in the adjusting scheme due to the intensified vertical velocity and hence vertical advection in the central-eastern equatorial Pacific Ocean. 相似文献
4.
A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k-ε equations in this model. The eddy viscosity is calculated from the k-ε turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field,and then applied to simulate the tidal flow in the Bohai Sea. 相似文献
5.
Mayumi Murakami 《Journal of Oceanography》1986,42(3):224-239
The dynamical structure of a two-dimensional (depth and axial directions) estuary is studied analytically. A set of governing equations describing the time-averaged velocity and salinity in the estuary is used, where all of the external parameters (depth, width, freshwater discharge and horizontal and vertical coefficients of eddy viscosity and eddy diffusivity) are assumed to be constant along the estuary.Two dynamical relations are taken into consideration in the theory. One of them is the dependence of the longitudinal scaleL
d on the balance of longitudinal salt transport, and the other is the relation between the vertical stratification of salinity and the the Prandtl numberP
r=Av/Kv, whereA
v andK
v denote the coeffcients of the vertical eddy viscosity and diffusivity, respectively. The two relations result in an extension of the parameter range in which the linear balance of momentum holds.A linear state of motion (LM-state) is defined as the state where the momentum balance is linear. The LM-state comprises the so-called vertically homogeneous and the so-called partially mixed state. Perturbation analysis is introduced and dynamical theory is developed in the LM-state. Since the LM-state covers a fairly wide regime with respect to the balance of salt transport, the state is subdivided into three stages: the diffusive, intermediate and advective stages. In the diffusive stage the upstream salt transport is mainly attributed to the horizontal diffusion, whereas in the advective stage it is attributed to advection caused by gravitational circulation. The salinity balance is also linear in the diffusive and intermediate stages, whereas the balance is nonlinear in the advective stage. The advective stage of the LM-state is regarded as a stage bordering the salt wedge state.The longitudinal distribution of depth-mean salinity is found to take an exponential form in the diffusive stage, a nearly linear form in the advective stage and an intermediate form between them in the intermediate stage. 相似文献
6.
Using a parallel-track approach to estimate geostrophic surface velocities, an estimate of the statistics of ocean geostrophic surface currents and momentum stresses is provided on a 10 km along-track resolution from the first 49 repeat cycles (16 months) of the Jason-TOPEX/Poseidon tandem altimetric sea surface height (SSH) data. Results are compared with estimates obtained in a traditional way from along-track SSH data at crossover points and with in situ, Acoustic Doppler Current Profiler (ADCP) measurements obtained on board the VOS Oleander along a nominal path connecting Bermuda with the U.S. mainland. Agreements with the Oleander data are reasonable when simultaneous (in space and time) sampling is available. However, amplitudes of parallel-track geostrophic velocity variances are about 25% lower as compared to Oleander measurements which represent geostrophic and ageostrophic flow components. Estimates of velocity variances show clear signs of an anisotropic eddy field in the vicinity of all major current systems. At the same time estimates of Reynolds stresses and eddy momentum fluxes show a convergence of eddy momentum in all those regions, suggesting a forcing of the mean flow by the eddy field there. 相似文献
7.
Using a parallel-track approach to estimate geostrophic surface velocities, an estimate of the statistics of ocean geostrophic surface currents and momentum stresses is provided on a 10 km along-track resolution from the first 49 repeat cycles (16 months) of the Jason-TOPEX/Poseidon tandem altimetric sea surface height (SSH) data. Results are compared with estimates obtained in a traditional way from along-track SSH data at crossover points and with in situ, Acoustic Doppler Current Profiler (ADCP) measurements obtained on board the VOS Oleander along a nominal path connecting Bermuda with the U.S. mainland. Agreements with the Oleander data are reasonable when simultaneous (in space and time) sampling is available. However, amplitudes of parallel-track geostrophic velocity variances are about 25% lower as compared to Oleander measurements which represent geostrophic and ageostrophic flow components. Estimates of velocity variances show clear signs of an anisotropic eddy field in the vicinity of all major current systems. At the same time estimates of Reynolds stresses and eddy momentum fluxes show a convergence of eddy momentum in all those regions, suggesting a forcing of the mean flow by the eddy field there. 相似文献
8.
The upper ocean response to a moving typhoon 总被引：1，自引：0，他引：1
The upper ocean response to the translation speed of typhoons is studied using a three-dimensional primitive equation model.
Similar models studied previously have applied stability criteria rather than the diffusion term to simulate the vertical
mixing process. This study retains the diffusion term and uses the level-2 turbulence closure scheme to estimate the vertical
eddy viscosity. The model results indicate that in the forced period, the mixed-layer temperature decrease is greater for
a slow-moving storm due to stronger upwelling caused by the longer residence time. A fast-moving storm can attain a similar
cooling intensity in the wake period if its residence time allows the wind to resonate with the current. The significant downward
momentum diffusion and advection in the first few inertial periods of these events leads to strong, persistent inertial pumping
throughout the upper ocean in the wake period. The mixed layer is further cooled by turbulent mixing supported by vertical
current shears. Meanwhile, the upper thermocline exhibits a compensating temperature increase. The vertical transfer magnitude
and penetration scale are smaller in the slow-moving case, when the inertial motion decays rapidly. The model results also
indicate that the dominant cooling process can be inferred from the non-dimensional storm speed. However, this value may be
misleading for rapidly moving storms in which the current response is so distant from the storm that little wind work is performed
on the ocean. 相似文献
9.
Nonlinear axisymmetric oscillations of a warm baroclinic eddy are considered within the framework of an reduced-gravity model of the dynamics of a multilayer ocean. A class of exact analytical solutions describing pure inertial oscillations of an eddy formation is found. The thicknesses of layers in the eddy vary according to a quadratic law, and the horizontal projections of the velocity in the layers depend linearly on the radial coordinate. Owing to a complicated structure of the eddy, weak limitations on the vertical distribution of density, and an explicit form of the solution, the latter can be treated as a generalization of the exact analytical solutions of this form that were previously obtained for homogeneous and baroclinic eddies in the ocean. 相似文献
10.
Existing models of the wave bottom boundary layer have focused on the vertical and temporal dynamics associated with monochromatic forcing. While these models have made significant advances, they do not address the more complicated dynamics of random wave forcing, commonly found in natural environments such as the surf zone. In the closed form solution presented here, the eddy viscosity is assumed to vary temporally with the bed shear velocity and linearly with depth, however, the solution technique is valid for any eddy viscosity which is separable in time and space. A transformation of the cross-shore velocity to a distorted spatial domain leads to time-independent boundary conditions, allowing for the derivation of an analytic expression for the temporal and vertical structure of the cross-shore velocity under an arbitrary wave field. The model is compared with two independent laboratory observations. Model calculations of the bed shear velocity are in good agreement with laboratory measurements made by Jonsson and Carlsen (1976, J. Hydraul. Res., 14, 45–60). A variety of monochromatic, skewed, and asymmetric wave forcing conditions, characteristic of those found in the surf zone, are used to evaluate the relative effects on the bed shear. Because the temporal variation of the eddy viscosity is assumed proportional to the bottom shear, a weakly nonlinear interaction is created, and a fraction of the input monochromatic wave energy is transferred to the odd harmonics. For a monochromatic input wave, the ratio of the third harmonic of velocity at the bed to the first is <10%. However, for a skewed and asymmetric input wave, this ratio can be as large as 30% and is shown to increase with increasing root-mean-square input wave acceleration. The work done by the fluid on the bed is shown to be a maximum under purely skewed waves and is directed onshore. Under purely asymmetric waves, the work done is significantly smaller and directed offshore. 相似文献
11.
A complication of finite-volume triangular C-grid methods is the numerical emergence of horizontal divergence errors that lead to grid-scale oscillations in vertical velocity. Nonlinear feedback via advection of momentum can lead to numerical instability in velocity modes via positive feedback with spurious vertical velocities induced by horizontal divergence truncation error. Existing strategies to mitigate divergence errors such as direct divergence averaging and increased diffusion do not completely mitigate horizontal vertical velocity oscillations. We present a novel elliptic filtering approach to mitigate this spurious error and more accurately represent vertical velocities via improved calculation of horizontal divergences. These results are applied to laminar curved channel flows, demonstrating the applicability of the method to reproduce secondary flow features. 相似文献
12.
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. 相似文献
13.
HUANG Rui Xin 《海洋学报(英文版)》2014,33(3):19-39
Study of oceanic circulation and climate requires models which can simulate tracer eddy diffusion and ad vection accurately. It is shown that the traditional Eulerian coordinates can introduce large artificial hori zontal diffusivity/viscosity due to the incorrect alignment of the axis. Therefore, such models can smear sharp fronts and introduce other numerical artifacts. For simulation with relatively low resolution, large lateral diffusion was explicitly used in models; therefore, such numerical diffusion may not be a problem. However, with the increase of horizontal resolution, the artificial diffusivity/viscosity associated with hori zontal advection in the commonly used Eulerian coordinates may become one of the most challenging ob stacles for modeling the ocean circulation accurately. Isopycnal eddy diffusion （mixing） has been widely used in numerical models. The common wisdom is that mixing along isopycnal is energy free. However, a careful examination reveals that this is not the case. In fact, eddy diffusion can be conceptually separated into two steps： stirring and subscale diffusion. Due to the thermobaric effect, stirring, or exchanging water masses, along isopycnal surface is associated with the change of GPE in the mean state. This is a new type of instability, called the thermobaric instability. In addition, due to cabbeling subscale diffusion of water parcels always leads to the release of GPE. The release of GPE due to isopycnal stirring and subscale diffusion may lead to the thermobaric instability. 相似文献
14.
再谈海冰边缘区域中尺度涡旋形成机制——非线性平流 总被引：1，自引：1，他引：0
利用三维海洋模式与二维海冰模式耦合,研究海冰边缘区域中尺度涡旋形成最重要的机制之一——非线性平流机制。二维海洋模型模拟结果表明,非线性平流机制在水深比较浅的时候更加重要。不同于把海洋考虑成一个正压流体的二维模型,三维海洋模型中海冰通过海-冰相互作用直接影响海洋表层。我们发现在三维海洋模型实验中,中尺度涡旋和海洋表面抬升都对水深变化敏感。海流速度的垂直结构表面,当海水变浅,各层海流都变得更快。相同风应力作用相同时间之后,表面抬升与海水深度成反比关系。同时我们还发现由于垂直运动,在三维海洋模型实验结果中,海面抬升非常小,只有二维海洋模型实验结果的1%。垂直运动是三维海洋模型和二维海洋模型实验结果不同的根本原因。 相似文献
15.
Shizuo Tsunogai 《Journal of Oceanography》1972,28(4):145-152
A simple advection-diffusion model is applied to the deep water of the North Pacific Ocean. The physical mixing parameter, i.e., the ratio of vertical advection velocity (W) to vertical eddy diffusivity (D), is obtained from the vertical distribution of a conservative property such as salinity. The rate of decomposition of organic matter is estimated from the oxygen consumption rate which is obtained from dissolved oxygen content. The calcium carbonate flux in the deep water is obtained from alkalinity. Using these values and the vertical distribution of a radioisotope,14C or226Ra, the vertical eddy diffusivity and the upwelling velocity are found to be 1.2 cm2/sec and 1.4 ×10–5 cm/sec, respectively, at the Geosecs 1969 station. The oxygen consumption rate at 3 km depth of the station is found to be 1.4×10–3ml/l/yr. 相似文献
16.
Energetics of lateral eddy diffusion/advection：Part Ⅰ. Thermodynamics and energetics of vertical eddy diffusion 总被引：2，自引：2，他引：0
HUANG Rui Xin 《海洋学报(英文版)》2014,33(3):1-18
Two important nonlinear properties of seawater thermodynamics linked to changes of water density, cabbeling and elasticity(compressibility), are discussed. Eddy diffusion and advection lead to changes in density; as a result, gravitational potential energy of the system is changed. Therefore, cabbeling and elasticity play key roles in the energetics of lateral eddy diffusion and advection. Vertical eddy diffusion is one of the key elements in the mechanical energy balance of the global oceans. Vertical eddy diffusion can be conceptually separated into two steps: stirring and subscale diffusion. Vertical eddy stirring pushes cold/dense water upward and warm/light water downward; thus, gravitational potential energy is increased. During the second steps, water masses from different places mix through subscale diffusion, and water density is increased due to cabbeling. Using WOA01 climatology and assuming the vertical eddy diffusivity is equal to a constant value of 2×103 Pa2/s, the total amount of gravitational potential energy increase due to vertical stirring in the world oceans is estimated at 263 GW. Cabbeling associated with vertical subscale diffusion is a sink of gravitational potential energy, and the total value of energy lost is estimated at 73 GW. Therefore, the net source of gravitational potential energy due to vertical eddy diffusion for the world oceans is estimated at 189 GW. 相似文献
17.
为完全拟合河口近海复杂岸线和工程结构以及有效局部加密,设计并建立了一个无结构三角形网格二维河口海岸水动力数值模式。空间离散主要基于有限体积法以保证守恒性,时间积分采用预估修正法以提高精度。水位在三角形网格中心通过连续方程求解;水平x方向和y方向的流速U和V均在网格边中点上通过动量方程求解。流速平流项的求解中采用了TVD格式。TVD流速平流通量为一个一阶迎风格式通量和一个二阶格式通量的组合,一阶格式通量和二阶格式通量根据流速的局部分布情况得出配比,最终组合得到TVD通量。TVD格式具有低耗散和无频散的优点,提高了模式的稳定性。应用实测资料验证建立的模式,结果显示水位、流速和流向的计算值与实测值均符合良好。 相似文献
18.
Three velocity components of subsurface flow, observed in a rectangular tank under the action of a constant wind speed, are measured systematically at mesh points distributed uniformly over a vertical cross-section of the tank. Measurements are carried out for two cases: 1) reference wind speedU
r
=7.5 m/s and fetchF=10 m; and 2)U
r
=10 m/s andF=25 m. A pair of Langmuir cells is observed for both cases; downwelling zones are found along both of the sidewalls and an upwelling zone in the centre of the tank. Near the water surface, the vertical momentum flux is dominated by the Reynolds stress resulting from small-scale turbulence, while over the entire cross-section except near the surface, the Reynolds stress due to the Langmuir cells dominates the vertical momentum flux. As the result of the occurrence of this Langmuir cells, the vertical momentum flux, which consists of both mean advection and small-scale turbulence, is markedly inhomogeneous in the spanwise direction; for example, the largest vertical flux of the order of the wind stress is observed in the downwelling zone near one sidewall, while at the centre of the tank, the vertical momentum flux occupies only 30% of the wind stress. This indicates that a pair of Langmuir cells plays more important role than small-scale turbulence in the mixing process in a greater part of the wind-wave tank.Address after April 1, 1992: Department of Civil Engineering, Hiroshima Institute of Technology, Miyake 2-1-1, Saeki-ku, Hiroshima 731-51, Japan. 相似文献
19.
HUANG Rui Xin 《海洋学报(英文版)》2014,33(3):58-81
Gravitational potential energy （GPE） source and sink due to stirring and cabbeling associated with sigma dif fusion/ advection is analyzed. It is shown that GPE source and sink is too big, and they are not closely linked to physical property distribution, such as temperature, salinity and velocity. Although the most frequently quoted advantage of sigma coordinate models are their capability of dealing with topography; the exces sive amount of GPE source and sink due to stirring and cabbeling associated with sigma diffusion/advec tion diagnosed from our analysis raises a very serious question whether the way lateral diffusion/advection simulated in the sigma coordinates model is physically acceptable. GPE source and sink in three coordinates is dramatically different in their magnitude and patterns. Overall, in terms of simulating lateral eddy diffu sion and advection isopycnal coordinates is the best choice and sigma coordinates is the worst. The physical reason of the excessive GPE source and sink in sigma coordinates is further explored in details. However, even in the isopycnal coordinates, simulation based on the Eulerian coordinates can be contaminated by the numerical errors associated with the advection terms. 相似文献
20.
CHEN Yongping 《中国海洋工程》2001,(2):241-252
Based on the fluid motion equations, the physical meaning of eddy viscosity coefficient and the rationality of the Boussinesq hypothesis are discussed in this paper. The effect of the coefficient on numerical stability is analyzed briefly. A semi-enclosed rectangular sea area, with an orthogonal spur dike, is applied in a 2-D numerical model to study the effect of horizontal eddy viscosity coefficient (AH). The computed result shows that AH has little influence on the tidal level and averaged flow velocity, but has obvious influence on the intensity and the range of return flow around near the spur dike. Correspondingly, a wind-driven current pool and an annular current are applied in a 3-D numerical model respectively to study the effect of vertical eddy viscosity coefficient (Av). The computed result shows that the absolute value of Av is inversely proportional to that of horizontal velocity, and the vertical gradient value of Av determines the vertical distribution of horizontal velocity. The distrib 相似文献