This study, which includes both experimental and numerical-modelling components, investigates the potentially dangerous conditions that can often occur when low-head dams (or weirs) are overtopped and ‘submerged’-type hydraulic jumps subsequently form downstream of them. The combination of high local turbulence levels, air entrainment, and strong surface currents associated with submerged jumps pose a significant risk to safety of boaters and swimmers. In this study, a wide range of flow regimes and different experimental conditions (i.e. crest length and downstream apron elevation) were considered. The experimental phase involved physical model testing to determine: (i) the hydraulic conditions that govern submerged jump formation, and (ii) the hydrodynamic characteristics of the submerged vortex. The numerical model, developed using OpenFOAM, was validated with the obtained experimental data. This research seeks to help develop improved guidelines for the design and safe operation of low-head dams. The experimental phase of the study involved physical model testing to determine: (i) the hydraulic conditions that govern submerged jump formation, and (ii) the hydrodynamic characteristics of the submerged vortex. The numerical modelling work involved using interFoam (OpenFOAM toolbox) for simulating the experimental results. InterFoam is an Eulerian 3-D solver for multiphase incompressible fluids that employs the Volume of Fluid approach (VOF) to capture the water-air interface. The developed numerical model was subsequently validated using the experimental data collected and processed by the author of this study. The range of tailwater depths associated with submerged hydraulic jump formation is dramatically reduced when a broad-crested weir is coupled with an elevated downstream apron, especially under high flow rate conditions. However sharp-crested weirs induced vortices which displayed reduced velocities and decreased spatial development, which were judged to be safer than broad crest lengths under the same discharge conditions. The classical formulation for the degree of submergence was not explicative when used to evaluate “how submerged the vortex was”. Consequently, a new normalized formulation which compares the local tailwater depth to the lower and upper tailwater limits for the submerged hydraulic jump is proposed. The numerical model developed for this study demonstrated the existence of residual turbulent kinetic energy at downstream sections located within the vortex’s extension, at instants coinciding with the presence of a fully formed roller. This turbulent energy is arguably responsible for the stationary nature of the vortex under constant flow conditions. Residual vertical and horizontal velocities at points located within the vortex’s domain are indicative of the existence of the free surface current.

Natural disasters usually occur without any warning. They can leave trail of destruction and cause much tragedy. We are at a time when we witness fast technological advances; hence, we need to apply the force of scientific advancements to decrease economic losses and the number of human deaths. Tsunami is one of the extreme environmental events that leaves nothing but a path of death and destruction, and as a result, it is essential to understand this phenomenon and identify the mitigation strategies. Several mitigation strategies have been proposed so far; however, more investigations are still required to achieve an acceptable solution. Researchers around the world are studying different aspects of this phenomenon. One of the proposed solutions that has received much attention is designing tsunami-resistant structures which can withstand the force of a tsunami bore. Various studies have been done so far to understand the base shear force of tsunami bore on structures. The focus of this thesis is to improve and better understand the characteristics of the tsunami base shear forces on structures. Hence, in this thesis, two numerical studies were proposed and performed with the main goal of estimating the total tsunami forces on structure under two different conditions. Those include structures with various cross sections, as well as positioning a mitigation wall at an appropriate location relative to the structure. The first study focused on developing a numerical model to study the relationship between tsunami forces and the geometry of the structure. The main goal of this study was to define a numerical model capable of simulating this case precisely. To ensure the accuracy of the model, a comparison was carried out between the results of the numerical model and experimental test performed at the NRC-CHC (National Research Council- Canadian Hydraulics Center) laboratory in Ottawa, Canada and Université Catholique de Louvain (UCL), Belgium, which revealed a very good agreement between the results of the experimental test and numerical model. Further, the validated model was applied to investigate the tsunami force on structures with various cross sections. The second study focus was on developing a numerical model for understanding the role of mitigation wall (a novel idea proposed as a mitigation strategy by the second author of technical paper 2) on reducing the exerted force of tsunami on structures. After developing various models and applying several turbulence models, a valuable result was obtained which demonstrated that a Large Eddy Simulation (LES) model seems to be an excellent approach for predicting the tsunami forces on the structure with a mitigation wall in the direction of the flow. The results of this study will be used to better estimate the tsunami forces exerted on coastal structures which will light the path to the main goal of designing tsunami resistant-structures.

