An experimental research was performed to study the effects of surface roughness on the characteristics of separated and reattached turbulent flows in an open channel. A backward facing step was used to induce flow separation. The rough surfaces comprised wire mesh grit-80 and sand grains of average diameter 1.5 mm. In each experiment, the Reynolds number based on the step height and freestream velocity of approach flow was fixed at 3240 and the Reynolds number based on the approach flow depth and freestream velocity was kept constant at 25130. Particle image velocimetry (PIV) technique was used to measure the flow velocity. The results showed that roughness effects on the mean and turbulent quantities are evident only in the recovery region. Moreover, roughness effects on the flow dynamics are dependent on the specific roughness element.

This thesis aimed at developing multi-scale models and computational methods for aerother-modynamics applications. The research on multi-scale models has focused on internal energy excitation and dissociation of molecular gases in atmospheric entry flows. The scope was two-fold: to gain insight into the dynamics of internal energy excitation and dissociation in the hydrodynamic regime and to develop reduced models for Computational Fluid Dynamics applications. The reduced models have been constructed by coarsening the resolution of a detailed rovibrational collisional model developed based on ab-initio data for the N2 (1Σ+g)-N (4Su) system provided by the Computational Quantum Chemistry Group at NASA Ames Research Center. Different mechanism reduction techniques have been proposed. Their appli-cation led to the formulation of conventional macroscopic multi-temperature models and vi-brational collisional models, and innovative energy bin models. The accuracy of the reduced models has been assessed by means of a systematic comparison with the predictions of the detailed rovibrational collisional model. Applications considered are inviscid flows behind normal shock waves, within converging-diverging nozzles and around axisymmetric bodies, and viscous flows along the stagnation-line of blunt bodies. The detailed rovibrational colli-sional model and the reduced models have been coupled to two flow solvers developed from scratch in FORTRAN 90 programming language (SHOCKING_F90 and SOLV-ER_FVMCC_F90). The results obtained have shown that the innovative energy bin models are able to reproduce the flow dynamics predicted by the detailed rovibrational collisional model with a noticeable benefit in terms of computing time. The energy bin models are also more accurate than the conventional multi-temperature and vibrational collisional models. The research on computational methods has focused on rarefied flows. The scope was to formu-late a deterministic numerical method for solving the Boltzmann equation in the case of multi-component gases with internal energy by accounting for both elastic and inelastic collisions. The numerical method, based on the weighted convolution structure of the Fourier trans-formed Boltzmann equation, is an extension of an existing spectral-Lagrangian method, valid for a mono-component gas without internal energy. During the development of the method, particular attention has been devoted to ensure the conservation of mass, momentum and en-ergy while evaluating the collision operators. Conservation is enforced through the solution of constrained optimization problems, formulated in a consistent manner with the collisional in-variants. The extended spectral-Lagrangian method has been implemented in a parallel com-putational tool (best; Boltzmann Equation Spectral Solver) written in C programming lan-guage. Applications considered are the time-evolution of an isochoric gaseous system initially set in a non-equilibrium state and the steady flow across a normal shock wave. The accuracy of the proposed numerical method has been assessed by comparing the moments extracted from the velocity distribution function with Direct Simulation Monte Carlo (DSMC) method predictions. In all the cases, an excellent agreement has been found. The computational results obtained for both space homogeneous and space inhomogeneous problems have also shown that the enforcement of conservation is mandatory for obtaining accurate numerical solutions.

