首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 296 毫秒
1.
A three-dimensional (3D) numerical model of fixed Oscillating Water Column system (OWC) is presented and validated. The steady-state potential flow boundary value problem due to regular wave interaction with the OWC is solved by a first order mixed distribution panel method. Ocean response predictions are derived using a deterministic statistical model based on a spectral analysis method. The model validation focusses on diffraction predictions and involves convergence tests and numerical comparisons with independent potential flow computations. Predictions of both regular and irregular wave responses are also compared against experimental results. Sample results including the yearly-averaged power conversion efficiency are presented in the final section to illustrate the method’s suitability to a 3D hydrodynamic design optimisation.  相似文献   

2.
A surface panel method is employed for the thin boundary layer calculation of heavily loaded marine propellers in steady state conditions. Employing the surface panel method, known as the “Morino Method”, the flow field around the propeller is represented by an unknown potential. The majority of the flow field is governed by the potential theory while the viscosity is assumed to be largely confined to thin shear layer on the propeller surface. The boundary layer calculations are performed by using Cebeci-Smith two dimensional model and the local skin friction coefficients and blowing velocities are obtained along the pre-computed on-body streamlines. It is shown that the prediction of torque of the propeller is improved when the boundary layer calculations are used instead of the boundary layer corrections based on the formulae established for the flat plates.  相似文献   

3.
A transformation method is presented by which current profiles (of tidal or wind-induced origin) can be extracted at any horizontal position and moment in time from a vertically integrated, two-dimensional, hydrodynamic numerical model. An arbitrary vertical variation of eddy viscosity can be included in the method, which can incorporate a no-slip bottom boundary condition. The technique assumes that the sea is homogeneous.The method is used to improve the representation of bottom stress within the two-dimensional model, whereby the bottom stress is no longer related simply to the depth-mean current as in the “conventional” two-dimensional, vertically integrated model.Idealized calculations for a range of eddy viscosity profiles, show that elevations, current profiles, and time series of current extracted from this “enhanced” two-dimensional numerical model are in good agreement with currents obtained from a full three-dimensional model.  相似文献   

4.
5.
Nan Xie  Dracos Vassalos   《Ocean Engineering》2007,34(8-9):1257-1264
The purpose of the present paper is to develop a potential-based panel method for determining the steady potential flow about three-dimensional hydrofoil under free surface. The method uses constant-strength doublets and source density distribution over the foil body surface and thereby Dirichlet-type boundary condition is used instead of Neumann-type condition. On the undisturbed free surface source density is used to meet the free surface condition that is linearised in terms of double-body model approach and is discretised by a one-side, upstream, four-point finite difference operator. After solving the doublets on the foil and sources on the free surface, the numerical results of pressure, lift and resistance coefficients and also wave profiles can then be calculated for different Froude number and depth of submergence to demonstrate the influence of free surface and aspect ratio effects on performance of the hydrofoil.  相似文献   

6.
ZHAO  Ming 《中国海洋工程》2002,16(4):513-523
A numerical model is developed for estimation of local scour around a large circular cylinder under vvave action. The model includes wave diffraction around structures, bed shear stress calculation inside the vvave boundary layer and topo-graphical change model. The vvave model is based on the improved Boussinesq equations for varying depth. The vvave boundary layer is calculaled by solving the integrated momentum equation over the boundary layer. The bed shear stress due to streaming, an important factor affecting the sediment transport around a large-scale cylinder, is calculated. The Lagrangian drift velocity is included in calculation of the suspended sediment transport rates. The model is implemented by a finite element method and the results from the present model, which agree well with experimental data, are com-pared vvith those from other methods.  相似文献   

