现场海域水平扩散系数的探讨 Ⅱ．黄河三角洲埕北海域的计算

根据Ｔａｙｌｏｒ的湍流扩散理论和Ｈａｙ、Ｐａｓｑｕｉｌｌ提出的拉格朗日流和欧拉流的相关函数的相似性，用黄河三角洲埕北海域的实测资料，按四种不同方案计算出该海域的水平扩散系数和谱密度分布。指出：该海区潮流扩散系数主要受潮流的影响，其分布特征与潮流的分布基本一致，所求的水平扩散系数主要是由周期性的潮运动给予潮流扰动所致。方案４是通过富氏变换滤波消去了周期性潮流扰动，本方案计算的水平扩散系数是由周期为６ｈ以下的小尺度涡提供能量。     

Influence of large-diameter pipe pile driving on surrounding marine soils and adjacent submarine pipelines

Abstract

With the large-scale development and utilization of ocean resources and space, it is inevitable to encounter existing submarine facilities in pile driving areas, which necessitates a safety assessment. In this article, by referring to a wharf renovation project as a reference, the surrounding soil response and buried pipe deformation during pile driving in a near-shore submarine environment are investigated by three-dimensional (3D) numerical models that consider the pore water effect. Numerical studies are carried out in two different series: one is a case of a single pile focusing on the effect of the minimum plane distance of the pile–pipe, and the other is a case of double piles focusing on the effect of the pile spacing.     

酸性污染物欧拉型输送转化模式的研究

数值模式是研究酸性污染及其长距离输送特性的一种重要方法。７０年代以来，国内外发展了许多欧拉型的输送转化模式，我国在欧拉型酸性污染物长距离输送模式的研究中也取得了重要进展。主要介绍了国内外的欧拉酸性污染物输送转化模式，并就模式的发展方向进行了讨论。     

Theoretical and experimental study of breaking wave on sloping bottoms

This paper studies the continuous evolution of breaking wave for the surface water waves propagating on a sloping beach. A Lagrangian asymptotic solution is derived. According to the solution coupled with the wave breaking criteria and the equations of water particles motion, the wave deformation and the continuous wave breaking processes for the progressive water waves propagating on a sloping bottom can be derived. A series of experiments are also conducted to compare with the theoretical solution. The results show that the present solution can reasonably describe the plunging or spilling wave breaking phenomenon.     

Numerical implementation of the installation/mooring line and application to analyzing comprehensive anchor behaviors

With the application of innovative anchor concepts and advanced technologies in deepwater moorings, anchor behaviors in the seabed are becoming more complicated and pose a great challenge to the analytical methods. In the present work, a large deformation finite element (FE) analysis employing the coupled Eulerian–Lagrangian technique is performed to simulate the installation/mooring line, and then is applied to analyzing comprehensive anchor behaviors in the seabed. By connecting cylindrical units with each other using connector elements, the installation/mooring line is constructed. With the constructed installation/mooring line, FE simulations are carried out to investigate comprehensive anchor behaviors in the seabed, including long-distance penetration of drag anchors, keying of suction embedded plate anchors and non-catastrophic behavior of gravity-installed anchors. Through comparative studies, the accuracy of the proposed method is well examined. A parametric study is also undertaken to quantify the effects of the frictional coefficient, initial embedment depth, and soil weight on the profile of the embedded anchor line and the shackle load. The present work demonstrates that the proposed FE model, which incorporates the installation/mooring line and the anchor, is effective in analyzing the comprehensive anchor behaviors in the seabed.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

Numerical study of sand scour with a modified Eulerian model based on incipient motion theory

Application of the standard Eulerian model to simulations of sand scour results in unrealistic phenomena. Therefore, the present work develops a modified Eulerian model based on sand incipient motion theory. The modified model is applied for simulating a two-dimensional single vertical jet and a moving planar jet. The simulation results generally demonstrate fairly good agreement with published results of scour profiles and the velocity contours of the water and sand phases. In addition, equations to describe self-similar scour profiles for the moving planar jet cases are given. The results demonstrate that the modified model efficiently and accurately simulates the two-dimensional sand scour produced by jets, particularly for the moving jet cases.     

