Three-dimensional liquid sloshing in a tank with baffles

A numerical model has been developed to study three-dimensional (3D) liquid sloshing in a tank with baffles. The numerical model solves the spatially averaged Navier-Stokes equations, which are constructed on a non-inertial reference frame having six degree-of-freedom (DOF) of motions. The large-eddy-simulation (LES) approach is employed to model turbulence by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. The baffles in the tank are modeled by the concept of virtual boundary force (VBF) method. The numerical model is first validated against the available analytical solution and experimental data for two-dimensional (2D) liquid sloshing in a tank without baffles. The 2D liquid sloshing in tanks with baffles is then investigated. The numerical results are compared with other results from available literatures. Good agreement is obtained. Finally, the model is used to study 3D liquid sloshing in a tank with vertical baffles. The effect of the baffle is investigated and discussed.     

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.     

Implementation of the European Water Framework Directive: Natural variability associated with the CARLIT method on the rocky shores of the Ligurian Sea (Italy)

The Water Framework Directive (WFD) 2000/60/EC, adopted by the European Community in 2000 with the goal of maintaining and improving the aquatic environments, requires that member states achieve and maintain a good ecological status of all water bodies by 2015. In the marine context, the ecological status has to be quantified applying indexes based on appropriate key biological elements, which allow the categorization of water bodies into five Ecological Status (ES) classes. The CARLIT index is a cartographic monitoring tool enabling the Ecological Quality Ratio (EQR) to be calculated using macroalgae in coastal hard bottoms as a key biological element; at present it is being applied in Spain, France and Italy. To detect actual changes of water quality and, consequently, rely on the application of indexes for the assessment of the ecological status in the marine environment, it is necessary to evaluate their sensitivity to natural variability at different temporal and spatial scales. In this study we present a first quantification of natural (spatial and temporal) variability of EQR‐CARLIT quality assessment in 2006 and 2007 along the Ligurian coast (North‐Western Mediterranean, Italy). The EQR‐CARLIT values recorded along the Ligurian coastline show that natural variability of EQR‐CARLIT is low and that it does not affect the attribution of a given stretch of coast to a particular ES class, in agreement with the major human pressures acting in the area (urbanization, river load, sea‐farming). A small‐scale variability was detected, strengthening the use of cartography of the whole rocky shore, whenever possible, or, alternatively, the assessment of the ecological status for long stretches of coast, to encompass the small‐scale variability due to local factors.     

Nonlinear numerical simulation on extreme-wave kinematics

A fully nonlinear numerical model based on a time-domain higher-order boundary element method (HOBEM) is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here.     

Detection of nighttime sea fog/stratus over the Huanghai Sea using MTSAT-1R IR data

A dual channel difference (DCD) method is applied to detect nighttime sea fog/stratus over the Huanghai Sea using the infrared (IR) data of shortwave (3.5–4.0 μm) and longwave (10.3–11.3 μm) channels from the Multi-functional Transport Satellite (MTSAT)-1R, i.e., shortwave minus longwave brightness temperature difference (SLTD). Twenty-four sea fog events over the Huanghai Sea during March to July of 2006 and 2007 are chosen to determine a suitable value of SLTD for nighttime sea fog/stratus detection, and ...     

Impacts of maintenance channel dredging in a northern Adriatic coastal lagoon. I: Effects on sediment properties,contamination and toxicity

Conservation and management of coastal lagoons envisage direct human intervention. To prevent siltation and to preserve the hydrodynamics features of the lagoon system, the inner channels undergo regular maintenance dredging.     

Winkler's method overestimates dissolved oxygen in seawater: Iodate interference and its oceanographic implications

Dissolved oxygen in seawater has been determined by using the Winkler's reaction scheme for decades. An interference in this reaction scheme that has been heretofore overlooked is the presence of naturally occurring iodate in seawater. Each mole of iodate can result in an apparent presence of 1.5 mol of dissolved oxygen. At the concentrations of iodate in the surface and deep open ocean, it can lead to an overestimation of 0.52 ± 0.15 and 0.63 ± 0.05 μmol kg^− 1 of oxygen in these waters respectively. In coastal and inshore waters, the effect is less predictable as the concentration of iodate is more variable. The solubility of oxygen in seawater was likely overestimated in data sources that were based on the Winkler's reaction scheme for the determination of oxygen. The solubility equation of García and Gordon [Garcia H.E., Gordon, L.I., 1992. Oxygen solubility in seawater: Better fitting equations. Limnol. Oceanogr. 37, 1307–1312] derived from the results of Benson and Krause [Benson, F.B., Krause, D. Jr., 1984. The concentration and isotopic fractionation of oxygen dissolved in freshwater and seawater in equilibrium with the atmosphere. Limnol. Oceanogr. 29, 620–632] is free from this source of error and is recommended for general use. By neglecting the presence of iodate, the average global super-saturation of oxygen in the surface oceans and the corresponding efflux of oxygen to the atmosphere both have been overestimated by about 8%. Regionally, in areas where the degree of super-saturation or under-saturation of oxygen in the surface water is small, such as in the tropical oceans, the net air–sea exchange flux can be grossly under- or overestimated. Even the estimated direction of the exchange can be reversed. Furthermore, the presence of iodate can lead to an overestimation of the saturation anomaly of oxygen in the upper ocean attributed to biological production by 0.23 ± 0.07%. AOU may have been underestimated by 0.52 ± 0.15 and 0.63 ± 0.05 μmol kg^− 1 in the surface mixed layer and deep water, while preformed phosphate and preformed nitrate may have been overestimated by 0.004 ± 0.001 and 0.06 ± 0.02 μmol kg^− 1 in the surface mixed layer, and 0.005 ± 0.0004 and 0.073 ± 0.006 μmol kg^− 1 in the deep water. These are small but not negligible corrections, especially in areas where the values of these parameters are small. At the increasing level of sophistication in the interpretation of oxygen data, this source of error should now be taken into account. Nevertheless, in order to avoid confusion, an internationally accepted standard needs to be adopted before these corrections can be applied.     

Modified Moving Particle Semi-implicit methods for the prediction of 2D wave impact pressure

As a gridless particle method, the MPS (Moving Particle Semi-implicit) method has proven useful in a wide variety of engineering applications including free-surface hydrodynamic flows. Despite its wide range of applicability, the MPS method suffers from some shortcomings such as non-conservation of momentum and spurious pressure fluctuation. By introducing new formulations for the pressure gradient and a new formulation of the source term of the Poisson Pressure Equation (PPE), and by allowing a slight compressibility, we have proposed modified MPS methods for the prediction of wave impact pressure on a coastal structure. The improved performance of the modified methods is shown through the simulation of numerous wave impact problems (including the impacts by a dam break flow, a flip-through and two cases of slightly-breaking waves) in comparison with the experimental data.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

Nonlinear random wave-induced drag force on a vegetation field

This paper provides a practical method by which the drag force on a vegetation field beneath nonlinear random waves can be estimated. This is achieved by using a simple drag formula together with an empirical drag coefficient given by Mendez et al. (Mendez, F.J., Losada, I.J., Losada, M.A., 1999. Hydrodynamics induced by wind waves in a vegetation field. J. Geophys. Res. 104 (C8), 18383–18396). Effects of nonlinear waves are included by using Stokes second order wave theory where the basic harmonic motion is assumed to be a stationary Gaussian narrow–band random process. An example of calculation is also presented.