σ-Coordinates hydrodynamic numerical model for coastal and ocean three-dimensional circulation

A three-dimensional finite-difference hydrodynamic model has been developed using σ-coordinate for the vertical dimension. An explicit scheme for temporal integration and a staggered grid for spatial discretization have been adopted. The model has been tested against analytical or literature cases for wind and tide induced circulation. Results are in good agreement both with analytical solutions under idealised conditions and with results from the model of Shankar et al. (1996).     

A three-dimensional hydrodynamic model of estuarine circulation with an application to Southampton Water, UK

A three-dimensional hydrodynamic model has been developed to simulate water mass circulation in estuarine systems. This model is based on the primitive equation in Cartesian coordinates with a terrain-following structure, coupled with a Mellor–Yamada 2.5 turbulence scheme. A fractional-step method is applied and the subset of equations is solved with finite volume and finite element methods. A dry–wet process simulates the presence of the tidal flat at low water. River inputs are introduced using a point-source method. The model was applied to a partially mixed, macrotidal, temperate estuary: Southampton Water, UK. The model is validated by comparisons with sea surface elevation, ADCP measurements and salinity data collected in 2001. The mean spring range 2(M₂ + S₂) and the mean neap range 2(M₂ − S₂) are modelled with an error relative to observation of 12 and 16%, respectively. The unique tidal regime of the system with the presence of the ‘young flood stand’ corresponding to the slackening conditions occurring at mid flood and ‘double high water’ corresponding to an extension of the slackening conditions at high tide is accurately reproduced in the model. The dynamics of the modelled mean surface and bottom velocity closely match the ADCP measurements during neap tides (rms of the difference is 0.09 and 0.01 m s⁻¹ at the bottom and at the surface, respectively), whereas at spring the difference is greater (rms of the difference is 0.25 and 0.20 m s⁻¹ at bottom and surface, respectively). The spatial and temporal variation of the degree of stratification as indicated by salinity distributions compares well with observations.     

A three-dimensional hydrodynamic and salinity transport model of estuarine circulation with an application to a macrotidal estuary

A three-dimensional semi-implicit finite volume numerical model has been developed and applied to study tidal circulation and salinity stratification in the region of Oujiang River Estuary, China. The model employs horizontally unstructured grids and boundary-fitted coordinate system in the vertical direction. Governing equations consisting of continuity, momentum, and transport equations are all solved in the integral form of the equations, which provides a better representation of the conservative laws for mass, momentum, and transport in the coastal region with complex geometry and bottom bathymetry. The model performance was firstly quantified with skill assessment statistics on the choice of different parameters and validated with observed tidal elevation, current velocity, direction and salinity data over a spring–neap tidal cycle collected in 2006. Numerical results show that the model with wetting–drying capability successfully simulated the tidal currents and salinity fields with a reasonable accuracy and indicate that the Oujiang River Estuary is a macrotidal estuary with strong tidal mixing. In addition, the model results also show that the Oujiang River Estuary is a well-mixed estuary during spring tide. Then, the numerical simulations were performed to compare the hydrodynamic process and salinity distribution before and after a river training, which was conducted by blocking the south branch of the Oujiang River mouth. The results reveal that with the only north access to the sea, the influence of the blocking project on the flood discharge capacity is limited and the incremental velocity is beneficial to the navigation channel maintenance, although it will cause some scour to the embankment. Furthermore, the redistribution of tidal prism passing in or out the north branch makes a little severe salinity intrusion during high tide or low tide. However, the salinity intrusion is still within acceptable range, although it can cause some adverse effect on water intaking of production and life. The variations of salinity levels in Yueqing Bay situated at the north of the river mouth are not obvious, so the blocking project will not bring damage to local aquiculture. However, significant changes of salinity happen inside or outside of the south branch, so enough attention need to be paid to the changes of environment caused by the salinity variation after the blocking project. Overall, by weighing advantages and disadvantages of the blocking project, it is feasible and the model can be considered as a tool for managing and studying estuarine circulation.     

A three-dimensional ocean general circulation model for mesoscale eddies——Ⅱ Diagnostic analysis

A three-dimensional density field associated with mesoscaie unstable waves generated by the 3-D, primitive-equation model (Wang and Ikeda, 1996) is provided to the quasi-geostrophic pressure tendency and ω-equations, and to the (ageostrophic) Q-vector equation. Diagnostic analyses, analogous to the approaches in meteorology: ω-equation and Q-vector method, are for the first time developed to examine the mesoscaie dynamical processes and mechanisms of the unstable waves propagating in the mid-latitude ocean. The weaknesses and strengths of these two diagnostic approaches are evaluated and compared to the model results. The Q-vector method is then recommended to diagnose the vertical motion associated with the mesoscaie dynamics from a hydrographic CTD (conductivity-temperature-depth) array, while the quasi-geostrophic equations produce some small-scale features (errors) in the diagnosed fields.     

An ADI technique for solving three-dimensional coastal circulation

