The topographical variations of the seas around Taiwan are great, which make the tides complicated. Taiwan is located in the juncture of the tropical and subtropical area. Geographically, it is located within the region of northwestern Pacific typhoon path. These seasonal and geographical situations causing Taiwan frequently threaten by typhoons during summer and autumn. In addition to natural disasters, the coastal area is over developed for the last few decades, which destroys the balance between nature and man. Storms and floods constantly threaten the lowland areas along the coast. An accurate and efficient storm surge model can be used to predict tides and storm surges. The model can be calibrated and verified with the field observations. Data measured by instruments at the tidal station constituting daily tidal variations and storm surge influences during typhoons. The model can offer both predictions to the management institutions and to the general public as pre-warning system and thus taking disaster-prevention measures. This study implements the numerical model, developed by Yu (1993) and Yu et al. (1994) to calculate the hydrodynamic in the seas around Taiwan. The main purpose of this study is to make a calibration and sensitivity analysis of the model parameters. Tidal gauge data around Taiwan coastal stations collected from June to October 2005 are used for the analysis and the comparison between the modeled data and the observations. Two steps have been taken for the model calibration and sensitivity analysis. First step is to calibrate the model for accurate prediction of the astronomical tide, and then the compound tide with meteorological influences. For the calibration of the astronomical tides, sensitivity analysis has been carried out by adjusting the horizontal diffusion coefficient and the bottom friction coefficients used in the model. The sensitivity of the time-step size used in the model and model grids fitted to coastlines are also checked. A depth dependent Chézy numbers are used in the model to describe bottom friction. The model has a better result when the Chézy value varied within 65 to 85. Modifying grids fitted to the coastline has improved the model results significantly. By improving the dynamic phenomenon brought about by the land features, the model calculation fits the real tidal phenomenon better. The analysis has shown that the model is less sensitive to the horizontal diffusion coefficient. Data from 22 tidal stations around Taiwan have been used for the comparisons. The maximum RMSE (root-mean-square error) is about 10 cm at WAi-Pu, whereas the minimum RMSE is about 1 cm for the stations along eastern coast. The calibration of the compound tide is divided into three cases. The first case is to calibrate the forecasted wind field. This has been done by comparing the forecasted wind field from the Central Weather Bureau with the satellite data obtained from QuikSCAT¡XLevel 3. The satellite wind speed has been applied to adjust the forecasted wind speed. The adjusted forecast wind field has shown improvement to the model predictions in the tidal stations south of Taichung, slightly improved in the eastern coast. The second case is tuning the drag coefficient on sea surface used by the hydrodynamic model. Several empirical formulas to describe the sea surface drag have been tested. The model result has shown little influence using various drag formulations. The third case is to single the influences by the meteo-inputs, i.e. the wind field and the atmospheric pressure. The tidal level is more sensitive to the variation of the atmospheric pressure through out the tests carried out during typhoon periods. The model simulation for 2006 using the best selected parameters has shown that the model is consisted with good stability and accuracy for both stormy and calm weather conditions.

