Bahia de la Ascension (BA) is a pristine, shallow, karst bay located in the Mexican Caribbean, a region experiencing rapid population growth stimulated by intense tourism development. The overall objective of this study was to address the natural hydrographic variability of this inherently vulnerable ecosystem and assess its influence on a key habitat, the seagrass. The chapters follow the three-branched nature of the study which tackled the connected ecosystem issues of coastal hydrology, physical dynamics of flow and circulation, and the ecological dynamics of the seagrass species Thalassia testudinum in BA. Freshwater input to BA is primarily by submerged groundwater discharges and surface runoff; both sources are derived from fissures in the aquifer but feature distinct water quality due to the interaction with adjacent wetlands. Hurricanes explain 36 percent of the interannual precipitation variability in the region. The water balance indicates a persistent net outflow from BA to the adjacent shelf, suggesting an intense exchange across inlets. Both diurnal and semidiurnal tidal frequencies are attenuated in the inner bay, where a meteorologically-induced subtidal water level increase may occur during four-day southeasterly winds. A clear SW-NE salinity gradient was established during dry and rainy seasons, with a strong tidally-driven marine influence throughout the central basin, and a perennial mesohaline ambient in the southwestern-most bay, where hydrodynamics are primarily controlled by wind stress. Thalassia testudinum is the dominant seagrass species in BA, occupying ~90 percent of the substrate, including the freshwater-influenced inner bay. High nutrient inputs, including phosphorus which might have limiting effects in karst environments, along with the wind-driven circulation controlling water residence times are associated with the successful development of T. testudinum (up to 1,461.23 g DW m-2) within the SW bay. Farthest into the central basin, Thalassia consistently exhibited an inverse correlation between abundance and density of shoots. This pattern was enhanced under exceptional precipitation and inputs of denuded organic matter resulting from hurricanes making landfall on this region. The relationship between nutrient distribution and the above/belowground ratio suggested that Thalassia growing in BA favors the development of the aerial component as nutrients availability increases. This study provides a basic understanding of the most important processes molding the patterns of variability exhibited by T. testudinum in Bahia de la Ascension. The salinity gradient and external nutrient supply, along with the hydrodynamic component, define the spatial scale at which the connectivity between the adjacent wetland, the bay, and the shelf may occur.