Recent tsunami disasters caused devastating damages to well-engineered coastal infrastructures. In fact, the current design guidelines are not able to provide realistic estimations of tsunami loads in order to design structures to withstand tsunamis. Tsunami hydrodynamic forces are estimated using the drag coefficient. This coefficient is traditionally calculated based on a steady flow analogy. However, tsunami bores behave like unsteady flows. The present work aims at investigating the tsunami forces for different structure geometries to provide realistic guidelines to estimate drag coefficients considering unsteady flows. In the present paper, the dam-break approach is used to investigate the tsunami-like bore interaction with structures. A three-dimensional multiphase numerical model is implemented to study the tsunami induced forces on rectangular shape structures with various aspect ratios (width/depth) and orientations. The numerical model results are validated using measured forces and bore surface elevations of the physical experiments. A scaled-up domain is modeled in order to eliminate the effects of domain sidewalls in the simulation results. The drag coefficient relations with structure geometries and bore depths are provided. The obtained hydrodynamic forces and drag coefficients are compared with existing data in the literature and design codes. For the second topic, a multi-phase three-dimensional numerical reproduction of a large scale laboratory experiment of tsunami-like bores interaction with a surface-piercing circular column is presented. The numerical simulation is conducted in OpenFOAM. The dam-break mechanism is implemented in order to generate tsunami-like bores. The numerical model is validated using the experimental results performed at Canadian Hydraulics Center of the National Research Council (NRC-CHC) in Ottawa. The unsteady Reynolds Averaged Navier-Stokes equations (RANS) are used in order to treat the turbulence effects. The Shear Stress Transport (SST) k-ω turbulence model showed high level of accuracy in replication of the bore-structure interaction. Further, a scaled-up domain is used to investigate the influence of the bed condition in terms of various downstream depths and roughness. Finally, a broad investigation on the bore propagation characteristics is performed. The resulting stream-wise forces exerted on the structural column as well as the bore velocity are compared and analyzed for smooth, rough, dry and wet beds with varying depths.

Vapor generation is central to the flow dynamics within fuel injector nozzles. Because the degree of atomization affects engine emissions and spray characteristics, quantification of phase change within diesel fuel injectors is a topic of design interest. Within the nozzle, the large pressure gradient between the upstream and downstream plena induce large velocities, creating separation and further pressure drop at the inlet corner. When local pressure in the throat drops below the fluid vapor pressure, phase change can occur with sufficient time. At the elevated temperatures present in diesel engines, this process can be dependent upon the degree of superheat, motivating the modeling of heat transfer from the wall. By modeling cavitation and flash boiling phenomena as a departure from equilibrium conditions, the HRMFoam model accurately reproduces canonical adiabatic flows. An experimentally determined relaxation time controls the rate at which vapor is generated, and includes model constants tuned for water and a diesel fuel surrogate. The model is shown to perform well for several benchmark experimental cases, including the work of Reitz, Lichtarowicz, and Nurick. With the implementation of the Farve-averaged energy equation, the present work examines and validates the transport of enthalpy through the fixed heat flux and fixed wall temperature boundary conditions. The pipe heat transfer experiments of Boelter and Allen are replicated using the kEpsilon, Realizable kEpsilon, and Spalart-Allmaras models. With proper turbulence model selection, Allen's heat transfer coefficient data is reproduced within 2.9%. Best-case bulk temperature rise prediction is within 0.05%. Boelter's bulk temperature rise is reproduced within 16.7%. Turbulent diffusivity is shown to determine radial enthalpy distribution.