The complex interaction between flow, vegetation and sediment drives the never settled changes of riverine system. Vegetation intercepts rainfall, adds resistance to surface flow, and facilitates infiltration. The magnitude and timing of flood flow are closely related to the watershed vegetation coverage. In the meantime, flood flow can transport a large amount of sediment resulting in bank erosion, channel degradation, and channel pattern change. As climate changes, future flood frequency will change with more intense rainfalls. However, the quantitative simulation of flood flow in vegetated channel and the influence of climate change on flood frequency, especially for the arid and semi-arid Southwest, remain challenges to engineers and scientists. Therefore, this research consists of two main parts: simulate unsteady flow and sediment transport in vegetated channel network, and quantify the impacts of climate change on flood frequency. A one-dimensional model for simulating flood routing and sediment transport over mobile alluvium in a vegetated channel network was developed. The modified St. Venant equations together with the governing equations for suspended sediment and bed load transport were solved simultaneously to obtain flow properties and sediment transport rate. The Godunov-type finite volume method is employed to discretize the governing equations. Then, the Exner equation was solved for bed elevation change. Since sediment transport is non-equilibrium when bed is degrading or aggrading, a recovery coefficient for suspended sediment and an adaptation length for bed load transport were used to quantify the differences between equilibrium and non-equilibrium sediment transport rate. The influence of vegetation on floodplain and main channel was accounted for by adjusting resistance terms in the momentum equations for flow field. A procedure to separate the grain resistance from the total resistance was proposed and implemented to calculate sediment transport rate. The model was tested by a flume experiment case and an unprecedented flood event occurred in the Santa Cruz River, Tucson, Arizona, in July 2006. Simulated results of flow discharge and bed elevation changes showed satisfactory agreements with the measurements. The impacts of vegetation density on sediment transport and significance of non-equilibrium sediment transport model were accounted for by the model. The two-dimensional surface flow model, called CHRE2D, was improved by considering the vegetation influence and then applied to Santa Cruz River Watershed (SCRW) in the Southern Arizona. The parameters in the CHRE2D model were calibrated by using the rainfall event in July 15th, 1999. Hourly precipitation data from a Regional Climate Model (RCM) called Weather Research and Forecasting model (WRF), for three periods, 1990-2000, 2031-2040 and 2071-2079, were used to quantify the impact of climate change on the magnitude and frequency of flood for the Santa Cruz River Watershed (SCRW) in the Southern Arizona. Precipitation outputs from RCM-WRF model were bias-corrected using observed gridded precipitation data for three periods before directly used in the watershed model. The watershed model was calibrated using the rainfall event in July 15th, 1999. The calibrated watershed model was applied to SCRW to simulate surface flow routing for the selected three periods. Simulated annual and daily maximum discharges are analyzed to obtain future flood frequency curves. Results indicate that flood discharges for different return periods are increased: the discharges of 100-year and 200-year return period are increased by 3,000 and 5,000 cfs, respectively.

<p><font face="TimesNewRoman"> <p align="left">Research findings in current study provide a new insight into the fractured rock aquifers in the TMG area. Some of the results will have wide implications on the groundwater management and forms a solid basis the further study of the TMG aquifers.</p> </font></p>

A fundamental stage in rocky-shore seaweed life history is the recruitment process involving external fertilisation and then settlement of the propagules on a suitable substrate. The ultimate step in this settlement stage is the crossing of the viscous sub-layer and attachment to the substrate. Given the extreme conditions met in the intertidal zone, propagules can be dislodged at any time before they secure a strong enough anchoring. Flow conditions and propagule properties are key to this process. The settlement process under turbulent conditions was recreated within a stirred benthic chamber for five different species. Whereas propagule properties (size, density) vary with species, and propagules are adapted to the different conditions in the intertidal, they exhibit the same settlement behaviour. They nevertheless exhibit different settling velocities and settlement thresholds. Several methods of characterisation of the tank flow from particle tracking velocimetry and acoustic Doppler velocimetry data are reviewed, as well as an analytical model. Turbulent settling was found to be independent of the well-mixed tank bulk flow and to depend only on the boundary-layer mechanics. A model of settlement threshold is presented from which propagule mucilage adhesiveness estimates are derived, leading to good correlations between adult plant exposure and the stickiness of its propagules and to the conclusion that settlement can only occur in calm conditions. To extend the work, computational fluid dynamic techniques are developed by extending the Gerris Flow Solver. A 2-D approach to tank modelling and a pilot study of expansion to 3-D is described. This extends the perspective given by the experiments, notably through output of the hydrodynamical forces experienced by the propagules. Finally, in the view of realizing direct numerical simulations of propagule behaviour in the viscous sub-layer, a new and unique 2-D/3-D fully conservative solid/fluid interaction model is developed and tested with success.

