波浪与抛石潜堤相互作用的MPS法数值模型研究

为了研究波浪与抛石潜堤相互作用过程中大自由表面变形和堤内渗流等强非线性紊流运动问题,利用改进的MPS法,建立了模拟波浪与抛石潜堤相互作用的MPS法数值计算模型。模型将抛石潜堤假定为均质多孔介质,采用Drew的二相流运动方程描述多孔介质内外的流体运动;通过在动量方程中增加非线性阻力项,并引入亚粒子尺度紊流模型,模拟波浪与可渗结构物相互作用过程中的紊流运动。选取“U”型管中多孔介质内渗流过程和孤立波与可渗潜堤相互作用两个典型的渗流问题,通过将数值计算结果与理论解和实测值的对比分析,对所提出的MPS法紊流渗流模型的模拟精度进行验证。结果表明:基于改进的MPS法构建的垂向二维紊流渗流模型可以很好地再现“U”型管中多孔介质内渗流以及波浪作用下可渗潜堤内外的复杂流场,显著缓解流-固界面处的压力震荡与粒子分布不均匀问题,实现了较高的模拟精度。     

完全非线性孤立波的直墙反射

报道了应用边界积分方法模拟完全非线性孤立波的传播与直墙反射,给出了波形演变过程。结果表明,本模型对计算孤立波的传播与直墙反射是有效的。三阶Boussinesq方程的孤立波解比低阶方程的孤立波解更接近完全非线性的数值解.当来波波高增大时,孤立波直墙反射的相位滞后变小。若考虑大波高孤立波的直墙反射或波——波相互作用,一阶理论预报的相位滞后往往低估实际情况。     

波浪与大孔隙多孔介质相互作用的耦合数学模型

建立了波浪与大孔隙多孔介质相互作用的耦合数学模型,波浪域的控制方程为雷诺时均方程和k-ε紊流模型。对于计算域的入射波采用推板式造波,它可以是线性波、椭圆余弦波和孤立波。采用PLIC-VOF法追踪波浪自由表面。对于多孔介质内的孔隙流场采用非线性Forchheimer方程,两区域共享连续方程,最后导出的波浪域与孔隙流域的压力修正方程具有完全相同的形式,利用这个方程能够同时而不是分别求解波浪场和孔隙流场,避免了在内部边界上给定匹配条件,实现了波浪场与孔隙流场的同步耦合。波浪与粗颗粒海床、平底床面上抛石潜堤及斜坡上抛石潜堤相互作用的验证计算结果表明该模型可用于研究波浪与大孔隙多孔介质相互作用的问题。     

孤立波对礁坪上直墙冲击试验和RANS数值模拟

在有、无礁冠两种典型岛礁剖面地形上开展孤立波冲击礁坪上直墙的二维水槽试验和数值模拟研究。通过OpenFOAM结合k-ε湍流模型进行RANS数值模拟,基于试验数据和验证后的RANS模拟结果,研究了岛礁剖面上孤立波传播的水动力特性及其对直墙的作用,并探讨了礁冠的存在对直墙所受冲击的影响。结果表明:RANS(k-ε)模型能够准确还原孤立波在岛礁地形上的主要传播特征,但在波浪破碎的情况下,RANS数值模拟不够精细,尤其对直墙所受瞬击荷载的捕捉能力欠缺;礁冠的存在,增强了波浪反射,使得波浪提前破碎,减小了直墙上的波浪爬升和最大动压。     

吕宋海峡内波传播演化的数值模拟

吕宋海峡生成的内波在往南海西部传播的过程中,具有明显的二维辐射特征,这是传统的一维KdV模式无法刻画的.基于弱二维的KP方程,并结合南中国海吕宋海峡附近内孤立波的观测资料,模拟了潮成内波的演化特征.数值结果较好地反应了内孤立波的二维特征.相比于一维的KdV方程,KP在内孤立波的仿真反演方面具有更大的优势.     

内孤立波波-波相互作用的数值模拟

基于弱二维的KP方程，并结合南中国海东沙群岛附近内孤立波的观测资料，模拟了内孤立波的波-波相互倌用0数值结果较好的反应了内孤立波的二维特征，同时体现两个内孤立波波-波相互作用的非线性特征，即两波相交处相速随振幅的增大而变大。相比于一维的KdV方程，KP在内孤立波的仿真反演方面具有更大的优势。     

孤立波与带窄缝双箱相互作用模拟研究

针对孤立波与带窄缝双箱的作用问题,应用时域高阶边界元方法建立了二维数值水槽。其中,自由水面满足完全非线性运动学和动力学边界条件,对瞬时自由表面流体质点采用混合欧拉-拉格朗日法追踪,采用四阶龙格库塔法对下一时刻的自由水面的速度势和波面升高进行更新。采用加速度势法求解物体湿表面的瞬时波浪力。采用推板方法生成孤立波。通过模拟孤立波在直墙上的爬高以及施加在直墙上的波浪力,并与已发表的实验和数值结果对比,验证本数值模型的准确性。通过数值模拟计算研究了窄缝宽度、方箱尺寸对波浪在箱体迎浪侧爬高,窄缝内波面升高,箱体背浪侧透射波高及箱体受波浪荷载的影响。同时研究了有一定时间间隔的双孤立波与带窄缝双箱系统作用问题。     