This thesis documents the studies conducted in deriving a mathematical model representing the dynamics of a maneuvering ship to be implemented as part of an interactive real-time simulation system, as well as the details and results of the implementation process itself. Different effects on the dynamics of ship motions are discussed separately, meaning that the effects are considered to be applied to the system one at a time and they are included in the model simply by the principle of superposition. The model is intended to include the hydrodynamic interactions between the ship hull and the ocean via added mass (added inertia), damping and restoring force concepts. In addition to these effects, which are derived considering no incident waves are present on the ocean, the environmental disturbances, such as wind, wave and ocean current are also taken into account for proposing a mathematical model governing the dynamics of the ship. Since the ultimate product of this thesis work is a running computer code that can be integrated into an available simulation software, the algorithm development and code implementation processes are also covered. Improvements made on the implementation to achieve &ldquo / better&rdquo / real-time performance are evaluated comparatively in reference to original runs conducted before the application of improvement under consideration. A new method to the computation of the wave model that allows faster calculation in real-time is presented. A modular programming approach is followed in the overall algorithm development process in order to make the integration of new program components into the software, such as a new hull or propulsion model or a different integrator type possible, easily and quickly.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Esta tese tem por objetivo propor uma estratégia de obtenção automática de parâmetros hidrodinâmicos e de transporte através da solução de problemas inversos. A obtenção dos parâmetros de um modelo físico representa um dos principais problemas em sua calibração, e isso se deve em grande parte à dificuldade na medição em campo desses parâmetros. Em particular na modelagem de rios e estuários, a altura da rugosidade e o coeficiente de difusão turbulenta representam dois dos parâmetros com maior dificuldade de medição. Nesta tese é apresentada uma técnica automatizada de estimação desses parâmetros através deum problema inverso aplicado a um modelo do estuário do rio Macaé, localizado no norte do Rio de Janeiro. Para este estudo foi utilizada a plataforma MOHID, desenvolvida na Universidade Técnica de Lisboa, e que tem tido ampla aplicação na simulação de corpos hídricos. Foi realizada uma análise de sensibilidade das respostas do modelo com relação aos parâmetros de interesse. Verificou-se que a salinidade é uma variável sensível a ambos parâmetros. O problema inverso foi então resolvido utilizando vários métodos de otimização através do acoplamento da plataforma MOHID a códigos de otimização implementados em Fortran. O acoplamento foi realizado de forma a não alterar o código fonte do MOHID, possibilitando a utilização da ferramenta computacional aqui desenvolvida em qualquer versão dessa plataforma, bem como a sua alteração para o uso com outros simuladores. Os testes realizados confirmam a eficiência da técnica e apontam as melhores abordagens para uma rápida e precisa estimação dos parâmetros. / This thesis presents a strategy for automatically obtaining hydrodynamic and transport parameters by means of the solution of inverse problems. Obtaining the parameters of a physical model represents a major problem in its calibration, and this is largely due to the difficulty associated to the field measurements of these parameters. In particular by modeling rivers and estuaries, the roughness height and the turbulent diffusion coefficient represent two of the most difficult parameters to measure or estimate. Here, an automated technique for estimation of these parameters through an inverse problem is applied to a model of the Macaé estuary, located in northern Rio de Janeiro. For this investigation, hydrodynamic and transport models were built in the MOHID platform, developed in the Technical University of Lisbon, which has had wide application in simulation of water bodies. A sensitivity analysis was performed regarding the model responses with respect to the parameters of interest. It was verified that salinity is a sensitive variable for both parameters. The inverse problem was then solved using various optimization methods by coupling the MOHID platform to optimization codes implemented in Fortran. The coupling was carried in a way to not change the MOHID source code, allowing the use of the computational tool here developed in any version of this platform, as well as its modification for use with other simulators. The tests confirm the efficiency of the technique and suggest the best approaches for a fast and accurate estimation of the parameters.

Thailand has been often affected by severe flood events over the past century. The 2011’s Thailand Flood Catastrophe was the costliest in country’s history, and it was ranked to be the second most damaging natural hazard in the world in terms of economic losses. The Chao Phraya River Basin was noted to be the most vulnerable area prone to flooding in Thailand. The dynamics of flood risk in the river basin have changed drastically over the past fifty years. In particular, flood exposure increased due to rapid urbanization and population growth. Since 2012, integrated flood risk management has been addressed to be the major framework of water-related disasters with the goal of losses and damage reductions. However, there is currently little research in Thailand on how to quantify flood risks and mitigate flood inundation damage on the relation between the occurrence of flood events and their consequential socio-economic implications. In this study, a tradition method in flood risk assessment is implemented by integrating 2D hydrodynamic modelling and the assessment of socio-economic impact of floods into the Chao Phraya River Basin. More specifically, the fully 2D version of the LISFLOOD-FP model code was used to model flood inundation processes. The output of the model was then used to map inundation depth and assess the levels of physical/environmental risk associated to flood hazards on multiple receptors/elements at risk. The European Flood Directive and the KULTURisk methodology were applied to quantify flood risks in monetary terms for residential, industrial, and agricultural sectors. The 2011 flood event was used for model calibration, while a hypothetical flood event with a return period of 100 years was simulated to identify the potential flood losses. Depth-damage functions comprising of JRC-ASIA, the Flemish, and JICA models were used to estimate potential damage for residential and industrial structures. The results showed that LISFLOOD-FP could satisfactorily reproduce the flood inundation extent obtained from satellite imagery in 2011. The model performance (Critical Success Index or F1) was of 56%, with a Bias of 112%. The latter meant the total inundated area was 12% larger than flood extent’s observation. Moreover, the model could simulate flood levels with overall Root Mean Square Error (RMSE) and Mean Absolute Error (MAE) of 2.03 m a.s.l. and 1.78 m a.s.l., respectively. For the estimation of flood damage and losses, the Flemish model showed the strongest agreement with the reported flood damage in the residential sector, while JICA-ASIA model underestimated flood damage for industrial sector by just 1%. The KULTURisk methodology also well-estimated crop losses in the 2011 event which an overestimation about 21% from the reported value. Apart from that, fully 2D numerical method could not perfectly represent 1-in-100 year flood inundation due to non-consideration of important features such as the precise river channel topography, hydraulic infrastructures, and flood protection schemes in the river basin. Lack of such features results in an overestimation of flood damage and losses for 1-in-100 year flood comparing to the national flood hazard map and damage assessment which are simulated and estimated by JICA’s study. Such features can be better handled by using a coupled 1D/2D numerical method in order to simulate flood inundation extent more realistically and estimate flood losses. This could help the Thai government to better prepare a budget for flood risk prevention. In addition, even if the Flemish model indicates a good representation of relative flood damage to housing structures, the government should establish depth-damage curves specific for Thailand.

