Simulating Everglades National Park hydrology and phosphorus transport under existing and future scenarios using numerical modeling

The Florida Everglades has a long history of anthropogenic changes which have impacted the quantity and quality of water entering the system. Since the construction of Tamiami Trail in the 1920's, overland flow to the Florida Everglades has decreased significantly, impacting ecosystems from the wetlands to the estuary. The MIKE Marsh Model of Everglades National Park (M3ENP) is a numerical model, which simulates Everglades National Park (ENP) hydrology using MIKE SHE/MIKE 11software. This model has been developed to determine the parameters that effect Everglades hydrology and understand the impact of specific flow changes on the hydrology of the system. As part of the effort to return flows to the historical levels, several changes to the existing water management infrastructure have been implemented or are in the design phase. Bridge construction scenarios were programed into the M3ENP model to review the effect of these structural changes and evaluate the potential impacts on water levels and hydroperiods in the receiving Northeast Shark Slough ecosystem. These scenarios have shown critical water level increases in an area which has been in decline due to low water levels. Results from this work may help guide future decisions for restoration designs. Excess phosphorus entering Everglades National Park in South Florida may promote the growth of more phosphorus-opportunistic species and alter the food chain from the bottom up. Two phosphorus transport methods were developed into the M3ENP hydrodynamic model to determine the factors affecting phosphorus transport and the impact of bridge construction on water quality. Results showed that while phosphorus concentrations in surface waters decreased overall, some areas within ENP interior may experience an increase in phosphorus loading which the addition of bridges to Tamiami Trail. Finally, phosphorus data and modeled water level data was used to evaluate the spectral response of Everglades vegetation to increasing phosphorus availability using Landsat imagery.

Performance assessment of a vertical axis turbine in a marine current flume tank and CFD modelling

As advances in numerical modelling techniques support the increased confidence in predictions from computer simulations, the need remains to have experimental verification built into the design process. This paper outlines the experimental investigation carried out on a shielded vertical axis turbine in a marine environment. The experiments consist of performance measurements and the use of particle image velocimetry on a small scale device in a marine current flume. The results demonstrate that the performance of the device can be modelled numerically; in particular, the results show that the numerical model used can correctly predict the increase in performance with Reynolds number.

Monitoring and modelling the morphodynamic evolution of a breached barrier beach system

Predicting the evolution of a coastal cell requires the identification of the key drivers of morphology. Soft coastlines are naturally dynamic but severe storm events and even human intervention can accelerate any changes that are occurring. However, when erosive events such as barrier breaching occur with no obvious contributory factors, a deeper understanding of the underlying coastal processes is required. Ideally conclusions on morphological drivers should be drawn from field data collection and remote sensing over a long period of time. Unfortunately, when the Rossbeigh barrier beach in Dingle Bay, County Kerry, began to erode rapidly in the early 2000’s, eventually leading to it breaching in 2008, no such baseline data existed. This thesis presents a study of the morphodynamic evolution of the Inner Dingle Bay coastal system. The study combines existing coastal zone analysis approaches with experimental field data collection techniques and a novel approach to long term morphodynamic modelling to predict the evolution of the barrier beach inlet system. A conceptual model describing the long term evolution of Inner Dingle Bay in 5 stages post breaching was developed. The dominant coastal processes driving the evolution of the coastal system were identified and quantified. A new methodology of long term process based numerical modelling approach to coastal evolution was developed. This method was used to predict over 20 years of coastal evolution in Inner Dingle Bay. On a broader context this thesis utilised several experimental coastal zone data collection and analysis methods such as ocean radar and grain size trend analysis. These were applied during the study and their suitability to a dynamic coastal system was assessed.

Modelling of a Stirling engine with parabolic dish for thermal to electric conversion of solar energy

Stirling engines with parabolic dish for thermal to electric conversion of solar energy is one of the most promising solutions of renewable energy technologies in order to reduce the dependency from fossil fuels in electricity generation. This paper addresses the modelling and simulation of a solar powered Stirling engine system with parabolic dish and electric generator aiming to determine its energy production and efficiency. The model includes the solar radiation concentration system, the heat transfer in the ther- mal receiver, the thermal cycle and the mechanical and electric energy conversion. The thermodynamic and energy transfer processes in the engine are modelled in detail, including all the main processes occur- ring in the compression, expansion and regenerator spaces. Starting from a particular configuration, an optimization of the concentration factor is also carried out and the results for both the transient and steady state regimes are presented. It was found that using a directly illuminated thermal receiver with- out cavity the engine efficiency is close to 23.8% corresponding to a global efficiency of 10.4%. The com- ponents to be optimized are identified in order to increase the global efficiency of the system and the trade-off between system complexity and efficiency is discussed.

Tidal effects on sea water intrusion in unconfined aquifers

