波浪作用下海床动力反应的数值分析

近海和离岸建筑物的基通常处于连续不断的小风浪作用之下，可将其变形视为弹性。基于二维广义Biot理论，提出了线性或非线性波浪作用下饱和弹性海床动力应应的时域有限元数值解法。静力平衡条件和Biot方程组成的边值方程可视为其特例，在比较算例中，数值计算得到的孔压和有效应力幅值沿海床深度的分布与解析解十分吻合。土骨架和孔隙流体的加速度对海床动力反应的影响很小。具体算例表明，线性波沿缓坡海床传播时，土层中超静孔压和有效应力幅值随之增大，有可能发生滑动坡坏。所提出的数值解法能够灵活地处理非线性波浪荷载，海底复杂地质条件和波浪沿缓坡传播等复杂情况。     

Modified mild slope equation and open boundary conditions

A numerical model to compute wave field is developed. It is based on the Berkhoff diffraction-refraction equation, in which an energy dissipation term is added, to take into account the breaking and the bottom friction phenomena. The energy dissipation function, by breaking and by bottom friction, is introduced in the Berkhoff equation to obtain a new equation of propagation.The resolution is done with the hybrid finite element method, where lagrangians elements are used.     

Stress and deformation of offshore piles under structural and wave loading

Various offshore structures, especially large structures such as Tension Leg Platforms (TLP), are usually supported by concrete piles as the foundation elements. The stress distribution within such a large structure is a dominant factor in the design procedure of an offshore pile. To provide a more accurate and effective design for offshore foundation systems under axial and lateral wave loads, a finite element model is employed herein to determine the stresses and displacements in a concrete pile under similar loading conditions. A parametric study is also performed to examine the effects of the stress distribution due to the changing loading conditions.     

Uncertainty in predictions of oil spill trajectories in open sea

A method is introduced to calculate and to account for the uncertainties in the predictions of oil spill trajectories using a classic oil spill model. The method considers the output of the oil spill model as a function of random variables, which are the input parameters, and calculates the standard deviation of the output results which gives a measure of the uncertainty of the model given the uncertainties of the input parameters.Instead of a single trajectory that is calculated by the oil spill model using the mean values of the parameters, a band of trajectories can be defined when various simulations are done taking into account the uncertainties of the input parameters. This band of trajectories defines envelopes of the trajectories that are likely to be followed by the spill given the uncertainties of the input.The method is applied to an oil spill that occurred in open sea near Madeira Islands, in the Atlantic Ocean, in December 1989. The simulations allow the understanding of how a change in the wind direction drove the spill towards the Islands.The envelope of likely trajectories that is obtained with the uncertainty modelling shows a band of trajectories that is in better agreement with the observations than the single trajectory simulated by the oil spill model, based on mean parameters.     

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

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.     

Boussinesq-type modelling using an unstructured finite element technique

A model for solving the two-dimensional enhanced Boussinesq equations is presented. The model equations are discretised in space using an unstructured finite element technique. The standard Galerkin method with mixed interpolation is applied. The time discretisation is performed using an explicit three-step Taylor–Galerkin method. The model is extended to the surf and swash zone by inclusion of wave breaking and a moving boundary at the shoreline. Breaking is treated by an existing surface roller model, but a new procedure for the detection of the roller thickness is devised. The model is verified using four test cases and the results are compared with experimental data and results from an existing finite difference Boussinesq model.     

Numerical study of the meridional overturning circulation with “mixing hotspots” in the Pacific Ocean

Using an idealized ocean general circulation model, we examine the effect of “mixing hotspots” (localized regions of intense diapycnal mixing) predicted based on internal wave-wave interaction theory (Hibiya et al., 2006) on the meridional overturning circulation of the Pacific Ocean. Although the assumed diapycnal diffusivity in the mixing hotspots is a little larger than the predicted value, the upwelling in the mixing hotspots is not sufficient to balance the deep-water production; out of 17 Sv of the downwelled water along the southern boundary, only 9.2 Sv is found to upwell in the mixing hotspots. The imbalance as much as 7.8 Sv is compensated by entrainment into the surface mixed layer in the vicinity of the downwelling region. As a result, the northward transport of the deep water crossing the equator is limited to 5.5 Sv, much less than estimated from previous current meter moorings and hydrographic surveys. One plausible explanation for this is that the magnitude of the meridional overturning circulation of the Pacific Ocean has been overestimated by these observations. We raise doubts about the validity of the previous ocean general circulation models where diapycnal diffusivity is assigned ad hoc to attain the current magnitude suggested from current meter moorings and hydrographic surveys.     

The Structure of A Turbulent Jet in A Crossflow-Effect of Jet-Crossflow Velocity

A comprehensive numerical study on the three-dimensional structure of a turbulent jet in crossflow is performed. The jet-to-crossflow velocity ratio (R) varies in the range of 2 - 16; both vertical jets and inclined jets without excess streamwise momentum are considered. The numerical results of the Standard two-equation k-ε model show that the turbulent structure can be broadly categorised according to the jet-to-crossflow velocity ratio. For strong to moderate jet discharges, i.e. R> 4, the jet is characterized by a longitudinal transition through a bent-over phase during which the jet becomes almost parallel with the main freestream, to a sectional vortex-pair flow with double concentration maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless distance of around 20-60. The similarity coefficients are only weakly dependent on R. The cross-section scalar field is kidney-shaped and bifurcated, vvith distinct double concentr