利用波浪数模试验优化设计方案

介绍葫芦岛港一期工程波浪数模试验结果及利用波浪数模试验结果进行的设计优化

福建沿海海域波浪能资源分析与评价

采用波浪模拟的方法,较准确计算得出福建沿海海域波浪能资源分布状况,并给出相应的分析和综合评价.主要结论如下:(1)福建沿海海域波浪能平均密度为2.6～7.3 kW/m,波浪能资源储量为2 210.45 MW,在我国沿海海域仅次于台湾和广东,是波浪能开发利用可以优先考虑的海区之一.(2)福建沿海海域波浪能资源储量的70%分布于平潭岛以北海域,其值达1 512.49 MW.其中,尤以北礵地区值最大,为378.80 MW.(3)以年平均波高为指标,福建沿海海域中东山区段为三类区,其他区段均为一类区和二类区,具有良好的开发前景.(4)福建沿海海域波浪能具有波功率密度低、资源分布广泛且不均匀、波功率密度随季节变化、能量具有多向性等分布特点.(5)基于福建波浪能的开发与利用现状,建议应优先着眼于解决边远海岛等特殊场所的用电问题.

基于FHDI-GNWM的某海域波浪时空分布特性分析

采用2015—2019年的FHDI-GNWM后报波浪资料,并结合波浪数学模型,对某海域的波浪分布特征进行了初步统计分析,得到如下结论:1)海区波向与风向基本一致,冬季主浪向为E向;夏季主浪向为SSW向;2)海区波浪以风浪为主,且表现出明显的季风特征,夏季波高较小,平均波高介于0.4~0.7 m;冬季波高较大,平均波高介于0.85~1.1 m。此外,基于某集装箱码头项目,结合船舶系泊试验,分析了波浪引起船舶装卸不可作业时间。文章成果对认识某近岸海域的风、浪特征以及对相关码头设计具有一定参考价值。

A VOF-based numerical model for breaking waves in surf zone

This paper introduces a numerical model for studying the evolution of a periodic wave train, shoaling, and breaking in surf zone. The model can solve the Reynolds averaged Navier-Stokes （RANS） equations for a mean flow, and the k-ε equations for turbulence kinetic energy k and turbulence dissipation rate ε. To track a free surface, the volume of fluid （VOF） function, satisfying the advection equation was introduced. In the numerical treatment, third-order upwind difference scheme was applied to the convection terms of the RANS equations in order to reduce the effect of numerical viscosity. The shoaling and breaking processes of a periodic wave train on gently sloping beaches were modeled. The computed wave heights of a sloping beach and the distribution of breaking wave pressure on a vertical wall were compared with laboratory data.

Simulation of double cold cores of the 35°N section in the YellowSea with a wave-tide-circulation coupled model

Based on the MASNUM wave-tide-circulation coupled numerical model, the temperature structure along 35°N in the Yellow Sea was simulated and compared with the observations. One of the notable features of the temperature structure along 35°N section is the double cold cores phenomena during spring and summer. The double cold cores refer to the two cold water centers located near 122°E and 125°E from the depth of 30m to bottom. The formation, maintenance and disappearance of the double cold cores are discussed. At least two reasons make the temperature in the center (near 123°E) of the section higher than that near the west and east shores in winter. One reason is that the water there is deeper than the west and east sides so its heat content is higher. The other is invasion of the warm water brought by the Yellow Sea Warm Current (YSWC) during winter.This temperature pattern of the lower layer (from 30m to bottom) is maintained through spring and summer when the upper layer (0 to 30m) is heated and strong thermocline is formed. Large zonal span of the 35°N section (about 600 km) makes the cold cores have more opportunity to survive. The double cold cores phenomena disappears in early autumn when the west cold core vanishes first with the dropping of the thermocline position.

Tidal effects on temperature iront in the Yellow Sea

Temperature front (TF) is one of the important features in the Yellow Sea, which forms in spring,thrives in summer, and fades in autumn as thermocline declines. TF intensity |ST| is defined to describe the distribution of TF. Based on the MASNUM wave-tide-circulation coupled model, temperature distribution in the Yellow Sea was simulated with and without tidal effects. Along 36°N, distribution of TF from the simulated results are compared with the observations, and a quantitative analysis is introduced to evaluate the tidal effects on the forming and maintaining processes of the TF. Tidal mixing and the circulation structure adapting to it are the main causes of the TF.

Hydrodynamic Characteristics of Three Rows of Vertical Slotted Wall Breakwaters