Devido à proximidade de grandes cidades, ecossistemas aquáticos são alvos de despejos industriais, esgotos domésticos e resíduos sólidos. Modelos matemáticos são largamente utilizados para entender os padrões de fluxo e o transporte de substâncias nestes ecossistemas. Neste trabalho, desenvolveu-se um modelo hidrodinâmico e de transporte de poluentes bidimensional para corpos d’água rasos que utiliza grades não estruturadas. O modelo foi desenvolvido em dois módulos: hidrodinâmico e de transporte de massa. O módulo hidrodinâmico utiliza as equações completas de Saint Venant para quantificação da circulação em corpos d’águas. No módulo hidrodinâmico aqui desenvolvido é feito um aperfeiçoamento de um esquema semi-implícito de diferenças finitas utilizado no modelo IPH-ECO para um esquema semi-implícito de volumes finitos. O módulo de transporte de massa, por sua vez, representa o transporte de um escalar conservativo sob influência dos processos físicos do meio. Dois esquemas foram empregados e comparados para a solução numérica da equação de transporte: um esquema de diferenças centrais e um esquema com limitador de fluxo. O modelo foi aplicado em um estudo de caso simplificado formado por uma lagoa circular e no rio Guaíba. O módulo hidrodinâmico foi avaliado através de testes de conservação de volume feitos comparando com os resultados encontrados nno modelo IPHECO na aplicação realizada na lagoa circular e na solução analítica no rio Guaíba. O módulo hidrodinâmico do modelo de grades não estruturadas (IPH-UnTRIM2D) mostrou-se mais conservativo que o modelo IPH-ECO em todos os cenários simulados na lagoa circular. A solução hidrodinâmica e de transporte do modelo de grades não estruturadas no estudo de caso idealizado no rio Guaíba representou meandros, sulcos, estreitamentos e alargamentos de canais com comprometimento com uma situação real. Quanto à simulação de transporte de escalares e poluentes, foi comprovada a eficiência do esquema com limitador de fluxo sobre os esquemas de diferenças centrais. Os testes de conservação de massa para diferentes cenários no rio Guaíba indicaram erros pouco significativos em ambos os esquemas quando comparados com a ordem de grandeza do ecossistema. / Aquatics ecossistems around large towns are always exposed waste domestic and industrial and solid residues. Mathematical models are widely used to understand flow patterns and transport of mass in these ecosystems. In this work a model was developed on unstructured two-dimensional complex dynamic model for shallow lakes. The model was divided in two modules: (a) hydrodynamic module, describing quantitative flows and water level; (b) a mass transport module. The hydrodynamic module solves Saint Venant equations to quantify water circulation in a water body. Futhermore, the hydrodynamic module was improved to semiimplicit finite volume scheme. The mass transport module represents transport of a conservative scalar under the physical processes influence. Two schemes were used and compared to the transport equation numerical solution : a central difference scheme and a flux limiter scheme. The model was applied in the circular lake and river Guaíba. The hydrodynamic module was validated by evaluating tests of volume conservation made comparing with results of the model IPH-ECO and analytical solution. It proved to be more conservative than the model IPH-ECO in all scenarios simulated in the circular lake. The hydrodynamic and transport solution in river Guaíba study case represented bends, grooves, and channels widening and narrowing. In the transport simulation, the flux limiter scheme was more efficient than central differences scheme. Futhermore, mass conservation tests for different scenarios in Guaíba indicated minor errors in both schemes when compared with ecosystem magnitude.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Esta tese tem por objetivo propor uma estratégia de obtenção automática de parâmetros hidrodinâmicos e de transporte através da solução de problemas inversos. A obtenção dos parâmetros de um modelo físico representa um dos principais problemas em sua calibração, e isso se deve em grande parte à dificuldade na medição em campo desses parâmetros. Em particular na modelagem de rios e estuários, a altura da rugosidade e o coeficiente de difusão turbulenta representam dois dos parâmetros com maior dificuldade de medição. Nesta tese é apresentada uma técnica automatizada de estimação desses parâmetros através deum problema inverso aplicado a um modelo do estuário do rio Macaé, localizado no norte do Rio de Janeiro. Para este estudo foi utilizada a plataforma MOHID, desenvolvida na Universidade Técnica de Lisboa, e que tem tido ampla aplicação na simulação de corpos hídricos. Foi realizada uma análise de sensibilidade das respostas do modelo com relação aos parâmetros de interesse. Verificou-se que a salinidade é uma variável sensível a ambos parâmetros. O problema inverso foi então resolvido utilizando vários métodos de otimização através do acoplamento da plataforma MOHID a códigos de otimização implementados em Fortran. O acoplamento foi realizado de forma a não alterar o código fonte do MOHID, possibilitando a utilização da ferramenta computacional aqui desenvolvida em qualquer versão dessa plataforma, bem como a sua alteração para o uso com outros simuladores. Os testes realizados confirmam a eficiência da técnica e apontam as melhores abordagens para uma rápida e precisa estimação dos parâmetros. / This thesis presents a strategy for automatically obtaining hydrodynamic and transport parameters by means of the solution of inverse problems. Obtaining the parameters of a physical model represents a major problem in its calibration, and this is largely due to the difficulty associated to the field measurements of these parameters. In particular by modeling rivers and estuaries, the roughness height and the turbulent diffusion coefficient represent two of the most difficult parameters to measure or estimate. Here, an automated technique for estimation of these parameters through an inverse problem is applied to a model of the Macaé estuary, located in northern Rio de Janeiro. For this investigation, hydrodynamic and transport models were built in the MOHID platform, developed in the Technical University of Lisbon, which has had wide application in simulation of water bodies. A sensitivity analysis was performed regarding the model responses with respect to the parameters of interest. It was verified that salinity is a sensitive variable for both parameters. The inverse problem was then solved using various optimization methods by coupling the MOHID platform to optimization codes implemented in Fortran. The coupling was carried in a way to not change the MOHID source code, allowing the use of the computational tool here developed in any version of this platform, as well as its modification for use with other simulators. The tests confirm the efficiency of the technique and suggest the best approaches for a fast and accurate estimation of the parameters.