Liquid wastes discharged from industrial outfalls have been researched for many years in the past. Majority of past studies, initiated in 1960s, were experimental studies mainly focused on basics of discharges such as key geometrical properties. Eventually, more robust experimental studies were performed to measure the mixing properties of effluent discharges with various jet configurations and ambient water conditions. Discharges could be as a means of submerged diffusers or surface channels and receiving water could vary from a homogenous calm ambient to a very complex stratified turbulent cross flow ambient. Depending on the bathymetric and economic situation around an outfall project, submerged discharges are preferred designs for most of ocean outfalls. It is the reason that majority of past studies have evaluated the mixing characteristics of submerged jets. Since early 1990s, the numerical modelling has emerged to support complex fluid mechanic problems. Later in 1990s and early in 2000s, the use of computational fluid dynamic (CFD) tools emerged in predicting the jet properties for the effluent discharges. Since then different numerical models have been developed for different applications. Similar to experimental studies, most of numerical studies have been focused on the submerged dense jet discharges. The current study intends to stay focused on the numerical modelling of such jets too; however, to cover the gaps in the literature. To achieve this, a thorough literature review was performed on the past CFD studies of over past 20 years to better understand what was done and what the gaps are. The results of this thorough review revealed that although there has been a great progress in the CFD studies in the field of effluent discharges, there are some applications that have not been investigated before, yet. It was found that there are some discharge inclinations that were not studied numerically before. Four discharge angles of 60°,75°, 80° and 85° were selected in this study, as previous studies mostly focused on 30° and 45°. The higher inclinations are more suitable for deep water outfalls where terminal rise height of the jet does not attach to the ambient water surface. The numerical model OpenFOAM was used in this study which is based on the Finite Volume Method (FVM) applying LRR turbulence model closure. LRR turbulence models was proved to be a capable choice for effluent discharge modelling. The second gap identified in the comprehensive literature review completed was the submerged dense effluent discharge into shallow water with surface attachment (for both inclined and vertical discharges). There was no previous numerical study of such jets identified. Three different regimes were identified: full submergence, plume contact and centerline impingement regimes (i.e. FSR, PCR and CIR). Key geometrical and dilution properties of these jets at surface contact (Xs, Ss) and return point (Xr, Sr) were extracted numerically and compared to those available from experiments. Two discharge angles (30° and 45°) were investigated based on the available experimental data. Five Reynolds-averaged Navier-Stokes (RANS) turbulence models were examined in this study: realizable k-ε and k-ω SST models (known as two-equation turbulence models), v2f (four equations to model anisotropic behavior) and LRR and SSG turbulence models (known as Reynolds stress models - six equations to model anisotropic behavior). Vertical dense effluent discharges are popular in the design of outfall systems. Vertical jets provide the opportunity to be efficient for a range of ambient currents, where the jet will be pushed away not to fall on itself. This research work investigates worst case scenario in terms of mixing and dilution of such jets: vertical dense effluent discharges with no ambient current and in shallow water where jet impacts the surface. This scenario provides a conservative design criteria for such outfall systems. The numerical modelling of such jets has not been studied before and this research work provides novel, though preliminary, insights in simulations of vertical dense effluent discharges in shallow waters. Turbulent vertical discharges with Froude numbers ranging from 9 to 24 were simulated using a Reynolds stress model (RSM), based on the results from inclined dense discharges to characterize the geometrical (i.e., maximum discharge rise Zm and lateral spread Rsp) and dilution μmin properties of such jets. Three flow regimes were reproduced numerically, based on the experimental data: deep, intermediate and impinging flow regimes.

Tato disertační práce se zabývá experimentálním a numerickým výzkumem kavitace vyvolané rotací. Pro potřeby tohoto výzkumu byla využita transparentní osově symetrická Venturiho dýza, díky čemuž bylo možné zkoumat dynamiku kavitujícího proudění pomocí analýzy vysokorychlostních nahrávek.

No description available.

Effective solvent extraction incorporating electromagnetic heating is a relatively new concept that relies on Radio Frequency heating and solvents to replace steam in current thermal processes for the purpose of extracting bitumen from oil rich sands. The work presented here will further the understanding of the near wellbore flow of this two phase system in order to better predict solvent vaporization dynamics and heat rates delivered to the pay zone. This numerical study details the aspects of phase change of immiscible, two component, liquid/vapor systems confined in porous media heated by electromagnetic radiation, approximated by a spatially dependent volumetric heat source term in the energy equation. The objective of this work is to utilize the numerical methodology presented herein to predict maximum solvent delivery rates to a heated isotropic porous matrix to avoid the over-saturation of the heated pay zone. The total liquid mass content and mean temperature in the domain are monitored to assess whether the liquid phase is fully vaporized prior to flowing across the numerical domain boundary. The distribution of the volumetric heat generation rate used to emulate the physics of electromagnetic heating in the domain decays away from the well bore. Some of the heat generated acts to superheat the already vaporized solvent away from the interface, requiring heat delivery rates that are many times greater than the energy required to turn the liquid solvent to vapor determined by an energy balance. Results of the parametric study from the pay zone simulations demonstrate the importance of the Darcian flow resistance forces added by the porous media to stabilize the flow being pulled away from the wellbore in the presence of gravity. For all cases involving an increase in solvent delivery rate with a constant heat rate, the permeability range required for full vaporization must decrease in order to balance the gravitational forces pulling the solvent from the heated region. For all conditions of permeability and solvent delivery rates, sufficiently increasing the heat rate results in complete vaporization of the liquid solvent. For the case of decreasing solvent delivery rate, a wider range of higher permeabilities for a given heat rate can be utilized while achieving full vaporization. A three dimensional surface outlining the transition from partially vaporized to fully vaporized regimes is constructed relating the solvent delivery rate, the permeability of the porous near wellbore zone and the heat rate supplied to the domain. For the range of permeabilities ~3000mD observed in these types of well bores, low solvent delivery rates and high heat rates must be utilized in order to achieve full vaporization.

