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.     

小振幅海洋内波的演变、破碎和所致混合

利用基于谱方法和MPI并行运算的数值模式SpectralModel,直接数值模拟了三维小振幅海洋内波的演变、破碎和所致湍流混合,指出导致其不稳定而破碎的为PSI(parametric subharmonic instability)机制;对于内波破碎所致的湍流混合过程,分析了跨等密度面扩散系数kρ、混合效率γ、浮力通量谱、动能谱以及势能谱等统计性质:内波破碎前,kρ和γ保持低值水平,浮力通量谱值为负,且集中在低波数段;内波破碎后,kρ和γ迅速增大,最大值分别约为0.9×103m2/s和0.18,浮力通量谱值在低波数段为负值,在高波数段为正值,这是因为层化湍流中势能向小尺度运动传递和动能向小尺度运动传递相比更为有效。在内波破碎、强湍流混合阶段,势能谱存在一谱段满足kz3律,P(kz)=0.2N2kz3。此外,与二维模拟结果相比较,导致内波不稳定而破碎的均为PSI机制,kρ、浮力通量谱、势能谱变化趋势大体一致;但三维数值实验中,内波破碎时间提前,湍流衰减加快;KE谱在高波数部分下降速度相对减小,更接近于kz3。     

伶仃洋海域潜在性富营养化评价与研究

根据2006年7月至2007年12月的调查资料,一年四季分析珠江口伶仃洋海域营养盐的空间分布特征及季节变化,该海域DIP表层年平均值为0.054mg/L,底层年平均值为0.53mg/L,DIN表层平均值为1.20mg/L,底层平均值为1.03mg/L;采用近海海域水体潜在性富营养化等级的划分标准和评价方法进行评价,评价结果表明,伶仃洋海域营养化水平基本上为磷限制潜在性富营养化至III级富营养水平。     

山区年、月平均温度推算方法的研究

本文提出了山区年、月平均温度空间分布的推算方法,并提供了相应的小地形温度订正值。该方法在福建省沙溪流域山区实际应用的结果表明,能有效地揭示山区年、月平均温度空间分布的宏观与微观变化,效果满意,精度较高。     

利用位场及其高阶导数的功率谱计算岩层的深度与厚度

本文导出了利用重力与磁异常及其高阶导数的功率谱计算密度体和磁性体上顶与下底埋深的近似公式,还提供了利用垂向或水平导数谱计算深度的方法。理论模型的试验结果表明,由近似公式计算的精度可以满足实际需要。本文最后根据一条重力剖面资料,计算了密度层上、下界面。     

浅源强震的判据和发震概率

根据确定潜在震源区的多判据原理,强震的发生不仅要满足准静态强度条件,而且必须具备突发释放足够能量的条件,而所有这些条件可以通过突发应力降判据、总能量判据和强度判据加以表达。本文建立了它们相应的概率表达式和估计发震总概率的方法,并针对我国28例地震进行了验算     

应用模式识别定量划分潜在震源区

本文以陕西关中及部分邻区为例,考虑地质构造、地震活动及地球物地场等因素,采用不同的模式识别方法和多种计算方案,以定量判定潜在震源区。其结果表明,模式识别方法有利于地质、地震活动及地球物理场等多种因素综合应用,并能选择和显示判定潜有震源区的主要特征。不同模式识别方法的比较和多种试验方案的综合,则可提高潜在震源区判定的可靠性和稳定性。     

垂直地震力作用下弹性半空间上扁壳基础的解析解

用布希涅斯克定义的弹性半空间内的垂直位移包括两项积分,除了积分号前面系数的差别之外,第一项积分是单层位势而第二项积分为双层位势。若扁壳基础是正高斯曲率的几何曲面,则壳底与半空间表面间的挤压强度就是半空间表面作用的分布垂直荷载。当越过边界时,双层势位的函数值和单层势位的法向导数值发生跳跃。利用这些性质,本文得出布希涅斯克积分的反演公式,从而避开要求解偏微分—积分方程组的巨大数学困难而易于得出解析解。以椭园抛物面扁壳为例说明本文方法的应用。