Devido à proximidade de grandes cidades, ecossistemas aquáticos são alvos de despejos industriais, esgotos domésticos e resíduos sólidos. Modelos matemáticos são largamente utilizados para entender os padrões de fluxo e o transporte de substâncias nestes ecossistemas. Neste trabalho, desenvolveu-se um modelo hidrodinâmico e de transporte de poluentes bidimensional para corpos d’água rasos que utiliza grades não estruturadas. O modelo foi desenvolvido em dois módulos: hidrodinâmico e de transporte de massa. O módulo hidrodinâmico utiliza as equações completas de Saint Venant para quantificação da circulação em corpos d’águas. No módulo hidrodinâmico aqui desenvolvido é feito um aperfeiçoamento de um esquema semi-implícito de diferenças finitas utilizado no modelo IPH-ECO para um esquema semi-implícito de volumes finitos. O módulo de transporte de massa, por sua vez, representa o transporte de um escalar conservativo sob influência dos processos físicos do meio. Dois esquemas foram empregados e comparados para a solução numérica da equação de transporte: um esquema de diferenças centrais e um esquema com limitador de fluxo. O modelo foi aplicado em um estudo de caso simplificado formado por uma lagoa circular e no rio Guaíba. O módulo hidrodinâmico foi avaliado através de testes de conservação de volume feitos comparando com os resultados encontrados nno modelo IPHECO na aplicação realizada na lagoa circular e na solução analítica no rio Guaíba. O módulo hidrodinâmico do modelo de grades não estruturadas (IPH-UnTRIM2D) mostrou-se mais conservativo que o modelo IPH-ECO em todos os cenários simulados na lagoa circular. A solução hidrodinâmica e de transporte do modelo de grades não estruturadas no estudo de caso idealizado no rio Guaíba representou meandros, sulcos, estreitamentos e alargamentos de canais com comprometimento com uma situação real. Quanto à simulação de transporte de escalares e poluentes, foi comprovada a eficiência do esquema com limitador de fluxo sobre os esquemas de diferenças centrais. Os testes de conservação de massa para diferentes cenários no rio Guaíba indicaram erros pouco significativos em ambos os esquemas quando comparados com a ordem de grandeza do ecossistema. / Aquatics ecossistems around large towns are always exposed waste domestic and industrial and solid residues. Mathematical models are widely used to understand flow patterns and transport of mass in these ecosystems. In this work a model was developed on unstructured two-dimensional complex dynamic model for shallow lakes. The model was divided in two modules: (a) hydrodynamic module, describing quantitative flows and water level; (b) a mass transport module. The hydrodynamic module solves Saint Venant equations to quantify water circulation in a water body. Futhermore, the hydrodynamic module was improved to semiimplicit finite volume scheme. The mass transport module represents transport of a conservative scalar under the physical processes influence. Two schemes were used and compared to the transport equation numerical solution : a central difference scheme and a flux limiter scheme. The model was applied in the circular lake and river Guaíba. The hydrodynamic module was validated by evaluating tests of volume conservation made comparing with results of the model IPH-ECO and analytical solution. It proved to be more conservative than the model IPH-ECO in all scenarios simulated in the circular lake. The hydrodynamic and transport solution in river Guaíba study case represented bends, grooves, and channels widening and narrowing. In the transport simulation, the flux limiter scheme was more efficient than central differences scheme. Futhermore, mass conservation tests for different scenarios in Guaíba indicated minor errors in both schemes when compared with ecosystem magnitude.