7.
X.T. Zhang  B.C. Khoo  J. Lou 《Ocean Engineering》2006,33(17-18):2310-2331
The problem of wave propagation in a fully nonlinear numerical wave tank is studied using desingularized boundary integral equation method coupled with mixed Eulerian–Lagrangian formulation. The present method is employed to solve the potential flow boundary value problem at each time step. The fourth-order predictor–corrector Adams–Bashforth–Moulton scheme is used for the time-stepping integration of the free surface boundary conditions. A damping layer near the end-wall of wave tank is added to absorb the outgoing waves with as little wave reflection back into the wave tank as possible. The saw-tooth instability is overcome via a five-point Chebyshev smoothing scheme. The model is applied to several wave propagations including solitary, irregular and random incident waves.  相似文献   

8.
The objective of the present work is to discuss the implementation of an active wave generating–absorbing boundary condition for a numerical model based on the Volume Of Fluid (VOF) method for tracking free surfaces. First an overview of the development of VOF type models with special emphasis in the field of coastal engineering is given. A new type of numerical boundary condition for combined wave generation and absorption in the numerical model VOFbreak2 is presented. The numerical boundary condition is based on an active wave absorption system that was first developed in the context of physical wave flume experiments, using a wave paddle. The method applies to regular and irregular waves. Velocities are measured at one location inside the computational domain. The reflected wave train is separated from the incident wave field in front of a structure by means of digital filtering and subsequent superposition of the measured velocity signals. The incident wave signal is corrected, so that the reflected wave is effectively absorbed at the boundary. The digital filters are derived theoretically and their practical design is discussed. The practical use of this numerical boundary condition is compared to the use of the absorption system in a physical wave flume. The effectiveness of the active wave generating–absorbing boundary condition finally is proved using analytical tests and numerical simulations with VOFbreak2.  相似文献   

9.
Sheet flow and suspension of sand in oscillatory boundary layers   总被引:1,自引:0,他引:1  
after revisionTime-dependent measurements of flow velocities and sediment concentrations were conducted in a large oscillating water tunnel. The measurements were aimed at the flow and sediment dynamics in and above an oscillatory boundary layer in plane bed and sheet-flow conditions. Two asymmetric waves and one sinusoidal wave were imposed using quartz sand with D50 = 0.21 mm. A new electro-resistance probe with a large resolving power was developed for the measurement of the large sediment concentrations in the sheet-flow layer. The measurements revealed a three layer transport system consisting of a pick-up/deposition layer, an upper sheet flow layer and a suspension layer.In the asymmetric wave cases the total net transport was directed “onshore” and was mainly concentrated in the thin sheet flow layer (< 0.5 cm) at the bed. A small net sediment flux was directed “offhore” in the upper suspension layer. The measured flow velocities, sediment concentrations and sedimenl fluxes showed a good qualitative agreement with the results of a (numerical) 1DV boundary-layer flow and transport model. Although the model did not describe all the observed processes in the sheet-flow and suspension layer, the computational results showed a reasonable agreement with measured net transport rates in a wide range of asymmetric wave conditions.  相似文献   

10.
A numerical solution is developed to investigate the generation and propagation of small-amplitude water waves in a semi-infinite rectangular wave basin. The three-dimensional wave field is produced by the prescribed “snake-like” motion of an array of segmented wave generators located along the wall at one end of the tank. The solution technique is based on the boundary element approach and uses an appropriate three-dimensional Green function which explicitly satisfies the tank-wall boundary conditions. The Green function and its derivatives which appear in the integral equation formulation can be shown to be slowly convergent when the source and field points are in close proximity. Therefore, when computing the velocity potentials on the wave generators, the source points are chosen outside the fluid domain, thereby ensuring the rapid convergence of these functions and rendering the integral equations non-singular. Numerical results are shown which illustrate the influence of the various wavemaker and basin parameters on the generated wave field. Finally, the complete wave field produced by the diffraction of oblique waves by a vertical circular cylinder in a basin is presented.  相似文献   

11.
X.T. Zhang  B.C. Khoo  J. Lou 《Ocean Engineering》2007,34(10):1449-1458
A numerical approach based on desingularized boundary element method and mixed Eulerian–Lagrangian formulation [Zhang et al., 2006. Wave propagation in a fully nonlinear numerical wave tank: a desingularized method. Ocean Engineering 33, 2310–2331] is extended to solve the water wave propagation over arbitrary topography in a three-dimensional wave tank. A robust damping layer applicable for regular and irregular incident waves is employed to minimize the outgoing wave reflection back into the wave tank. Numerical results on the propagation of regular and irregular incident waves over the flat bottom and linear incident waves over an elliptical shoal show good concurrence with the corresponding analytical solutions and experimental data.  相似文献   

