Errors in dynamical fields inferred from oceanographic cruise data: Part II. The impact of the lack of synopticity

Diagnostic studies of ocean dynamics based on the analysis of oceanographic cruise data are usually quite sensitive to observation errors, to the station distribution and to the synopticity of the sampling. The first two sources have been evaluated in the Part I of this work. Here we evaluate synopticity errors for different sampling strategies applied to simulated unstable baroclinic waves. As suggested in previous studies, downstream and upstream cross-front samplings produce larger errors than along-front samplings. In our particular case study, the along-front sampling results in fractional errors (rms error divided by the standard deviation of the field) of about 15% for dynamic height and more than 50% for relative vorticity and vertical velocity. These values are significantly higher than those obtained in Part I for typical observation errors and sampling limitations (about 6% for dynamic height and between 15 and 30% for geostrophic vorticity and vertical velocity).We also propose and test two methods aimed at reducing the impact of the lack of synopticity. The first one corrects the observations using the quasi-geostrophic tendency equation. The second method combines the relocation of stations (based on a system velocity) and the correction of observations (through the estimation of a growth rate). For the fields simulated in this work, the second method gives better results than the first, being able to eliminate practically all synopticity errors in the case of the along-front sampling. In practice, the error reduction is likely to be less effective, since actual fields cannot be expected to have a system velocity as homogeneous as for the single-mode waves simulated in this work.     

河道三维水流数值模型及应用

天然河道水流、泥沙运动在空间上属三维问题,故提出一种三维明渠水流数学模型,将该模型应用到黄河贵德水文站下游某河段,模拟的流速、水位与实测资料吻合较好;流量沿程分布结果表明模型具有良好的物理量守恒性;给出的流速垂向分布、紊动能量等参数亦是合理的,为河道三维水流结构分析提供了依据。     

Hydrodynamic damping of a circular cylinder at low KC: Experiments and an associated model

This paper investigates the hydrodynamic damping of a smooth circular cylinder undergoing forced oscillations at Keulegan-Carpenter (KC) numbers smaller than 5 and Reynolds (Re) numbers from 10³–10⁵ with and without background steady currents. A series of experiments are conducted with a circular cylinder oscillating in still water, in-line currents and cross currents. The measured drag coefficients of the smooth cylinder in the still water condition match with the well-published results and the theoretical solution of Stokes and Wang at very small KC numbers. The hydrodynamic damping increases with the in-line steady current whereas it remains almost constant at small transverse velocities and increases notably when the latter becomes large. To predict the hydrodynamic damping in in-line steady currents, the performance of the Morison equation based on relative velocity and independent velocity is explored, respectively. The latter model, by separating the drag into two independent parts, leads to a better fit of the drag force than the former, which is not surprising. However, the former is still a preferable option for engineering design due to its simplicity. The experimental data suggest that the existing design guidelines such as ISO-19902 or DNVGL-RP-C205 should be used with caution for KC < 5.     

Numerical analysis of surface piercing propeller in unsteady conditions and cupped effect on ventilation pattern of blade cross-section

The aim of this study is to calculate hydrodynamic performance and ventilation flow around wedge, 2D blade and 3D surface piercing propeller (SPP), using computational fluid dynamic based on Reynolds-averaged Navier–Stokes method. First, numerical analyses for two-phase fluid flow around the wedge and 2D blade section (cupped and non-cupped) are presented. Flow ventilation, pressure distribution and forces are determined and compared with experimental data. Then, the method is extended to predict the hydrodynamic performance of propeller SPP-841B. The propeller exhibits a cupped blade. In the simulated configuration, SPP is one-third submerged (I = h/D = 0.33) and is working at various loadings with full ventilation occurring at low advance coefficient (J). The open water performance, pressure distribution, forces/moments and ventilation pattern on the SPP-841B model are obtained and compared with experimental data. The numerical results are in good agreement with experimental measurements, especially at high advance coefficient.     

A numerical transport model for predicting the distributions of Cd, Cu, Ni, Pb and Zn in the southern North Sea: the sensitivity of model results to the uncertainties in the magnitudes of metal inputs