This thesis presents a two-phase liquid-solid numerical model using Smoothed Particle Hydrodynamics (SPH). The scheme is developed for multi-phase flows in industrial tanks containing sediment used in the nuclear industry for decommissioning. These two-phase liquid-sediments flows feature a changing interfacial profile, large deformations and fragmentation of the interface with internal jets generating resuspension of the solid phase. SPH is a meshless Lagrangian discretization scheme whose major advantage is the absence of a mesh making the method ideal for interfacial and highly non-linear flows with fragmentation and resuspension. Emphasis has been given to the yield profile and rheological characteristics of the sediment solid phase using a yielding, shear and suspension layer which is needed to predict accurately the erosion phenomena. The numerical SPH scheme is based on the explicit treatment of both phases using Newtonian and non-Newtonian Bingham-type constitutive models. This is supplemented by a yield criterion to predict the onset of yielding of the sediment surface and a suspension model at low volumetric concentrations of sediment solid. The multi-phase model has been compared with experimental and 2-D reference numerical models for scour following a dry-bed dam break yielding satisfactory results and improvements over well-known SPH multi-phase models. A 3-D case using more than 4 million particles, that is to the author’s best knowledge one of the largest liquid-sediment SPH simulations, is presented for the first time. The numerical model is accelerated with the use of Graphic Processing Units (GPUs), with massively parallel capabilities. With the adoption of a multi-phase model the computational requirements increase due to extra arithmetic operations required to resolve both phases and the additional memory requirements for storing a second phase in the device memory. The open source weakly compressible SPH solver DualSPHysics was chosen as the platform for both CPU and GPU implementations. The implementation and optimisation of the multi-phase GPU code achieved a speed up of over 50 compared to a single thread serial code. Prior to this thesis, large resolution liquid-solid simulations were prohibitive and 3-D simulations with millions of particles were unfeasible unless variable particle resolution was employed. Finally, the thesis addresses the challenging problem of enforcing wall boundary conditions in SPH with a novel extension of an existing Modified Virtual Boundary Particle (MVBP) technique. In contrast to the MVBP method, the extended MVBP (eMVBP) boundary condition guarantees that arbitrarily complex domains can be readily discretized ensuring approximate zeroth and first order consistency for all particles whose smoothing kernel support overlaps the boundary. The 2-D eMVBP method has also been extended to 3-D using boundary surfaces discretized into sets of triangular planes to represent the solid wall. Boundary particles are then obtained by translating a full uniform stencil according to the fluid particle position and applying an efficient ray casting algorithm to select particles inside the fluid domain. No special treatment for corners and low computational cost make the method ideal for GPU parallelization. The models are validated for a number of 2-D and 3-D cases, where significantly improved behaviour is obtained in comparison with the conventional boundary techniques. Finally the capability of the numerical scheme to simulate a dam break simulation is also shown in 2-D and 3-D.