Este trabalho teve como objetivo analisar os possíveis efeitos de um aumento do nível do mar no padrão de correntes no Sistema Estuarino de Santos (SP), através do uso do modelo hidrodinâmico do Delft3D-FLOW, considerando diferentes taxas de aumento de acordo com o quinto relatório do Intergovernmental Panel on Climate Change (AR5). Os resultados indicam alteração na velocidade da corrente ao longo do domínio implementado, mais evidente durante a sizígia. Observa-se um aumento da magnitude da corrente em praticamente todo o Canal de São Vicente e redução da magnitude nos canais do Porto de Santos, canal de Bertioga e na região interna do Estuário de Santos. Os resultados do modelo numérico foram comparados com diversos dados coletados na região de estudo e verificou-se a correta representação dos padrões de circulação da área. Este estudo comprovou a importância de utilizarmos resultados de um modelo de maior escala (HYCOM) como forçantes, através da condição de fronteira denominada Riemann, permitindo que o modelo represente fenômenos de baixa frequência, que geram circulações e variações do nível do mar com a mesma ordem de grandeza dos efeitos da maré astronômica na área de estudo. / This study aimed to analyze the possible effects of a sea level rise in Sistema Estuarino de Santos (SP) circulation pattern, through the use of Delft3D-FLOW Hydrodynamic Model considering different elevation rates according to Intergovernmental Panel on Climate Change fifth report (AR5). The results indicate change in the current\'s magnitude, being more evident during the spring tide. There is an increase in the velocity practically all over São Vicente Channel and the reduction of magnitude in the Porto de Santos and Bertioga channels and in the inner region of the estuary. The numerical model results were compared with several measured data and it was verified that the model correctly representats the circulation patterns of the area. This study proved the importance of using results of a larger scale model (HYCOM), forced through a boundary condition called Riemann, allowing the model to represent low frequency phenomena, which generate sea level variations with the same order of magnitude of the effects of the astronomical tide in the study area.