ＡｎＡＤＩｔｅｃｈｎｉｑｕｅｆｏｒｓｏｌｖｉｎｇｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｃｏａｓｔａｌｃｉｒｃｕｌａｔｉｏｎＹａｎｇＬｉａｎｗｕ１，＊ａｎｄＪｏｓｅＯｚｅｒ２（ＲｅｃｅｉｖｅｄＮｏｖｅｍｂｅｒ２７，１９９６；ａｃｃｅｐｔｅｄＤｅｃｅｍｂｅｒ...     

Comparison of two three-dimensional hydrodynamic modeling systems for coastal tidal motion

A comparison of two three-dimensional numerical modeling systems for tidal elevations and velocities in the coastal waters is presented. The two modeling systems are: (1) the Princeton Ocean Model (POM) and (2) the MIKE 3 flow model. The model performance results for Singapore's coastal waters show that the predicted tidal elevations from the two hydrodynamic modeling systems are almost identical and are in very good agreement with field measurement data. The simulated tidal current velocities match well with field measurement data at the selected stations, but it seems that the POM provides the slightly better simulation, compared to the MIKE 3 flow model. The depth profiles of the velocities obtained from the two modeling systems may be greatly different at some time, due to the vertical diffusion coefficient calculated from different turbulent sub-models in the two modeling systems. The POM generally predicts larger peak tidal velocities. The maximum speed differences for the model results from the two modeling systems occur in the top and differ from time to time and from location to location, reaching up to 20%.     

A three-dimensional model of the breeze-induced circulation of waters in the Kerch strait

We generalize the problem of evaluation of currents caused by the action of breeze winds to the three-dimensional case. In the “solid-lid” approximation, the problem is reduced to the numerical solution of a two-dimensional equation for the integral current function (with complex-valued coefficients) with subsequent evaluation of the components of current velocity according to analytic formulas. The breeze is specified as acting in a narrow coastal strip in the form of a zonal wind. The three-dimensional structure and the variations of currents near the west boundary of the Kerch Strait as functions of time are studied in detail.     

Implementation of baroclinic terms in a three-dimensional hydrodynamic model and its application to Anzali Port, Iran

Baroclinic terms have been implemented in a three-dimensional fully hydrodynamic model developed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] modifying its momentum equations to account for density gradients and utilizing the scalar (salinity, temperature, etc.) conservation equation (SCE) and a state equation for the calculation of density. In the solution of advection-diffusion terms of the governing Navier-Stokes equations (NSE) and SCE, a symmetric splitting method was applied to ensure the long-term stability of simulations. Correction terms proposed by Ruddic et al. (1995) were applied to SCE to ensure the conservation of the scalar quantity. In the presence of baroclinic terms, the zero gradient pressure in the vertical direction in the vicinity of surface and bottom boundaries assumed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] created spurious currents. This problem was solved by assuming a hydrostatic pressure variation at those boundaries. The ability of extended model was validated by comparing its results with an experimental test case. The simulation of hydrodynamic and salt intrusion at Anzali Port located at the southern coasts of Caspian Sea in Iran was carried out by the model with both barotropic and baroclinic modes. The simulated results with baroclinic mode show a better agreement with measured data as compared to the results of barotropic mode that clearly demonstrate the significance of baroclinic terms in the simulation of cyclic intrusion of salt wedge into the Port Basin.     

A three-dimensional ocean general circulation model for mesoscale eddies——Ⅰ. Meander simulation and linear growth rate

Ａｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｏｃｅａｎｇｅｎｅｒａｌｃｉｒｃｕｌａｔｉｏｎｍｏｄｅｌｆｏｒｍｅｓｏｓｃａｌｅｅｄｄｉｅｓ──ⅠＭｅａｎｄｅｒｓｉｍｕｌａｔｉｏｎａｎｄｌｉｎｅａｒｇｒｏｗｔｈｒａｔｅ￥ＷａｎｇＪｉａａｎｄＭｏｔｏＩｋｅｄａ（Ｒ...     

Quasi-isopycnic layer model for large-scale ocean circulation

A hydrothermodynamic model of a multilayer ocean, incorporating the upper mixed layer (UML) is described. The model is based on a system of primitive equations integrated within each layer. All layers are assumed to be horizontally-inhomogeneous, however, the density in each thermocline layer changes within the limits determineda priori by the prescribed basic stratification. It is assumed that the layers may outcrop. Results of the numerical experiment on subduction simulation (downwelling of UML waters toward the main thermocline layers) are given.Translated by V. Puchkin.     

OpenMP在水动力数学模型并行计算中的应用

采用PGI Fortran7.1-2(Portland Group)的OpenMP技术对二维水动力数学模型进行了并行优化试验。并行后,数学模型运行时间明显减少,在对同一算例的水流计算模拟中,串行和并行所需运行时间分别为5 336.781 s和3 454.296 s,平均加速比为1.56,平均并行计算节省时间36%,明显提高了水动力学数学模型的运算速度。     

Three-dimensional hydrodynamic model of Xiamen waters

