Computation of solitary waves during propagation and runup on a slope

A numerical time-simulation algorithm for analysing highly nonlinear solitary waves interacting with plane gentle and steep slopes is described by employing a mixed Eulerian–Lagrangian method. The full nonlinear free surface conditions are considered here in a Lagrangian frame of reference without any analytical approximations, and thus the method is valid for very steep waves including overturning. It is found that the runup height is crucially dependent on the wave steepness and the slope of the plane. Pressures and forces exerted on impermeable walls of different inclinations (slopes) by progressive shallow water solitary waves are studied. Strong nonlinear features in the form of pronounced double peaks are visible in the time history of pressure and force signals with increasing heights of the oncoming solitary waves. The effect of nonlinearity is less pronounced as the inclination of the wall decreases with respect to the bottom surface.     

Pore characterization of organic-rich Lower Cambrian shales in Qiannan Depression of Guizhou Province,Southwestern China

Nine organic-rich shale samples of Lower Cambrian black shales were collected from a recently drilled well in the Qiannan Depression, Guizhou Province where they are widely distributed with shallower burial depth than in Sichuan Basin, and their geochemistry and pore characterization were investigated. The results show that the Lower Cambrian shales in Qiannan Depression are organic rich with TOC content ranging from 2.81% to 12.9%, thermally overmature with equivalent vitrinite reflectance values in the range of 2.92–3.25%, and clay contents are high and range from 32.4% to 53.2%. The samples have a total helium porosity ranging from 2.46% to 4.13% and total surface area in the range of 9.08–37.19 m²/g. The estimated porosity in organic matters (defined as the ratio of organic pores to the volume of total organic matters) based on the plot of TOC vs helium porosity is about 10% for the Lower Cambrian shales in Qiannan Depression and is far lower than that of the Lower Silurian shales (36%) in and around Sichan Basin. This indicates that either the organic pores in the Lower Cambrian shale samples have been more severely compacted than or they did not develop organic pores as abundantly as the Lower Silurian shales. Our studies also reveal that the micropore volumes determined by Dubinin–Radushkevich (DR) equation is usually overestimated and this overestimation is closely related to the non-micropore surface area of shales (i.e. the surface area of meso- and macro-pores). However, the modified BET equation can remove this overestimation and be conveniently used to evaluate the micropore volumes/surface area and the non-micropore surface areas of micropore-rich shales.     

A modified method to estimate chain inverse catenary profile in clay based on chain equation and chain yield envelope

Offshore floating facilities are fixed by anchoring systems embedded in seabed soils through chains or ropes. The chain inverse catenary profile embedded in soils influences both the anchor failure mechanism and the anchor holding capacity. The chain mobilizes varying soil normal and tangential resistances during motion, hence it is with difficulty to depict the chain profile. The present work proposed a modified method to estimate the chain inverse catenary profile with high accuracy based on the chain equations and the chain yield envelope. A testing arrangement with three load cells and two MEMS (Micro-electromechanical systems) accelerometers included was designed in model tests. By model tests, the loading combinations of the soil tangential and normal resistances on the chain were obtained and the yield envelopes for both chain and rope were determined. In addition, supplemental model tests were performed to validate the modified method proposed in this study, and the testing results indicated that the estimated chain inverse catenary profile was in good agreement with the actual one. Moreover, the testing arrangement is beneficial in investigating the chain-soil-anchor interaction.     

An energy-controlling boundary condition for partial wave reflections in the mild slope equation

An energy-controlling technique to actively manage the reflective property of waves from solid boundary is presented. As linear waves propagate through an energy-controlling area, a reduction in wave heights occurs due to energy dissipation, which can be placed under direct control through the imaginary part of the wavenumber and phase velocity. Based on this relationship, the present study investigates a new method to control reflected waves with desired heights in the mild slope equation model. The method is validated through numerical tests for various reflection coefficients and the results confirm the promising use of energy-controlling boundary condition for partial wave reflections.     

Explicit stress–strain equations for modeling frictional materials

In this paper, a comprehensive study on simulating the shearing behavior of frictional materials is performed. A set of two explicit equations, describing the relationship among the shear stress ratio and the distortional strain and the volumetric strain, are formulated independently. The equations contain three stress parameters and three strain parameters and another parameter representing the nonuniformity of stress and strain during softening. All the parameters have clear physical significance and can be determined experimentally. It is demonstrated that the proposed equations have the capacity of simulating the complicated shearing behavior of many types of frictional materials including geomaterials. The proposed equations are used to simulate the stress–strain behavior for 27 frictional materials with 98 tests. These materials include soft and stiff clays in both reconstituted and structured states, silicon sands and calcareous sands, silts, compacted fill materials, volcanic soils, decomposed granite soils, cemented soils (both artificially and naturally cemented), partially saturated soils, ballast, rocks, reinforced soils, tire chips, sugar, wheat, and rapeseed. It has been demonstrated that the proposed explicit constitutive equations have the capacity to capture accurately the shearing behavior of frictional materials both qualitatively and quantitatively. A study on model parameters has been performed.     

垂荡式波浪能装置运动模型数值分析

垂荡式波浪能装置在海洋可再生能源开发中被广泛应用,通过浮子与摇杆在垂荡方向的相对运动吸收波浪能。在以往相关的运动预报数值分析中,通常基于微幅波假设,仅考虑浮子与摇杆在垂荡方向的运动,忽略摇杆其他自由度运动。建立并求解了垂荡式波浪能装置的非线性联合运动方程组,分析垂荡式波浪能装置的波浪载荷、浮子与平台连接处的受力情况。数值计算系统的运动响应,并将计算结果与已有的试验数据进行比较验证,结果表明数值模拟的垂荡式波能装置的运动响应与试验结果相符合。最后,应用本计算方法分析PTO(power take-off)参数对波能装置发电性能的影响。     

MPI环境下基于边界校正方法的海洋模型的并行化

基于区域分解技术,设计了一种针对海洋模型的边界校正方法,该方法不需要对原有海洋模型进行较大改动,只需要利用消息传递接口(message passing interface,MPI)编写子区域间的数据传输模块,即可实现原有海洋模型的并行化。相对于重新开发并行系统,该方法能够较快把原有模型的串行解决方案转化为并行解决方案,并且并行化后的海洋模型能够获得较高的加速比。     

浅海内波垂直模式方程的一种解析解

给出了一些海区实测浮力频率剖面的数学模型.利用超几何函数对内波垂直模式方程进行了解析求解,得到内波模式方程的解析解和频散方程.利用得到的解析解给出了一种理想条件下内波前三阶简正模式的结构.     

Numerical modelling of deformation of multi-directional random waves over a varying topography

Two numerical formulations of the breaking phenomenon were implemented in a numerical model for random wave propagation based on the elliptic formulation of the mild-slope equation. The randomness of the wave field was simulated based on a spectral component method, in which the 3-D spectrum is discretised in components of equal energy. One of the breaking process formulations is based on the concept of breaking each independent spectral component. The other is based on the distribution of the local amount of energy dissipated through the independent spectral components. The model based on the concept of breaking each independent spectral component produces the best estimates of the wave field, when the numerical results are compared with laboratory data.     

基于扩散滤波的多尺度三维变分研究

将扩散方程引入三维变分分析,揭示了传统3D-VAR不能有效提取多尺度观测信息的根本原因,即观测导致的目标函数梯度在空间的不连续分布.将扩散滤波融入基于梯度的最优化算法,发展了基于扩散滤波的多尺度3D-VAR.海表面温度数据同化试验结果表明,新方法可从长波至短波有效地提取多尺度的观测信息.