A new transport model for metals (named NOSTRADAMUS) has been developed to predict concentrations and distributions of Cd, Cu, Ni, Pb and Zn in the southern North Sea. NOSTRADAMUS is comprised of components for water, inorganic and organic suspended particulate matter transport; a primary production module contributes to the latter component. Metal exchange between dissolved (water) and total suspended particulate matter (inorganic + organic) phases is driven by distribution coefficients. Transport is based on an existent 2-D vertically integrated model, incorporating a 35 × 35 km grid. NOSTRADAMUS is largely driven by data obtained during the Natural Environment Research Council North Sea Project (NERC NSP). The sensitivity of model predictions to uncertainties in the magnitudes of metal inputs has been tested. Results are reported for a winter period (January 1989) when plankton production was low. Simulated ranges in concentrations in regions influenced by the largest inflows, i.e. the NE English coast and the Southern Bight, are similar to the ranges in the errors of the concentrations estimated at the northern and southern open sea boundaries of the model. Inclusion of uncertainties with respect to atmospheric (up to ± 54%) and riverine (± 30%) inputs makes little difference to the calculated concentrations of both dissolved and particulate fractions within the southern North Sea. When all the errors associated with the inputs are included there is good agreement between computed and observed concentrations, and that for dissolved and particulate Cd, Cu and Zn, and dissolved Ni and Pb, many of the observations fall within, or are close to, the range of values generated by the model. For particulate Pb, model simulations predict concentrations of the right order, but do not reproduce the large scatter in actual concentrations, with simulated concentrations showing a bias towards lower values compared to those observed. A factor which could have contributed to observed concentrations, and which is not included in the model, is considered to be a substantial benthic input of dissolved lead during this winter period, coupled to a rapid and extensive scavenging of the dissolved lead to particles. Significant reductions in riverine and aeolian inputs of total Cd and Cu of 70% and 50%, respectively, consistent with aims of North Sea Conferences, are predicted to lead to minor decreases (~ 10%) in water column concentrations of dissolved and particulate Cd and Cu, except near river sources, where maximum reductions of ~ 30–40% may occur.     

Lagrangian data in a high-resolution numerical simulation of the North Atlantic: I. Comparison with in situ drifter data

A model/data comparison was performed between simulated drifters from a high-resolution numerical simulation of the North Atlantic and a data set from in situ surface drifters. The comparison makes use of pseudo-Eulerian statistics such as mean velocity and eddy kinetic energy, and Lagrangian statistics such as integral time scales. The space and time distribution of the two data sets differ in the sense that the in situ drifters were released inhomogeneously in space and time while the simulated drifters were homogeneously seeded at the same time over a regular 1° grid. Despite this difference, the total data distributions computed over the complete data sets show some similarities that are mostly related to the large-scale pattern of Ekman divergence/convergence.Comparisons of eddy kinetic energy and root mean square velocity indicate that the numerical model underestimates the eddy kinetic energy in the Gulf Stream extension and in the ocean interior. In addition, the model Lagrangian time scales are longer in the interior than the in situ time scales by approximately a factor of 2. It is suggested that this is primarily due to the lack of high-frequency winds in the model forcing, which causes an underestimation of the directly forced eddy variability. Regarding the mean flow, the comparison has been performed both qualitatively and quantitatively using James' statistical test. The results indicate that over most of the domain, the differences between model and in situ estimates are not significant. However, some areas of significant differences exist, close to high-energy regions, notably around the Gulf Stream path, which in the model lies slightly north of the observed path, although its strength and structure are well represented overall. Mean currents close to the buffer zones, primarily the Azores Current, also exhibit significant differences between model results and in situ estimates. Possibilities for model improvement are discussed in terms of forcings, buffer zone implementations, turbulence and mixed layer parameterizations, in light of our model/data comparison.     

三峡蓄水以来青岩子河段的演变规律*

通过大量的地形、水文、流速流向资料,详细分析了三峡蓄水以来变动回水区下段青岩子河段的泥沙淤积特点,得到了几点认识:1)三峡蓄水后青岩子河段汛期水位抬高较少,流速约有减缓,以致未大幅淤积,这与模型试验论证成果略有差别;2)青岩子河段淤积时间主要发生在9-10月间坝前水位抬升时,淤积部位主要发生在金川碛尾、五羊溪、牛屎碛深槽及桌子角;3)淤积区域与所处的河型有关;4)随着三峡水库的继续运行,青岩子河段将会发生河型转化.     

URANS simulations of catamaran interference in shallow water

