In the field of ocean wave energy converters (WECs), active areas of research are on a priori or in situ methods for power production estimates and on control system design. Linear potential flow theory modelling techniques often underpin these studies; however, such models rely upon small wave and body motion amplitude assumptions and therefore cannot be applied to all wave conditions. Nonlinear extensions can be applied to the fluid loads upon the structure to extend the range of wave conditions for which these models can provide accurate predictions. However, careful consideration of the thresholds of wave height and periods to which these models can be applied is still required. Experimental modelling in wave tank facilities can be used for this purpose by comparing experimental observations to numerical predictions using the experimental wave field as an input. This study establishes a recommended time domain numerical modeling approach for power production assessments of oscillating wave surge converters (OWSCs), a class of WEC designed to operate in shallow and intermediate water depths. Three candidate models were developed based on nonlinear numerical modelling techniques in literature, each with varying levels of complexity. Numerical predictions provided by each model were found to be very similar for small wave amplitudes, but divergence between the models was observed as wave height increased. Experimental data collected with a scale model OWSC for a variety of wave conditions was used to evaluate the accuracy of the candidate models. These experiments were conducted in a small-scale wave flume at the University of Victoria. A challenge with this experimental work was managing wave reflections from the boundaries of the tank, which were significant and impacted the dynamics of the scale model OWSC. To resolve this challenge, a modified reflection algorithm based upon the Mansard and Funke method was created to identify the incident and reflected wave amplitudes while the OWSC model is in the tank. Both incident and reflected wave amplitudes are then input to the candidate models to compare numerical predictions with experimental observations. The candidate models agreed reasonably well with the experimental data, and demonstrated the utility of the modified wave reflection algorithm for future experiments. However, the maximum wave height generated in the wave tank was found to be limited by the stroke length of the wavemaker. As a result, no significant divergence of the candidate model predictions from the experimental data could be observed for the limited range of wave conditions, and therefore a recommended model could not be selected based solely on the experimental/numerical model comparisons. Preliminary assessments of the annual power production (APP) for the OWSC were obtained for a potential deployment site on the west coast of Vancouver Island. Optimal power take-off (PTO) settings for the candidate models were identified using a least-squares optimization to maximize power production for a given set of wave conditions. The power production of the OWSC at full scale was then simulated for each bin of a wave histogram representing one year of sea states at the deployment site. Of the three candidate models, APP estimates were only obtained for Model 1, which has the lowest computational requirements, and Model 3, which implements the most accurate algorithm for computing the fluid loads upon the OWSC device. Model 2 was not considered as it provides neither advantages of Models 1 and 3. The APP estimates from Models 1 and 3 were 337 and 361 MWh per year. For future power production assessments, Model 3 is recommended due to its more accurate model of the fluid loads upon the OWSC. However, if the high computational requirements of Model 3 are problematic, then Model 1 can be used to obtain a slightly conservative estimate of APP with a much lower computational effort. / Graduate

