The challenge of tackling global climate change and our increasing reliance on power means that new and diverse renewable energy generation technologies are a necessity for the future. From a number of technologies reviewed at the outset, the cross-flow tidal turbine was chosen as the focus of the research. The numerical investigation begins by choosing to model flow around a circular cylinder as a challenging benchmarking and evaluation case to compare two potential solvers for the ongoing research, ANSYS CFX and OpenFOAM. A number of meshing strategies and solver limitations are extracted, forming a detailed guide on the topic of cylinder lift, drag and Strouhal frequency prediction in its own right. An introduction to cross-flow turbines follows, setting out turbine performance coefficients and a strategy to develop a robust numerical modelling environment with which to capture and evaluate hydrodynamic phenomena. The validation of a numerical model is undertaken by comparison with an experimentally tested lab scale turbine. The resultant numerical model is used to explore turbine performance with varying Reynolds number, concluding with a recommended minimum value for development purposes of Re = 350 × 103 to avoid scalability errors. Based on this limit a large scale numerical simulation of the turbine isconducted and evaluated in detail, in particular, a local flow sampling method is proposed and presented. The method captures flow conditions ahead of the turbine blade at all positions of motion allowing local velocities and angles of attack to be interrogated. The sampled flow conditions are used in the final chapter to construct a novel blade pitching strategy. The result is a highly effective optimisation method which increases peak turbine power coefficient by 20% for only two further case iterations of the numerical solution.

Long linear earthwork assets constructed in high-plasticity overconsolidated clay are known to be deteriorating due to long-term effects of wetting and drying stress cycles as a result of seasonal weather patterns. These stress cycles can lead to shallow first-time failures due to the mobilisation of post-peak strength and progressive failure. Design requirements of new earthworks and management of existing assets requires improved understanding of this critical mechanism; seasonal ratcheting. Incremental model development and validation to allow investigation of multiple inter-related strength deterioration mechanisms of cut slope behaviour in high-plasticity overconsolidated clay slopes has been presented. Initially, the mechanism of seasonal ratcheting has been considered independently and a numerical modelling approach considering unsaturated behaviour has been validated against physical modelling data. Using the validated model, the effects of slope geometry, design parameter selection and design life have been considered. Following this, an approach to allow undrained unloading of soil, stress relief, excess pore water pressure dissipation, seasonal ratcheting and progressive failure with wetting and drying boundary conditions has been considered. Hydrogeological property deterioration and the potential implications of climate change have been explored using the model. In both cases the serviceable life of cut slopes is shown to reduce significantly in the numerical analyses. Finally, a model capable of capturing hydrogeological behaviour of a real cut slope in London Clay has been developed and validated against long-term field monitored data. Using the validated model, a climate change impact assessment for the case study slope has been performed. The numerical analyses performed have indicated that seasonal ratcheting can explain shallow first-time failures in high-plasticity overconsolidated clay slopes and that the rate of deterioration of such assets will accelerate if current climate change projections are representative of future weather.

Este trabalho aborda a influência da hidrodinâmica da Lagoa dos Patos (RS) no comportamento de variáveis indicadoras da qualidade da água, como o material em suspensão, o oxigênio dissolvido, os nutrientes inorgânicos dissolvidos e a clorofila-a. Adicionalmente, o trabalho apresenta uma avaliação do efeito do lançamento de efluentes na magnitude das concentrações destas variáveis e nas suas inter-relações. Para o desenvolvimento deste estudo aplicou-se o modelo numérico MOHID 2D, utilizado com êxito na gestão ambiental de regiões costeiras. Os resultados indicam que a Lagoa dos Patos é um corpo lagunar heterogêneo no qual sua hidrodinâmica é regida principalmente pela atuação dos ventos e da vazão fluvial. Esses fatores influenciam diretamente o comportamento das variáveis estudadas. A presença de efluentes altera significativamente a magnitude das concentrações dos nutrientes e da clorofila-a no ambiente, principalmente nas regiões próximas aos locais de lançamento. A modelagem numérica é uma importante ferramenta para a gestão ambiental de sistemas dinâmicos e que variáveis ecológicas são condicionadas pelas suas inter-relações e pela circulação do sistema. / This work assesses the influence of Patos Lagoon- RS hydrodynamics in the behavior of water quality indicator variables as like suspended material, dissolved oxygen, dissolved inorganic nutrients and chlorophyll-a. Also gives an assessment of the effect of effluent discharge in the magnitude of the concentrations of these variables and in their interrelationships. To develop this study we applied the numerical model MOHID 2D, successfully applied in environmental management of coastal areas. The results indicate that Patos Lagoon is a heterogeneous system in which its hydrodynamic is governed mainly by the action of wind and river flow, which directly influence the behavior of water quality variables. Furthermore, the presence of effluent significantly changes the magnitude of the nutrients and chlorophyll-a concentration, particularly in areas close to the sewage launch sites. We concluded that the numerical modeling is an important tool for the dynamics systems management and that ecological variables are not only conditioned by their interrelations, but also by system circulation.