A variable-density groundwater model is used to analyse the effects of tidal fluctuations on sea-water intrusion in an unconfined aquifer. It is shown that the tidal activity forces the sea-water to intrude further inland and it also creates a thicker interface than would occur without tidal effects. Moreover, the configuration of the interface is radically changed when the tidal fluctuations are included. This is because of the dramatic changes in the flow pattern and velocity of the groundwater near the shoreline. For aquifer depths much larger than tidal amplitudes, the tidal fluctuation does not have much effect on how far the sea-water intrudes into the aquifer; nevertheless, a significant change in the configuration of concentration contours because of the effect of tidal fluctuations is observed. This change is more noticeable at the top of the aquifer, near the water table, than at the bottom of the aquifer, and is caused by the infiltration of salt water into the top of the aquifer at higher tidal levels. A flatter beach slope, therefore, intensifies this phenomenon. The interface configurations do not change noticeably over the course of a tidal cycle. Neglecting tidal fluctuation effects results in an inaccurate evaluation of the water table elevation at the land end of the aquifer, although no distinguishable difference is seen between the water tables near the shoreline. Where the landward boundary condition is a constant head, the effects of tidal fluctuations on sea-water intrusion are more pronounced than for cases where the landward boundary condition is a specified flux. Also it is shown that the effects of tidal fluctuations are more significant for a sloping beach than for a vertical shoreline and the salt water intrudes further inland for the sloping case. (C) 1999 Elsevier Science B.V. All rights reserved.

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

Quantifying stresses and support requirements in the undercut and production level drifts of block and panel caving mines

A numerical modelling strategy has been developed in order to quantify the magnitude of induced stresses at the boundaries of production level and undercut level drifts for various in situ stress environments and undercut scenarios. The results of the stress modelling were in line with qualitative experiential guidelines and a limited number of induced stress measurements documented from caving sites. A number of stress charts were developed which quantify the maximum boundary stresses in drift roofs for varying in situ stress regimes, depths and undercut scenarios. This enabled many of the experiential guidelines to be quantified and bounded. A limited number of case histories of support and support performance in cave mine drifts were compared to support recommendations using the NGI classification system, The stress charts were used to estimate the Stress Reduction Factor for this system. The back-analyses suggested that the NGI classification system might be able to give preliminary estimates of support requirements in caving mines with modifications relating to rock bolt length and the support of production level intersections. (C) 2002 Elsevier Science Ltd. All rights reserved.

Experimental investigation of contaminant transport in coastal groundwater

Contaminant transport in coastal aquifers is of increasing interest since, with the development of coastal areas, contaminants from surface sources may enter coastal aquifers and pollute the groundwater flow. Coastal groundwater flow is complicated because of the presence of a freshwater-saltwater diffusion zone and the tidal variation of sea level at the seaward end. This paper investigates experimentally the behaviour of contaminant plumes with different densities in an unconfined coastal aquifer. Experiments were performed in a flow tank filled with glass beads as the porous medium. Results show that the dense contaminant has a more diffusive front than the less dense one in the seaward direction towards the coastline. The plume becomes more diffusive when it travels closer to the saltwater interface. On the contrary, the less dense contaminant presents a relatively sharp outline. It tends to migrate in the upper portion of the aquifer and exits in a concentrated manner over a small discharge area at the coastline, not further seaward under the sea. Non-dimensional parameters show that instabilities occur in our experiments for a density difference of 1.2% or larger between the contaminant and the ambient water. The experimental results provide guidance for field monitoring and numerical modelling. (C) 2002 Elsevier Science Ltd. All rights reserved.

Predicting cave initiation and propagation

The most widely used method for predicting the onset of continuous caving is Laubscher's caving chart. A detailed examination of this method was undertaken which concluded that it had limitations which may impact on results, particularly when dealing with stronger rock masses that are outside current experience. These limitations relate to inadequate guidelines for adjustment factors to rock mass rating (RMR), concerns about the position on the chart of critical case history data, undocumented changes to the method and an inadequate number of data points to be confident of stability boundaries. A review was undertaken on the application and reliability of a numerical method of assessing cavability. The review highlighted a number of issues, which at this stage, make numerical continuum methods problematic for predicting cavability. This is in particular reference to sensitivity to input parameters that are difficult to determine accurately and mesh dependency. An extended version of the Mathews method for open stope design was developed as an alternative method of predicting the onset of continuous caving. A number of caving case histories were collected and analyzed and a caving boundary delineated statistically on the Mathews stability graph. The definition of the caving boundary was aided by the existence of a large and wide-ranging stability database from non-caving mines. A caving rate model was extrapolated from the extended Mathews stability graph but could only be partially validated due to a lack of reliable data.

Modelagem da dispersão de odores em ambientes urbanos