L'utilisation extensive des matériaux composites dans l'industrie aéronautique nécessite une compréhension étendue de leur comportement mécanique et mécanismes d'endommagement. La prise en compte des effets environnementaux, température et vieillissement humide, dans la prévision du comportement à l'impact, des endommagements et des mécanismes associés reste un sujet encore peu traité actuellement. Dans ce travail différentes conditions de vieillissement humide (dénommées ci-après état neuf, 65% VH, et 90%VH) ont été appliquées au matériau de l'étude. La caractérisation des propriétés du pli et de l'interface ont été effectuée à température ambiante et 80°e. Une attention particulière a été portée à l'identification des propriétés en cisaillement plan. Un modèle d'endommagement progressif basé sur la mécanique de l'endommagement dans les milieux continus est proposé afin de décrire au mieux le comportement du pli élémentaire en prenant en compte les comportements fortement non-linéaires identifiés pour le cisaillement plan de ce matériau. Ce modèle matériau de comportement du pli est combiné avec des lois cohésives d'interface inter-laminaire utilisés dans des simulations numériques d'impact à faible vitesse afin de prédire la réponse mécanique à l'impact et les dommages. Le modèle de matériau a été implémenté dans le code de calcul explicite LS-Dyna. Des essais d'impact à différentes énergies ont été effectués afin de valider par comparaison essais-calcul l'approche utilisée. Une bonne adéquation est observée entre les résultats numériques et expérimentaux du matériau non-vieilli et testé à température ambiante pour ce qui concerne l'historique de la réponse globale et les dommages internes (rupture de fibre, fissures de matrice et délaminages). / The extensive application of composite laminates in aeronautical industries demands a thorough understanding of their mechanical behaviours and damage. The prediction of impact behaviour, impact damage and the associated mechanisms under complex environmental effects remain still an open issue. The effects of various humid ageing conditions (No ageing, 65% relative humidity and 90% relative humidity ageing) and environ mental temperatures (room temperature and 80°C) on the intra-Iaminar and inter-Iaminar properties are characterized. A particular emphasis is laid on the study of different non linear behaviours of in plane shear. A progressive damage model based on continuum damage mechanics is proposed to describe intra laminar mechanical behaviour by taking into account the identified in-plane shear nonlinear phenomenon. This model is then applied in numerical simulation on low velocity impact combined with cohesive law in interface to predict the impact response and associated damage (ie. Fibre failure, matrix cracks and delaminations). The damage model is implemented in LS Dyna as user defined material laws.lmpact tests on laminate coupons with different stratifications are performed to validate the numerical simulation. Good agreement is achieved especially for nonaged material under room temperature between numerical and experimental results in terms of contact force history, platedeflection and impact duration. By means of non-destructive inspection, the impact damage is coherent with predicted results concerning fiber damage location, matrix cracks especially the through thickness delamination.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / In this thesis, we presents a development and validation of a 3D numerical model for the advection-diffusion equation. Models of this kind has been developed for scientific investigations and to support atmospheric emissions control and environmental policy decisions. To develop this model, we used the computational implementation of an explicit numerical scheme for the discretization of the envolved equations. During this procedures, exaustive tests were performed to ensure that the used implementations agrees to the stability, consistence and convergence criterias. As a way for minimizing one of the main deficiences found in almost the major atmospheric dispersion models, i.e. imprecisions in the meteorological input data for initializing this models, we used a realistic atmospheric flow field generated by mesoscale circulation model. As the mesoscale model gives information at scale larger than the necessary for description of a plume trajectory, a weighted linear average proper interpolation was developed for intermediate these distances. Our model considers the assumption that atmospheric turbulence is not isotropic, where diffusion coefficients are variables in time and space and are different for lateral and vertical directions. In our model we estimate this coefficients by the atmospheric boudary layer parameterizations proposed byMoraes (2000). For validation of themodel, we used experimental datasets from field experiment carried near a thermoelectric power plant presidente Médici, in the city of Candiota/RS. These datasets contains surface SO2 concentrations, surface wind velocity measured in meteorological towers as well as turbulence data measured in micrometeorological towers. The results of the validation indicates that the model works well, at least for the source and the terrain were it is located. i.e. continuous emission and homogeneous topography. / Nesta tese, apresentamos o processo de desenvolvimento e validação de um modelo numérico 3D para a equação de difusão-advecção. Modelos como este têm sido desenvolvidos objetivando a investigação científica para dar suporte ao controle de emissões de poluentes atmosféricos e à tomada de decisões no desenvolvimento de políticas ambientais. Para a elaboração deste modelo utilizamos a implementação de um esquema numérico explícito na discretização das equações envolvidas. Durante este processo, exaustivos testes foram realizados para que a implementação garantisse estabilidade, consistência e convergência. Como forma de minimizar uma das principais deficiências encontradas na maioria dos modelos de dispersão atmosféricos (MDA), qual seja, a imprecisão da entrada de dados meteorológicos para a inicialização destes modelos, utilizamos um campo de vento realístico gerado por um modelo numérico de circulação de mesoescala. No entanto, como este modelo de mesoescala fornece informações com uma escala maior do que a necessária para a descrição da trajetória de uma pluma, foi preciso desenvolver um método de interpolação apropriado para intermediar estas distâncias. Nosso modelo contempla, ainda, a hipótese da turbulência atmosférica não ser isotrópica, onde os coeficientes de difusão turbulenta são variáveis no tempo e espaço e, ainda, diferentes nas direções lateral e vertical. Em nosso modelo, estes coeficientes foram estimados utilizando a parametrização de camada limite proposta por Moraes (2000). Para a validação do modelo, utilizamos dados experimentais obtidos em um experimento realizado próximo à usina termelétrica Presidente Médici, no município de Candiota/RS. Estes dados compreendem medidas de concentração superficial de SO2, velocidade do vento em superfície medidos em torres meteorológicas, bem como dados de turbulência medido em torre micrometeorológica. Os resultados da validação mostraram que o modelo funciona bem, ao menos para o tipo de fonte e topografia onde está localizada, ou seja, emissão contínua e topografia homogênea.

