Bubbly flows around ships have been studied for years, mostly in relation with ship acoustic signatures. Bubbles are generated at the bow and shoulder breaking waves, at the hull/free surface contact line, the propeller and the highly turbulent stern flow. These bubbles are further transported downstream by the flow forming a two-phase mixture in the wake that can be kilometers long. The presence of bubbles in the wake of a ship significantly affects the acoustic response of the medium and can be detected by measuring acoustic attenuation and backscattering making a ship vulnerable to detection. Additionally, the bubbly wake shows at the surface as a characteristic signature of white water, and given the length of the bubbly wake, it makes a ship visible from satellites. Therefore, the bubbly wake can be used to detect and identify surface ships. Bubbly flows do not scale to model scale experiments, and experiments on full scale ships are scarce mostly due to difficult access areas and the high speeds involved. It is therefore of interest to simulate the bubbly flow around ships to provide information difficult, if not impossible, to obtain with experiments. This work presents the development of a code for the simulation of polydispersed bubbly flows with a focus on ship hydrodynamics. The mathematical model implemented is based on a two-fluid formulation coupled with a Boltzmann-like transport equation describing the bubbly phase. The tool developed attempts to include most of the relevant physics of the problem to represent better the conditions of real scenarios. The resulting code allows the simulation of polydispersed bubbly flows in situations including free surface and air entrainment, high void fraction levels and moving control surfaces and propulsors. The code is two-way coupled, with a strong coupling between the two phases and between the bubble sizes. The complexity of the problems tackled in this research required the development of novel numerical methods solving issues never identified before or simply neglected. These methods play an essential role in the accuracy, robustness and efficiency of the code and include: a two-phase projection method that not only couples pressure and velocity but also implicitly couples void fraction, a time splitting marching scheme to solve separately coupling in space and in bubble sizes, and a stable numerical method to integrate the strong coupling introduced by collision forces. The implemented code is applied to the simulation of the bubbly flow around a full scale ship using the latest available models and computational techniques. A study is performed on the influence of several mechanisms on the predicted bubbly wake and comparisons with available experimental data are presented. The influence of breakup in the boundary layer is analyzed in detail as well. In addition, this work identifies several modeling and implementations issues and attempts to provide a path for future studies. To illustrate the flexibility and robustness of the code, a final demonstration case is presented that includes rotating propellers. The computation is performed at full scale, with the fully appended geometry of the vessel and includes incoming waves, oceanic background and rectified diffusion models. Many of these features are unique to this computation and make it the first of its kind.

A high fidelity computational fluid dynamics approach to perform direct simulations of ship maneuvers is presented in this thesis. The approach uses dynamic overset grids with a hierarchy of bodies to enable arbitrary motions between objects, and overcome the difficulties in simulation of the moving rudder and rotating propeller. To better resolve propeller/rudder interaction a Delayed Detached Eddy Simulation turbulence model based on Menter’s SST is used. The methodology was implemented in the general purpose RANS/DES/DDES research code REX, and is applied to the KRISO Container Ship (KCS) with moving rudder and rotating propeller in deep and shallow water. For the first time, a grid study is conducted for the self-propulsion condition for the propeller RPM, thrust, torque and lateral force, and for the roll and pitch motions, using grids of 8.7 (coarse), 24.6 (medium) and 71.3 (fine) million points. A grid study is also performed for the zigzag maneuver evaluating the maximum and minimum values of propeller thrust, torque and lateral force roll, pitch, yaw, roll rate, yaw rate and drift throughout the maneuver. An extensive comparison between predicted motions and forces of the direct simulations and the experimental data collected by Schiffbau-Versuchsanstalt Potsdam GmbH (SVA) and Flanders Hydraulics Research (FHR) are presented. While the results and comparisons with experimental data show that using direct CFD to compute modified and standard maneuvers with moving rudder and rotating discretized propeller is feasible, computational cost remains an impediment for many practical applications. Coupling a dynamic overset CFD solver with a potential propeller code can dramatically reduce the computational time to perform maneuvering simulations by using one order of magnitude larger time step than direct simulation. This thesis investigates the ability of a coupled CFD/potential propeller code approach to simulate maneuvers in ships, where the rudder is located downstream of the propeller. While the approach has been successfully applied to submarine maneuvers, in which the propeller wake is free of interference, the concept had not been evaluated before for cases where an object (the rudder) is immersed in the wake. The study is performed using the CFD code REX and the propeller code PUF-14. Performance of the coupled REX/PUF-14 approach is first tested studying propeller/rudder interaction, evaluating influence of the propeller/rudder gap size and rudder deflection on propeller performance curves and rudder forces, comparing against DDES simulations with a discretized rotating propeller. A grid study was performed for advance coefficient J=0.6 and a rudder angle δ=20 degrees for a propeller rudder gap of 0.2 times the rudder radius, with the resulting grid uncertainties for propeller thrust and torque coefficients suggesting that the effects of the grid changes are small for the present range of grid sizes. A 15/1 zigzag maneuver for the KCS container ship, in which case the rudder is very close downstream of the propeller, is then analyzed, and compared against discretized propeller simulations and experimental data. Self-propulsion coupled REX/PUF-14 results agree very well with experiments and discretized propeller simulations. Prediction of motions, forces and moments, and mean flow field with the coupled REX/PUF-14 approach are comparable to results obtained with discretized propeller simulations and agree with experiments well, though as implemented the coupled approach is unable to resolve tip vortices and other flow structures that interact with the rudder, potentially affecting prediction of flow separation. It can be concluded that coupled CFD/potential flow propeller approaches are an effective and economical way to perform direct simulation of surface ship maneuvers with CFD.