Thesis (MEng)--Stellenbosch University, 2015. / ENGLISH ABSTRACT: In the present study, two-phase natural circulation flow in a multi-parallel channel system was investigated using experimental and numerical modelling. The investigation was carried out under different power excitations and various system operations (open system, closed system and heat pipe mode). The multichannel system was equipped at the upper end with a condenser enclosed within a steam drum, while the lower portion of each channel was heated to heat the system. For the numerical modelling, transient one-dimensional conservation equations were derived from first principle for both single- and two-phase fluids and used to computer program the system’s discretised simulation model. Temperatures and mass flow rates of the fluid responses as a result of different power excitations and operations were obtained for both the experimental and numerical modelling. It was observed from the results that the fluid experiences a start-up transient before accomplishing steady-state conditions. It was further noted that the transient duration varies with power excitations and system operation modes and hence with the stability of the system. A rise in power proved not to necessarily increase the fluid mass flow rate, but invited oscillations with higher amplitudes, depending on the system’s mode of operation. Type I instability and low-quality steam oscillations were witnessed at low power and open system operation mode (system open to the atmosphere). Type II instabilities and flashing instability were observed to be associated with medium and high-power excitations for the open system mode of operation. The fluid flow became more stable and less oscillatory at all power excitations for the closed system operation mode (system not open to the atmosphere). However, a sub-cooling effect was evident at higher power, where the two-phase fluid temperatures oscillated in a sinusoidal manner. However, the mass flow rates oscillated with high amplitudes in the forward direction in some channels and assumed a unidirectional flow in other channels. In general, steady-state conditions were obtained earlier when the system was operated as a closed system. For the heat pipe mode of operation, the system transient response in all channels exhibited a geysering instability followed by flashing-induced boiling. In-phase (flow in channels exhibiting the same behaviour) and out-of-phase (flow in channels exhibiting contradictory conduct) behaviour between adjacent channels were observed at all power excitations and system operation modes. Flow reversal in heated channels of a natural circulation system were proven to exists even under equal power excitations. / AFRIKAANSE OPSOMMING: In hierdie studie is tweefasige natuurlike sirkulasievloei in ’n parallelle multikanaalstelsel ondersoek deur middel van eksperimentele en numeriese modellering. Die ondersoek is onder verskillende kragopwekkings en verskeie stelselwerkings (oop stelsel, toe stelsel en hittepypmodus). Die multikanaalstelsel is aan die bopunt met ’n kondensor binne ’n stoomdrom toegerus, terwyl die laer gedeelte van elke kanaal verhit is om die stelsel te verhit. Vir die numeriese modellering, is oorgangseendimensionele behoundsvergelykings vanaf die eerste beginsel vir beide een- en tweefasige vloeistowwe afgelei en dit is gebruik om die stelsel se gediskretiseerde simulasiemodel vir ’n rekenaar te programmeer. Temperature en massavloeitempo’s van die vloeistofrespons as gevolg van verskillende kragopwekkings en -werkings is vir beide die eksperimentele en die numeriese modellering verkry. Dit is in die resultate waargeneem dat die vloeistof ’n aansitoorgang ervaar voor dit vloeiewewigstoestande bereik. Daar is verder waargeneem dat die duur van die oorgang wissel volgens kragopwekkings en stelselwerkingsmodusse en dus op grond van die stabiliteit van die stelsel. ’n Toename in krag het nie noodwendig die vloeitempo van die vloeistofmassa verhoog nie, maar het aanleiding gegee tot ossillasies met groter amplitudes, afhangende van die stelsel se metode van werking. Tipe I-onstabiliteit en stoom-ossillasies van ’n lae intensiteit is teen lae krag en oop stelselwerkingsmodus waargeneem (stelsel oop aan die atmosfeer). Tipe II-onstabiliteit en flitsingsonstabiliteit (flashing instability) is met medium- en hoë kragopwekkings vir die oop stelsel modus van werking waargeneem. Die vloeistofvloei het meer stabiel en minder ossillerend geraak by alle kragopwekkings in die geslote stelsel van werking (stelsel nie oop na die atmosfeer nie). ’n Subverkoelingseffek was egter teen hoër krag duidelik, waar die tweefasige vloeistof se temperature sinusvormig geossilleer het. Die massavloeitempo’s het egter met hoë amplitudes in die vorentoe rigting in sommige kanale gevloei en eenrigtingvloei in ander kanale vertoon. Oor die algemeen is vloei-ewewigstoestande vroeër verkry toe die stelsel as ’n geslote stelsel bedryf is. Vir die hittepypmodus van werking het die stelsel se oorgangsweergawe in alle kanale ’n geysering onstabiliteit getoon, gevolg deur flitsinggeïnduseerde (flashing induced) kook. Gelykfasige gedrag (vloei in kanale vertoon dieselfde gedrag) en ongelykfasige gedrag (vloei in kanale vertoon teenstrydige gedrag) tussen langsliggende kanale is met al die kragopwekkings en stelselwerkingsmodusse waargeneem. Vloei-omkering in die verhitte kanale van ’n natuurlike sirkulasiestelsel is bewys om selfs onder gelyke kragopwekkings te bestaan.

Tidal stream energy is a predictable source of renewable energy. Tidal stream turbines have been proposed as a way to extract useful energy from the tide. Many arrays of such devices will need to be installed to extract significant amounts of energy. The presence of an array of turbines within a tidal flow will impact the flowfield, as complex fluid interactions occur across multiple scales. This thesis is concerned with the behaviour of tidal turbines arrayed across channels. Experimental and analytical work is carried out to investigate array behaviour and to create new modelling tools to replicate this behaviour. Linear Momentum Actuator Disc Theory (LMADT) is employed to develop a new analytical model for a long row array of tidal turbines split into multiple smaller, co- linear row arrays. An argument of separation of scales is used to facilitate this model. It is found that increases in power extraction beyond that of a single continuous row array are possible. Experimental work is carried out on a row array of eight porous discs, simulating a short row array of tidal turbines. Disc porosity and spacing are varied to investigate thrust on the array, flow behaviour behind the array and an 'inferred' power removed from the flow. The results are compared to previously developed theoretical models. Good agreement is found with the trends of the analytical model, for example that there is a peak power coefficient which can be reached through appropriate selection of spacing and disc resistance. Differences from theory are found in the total thrust and power measurements, as well as in some aspects of the flow behaviour in the array wake. Reductions in thrust and power towards the ends of the array are also identified as 'end effects' which are not included in the analytical model. Based on these results a new semi-empirical model is proposed, using LMADT with experimental data closure. This model allows variation of the disc resistance across a row array. Values from the experimental work are used as inputs to the model, and the results compared to experimental measurements of flowspeed, thrust and power. Although agreement with experimental results is found in some areas, there are still some discrepancies between the analytical model and the experimental results. This indicates that there are additional factors that contribute to end effects on a short row array.