共查询到20条相似文献，搜索用时 218 毫秒
1.
基于同量阶迭代法，在保留同阶面的前提下，对林建国等（1998a)得到的二阶Boussinesq类方程进行了求解，得到了与其量阶相对应的取立波解，并春与Euler方程的二阶孤立波解进行了比较，结果显示，本文解比传统Boussinesq方程的孤立波解有明显的改善，扩大了孤立的适用范围。 相似文献
2.
B. V. Khar’kov 《Oceanology》2007,47(2):161-171
The nonlinear dynamics of the low-frequency variability of a mid-latitude ocean are studied. The mechanism of the separation of the western boundary current from the western wall, as well as the meridional displacements of the separation point and the separated eastward jet, is analyzed. A regional barotropic quasigeostrophic eddy-resolving numerical model is used for the analysis. The flow in a rectangular domain is simulated by the constant inflow and outflow of fluid through the boundaries. A regime when the nonlinearity prevails over the dissipation and the advection and β terms are of the same order of magnitude is considered, which is characteristic of the actual ocean. When the nonlinearity exceeds the threshold value, a periodic solution is obtained. The solution is determined by the nonslip boundary condition at the western wall. The solution obtained is studied in detail. The meridional displacements of the western boundary current separation point with respect to the western wall and the separated eastward jet can reach a few hundred kilometers. Their intensities and the intensity of the recirculation gyre in the western boundary current are found to oscillate with a period of about five years. 相似文献
3.
4.
利用2003—2008年在渤海所测的212组气溶胶数据,研究了气溶胶光学厚度(AOT)、浑浊度系数(β)和ngstrm指数α的时空分布规律,并初步探讨了气溶胶特征参数之间的关系。结果表明,在时间分布方面,渤海气溶胶光学厚度、浓度和粒径尺度在一天中的变化比较小,中午AOT较大,10点左右AOT较小。从6月份到9月份,AOT、β和α的范围和均值均呈下降趋势;3月份AOT、β和α均明显小于六、8、9月份的对应值。离岸25km之内的大部分地区,AOT随离岸距离的增大而减小;但25km以外,气溶胶变化没有规律。比较而言,气溶胶AOT和β变化比较一致。浑浊度系数(浓度因子)β比ngstrm指数(尺度因子)α对气溶胶光学厚度AOT的影响更明显,AOT和β的关系可采用一阶线性回归方程表示。 相似文献
5.
任意水深变化Boussinesq型方程非线性波数值计算 总被引：6，自引：2，他引：4
首先从一个含有耗散项的高阶非线性和频散性波浪理论模型出发,建立了适用于任意底坡变化,相对水深h／Lo≤1的非线性波数学模型。应用全隐式交错网格和二阶精度中心差分法。得到离散方程组。进一步对其一阶导数项进行修正,达到与方程高阶项同量阶精度。精度检验表明本文计算结果与理论解和物理模型结果符合良好。 相似文献
6.
In this paper our previously developed advanced system identification technique [1] has been applied to extract the frequency dependent roll damping from a series of model tests run in irregular (random) waves. It is shown that this methodology accurately models the roll damping which can then be used to produce accurate predictions of the ships roll motion. These roll motion predictions are not only more accurate than the potential flow predictions but more accurate than potential flow models corrected using either empirical prediction methods [2] and even those corrected using roll damping obtained from free decay sallying experiments. This methodology has the potential to significantly improve roll motion prediction during full scale at sea trails of vessels in order to dramatically improve safety of critical operations such as helicopter landing or ship to ship cargo transfer. 相似文献
7.
This paper presents a potential based boundary element method for solving a nonlinear free surface flow problem for a ship moving with a uniform speed in finite depth of water. The free surface boundary condition is linearized by the systematic method of perturbation in terms of a small parameter up to third order. The surfaces are discretized into flat quadrilateral elements and the influence coefficients are calculated by Morino's analytical formula. Dawson's upstream finite difference operator is used in order to satisfy the radiation condition. The second order solution gives better result than the first or third order solution. So the present method with the second order solution can be adopted as a powerful tool for the hydrodynamic analysis of the thin ship in finite depth of water. 相似文献
8.
Hong Guangwen 《海洋学报(英文版)》1992,11(2):287-303
An exact analytic solution for wave diffraction by wedge or corner with arbitrary angle (rπ) and reflection coefficients (R0 and Rr) is presented in this paper. It is expressed in two forms-series and integral representations, corresponding recurrence relation and asymptotic expressions are also derived. The solution is simplified for some special cases of rπ. For Rr= R0,r= 1/N and Rr≠R0,r = 1/2N, the solution can be reduced to linear superpositions of incident and partially reflected waves, hence a nonlinear solution of forth order for this problem can be obtained by using the author's theory of nonlinear interaction among gravity surface waves. The given solution is related to inhomogeneous Robin boundary conditions, which include the Neumann boundary conditions usually accepted in wave diffraction theory. 相似文献
9.
Y. Chang D. Hammond A.C. Haza P. Hogan H.S. Huntley A.D. Kirwan B.L. Lipphardt V. Taillandier A. Griffa T.M. Özgökmen 《Ocean Modelling》2011,36(3-4):179-197
An investigation to improve trajectory prediction using Lagrangian data is presented. The velocity field of a data assimilating model, EAS-16, is corrected using drifter observations taken during an experiment off Taiwan. The results are tested using another independent Lagrangian data set provided by sonobuoys launched in the same area. The latter have instrument chains that extend well into the water column. Consequently the corrected model velocities were projected into the water column in order to calculate sonobuoy trajectories. The drifter and sonobuoy trajectories both show two distinct regimes in the considered area of approximately 1/2° square. One regime is dominated by shelf dynamics, the other by meandering of the Kuroshio, with a sharp boundary dividing the two. These two regimes are not reproduced by the trajectories of the EAS-16 model. When the drifter data are blended with the model velocities, synthetic sonobuoy trajectories track the observed ones much better, and the two regimes are clearly depicted. Two different methods for the velocity reconstruction are tested. One is based on a variational approach and the other on a normal mode decomposition. Both methods show qualitatively similar improvements in the prediction of sonobuoys trajectories, with a quantitative improvement in the total rms error of approximately 50% and 25%, respectively. 相似文献
10.
Ｏｎｈｙｄｒａｕｌｉｃｆａｌｌｓｏｆｔｗｏ－ｌａｙｅｒｆｌｏｗ￥ＸｕＺｈａｏｔｉｎｇ;ＬｏｕＳｈｕｎｌｉ;ＴｉａｎＪｉｗｅｉａｎｄＳａｍｕｅｌＳｈａｎｐｎＳｈｅｎ（ＩｎｓｔｉｔｕｔｅｏｆＰｈｙｓｉｃａｌＤｅｅａｎｏｐａｐｈｙａｎｄＰｈｙｓｉｃａｌｏｃ... 相似文献
11.
Study on Wave Dissipation of the Structure Combined by Baffle and Submerged Breakwater 总被引：1，自引：0，他引：1
This paper proposes a structure combined by baffle and submerged breakwater (abbreviated to SCBSB in the following texts). Such a combined structure is conducive to the water exchange in the harbor, and has strong capability on wave dissipation. Our paper focuses on the discussion of two typical structures, i.e., the submerged baffle and rectangular breakwater combined with the upper baffle respectively, which are named as SCBSB1 and SCBSB2 for short. The eigenfunction method corrected by experimental results is used to investigate the wave dissipation characteristics. It shows that the calculated results agree well with the experimental data and the minimum value of the wave transmission coefficient can be obtained when the distance between the front and rear structures is from 1/4 to 1/2 of the incident wave length. 相似文献
12.
A solution is presented for the wave induced drift forces acting on a submerged sphere in a finite water depth based on linearised velocity potential theory. In order to obtain the velocity potential, use has been made of multipole expansions in terms of an infinite series of Legendre functions with unknown coefficients. The series expression for the second order mean forces (drift forces) is provided by integrating the fluid pressure over the body surface. The horizontal drift force is also expressed by a series solution obtained using the far-field method. 相似文献
13.
14.
The recovery of quantities related to the gravity field (i.e., geoid heights and gravity anomalies) is carried out in a test area of the central Mediterranean Sea using 5' × 5' marine gravity data and satellite altimeter data from the Geodetic Mission (GM) of ERS‐J. The optimal combination of the two heterogeneous data sources is performed using (1) the space‐domain least‐squares collocation (LSC) method, and (2) the frequency‐domain input‐output system theory (IOST). The results derived by these methods agree at the level of 2 cm in terms of standard deviation in the case of the geoid height prediction. The gravity anomaly prediction results by the same methods vary between 2.18 and 2.54 mGal in terms of standard deviation. In all cases, the spectral techniques have a much higher computational efficiency than the collocation procedure. In order to investigate the importance of satellite altimetry for gravity field modeling, a pure gravimetric geoid solution, carried out in a previous study for our lest area by the fast collocation approach (FCOL), is used in comparison with the combined geoid models. The combined solutions give more accurate results, at the level of about 15 cm in terms of standard deviation, than the gravimetric geoid solution, when the geoid heights derived by each method are compared with TOPEX altimeter sea surface heights (SSHs). Moreover, nonisotropic power spectral density functions (PSDs) can be easily used by IOST, while LSC requires isotropic covariance functions. The results show that higher prediction accuracies are always obtained when using a priori nonisotropic information instead of isotropic information. 相似文献
15.
M. Taylan 《Ocean Engineering》1999,26(11):255
Mathematical modeling of the nonlinear roll motion of ships is one subject widely dealt with in nonlinear ship dynamics. This paper investigates setting up a form of nonlinear roll motion model and developing its periodic solution by the generalized Krylov–Bogoliubov asymptotic method in the time domain. In this model, nonlinearities are introduced through damping and restoring terms. The restoring term is approximated as a third-order odd polynomial whereas the quadratic term is favored to represent the nonlinear damping. The ship is assumed to be under the influence of a sinusoidal exciting force. Although the method is expressible to contain any order of the perturbing term, a single degree is chosen to avoid cumbersome mathematical complexity. In order to improve the solution a first-order correction term is also included. Moreover, a numerical example is carried out for a small vessel in order to validate the solution scheme. 相似文献
16.
Abstract We have obtained the residual depth and residual geoid anomalies of the South Atlantic Ocean and interpreted them in terms of upper mantle processes. Starting from the 5’ X 5’ SYNBAPS data, we computed the 1° X 1° mean water depth. We digitized a recent sediment thickness map of the South Atlantic and corrected for the isostatic sediment loading effects. A plot of the corrected basement depth against crustal ages shows a good match to the depth‐age curve of the plate model. We then subtracted the predicted plate model depths from the corrected basement depths and obtained the 1° X 1° residual depth. The residual depth anomalies have positive values over the topographic highs and relative negative values over the ocean basins. A prominent asymmetry is observed between the residual depth over the Argentine Basin and that over the Cape Basin. We have obtained the 1° X 1° residual geoid of the South Atlantic by subtracting both the long wavelength features and the geoid variations due to the plate cooling from the 1° X 1° Seasat altimeter derived geoid. The long wavelength features are represented by the degree and order 10 geoid of GEM1OC, and the geoid variations due to the plate cooling effects are predicted by the plate model geoid‐age relationship. The residual geoid anomalies also show an asymmetry although weaker than the case of the residual depth, between the Argentine and Cape basins. By taking the 5° X 5° averages, we removed possible plate bending effects on the depth and geoid anomalies. We then compared those two data sets with respect to the reported hot spot tracks in the South Atlantic. The residual depth generally shows positive values over the hot spot tracks, whereas the residual geoid does not show any sign of them. A prominent asymmetric feature of depth and geoid anomalies is observed over the Argentine and Cape basins. This asymmetry is probably caused by hotter and less dense materials beneath the Cape Basin. Hot spots or other mantle upwellings may be the heat sources. 相似文献
17.
In this paper we propose an integral form of the fully non-linear Boussinesq equations in contravariant formulation, in which Christoffel symbols are avoided, in order to simulate wave transformation phenomena, wave breaking and nearshore currents in computational domains representing the complex morphology of real coastal regions. Following the approach proposed by Chen (2006), the motion equations retain the term related to the approximation to the second order of the vertical vorticity. A new Upwind Weighted Essentially Non-Oscillatory scheme for the solution of the fully non-linear Boussinesq equations on generalised curvilinear coordinate systems is proposed. The equations are rearranged in order to solve them by a high resolution hybrid finite volume–finite difference scheme. The conservative part of the above-mentioned equations, consisting of the convective terms and the terms related to the free surface elevation, is discretised by a high-order shock-capturing finite volume scheme in which an exact Riemann solver is involved; dispersive terms and the term related to the approximation to the second order of the vertical vorticity are discretised by a cell-centred finite difference scheme. The shock-capturing method makes it possible to intrinsically model the wave breaking, therefore no additional terms are needed to take into account the breaking related energy dissipation in the surf zone. The model is verified against several benchmark tests, and the results are compared with experimental, theoretical and alternative numerical solutions. 相似文献
18.
Correction method for full-depth current velocity with lowered acoustic Doppler current profiler (LADCP) 总被引：1，自引：0，他引：1
A new method is presented to process and correct full-depth current velocity data obtained from a lowered acoustic Doppler
current profiler (LADCP). The analysis shows that, except near the surface, the echo intensity of a reflected sound pulse
is closely correlated with the magnitude of the difference in vertical shear of velocity between downcast and upcast, indicating
an error in velocity shear. The present method features the use of echo intensity for the correction of velocity shear. The
correction values are determined as to fit LADCP velocity to shipboard ADCP (SADCP) and LADCP bottom-tracked velocities. The
method is as follows. Initially, a profile of velocity relative to the sea surface is obtained by integrating vertical shears
of velocity after low-quality data are rejected. Second, the relative velocity is fitted to the velocity at 100–800 dbar measured
by SADCP to obtain an “absolute” velocity profile. Third, the velocity shear is corrected using the relationship between the
errors in velocity shears and echo intensity, in order to adjust the velocity at sea bottom to the bottom-tracked velocity
measured by LADCP. Finally, the velocity profile is obtained from the SADCP-fitted velocity at depths less than 800 dbar and
the corrected velocity shear at depths greater than 800 dbar. This method is valid for a full-depth LADCP cast throughout
which the echo intensity is relatively high (greater than 75 dB in the present analysis). Although the processed velocity
may include errors of 1–2 cm s−1, this method produced qualitatively good current structures in the Northeast Pacific Basin that were consistent with the
deep current structures inferred from silicate distribution, and the averaged velocities were significantly different from
those calculated by the Visbeck (2002) method. 相似文献
19.
The principal purpose of this paper is to extract entire sea surface wind＇s information from spaceborne lidar, and particularly to utilize a appropriate algorithm for removing the interference information due to white caps and subsurface water. Wind speeds are obtained through empirical relationship with sea surface mean square slopes. Wind directions are derived from relationship between wind speeds and wind directions im plied in CMOD5n geophysical models function （GMF）. Whitecaps backscattering signals were distinguished with the help of lidar depolarization ratio measurements and rectified by whitecaps coverage equation. Subsurface water backscattering signals were corrected by means of inverse distance weighted （IDW） from neighborhood non-singular data with optimal subsurface water backscattering calibration parameters. To verify the algorithm reliably, it selected NDBC＇s TAO buoy-laying area as survey region in camparison with buoys＇ wind field data and METOP satellite ASCAT of 25 km single orbit wind field data after temporal-spa tial matching. Validation results showed that the retrieval algorithm works well in terms of root mean square error （RMSE） less than 2m/s and wind direction＇s RMSE less than 21 degree. 相似文献
20.
孙孚 《海洋学报(英文版)》1987,6(1):8-19
Two-dimensional non-linear hydrodynamical equations are solved by using perturbation method and treating slopping beaches as bottom boundary conditions so that a kind of solution for nonlinear progressing waves is obtained. The first order of approximation is the same potential function as used by Biesel, and the second order is calculated numerically. Based on the solution, wave characteristics before breaking, especially the wave set-down, are discussed. It turns out that for the whole course of waves propagating from deep to shallow waters the theory proposed in this paper has a wider valid range of application than others. 相似文献