波浪与多孔介质结构相互作用SPH模拟

应用经CSPM法和黎曼解修正后的光滑粒子流体动力学(SPH)方法,建立了主动吸收无反射数值波浪水槽,研究波浪作用下多孔介质结构的水动力特性。流体运动控制方程采用N-S方程,多孔渗水结构内流体的运动控制方程考虑渗流力的影响。数值计算结果给出了水槽内不同位置测点的波面历时曲线和越浪量随时间变化曲线,并同试验结果和Philip Liu的数值计算结果进行了比较。并对一个波浪周期内斜坡堤多孔介质结构内外的速度场和压力分布进行了讨论分析。计算分析表明,数值计算波面较Philip Liu的计算结果与试验结果吻合更好。说明应用SPH方法建立的二维数学模型能够较好地模拟破碎波在多孔渗水斜坡上的爬坡和越浪。     

内孤立波数值模拟方法研究

本文选用RNGk-ε湍流模型封闭二维不可压缩N-S方程组,选用KdV、mKdV理论作为输入条件进行内孤立波造波模拟。运用流体模拟软件Fluent建立二维内波数值水槽,用动网格手段实现网格变形模拟物理边界造波法,用VOF方法进行密度分层流体的内界面追踪,对不同非线性强度的二相流内孤立波进行了数值模拟。并将数值模拟结果与KdV、mKdV理论波形以及内波试验水槽所造内孤立波进行对比,验证了此套造波方法在Fluent中进行内孤立波模拟的可行性。     

孤立波爬坡的二维数值模拟

波浪在滩地上以及遇海岸工程后传播发生变形、爬坡等现象,对其进行数值模拟具有广阔的工程应用背景.应用基于Boltzmann方程的KFVS(kinetic flux vector splitting)格式求解二维浅水方程,同时采用干底Riemann解模拟动边界问题.模型模拟了孤立波在滩地上传播变形、爬坡的过程,以及孤立波在滩地上遇圆柱后绕射、变形和爬坡的过程.计算结果与实验结果非常吻合,表明模型具有较大的推广应用价值.     

Microbial Biomass and Organic Nutrients in the Deep-Sea Sediments of the Central Indian Ocean Basin

In order to assess the impact of deep-sea mining on the in situ benthic life, we measured the microbial standing stock and concentration of organic nutrients in the deep-sea sediments of the Central Indian Ocean Basin in the Indian pioneer area. Sediments were collected using box core and grab samples during September 1996. The total bacterial numbers ranged from 10 10 -10 11 cells per g -1 dry weight sediment. There was a marginal decrease in the number of bacteria from surface to 30 cm depth, though the subsurface section registered a higher number than did the surface. The highest numbers were encountered at depths of 4-8 cm. The retrievable number of bacteria were two orders less in comparison with the direct total counts of bacteria. An almost homogeneous distribution of bacteria, total organic carbon, living biomass, and lipids throughout the depth of cores indicates active microbial and benthic processes in the deep sea sediments. On the other hand, a uniform distribution of total counts of bacteria, carbohydrates, and total organic carbon in all the cores indicates their stable nature and suggests that they can serve as useful parameters for long-term monitoring of the area after the benthic disturbance. Further studies on temporal variability in this region would not only verify the observed norms of distribution of these variables but would also help to understand restabilization processes after the simulated benthic disturbance.     

A Review of Some Environmental Issues Affecting Marine Mining

This article reviews information recently available from existing marine and coastal mining for responses to environmental issues affecting marine mining at different depths. It is particularly but not exclusively concerned with those issues affecting seabed biodiversity impact and recovery. Much information has been gathered in the past 10 years from shallow mining operations for construction aggregate, diamonds, and gold, from coastal mines discharging tailings to shallow and deep water, and from experimental deep mining tests. The responses to issues identified are summarized in a series of eight tables intended to facilitate site-specific consideration. Since impacts can spread widely in the surface mixing layer SML, and can affect the biologically productive euphotic zone, the main issues considered arise from the depth of mining relative to the SML of the sea. Where mining is below the SML, the issue is whether it is environmentally better to bring the extraction products to the surface vessel for processing (and waste discharge), or to process the extraction products as much as possible on the seabed. Responses to the issues need to be sitespecific, and dependent on adequate preoperational environmental impact and recovery prediction. For deep tailings disposal from a surface vessel, there are four important environmental unknowns: (1) the possible growth of "marine snow" (bacterial flocs) utilizing the enormous quantities of fine tailings particles (hundreds or thousands of metric tons per day) as nuclei for growth, (2) the possibility that local keystone plankton and nekton species may migrate diurnally down to and beyond the depth of deep discharge and hence be subjected to tailings impact at depth, (3) the burrow-up capability of deep benthos and their ability to survive high rates of tailings deposition, and (4) the pattern and rate of dispersion of a tailings density current through the deep water column from discharge point to seabed. Actions to obtain relevant information in general and site-specifically are suggested.     