The Upper Mississippi River System (UMRS) is a diverse and dynamic ecosystem that includes the main stem river channel, side channels, backwater floodplains and lakes, islands, wetlands, grasslands, and floodplain forests. The hydrology of this rich ecosystem is one of the key drivers for physical, chemical and biological processes. However, the hydrology and hydraulics of the UMRS has been drastically altered from its natural state as a result of the construction of the locks and dams in the 1930s. Beginning with the Water Resources Development Act of 1986, biologists, ecologists, and engineers have been working to restore the river to a more natural state within the current constraints imposed by the lock and dam system. In an effort to restore rivers to a more natural state, the determination of a hydraulic reference condition is essential to understanding the "why and how" of historical river system function. Understanding the fundamental processes of historical conditions will help prioritize resources and better quantify possible outcomes for riverine restoration. The main goal of this study was to construct a hydrodynamic reference condition for Pool 8 of the Upper Mississippi River System using hydrodynamic computational fluid dynamic (CFD) modeling. The CFD model will provide a better understanding of pre-impoundment flow conditions as compared to post-impoundment conditions today. The numerical model was constructed and developed primarily from a pre-impoundment 1890s topographic map with bathymetric cross-sections in the channels. The 1890s map and other sources from the U.S. Army Corps of Engineers provided historic elevation and hydraulic reference data for model calibration. The calibrated historic model was then compared with a current model of similar scale representing post-impoundment conditions, allowing for quantitative analysis of the differences between the two conditions. Model results indicated large changes in average depth and average velocity between historic and current conditions in certain parts of the pool, while others remained relatively unchanged. For example, velocities decreased in main channel aquatic areas in the lower part of Pool 8 from an average of 0.6 m/s (2.0 ft/s) under historic conditions to 0.1 m/s (0.3 ft/s) under current conditions. In the same part of the pool, however, velocities in contiguous backwater areas remained relatively constant, with most remaining less than 0.25 m/s (0.82 ft/s). Additionally, in the lower part of the pool, discharge distribution between the floodplain areas and the main channel was historically much more dynamic, with flow concentrated in the main and secondary channels at discharges less than 2265 m3/s and in the floodplains at greater than 2265 m3/s. Under current conditions, discharge distribution is much less dynamic, with approximately 2/3 of the total discharge conveyed on the floodplain for all discharges modeled (283 m3/s to 2832 m3/s or 10,000 ft3/s to 100,000 ft3/s).

High-fidelity CFD-MBD FSI (Computational Fluid Dynamics - Multi Body Dynamics Fluid-Structure Interaction) code development and validation by full-scale experiments is presented, for a novel hull form, WAM-V (Wave Adaptive Modular Vessel). FSI validation experiments include cylinder drop with suspended mass and 33 ft WAM-V sea-trials. Calm water and single-wave sea-trails were with the original suspension, while the rough-water testing was with a second generation suspension. CFDShip-Iowa is used as CFD solver, and is coupled to Matlab Simulink MBD models for cylinder drop and second generation WAM-V suspension. For 1DOF cylinder drop, CFD verification and validation (V&V) studies are carried out including grid and time-step convergence. CFD-MBD results for 2DOF cylinder drop show that 2-way coupling is required to capture coupled physics. Overall, 2-way results are validated with an overall average error value of E=5.6%DR for 2DOF cylinder drop. For WAM-V in calm water, CFD-MBD 2-way results for relative pod angle are validated with E=14.2%DR. For single-wave, CFD-MBD results show that 2-way coupling significantly improves the prediction of the peak amplitude in pontoon motions, while the trough amplitudes in suspension motions are under-predicted. The current CFD-MBD 2-way results for single-wave are validated with E=17%DR. For rough-water, simulations are carried out in regular head waves representative of the irregular seas. CFD-MBD 2-way results are validation with E=23%D for statistical values and the Fourier analysis results, which is reasonable given the differences between simulation waves and experiments.