O presente trabalho investigou o problema da modelagem da dispersão de compostos odorantes em presença de obstáculos (cúbicos e com forma complexa) sob condição de estabilidade atmosférica neutra. Foi empregada modelagem numérica baseada nas equações de transporte (CFD1) bem como em modelos algébricos baseados na pluma Gausseana (AERMOD2, CALPUFF3 e FPM4). Para a validação dos resultados dos modelos e a avaliação do seu desempenho foram empregados dados de experimentos em túnel de vento e em campo. A fim de incluir os efeitos da turbulência atmosférica na dispersão, dois diferentes modelos de sub-malha associados à Simulação das Grandes Escalas (LES5) foram investigados (Smagorinsky dinâmico e WALE6) e, para a inclusão dos efeitos de obstáculos na dispersão nos modelos Gausseanos, foi empregado o modelo PRIME7. O uso do PRIME também foi proposto para o FPM como uma inovação. De forma geral, os resultados indicam que o uso de CFD/LES é uma ferramenta útil para a investigação da dispersão e o impacto de compostos odorantes em presença de obstáculos e também para desenvolvimento dos modelos Gausseanos. Os resultados também indicam que o modelo FPM proposto, com a inclusão dos efeitos do obstáculo baseado no PRIME também é uma ferramenta muito útil em modelagem da dispersão de odores devido à sua simplicidade e fácil configuração quando comparado a modelos mais complexos como CFD e mesmo os modelos regulatórios AERMOD e CALPUFF. A grande vantagem do FPM é a possibilidade de estimar-se o fator de intermitência e a relação pico-média (P/M), parâmetros úteis para a avaliação do impacto de odores. Os resultados obtidos no presente trabalho indicam que a determinação dos parâmetros de dispersão para os segmentos de pluma, bem como os parâmetros de tempo longo nas proximidades da fonte e do obstáculo no modelo FPM pode ser melhorada e simulações CFD podem ser usadas como uma ferramenta de desenvolvimento para este propósito. Palavras chave: controle de odor, dispersão, fluidodinâmica computacional, modelagem matemática, modelagem gaussiana de pluma flutuante, simulação de grandes vórtices (LES).

Estudo de falha em placas laminadas usando métodos analíticos e numéricos

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia de Manutenção

Cartografia e avaliação geotécnica de blocos rochosos em plataformas costeiras (Gaia, NW Portugal)

Mestrado em Engenharia Geotécnica e Geoambiente

Study of AC losses in medium-sized high temperature superconducting coils

The study of AC losses in superconducting pancake coils is of utmost importance for the development of superconducting devices. Due to different technical difficulties this study is usually performed considering one of two approaches: considering superconducting coils of few turns and studying AC losses in a large frequency range vs. superconducting coils with a large number of turns but measuring AC losses only in low frequencies. In this work, a study of AC losses in 128 turn superconducting coils is performed, considering frequencies ranging from 50 Hz till 1152 Hz and currents ranging from zero till the critical current of the coils. Moreover, the study of AC losses considering two different simultaneous harmonic components is also performed and results are compared to the behaviour presented by the coils when operating in a single frequency regime. Different electrical methods are used to verify the total amount of AC losses in the coil and a simple calorimetric method is presented, in order to measure AC losses in a multi-harmonic context. Different analytical and numerical methods are implemented and/or used, to design the superconducting coils and to compute the total amount of AC losses in the superconducting system and a comparison is performed to verify the advantages and drawbacks of each method.

Análisis mecánico de albañilería arqueológica de adobe bajo cargas de compresión uniaxial: El caso de Huaca de la Luna

This article presents an experimental and numerical study for the mechanical characterization under uniaxial compressive loading of the adobe masonry of one of the most emblematic archaeological complex in Peru, 'Huaca de la Luna' (100-650AD). Compression tests of prisms were carried out with original material brought to the laboratory. For measuring local deformations in the tests, displacement transducers were used which were complemented by a digital image correlation system which allowed a better understanding of the failure mechanism. The tests were then numerically simulated by modelling the masonry as a continuum media. Several approaches were considered concerning the geometrical modelling, namely 2D and 3D simplified models, and 3D refined models based on a photogrammetric reconstruction. The results showed a good approximation between the numerical prediction and the experimental response in all cases. However, the 3D models with irregular geometries seem to reproduce better the cracking pattern observed in the tests.

A FEM-based model to study the behaviour of corroded RC beams shear repaired by NSM CFRP rods technique

This paper presents the main features of finite element FE numerical model developed using the computer code FEMIX to predict the near-surface mounted NSM carbon-fiber-reinforced polymer CFRP rods shear repair contribution to corroded reinforced concrete RC beams. In the RC beams shear repaired with NSM technique, the Carbon Fibre Reinforced Polymer (CFRP) rods are placed inside pre-cut grooves onto the concrete cover of the RC beam’s lateral faces and are bonded to the concrete with high epoxy adhesive. Experimental and 3D numerical modelling results are presented in this paper in terms of load-deflection curves, and failure modes for 4 short corroded beams: two corroded beams (A1CL3-B and A1CL3-SB) and two control beams (A1T-B and A1T-SB), the beams noted with B were let repaired in bending only with NSM CFRP rods while the ones noted with SB were repaired in both bending and shear with NSM technique. The corrosion of the tensile steel bars and its effect on the shear capacity of the RC beams was discussed. Results showed that the FE model was able to capture the main aspects of the experimental load-deflection curves of the RC beams, moreover it has presented the experimental failure modes and FE numerical modelling crack patterns and both gave similar results for non-shear repaired beams which failed in diagonal tension mode of failure and for shear-repaired beams which failed due to large flexural crack at the middle of the beams along with the concrete crushing, three dimensional crack patterns were produced for shear-repaired beams in order to investigate the splitting cracks occurred at the middle of the beams and near the support.