Interaction mechanisms among waves,currents and a submerged plate

Wave-induced loads on a submerged plate, representative of submerged breakwater, coastal-bridge deck and a certain type of wave energy converter, in a uniform current are investigated in this study using fully nonlinear numerical wave tanks (NWTs) based on potential flow theory. The coupling effect of wave and current is explored, and the underlying interaction mechanisms of the hydrodynamic forces are described. The presence of a background current modifies the frequency dispersion. It produces changes of the water-surface elevation, and also has an effect on wave-induced loads. Depending on the nonlinearity, higher harmonic wave components are generated above the submerged plate. These contribute to the wave forces. It is found that the horizontal and the vertical force, hence the moment, are affected in the opposite way by the currents. The Doppler shifted effect dominates the vertical force and the moment on the plate. Whereas, the Doppler shifted effect and the generation of higher wave harmonics play opposite roles on the horizontal forces. The contribution of 2^nd order harmonics is found to be up to 30% of the linear component. The current-induced drag force, represented by the advection term ρU∂φ/∂x in the pressure equation, is found to lead to a decrease in the moment for the most range of wavelengths considered, and an increase in the moment for a small range of longer waves.     

Changes in fishing behaviour of two fleets under fully documented catch quota management: Same rules,different outcomes

A Dutch pilot study of fully documented fisheries provided the opportunity to observe actual changes in fishing behaviour under catch quota management (CQM). Interviews with fishers in the pilot study aided in interpreting the results and giving insight in the decision making process and reasoning of fishers. The CQM pilot study entailed a fleet of small and large demersal vessels. For these vessels, all cod catches were counted against quota, including catches of individuals below minimum landings size. To obtain reliable catch data all vessels were equipped with electronic monitoring (EM) systems. These systems recorded videos of all fishing and processing activities on board. In return, fishers received a 30% quota bonus for cod and were compensated with more flexibility on effort regulations. It was hypothesized that vessels in the CQM will (i) increase their landings by 30% according to their quota bonus, (ii) increase the use of gear with large mesh size, and (iii) change effort towards fishing locations with high catch rates of large cod and avoid areas with high catch rates of undersized cod. The results showed that CQM had no effect on fishing behaviour of the small vessels. In contrast, large vessels significantly increased their cod landings (216%) and avoided undersized cod. This difference in response of different fleets suggested that implementation of CQM, for instance in the context of the European Common Fisheries Policy, should consider fleet characteristics. It seemed that larger vessels in this study more easily adapted their behaviour to new management regimes and that the quota bonus opened up new fishing strategies, that were not envisaged during the implementation.     

基于离散元模拟的巨厚煤层分层开采覆岩裂隙演化分形特征

为了探究左权五里堠煤业有限公司上组煤5号煤层老空水是否对下组煤15号煤层首采面综放开采造成突水水害威胁,以下组煤首采面为背景,系统分析了研究区地质及水文地质条件,建立了两煤层开采后顶底板变形破坏空间分布特征模型,提出了老空水是否威胁下组煤综放开采的3个可能条件。在此基础上,采用具有针对性的经验公式和数值模拟相结合的方法,综合获得了下组煤15号煤层首采面最大导水裂隙带高度和上组煤5号煤层底板最大破坏深度量化参数,并与下组煤首采面附近施工的地面钻孔揭露这两层煤间距进行对比。认为下组煤15号煤层开采后,通过导水裂隙带将会波及到上组煤采动底板破坏带内部,在正常条件下可能与上组煤5号煤层老空水导通,建议开采前要进行提前探查和疏放。研究结果对下组煤首采工作面顶板老空水害防治具有很好的参考价值。     

A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalised associated Legendre functions

 Spherical harmonic expansions form partial sums of fully normalised associated Legendre functions (ALFs). However, when evaluated increasingly close to the poles, the ultra-high degree and order (e.g. 2700) ALFs range over thousands of orders of magnitude. This causes existing recursion techniques for computing values of individual ALFs and their derivatives to fail. A common solution in geodesy is to evaluate these expansions using Clenshaw's method, which does not compute individual ALFs or their derivatives. Straightforward numerical principles govern the stability of this technique. Elementary algebra is employed to illustrate how these principles are implemented in Clenshaw's method. It is also demonstrated how existing recursion algorithms for computing ALFs and their first derivatives are easily modified to incorporate these same numerical principles. These modified recursions yield scaled ALFs and first derivatives, which can then be combined using Horner's scheme to compute partial sums, complete to degree and order 2700, for all latitudes (except at the poles for first derivatives). This exceeds any previously published result. Numerical tests suggest that this new approach is at least as precise and efficient as Clenshaw's method. However, the principal strength of the new techniques lies in their simplicity of formulation and implementation, since this quality should simplify the task of extending the approach to other uses, such as spherical harmonic analysis. Received: 30 June 2000 / Accepted: 12 June 2001     

Finite volume coupling strategies for the solution of a Biot consolidation model

