共查询到20条相似文献，搜索用时 15 毫秒
1.
修正型缓坡方程的有限元模型 总被引：1，自引：1，他引：0
与缓坡方程相比,修正型缓坡方程增加了地形曲率项和坡度平方项,从而提高了数值求解的复杂性。本文将计算域划分为内域和外域,内域为水深变化区域,使用修正型缓坡方程,其中的地形曲率项和坡度平方项可用有限单元各节点的水深信息和单元插值函数表示,外域为水深恒定区,速度势满足Helmholtz方程,通过内外域的边界匹配建立有限元方程,并用高斯消去法求解。进而分别模拟了波浪传过Homma岛和圆形浅滩的变形,其结果与相关的解析解和实验数据吻合良好,证明了本文有限元模型的正确性。同时,通过与实验数据的对比也明显看出,在地形坡度较陡的情况下,修正型缓坡方程较缓坡方程具有更高的计算精度。 相似文献
2.
Lin Mingchung Hsiao Sungshan Hsu Yungcheng
Professor Dept. of Naval Architecture Ocean Engineering National Taiwan University Taiwan. Doctor Course Student Dept. of Naval Architecture Ocean Engineering National Taiwan University Taiwan 《中国海洋工程》1994,(3)
-Wave refraction-diffraction due to a large ocean structure and topography in the presence of a 'current are studied numerically. The mathematical model is the mild-slope equation developed by Kirby (1984). This equation is solved using a finite and boundary element method. The physical domain is devid-ed into two regions: a slowly varying topography region and a constant water depth region. For waves propagating in the constant water depth region, without current interfering, the mild- slope equation is then reduced to the Helmholtz equation which is solved by boundary element method. In varying topography region, this equation will be solved by finite element method. Conservation of mass and energy flux of the fluid between these two regions is required for composition of these two numerical methods. The numerical scheme proposed here is capable of dealing with water wave problems of different water depths with the main characters of these two methods. 相似文献
3.
为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析，建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域，每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解；而在没有结构物的外域，波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合，从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟，数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性；然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用，并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析，进一步证明了耦合模型的正确性与高效性。 相似文献
4.
A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable. 相似文献
5.
Second-order springing on an elastic body with forward speed is analyzed by numerical simulations. The boundary-value problem for the velocity potential is solved by means of the direct time-domain higher-order boundary element method (HOBEM). The free-surface boundary condition in the boundary-value problem is approximated on the mean surface up to second order by use of perturbation and Taylor-series expansion methods. The body boundary condition for an elastic body is derived with various quantities which are redefined in the generalized mode. These variables such as mode shape, normal vector, etc. are obtained by using directional derivative and continuum mechanics, and the same mathematical expressions are used to obtain several second-order generalized forces. To validate the numerical results, the second-order hydrodynamic force on the bottom-mounted rigid/elastic cylinders without forward speed is compared with other semi-analytic results. The property of second-order forces on an elastic ship is studied by changing the flexural rigidity and forward speed with elastic response. It is confirmed that the second-order velocity potential is important for a body with forward speed and investigation should be made more on numerical methods for accurate computation of the second-order velocity-potential force with forward speed. 相似文献
6.
系泊船非线性波浪力时域计算: 二维模型 总被引：6，自引：1，他引：6
为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算. 相似文献
7.
8.
波浪与外壁透空双方形沉箱相互作用 总被引：1，自引：1，他引：0
在线性波浪理论下,利用复合边界元素法(composite BEM)数值解析在等水深、规则波浪入射二种外壁透空双方形沉箱的无因次波力及Kd绕射分布图,并分别与其他研究者所作双圆筒内、外圆柱无因次波力及绕射分布图的计算结果进行比较,都说明本模式的合理性和可行性。在考虑不同透水参数下,分别计算波浪作用在内、外结构物的波力及沉箱四周绕射系数大小分布图。计算结果显示:波浪作用于外壁透空全透水双方形沉箱,外方柱无因次最大波力值会随着透水参数增加而降低;内方柱无因次最大波力值则随着透水参数增加而增加。本研究结果可供设计外海透水方形沉箱结构交互作用参考。 相似文献
9.
10.
A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropriate wave theory. At the downstream boundary, an artificial damping zone is used to prevent wave reflection back into the computational domain. Using the image Green function in the whole fluid domain, the integrations on the two lateral walls and bottom are excluded. The simulation results on extreme wave elevations in finite and infinite water-depths are compared with experimental results and second-order analytical solutions respectively. The wave kinematics is also discussed in the present study. 相似文献
11.
12.
The performance of coastal vertical seawalls in extreme weather events is studied numerically, aiming to provide guidance in designing and reassessing coastal structures with vertical wall. The extreme wave run-up and the pressure on the vertical seawall are investigated extensively. A time-domain higher-order boundary element method (HOBEM) is coupled with a mixed Eulerian-Lagrangian technique as a time marching technique. Focused wave groups are generated by a piston wave-maker in the numerical wave tank using a wave focusing technique for accurately reproducing extreme sea states. An acceleration-potential scheme is used to calculate the transient wave loads. Comparisons with experimental data show that the extended numerical model is able to accurately predict extreme wave run-ups and pressures on a vertical seawall. The effects of the wave spectrum bandwidth, the wall position and the wave nonlinearity on the wave run-up and the maximum wave load on the vertical seawall are investigated by doing parametric studies. 相似文献
13.
The resonance phenomenon of fluid motions in the gap between ship section, seabed and vertical quay wall is studied numerically and experimentally. The natural frequency of the fluid motions in the gap is derived. A two-dimensional time-domain coupled numerical model is developed to calculate the non-linear wave forces acting on a ship section against vertical quay in a harbor. The fluid domain is divided into an inner domain and an outer domain. The outer domain is the area between the left side of ship section and the incident boundary, where flow is expressed by Boussinesq equations. The rest area is the inner domain, which is the domain beneath the ship section plus the domain between the right side of ship section and vertical quay wall. The flow in the inner domain is expressed by Newton's Second Law. Matching conditions on the interface between the inner domain and the outer domain are the continuation of volume flux and the equality of pressures. The numerical results are validated by experimental data. 相似文献
14.
A 3-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear waves acting on a box-shaped ship fixed in a harbor. The domain is divided into the inner domain and the outer domain. The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations. The remaining area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions. Along the interface boundaries between the inner domain and the outer domain, the volume flux is assumed to be continuous and the wave pressures are equal. Relevant physical experiment is conducted to validate the present model and it is shown that the numerical results agree with the experimental data. Compared the coupled model with the flow in the inner domain governed by the Laplace equation, the present coupled model is more efficient and its solution procedure is simpler, which is particularly useful for the study on the effect of the nonlinear waves acting on a fixed box-shaped ship in a large harbor. 相似文献
15.
-In this paper, an analytical solution in the outer region of finite water depth is derived for the second-order diffraction potential, which gives a clear physical meaning of the wave transmission and reflection characteristics in the far field. A numerical method-simple Green's function technique-for calculating the second-order diffraction potential in the inner region is also described. Numerical results are provided for the second-order wave forces on a semi-submerged cylinder. It is found that the contribution of second-order diffraction potential to second-order wave forces is important. The effect of water depth and submerged depth on the wave force is also discussed. 相似文献
16.
Numerical and experimental study on hydrodynamic interaction of side-by-side moored multiple vessels 总被引：2，自引：0，他引：2
This paper aims to investigate the basic interaction characteristics of side-by-side moored vessels both numerically and experimentally. A higher-order boundary element method (HOBEM) combined with generalized mode approach is applied to analysis of motion and drift force of side-by-side moored multiple vessels (LNG FPSO, LNGC and shuttle tankers). Model tests were carried out for the same floating bodies investigated in the numerical study in regular and irregular waves. Global and local motion responses and drift forces of three vessels are compared with those of calculations. Discussions is highlighted on applicability of numerical method to prediction of sophisticated multi-body interaction problem of which motion behavior is very important to analysis of mooring dynamics of deep sea floating bodies. 相似文献
17.
A time-domain higher-order boundary element scheme was utilized to simulate the linear and nonlinear irregular waves and diffractions due to a structure. Upon the second-order irregular waves with four Airy wave components being fed through the inflow boundary, the fully nonlinear boundary problem was solved in a time-marching scheme. The open boundary was modeled by combining an absorbing beach and the stretching technique. The proposed numerical scheme was verified by simulating the linear regular and irregular waves. The scheme was further applied to compute the linear and nonlinear irregular wave diffraction forces acting on a vertical truncated circular cylinder. The nonlinear results were also verified by checking the accuracy of the nonlinear simulation. 相似文献
18.
采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究，对自由表面边界条件进行Taylor级数展开，应用摄动展开可以建立相应的边值问题，而且此边值问题的计算域不随时间变化，运用基于B－样条的边界元方法求解每一时刻的波浪场，二阶自由表面边界条件在时间上进行数值积分，在自由表面加了一个人工阻尼层以避免波浪的反射，速度势分解为已知的入射势和未知的散射势，初始条件采用二阶Stokes波浪场，通过加入物体表面边界条件，得到散射势在时间和空间上的发展，本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算，数值结果表明此理论计算准确，效率高，数值稳定。 相似文献
19.
The mooring of offshore floating structures, such as offshore platforms, in large waves against drift forces and rotational moments is a challenging problem in offshore engineering. To accurately investigate such problems, called positioning problems, the time-averaged steady forces of the second order known as the wave drift forces must be taken into account. Fortunately, a cloaking phenomenon occurs under certain conditions and dramatically reduces the wave drift force acting on such a floating body, as previously reported by several researchers. In the diffraction problem of water waves, cloaking refers to the condition where there is no scattering in the form of radial outgoing waves. The reduction of wave drift force on a truncated cylinder with the occurrence of cloaking phenomenon has been numerically and experimentally confirmed. In this paper, the arrangement of several small circular cylinders at regular intervals in a circle concentric with a fixed floating body is considered as an effective means of reducing the wave drift force. Using a combination of a higher-order boundary element method (HOBEM) and wave interaction theory, the influences of the geometric parameters of the outer surrounding cylinders on the wave drift force and the total scattered-wave energy are systematically investigated and discussed. A quasi-cloaking phenomenon is first found and reported in the present study, which is beneficial and flexible for application in practical engineering. More than one quasi-cloaking trigger (where a trigger is an occurrence condition) can be found simply by varying the distance between the inner and outer floating bodies. 相似文献
20.
Higher Harmonics Induced by Waves Propagating over A Submerged Obstacle in the Presence of Uniform Current 总被引：1，自引：0，他引：1
To investigate higher harmonics induced by a submerged obstacle in the presence of uniform current, a 2D fully nonlinear numerical wave flume（NWF） is developed by use of a time-domain higher-order boundary element method（HOBEM） based on potential flow theory. A four-point method is developed to decompose higher bound and free harmonic waves propagating upstream and downstream around the obstacle. The model predictions are in good agreement with the experimental data for free harmonics induced by a submerged horizontal cylinder in the absence of currents. This serves as a benchmark to reveal the current effects on higher harmonic waves. The peak value of non-dimensional second free harmonic amplitude is shifted upstream for the opposing current relative to that for zero current with the variation of current-free incident wave amplitude, and it is vice versa for the following current. The second-order analysis shows a resonant behavior which is related to the ratio of the cylinder diameter to the second bound mode wavelength over the cylinder. The second-order resonant position slightly downshifted for the opposing current and upshifted for the following current. 相似文献