Three-dimensional liquid sloshing in a tank with baffles

A numerical model has been developed to study three-dimensional (3D) liquid sloshing in a tank with baffles. The numerical model solves the spatially averaged Navier-Stokes equations, which are constructed on a non-inertial reference frame having six degree-of-freedom (DOF) of motions. The large-eddy-simulation (LES) approach is employed to model turbulence by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. The baffles in the tank are modeled by the concept of virtual boundary force (VBF) method. The numerical model is first validated against the available analytical solution and experimental data for two-dimensional (2D) liquid sloshing in a tank without baffles. The 2D liquid sloshing in tanks with baffles is then investigated. The numerical results are compared with other results from available literatures. Good agreement is obtained. Finally, the model is used to study 3D liquid sloshing in a tank with vertical baffles. The effect of the baffle is investigated and discussed.     

Adaptation of the vertical resolution in the mixed layer for HYCOM

This paper focuses on the dynamics of the mixed layer. When the mixed layer depth increases, the vertical discretisation eventually becomes too sparse at the bottom of this layer to accurately resolve its evolution and strong numerical errors can appear. This is linked to the fact that the vertical resolution is concentrated in the upper part of the ocean and does not adapt to the deepening of the mixed layer.Knowing that the HYbrid Coordinate Ocean Model (HYCOM) is able to modify the distribution of the vertical levels, we suggest in this paper a method to adapt the resolution to the mixed layer extension. This method is tested in 1-D configurations for two academic atmospheric forcing conditions (strong convection and wind-mixing) and a realistic forcing extending over one year, with seasonal restratification following strong winter convection. The new method improves the results in all cases, and in particular when the mixed layer reaches deep layers.     

Application of Line Boundary Technique to 2D Tidal Current Simulation

To deal with the problems concerning the shore boundary,moving boundary and engineeringboundary which are encountered frequently in 2D tidal current simulation by the finite difference method,theconcept of line boundary is introduced and studied here,and then the line boundary technique in common useis proposed in this paper.Analysis of some calculation cases shows that this technique is practical,effective,and simple in 2D tidal current simulation involving different boundaries.     

基于新型高频声学多普勒流速剖面仪的河口近底部边界层观测初探

河口物质输运、能量交换与底边界层内的水动力过程密切相关,底边界层参数(如切应力、拖曳系数)的确定至关重要。挪威Nortek公司生产的新型声学多普勒流速剖面仪AD2CP相比传统ADCP具有高频、低噪的优点,可用于高频(16Hz)流速剖面观测,而被广泛应用于底边界层观测的ADV只能测量单点高频流速。本文采用AD2CP在长江口南槽最大浑浊带区域进行座底式观测,并与同步近底部三脚架上ADV的观测结果进行对比。结果表明,使用AD2CP测得的近底部平均流速与ADV的测量结果吻合良好;使用惯性耗散法计算了底切应力,基于ADV的单点高频流速数据计算结果为2.16×10~(-2)～5.69×10~(-1)N/m~2,基于AD2CP的结果为2.09×10~(-2)～4.26×10~(-1)N/m~2,二者范围大致相当。在此基础上,基于AD2CP数据计算出摩阻流速为4.55×10~(-3)～2.06×10~(-2)m/s、底拖曳系数范围为1.84×10~(-4)～2.49×10~(-3),与ADV的计算结果基本一致。此外,由于AD2CP可以获得高频的流速剖面数据,优于单点ADV,具备观测近底部边界层参数和边界层内湍流剖面的潜力。     

Scour depth under pipelines placed on weakly cohesive soils

We here study the scouring processes that evolve around a submarine pipeline placed on a weakly cohesive seabed. We first analyze some laboratory tests carried out by Vijaya Kumar et al. [21], Xu et al. [25] and Zhou et al. [28] that focused on the scouring around a horizontal cylinder lying on a cohesive bed, subject to waves and currents. The specific purpose is that of finding a new formula for the prediction of the equilibrium scour depth under submarine pipelines. After a theoretical analysis of the main parameters, the sought formula has been found to be a function of: (i) the hydrodynamic forces acting on the cylinder (through the Keulegan–Carpenter parameter KC), (ii) the clay content of the soil C_c, and (iii) the burial depth e₀/D. In the presence of small amounts of clay (C_c< 5 %), the scour depth depends directly on KC (as confirmed by many literature works for pipelines lying on sandy soils, e.g. [18]) and inversely on C_c (as already seen for bridge abutments on cohesive soils, e.g. [1]), the best-fit law being characterized by a coefficient of determination R² = 0.62. If some burial depth is accounted for, this being a novelty of the present work, a more general formulation can be used, valid in the presence of weakly-cohesive soils and with burial depths of the pipe smaller than 0.5 (R² = 0.79). For large clay-content ranges (2% < C_c < 75 %), the scour depth depends directly on both KC and C_c, this giving R² = 0.79 (no burial depth) and 0.91 (some burial depth). However, this finding is at odds with the main literature, because, for large amounts of clay, it is fundamental to consider the liquidity index LI, which accounts for some important clay properties, like the plasticity. We argue that the absence of LI is balanced by the direct dependence of the scour depth on C_c. Notwithstanding the small number of available data, a formula for the prediction of the scour depth under pipelines lying on cohesive soils is fundamental for several engineering applications. The present contribution represents the first attempt to build such a formula, when the pipeline is subject to the wave-current forcing and the seabed is characterized by a relatively small clay content.     