L’objectif général de cette thèse est de contribuer à l’avancement de la connaissance de la variabilité du Courant Nord Méditerranéen (CN) et de ses interactions avec la dynamique côtière, en s’appuyant principalement sur une configuration numérique réaliste à haute résolution de la façade méditerranéenne française, basée sur le modèle de circulation océanique NEMO et nommée GLAZUR64.La validation de cette configuration avec toutes les observations disponibles sur la période d'étude (CTD, gliders, ARGO, radar HF, ADCP, altimétrie et SST satellite) a permis d'évaluer le réalisme des simulations et leur paramétrisation, et de montrer l'apport de la haute-résolution par rapport aux configurations de bassin au 1/12° utilisées aux frontières de GLAZUR64.Enfin, l'utilisation d'un forçage océanique opérationnel a permis d'utiliser une simulation en complément des données d'une campagne en mer, pour l'étude ciblée d'un tourbillon anticyclonique associé à un méandre du CN au large de Toulon, en avril 2011 [Guihou et al., 2013]. / The main objective of this work is to improve our knowledge of the Northern Mediterrranean Current (NC) variability and its interactions with coastal dynamics, using high-resolution modelling of the North-Western Mediterranean Sea, based on the NEMO code and nammed GLAZUR64.The configuration has been validated using all available observations for the period of study (CTD, gliders, ARGO, HF radar, ADCP, satellite altimetry and SST), in order to assess the realism of the simulations and their parameterisations. The contribution of the high-resolution for the simulation of the mesoscale dynamics, compared to 1/12° configurations used at the boundaries is shown.Finally, the use of an operationnal oceanic boundary forcing allowed us to use a combined data-model approach to study the generation and advection of an anticyclonic eddy trapped at the coast and associated to a NC meander, flowing off Toulon in April 2011 [Guihou et al., 2013].

A numerical study is made of the wintertime circulation within the Southeast Indian Ocean (SEIO). The downwelling favourable winds result in a continuous eastward Coastal Current (CC) extending from Cape Leeuwin to the eastern coast of Tasmania, where it forms a confluence with the south branch of the East Australian Current. An additional forcing mechanism for the CC is the Leeuwin Current in the western part of the domain. The study here is divided in two parts: (1) available data and the wintertime averaged results from the Ocean Circulation and Climate Advanced Model (OCCAM) are analysed to provide a first order description of the large-scale circulation; (2) a high resolution model (Princeton Ocean Model) is nested within OCCAM to examine the shelf-slope circulation within the eastern SEIO. The nested model is forced with climatological monthly average winds and several experiments were run to simulate the effects of surface fluxes of density, enhanced bottom friction and stronger winds. In summary, the shelf-slope circulation is governed by a surface south-eastward CC that carries around 2 Sv and reaches velocities of up to 50 cm/s, where the shelf is narrowest. The core of the current is generally constrained to the shelf-break region. Zonal winds and geostrophic control of the CC lead to a transport of 1 Sv through Bass Strait and a north-eastward jet that is directed into the strait between King Is. and Tasmania. Further south, the CC is poleward and known as the Zeehan Current (ZC). Between Cape Leeuwin and Tasmania and over the slope region, a westward current (the Flinders Current) is found at depths of 500-1000 m and has an associated transport of 5-7 Sv. The current is shown to result from a northward Sverdrup transport in the deep ocean. Meso-scale eddies are shown to result from baroclinic instability and have wavelengths of around 250 km and transports of 3-4 Sv, and can dominate the slope circulation. A comparison of the numerical results is also made with two current meter data sets and results show an interannual variability in the ZC strength, that is probably related to ENSO.

The unsteady nature of coastal hydrodynamics is associated with complex boundary layer dynamics and hence engineering predictions of shear stresses and sediment transport are difficult. This thesis explores some of the complex hydrodynamic problems and boundary layer behaviour in the coastal zone and seeks to provide new and improved modelling approaches. The latest experimental results are used to inform the model development process. New laboratory experiments carried out as part of this thesis illustrate the value of convolution integral calculations for both pressure and skin friction forces on particles on the bed. The experiments also highlight the importance of the phase differences between free stream velocity and boundary layer shear stresses. The use of a “bed” shear stress as a model input is found to be problematic whenever there is a large vertical gradient in the boundary layer shear stress. New experimental and modelling work has helped to improve our understanding of sheet flow boundary layer dynamics. This thesis builds on some of these new discoveries to propose a new simplified model framework for sheet flow sediment transport prediction using convolution integrals. This time domain technique has the advantage of simplicity while incorporating the most important physical processes from more detailed models. The new model framework could be incorporated into any depth averaged coastal hydrodynamic modelling software package. Boundary layer analysis techniques presented in the thesis provide an improved understanding of the effective roughness of mobile beds and can be used to calculate instantaneous shear stress profiles throughout the mobile bed boundary layer. New solutions for swash zone hydrodynamics are presented which illustrate the limitations of the previous benchmark analytical model for swash hydrodynamics. It is shown that real swash necessarily involves a much larger influx of mass and momentum than the analytical solution which was previously used by many in the swash sediment transport research community. Models for swash boundary layer development are also presented.