This paper investigates the interference effects of wave systems on a multi-hull vessel in shallow water. A numerical analysis is made using the URANS code CFDSHIP-Iowa V.4 on the DELFT Catamaran model 372. The test matrix for numerical computations includes two separation distances (s = 0.17; 0.23) and the depth values of h/T = 8.2, 2.5 and 2, at several speeds ranging within Fr _H = 0.775–1.739. Numerical results are compared with the experimental data of the Bulgarian Ship Hydrodynamic Center, and verification and validation for resistance, sinkage and trim are also performed. Results show that, at critical speed (Fr _H ≈ 1), the presence of a finite depth significantly affects the catamaran total resistance, which, in shallower water, increases considerably with respect to deep water. At low h/T, small effects of the water depth on resistance occur at subcritical and supercritical speeds. The interference effects seem to be more relevant in shallow, rather than in deep water, with maximum IF values registered at critical speeds (Fr _H ≈ 1). Similarly to deep water, the lower the separation distance the greater the interference value. Moreover, in shallow water some negative interference is observed at Fr > 0.5. Wave patterns and wave profiles are analyzed and a comparison is made between several configurations of catamaran and a mono-hull vessel, in order to analyze how water depth and separation distance determine resistance and interference. Finally, a vortex instability study is also included.     

长江口海域波浪场模拟研究

采用MIKE21 SW波浪谱模型对长江口海域在凤凰台风期间的波浪场进行了模拟计算,考虑了水位和流速的影响,模拟结果和实测吻合较好,能够反映长江口海域的波浪分布。结果显示出水位对波高的影响是比较明显的,高水位条件下深水航道两侧和坝田区波高达到1 m以上,为泥沙的再悬浮提供了动力。     

Turn and zigzag maneuvers of a surface combatant using a URANS approach with dynamic overset grids

Unsteady Reynolds averaged Navier–Stokes (URANS) computations of standard maneuvers are performed for a surface combatant at model and full scale. The computations are performed using CFDShip-Iowa v4, a free surface solver designed for 6DOF motions in free and semi-captive problems. Overset grids and a hierarchy of bodies allow the deflection of the rudders while the ship undergoes 6DOF motions. Two types of maneuvers are simulated: steady turn and zigzag. Simulations of steady turn at 35° rudder deflection and zigzag 20/20 maneuvers for Fr = 0.25 and 0.41 using constant RPM propulsion are benchmarked against experimental time histories of yaw, yaw rate and roll, and trajectories, and also compared against available integral variables. Differences between CFD and experiments are mostly within 10 % for both maneuvers, highly satisfactory given the degree of complexity of these computations. Simulations are performed also with waves, and with propulsion at either constant RPM or torque. 20/20 zigzag maneuvers are simulated at model and full scale for Fr = 0.41. The full scale case produces a thinner boundary layer profile compared to the model scale with different reaction times and handling needed for maneuvering. Results indicate that URANS computations of maneuvers are feasible, though issues regarding adequate modeling of propellers remain to be solved.     

Recovery of North-East Atlantic temperature fields from profiling floats: Determination of the optimal float number from sampling and instrumental error analysis

Argo is an international project that is deploying an array of temperature and salinity profiling floats over the global ocean. Here we use the error formulation derived from Optimal Statistical Interpolation to estimate statistical errors associated with the recovery of the temperature field in the North-East Atlantic ocean. Results indicate that with the present distribution of floats (119 in the considered domain), scales of wavelength larger than 500 km can be recovered with a relative uncertainty (rms error relative to the standard deviation of the field) of about 7% at 50 m, 8% at 200 m and 10% at 1000 m. This corresponds to mean absolute errors of 0.111 °C at 50 m, 0.104 °C at 200 m and 0.073 °C at 1000 m.The splitting of total errors into instrumental and sampling contributions reveals that, in the present scenario, errors are more due to the small number of floats than to instrumental errors, especially at upper levels. For scales larger than 500 km this will hold true until 200–250 floats are deployed (less than 200 for deep levels). In such a simulated scenario, the number of observations and the technology become approximately equally limiting factors for the accuracy of the temperature field mapping, with total relative errors of less than 2% at upper levels and about 3% at 1000 m.     

PIV study of the effect of piston position on the in-cylinder swirling flow during the scavenging process in large two-stroke marine diesel engines