In this paper a finite volume (FV) numerical method is implemented to solve a Biot consolidation model with discontinuous coefficients. Our studies show that the FV scheme leads to a locally mass conservative approach which removes pressure oscillations especially along the interface between materials with different properties and yields higher accuracy for the flow and mechanics parameters. Then this numerical discretization is utilized to investigate different sequential strategies with various degrees of coupling including: iteratively, explicitly and loosely coupled methods. A comprehensive study is performed on the stability, accuracy and rate of convergence of all of these sequential methods. In the iterative and explicit solutions four splits of drained, undrained, fixed-stress and fixed-strain are studied. In loosely coupled methods three techniques of the local error method, the pore pressure method, and constant step size are considered and results are compared with other types of coupling methods. It is shown that the fixed-stress method is the best operator split in comparison with other sequential methods because of its unconditional stability, accuracy and the rate of convergence. Among loosely coupled schemes, the pore pressure and local error methods which are, respectively, based on variation of pressure and displacement, show consistency with the physics of the problem. In these methods with low number of total mechanical iterations, errors within acceptance range can be achieved. As in the pore pressure method mechanics time step increases more uniformly, this method would be less costly in comparison with the local error method. These results are likely to be useful in decision making regarding choice of solution schemes. Moreover, the stability of the FV method in multilayered media is verified using a numerical example.     

上海市浦西地区地下水三维数值模拟

以上海市浦西地区第二承压含水层和第三承压含水层为例,利用实际资料以及部分插值的方法,建立了该区地下水系统完整三维数学模型,较好地解决了相对隔水层缺失区越流系数无法调试的问题,为上海市整个地下水含水系统建立完整三维水流模型以及沉降模型奠定了基础.     

Fully coupled finite element model for dynamics of partially saturated soils

Understanding the response of partially saturated earth structures under various static and dynamic loads is important for the design and construction of economical and safe geotechnical engineering structures. In this study, the numerical approach is used to understand the dynamics of partially saturated soils. The mathematical equations governing the dynamics of partially saturated soils are derived based on the theory of mixtures and implemented within a finite element framework. The stress–strain behavior of the soil is represented by an elasto-plastic constitutive model for unsaturated soil based on bounding surface concept and the moisture-suction behavior is modeled using van Genuchten model. Fully coupled finite element simulations are performed to study the response of partially saturated soil embankment under earthquake loading and validated with centrifuge test results available in the literature. The predicted displacement responses are in good agreement with the measured responses. The pore water pressure, pore air pressure, matric suction, the degree of saturation in various elements and the response of the embankment under different initial moisture content are also discussed.     

Fully nonlinear wave-body interactions with surface-piercing bodies

Fully nonlinear wave-body interactions for stationary surface-piercing single and double bodies are studied by a potential-theory-based fully nonlinear 2D numerical wave tank (NWT). The NWT was developed in time domain by using boundary element method (BEM) with constant panels. MEL free surface treatment and Runge–Kutta fourth-order time integration with smoothing scheme was used for free-surface time simulation. The acceleration-potential scheme is employed to obtain accurate time derivative of velocity potential. Using the steady part of nonlinear force time histories, mean and a series of higher-harmonic force components are calculated and compared with the experimental and numerical results of other researchers. The slow-decaying second-harmonic vertical forces are investigated with particle velocities and corresponding body pressure. Typical patterns of two-body interactions, shielding effect, and the pumping/sloshing modes of water column in various gap distances are investigated. The pumping mode in low frequencies is demonstrated by the comparison of velocity magnitudes.     

Fully nonlinear 3D interaction of bubble dynamics and a submerged or floating structure

This paper is concerned with the interaction of bubbles, a submerged or floating structure, and free surface waves. A three-dimensional fully nonlinear model has been developed based on the coupling of the boundary integral method (BIM) for bubble dynamics and free surface waves and the finite element method for structure deformation. The present method is well validated by comparing the numerical results with the experimental data. Three structure characteristics, including fixed, rigidly moving and flexible, are investigated separately to determine their influence on bubble dynamics. For a free-floating structure, the free surface causes not only a larger reduction in peak pressure for a rigid structure compared with a fixed body but also the modification of the bubble period and structural response. The interaction between a bubble and a flexible structure, in the absence of a free surface, is simulated. Both the rigid motion and the deformation at the local structure appear in the simulation. The effect of the structural thickness on the reduction in peak pressure is also considered.     

Time domain three-dimensional fully nonlinear computations of steady body–wave interaction problem

Three-dimensional fully nonlinear waves generated by moving disturbances with steady forward speed without motions are solved using a mixed Eulerian–Lagrangian method in terms of an indirect boundary integral method and a Runge–Kutta time marching approach which integrates the fully nonlinear free surface boundary conditions with respect to time.A moving computational window is used in the computations by truncating the fluid domain (the free surface) into a computational domain. The computational window maintains the computational domain and tracks the free surface profile by a node-shifting scheme applied within it. An implicit implement of far field condition is enforced automatically at the truncation boundary of the computational window.Numerical computations are applied to free surface waves generated by Wigley and Series 60 hulls for the steady problem. The present numerical results are presented and compared with existing linear theory, experimental measurements, and other numerical nonlinear computations. The comparisons show satisfactory agreements for these hydrodynamic problems.