Submitted by Wagner Maciel Castilho (castilho.w@gmail.com) on 2017-12-05T10:53:57Z No. of bitstreams: 1 TeseWagner.pdf: 4189219 bytes, checksum: dbc77f46b31b8872c5aafe68d8a5b7d0 (MD5) / Approved for entry into archive by Pamella Benevides Gonçalves null (pamella@feg.unesp.br) on 2017-12-05T12:38:21Z (GMT) No. of bitstreams: 1 catilho_wm_dr_guara.pdf: 4189219 bytes, checksum: dbc77f46b31b8872c5aafe68d8a5b7d0 (MD5) / Made available in DSpace on 2017-12-05T12:38:21Z (GMT). No. of bitstreams: 1 catilho_wm_dr_guara.pdf: 4189219 bytes, checksum: dbc77f46b31b8872c5aafe68d8a5b7d0 (MD5) Previous issue date: 2017-11-10 / Os fenômenos físicos envolvidos em colisões relativísticas de íons pesados foram estudados utilizando o código hidrodinâmico NeXSPheRIO, com condições iniciais flutuantes evento a evento, caracterizado por tubos de fluxo longitudinais de altas energias. Foram realizadas análises para o espectro de hádrons em função dos momentos transversais pelos cálculos das distribuições de partículas em sete janelas de centralidade de colisão, de acordo com dados publicados pela Colaboração STAR. A coletividade da matéria nuclear formada durante as colisões foi investigada em termos de correlações azimutais entre hádrons. Verificou-se que as estruturas da correlação dependem do momento transversal e da centralidade. Tais estruturas são atribuídas à interação entre um fundo dependente da centralidade e as flutuações de evento a evento na multiplicidade. Como resultado, as correlações observadas são apresentadas por um duplo pico na região do away-side que mudam continuamente para um pico único sob certas circunstâncias. Em termos do modelo do tubo periférico, o fundo é dominado pelo fluxo elíptico que aumenta com a descentralização. Pelos cálculos numéricos, a anisotropia do sistema é estudada quanto aos coeficientes de fluxo harmônico até a quarta ordem, de acordo com os dados da colaboração PHENIX. A anisotropia do fundo contribui para as correlações entre hádrons e pôde ser avaliada pelo método ZYAM. Após sua subtração na produção do sinal de correlação próprio, a estrutura de correlação resultante é uma função da centralidade, plano de evento e do momento transverso de acordo com os dados das colaborações PHENIX e STAR. Mostrou-se que a estrutura do duplo pico no away-side muda continuamente para um pico único com a redução na centralidade. Além disso, em relação à dependência do plano de eventos, a estrutura de duplo pico no away-side evolui para um pico único à medida que o ângulo azimutal das partículas gatilho se move do plano para o da direção fora do plano de eventos. Também foi mostrado que a estrutura resultante do away-side pode ser sensível ao esquema detalhado pelo método ZYAM, particularmente, a subtração do terceiro coeficiente harmônico. Por último, mas não menos importante, a magnitude da correlação diminuiu com o aumento do momento transversal das partículas associadas. Nossas análises apresentaram um acordo razoável com os dados experimentais publicados pela colaboração PHENIX e STAR do RHIC / The physical phenomena involved in relativistic heavy-ion collisions were studied using a hydrodynamic code, NeXSPheRIO, with event-by-event fluctuating initial conditions featured by high-energy longitudinal flux tubes. Analyses were carried out for hadron spectrum as a function of the transverse momentum by the calculations of particle distributions in seven collision centrality windows, according to the data published by the STAR Collaboration. Also, the collectivity of the nuclear matter formed during the collisions was investigated in terms of dihadron azimuthal correlations. The correlation structures were found to be dependent on the transverse momentum and centrality. They are attributed to the interplay between the centrality dependent background and the eventby-event multiplicity fluctuations. As a result, the observed correlations are featured by a double-peak in the away-side region that continuously change to a single peak under certain circumstances. In terms of the peripheral tube model, the background is dominated by the elliptic flow which increases with decreasing centrality. By numerical calculations, the anisotropy of the system is studied regarding the harmonic flow coefficients up to fourth order, in accordance with the existing RHIC data. The background anisotropy contributes to the dihadron correlations, and can be evaluated by the ZYAM method. After it is subtracted from the proper correlation yields, the resultant correlation structure is a function of centrality, event plane and transverse moment, in accordance with the data from the PHENIX and STAR Collaborations. It was shown that the double-peak structure in the away-side continually changes to a single peak as centrality decreases. In addition, concerning the event plane dependence, the double-peak structure in the away-side evolves into a single peak as the azimuthal angle of the trigger particles moves from the in-plane to the out-of-plane direction. It was also shown that the resultant away-side structure can be sensitive to the detailed schemes of the ZYAM method, particularly, to the subtraction of the third harmonic coefficient. Last but not least, the magnitude of the correlation was found to decrease with increasing transverse moment of the associated particles. Our analyses manifestly presented a reasonable agreement with the experimental data published by the PHENIX and STAR Collaborations