A presente tese versa sobre ampla investigação em torno da distribuição das ações verticais entre paredes de alvenaria estrutural dentro de um edifício. O edifício foi caracterizado como pertencente ao sistema complexo de paredes, o qual é de uso mais comum no Brasil. Desenvolveu-se análise numérica no edifício, com objetivo de predizer o comportamento na análise experimental, assim como na comparação entre os próprios modelos numéricos. Serviu de modelo experimental um edifício em escala real, o qual foi monitorado, durante sua execução e depois dela, com ensaio de prova de carga. Realizou-se monitoramento mediante uso de unidades especiais de alvenaria instrumentadas para se converterem em células de carga, nas quais as tensões em cada ponto de análise puderam ser determinadas. Tabelas e gráficos foram montados para tornarem possível a verificação das diferenças entre todos os modelos analisados. Como principais conclusões uma se mostrou de grande importância e foi esta: as paredes sofrem efeito da flexão e a causa são as rotações nas lajes apoiadas por elas. / The present thesis contemplated a wide investigation about the vertical action distribution among structural masonry walls within a building. The building was characterized as related to the complex system of walls, which practice is the most customary in Brazil. It was developed a numerical analysis in the building aiming to predict the structural behavior on the experimental analysis, as well as a comparison among numerical models themselves. The experimental model consisted on a real scale building, which was monitored during its execution and after that, with a live load test. The monitoring was done with the usage of special masonry units, which were instrumented to convert itself in a load cell, on which the stress in each point of analysis could be determined. Tables and graphics were mounted to make feasible to verify the differences among all models analyzed. As main conclusions observed, one showed to be of great importance, which was the effect of bending on the walls caused by the rotations on the slabs supported by them.

A utilização de concreto projetado como suporte de túneis é uma prática amplamente difundida no mundo inteiro. Este tipo material possui a característica de começar a agir estruturalmente desde pequenas idades. Apesar disso, os correntes processos de dimensionamento de suportes negligenciam o desenvolvimento de suas propriedades com o tempo, em acoplamento aos efeitos tridimensionais da região onde se localiza a frente de escavação. O presente trabalho tem a finalidade de relatar os procedimentos utilizados na análise da influência de alguns parâmetros da interação maciço - suporte, sobre os esforços solicitantes e os deslocamentos radiais finais do suporte de um túnel, tanto para o caso de concreto projetado, com suas propriedades dependentes do tempo, quanto para materiais com propriedades constantes. São elaboradas soluções adimensionais para o problema da quantificação de esforços solicitantes no suporte e de deslocamentos radiais na interface entre maciço e suporte. Além disso, é proposta uma forma de se determinar, através destas soluções adimensionais, coeficientes de alívio de tensões que auxiliam em simulações bidimensionais de escavações subterrâneas. / The utilization of shotcrete as tunnel support is a widely diffused practice in the whole world. This kind of material has the feature of beginning to act structurally since early ages. Nevertheless, the current processes of support design neglect the development of its properties after some time in connection to the 3D effects of the region where the face of the tunnel is located. The present work relates the procedures adopted in analyzing the influence of some ground-support interaction parameters on the support internal forces and interface radial displacements of a tunnel, not only when shotcrete is used, with its time-dependent properties, but for materials with constant properties as well. Dimensionless solutions are developed for the support thrust and radial displacement quantification problem. Moreover, through this dimensionless solution, a way of quantifying stress relief factors, which are intermediate steps in 2D simulations of underground excavations, is proposed.

Orientador: Luiz Felipe Mendes de Moura / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-12T21:08:34Z (GMT). No. of bitstreams: 1 Brandao_MatheusMarcelino_M.pdf: 452664 bytes, checksum: 124f6ab109716c77cffe6f192c006b30 (MD5) Previous issue date: 2005 / Resumo: No presente trabalho se realiza o estudo e desenvolvimento de um modelo numérico computacional com a finalidade da simulação de sistemas de refrigeração por compressão de vapor. Em um primeiro momento, foram analisados e discutidos os principais tipos de modelamento para cada componente do sistema de refrigeração (compressor, condensador, válvula de expansão e evaporador) encontrados na literatura. Após esta primeira análise, foram escolhidos os modelos utilizados para o desenvolvimento do programa computacional. O programa foi validado a partir da comparação dos resultados obtidos numericamente com os obtidos em ensaios experimentais. Após a validação, foram feitos dois estudos de casos, para a verificação da influência dos parâmetros sobre os sistemas simulados e otimização do desempenho de tais sistemas / Abstract: This work performs the study and development of a numerical computational model in order to simulate refrigeration systems by steam compression. At first, the main types of modelation for each component of the refrigeration system (compressor, condenser, expansion valve and evaporator) found in literature are analysed and discussed. After this first analysis, the models for the development of a computational program are chosen. The program is validated by means of the comparison of the results obtained numerically with those acquired in experimental tests. Just after validation, two studies of cases are made to verify the influence of the parameters on the simulated systems and to optimize the performance of those systems / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica