Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

分潮波模型多资料数值计算方法研究 Ⅱ.有限差分最小二乘方法

对含多资料的分潮波模型的超定问题进行探讨，建立一种有限差分—最小二乘方法，该方法使得岸边资料直接进入了计算方程，并保证岸界法向流速为零条件。通过对渤海的四大分潮的试算，证实本方法可行。     

Simulation of wave power devices

Classical frequency and time domain models of a single degree of freedom wave power device are presented. In the time domain, a convolution integral is conventionally used to represent the fluid dynamic radiation force, characterised by added mass and damping in the frequency domain. This integral is replaced by an approximate ordinary differential equation (ODE) model which is faster and more convenient in simulations. A time domain model of the fluid dynamics of an oscillating water column (OWC) device is derived to illustrate the technique. Digital simulations of the OWC are used to compare the accuracy of the classical and ODE models. The simulation of the ODE model runs about six times as fast as the classical model based on convolution, yet characterises the fluid dynamics accurately.     

Numerical Study on Instantaneous Discharge of Unsorted Particle Cloud in Cross Flow

The mixing characteristics of particles such as dredged sediment of variable size discharged into cross flow are studied by a 3D numerical model, which is developed to model the particle-fluid two-phase flow. The Eulerian method with the modified k-ε parameterization of turbulence for the fluid phase is used to solve fluid phase, while a Lagrangian method for the solid phase (particles), both the processes are coupled through the momentum sources. In the model the wake turbulence induced by particles has been included as additional source term in the k-ε model; and the variable drift velocities of the particles are treated efficiently by the Lagrangian method in which the particles are tracked explicitly and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumping a cloud of particles is governed by the total buoyancy of the cloud, the drag force on each particle and the velocity of cross-flow. The computed results show a roughly linear relationship between the displacement of the frontal position and the longitudinal width of the particle cloud. The particle size in the cloud and the velocity of cross flow dominate the flow behavior. The computed results are compared with the results of laboratory experiments and satisfactory agreement is obtained.     

Ensuring confidence in radionuclide-based sediment chronologies and bioturbation rates

Sedimentary records of naturally occurring and fallout-derived radionuclides are widely used as tools for estimating both the ages of recent sediments and rates of sedimentation and bioturbation. Developing these records to the point of data interpretation requires careful sample collection, processing, analysis and data modeling. In this work, we document a number of potential pitfalls that can impact sediment core records and their interpretation. This paper is not intended as an exhaustive treatment of these potential problems. Rather, the emphasis is on potential problems that are not well documented in the literature, as follows: (1) the mere sampling of sediment cores at a resolution that is too coarse can result in an apparent diffusive mixing of the sedimentary record at rates comparable to diffusive bioturbation rates observed in many locations; (2) ²¹⁰Pb profiles in slowly accumulating sediments can easily be misinterpreted to be driven by sedimentation, when in fact bioturbation is the dominant control. Multiple isotopes of different half lives and/or origin may help to distinguish between these two possible interpretations; (3) apparent mixing can occur due simply to numerical artifacts inherent in the finite difference approximations of the advection diffusion equation used to model sedimentation and bioturbation. Model users need to be aware of this potential problem. Solutions to each of these potential pitfalls are offered to ensure the best possible sediment age estimates and/or sedimentation and bioturbation rates can be obtained.     

黄、东海环流的数值研究Ⅲ——正压环流的数值模拟

依据黄、东海环流的的动力学模型 ,运用“流速分解法”对黄、东海正压环流进行了数值模拟。计算结果表明冬季黄海正压环流主要受风应力影响 ,基本形态为黄海暖流由济州岛西南进入南黄海中部 ,其东西两侧分别为两支向南流动的沿岸流 ;夏季主要受到潮致体力的影响 ,为一逆时针涡旋。东海环流主要是边界力作用驱动的结果 ,东海黑潮、台湾暖流和对马暖流较稳定。冬季风应力对东海环流表层流场有消弱作用 ,在夏季则有一定增强作用。     

High-resolution sediment echosounding off Peru: Late Quaternary depositional sequences and sedimentary structures of a current-dominated shelf

About 6000 km of both bathymetric and high-resolution acoustic profiles were acquired on the shelf and upper slope offshore Peru between 9° S and 14° S. Two new sediment echosounder systems – SEL-96 and SES-2000DS – provided details of the sedimentary structures of the Quaternary sequences within the Sechura-Salaverry, Huacho and Pisco Basins. To a great extent, the poleward undercurrent determines the distribution of sediments. The undercurrent has generated numerous erosional unconformities, it has winnowed hardgrounds and has created mudwaves common between 250 m and 400 m water depth. Distinct subbottom reflectors within sedimentary units represent hiatuses due to periods of intensified winnowing or non-deposition. Erosional unconformities usually marked by pronounced reflectors suggest shifts of the undercurrent system related to climatic changes and eustatic variations of sea level. On a larger scale, the stacked prograding depositional sequences reflect the sea-level cycles of the Middle Pleistocene to the Holocene. Based on the stratigraphy of our piston cores and that of Ocean Drilling Program (ODP) Site 680, the depositional sequences limited by extended unconformities were assigned to oxygen isotope stages 1 to 7. Other sedimentary structures are small straight channels that were conduits for downslope sediment transport. Deformed sediments associated with synsedimentary normal faults result from creep movements indicating beginning slope failure.