Devido à proximidade de grandes cidades, ecossistemas aquáticos são alvos de despejos industriais, esgotos domésticos e resíduos sólidos. Modelos matemáticos são largamente utilizados para entender os padrões de fluxo e o transporte de substâncias nestes ecossistemas. Neste trabalho, desenvolveu-se um modelo hidrodinâmico e de transporte de poluentes bidimensional para corpos d’água rasos que utiliza grades não estruturadas. O modelo foi desenvolvido em dois módulos: hidrodinâmico e de transporte de massa. O módulo hidrodinâmico utiliza as equações completas de Saint Venant para quantificação da circulação em corpos d’águas. No módulo hidrodinâmico aqui desenvolvido é feito um aperfeiçoamento de um esquema semi-implícito de diferenças finitas utilizado no modelo IPH-ECO para um esquema semi-implícito de volumes finitos. O módulo de transporte de massa, por sua vez, representa o transporte de um escalar conservativo sob influência dos processos físicos do meio. Dois esquemas foram empregados e comparados para a solução numérica da equação de transporte: um esquema de diferenças centrais e um esquema com limitador de fluxo. O modelo foi aplicado em um estudo de caso simplificado formado por uma lagoa circular e no rio Guaíba. O módulo hidrodinâmico foi avaliado através de testes de conservação de volume feitos comparando com os resultados encontrados nno modelo IPHECO na aplicação realizada na lagoa circular e na solução analítica no rio Guaíba. O módulo hidrodinâmico do modelo de grades não estruturadas (IPH-UnTRIM2D) mostrou-se mais conservativo que o modelo IPH-ECO em todos os cenários simulados na lagoa circular. A solução hidrodinâmica e de transporte do modelo de grades não estruturadas no estudo de caso idealizado no rio Guaíba representou meandros, sulcos, estreitamentos e alargamentos de canais com comprometimento com uma situação real. Quanto à simulação de transporte de escalares e poluentes, foi comprovada a eficiência do esquema com limitador de fluxo sobre os esquemas de diferenças centrais. Os testes de conservação de massa para diferentes cenários no rio Guaíba indicaram erros pouco significativos em ambos os esquemas quando comparados com a ordem de grandeza do ecossistema. / Aquatics ecossistems around large towns are always exposed waste domestic and industrial and solid residues. Mathematical models are widely used to understand flow patterns and transport of mass in these ecosystems. In this work a model was developed on unstructured two-dimensional complex dynamic model for shallow lakes. The model was divided in two modules: (a) hydrodynamic module, describing quantitative flows and water level; (b) a mass transport module. The hydrodynamic module solves Saint Venant equations to quantify water circulation in a water body. Futhermore, the hydrodynamic module was improved to semiimplicit finite volume scheme. The mass transport module represents transport of a conservative scalar under the physical processes influence. Two schemes were used and compared to the transport equation numerical solution : a central difference scheme and a flux limiter scheme. The model was applied in the circular lake and river Guaíba. The hydrodynamic module was validated by evaluating tests of volume conservation made comparing with results of the model IPH-ECO and analytical solution. It proved to be more conservative than the model IPH-ECO in all scenarios simulated in the circular lake. The hydrodynamic and transport solution in river Guaíba study case represented bends, grooves, and channels widening and narrowing. In the transport simulation, the flux limiter scheme was more efficient than central differences scheme. Futhermore, mass conservation tests for different scenarios in Guaíba indicated minor errors in both schemes when compared with ecosystem magnitude.