In this study, we examine the hydrodynamic characteristics of three rows of vertical slotted wall breakwaters in which the front and middle walls are permeable and partially immersed in a water channel of constant depth, whereas the third wall is impermeable. The wave–structure interaction and flow behavior of this type of breakwater arrangement are complicated and must be analyzed before breakwaters can be appropriately designed. To study the hydrodynamic breakwater performance, we developed a mathematical model based on the eigenfunction expansion method and a least squares technique for predicting wave interaction with three rows of vertical slotted wall breakwaters. We theoretically examined the wave transmission, reflection, energy loss, wave runup, and wave force under normal regular waves. Comparisons with experimental measurements show that the mathematical model results adequately reproduce most of the important features. The results of this investigation provide a better understanding of the hydrodynamic performance of triple-row vertical slotted wall breakwaters.

Wave-induced mixing in the Yellow Sea

Vertical wave-induced mixing parameter Bv expressed in wave number spectrum was estimated in the Yellow Sea. The spatial distributions of Bv averaged over upper 20 m in 4 seasons were analyzed. It is the strongest in winter because of winter monsoon, and the weakest in spring. Since in summer it plays an important role for circulation of upper layers, its vertical structure was also discussed. Two simulations with and without wave-induced mixing in this season were performed to evaluate its effect on temperature distribution. Numerical results indicate that wave-induced mixing could increase the mixed layer thickness greatly.

Numerical Simulation of Nonlinear Wave Force on a Quasi-ellipse Caisson

A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid （VOF） method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force.

Comparison of TOPEX/POSEIDON altimeter derived wave period with ocean buoy data in the East China Sea and South China Sea

Altimeter wave period data obtained from continental shelf seas are analyzed in this paper. Empirical models are introduced for zero up-crossing and peak wave period calculation with TOPEX/POSEIDON data. Their performances are assessed using independent validation dataset in four sites in the open ocean of China. To provide more accurate wave period estimation, new coefficients are applied to reliable in situ data. Comparison of our estimated the wave periods with new linear calibrations based on independent data of Seapac 2100 deployed in the East China Sea and South China Sea showed that the accuracy was improved over estimates determined from earlier empirical models. Regional analysis indicated that the wave period model works better under wind sea condition.

RANSE Simulation of High-speed Planning Craft in Regular Waves

This paper presents a study on the numerical simulation of planing crafts sailing in regular waves. This allows an accurate estimate of the seas keeping performance of the high speed craft. The simulation set in six-degree of freedom motions is based on the Reynolds averaged Navier Stokes equations volume of fluid （RANSE VOF） solver. The trimming mesh technique and integral dynamic mesh method are used to guarantee the good accuracy of the hydrodynamic force and high efficiency of the numerical simulation. Incident head waves, oblique waves and beam waves are generated in the simulation with three different velocities （Fn =1.0, 1.5, 2.0）. The motions and sea keeping performance of the planing craft with waves coming from different directions are indicated in the flow solver. The ship designer placed an emphasis on the effects of waves on sailing amplitude and pressure distribution of planing craft in the configuration of building high speed crafts.

Numerical simulation of scatterometer assimilated wind and ocean wave in eastern China seas and adjacent waters

Using the latest version of Mesoscale Modeling System （MMSv3）, we assimilated wind data from the scatterometer and built a model to assimilate the wind field over eastern China seas and adjacent waters and applied the wave model WAVEWATCH-Ⅲ to test the sea area with assimilative wind and blended wind of QSCAT and NCEP as driving forces. High precision and resolution numerical wave results were obtained. Analysis indicated that if we replace the model wind result with the blended wind, better sea surface wind results and wave results could be obtained.

Study on Unsteady Hydrodynamic Performance of Propeller in Waves

The speed of a ship sailing in waves always slows down due to the decrease in efficiency of the propeller. So it is necessary and essential to analyze the unsteady hydrodynamic performance of propeller in waves. This paper is based on the numerical simulation and experimental research of hydrodynamics performance when the propeller is under wave conditions. Open-water propeller performance in calm water is calculated by commercial codes and the results are compared to experimental values to evaluate the accuracy of the numerical simulation method. The first-order Volume of Fluid(VOF) wave method in STAR CCM+ is utilized to simulate the three-dimensional numerical wave. According to the above prerequisite, the numerical calculation of hydrodynamic performance of the propeller under wave conditions is conducted, and the results reveal that both thrust and torque of the propeller under wave conditions reveal intense unsteady behavior. With the periodic variation of waves, ventilation, and even an effluent phenomenon appears on the propeller. Calculation results indicate, when ventilation or effluent appears, the numerical calculation model can capture the dynamic characteristics of the propeller accurately, thus providing a significant theory foundation forfurther studying the hydrodynamic performance of a propeller in waves.

3D Numerical Modeling of Wave Forces on Tandem Fixed Cylinders Using the BEM