Abstract-A semi-implicit and Eulerian-Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-momentturbulence closure model and transport equations of salinity and temperature. The simulation for flood-ing and drying of mudflats has been improved. The model is applied to Xiamen waters. Based on exten-sive survey data, water level elevation, temperature and salinity field along the eastern open boundaryand at the Jiulong River inlets and runoffs are analyzed, specified and calibrated. The computed resultsshow good agreement with the measured data, reproduce flooding, emergence of large and complexmudflat region.     

An integrated geochemical and hydrodynamic model for tidal coastal environments

In this paper, the design, calibration and application of an integrated geochemical–hydrodynamic model are described. The model comprises three parts: A hydrodynamic submodel that was adopted from a depth-averaged, semi-implicit hydrodynamic model, a geochemical submodel based on equilibrium partitioning of chemicals between aqueous and particulate phases, and a particle dynamic submodel that simulates resuspension, transport and settling of suspended particulate matter (SPM). The integrated model was implemented in San Diego Bay (SDB), a heavily urbanized, semi-closed mesotidal embayment. A series of model calibrations were carried out based on observations on salinity, polychlorinated biphenyls (PCBs) and SPM. Salinity calibrations indicated that only 15% of precipitation in the drainage area of SDB could reach the bay, presumably due to the dams on the tributary rivers. Steady-state calibrations of PCBs based on fixed concentrations at known ‘hot spots’ have reproduced observed PCB concentrations in both the dissolved and particulate phases. SPM calibrations showed that shipping-induced resuspension produce more SPM than natural processes. Based on the calibrated model, the annual transport of PCBs out of SDB was estimated to be 3.85 kg (3.5 kg and 0.35 kg in the dissolved and particulate phases, respectively), much higher than the previous estimates based on steady-state assumptions. It was also found out that only a small portion of the fine sediment exported from SDB was derived from riverine input. This model can be applied to the studies of the transport and fate of other chemical species. It can be transplanted to other coastal areas as well. The integrated model represents a novel framework in which geochemical processes in coastal environments can be investigated on a truly dynamic basis.     

一种新型三维水流数值模型

以不可压缩流体的N-S方程为基本控制方程,用快速粒子level set方法(FPLS)追踪自由表面,提出了一种新的三维水流数值模型。在自由表面处应用虚拟压力法来封闭压力泊松方程,同时用速度等值外插的方法构造自由表面外侧的虚拟速度分布。通过模拟水波振荡、水柱崩塌、水滴滴落和空箱注水过程证明了模型的有效性。     

Diagnoses of vertical transport in a three-dimensional finite element model of the tidal circulation around an island

A three-dimensional finite element model is used to investigate the formation of shallow-water eddies in the wake of Rattray Island (Great Barrier Reef, Australia). Field measurements and visual observations show that stable eddies develop in the lee of the island at rising and falling tides. The water turbidity downstream of the island suggests the existence of strong upwelling that would be responsible for carrying bed sediments up to the sea surface. We first propose to look at the upwelling velocity and then use the theory of the age to diagnose vertical transport. The water age is defined as the time elapsed since particles of water left the sea bottom, where the age is prescribed to be zero. Two versions of this diagnosis are considered. Although the model predicts upwelling within the eddies, it is not sufficiently intense to account for vertical transport throughout the water column during the life span of the eddies. As mesh resolution increases, this upwelling does not intensify. However, strong upwelling is then resolved off the island's tips, which is confirmed by the results obtained with the age. This study also shows that the finite element method, together with unstructured meshes, performs well for representing three-dimensional flow past an island.     

Multilevel hydrothermodynamic model of a deep basin

A hydrothermodynamic model based on the traditional system of differential equations is discussed. The model uses conservative finite-difference schemes based on the methods of identifying the barotropic and baroclinic velocity components and on the complete inversion of the dynamic operator. Test computations for the Black Sea basin have been conducted.Translated by Vladimir A. Puchkin.     

Wake model of hydrodynamic forces on pipelines

The hydrodynamic force model for pipelines presented includes flow history effects (wake effects) and time dependence in the force coefficients. These two features in the model were necessary to obtain satisfactory agreement between model predictions and full scale field measurements of pipeline forces. Conventional force models which represent adaptations of Morison's equation with ambient velocity and constant coefficients give predictions which for the lift component of the force in particular are in very poor agreement with the measurements. The parameters in the new model have been estimated on the basis of the full scale measurements and reflect a wide range of flow conditions. The model can be used in pipeline on-bottom stability design calculations for regular or irregular waves.     

Boussinesq相位解析的海岸水动力学数学模型研究进展

由于在平衡计算效率和精度上具有优势,Boussinesq相位解析数学模型研究不断取得突破,已成为波浪和水流精细化模拟的较优解析方式,为海岸工程、环境、地质等问题提供了实用和高效的研究手段。本文对已有Boussinesq类模型的研究进行了评述,深入探讨其重要发展、实际应用和理论瓶颈,从高阶非静压修正、GPU准三维高性能算法编译、波浪破碎和泥沙运移沉积等4个方面提出未来可能的科学突破方向。