A simplified model of a low speed large two-stroke marine diesel engine cylinder is developed. The effect of piston position on the in-cylinder swirling flow during the scavenging process is studied using the stereoscopic particle image velocimetry technique. The measurements are conducted at different cross-sectional planes along the cylinder length and at piston positions covering the air intake port by 0, 25, 50 and 75 %. When the intake port is fully open, the tangential velocity profile is similar to a Burgers vortex, whereas the axial velocity has a wake-like profile. Due to internal wall friction, the swirl decays downstream, and the size of the vortex core increases. For increasing port closures, the tangential velocity profile changes from a Burgers vortex to a forced vortex, and the axial velocity changes correspondingly from a wake-like profile to a jet-like profile. For piston position with 75 % intake port closure, the jet-like axial velocity profile at a cross-sectional plane close to the intake port changes back to a wake-like profile at the adjacent downstream cross-sectional plane. This is characteristic of a vortex breakdown. The non-dimensional velocity profiles show no significant variation with the variation in Reynolds number.     

浅水中低速斜航船舶水动力预报及验证与确认分析

浅水中的斜航船舶受到浅水阻塞效应和不对称流的综合影响。为预报该运动中的船舶水动力,文章采用基于定常雷诺平均纳维—斯托克斯方程的计算流体动力学方法,对浅水中做斜航运动的船舶粘性绕流场进行数值模拟。考虑低航速运动的特点,忽略航速影响下的自由面兴波,由数值计算得到水动力系数在漂角影响下的变化规律。针对计算精度问题,在数值模拟中从验证和确认角度分析和评估计算结果：通过网格收敛性分析分析数值误差与不确定度;结合试验数据考察计算模型的误差。此外,从计算区域尺度、湍流模型、边界条件、船体下沉和纵倾作用方面对模型误差的影响因素进行探讨,可为改进计算模型、提高数值模拟精度提供参考依据。     

Numerical simulation of the water entry of a structure in free fall motion

To solve the problems concerning water entry of a structure, the RANS equations and volume of fluid （VOF） method are used. Combining the user-defined function （UDF） procedure with dynamic grids, the water impact on a structure in free fall is simulated, and the velocity, displacement and the pressure distribution on the structure are investigated. The results of the numerical simulation were compared with the experimental data, and solidly consistent results have been achieved, which validates the numerical model. Therefore, this method can be used to study the water impact problems of a structure.     

航电枢纽调节山区航道设计最高通航水位的确定

航电枢纽人工调节将引发下游航道水位和流速的变化,对船舶通航条件带来较大影响。针对航电枢纽调节下典型山区航道——兰江航道,采用数值模拟手段,计算航道上游规划姚家枢纽和灵马枢纽建成后,不同重现期洪水下的水位、流速分布等水力特性,最终得到航道设计最高通航水位。计算结果表明:10 a和5 a重现期下航道特征水位超过警戒水位较多,不利于河道防洪和航行安全;3 a重现期洪水对应水位的局部分航段流速较大,影响航道运行安全和航道运输的综合效益;结合航道建设和运行的经济效益、航道以及行船安全,认为以2 m/s适航流速作为参考条件、经试算获得的90%洪水频率对应的特征水位作为兰江航道设计最高通航水位较为合理。     

Experimental uncertainty of a physical model of a tanker moored to a terminal in a port

This paper presents the uncertainty modelling of experimental results for a physical model of a tanker moored to a terminal inside a port. The physical model was built for an oil terminal at the port of Leixões in Portugal. The model incorporates the new modified port layout, as well as a future 300 m extension of the port outer north breakwater to enhance operational conditions. The physical model tests were performed on a scale of 1:80 in the Portuguese Civil Engineering Laboratory (LNEC). A generic mooring system of four mooring lines and two fenders is simulated using a nonlinear spring system. Decay tests are carried out to evaluate the natural periods of the moored model. Then, tests are carried out for the moored model in waves. The major aim of the experimental study is to obtain novel results for the wave elevation and direction at various locations inside the port, the ship motions at six degrees of freedom, and loads on mooring lines and fenders including the modified port layout. As the physical model measurements are subjected to different types of uncertainties, a systematic uncertainty analysis is carried out here, following ITTC guidelines and recommendations, to quantify all possible sources of uncertainties. The results are discussed, and several conclusions are reached. Based on the experimental results, the presented physical model study may replicate the results for waves and motions with uncertainties less than 9% of the significant amplitudes. The assessment of the applied nonlinear spring model reveals load predictions on the moorings, with uncertainties less than 4% of the maximum mooring loads.     

