Numerical assessment of tsunami attack on a rubble mound breakwater using OpenFOAM®

A numerical assessment study of tsunami attack on the rubble mound breakwater of Haydarpasa Port, located at the southern entrance of the Istanbul Bosphorus Strait in the Sea of Marmara, Turkey, is carried out in this study using a Volume-Averaged Reynolds-Averaged Navier-Stokes solver, IHFOAM, developed in OpenFOAM^® environment. The numerical model is calibrated with and validated against the data from solitary wave and tsunami overflow experiments representing tsunami attack. Furthermore, attack of a potential tsunami near Haydarpasa Port is simulated to investigate effects of a more realistic tsunami that cannot be generated in a wave flume with the present state of the art technology. Discussions on practical engineering applications of this type of numerical modeling studies are given focusing on pressure distributions around the crown-wall of the rubble mound breakwater, and the forces acting on the single stone located behind the crown-wall at the rear side of the breakwater. Numerical modeling of stability/failure mechanism of the overall cross-section is studied throughout the paper.The present study shows that hydrodynamics along the wave flume and over the breakwater can be simulated properly for both solitary wave and tsunami overflow experiments. Stability of the overall cross-section can only be simulated qualitatively for solitary wave cases; on the other hand, the effect of the time elapsed during tsunami overflow cannot be reflected in the simulations using the present numerical tool. However, the stability of the overall cross-section under tsunami overflow is assessed by evaluating forces acting on the rear side armor unit supporting the crown-wall of the rubble mound breakwater as a practical engineering application in the present paper. Furthermore, two non-dimensional parameters are derived to discuss the stability of this armor unit; and thus, the stability condition of the overall cross-section. Approximate threshold values for these non-dimensional parameters are presented comparing experimental and numerical results as a starting point for engineers in practice. Finally, investigations on the solitary wave and tsunami overflow experiments/simulations are extended to the potential tsunami simulation in the scope of both representation of a realistic tsunami in a wave flume and stability of the rubble mound breakwater.     

Quadrature-based approach for the efficient evaluation of surge hazard

The Joint Probability Method (JPM) has been used for hurricane surge frequency analysis for over three decades, and remains the method of choice owing to the limitations of more direct historical methods. However, use of the JPM approach in conjunction with the modern generation of complex high-resolution numerical models (used to describe winds, waves, and surge) has become highly inefficient, owing to the large number of costly storm simulations that are typically required. This paper describes a new approach to the selection of the storm simulation set that permits reduction of the JPM computational effort by about an order of magnitude (compared to a more conventional approach) while maintaining good accuracy. The method uses an integration scheme called Bayesian or Gaussian-process quadrature (together with conventional integration methods) to evaluate the multi-dimensional joint probability integral over the space of storm parameters (pressure, radius, speed, heading, and any others found to be important) as a weighted summation over a relatively small set of optimally selected nodes (synthetic storms). Examples of an application of the method are shown, drawn from the recent post-Katrina study of coastal Mississippi.     

Balancing deformation in NW Borneo: Quantifying plate-scale vs. gravitational tectonics in a delta and deepwater fold-thrust belt system

Recent GPS measurements demonstrate that NW Borneo undergoes 4–6 mm of plate-scale shortening a year, which is not accommodated by plate-scale structures. The only geological structure in NW Borneo described to accommodate on-going shortening is the Baram Delta System located on the outer shelf to basin floor. Delta toe fold-thrust belts are commonly thought to be caused by margin-normal compressional stresses generated by margin-parallel upslope gravitational extension.     

Experimental and numerical investigation of local scour around a submerged vertical circular cylinder in steady currents

Local scour around a submerged vertical circular cylinder in steady currents was studied both experimentally and numerically. The physical experiments were conducted for two different cylinder diameters with a range of cylinder height-to-diameter ratios. Transient scour depth at the stagnation point (upstream edge) of the cylinder was measured using the so-called conductivity scour probes. Three-dimensional (3D) seabed topography around each model cylinder was measured using a laser profiler. The effect of the height-to-diameter ratio on the scour depth was investigated. The experimental results show that the scour depth at the stagnation point is independent on cylinder height-to-diameter ratio when the later is smaller than 2. The increase rate of equilibrium scour depth with cylinder height increases with an increase in Shields parameter.     

Modelling the effect of wave overtopping on nearshore hydrodynamics and morphodynamics around shore-parallel breakwaters