On the meridional circulation and balance of momentum in the Southern Ocean of POP

 The circulation of the Southern Ocean is studied in the eddy-resolving model POP (Parallel Ocean Program) by an analysis of zonally integrated balances. The TEM formalism (Transformed Eulerian Mean) is extended to include topography and continental boundaries, thus deviations from a zonally integrated state involve transient and standing eddies. The meridional circulation is presented in terms of the Eulerian, eddy-induced, and residual streamfunctions. It is shown that the splitting of the meridional circulation into Ekman and geostrophic transports and the component induced by subgrid and Reynolds stresses is identical to a particular form of the zonally integrated balance of zonal momentum. In this balance, the eddy-induced streamfunctions represent the interfacial form stresses by transient and standing eddies and the residual streamfunction represents the acceleration of the zonal current by density fluxes in a zonally integrated frame. The latter acceleration term is directly related to the surface flux of density and interior fluxes due to the resolved and unresolved eddies. The eddy-induced circulation is extremely vigorous in POP. In the upper ocean a shallow circulation, reversed in comparison to the Deacon cell and mainly due to standing eddies, appears to the north of Drake Passage latitudes, and in the Drake Passage belt of latitudes a deep-reaching cell is induced by transient eddies. In the resulting residual circulation the Deacon cell is largely cancelled and the residual advection of the zonal mean potential density is balanced by diapycnal eddy and subgrid fluxes which are strong in the upper few hundred meters but small in the ocean interior. The balance of zonal momentum is consistent with other eddy-resolving models; a new aspect is the clear identification of density effects in the zonally integrated balance. We show that the wind stress and the stress induced by the residual circulation drive the eastward current, whereas both eddy species result in a braking. Finally, we extend the Johnson–Bryden model of zonal transport to incorporate all relevant terms from the zonal momentum balance. It is shown that wind stress and induction by the residual circulation carry an eastward transport while bottom form stress and the stress induced by standing eddies yield westward components of transport. Received: 26 June 2001 / Accepted: 2 November 2001     

Turbulence Scale Dependence of the Richardson Constant in Lagrangian Stochastic Models

We investigate the relative dispersion properties of the well-mixed class of Lagrangian stochastic models. Dimensional analysis shows that, given a model in the class, its properties depend solely on a non-dimensional parameter, which measures the relative weight of Lagrangian-to-Eulerian scales. This parameter is formulated in terms of Kolmogorov constants, and model properties are then studied by modifying its value in a range that contains the experimental variability. Large variations are found for the quantity, g* = 2gC₀^{− 1}, where g is the Richardson constant.     

Energetics of lateral eddy diffusion/advection：Part Ⅱ. Numerical diffusion/diffusivity and gravitational potential energy change due to isopycnal diffusion

Study of oceanic circulation and climate requires models which can simulate tracer eddy diffusion and ad vection accurately. It is shown that the traditional Eulerian coordinates can introduce large artificial hori zontal diffusivity/viscosity due to the incorrect alignment of the axis. Therefore, such models can smear sharp fronts and introduce other numerical artifacts. For simulation with relatively low resolution, large lateral diffusion was explicitly used in models; therefore, such numerical diffusion may not be a problem. However, with the increase of horizontal resolution, the artificial diffusivity/viscosity associated with hori zontal advection in the commonly used Eulerian coordinates may become one of the most challenging ob stacles for modeling the ocean circulation accurately. Isopycnal eddy diffusion （mixing） has been widely used in numerical models. The common wisdom is that mixing along isopycnal is energy free. However, a careful examination reveals that this is not the case. In fact, eddy diffusion can be conceptually separated into two steps： stirring and subscale diffusion. Due to the thermobaric effect, stirring, or exchanging water masses, along isopycnal surface is associated with the change of GPE in the mean state. This is a new type of instability, called the thermobaric instability. In addition, due to cabbeling subscale diffusion of water parcels always leads to the release of GPE. The release of GPE due to isopycnal stirring and subscale diffusion may lead to the thermobaric instability.