The main objective of this thesis aims at developing a new generation of software products, in order to obtain a fully automatic simulator predicting the entire Czochralski process while handling correctly the switches between the different growth stages. First of all, new efficient, robust and high-quality automatic mesh generation algorithms with enough flexibility for any complex geometry were implemented, including a 1D mesh generator by global grade-adaptive method, a 2D initial triangulation algorithm by improved sweep line technique and an automatic 2D shape-quality unstructured mesh generator by modified incremental Delaunay refinement technique. Secondly, a Finite Element Navier-Stokes solver based on unstructured meshes was developed and validated. Enhanced turbulence models based on the classical mixing-length or k-l model, together with a generic transformation method to avoid negative k when solving the turbulent kinetic energy equation by the Newton-Raphson iterative method were introduced and implemented. Moreover, laminar and turbulent mathematical models governing the gas convection, thermal distribution and oxygen concentration were developed, and Finite Element numerical methods to solve these governing equations on unstructured meshes were implemented, while appropriate numerical approaches to capture the wall shear stress exerted by the gas flow and experienced by the silicon melt were investigated. Finally, a series of numerical experiments devoted to investigate the industrial Czochralski crystal growth process under various growth conditions are presented based on all the developments implemented. Comparisons of the simulation results with literature and available experimental observations are also presented, and conclusions are drawn based on these simulation results and observations.

Thermoelectric devices are capable of providing both localized active cooling and waste heat power generation. This work will explore the possibility of embedding thermoelectric devices within electronic packaging in order to achieve better system performance. Intel and Nextreme, Inc. have produced thin-film superlattice thermoelectric devices that have above average performance for thermoelectrics and are much thinner than most devices on the market currently. This allows them to be packaged inside of the electronic package where the thermoelectric devices can take advantage of the increased temperatures and decreased thermal lag as compared to the devices being planted on the outside of the package. This work uses the numerical CFD solver FLUENT and the analog electronic circuit simulator SPICE to simulate activity of thermoelectric devices within an electronics package.

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

The mechanical response of geomaterials is highly influenced by geometrical and material heterogeneity. To date, most modelling practices consider heterogeneity qualitatively and the choice of input parameters can be subjective. In this study, a novel approach to combine detailed micro-scale characterization with modelling of heterogeneous geomaterials is presented. The influence of micro-scale heterogeneity and microcracks on the mechanical response and brittle fracture of a crystalline rock was studied using numerical and experimental tools. An existing Combined Finite-Discrete element (FEM/DEM) code was extended to suit heterogeneous, discontinuous, brittle rocks. By conducting grid micro-indentation and micro-scratch tests, the Young's modulus and fracture toughness of the constituent phases of the rock were obtained and used as accurate input parameters for the numerical models. The models incorporated the exact phase mapping obtained from a MicroCT-scanned specimen and the existing microcrack density obtained from thin section analysis. The results illustrated that by incorporating accurate micromechanical input parameters and the intrinsic rock geometric features, the numerical simulations could more accurately predict the mechanical response of the specimen, including the fracture patterns and tensile strength. The numerical simulations illustrated that microstructural flaws such as microcracks should be included in the models to more accurately reproduce the rock strength. In addition, the differential elastic deformations caused by rock heterogeneity altered the stress distribution in the specimen, creating zones of local tensile stresses, in particular, on the boundaries between different mineral phases. As a result, heterogeneous models exhibited rougher fracture surfaces. MicroCT observations emphasized the influence of heterogeneity and, in particular, biotite grains on the fracture trajectories in the specimens. Favourably oriented biotite flakes and cleavage splitting significantly deviated the cracks. The interaction of the main crack with perpendicular cleavage planes of biotite caused strong crack deviation and termination. Considering heterogeneity and the strength degradation caused by microcracks, the simulations captured reasonably accurate mechanical responses and failure mechanisms for the rock, namely, the nonlinear stress-strain relationships. The insights presented in this study improve the understanding of the role of heterogeneity and microstructure on damage and mechanical behaviour of brittle rock.