Self propulsion modelling is important in order to accurately simulate ships and submarinesusing Computational Fluid Dynamics (CFD). However, fully resolved simulations of hull andpropeller geometries are computationally heavy and time consuming. As such there is a greatinterest in lower order CFD models of propellers. This work investigates three lower ordermodels of a non-cavitating generic submarine propeller (INSEAN E1619) in OpenFOAM. Themodels investigated are Actuator Disk (AD). Rotor Disk (RD) and Actuator Line Model (ALM).The AD model applies a momentum change based on propeller performance coefficients overa disc cell set. The RD uses Blade Element Method (BEM) to calculate a more realistic thrustdistribution over the disk. Finally the ALM applies BEM over seven rotating lines within the cellset disc. The source code to the RD model was modified according to suggestions provided fromearlier studies on the model. The ALM used was originally designed for turbines which wasrectified by changing the force projection vectors in the source code to model propellers instead.There was not enough published data to directly utilize BEM on the E1619 propeller, thus thedata was generated by conducting 2D simulations on every element. The simulations were setup to replicate results provided in earlier works with higher order models in order to compareboth quantitative and qualitative results. It was found the ALM matched the reference databest out of the models tested in this work. The RD was qualitatively similar to the time averageof the ALM fields but numerically inaccurate. The AD results were poor, both quantitativelyand qualitatively.

In the field of heavy-duty vehicles, fuel efficiency and environmental protection are factors that need to be focused on, while the aerodynamic drag generated during vehicle travelling is one of the most influential aspects. This thesis delves into the aerodynamic simulation of Scania trucks using the open-source Computational Fluid Dynamics (CFD) tool, OpenFOAM v2206. This study rigorously investigates the aerodynamics of two Scania truck models under different operating conditions, including scenarios with different crosswind environments at high speeds.The core of this study is to compare and analyse the computational results of OpenFOAM v2206 and its predecessor OpenFOAM v3.0+ in a number of aspects, in order to elucidate the evolution and improvement of CFD techniques and their practical impact on vehicle simulation performance. In order to save computational resources, the RANS method was used for the steady-state simulations. Preliminary comparisons were also made with results from PowerFLOW, another CFD software widely used within the Scania group.Another important part of this thesis is the exploration of an alternative meshing method (ANSA Hextreme Mesh) in CFD simulations. As a widely used pre-processing software in the Scania group today, analysing and comparing the advantages and disadvantages of ANSA and OpenFOAM in terms of meshing, such as the ease of meshing and the accuracy of aerodynamic predictions, can help to provide valuable guidance for the application of truck shape design and aerodynamic simulation.The results indicate that OpenFOAM v2206 excels in predicting aerodynamics and has utility in optimising truck design. Compared to OpenFOAM v3.0+, OpenFOAM v2206 shows smaller discrepancies in results with PowerFLOW. Further exploration is required regarding transient simulations using OpenFOAM. In terms of meshing methods, a simplified model (Allan Body) was investigated, and there is further research to be done on meshing the complete truck.In conclusion, this thesis presents a comprehensive and in-depth exploration of truck aerodynamics using advanced CFD tools. The results not only deepen the understanding of airflow dynamics around heavy vehicles, but also pave the way for the development of more aerodynamically efficient and environmentally friendly truck designs.