In this paper a 3D numerical model was developed to study the complicated interaction between waves and a set of tandem fixed cylinders.The fluid was considered to be inviscid and irrotational.Therefore,the Helmholtz equation was used as a governing equation.The boundary element method（BEM） was adopted to discretize the relevant equations.Open boundaries were used in far fields of the study domain.Linear waves were generated and propagated towards tandem fixed cylinders to estimate the forces applied on them.Special attention was paid to consideration of the effect on varying non-dimensional cylinder radius and distance between cylinders,ka and kd on forces and trapped modes.The middle cylinder wave forces and trapped modes in a set of nine tandem cylinders were validated utilizing analytical data.The comparisons confirm the accuracy of the model.The results of the inline wave force estimation on n tandem cylinders show that the critical cylinder in the row is the middle one for odd numbers of cylinders.Furthermore the results show that the critical trapped mode effect occurs for normalized cylinder radiuses close to 0.5 and 1.0.Finally the force estimation for n tandem cylinders confirms that force amplitude of the middle cylinder versus normalized separation distance fluctuates about that of a single cylinder.

Study of algorithms of new slender ship theory of wave resistance

In this paper, Noblesse's New Slender-Ship Wave-Making Theory was investigated numerically. Detailed expressions of zeroth and lst order wave resistance have been derived and calculation programs have also been compiled. In the single and double integral terms of Green function, the kernel function of wave resistance expression, special function expansion method and Chebyshev polynomials approach have been adopted respectively, which greatly simplify the calculation and increase the convergence speed.

Analytical Solution of Boundary Integral Equations for 2-D Steady Linear Wave Problems

Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed,free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigensolution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore,an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.

Numerical Simulation of Water Wave Propagation and Transformation

A numerical approach was performed to predict the propagation and transformation of nonlinear water waves. A numerical wave flume was developed based on the non-periodic high-order spectral (HOS) method. The flume was applied to analyze the effect of wave steepness and wavelength on the propagation of nonlinear waves. The results show that for waves of low steepness, the wave profile and the wave energy spectrum are stable, and that the propagation can be predicted by the linear wave theory. For waves of moderate steepness and steep waves, the effects associated with the interactions between waves in a wave group become significant and a train of initially sinusoidal waves may drastically change its form within a short distance from its original position.

Numerical Simulation of Motion Response of an Offshore Observation Platform in Waves

Offshore observation platforms are required to have great ability to resist waves when they are operating at sea. Investigation on the motion characteristics of the platforms in the sea can provide significant reference values during the platform design procedure. In this paper, a series of numerical simulation on the interaction of a triple-hulled offshore observation platform with different incident waves is carried out. All of the simulations are implemented utilizing our own solver naoe-FOAM-SJTU, which is based and developed on the open source tools of OpenFOAM. Duration curves of motion characteristics and loads acting on the platform are obtained, and a comparison between the results of the amplitude in different incident waves is presented. The results show that the solver is competent in the simulation of motion response of platforms in waves.

Numerical Simulation of Multi-Directional Random Wave Transformation in a Yacht Port

This paper extends a prediction model for multi-directional random wave transformation based on an energy balance equation by Mase with the consideration of wave shoaling, refraction, diffraction, reflection and breaking. This numerical model is improved by 1) introducing Wen's frequency spectrum and Mitsuyasu's directional function, which are more suitable to the coastal area of China; 2) considering energy dissipation caused by bottom friction, which ensures more accurate results for large-scale and shallow water areas; 3) taking into account a non-linear dispersion relation. Predictions using the extended wave model are carried out to study the feasibility of constructing the Ai Hua yacht port in Qingdao, China, with a comparison between two port layouts in design. Wave fields inside the port for different incident wave directions, water levels and return periods are simulated, and then two kinds of parameters are calculated to evaluate the wave conditions for the two layouts. Analyses show that Layout I is better than Layout II. Calculation results also show that the harbor will be calm for different wave directions under the design water level. On the contrary, the wave conditions do not wholly meet the requirements of a yacht port for ship berthing under the extreme water level. For safety consideration, the elevation of the breakwater might need to be properly increased to prevent wave overtopping under such water level. The extended numerical simulation model may provide an effective approach to computing wave heights in a harbor.

Study on the Cause and Dynamics of Seawater Intrusion in Laizhou Bay

In this paper, we conduct research on the cause and dynamics of seawater intrusion in the Laizhou bay. To consider modelling the Laizhou bay, we should then take the listed tools into consideration. In actual engineering, wave is a complicated stochastic process. If the numerical value of wave model is a little too large, the wave attenuation should be considered sufficiently, such as wave refraction caused by terrain change, wave attenuation due to submarine friction, wave into the very shallow waters of the broken and small wind area. We integrate the current conition of the Laizhou bay and the further model the seawater intrusion steps. The uncertainty research of ground water solute migration is at present a ground water research area quite popular topic, but the seawater intrusion is the variable density class solute migration issue, applies the uncertainty analysis of seawater intrusion simulation also to have certain difficulty the existing research results. In the future, more simulation will be down to verify the effectiveness.