Wave overtopping nearshore coastal structures, such as shore-parallel breakwaters, can significantly alter the current circulation and sediment transport patterns around the structures, which in turn affects the formation of tombolos and salients in the nearshore area. This paper describes the implementation of a wave overtopping module into an existing depth-averaged coastal morphological mode: COAST2D and model applications to investigate the effect of wave overtopping on the hydrodynamics and morphodynamics around a group of shore-parallel breakwaters. The hydrodynamic aspects of the model were validated against a series of laboratory conditions. The model was then applied to a study site at Sea Palling, Norfolk, UK, where 9 shore-parallel segmented breakwaters including 4 surface-piercing and 5 low-crested breakwaters are present, for the storm conditions in Nov 2006. The model results were compared with laboratory data and field measurements, showing a good agreement on both hydrodynamics and morphological changes. Further analysis of wave overtopping effect on the nearshore hydrodynamics and morphodynamics reveals that wave overtopping has significant impacts on the nearshore circulation, sediment transport and the resulting morphological changes within such a complex breakwater scheme under the storm and macro-tide conditions. The results indicate the importance of including the wave overtopping in modelling nearshore morphodynamics with the presence of coastal structures.     

库车坳陷克拉2气田异常流体压力演化史

根据克拉2气田现今实测地层压力分析了异常高压发育的特征,并通过流体包裹体测试求取包裹体形成时期的古压力。计算结果表明:下白垩统和古近系库姆格列木群砂岩储层在距今8～6Ma(康村组沉积期末)时,古压力系数大致为1.09～1.18,发育过剩压力为4.6～9.6MPa的异常高压。在此基础上,以古、今压力数据为约束条件,利用数值模拟技术定量恢复出下白垩统与侏罗系异常流体压力的演化历史。研究发现:下白垩统储层流体压力经历了沉积型异常流体压力形成发展(白垩纪-库车组沉积期末)和沉积型异常流体压力萎缩—构造挤压型异常流体压力孕育(库车组沉积期末至今)两个阶段。康村组沉积早期沉积型超压开始出现,库车组沉积期末达到顶峰,过剩压力达40～50MPa;此后,在沉积型超压萎缩—构造挤压型超压孕育过程的耦合下,过剩压力曾降至25～30MPa,更新世起超压略有升高。侏罗系流体压力演化较为简单,吉迪克组沉积期开始出现超压,至库车组沉积期末达到高峰,此后异常高压逐渐萎缩。在克拉2气田的主要成藏时期(距今3～1Ma),烃源岩较储层具有更高的异常压力,源—储压力差为10～30MPa。     

Numerical and experimental study of a flapping foil generator

The energy extraction performance of a flapping foil generator is studied through experiment and numerical simulation. A practical flapping foil generator has been proposed. The heave motion and the pitch motion of the foil are adjusted through a crankshaft-like structure. The heave and pitch motions of the foil are transferred to the rotational motion of the main shaft. A pair of gears is adopted to increase the pitch angle. A prototype with pitch amplitude θ₀ = 60^∘ has been built and the experiment is carried out in a tunnel. The overall performance of the mechanism has been analysed. Good agreement of numerical results and experiment data has been found. Further simulations with larger pitch amplitudes are carried out. It is found that higher efficiency can be achieved with larger pitch amplitude at medium frequency.     

A survey of numerical methods for stochastic differential equations

The development of numerical methods for stochastic differential equations has intensified over the past decade. The earliest methods were usually heuristic adaptations of deterministic methods, but were found to have limited accuracy regardless of the order of the original scheme. A stochastic counterpart of the Taylor formula now provides a framework for the systematic investigation of numerical methods for stochastic differential equations. It suggests numerical schemes, which involve multiple stochastic integrals, of higher order of convergence. We shall survey the literature on these and on the earlier schemes in this paper. Our discussion will focus on diffusion processes, but we shall also indicate the extensions needed to handle processes with jump components. In particular, we shall classify the schemes according to strong or weak convergence criteria, depending on whether the approximation of the sample paths or of the probability distribution is of main interest.     

青藏高原对冷空气爆发影响的个例分析(一)

本文通过对1979年4月10—13日在我国的一次寒潮爆发过程进行了有、无地形及修改地形的数值试验和对比分析,发现高原地形对冷空气活动强度及路径等有明显的影响:(1)由于高原地形的存在使冷空气快速沿地形边缘向南推移,冷锋偏南,给我国平原地区带来较大的降温及较多的降雨。(2)这种冷锋的移动快而达到很偏南的纬度主要是低层冷空气沿高原地形边缘绕流所产生的。(3)如果高原大地形不存在,这时冷空气直侵印度,给印度北方带来强降温。     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。