Calculation of transverse hydrodynamic coefficients using computational fluid dynamic approach

Computational Fluid Dynamic (CFD) based on Reynolds Averaged Navier–Stokes equation is used for determining the transverse hydrodynamic damping force and moment coefficients that are needed in the maneuverability study of marine vehicles. Computations are performed for two geometrical shapes representing typical AUVs presently in use. Results are compared with available data on similar geometries and from some of the available semi-empirical relations. It is found that the CFD predictions compares reasonable well with these results. In particular, the CFD predictions of forces and moments are found to be nonlinear with respect to the transverse velocity, and therefore both linear and nonlinear coefficients can be derived. A discussion on the sources of the component forces reveal that the total force and moment variations should in fact be nonlinear.     

Computational study of sail performance in upwind condition

Airflow around yacht sails with imposed final geometry is simulated using a CFD code, reproducing experimental tests carried in a wind tunnel. Two configurations are considered: one, mast–main sail and the other, mast–jib–main sail. Both configurations were studied in the same flow conditions of air velocity and incidence angle. The grid is structured-like next to the sails and unstructured in the rest of the domain. The turbulence model used is Shear Stress Transport. The results are compared against experimental and numerical results.     

Investigation of a longfin inshore squid's swimming characteristics and an underwater locomotion during acceleration

In the present study, locomotion of a real longfin inshore squid (Doryteuthis pealeii) was numerically investigated. Geometry of a real squid was obtained using computed tomography (CT) images. In addition to a two-dimensional axisymmetric squid model, a modified squid model with no cavities around her head and an ellipse shaped model were generated with a fineness ratio (the ratio of body length to maximum body diameter) of 7.56. These numerical models were exposed to an acceleration with two different velocity programs. Added mass and basset forces on bodies were calculated during acceleration of the squid models starting from rest. Pressure and viscous drag forces were also calculated due to pressure variation along the squid models and friction on the surface of the models. The effect of a nozzle diameter on jet velocities and propulsive efficiency at all bodies were evaluated when time dependent velocity profiles (from 0 to 10 m/s in 0.5 and 1 s time durations) were set for the inlet of computational domain. The modified squid model required least thrust force during acceleration phase of time dependent velocity profile compared to the other models while the 0.02 m nozzle diameter provided largest propulsive efficiency for all models.     

Direct scantling assessment of propeller blades

Traditionally, propeller design has been focused on all activities necessary to obtain a propeller featuring a high efficiency, avoiding erosive cavitation for given operating conditions and having adequate structural strength. In recent years, more and more challenging requirements have been imposed, such as the reduction of radiated noise and pressures pulses, requiring more precise analyses and methods in the optimization of the propeller performance. On the other hand, the evaluation of the propeller strength still relies on simplified methods, which basically consider the blade as a cantilever beam subjected to characteristic static forces. Since the loads acting on a blade are variable in the blade revolution and in different operating conditions throughout the ship life, a procedure to account for the influence of fatigue phenomena is proposed. The fatigue assessment could reduce the safety factor in the propeller scantling rules and allow improving the quality of propeller design (e.g. obtaining higher efficiency, margin on cavitation phenomena, less noise).     

CFD investigations of OXYFLUX device,an innovative wave pump technology for artificial downwelling of surface water

No other environmental variable of such ecological importance to estuarine and coastal marine ecosystems around the world has changed so drastically, in such a short period of time, as dissolved oxygen in coastal waters. The prevalent methods for counteracting anoxic sea events are indirect measures which aim to cut-down anthropic loads introduced in river and marine environments. To date, no direct approaches, like artificial devices have been investigated except the WEBAP and OXYFLUX devices. The present paper adopts a numerical approach to the analysis of the pumped surface water as well as the analysis of the dynamic response of the OXYFLUX device. By means of a CFD-RANS code and through the application of overset grid method, the 1/16 OXYFLUX model’s dynamic response and pumping performance are evaluated. The appropriate grid is selected after an extensive sensitivity analysis carried out on 9 different grids. The CFD model is validated by comparing numerical and physical results of heave decay test, heave response, and surface water discharge under the action of regular waves. The extensive comparison with experimental results shows consistently accurate predictions. The main findings of the study show that nonlinear effects remarkable reduce the dynamic behaviour of the OXYFLUX and generate an unexpected second harmonic for pitch response intensifying the overtopping discharge also for small waves caused by the summer's low intensity winds.     

Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of random variability in bed elevation

Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

CFD simulation of airflow over a regular array of cubes. Part II: analysis of spatial average properties

In the first part of this study, results of a computational fluid dynamics simulation over an array of cubes have been validated against a set of wind-tunnel measurements. In Part II, such numerical results are used to investigate spatially-averaged properties of the flow and passive tracer dispersion that are of interest for high resolution urban mesoscale modelling (e.g. non resolved obstacle approaches). The results show that vertical profiles of mean horizontal wind are linear within the canopy and logarithmic above. The drag coefficient, derived from the numerical results using the classical formula for the drag force, is height dependent (it decreases with height). However, a modification of the formula is proposed (accounting for subgrid velocity scales) that makes the drag coefficient constant with height. Results also show that the dispersive fluxes are similar in magnitude to the turbulent fluxes, and that they play a very important role within the canopy. Vertical profiles of turbulent length scales (to be used in k–l closure schemes, where k is the turbulent kinetic energy and l a turbulent length scale) are also derived. Finally the distribution of the values around the mean over the reference volumes are analysed for wind and tracer concentrations.     

Computational fluid dynamics modeling of abutment scour under steady current using the level set method

The scour and deposition pattern around an abutment under constant discharge condition is calculated using a three dimensional (3D) Computational Fluid Dynamics (CFD) model. The Reynolds-Averaged Navier Stokes (RANS) equations are solved in three dimensions using a CFD model. The Level Set Method (LSM) is used for calculation of both free surface and bed topography. The two-equation turbulence model (k-ε and k-ω) is used to calculate the eddy viscosity in the RANS equations. The pressure term in the RANS equations on a staggered grid is modeled using the Chorin's projection method. The 5th order Weighted Essentially Non-Oscillatory (WENO) scheme discretizes the convective term of the RANS equations. The Kovacs and Parker and Dey formulations are used for the reduction in bed shear stress on the sloping bed. The model also used the sandslide algorithm which limits bed shear stress reduction during the erosion process. The numerical model solution is validated against experimental results collected at the Politecnico di Milano, Milan, Italy. Further, the numerical model is tested for performance by varying the grid sizes and key parameters like the space and time discretization schemes. The effect of varying bed porosity has been evaluated. Overall, the free surface is well represented in a realistic manner and bed topography is well predicted using the Level Set Method (LSM).