Floc size variability under strong turbulence: Observations and artificial neural network modeling

The flocculation of cohesive sediment in the presence of waves is investigated using high-resolution field observations and a newly-developed flocculation model based on artificial neural networks. Vertical profiles of suspended sediment concentration and turbulent intensity are estimated using measurements of current profile and acoustic backscatter. The vertical distribution of floc size is estimated using an artificial neural network (ANN) that is trained and validated using floc size measurements at one vertical level. Data analysis suggests a linear correlation between suspended sediment concentration and turbulence intensity. Observations and numerical simulations show that floc size is inversely related to sediment concentration, turbulence intensity and water temperature. The numerical results indicate that floc growth is supported by low concentration and low turbulence. In the vertical direction, mean size of flocs decreases toward the bottom, suggesting floc breakage due to increasing turbulence intensity toward the bed. A significant decrease in turbulent shear could occur within the bottom few-cm, related to increased damping of turbulence by sediment induced density stratification. The results of the numerical simulations presented here are consistent with the concept of a cohesive sediment particle undergoing aggregation-fragmentation processes, and suggest that the ANN can be a precise tool to study flocculation processes.     

南海冬季之西边界流：1998年2月的水文结构与体积输运

The unique survey in December 1998 mapped the entire western boundary area of the South China Sea(SCS),which reveals the three-dimensional structure and huge volume transport of the swift and narrow winter western boundary current of the SCS(SCSwwbc) in full scale. The current is found to flow all the way from the shelf edge off Hong Kong to the Sunda Shelf with a width around 100 km and a vertical scale of about 400 m. It appears to be the strongest off the Indo-China Peninsula, where its volume transport reached over 20×10~6 m~3/s. The current is weaker upstream in the northern SCS to the west of Hong Kong. A Kuroshio loop or detached eddy intruded through the Luzon Strait is observed farther east where the SCSwwbc no more exists. The results suggest that during the survey the SCSwwbc was fed primarily by the interior recirculation of the SCS rather than by the"branching" of the Kuroshio from the Luzon Strait as indicated by surface drifters, which is likely a near-surface phenomenon and only contributes a minor part to the total transport of the SCSwwbc. Several topics related to the SCSwwbc are also discussed.     

南黄海盆地二叠系地震地层特征与识别

南黄海盆地钻遇二叠系的井仅有4口,对大面积无井区(如崂山隆起)的地层特征认识不统一。从二叠系岩性特征和海陆钻井地层对比为出发点,结合二叠系内部及其上下地层的地震波组和速度特征,认为南黄海盆地二叠系分上下2段。上段为大隆—龙潭组含煤碎屑岩地层,煤系地层较薄且分布不稳定,存在相变特征,使得该组地层在地震反射特征上的规律性变差,地震上表现为较连续、中—高频反射和低速特征;下段为栖霞组灰岩夹泥页岩地层,地震上表现为较连续、中—低频反射和高速特征。在无井区识别二叠系时,抓住二叠系地层界面极性为"2负+1正"的强反射特征。下二叠统—中、上石炭统地震波组特征相似,由4~5个同相轴组成,厚度稳定(时间约为200 ms),地震波极性表现为"上正下负"的地层界面反射特征,可作为全区追踪的标准层。     

淤泥质、粉沙质及沙质海岸航道回淤统一计算方法

根据20世纪80年代初提出的淤泥质海岸航道回淤计算方法的基础上,开展了拓展研究,使之统一适用于淤泥质海岸、粉沙质海岸及沙质海岸。主要研究内容有,在波浪和潮流综合作用下挟沙力含沙量研究中引入了特定的泥沙因子F1/F,从而挟沙力含沙量公式不仅适用于淤泥质泥沙,也适用于非淤泥质泥沙;在动力因素方面,除了一般寻常潮和波浪动力外,进一步考虑了风暴潮和破波的巨大掀沙能力造成航道骤淤的可能性。并结合连云港30万吨级主航道扩建及徐圩港区10万吨级航道的回淤问题(包括骤淤可能性问题),京唐港外航道和黄骅港外航道的骤淤问题,进行预测计算及校验计算,结果都是良好的。此外,对杭州湾强潮和涌潮情况下的挟沙力含沙量也进行了校验计算,结果也非常满意。