Error quantification of a high-resolution coupled hydrodynamic-ecosystem coastal-ocean model: Part 2. Chlorophyll-a, nutrients and SPM

Marine systems models are becoming increasingly complex and sophisticated, but far too little attention has been paid to model errors and the extent to which model outputs actually relate to ecosystem processes. Here we describe the application of summary error statistics to a complex 3D model (POLCOMS-ERSEM) run for the period 1988–1989 in the southern North Sea utilising information from the North Sea Project, which collected a wealth of observational data. We demonstrate that to understand model data misfit and the mechanisms creating errors, we need to use a hierarchy of techniques, including simple correlations, model bias, model efficiency, binary discriminator analysis and the distribution of model errors to assess model errors spatially and temporally. We also demonstrate that a linear cost function is an inappropriate measure of misfit. This analysis indicates that the model has some skill for all variables analysed. A summary plot of model performance indicates that model performance deteriorates as we move through the ecosystem from the physics, to the nutrients and plankton.     

基于万吨溢油回收船的平面被动式减摇水舱优化和模型试验

为进一步掌握平面被动式减摇水舱内水位的变化对减摇效率的影响规律,以某万吨级溢油回收船为研究对象,选择5档水舱水位和无水状态分别进行静水横摇衰减模型试验、规则波模型试验和4级海况下的不规则波模型试验,对减摇水舱内水位的变化对目标船固有周期、阻尼特性和减摇效率的影响规律进行对比分析。研究结果表明:水舱内的水位增加对目标船的横摇固有周期和阻尼系数都有影响;当水舱内的水位在水舱总高的50%以内时,减摇效率随水位的增加而提高;当水舱内的水位在水舱总高的50%~60%范围内时,减摇效率最高;当水舱内的水位在水舱总高的60%以上时,减摇效率有下降的趋势。     

基于MRF方法和滑移网格的螺旋桨水动力性能研究

本文基于计算流体力学(CFD)方法,对多重参考系模型(MRF)及滑移网格模型(SM)在计算螺旋桨水动力性能时的差异进行了探讨。将以上两种模型应用到4381螺旋桨的水动力性能计算中,首先将计算得到的推力系数及转矩系数与试验数据进行了对比,考察了两种计算模型对螺旋桨的敞水性能的预测情况,并进一步对两种模型计算得到的螺旋桨盘面的速度场、桨叶的压力分布、桨后涡量云图等进行了对比分析。计算结果表明,滑移网格模型相较于多重参考系模型,对螺旋桨的推力系数的模拟结果误差更小,扭矩系数方面,两种模型的模拟结果相差不大;对于进速系数较大时,两种模型模拟得到的压力分布及速度分布较为相似,但对于高负荷情况,滑移网格模型可以更好地捕捉桨叶的压力分布及桨盘面处的速度分布情况;进速系数较小时,多重参考系模型可以模拟出涡结构的发散现象,而滑移网格模型可以更好的在高进速系数情况下捕捉到梢涡结构。     

Optimization of integrated scheduling of handling and storage operations at automated container terminals

Increasing demand for containerization compels container terminals to improve their performance. Uncoordinated scheduling of operations is one of the main factors accounting for poor performance at automated container terminals (ACTs). To increase land utilization efficiency and lower operational times, a new storage system called the split-platform automated storage/retrieval system (SP-AS/RS) has been introduced for temporary storage of containers. This paper describes a multi-objective mixed-integer programming (MIP) model that is based on a combination of multiple interacting sub-tasks. It is aimed at optimizing the integrated scheduling of handling and storage operations in ACTs. The MIP model objective function is to minimize delays in the loading/unloading tasks of the cranes and the travel time of vehicles and platforms in the SP-AS/RS. At the same time, a simulated annealing algorithm (SAA) that provides near-optimal solutions for the problem in a reasonable computation time is appraised. The results of this study show that the objective function of the MIP model is, on average, 58 % lower than that of the non-integrated scheduling method. On the other hand, the best objective function values obtained by the SAA indicate only a 3.7 % disadvantage in comparison with optimal values determined by the MIP model, demonstrating that the SAA is able to provide near-optimal solutions for the integrated scheduling of handling and storage operations.