A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

Influence of large-diameter pipe pile driving on surrounding marine soils and adjacent submarine pipelines

Abstract

With the large-scale development and utilization of ocean resources and space, it is inevitable to encounter existing submarine facilities in pile driving areas, which necessitates a safety assessment. In this article, by referring to a wharf renovation project as a reference, the surrounding soil response and buried pipe deformation during pile driving in a near-shore submarine environment are investigated by three-dimensional (3D) numerical models that consider the pore water effect. Numerical studies are carried out in two different series: one is a case of a single pile focusing on the effect of the minimum plane distance of the pile–pipe, and the other is a case of double piles focusing on the effect of the pile spacing.     

Estimating hydraulic conductivity using grain-size analyses, aquifer tests, and numerical modeling in a riverside alluvial system in South Korea

Hydraulic conductivity (K) for an alluvial system in a riverbank filtration area in Changwon City, South Korea, has been studied using grain-size distribution, pumping and slug tests, and numerical modeling. The alluvial system is composed of layers: upper fine sand, medium sand, lower fine sand, and a highly conductive sand/gravel layer at the base. The geometric mean of K for the sand/gravel layer (9.89?×?10^?4 m s^?1), as determined by grain-size analyses, was 3.33 times greater than the geometric mean obtained from pumping tests (2.97?×?10^?4 m s^?1). The geometric mean of K estimates obtained from slug tests (3.08?×?10^?6 m s^?1) was one to two orders of magnitude lower than that from pumping tests and grain-size analyses. K estimates derived from a numerical model were compared to those derived from the grain-size methods, slug tests and pumping tests in order to determine the degree of deviation from the numerical model. It is considered that the K estimates determined by the slug tests resemble the uppermost part of the alluvial deposit, whereas the K estimates obtained by grain-size analyses and pumping tests are similar to those from the numerical model for the sand/gravel layer of the riverside alluvial system.     

含煤层地质环境下地震波场的数值模拟

本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。     

Testing earthquake prediction methods: «The West Pacific short-term forecast of earthquakes with magnitude MwHRV ≥ 5.8»

Analyzing the tables and probability maps posted by Yan Y. Kagan and David D. Jackson in April 2002–September 2004 at http://scec.ess.ucla.edu/~ykagan/predictions_index.html and the catalog of earthquakes for the same period, the conclusion is drawn that the underlying method could be used for prediction of aftershocks, while it does not outscore random guessing when main shocks are considered.     

A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping

In this paper, a multi-method approach for the assessment of the stability of natural slopes and landslide hazard mapping applied to the Dakar coastal region is presented. This approach is based on the effective combination of geotechnical field and laboratory works, of GIS, and of mechanical (deterministic and numerical) stability analysis. By using this approach, valuable results were gained regarding instability factors, landslide kinematics, simulation of slope failure and coastal erosion. This led to a thorough assessment and strong reduction in the subjectivity of the slope stability and hazard assessment and to the development of an objective landslide danger map of the SW coast of Dakar. Analysis of the results shows that the slides were influenced by the geotechnical properties of the soil, the weathering, the hydrogeological situation, and the erosion by waves. The landslide susceptibility assessment based on this methodological approach has allowed for an appropriate and adequate consideration of the multiple factors affecting the stability and the optimization of planning and investment for land development in the city.     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.