共查询到20条相似文献，搜索用时 162 毫秒
1.
An approach is developed to simulate wave–wave interactions using nonlinear elliptic mild-slope equation in domains where wave reflection, refraction, diffraction and breaking effects must also be considered. This involves the construction of an efficient solution procedure including effective boundary treatment, modification of the nonlinear equation to resolve convergence issues, and validation of the overall approach. For solving the second-order boundary-value problem, the Alternating Direction Implicit (ADI) scheme is employed, and the use of approximate boundary conditions is supplemented, for improved accuracy, with internal wave generation method and dissipative sponge layers. The performance of the nonlinear model is investigated for a range of practical wave conditions involving reflection, diffraction and shoaling in the presence of nonlinear wave–wave interactions. In addition, the transformation of a wave spectrum due to nonlinear shoaling and breaking, and nonlinear resonance inside a rectangular harbor are simulated. Numerical calculations are compared with the results from other relevant nonlinear models and experimental data available in literature. Results show that the approach developed here performs reasonably well, and has thus improved the applicability of this class of wave transformation models. 相似文献
2.
The problem of static three-dimensional, nonlinear, large deformation response of a marine riser is formulated within small strain theory and solved numerically. This type of analysis is necessary, for the new generation of drilling and production risers. The mathematical model takes properly into account the effects of internal and external pressure and the complete nonlinear boundary conditions, without linearizing the follower forces. The extensibility or inextensibility condition is used as the constitutive relation in the tangential direction. Torsion and bending are coupled. The external load and the boundary conditions are deformation dependent. A solution method is developed based on an incremental finite element algorithm, which involves a prediction-correction scheme. In the correction phase deformation dependent quantities are updated. The extensibility or inextensibility condition is used to reduce the degrees of freedom of the system. The numerical results of the developed computer code compare very well with available semi-analytical and numerical solutions. Three numerical applications are used to demonstrate the importance of large deformation, nonlinear and three-dimensional analyses. 相似文献
3.
In this work, transient resonant motions excited in linearised scattering interactions are approximated using a simple damped harmonic oscillator model. The scattering interactions considered involve the diffraction of an incident wave-train with a dominant regular time-harmonic component by structures which enclose a portion of the free surface. Provided a single resonant mode only is excited significantly, the fluid oscillation in the vicinity of the structure will primarily be composed of the resonant mode and incident wave mode contributions. The forced damped harmonic oscillator equation is used to predict the fluid motion and, in particular, the elevation of the internal free surface. The predictions are compared to the results from a numerical time-domain solver based on the linearised water-wave equations. It is shown that, given a good estimate of the location of the resonance in the complex frequency plane and a priori knowledge of the form of the incident wave, the model can successfully predict the time-dependent behaviour of the motion. Both two- and three-dimensional scattering problems are considered with a variety of scatterers in each case. 相似文献
4.
C.L. Kirk 《Applied Ocean Research》1985,7(1):2-13
This paper uses the Galerkin method in the solution of the marine riser differential equation and compares the dynamic bending stresses in a tension-leg-platform riser calculated by the linearised single wave and linearised spectral analysis methods.The results show that it is possible to make some rational assessment of approximate peak values of bending stress in the spectral method.The analysis methods are applicable to any riser system subject to horizontal exciting forces due to fluid and vessel motion. 相似文献
5.
斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式 总被引：1，自引：0，他引：1
依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍 相似文献
6.
将大气和海洋中f-平面上中尺度地形的扰动问题统一起来作为一个地球物理流体力学问题，应用涡度、能量和经圈动量守恒条件，将原来高度非线性方程变成一个二阶椭园型非线性方程，用数值模拟方法分析了大气中的过山运动和海洋中的沿岸上升流，给出了与观测接近的模拟结果以及运动过程对内部物理参数和外界条件的敏感性分析。所提方法可供进一步的动力学分析和数值研究参考。 相似文献
7.
The covariance matrix of sound-speed variations is determined from yo-yo CTD data collected during the SWARM 95 experiment at a fixed station. The data covered approximately 2 h and were collected during a period when nonlinear solitary internal waves were absent or negligible. The method of empirical orthogonal functions (EOF) is applied to the sound-speed covariance matrix assuming that the internal wave modes are uncorrelated. The first five eigenvectors are found to agree well with the theoretically modeled eigenfunctions based on the measured buoyancy frequency and the internal wave eigenmode equation. The mode amplitudes for the first five modes are estimated from the corresponding eigenvalues. They agree with the Garrett-Munk model if j*=1 is used instead of j*=3. A second method is used to deduce the mode amplitudes and mode frequency spectra by projecting the sound-speed variation (as a function of time) onto the theoretical mode depth functions. The mode amplitudes estimated with this method are in agreement with the EOF results. A modified Garrett-Munk model is proposed to fit the frequency spectrum of linear internal waves in shallow water 相似文献
8.
《Oceanic Engineering, IEEE Journal of》2004,29(4):1182-1199
In a recent study, satellite images have shown that internal solitons are active in the northern South China Sea (SCS). During the Asian Seas International Acoustic Experiment (ASIAEX) pilot studies, current profiler and thermistor chain moorings were deployed in the spring of 1999 and 2000 to investigate internal solitons northeast of Tung-Sha Island on the continental slope of the northern SCS. Most of the observed internal solitons were first baroclinic mode depression waves. The largest horizontal current velocity, vertical displacement, and temperature variation induced by the internal solitons were around 240 cm/s, 106 m, and 11/spl deg/C, respectively, while the estimated nonlinear phase speed was primarily westward at 152 /spl plusmn/ 4 cm/s. The observed internal solitons could be categorized as four types. The first type is the incoming wave from deep water and can be described reasonably well with the KdV equation. The second and third types are in the transition zone before and close to the turning point (where the upper and lower layer depths are equal), respectively. These two types of solitons were generally near the wave-breaking stage. The fourth type of soliton is a second baroclinic mode and probably was locally generated. The time evolutions are asymmetric, especially at the middle depths. A temperature kink following the main pulse of the soliton is often seen. Higher order nonlinear and shallow topographic effects could be the primary cause for these features. The appearance/disappearance of internal solitons coincides mostly with spring/neap tide. The internal soliton is irregularly seen during the neap tide period and its amplitude is generally small. The time interval between two leading solitons is generally around 12 h. The first baroclinic mode of the semidiurnal tide has a larger amplitude than the diurnal tide and could redistribute its energy into the soliton. 相似文献
9.
Excited by ocean currents, random wave and vessel motion, deepwater drilling risers exhibit significant dynamic response. In time domain, a method is proposed to calculate the nonlinear dynamic response of deepwater drilling risers subjected to random wave and dynamic large displacement vessel motion boundary condition. Structural and functional loads, external and internal pressure, free surface effect of irregular wave, hydrodynamic forces induced by current and wave, as well as wave and low frequency (drift) motion of the drilling vessel are all accounted for. An example is presented which illustrates the application of the proposed method. The study shows that long term drift motion of the vessel has profound effect on the envelopes of bending stress and lateral displacement, as well as the range of lower flex joint angle of the deepwater riser. It can also be concluded that vessel motion is the principal dynamic loading of nonlinear dynamic response for the deepwater risers rather than wave force. 相似文献
10.
In this paper, a beam without contact with water is called the “dry” beam and the one in contact with water is called the “wet” beam. For a partially (or completely) immersed uniform beam carrying an eccentric tip mass possessing rotary inertia, the conventional analytical (closed-form) solution is achieved by considering the inertial forces and moments of the tip mass and rotary inertia as the boundary conditions at the tip end of the beam. However, it has been found that the approximate solution for the last problem may be achieved by two techniques: Method 1 and Method 2. In Method 1, the basic concept is the same as the conventional analytical method; but in Method 2, the tip end of the beam is considered as a free end, while the inertial forces and moments induced by the tip mass and rotary inertia are considered as the external loads applied at the tip end of the beam. The main differences between the formulation of Method 1 and that of Method 2 are: In Method 1, the “normal” shapes of the “dry” beam are functions of the frequency-dependent boundary conditions but the external loads at the tip end are equal to zero; On the contrary, in Method 2, the “normal” mode shapes of the “dry” beam are determined based on the zero boundary conditions at the tip end of the beam but the external loads at the tip end due to the inertial effects of the tip mass and rotary inertia must be taken into consideration for the free vibration analysis of the “wet” beam. Numerical results reveal that the approximate solution obtained from Method 2 are very close to that from Method 1 if the tip mass moment of inertia is negligible. Besides, the two approximate solutions are also very close to the associated analytical (closed-form) solution or the finite element solution. In general, it is hoped that there exist several methods for tackling the same problem so that one may have more choices to incorporate with the specified cases. It is believed that the two approximate methods presented in this paper will be significant from this point of view. 相似文献
11.
在中国近海的不少海域曾经发现相当强的内潮现象,例如赵俊生(1992)报道了渤海海峡的内潮现象, Yamashiro(1988)在东海陆坡处亦发现内潮的存在。特别是Amoco石油公司在南海陆坡处进行了一年多的海流测量,测得最大的全日潮流振幅可超过60cm/s,这里正压潮流小于10cm/s,故观测到的潮流主要与内潮有关。迄今关于内潮的数值模式大多将表面潮作为已知输入条件、从而使模式的应用受到较多限制(Jiang and Fang，1992)。Heaps(1983)曾发展了一种自由海面多层模式,但他们忽略了非线性平流项,这无疑也会损失许多有用信息。本文将叙述一种自由海面的多层模式,它能够同时模拟表面潮和内潮,并包含了非线性平流项。 相似文献
12.
研究了一类非匹配不确定离散广义系统的准滑模控制问题。给出了带有扰动补偿的离散广义趋近律,消除了常规滑模控制中不确定项必须有界的限制,且不必满足匹配条件。所设计的滑模控制在有限时间内可达切换面,减小了准滑动模态带宽,有效地削弱了系统抖振,改善了系统动态品质。数值算例验证了该方法的可行性与有效性。 相似文献
13.
Risers/pipes conveying fluid are a typical kind of slender structures commonly used in marine engineering. It is of great academic significance and application value for us to evaluate and understand the vibration characteristics and nonlinear responses of these risers under the combined action of internal and external fluid flows. In this paper, the nonplanar vibrations and multi-modal responses of pinned-pinned risers in shear cross flow are numerically studied. With this objective in mind, the van der Pol wake oscillators are used to simulate the dynamical behavior of the vortex shedding in the wake. Two nonlinear equations of motion of the riser are proposed to govern the lateral responses of the riser structure. The nonplanar nonlinear equations for the riser and wake are then discretized by employing Galerkin's method and solved by using a fourth-order Runge–Kutta integration algorithm. Theoretical results show that the coupled frequencies for cross-flow (CF) and in-line (IL) motions and the corresponding coupled damping ratio could be influenced by the external and/or internal fluid velocities. Based on extensive calculations, the dynamical behavior of the riser with various internal and external flow velocities are presented in the form of bifurcation diagrams, time traces, phase portraits, oscillation trajectories and response spectrum curves. It is shown that some interesting dynamical phenomena, such as ‘lock-in’ state, ‘figure-of-eight’ trajectory and quasi-periodic oscillation, could occur in such a fluid-structure interaction system. Our results also demonstrate that the shear parameter can significantly affect the dynamic responses of the riser. When the shear parameter of the cross flow is large, multi-modal quasi-periodic responses of the riser can be excited, showing some new features undetected in the system of fluid-conveying risers in uniform cross flow. 相似文献
14.
A computational case study of coupled-mode 400-Hz acoustic propagation over the distance 27 km on the continental shelf is presented. The mode coupling reported here is caused by lateral gradients of sound-speed within packets of nonlinear internal waves, often referred to as solitary wave packets. In a waveguide having unequal attenuation of modes, directional exchange of energy between low- and high-loss modes, via mode coupling, can become time dependent by the movement of waves and can cause temporally variable loss of acoustic energy into the bottom. Here, that bottom interaction effect is shown to be sensitive to stratification conditions, which determine waveguide properties and, in turn, determine modal attenuation coefficients. In particular, time-dependent energy loss due to the presence of moving internal wave packets is compared for waveguides with and without a frontal feature similar to that found at the shelfbreak south of New England. The mean and variability of acoustic energy level 27 km distant from a source are shown to be altered in a first order way by the presence of the frontal feature. The effects of the front are also shown to be functions of source depth. 相似文献
15.
A coupled-mode model is developed for treating the wave–current–seabed interaction problem, with application to wave scattering by non-homogeneous, steady current over general bottom topography. The vertical distribution of the scattered wave potential is represented by a series of local vertical modes containing the propagating mode and all evanescent modes, plus additional terms accounting for the satisfaction of the free-surface and bottom boundary conditions. Using the above representation, in conjunction with unconstrained variational principle, an improved coupled system of differential equations on the horizontal plane, with respect to the modal amplitudes, is derived. In the case of small-amplitude waves, a linearised version of the above coupled-mode system is obtained, generalizing previous results by Athanassoulis and Belibassakis [J Fluid Mech 1999;389:275–301] for the propagation of small-amplitude water waves over variable bathymetry regions. Keeping only the propagating mode in the vertical expansion of the wave potential, the present system reduces to an one-equation model, that is shown to be compatible with mild-slope model concerning wave–current interaction over slowly varying topography, and in the case of no current it exactly reduces to the modified mild-slope equation. The present coupled-mode system is discretized on the horizontal plane by using second-order finite differences and numerically solved by iterations. Results are presented for various representative test cases demonstrating the usefulness of the model, as well as the importance of the first evanescent modes and the additional sloping-bottom mode when the bottom slope is not negligible. The analytical structure of the present model facilitates its extension to fully non-linear waves, and to wave scattering by currents with more general structure. 相似文献
16.
HAN Lei 《海洋学报(英文版)》2014,33(7):12-20
A time splitting technique is common to many free surface ocean models. The different truncation errors in the equations of the internal and external modes require a numerical adjustment to make sure that algorithms correctly satisfy continuity equations and conserve tracers quantities. The princeton ocean model （POM） has applied a simple method of adjusting the vertical mean of internal velocities to external velocities at each internal time step. However, due to the Asselin time filter method adopted to prevent the numerical instability, the method of velocity adjustment used in POM can no longer guarantee the satisfaction of the continuity equation in the internal mode, though a special treatment is used to relate the surface elevation of the internal mode with that of the external mode. The error is proved to be a second-order term of the coefficient in the Asselin filter. One influence of this error in the numerical model is the failure of the kinetic boundary condition at the sea floor. By a regional experiment and a quasi-global experiment, the magni- tudes of this error are evaluated, and several sensitivity tests of this error are performed. The characteristic of this error is analyzed and two alternative algorithms are suggested to reduce the error. 相似文献
17.
张力腿平台内孤立波作用特性数值模拟 总被引：1，自引：1，他引：0
依据三类内孤立波理论KdV、eK dV和MCC的适用性条件,采用Navier-Stokes方程为流场控制方程,以内孤立波诱导上下层深度平均水平速度作为入口边界条件,建立了两层流体中内孤立波对张力腿平台强非线性作用的数值模拟方法。结果表明,数值模拟所得内孤立波波形及其振幅与相应理论和实验结果一致,并且在内孤立波作用下张力腿平台水平力、垂向力及力矩数值模拟结果与实验结果吻合。研究同时表明,张力腿平台内孤立波载荷由波浪压差力、粘性压差力和摩擦力构成,其中摩擦力很小,可以忽略;水平力的主要成分为波浪压差力和粘性压差力,粘性压差力与波浪压差力相比较小却不可忽略,流体粘性的影响较小;垂向力中粘性压差力很小,流体粘性影响可以忽略。 相似文献
18.
The pycnocline in a closed domain is tilted by external wind forcing and tends to restore to a level posi- tion when the wind falls. An internal seiche oscillation exhibits if the forcing is weak, otherwise internal surge and internal solitary waves emerge, which serve as a link to cascade energy to small-scale processes. A two-dimensional non-hydrostatic code with a turbulence closure model is constructed to extend previous laboratory studies. The model could reproduce all the key phenomena observed in the corresponding labo- ratory experiments. The model results further serve as a comprehensive and reliable data set for an in-depth understanding of the related dynamical process. The comparative analyses indicate that nonlinear term favors the generation of internal surge and subsequent internal solitary waves, and the linear model predicts the general trend reasonably well. The vertical boundary can approximately reflect all the incoming waves, while the slope boundary serves as an area for small-scale internal wave breaking and energy dissipation. The temporal evolutions of domain integrated kinetic and potential energy are also analyzed, and the results indicate that about 20% of the initial available potential energy is lost during the first internal wave breaking process. Some numerical tactics such as grid topology and model initialization are also briefly discussed. 相似文献
19.
In this paper, motion response of a moored floating structure interacting with a large amplitude and steep incident wave field is studied using a coupled time domain solution scheme. Solution of the hydrodynamic boundary value problem is achieved using a three-dimensional numerical wave tank (3D NWT) approach based upon a form of Mixed-Eulerian–Lagrangian (MEL) scheme. In the developed method, nonlinearity arising due to incident wave as well as nonlinear hydrostatics is completely captured while the hydrodynamic interactions of radiation and diffraction are determined at every time step based on certain simplifying approximations. Mooring lines are modelled as linear as well as nonlinear springs. The horizontal tension for each individual mooring line is obtained from the nonlinear load-excursion plot of the lines computed using catenary theory, from which the linear and nonlinear line stiffness are determined. Motions of three realistic floating structures with different mooring systems are analyzed considering various combinations of linear and approximate nonlinear hydrodynamic load computations and linear/nonlinear mooring line stiffness. Results are discussed to bring out the influence and need for consideration of nonlinearities in the hydrodynamics and hydrostatics as well as the nonlinear modelling of the line stiffness. 相似文献
20.
Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｔｉｄａｌｃｕｒｒｅｎｔｉｎｔｈｅＢｏｂａｌＳｅａ￥ＤｅｎＺｈｅｎｘｉｎｇ;ＹａｎｇＬｉａｎｗｕ;ａｎｄＪｏｓｅＯｚｅｒ（ＲｅｃｅｉｖｅｄＭａｙ１,１９９３;ａｃｃｅｐｔｅｄ... 相似文献