The delivery of nutrients to the surface of marine algae can be controlled by the local hydrodynamic regime: in higher flow velocities, the Diffusive Boundary Layer (DBL) at the uptake surface is thinner, which can increase the flux of dissolved chemicals to the algal surface. If the primary productivity of an alga is controlled by the availability of a dissolved chemical, increased water flow should result in greater primary productivity due to increased chemical flux. To test the hypothesis that increased water flow will increase Photosystem II kinematics (PSII) in the green alga Caulerpa we used a Diving Pam Fluorometer to measure the maximum relative electron transport rate (Pmax), Saturation Irradiance (Ik), Non-photochemical quenching (NPQ), the light limited slope of photosynthesis vs. irradiance curve (α) and photo-chemical quenching (qP) and compared these measured values among treatments of varying flow speeds in a portable laboratory flume. We also measured the influence of water flow on values of Pmax, Ik, α , qP and NPQ in the field. Results showed that in C. racemosa collected from Tampa bay, and tested in a laboratory flume, values of Pmax and Ik were positively correlated to increase water flow, possibly indicating mass-transfer limitation. C. mexicana, collected from the Florida Keys, showed a decrease in values of Pmax, and Ik with increasing water velocity in flume experiments, indicating that the increased flow was resulting in physiological stress. This result was supported with field measurements for C. sertularioides, which showed a negative correlation between Pmax and flow velocity and Ik and flow velocity.

Three-dimensional scour propagation along offshore pipelines is a major reason to pipeline failures in an offshore environment. Although the research on scour in both numerical and experimental aspect has been extensive over the last three decades, the focus of the investigation has been limited to the two-dimensional aspect. The knowledge on three-dimensional scour is still limited. This dissertation presents the results of an experimental investigation on threedimensional scour along offshore pipelines in (1) steady currents (2) waves only and (3) combined waves and current. The major emphasis of the investigation is to investigate the propagation of the scour hole along the pipeline after the initiation of scour. Physical experiments conducted were used to quantify the effects of various parameters on scour propagation velocities along the pipeline. The problem of monitoring real time scour below a pipeline was solved by using specifically developed conductivity scour probes. Effects of various parameters such as pipeline embedment depth, incoming flow Shields parameter, Keuglegan- Carpenter (KC) number and flow incident angle to the pipeline on scour propagation velocities along the pipeline were investigated. The investigations clearly reveal that scour propagation velocities generally increase with the increase of flow but decrease with the increase of the pipeline embedment depth. A general predictive formula for scour propagation velocities is proposed and validated against the experimental results. There are still some common issues related to pipeline scour that is lacking in the literature to date. One of these issues is the effects of Reynolds number on two-dimensional scour beneath pipelines. A numerical approach was adopted to investigate the Reynolds-number dependence of two-dimensional scour beneath offshore pipelines in steady currents. A novel wall function is proposed in calculating the suspended sediment transport rate in the model. The effects of Reynolds number were investigated by simulating the same undisturbed Shields parameters in both model and prototype but with different values of Reynolds number in two separate calculations. The results revealed that scour depths for prototype pipelines are about 10~15% smaller than those for model pipelines. The normalized time scales was found to be approximately the same, and the simulated scour profiles for the model pipelines agree well with the experimental results from an independent study. The backfilling of pipeline trenches is also an important issue to the design and management of offshore pipelines. A numerical model is developed to simulate the self-burial of a pipeline trench. Morphological evolutions of a pipeline trench under steady-current or oscillatory-flow conditions are simulated with/without a pipeline inside the trench. The two-dimensional Reynolds-averaged continuity and Navier-Stokes equations with the standard k-e turbulence closure, as well as the sediment transport equations, are solved using finite difference method in a curvilinear coordinate system. Different time-marching schemes are employed for the morphological computation under unidirectional and oscillatory conditions. It is found that vortex motions within the trench play an important role in the trench development.

Abstract not available

Abstract not available

Abstract not available

A weakly nonlinear stability analysis was conducted for the flow induced in an incompressible, Newtonian, viscous fluid lying between two infinite parallel plates which form a channel. The plates are oscillating synchronously in simple harmonic motion. The disturbed velocity of the flow was written in the form of a series in powers of a parameter which is a measure of the distance away from the linear theory neutral conditions. The individual terms of this series were decomposed using Floquet theory and Fourier series in time. The equations at second order and third order in were derived, and solutions for the Fourier coefficients were found using pseudospectral methods for the spatial variables. Various alternative methods of computation were applied to check the validity of the results obtained. The Landau equation for the amplitude of the disturbance was obtained, and the existence of equilibrium amplitude solutions inferred. The values of the coefficients in the Landau equation were calculated for the nondimensional channel half-widths h for the cases h = 5, 8, 10, 12, 14 and 16. It was found that equilibrium amplitude solutions exist for points in wavenumber Reynolds number space above the smooth portion of the previously determined linear stability neutral curve in all the cases examined. Similarly, Landau coefficients were calculated on a special feature of the neutral curve (called a ???finger???) for the case h = 12. Equilibrium amplitude solutions were found to exist at points inside the finger, and in a particular region outside near the top of the finger. Traces of the x-components of the disturbance velocities have been presented for a range of positions across the channel, together with the size of the equilibrium amplitude at these positions. As well, traces of the x-component of the velocity of the disturbed flow and traces of the velocity of the basic flow have been given for comparison at a particular position in the channel.