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.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depth

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropriate wave theory. At the downstream boundary, an artificial damping zone is used to prevent wave reflection back into the computational domain. Using the image Green function in the whole fluid domain, the integrations on the two lateral walls and bottom are excluded. The simulation results on extreme wave elevations in finite and infinite water-depths are compared with experimental results and second-order analytical solutions respectively. The wave kinematics is also discussed in the present study.     

The connection between Labrador Sea buoyancy loss, deep western boundary current strength, and Gulf Stream path in an ocean circulation model

The sensitivity of the North Atlantic gyre circulation to high latitude buoyancy forcing is explored in a global, non-eddy resolving ocean general circulation model. Increased buoyancy forcing strengthens the deep western boundary current, the northern recirculation gyre, and the North Atlantic Current, which leads to a more realistic Gulf Stream path. High latitude density fluxes and surface water mass transformation are strongly dependent on the choice of sea ice and salinity restoring boundary conditions. Coupling the ocean model to a prognostic sea ice model results in much greater buoyancy loss in the Labrador Sea compared to simulations in which the ocean is forced by prescribed sea ice boundary conditions. A comparison of bulk flux forced hindcast simulations which differ only in their sea ice and salinity restoring forcings reveals the effects of a mixed thermohaline boundary condition transport feedback whereby small, positive temperature and salinity anomalies in subpolar regions are amplified when the gyre spins up as a result of increased buoyancy loss and convection. The primary buoyancy flux effects of the sea ice which cause the simulations to diverge are ice melt, which is less physical in the diagnostic sea ice model, and insulation of the ocean, which is less physical with the prognostic sea ice model. Increased salinity restoring ensures a more realistic net winter buoyancy loss in the Labrador Sea, but it is found that improvements in the Gulf Stream simulation can only be achieved with the excessive buoyancy loss associated with weak salinity restoring.     

Coupling winds to ocean surface currents over the global ocean

A Wind stress–Current Coupled System (WCCS) consisting of the HYbrid Coordinate Ocean Model (HYCOM) and an improved wind stress algorithm based on Donelan et al. [Donelan, W.M., Drennan, Katsaros, K.B., 1997. The air–sea momentum flux in mixed wind sea and swell conditions. J. Phys. Oceanogr. 27, 2087–2099] is developed by using the Earth System Modeling Framework (ESMF). The WCCS is applied to the global ocean to study the interactions between the wind stress and the ocean surface currents. In this study, the ocean surface current velocity is taken into consideration in the wind stress calculation and air–sea heat flux calculation. The wind stress that contains the effect of ocean surface current velocity will be used to force the HYCOM. The results indicate that the ocean surface velocity exerts an important influence on the wind stress, which, in turn, significantly affects the global ocean surface currents, air–sea heat fluxes, and the thickness of ocean surface boundary layer. Comparison with the TOGA TAO buoy data, the sea surface temperature from the wind–current coupled simulation showed noticeable improvement over the stand-alone HYCOM simulation.     

1998年季风爆发期南海大气边界层的日变化

在南海尚未有关于海洋大气边界层结构日变化方面的研究.研究海洋大气边界层,有利于我们更好地研究海洋表层结构变化影响机制,对发生在海洋大气边界层上的过程进行预测预报.因此了解南海海洋大气边界层的结构特征,对于我们进一步了解南海天气尺度海气相互作用有非常重要的意义.通过对1998年在南海南部(6°15′N,110°E)和南海北部(20°29′39″N,116°57′48″E)的南海季风试验中定点科考船释放一天四次的探空观测资料分析得出,季风爆发前海洋大气边界层存在规则的日变化,在中午达到深厚.季风爆发后南海北部大气边界层日变化消失,南海南部依然明显.分析表明对其日变化起重要影响的是短波辐射;潜热输送对大气边界层高度日变化影响不大.与大气边界层厚度日变化相对应,南海南部边界层内水汽日变化明显,而南海北部较不明显.     

极限波浪运动特性的非线性数值模拟

利用时域高阶边界元方法建立了模拟极限波浪运动的完全非线性数值模型,其中自由水面满足完全非线性自由水面条件.采用半混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶Runge-Kutta方法更新下一时间步的波面和速度势,同时应用镜像格林函数消除水槽两个侧面和底面上的积分.研究中利用波浪聚焦的方法产生极限波浪,并且在水槽中开展了物理模型实验,将测点试验数据与数值结果进行了对比,两者吻合得很好.对极限波浪运动的非线性和流域内速度分布进行了研究.     

Simulation of nonlinear wave run-up with a high-order Boussinesq model

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions. As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability to accurately simulate a moving wet–dry boundary is of considerable practical importance within coastal engineering, and the extension described in this work significantly improves the nearshore versatility of the present high-order Boussinesq approach.     

港域波浪数学模型的改进与验证

通过物理模型对改进的港内波浪传播变形数学模型进行验证。该数学模型以推广的时变缓坡方程为控制方程,采用含松弛因子的ADI法求解,并对波浪反射和透射边界模拟方法进行改进。先通过物理模型试验确定斜向浪入射条件下抛石防波堤前的波浪反射系数,作为数学模型中部分反射边界模拟的依据。然后进行了一个典型港口内波浪折射、绕射和反射的模型试验,测量港内波浪分布。对比模型试验和数学模型计算的结果表明,数学模型可较好地模拟港内复杂地形和边界条件下规则波和不规则波的传播变形。     

Nonlinear numerical simulation on extreme-wave kinematics

A fully nonlinear numerical model based on a time-domain higher-order boundary element method (HOBEM) is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here.     

天然气水合物的上界面——硫酸盐还原-甲烷厌氧氧化界面

定义了海域天然气水合物成矿带的上界面。指出在地球深部存在最原始的、从根本上不依靠光合作用来生存的生命系统。根据对ODP岩心样品中微生物数量的统计,海底以下沉积层中的生物数量可能占据全球原核生物总量的70%,其生物总碳量和地球表面所有植物的碳总量相当。地球内部如此巨大的生物总量应该在地壳中的气体分布等方面起着重要作用。甲烷在地壳层中广泛存在,并主要是微生物成因的。微生物产甲烷的途径主要有两个,一是二氧化碳还原,另一个是醋酸盐发酵。相应地,参与产甲烷的微生物菌群主要是产甲烷菌和食醋酸菌。甲烷在沉积层中的厌氧氧化是一个不争的事实。该过程发生在海底以下一个非常局限的区带,称为硫酸盐还原-甲烷厌氧氧化区带。通常,这个区带很窄,仅为一个面,因此,硫酸盐还原-甲烷厌氧氧化区带又称硫酸盐还原-甲烷厌氧氧化界面。这是一个基本的生物地球化学界面,在功能上它起到屏蔽甲烷向海底和大气逸散的作用,是一个巨大的甲烷汇。甲烷的厌氧氧化同样是一个由微生物介导的过程,参与此过程的微生物主要是食甲烷古菌和硫酸盐还原菌。硫酸盐还原-甲烷厌氧氧化界面在海洋沉积层中一般深可达海底以下上百米,浅可至海底。此界面为天然气水合物的上界面,该界面以上没有甲烷...