12.
The hybrid Cartesian/immersed boundary method is applied to simulate effects of flexibility on propulsive force acting on a heaving foil in a viscous flow. Immersed boundary nodes are distributed inside an instantaneous fluid domain. Velocity vector is reconstructed at the immersed boundary node based on an interpolation along a local normal line. Using the staggered/non-staggered grid method, the demand for pressure at boundary nodes is removed. Elastic deformation of the flexible foil is modelled based on the dynamic thin-plate mechanics. The developed code is validated through comparisons with other computations on flow fields around a flapping foil. The generation of the reverse Karman vortex street is investigated. Forces acting on heaving foils are compared for flexible and rigid cases and the increased thrust of the flexible foil is attributed to the deformed configuration near the tip. The flexibility of the heaving foil decreases vertical force and improves propulsion efficiency. The variations of force and deformation are investigated according to bending stiffness of the foil.  相似文献   

13.
14.
The seakeeping characteristics of various boom geometries in irregular waves and currents are investigated. The response of a floating boom section on the open sea is a function of a number of parameters, such as boom geometry, distribution of mass, buoyancy/weight ratio, and wave and current characteristics. To understand the relationship between these design parameters more clearly, a series of regular and irregular wave tests were conducted with six different 1:4 scale models for three current velocities and six different wave conditions. To simplify the problem, only rigid boom sections consisting of a buoyancy cylinder and vertical skirt were used. In parallel with this experimental program, a numerical model for the responses of two-dimensional floating boom sections in small-amplitude waves is also developed. The numerical results are compared with our large-scale experimental results. The boom effectiveness on the open sea is evaluated based on the concept of “effective draft” and “effective freeboard” assuming that drainage and oversplashing failures are the prime mechanisms of containment failure. Using the present results, a guideline for the optimum design/selection of future booms is developed.  相似文献   

15.
柱体绕流问题是流体力学领域一个非常经典的问题。当流体流经柱体时,由于黏性的存在,会发生许多复杂的流动现象,如流动分离、涡旋周期性生成与脱落等,经常被作为标准验证算例。同时,柱体绕流广泛存在于实际工程中,并在一定工况下可能对工程产生巨大危害,因此对柱体绕流进行深入研究具有重要意义。研究中,拟将一种无网格类方法——半隐式移动粒子方法(moving particle semi-implicit method,简称MPS)引入到柱体绕流问题的数值研究中,并对不同雷诺数下二维方柱绕流问题进行数值模拟。首先,使用基于MPS方法自主开发的MLParticle-SJTU求解器,结合入口边界条件和出口边界条件,模拟了雷诺数Re分别为40、200和1 000时均匀来流条件下的方柱绕流。随后,将模拟的绕流结果与文献中试验和数值计算结果进行了对比,结果吻合较好,并且在雷诺数为200和1 000时,可以清晰地捕捉到方柱尾流中的卡门涡街现象,验证了MPS方法在柱体绕流问题模拟上的有效性和适用性。  相似文献   

16.
A three-dimensional numerical model is established to simulate the turbulent oscillatory boundary layer over a fixed and rough bed composed by randomly arrayed solid spheres based on the lattice Boltzmann method and the large eddy simulation model.The equivalent roughness height,the location of the theoretical bed and the time variation of the friction velocity are investigated using the log-fit method.The time series of turbulent intensity and Reynolds stress are also investigated.The equivalent roughness height of cases with Reynolds numbers of 1×10~4–6×10~4 is approximately 2.81 d(grain size).The time variation of the friction velocity in an oscillatory cycle exhibits sinusoidal-like behavior.The friction factor depends on the relative roughness in the rough turbulent regime,and the pattern of solid particles arrayed as the rough bed in the numerical simulations has no obvious effect on the friction factor.  相似文献   

