共查询到18条相似文献，搜索用时 375 毫秒
1.
一个两时间层分裂显格式海洋环流模式(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. 相似文献
2.
Surface currents measured by high frequency （HF） radar arrays are assimilated into a regional ocean model over Qingdao coastal waters based on Kalman filter method. A series of numerical experiments are per- formed to evaluate the performance of the data assimilation schemes. In order to optimize the analysis pro- cedure in the traditional ensemble Kalman filter （ENKF）, a different analysis scheme called quasiensemble Kaman filter （QENKF） is proposed. The comparisons between the ENKF and the QENKF suggest that both them can improve the simulated error and the spatial structure. The estimations of the background error covariance （BEC） are also assessed by comparing three different methods： Monte Carlo method; Canadian quick covariance （CQC） method and data uncertainty engine （DUE） method. A significant reduction of the root-mean-square （RMS） errors between model results and the observations shows that the CQC method is able to better reproduce the error statistics for this coastal ocean model and the corresponding external forcing. In addition, the sensibility of the data assimilation system to the ensemble size is also analyzed by means of different scales of the ensemble size used in the experiments. It is found that given the balance of the computational cost and the forecasting accuracy, the ensemble size of 50 will be an appropriate choice in the Qingdao coastal waters. 相似文献
3.
Diapycnal mixing plays an important role in the ocean circulation.Internal waves are a kind of bridge relating the diapycnal mixing to external sources of mechanical energy.Difficulty in obtaining eigen solutions of internal waves over curved topography is a limitation for further theoretical study on the generation problem and scattering process.In this study,a kind of transform method is put forward to derive the eigen solutions of internal waves over subcritical topography in twodimensional and linear framework.The transform converts the curved topography in physical space to flat bottom in transform space while the governing equation of internal waves is still hyperbolic if proper transform function is selected.Thus,one can obtain eigen solutions of internal waves in the transform space.Several examples of transform functions,which convert the linear slope,the convex slope,and the concave slope to flat bottom,and the corresponding eigen solutions are illustrated.A method,using a polynomial to approximate the transform function and least squares method to estimate the undetermined coefficients in the polynomial,is introduced to calculate the approximate expression of the transform function for the given subcritical topography. 相似文献
4.
Internal tide is one of the major oceanic phenomena. Determination of internal tide is important for theoretical study and for ocean engineering research. As an inverse problem, extraction of internal tidal currenls from sea currents is diffi-cult. In this paper, a method is developed to extract internal tidal currents from a portion of the sea current profile based on the fact that the directions of internal tidal currents above and below the thermocline are inverse. Sea current data col-lected from the South China Sea is processed with this method. The internal tidal currents and the depth of the thermocline are successfully extracted. The depth of the thermocline determined is in good agreement with that measured in 1959. 相似文献
5.
Numerical study about vortex-induced vibration(VIV) related to a flexible riser model in consideration of internal flow progressing inside has been performed.The main objective of this work is to investigate the coupled fluid-structure interaction(FSI) taking place between tensioned riser model,external shear current and upward-progressing internal flow(from ocean bottom to surface).A CAE technology behind the current research which combines structural software with the CFD technology has been proposed.According to the result from dynamic analysis,it has been found that the existence of upward-progressing internal flow does play an important role in determining the vibration mode(/dominant frequency),vibration intensity and the magnitude of instantaneous vibration amplitude,when the velocity ratio of internal flow against external current is relatively high.As a rule,the larger the velocity of internal flow is,the more it contributes to the dynamic vibration response of the flexible riser model.In addition,multi-modal vibration phenomenon has been widely observed,for asymmetric curvature along the riser span emerges in the case of external shear current being imposed. 相似文献
6.
Yu Ke Su Weidong Postgraduate River Harbour Department Nanjing Hydraulic ResearchInstitute Nanjing Engineer River Harbour Department Nanjing Hydraulic Research Institute Nanjing 《中国海洋工程》1992,(1)
Generally, the reliability of ami -sliding and anti-overturning stability of an isolated gravity cylinder in a certain working period can be evaluated only when the statistical properties of short term stability are given first. The authors used numerical method to simulate the stability state function of a cylinder in short-crested sea, and further to get the probabilitical characteristics of the structure's stability by time domain analysis. The external loads appeared in the state functions include horizontal wave force, lift force and the respective moments, and the loads are correlated by co- spectrum. The numerical method presented in this paper can be used not only to solve short term reliability problem directly, but to calculate and analyse the long term reliability problem as well. For circular cylinders, an example of simulation and analysis is displayed in this paper. 相似文献
7.
Li Dejun Shen Guoguang Wang Rixin Cheng Hong Professor Dept. of Ocean Architecture Engineering Tianjin University Tianjin Lecturer Dept. of Ocean Architecture Engineering Tianjin University Tianjin Engineer Naval Research Centre P. L. A. of China Beijing 《中国海洋工程》1995,(4)
Using the axial symmetry results of marker and cell (MAC) method as initial value in this paper, two numerical calculating methods are presented for the late wavemaking response induced by explosion in harbour. One of the methods is the superposition method of the vibration mode based on fluid slosh in container. Another one is the joining method of the MAC results with the shallow wave theory calculation in time domain. As a practical example, it is conducted to the numerical calculation about 1000 ton TNT equivalent explosion within touch of water surface. The results show that it can be rationally described with the methods to the wavemaking progress and character. The numerical results are identical with the observed scene on the spot experiment. The methods are simple and applicable in the engineering design. 相似文献
8.
Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieval,which is widely used in spaceborne and airborne radars. While the effect of the ocean currents and waves is interactional. It is impossible to retrieve the ocean surface current speed from Doppler frequency shift directly. In order to study the relationship between the ocean surface current speed and the Doppler frequency shift, a numerical ocean surface Doppler spectrum model is established and validated with a reference. The input parameters of ocean Doppler spectrum include an ocean wave elevation model, a directional distribution function, and wind speed and direction. The suitable ocean wave elevation spectrum and the directional distribution function are selected by comparing the ocean Doppler spectrum in C band with an empirical geophysical model function(CDOP). What is more, the error sensitivities of ocean surface current speed to the wind speed and direction are analyzed. All these simulations are in Ku band. The simulation results show that the ocean surface current speed error is sensitive to the wind speed and direction errors. With VV polarization, the ocean surface current speed error is about 0.15 m/s when the wind speed error is 2 m/s, and the ocean surface current speed error is smaller than 0.3 m/s when the wind direction error is within 20° in the cross wind direction. 相似文献
9.
《中国海洋工程》1989,(3)
A preliminary numerical model of the buoyant drum by wave action is presented, which is based on The Finite Boundary Element Method of the source distribution. A series formula of Green's function is used for the derivation of the model. In the paper the calculated results of the examples for a sphere and a buoyant drum are close to the data in some references. A subject of great interest to ocean engineers is the response a floating body suffers in ocean waves. In this paper, the oscilatory motion of a buoyant drum is analysed and a calculating method is presented. Applying this calculating method, the hydrodynamic coefficient for any kind of floating bodies can be obtained. 相似文献
10.
Modified PIC Method for Sea Ice Dynamics 总被引：1，自引：0，他引：1
The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation. 相似文献
11.
斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式 总被引：1，自引：0，他引：1
依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍 相似文献
12.
The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model 总被引：54，自引：1，他引：54
The purpose of this study is to find a combination of optimal numerical algorithms for time-stepping and mode-splitting suitable for a high-resolution, free-surface, terrain-following coordinate oceanic model. Due to mathematical feedback between the baroclinic momentum and tracer equations and, similarly, between the barotropic momentum and continuity equations, it is advantageous to treat both modes so that, after a time step for the momentum equation, the computed velocities participate immediately in the computation of tracers and continuity, and vice versa, rather than advancing all equations for one time step simultaneously. This leads to a new family of time-stepping algorithms that combine forward–backward feedback with the best known synchronous algorithms, allowing an increased time step due to the enhanced internal stability without sacrificing its accuracy. Based on these algorithms we design a split-explicit hydrodynamic kernel for a realistic oceanic model, which addresses multiple numerical issues associated with mode splitting. This kernel utilizes consistent temporal averaging of the barotropic mode via a specially designed filter function to guarantee both exact conservation and constancy preservation properties for tracers and yields more accurate (up to second-order), resolved barotropic processes, while preventing aliasing of unresolved barotropic signals into the slow baroclinic motions. It has a more accurate mode-splitting due to redefined barotropic pressure-gradient terms to account for the local variations in density field, while maintaining the computational efficiency of a split model. It is naturally compatible with a variety of centered and upstream-biased high-order advection algorithms, and helps to mitigate computational cost of expensive physical parameterization of mixing processes and submodels. 相似文献
13.
《Ocean Modelling》2008,20(1):61-89
This paper focuses on the energy conservation properties of a hydrostatic, Boussinesq, coastal ocean model using a classic finite difference method. It is shown that the leapfrog time-stepping scheme, combined with the sigma-coordinate formalism and the motions of the free surface, prevents the momentum advection from exactly conserving energy. Because of the leapfrog scheme, the discrete form of the kinetic energy depends on the product of velocities at odd and even time steps and thus appears to be possibly negative when high-frequency modes develop. Besides, the study of the energy balance clarifies the numerical choices made for the computation of mixing processes. The time-splitting technique used to reduce the computation costs associated to the resolution of surface waves leads to the well-known external and internal mode equations. We show that these equations do not conserve energy if the coupling of these two modes is forward in time. Even if non-linear terms are negligible, this shortcoming can be significant regarding the pressure gradient term ‘frozen’ over a baroclinic time step. An alternative energy-conserving time-splitting technique is proposed in this paper. Discussion and conclusions are conducted in the light of a set of numerical experiments dedicated to surface and internal gravity waves. 相似文献
14.
《Ocean Modelling》2002,4(3-4):249-267
During the course of developing new numerical algorithms for a terrain-following ocean modeling system (TOMS), different numerical aspects have been evaluated through a comparison between two widely used community ocean models, the Princeton ocean model (POM) and the regional ocean modeling system (ROMS). While both models aim at modeling coastal to basin-scale problems using similar grids, their numerical algorithms, code structure, and parameterization options are very different. Sensitivity studies with an idealized channel flow and a steep seamount configuration demonstrate how different algorithms in the two models may affect numerical errors, the stability of the code and the computational efficiency. For example, new pressure gradient schemes using polynomial fits and new time stepping algorithms may reduce numerical errors and allow using longer time steps than standard schemes do. However, the new schemes may require more careful choices of time steps and the use of higher order advection schemes to maintain numerical stability. 相似文献
15.
数据同化利用观测信息对模型状态场调整的同时也可以对数值模型中的不确定参数进行估计,从而改进数值模型,提高数值模拟的精度。本文基于集合调整卡尔曼滤波方法,采用广义坐标系统的美国普林斯顿大学海洋模式的外模式开展了渤海和部分黄海海域M2分潮模拟中的水深估计研究。理想数据同化试验结果表明,集合调整卡尔曼滤波方法能很好地降低模式模拟的水位误差并反演出“真实”的水深参数。而在NAO.99Jb和验潮站数据的实际数据同化试验中,与验潮站数据相比较,水深参数估计后,模式模拟的M2分潮振幅与迟角误差分别降低了40.27%和49.19%。 相似文献
16.
在中国近海的不少海域曾经发现相当强的内潮现象,例如赵俊生(1992)报道了渤海海峡的内潮现象, Yamashiro(1988)在东海陆坡处亦发现内潮的存在。特别是Amoco石油公司在南海陆坡处进行了一年多的海流测量,测得最大的全日潮流振幅可超过60cm/s,这里正压潮流小于10cm/s,故观测到的潮流主要与内潮有关。迄今关于内潮的数值模式大多将表面潮作为已知输入条件、从而使模式的应用受到较多限制(Jiang and Fang，1992)。Heaps(1983)曾发展了一种自由海面多层模式,但他们忽略了非线性平流项,这无疑也会损失许多有用信息。本文将叙述一种自由海面的多层模式,它能够同时模拟表面潮和内潮,并包含了非线性平流项。 相似文献
17.
将大气和海洋中f-平面上中尺度地形的扰动问题统一起来作为一个地球物理流体力学问题，应用涡度、能量和经圈动量守恒条件，将原来高度非线性方程变成一个二阶椭园型非线性方程，用数值模拟方法分析了大气中的过山运动和海洋中的沿岸上升流，给出了与观测接近的模拟结果以及运动过程对内部物理参数和外界条件的敏感性分析。所提方法可供进一步的动力学分析和数值研究参考。 相似文献
18.
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. 相似文献