Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.     

桶形基础平台三维有限元稳定性分析

以我国第一座桶形基础采油平台工程实例为原型，采用三维有限元计算方法对桶形基础稳定性分析。平台在波浪、风等荷载作用下，其稳定性不仅与平台的整体结构有关，而且与桶形基础及其周围土体相互作用相关。在计算中，利用薄层摩擦单元模拟土与结构间的相互作用，取得了很好的效果。分析结果表明桶形基础采油平台在设计条件下运行是较为稳定的。     

Variation in Water Masses on Shelf of Northern Bering Sea in September 2016-2018

Abstract

Based on hydrological data obtained during the 7th to 9th Chinese National Arctic Research Expeditions in the summers of 2016–2018, the main water structure on the shelf of the northern Bering Sea and the volume and heat fluxes of the Bering Strait throughflow were analyzed. Distinct variability was identified in the three Pacific water masses feeding the strait - Anadyr Water (AW), Bering Shelf Water (BSW) and Alaskan coastal water (ACW). Overall, the temperature and salinity of the entire section increased each year, with 2018 showing significant anomalies, i.e., a temperature anomaly of up to 1?°C and a maximum salinity anomaly of 2. From 2016 to 2018, the extent of the ACW gradually narrowed in the eastern part of this section, while the AW expanded eastward each year. The net volume transport through each of the three sections increased poleward from 1.65?Sv to 2.76?Sv, with the AW increasing from 0?Sv to 1.03?Sv, the BSW varying between 0.52–1.65?Sv, and the ACW gradually decreasing from 1.04?Sv to disappearing completely. The net heat fluxes were also poleward, varying between 25.77 TW and 61.50 TW, and showing a significant increase. Significant variations in magnitude and extent were observed in each water mass of the Bering Strait throughflow, which could produce widespread effects in the Arctic Ocean and the global ocean beyond.     

Numerical simulation of solitary wave scattering by a circular cylinder array

The interaction of a solitary wave with an array of surface-piercing vertical circular cylinders is investigated numerically. The wave motion is modeled by a set of generalized Boussinesq equations. The governing equations are discretized using a finite element method. The numerical model is validated against the experimental data of solitary wave reflection from a vertical wall and solitary wave scattering by a vertical circular cylinder respectively. The predicted wave surface elevation and the wave forces on the cylinder agree well with the experimental data. The numerical model is then employed to study solitary wave scattering by arrays of two circular cylinders and four circular cylinders respectively. The effect of wave direction on the wave forces and the wave runup on the cylinders is quantified.     

Numerical simulation of the hydrodynamic behaviour of submerged plane nets in current

The aim of this paper is to investigate the shape and tension distribution of fishing nets in current. A numerical model is developed, based on lumped mass method to simplify the net. The motion equation is set up for each lumped mass. The Runge–Kutta–Verner fifth-order and sixth-order method is used to solve these simultaneous equations, and then the displacement and tension of each lumped mass are obtained. In order to verify the validity of the numerical method, model tests have been carried out. The results by the numerical simulation agree well with the experimental data.     

黄海、渤海底层温度数值预报方法的研究

根据文献［ｌ］建立的底层温度（ＴＨ）与其水柱垂向平均温度（）的经验关系，结合流体动力学方程和（垂向平均）热传导方程，发展了以水气温差和风速为已知量的底层温度二维数值预报模式。该模式避开了海面热量和动量输入在垂直水柱中分配的复杂物理过程而直接报出底层水温场，具有较好的实用性；此外，从试报结果看，效果令人满意。     

Wave interaction with a perforated wall breakwater with a submerged horizontal porous plate

This study examines the hydrodynamic performance of a new perforated-wall breakwater. The breakwater consists of a perforated front wall, a solid back wall and a submerged horizontal porous plate installed between them. The horizontal porous plate enhances the stability and wave-absorbing capacity of the structure. An analytical solution based on linear potential theory is developed for the interaction of water waves with the new proposed breakwater. According to the division of the structure, the whole fluid domain is divided into three sub-domains, and the velocity potential in each domain is obtained using the matched eigenfunction method. Then the reflection coefficient and the wave forces and moments on the perforated front wall and the submerged horizontal porous plate are calculated. The numerical results obtained for limiting cases are exactly the same as previous predictions for a perforated-wall breakwater with a submerged horizontal solid plate [Yip, T.L., Chwang, A.T., 2000. Perforated wall breakwater with internal horiontal plate. Journal of Engineering Mechanics ASCE 126 (5), 533–538] and a vertical wall with a submerged horizontal porous plate [Wu, J.H., Wan, Z.P., Fang, Y., 1998. Wave reflection by a vertical wall with a horizontal submerged porous plate. Ocean Engineering 25 (9), 767–779]. Numerical results show that with suitable geometric porosity of the front wall and horizontal plate, the reflection coefficient will be always rather small if the relative wave absorbing chamber width (distance between the front and back walls versus incident wavelength) exceeds a certain small value. In addition, the wave force and moment on the horizontal plate decrease significantly with the increase of the plate porosity.     

The Working Mechanism and An Analytic Method of the Deep Embedded Large Cylinder Structure

The large cylinder is a new-type structure that has been applied to harbor and offshore engineering. An analytic method of the relationship between loads and the structure displacement is developed based on the failure mode of deep embedded large cylinder structures. It can be used to calculate directly the soil resistance and the ultirnate bearing capacity of the structure under usage. A new criterion of the large cylinder structure, which discriminates the deep embedded cylinder from the shallow embedded cylinder, is defined. Model tests prove that the proposed method is feasible for the analysis of deep embedded large cylinder structures.     

Parabolic Wave Propagation Modelling in Orthogonal Coordinate Systems

This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild-slope equation in the orthogonal coordinate system. Two types of coordinate systems are studied: (a) a general form of orthogonal coordinate system and (b) the conformal system, a special form of orthogonal coordinate system. Two typical examples, namely, expanded breakwaters and a circular channel, are studied to validate the model. First, the examples are studied by use of the general orthogonal coordinates. Then the same examples are computed by use of the confonnal system. The computational results show that the confonnal coordinate system generally gives better predictions than the general orthogonal system. A numerical technique for generating the conformal grid is combined with the numerical model to improve the practicability of the model. The comparison between the result from the numerical grid system and that from the analytical grid system shows that reliable computational results can be obtained by use of the numerical confonnal grid system.     

双层局部开孔板沉箱对波浪反射的理论研究

提出了一种用于研究由双层开孔板和一个不透水后板的开孔结构对斜向波反射率的理论分析方法。整个流域被分成三个子域，在每个子域内应用特征函数展开法以得到该域内包含未知展开系数的势函数的表达式，在速度势的展开中，考虑了非传播模态波浪的影响。通过匹配开孔板处的边界条件可以求解待定的展开系数，继而求解双层开孔板防波堤结构对斜向波的反射率。数值计算结果与试验结果进行了比较，符合较好。并进一步讨论了几个重要因素对反射系数的影响。