17.
A numerical model is developed to predict the onset of local scour below offshore pipelines in steady currents and waves. The scour is assumed to start when the pressure gradient underneath the pipeline exceeds the floatation gradient of the sediments. In this model, the water flow field above the bed is determined by solving the two-dimensional (2-D) Reynolds-averaged Navier–Stokes equations with a k-ω turbulence closure. The seepage flow below the seabed is calculated by solving the Darcy's law (Laplace's equation) with known pressure distribution along the common boundaries of the flow domains-seabed. The numerical method used for both the turbulent flow around the pipeline and Darcy's flow in the seabed is a fractional finite element method. The average pressure gradient along the buried pipe surface is employed in the evaluation of onset condition with a calibration coefficient. The numerical model is validated against experimental data available in literature. A unified onset condition for steady currents and waves is proposed. Influences of flow parameters, including water depth, embedment depth, boundary layer thickness, Reynolds number (Re) and Keuleagan–Carpenter (KC) number, on the pressure drop coefficient over the pipeline are studied systematically.  相似文献   

18.
A large number of studies have been done dealing with sinusoidal wave boundary layers in the past. However, ocean waves often have a strong asymmetric shape especially in shallow water, and net of sediment movement occurs. It is envisaged that bottom shear stress and sediment transport behaviors influenced by the effect of asymmetry are different from those in sinusoidal waves. Characteristics of the turbulent boundary layer under breaking waves (saw-tooth) are investigated and described through both laboratory and numerical experiments. A new calculation method for bottom shear stress based on velocity and acceleration terms, theoretical phase difference, φ and the acceleration coefficient, ac expressing the wave skew-ness effect for saw-tooth waves is proposed. The acceleration coefficient was determined empirically from both experimental and baseline kω model results. The new calculation has shown better agreement with the experimental data along a wave cycle for all saw-tooth wave cases compared by other existing methods. It was further applied into sediment transport rate calculation induced by skew waves. Sediment transport rate was formulated by using the existing sheet flow sediment transport rate data under skew waves by Watanabe and Sato [Watanabe, A. and Sato, S., 2004. A sheet-flow transport rate formula for asymmetric, forward-leaning waves and currents. Proc. of 29th ICCE, ASCE, pp. 1703–1714.]. Moreover, the characteristics of the net sediment transport were also examined and a good agreement between the proposed method and experimental data has been found.  相似文献   

19.
Pressure variations and three-dimensional effects on liquid sloshing loads in a moving partially filled rectangular tank have been carried out numerically and experimentally. A numerical algorithm based on the volume of fluid (VOF) technique is used to study the non-linear behavior and damping characteristics of liquid sloshing. A moving coordinate system is used to include the non-linearity and avoid the complex boundary conditions of moving walls. The numerical model solves the complete Navier–Stokes equations in primitive variables by using of the finite difference approximations. In order to mitigate a series of discrete impacts, the signal computed is averaged over several time steps. In order to assess the accuracy of the method used, computations are compared with the experimental results. Several configurations of both baffled and unbaffled tanks are studied. Comparisons show good agreement for both impact and non- impact type slosh loads in the cases investigated.  相似文献   

20.
A two-layered model is considered in which the upper layer is continuously stratified and the lower layer is homogeneous. The system is driven by atmospheric forces. Bottom stress and topography are included in the model. The linear three-dimensional hydrodynamical equations are used to describe the system. Taking the eddy viscosity in the upper layer as inversely proportional to the static stability, the dependent variables are expanded in terms of continuous functions in the vertical (eigenfunctions). Using this method it is possible to compute currents and internal displacements at any depth in the upper layer. The three-dimensional structure of the lower layer is not considered in this model. The equations describing the lower layer are integrated over depth to give depth mean currents. Using a staggered finite-difference grid in the horizontal and a forward time-stepping procedure, numerical test experiments are carried out for a cross section and for a closed rectangular basin.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号