Spatial interaction models of the gravity type are used in conjunction with sample data on flows between origin and destination locations to analyse international and interregional trade, commodity, migration and commuting patterns. The focus is on the classical log-normal model version and spatial econometric extensions that have recently appeared in the literature. These new models replace the conventional assumption of independence between origin-destination flows with formal approaches that allow for spatial dependence in flow magnitudes. The paper also discusses problems that arise in applied practice when estimating (log-normal) spatial interaction models. (authors' abstract)

It is well established that non-equilibrium flows cannot properly be described by traditional hydrodynamics, namely, the Navier-Stokes-Fourier (NSF) equations. Such flows occur, for example, in micro-electro-mechanical systems (MEMS), and ultra vacuum systems, where the dimensions of the devices are comparable to the mean free path of a gas molecule. Therefore, the study of non-equilibrium effects in gas flows is extremely important. The general interest of the present study is to explore boundary value problems for moderately rarefied gas flows, with an emphasis on numerical solutions of the regularized 13--moment equations (R13). Boundary conditions for the moment equations are derived based on either phenomenological principles or on microscopic gas-surface scattering models, e.g., Maxwell's accommodation model and the isotropic scattering model. Using asymptotic analysis, several non-linear terms in the R13 equations are transformed into algebraic terms. The reduced equations allow us to obtain numerical solutions for multidimensional boundary value problems, with the same set of boundary conditions for the linearized and fully non-linear equations. Some basic flow configurations are employed to investigate steady and unsteady rarefaction effects in rarefied gas flows, namely, planar and cylindrical Couette flow, stationary heat transfer between two plates, unsteady and oscillatory Couette flow. A comparison with the corresponding results obtained previously by the DSMC method is performed. The influence of rarefaction effects in the lid driven cavity problem is investigated. Solutions obtained from several macroscopic models, in particular the classical NSF equations with jump and slip boundary conditions, and the R13--moment equations are compared. The R13 results compare well with those obtained from more costly solvers for rarefied gas dynamics, such as the Direct Simulation Monte Carlo (DSMC) method. Flow and heat transfer in a bottom heated square cavity in a moderately rarefied gas are investigated using the R13 and NSF equations. The results obtained are compared with those from the DSMC method with emphasis on understanding thermal flow characteristics from the slip flow to the early transition regime. The R13 theory gives satisfying results including flow patterns in fair agreement with DSMC in the transition regime, which the conventional Navier-Stokes-Fourier equations are not able to capture. / Graduate / 0548 / anirudh@uvic.ca

The present research focuses on the development and computer implementation of a novel threedimensional, anisotropic turbulence model not only capable of handling complex geometries but also the turbulence driven secondary currents. The model equations comprise advanced algebraic Reynolds stress models in conjunction with Reynolds Averaged Navier-Stokes equations. In order to tackle the complex geometry of compound meandering channels, the body-fitted orthogonal coordinate system is used. The finite volume method with collocated arrangement of variables is used for discretization of the governing equations. Pressurevelocity coupling is achieved by the standard iterative SIMPLE algorithm. A central differencing scheme and upwind differencing scheme are implemented for approximation of diffusive and convective fluxes on the control volume faces respectively. A set of algebraic equations, derived after discretization, are solved with help of Stones implicit matrix solver. The model is validated against standard benchmarks on simple and compound straight channels. For the case of compound meandering channels with varying sinuosity and floodplain height, the model results are compared with the published experimental data. It is found that the present method is able to predict the mean velocity distribution, pressure and secondary flow circulations with reasonably good accuracy. In terms of engineering applications, the model is also tested to understand the importance of turbulence driven secondary currents in slightly curved channel. The development of this unique model has opened many avenues of future research such as flood risk management, the effects of trees near the bank on the flow mechanisms and prediction of pollutant transport.

Gradient flows of energy functionals on the space of probability measures with Wasserstein metric has proved to be a strong tool in studying certain mass conserving evolution equations. Such gradient flows provide an alternate formulation for the solutions of the corresponding evolution equations. An important condition, which is known to guarantees existence, uniqueness, and continuous dependence on initial data is that the corresponding energy functional be displacement convex. We introduce a relaxed notion of displacement convexity and we show that it still guarantees short time existence and uniqueness of Wasserstein gradient flows for higher order energy functionals which are not displacement convex in the standard sense. This extends the applicability of the gradient flow approach to larger family of energies. As an application, local and global well-posedness of different higher order non-linear evolution equations are derived. Examples include the thin-film equation and the quantum drift diffusion equation in one spatial variable.

This thesis investigates why and how macroeconomic volatility in source countries interacts with their FDI outflows. The study focuses on FDI flowing out from OECD countries to less developed countries in the ASEAN region. Using a panel data encompassing 52 country-pairs over the period 1996-2011, I find a negative correlation between FDI outflows and macroeconomic volatility in source countries. More specifically the empirical results suggest an adverse relationship between inflation and output volatility (business cycles fluctuations) and FDI flows – the more macroeconomic volatility in developed economies the lesser FDI flows to less developed economies, which is explained by Keynesian theories. These findings derive from a gravity model approach, which enabled me to control for host country determinants. In order to estimate these relationships I adopted a random effects model and a tobit model. The reason behind the use of these two models derives from the different views within this branch of research because of censored FDI statistics. The thesis is inspired by Éric Rougier’s et al. work on how macroeconomic volatility in European countries interacts with FDI flows to the MENA region (2012).

Many problems in fluid mechanics and material sciences deal with liquid-vapour flows. In these flows, the ideal gas assumption is not accurate and the van der Waals equation of state is usually used. This equation of state is non-convex and causes the solution domain to have two hyperbolic regions separated by an elliptic region. Therefore, the governing equations of these flows have a mixed elliptic-hyperbolic nature. Numerical oscillations usually appear with standard finite-difference space discretization schemes, and they persist when the order of accuracy of the semi-discrete scheme is increased. In this study, we propose to use a Chebyshev pseudospectral method for solving the governing equations. A comparison of the results of this method with very high-order (up to tenth-order accurate) finite difference schemes is presented, which shows that the proposed method leads to a lower level of numerical oscillations than other high-order finite difference schemes, and also does not exhibit fast-traveling packages of short waves which are usually observed in high-order finite difference methods. The proposed method can thus successfully capture various complex regimes of waves and phase transitions in both elliptic and hyperbolic regimes

I studied how geology, land use and rainfall, correlated with peak flow responses in 110 headwater stream sites during a drought year. Highest discharges were observed in the most developed catchments and in the most poorly drained soils, but specific responses were variable depending on both geology and land disturbance. Redundancy analysis indicated that both surficial geology and land disturbance were important predictors of discharge and that rainfall was in general a poor predictor of discharge. I conclude that responses of headwater streams to individual storms are unpredictable from data generated using GIS, but increased peak flows occur associated with human development, mitigated by surficial geology. The headwater streams that are most vulnerable to flow alterations occur on poorly drained soils, and where urbanization tends to concentrate. Much greater attention to managing water is required if further degradation of stream ecosystems is to be prevented from our future land use.

In this thesis, chemically reacting flows are studied mainly for detonation problems under 1D, cylindrical and spherical symmetry conditions. The mathematical formulation of chemically reacting, inviscid, unsteady flows with species conservation equations and finite-rate chemistry is described. The Euler equations with finite-rate chemistry are discretized by Finite-Volume method and solved implicitly by using a time-spliting method. Inviscid fluxes are computed using Roe Flux Difference Splitting Model. The numerical solution is implemented in parallel using domain decomposition and PVM library routines for inter-process communication. The solution algorithm is validated first against the numerical and experimental data for a shock tube problem with and without chemical reactions and for a cylindrical and spherical propagation of a shock wave. 1D, cylindrically and spherically symmetric detonations of H2:O2:Ar mixture are studied next.