Increased Particle Fluxes at the INDEX Site Attributable to Simulated Benthic Disturbance

Particle fluxes were measured 7 m above the sea bottom during the predisturbance, disturbance, and postdisturbance periods by using time series sediment traps attached to seven deep-sea moorings deployed in the INDEX experiment site in the Central Indian Basin. The predisturbance particle fluxes varied between 22.3 to 55.1 mg m -2 day -1 . Increased and variable particle fluxes were recorded by the sediment traps during the disturbance period. The increase observed was 0.5 to 4 times more than the background predisturbance fluxes. The increases in particle fluxes (~4 times) recorded by the sediment trap located in the southwestern direction (DMS-1) were the greatest, which could be the result of preferential movement of resuspended particles generated during the deep-sea benthic disturbance along the general current direction prevailing in this area during the experimental period. Also, the traps located closer to the disturbance area recorded greater fluxes than did the traps far away, across the Deep Sea Sediment Resuspension System path. This variability in recorded particle fluxes by the traps around the disturbance area clearly indicates that physical characteristics such as grain size and density of the resuspended particles produced during the disturbance had an important effect on particle movement. The postdisturbance measurements during ~5 days showed a reduction in particle fluxes of ~50%, indicating rapid particle settlement.     

Assessing the Validity of Seismo-Acoustic Predictor Equations for Obtaining Seabed and Sub-Surface Sediment Physical Properties

The interdependence between the seismo-acoustic properties of a marine sediment and its geotechnical/physical parameters has been known for many years, and it has been postulated that this should allow the extraction of geotechnical information from seismic data. Though in the literature many correlations have been published for the surficial layer, there is a lack of information for greater sediment depths. In this article, a desktop study on a synthetic seafloor model illustrates how the application of published near-surface prediction equations to subsurface sediments (up to several tens of meters burial depth) can lead to spurious predictions. To test this further, acoustic and geotechnical properties were measured on a number of sediment core samples, some of which were subjected to loading in acoustically-equipped consolidation cells (oedometers) to simulate greater burial depth conditions. For low effective pressures (representing small burial depths extending to around 10 meters subsurface), the general applicability of established relationships was confirmed: the prediction of porosity, bulk density, and mean grain size from acoustic velocity and impedance appears generally possible for the investigated sedimentary environments. As effective pressure increases through, the observed relationships deviate more and more from the established ones for the near-surface area. For the samples tested in this study, in some instances increasing pressure even resulted in decreasing velocities. There are several possible explanations for this abnormal behavior, including the presence of gas, overconsolidation, or bimodal grain size distribution. The results indicate that an appropriate depth correction must be introduced into the published prediction equations in order to obtain reliable estimates of physical sediment properties for greater subsurface depths.     

Geoacoustic Inversion via Genetic Algorithm and Its Application to Manganese Sediment Identification

An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-waveletconvolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L 2 norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.     

海上大直径钢管桩打桩过程中桩周土强度弱化研究

海上大直径钢管桩打桩过程中,桩周土体受到强烈扰动而发生强度弱化,掌握桩周土体强度弱化规律对于准确预测打桩过程、保证工程安全具有重要意义。为研究土体强度弱化规律,开展了环剪试验模拟打桩对桩周土体的扰动,测试土体强度随剪切速率的变化规律,建立了描述土体强度弱化规律的拟合公式,引入到打桩分析软件中。研究结果表明：土体的强度折减程度不仅与土体本身的性质有关还受到土体的埋深和剪切速率的影响,埋深越深土体强度折减程度越低,剪切速率越高土体强度折减越高,在打桩分析中可采用这里推荐的线性折减方法来模拟不同深度处土体强度的折减规律。     

Determination of In Situ Sediment Shear Strength from Advanced Piston Corer Pullout Forces

The advanced piston corer (APC) has been used by the Ocean Drilling Program since 1985 for recovering soft sediments from the ocean floor. The pullout force measured on extracting the core barrel from the sediment is shown to correlate with the average shear strength of the sediment core measured in the ship's laboratory. A simple rule of thumb is derived relating the shear strength of the sediment to the pullout force. Multiple APC holes at individual sites allow the consistency of the pullout measurements to be assessed. The effects of different operational procedures during APC coring are also explored. Although generally applicable, the correlation between pullout force and laboratory measurements of shear strength breaks down for some APC holes, possibly because of the disturbance of some sediment types during the APC coring process. A better understanding of the physical process of APC coring, and its effect on the properties of the sediment both inside and immediately outside the core barrel, would indicate what confidence can be put on the measurement of pullout force as a way of evaluating the in situ shear strength of deep sea sediments.     

Interannual Sea Level Variations and Annual Tides in the Northwestern Pacific

Specific properties of the interannual sea level variations and annual tides in the Northwestern Pacific were studied. Several tide stations were monitored. The monthly mean sea level for the year of 1995 was analyzed at each tide station. A seismic event in 1995, some tectonic activity around the subject area, and the Kuroshio (the oceanic western boundary current) may possibly contaminate results which would have occurred from the astronomical annual tide alone.     

《Acta Oceanologica Sinica》Vol. 33(2014) contents

<正>~~