Artigos de Revista
Info
Physical modelling as a fundamental tool for the design of harbours and maritime structures
In this paper, a brief description of the role of physical modelling to support the design of harbours and maritime structures is presented, followed by an example of the tests carried out for the studies of the expansion of Rabo de Peixe fi shing harbour, at São Miguel
Ano: 2015
Número Páginas:
135-145pp.
Autor(es): Reis, M. T.; Silva, L. G.; Neves, M. G.; Lemos, R.; Capitão, R.; Fortes, C. J. E. M.
Revista: PIANC Yearbook 2014
Keywords: Rabo de peixe harbour; Maritime structures; Harbours; Physical modelling
Info
Propuesta de estudio de la morfología determinada por el oleaje en playas fluviales. "La Concordia" en el río Uruguay, un caso práctico
Tanto el perfil como la forma en planta de las playas fluviales vienen establecidos por numerosos factoresentre los que cabría destacar las características del sedimento y la dirección y módulo de la velocidad de la corriente en el río. Sin embargo, en ciertos casos, lejos incluso de su estuario o desembocadura, esta morfología puede estar determinada principalmente por la incidencia del oleaje pese a no tratarse de un entorno marítimo.Se presenta aquí una metodología cualitativa que nos permitirá reconocer los aspectos que influyen en la existencia de este tipo de morfologías. Para ello, se estudiará un caso particular, el de la playa de La Concordia ubicada en la margen izquierda del río Uruguay a 70 km de su desembocadura en el estuario del Rio de la Plata. Asimismo se analizará la forma en planta y se verificará que la misma está determinada por la acción del oleaje.
Ano: 2015
Número Páginas:
49
Autor(es): Muñoz-Perez, J.; Solari, S.; Teixeira, L; Alonso, R.; Neves, M. G.
Revista: Geo-temas
Volume:
15.
Keywords: Oleaje; Difracción; Refracción; Playas fluviales
Info
Sinking of concrete modules into a sandy seabed: a case study
Three submerged coastal structures on the sandy seabed at SantaMaria delMar (SMM) Beach in Southwest Spain were monitored during the sixmonths after installation, starting inNovember of 2005. The monitoring evaluated the self-burial phenomenon of the 40 precast modular concrete elements of the three structures. Assessment included scouring around the structures, their vertical movement, the sinking velocity of the modules and the resulting beach profiles. Pressure sensors and precision topographies permitted continuous monitoring of the self-burial process for the first time in a full-scale case. Comparison of two bathymetries (one performed immediately prior to installation and another performed six months after installation) indicated slight accretion of the seabed,whereas intermediate topographic surveys indicated extensive scouring around the structures. The three structures began sinking into the sandy bottomplatform immediately after placement and continued until they reached the rocky bottom. This study implemented a new methodology using pressure sensors attached to the concrete modules to monitor the structural sinking. Unexpectedly, the average sinking speed was extremely rapid at approximately 3
Ano: 2015
Número Páginas:
26
Autor(es): Muñoz-Perez, J.; Khan-Mozahedy, A.B.M.; Neves, M. G.; Tejedor, B.; Gomez-Pina, G.; Campo, J.; Negro, V.
Revista: Coastal Engineering
Volume:
99.
Keywords: Self-buried; Low crested breakwaters; Liquefaction; Scouring; Sandy seabed
Info
Stability formula and failure probability analysis of wave-dissipating blocks considering wave breaking
The effects ofwave steepness andwave breaking on the stability of wave-dissipating blocks placed on caisson breakwaters are examined by hydraulic experiments under conditions with and without wave breaking. It is confirmed that there is a clear dependence of the block stability on the wave steepness through the wave-breaking type, such as plunging and spilling breakers. The stability parameter in thespilling breaker case is larger than that in the plunging breaker case. For non-breaking wave conditions, there is no effect of the wave steepness. An empirical stability formula is proposed taking into account the effects of the wave steepness and the wave-breaking type. Finally, by using the stability formula, failure probabilities of the wave-dissipating blocks having a nominal diameter are evaluated by the first-orderreliability method of the Level II reliability analysis.
Ano: 2015
Número Páginas:
45-54pp.
Autor(es): Mase, H.; Yasuda, T.; Reis, M. T.; Karunarathna , H.; Yang , J.
Revista: Journal of Ocean Engineering and Marine Energy
Editor: Springer
Volume:
Vol. 1(1).
Keywords: Failure probability; Reliability analysis; Wave-breaking effect; Stability formula; Wave-dissipating block stability
Info
Tools for forecasting and evaluation of structural risk on maritime works
This paper describes recent advances on maritime structures
Ano: 2015
Número Páginas:
146-158pp.
Autor(es): Fortes, C. J. E. M.; Lemos, R.; Reis, M. T.; Capitão, R.; Neves, M. G.; Silva, L. G.; Poseiro, P.
Revista: PIANC Yearbook 2014
Editor: PIANC
Keywords: ANOSOM; Maritime inspections; Warning system; Risk evaluation; Maritime structures
Info
Wave propagation on a flume: numerical simulation
This paper presents the numerical simulations done by using the waves2Foam, an OpenFOAM® library, to simulate the propagation of regular waves without breaking in a three-dimensional flume. The numerical code solves the unsteady Navier-Stokes equations and uses a Volume-of-Fluid (VoF) method to identify the free surface. A regular incident wave with a 1.5s period and 0.1m wave height, was considered. This is one of the conditions from the wide range of wave flume tests conducted at the National Laboratory for Civil Engineering (LNEC) whose objective was to analyze the hydrodynamics of wave transformation and wave breaking for different incident conditions over a variable bathymetry. Comparisons are made between the numerical and the experimental results. These comparisons include time-series of wave-gauges records at several locations along the flume and the corresponding amplitude spectra; significant wave height and average period evolution along the flume; time-series of the velocity components at one section of the flume, measured at the middle of the water column; and hodograph representation of the velocity components, in the middle of the water column, in the xy, xz, and yz planes along the flume. It was found that the numerical results obtained are close to the experimental data. The observed differences are attributable to numerical inaccuracies as well as the differences between the wave generation method in the numerical and experimental tests.
Ano: 2015
Número Páginas:
95-102pp.
Autor(es): Conde, J. M.; Roberto, P.T.S; Fortes, C. J. E. M.
Revista: Engenharia Térmica (Thermal Engineering)
Volume:
Vol 14 - nº 1.
Keywords: Wave propagation; Free-surface; OpenFOAM; Navier
Info
A modelação física no apoio ao projeto de obras marítimas
Neste artigo, apresenta-se uma breve descrição do papel da modelação física no apoio ao projeto de obras marítimas, seguida de um exemplo de ensaios realizados no âmbito da ampliação do porto de pesca de Rabo de Peixe, São Miguel
Ano: 2014
Número Páginas:
26-30pp.
Autor(es): Fortes, C. J. E. M.; Reis, M. T.; Neves, M. G.; Silva, L. G.; Capitão, R.; Lemos, R.
Revista: Construção Magazine
Volume:
Nº 62.
Info
Avaliação da probabilidade de ocorrência de galgamentos não admissíveis no posto 2 do Terminal de granéis líquidos do porto de Sines
Descrevem-se os procedimentos para avaliar o caudal médio galgado sobre um trecho de um quebra-mar de taludesa partir das caraterísticas da agitação marítima incidente no mesmo trecho. Apresenta-se uma metodologia paraestimação de valores extremos daquele caudal, partindo das séries temporais da mesma variável, e ilustram-se essesprocedimentos num trecho do molhe oeste do porto de Sines que protege o posto 2 do terminal de granéis líquidos daação da agitação marítima. Este regime de extremos do caudal médio galgado é fundamental para a avaliação do riscode ocorrência de galgamentos não admissíveis no posto 2.
Ano: 2014
Número Páginas:
95-113pp.
Autor(es): Reis, M. T.; Santos, J. A.; Capitão, R.; Neves, D. R.; Fortes, C. J. E. M.
Revista: Territorium
Volume:
21.
Keywords: Porto de sines; Caudais médios extremos; Redes neuronais; Galgamento
Info
Interação entre ondas oceânicas e fundo marinho: Resultados na plataforma continental setentrional do Rio Grande do Norte
A plataforma continental setentrional do Rio Grande do Norte caracteriza-se por apresentar uma variedade de feições de fundo marinho que variam desde 10 cm até quilômetros, sobretudo submersos nas porções mais rasas da plataforma e indicam processos de retrabalhamento. Além do caráter singular dos processos, formas efáceis de fundo diferenciados a qual esta plataforma está submetida, compreendese a importância do conhecimento da dinâmica das ondas, uma vez que estas são consideradas o agente forçador dominante na hidro e morfodinâmica costeira. O objetivo principal deste trabalho é caracterizar o clima de ondas por meio da modelagem numérica, através do modelo de ondas de superfície do mar, SWAN, e analisar sua interação com o fundo marinho raso da plataforma continentalsetentrional do Rio Grande do Norte. Foram realizadas várias medições de campo entre 2011 e 2012, com o objetivo de se obter os parâmetros de ondas: altura signifi cativa (HS), período médio (Tmed), direção (DIR) para servir de comparação e validação da modelagem.Os resultados apontam que o SWAN conseguiu, no geral, reproduzir adequadamente o clima de ondas para este litoral. Entretanto, ao analisar os resultados, confi rmou-se que as alturas signifi cativas apresentam resultados mais aceitáveis em relação aos demais parâmetros. A interação da onda com o fundo marinho revela alteração das alturas significativas e direções ao propagarem-se sobre as formas de fundos presentes, como as rochas praiais, os campos de dunas longitudinais, transversais e subaquáticos. Tais resultados mostram que nesta plataforma continental o mecanismo de dissipação dominante à propagação das ondas nestas águas rasas é o atrito de fundo.
Ano: 2014
Número Páginas:
371-391pp.
Autor(es): Matos, M. F.; Amaro, V. E.; Fortes, C. J. E. M.; Scudelari, A. C.
Revista: Revista Brasileira de Geomorfologia
Volume:
v. 15, nº 3.
Keywords: Rio grande do norte; Fundo marinho raso; Modelagem de ondas
Info
INTERAÇÃO ENTRE ONDAS OCEÂNICAS E FUNDO MARINHO: RESULTADOS NA PLATAFORMA CONTINENTAL SETENTRIONAL DO RIO GRANDE DO NORTE
A plataforma continental setentrional do Rio Grande do Norte caracteriza-se porapresentar uma variedade de feições de fundo marinho que variam desde 10 cm atéquilômetros, sobretudo submersos nas porções mais rasas da plataforma e indicamprocessos de retrabalhamento. Além do caráter singular dos processos, formas efáceis de fundo diferenciados a qual esta plataforma está submetida, compreendesea importância do conhecimento da dinâmica das ondas, uma vez que estas sãoconsideradas o agente forçador dominante na hidro e morfodinâmica costeira.O objetivo principal deste trabalho é caracterizar o clima de ondas por meio damodelagem numérica, através do modelo de ondas de superfície do mar, SWAN,e analisar sua interação com o fundo marinho raso da plataforma continentalsetentrional do Rio Grande do Norte. Foram realizadas várias medições de campoentre 2011 e 2012, com o objetivo de se obter os parâmetros de ondas: altura signifi cativa (HS), período médio (Tmed), direção (DIR) para servir de comparação e validação da modelagem.Os resultados apontam que o SWAN conseguiu, no geral, reproduzir adequadamente o clima de ondas para estelitoral. Entretanto, ao analisar os resultados, confi rmou-se que as alturas signifi cativas apresentam resultados maisaceitáveis em relação aos demais parâmetros. A interação da onda com o fundo marinho revela alteração das alturassignifi cativas e direções ao propagarem-se sobre as formas de fundos presentes, como as rochas praiais, os camposde dunas longitudinais, transversais e subaquáticos. Tais resultados mostram que nesta plataforma continental omecanismo de dissipação dominante à propagação das ondas nestas águas rasas é o atrito de fundo.
Ano: 2014
Número Páginas:
371-391pp.
Autor(es): Fortes, C. J. E. M.; Scudelari, A. C.; Matos, M. F.; Amaro, V. E.
Revista: Revista Brasileira de Geomorfologia
Volume:
v. 15, nº 3.
Keywords: Rio grande do norte; Fundo marinho raso; Modelagem de ondas
|
Comunicações
Info
Experimental characterization of ship motions induced by passing ships
Using physical modelling, the movements of an otherwise motionless ship induced by the wake waves generated by another ship that sails in its vicinity are analyzed.The tests were carried out at one of the wave tanks of the Portuguese Civil Engineering Laboratory (LNEC). The moving ship is a self-propelled scale model of the
Ano: 2017
Número Páginas:
12p.
Autor(es): Pedro, F.; Santos, J. A.; Hinostroza , M.; Pinheiro, L.; Fortes, C. J. E. M.
Editor: Artes Gráficas J. Martínez, S.L
Keywords: Ship motions; Interaction between ships; Physical model tests; Wake waves
Info
Forecast, risk assessment and alert system for ships in ports
This work presents the development of a numerical integrated system to help prevent emergency situations and minimize risks associated with sea waves
Ano: 2017
Número Páginas:
13p.
Autor(es): Pinheiro, L.; Fortes, C. J. E. M.; Santos, J. A.; Lopes, P.; Poseiro, P.
Editor: Artes Gráficas J. Martínez, S.L
Keywords: Waves; Alert system; Risk assessment; Moored ships
Info
Galgamento em estruturas verticais: comparação de dados de ensaios com resultados do modelo numérico AMAZON e de fórmulas empíricas
Neste trabalho estuda-se o galgamento em estruturas marítimas verticais impermeáveis através da aplicação do modelo numérico AMAZON, que é um modelo bidimensional, baseado na resolução das equações não lineares de águas pouco profundas. Para validar a sua aplicação, os resultados de caudal médio galgado obtidos com o modelo foram comparados com os resultados de ensaios em modelo físico bidimensional realizados anteriormente e disponíveis na base de dados do projeto europeu CLASH. Os resultados foram ainda comparados com os das fórmulas empíricas / semi-empíricas de Van der Meer & Bruce (2014), Goda (2008) e de Hedges & Reis (Reis et al., 2013). A fórmula de Van der Meer & Bruce apresenta melhores aproximações aos caudais médios medidos que constam na base de dados do CLASH, seguida da fórmula de Goda e dos resultados do AMAZON, embora estes dois últimos apresentem valores bastante próximos. Assim, enquanto não se fizerem melhorias no modelo e/ou não se estender a outros ensaios a análise feita no âmbito deste trabalho, os resultados sugerem que não há vantagem na sua utilização face às formulações disponíveis, tendo também em conta o tempo computacional significativo exigido para executar o AMAZON.
Ano: 2017
Número Páginas:
17p.
Autor(es): Reis, M. T. L. G. V.; Neves, M. G.; Gaspar, J.J.J.; Hu, K.
Editor: PIANC
Keywords: Modelo numérico AMAZON; Estruturas verticais; Galgamento
Info
Galgamento num quebra-mar de talude em cenários de alterações climáticas: resultados obtidos em modelo físico e com fórmulas empíricas
O projeto HYDRALAB+ (H2020-INFRAIA-2014-2015) reúne uma rede Europeia de institutos de hidráulica para o estudo de estratégias de adaptação às alterações climáticas através de ensaios experimentais e cujos resultados são essenciais para a calibração e validação de modelos numéricos, fórmulas empíricas, redes neuronais, etc. Os ensaios experimentais desenvolvidos no LNEC, no âmbito do HYDRALAB+, Tarefa 8.2 do RECIPE tem como objetivo a análise do dano num quebra-mar de talude, do galgamento e do espraiamento para quatro cenários de temporal correspondentes a situações extremas associadas a cenários de alterações climáticas. Os ensaios foram realizados num canal de ondas irregulares no qual foiimplementada uma secção transversal do quebra-mar, com um manto constituído por enrocamento.No presente estudo comparam-se os valores do galgamento medidos nos ensaios em modelo físico com resultados estimados por fórmulas empíricas. Testou-se o efeito de um temporal num cenário de alterações climático, através do aumento da altura significativa de onda com o aumento do nível de água. Foi testado um conjunto de condições de agitação incidentes em combinação com dois níveis de maré.A análise dos resultados é feita com o intuito de estimar com precisão a ocorrência de galgamento em quebra-mares, um fenómeno que ainda se torna mais relevante em cenários de alterações climáticas. Com efeito, nestes cenários é previsível a subida do nível médio do mar, bem como o aumento da intensidade e frequência de temporais, o que pode dar origem a fenómenos de espraiamento/galgamento cada vez mais frequentes, tornando assim cada vez mais vulneráveis as estruturas de proteção costeiras e portuárias, bem como as zonas por elas protegidas.
Ano: 2017
Número Páginas:
13p.
Autor(es): Mendonça, A. C.; Fortes, C. J. E. M.; Reis, M. T. L. G. V.; Lemos, R.; Neves, M. G.; Ramos, A.; Capitão, R.
Editor: PIANC
Keywords: Fórmulas empíricas; Modelação física; Alterações climáticas; Galgamento; Quebra-mar
Info
HIDRALERTA: o sistema e as ferramentas de validação
Este trabalho pretende mostrar os procedimentos utilizados pelo sistema HIDRALERTA para a previsão dos galgamentos costeiros, bem como as ferramentas que têm sido utilizadas para a sua validação. Foram reproduzidas com o sistema as situações de tempestades passadas, para as quais houve relatos documentados de galgamentos significativos e que causaram danos. Esta simulação permitiu perceber a fiabilidade global do sistema e ajustar os limites pré- estabelecimentos para o acionamento de alertas. Com o sistema operacional desde 2015, a documentação que tem sido levantada, através de fotografias e vídeos, pelas autoridades e comunidade locais aquando da previsão da ocorrência de tempestades, têm sido um grande contributo para o melhoramento do sistema. Foi instalado no campo, um sistema de videovigilância que permite detetar e captar automaticamente imagens de galgamentos para comparar com os resultados do sistema. De forma complementar, para avaliação da agitação marítima, utilizam-se as medições da boia ondógrafo instalada junto ao porto e para a validação do fenómeno do galgamento, compararam-se as medições obtidas nos ensaios em modelo físico do perfil tipo do quebra-mar sul, com os resultados da ferramenta neuronal NN_OVERTOPPING2 (Coeveld et al., 2005) (atualmente incorporada no sistema) e com o modelo numérico SWASH (Zijlema et al., 2011) (que se encontra em processo de implementação).
Ano: 2017
Número Páginas:
20p.
Autor(es): Poseiro, P.
Editor: PIANC
Keywords: Porto da Praia da Vitória; Previsões; Eventos históricos; Validação; Sistema HIDRALERTA
Info
Measuring damage in physical model tests of rubble mounds
This paper studies novel ways to evaluate armour damage in physical models of coastal structures. High-esolution damage data for reference rubble mound breakwaters obtained under the HYDRALAB+ joint-esearch project are analysed and discussed. These tests are used to analyse the way to describe damage, the influence of the sequence of testing, and touches on the possible influence of sea level rise. Results of two test programmes were used. Firstly, 3D physical model tests carried out at the University of Porto, in cooperation with Deltares, were used. Here a wide breakwater trunk was used for statistical reasons. Additionally, 2D test results from LNEC were analysed. Tests for a sea level rise scenario resulted in less damage to the seaside slope. In addition, clear differences between
Ano: 2017
Número Páginas:
10p.
Autor(es): Hofland, B.; Rosa-Santos, P.; Taveira-Pinto, F.; Almeida, E.; Lemos, R.; Mendonça, A. C.; Fortes, C. J. E. M.
Editor: ICE - Institution of Civil Engineers
Keywords: Damage depth; Damage parameter; Breakwaters; Rubble-mound
Info
Modelação física e numérica da interação hidrodinâmica entre dois navios
Nesta comunicação, analisa-se, recorrendo à modelação física e à modelação numérica, a interação hidrodinâmica entre dois navios, um dos quais parado, o navio
Ano: 2017
Número Páginas:
20p.
Autor(es): Pedro, F.; Pinheiro, L.; Santos, J. A.; Fortes, C. J. E. M.; Hinostroza , M.
Editor: PIANC
Keywords: Interação navio-navio; Modelo numérico; Modelo físico; Movimentos dos navios; Ondas de esteira
Info
Numerical analysis of the distribution of incident wave energy on an onshore oscillating water column wave energy converter
In this study, a numerical analysis to evaluate the energy efficiency of an Oscillating Water Column (OWC) wave energy converter is carried out. The device geometry is based on the pilot plant installed on Pico Island, Azores, Portugal. The methodology consists in calculating the percentage of the incident energy that is distributed along the wave energy conversion. The FLUENT® numerical model, which is based on the Reynolds-Averaged Navier-Stokes (RANS) equations, is used. The Volume of Fluid (VoF) method is employed to take into account the water and air fluids. The case study comprises a flume 10 m deep with an onshore OWC at its end. The distribution of the incident wave energy (1 m wave height with a 9 s period) is determined in terms of reflected wave energy, water column oscillation potential energy, pneumatic energy and energy losses due to viscosity and turbulence dissipations. Results show that the incident wave energy is distributed around 37% for the sum of reflected wave energy and viscosity energy loss, 52% for pneumatic energy and 11% for water column oscillation potential energy.
Ano: 2017
Número Páginas:
13p.
Autor(es): Gaspar, L. A.; Teixeira, P.; Didier, E.; Neves, M. G.
Editor: Não disponível
Keywords: Volume of fluid; Oscillating water column device; Finite volume method; Numerical simulation; Wave energy
Info
Numerical analysis of the influence of air compressibility effects on the oscillating water column water energy converter chamber
The most studied device used for extracting wave energy is the Oscillating Water Column (OWC). In general, numerical simulations of these cases by means of models based on Reynolds-Averaged-Navier-Stokes (RANS) equations adopt the Volume of Fluid (VoF) method to take into account the water and air whose flows are incompressible. The aim of this study is to investigate the influence of the compressibility effect on the air inside the chamber by means of FLUENT® numerical model. Both water and air flows are still considered incompressible, but, at every instant, a pressure based on an analytical equation that takes into account the turbine characteristic relation (kt) of a Wells turbine and the isentropic transformation of the air is imposed on the top boundary of the chamber. Incident waves with period of 9 s and 1 and 2 m high are imposed in a flume 10 m deep with an onshore OWC at its end (kt = 240 Pa s m-3). Results show differences between incompressible and compressible air around 12% in the air pressure amplitude inside the chamber and about 21% in the average pneumatic power.
Ano: 2017
Número Páginas:
14p.
Autor(es): Gonçalves, R.; Teixeira, P.; Didier, E.; Torres, F.
Editor: Não disponível
Keywords: Volume of fluid; Finite volume method; Numerical simulation; Oscillating water column device; Wave energy
Info
Overtopping events in breakwaters: comparison of 2D physical experiments and empirical formulae
Reliable prediction of wave run-up/overtopping and structure damage is a key task in the design and safety assessment of coastal and harbor structures. Run-up/overtopping and damage must be below acceptable limits, both in extreme and in normal operating conditions, to guarantee the stability of the structure as well as the safety of people and assets on and behind the structure. The mean-sea-level rise caused by climate change and its effects on wave climate may increase run-up/overtopping events and make the existing coastal/harbour structures more vulnerable to damage. Accurate estimates, through physical modelling, of the statistics of overtopping waves for a set of climate change conditions, are needed. Project HYDRALAB+ (H2020-INFRAIA-2014-2015) gathers an advanced network of environmental hydraulic institutes in Europe, which provides access to a suite of environmental hydraulic facilities. They play a vital role in the development of climate change adaptation strategies, by allowing the direct testing of adaptation measures and by providing data for numerical model calibration and validation. The use of physical (scale) models allows the simulation of extreme events as they are now, and as they are projected to be under different climate change scenarios. The experimental work developed at LNEC within HYDRALAB+, TASK 8.2 of RECIPE, entitled
Ano: 2017
Número Páginas:
3749-3757pp.
Autor(es): Mendonça, A. C.; Lemos, R.; Reis, M. T. L. G. V.; Fortes, C. J. E. M.; Neves, M. G.; Ramos, A.; Capitão, R.
Editor: IAHR & USAINS HOLDING SDN BHD
Keywords: Hydralab; Empirical formulae; Physical experiments; Breakwaters; Overtopping
|
|
Capítulos de Livros
Info
Multi-functional artificial reefs for coastal protection
N/A
Ano: 2009
Número Páginas:
153-210 pp.
Autor(es): Voorde, M.; Antunes do Carmo, J.; Neves, M. G.
Editor: Ch. A. Hudspeth, T. E. Reeve
|
Teses de Doutoramento
Info
Physical and Numerical Modelling of Wave Propagation Over Vegetation
Os efeitos das alterações climáticas aumentaram a necessidade de investigar a capacidade da vegetação em atenuar a força das ondas, nas zonas costeiras. Apoiando-se em modelação física, analítica e numérica da propagação de ondas sobre vegetação, este estudo teve como objetivos: i) analisar a diferença na dissipação de ondas sobre campos de vegetação rígida e flexível, ii) caracterizar o movimento de vegetação flexível, iii) estimar o coeficiente de arrasto de vegetação rígida e vegetação flexível, e iv) analisar a influência da flexibilidade (movimentoda vegetação) e da força de inércia no coeficiente de arrasto.
Ano: 2022
Autor(es): Reis, R.
Editor: IST
Keywords: Propagação/Dissipação de ondas; Análise de imagem; Movimento da vegetação; Coeficiente de arrasto; Vegetação flexível
Info
Modelação numérica da propagação de ondas e do comportamento de navios amarrados no interior de portos
Na presente dissertação, desenvolveu-se uma metodologia para caracterização do movimento de um navio amarrado no interior de portos na presença de um campo de ondas. Esta metodologia consiste em utilizar vários modelos numéricos para caracterizar convenientemente: a) a propagação das ondas desde o largo até ao local de acostagem; b) a modificação das características das ondas no interior do porto; c) a ação das ondas sobre o navio; d) a resposta do navio; e) as forças resultantes no sistema de marração. Os modelos utilizados são:· SWAN: modelo não linear espectral para simular a propagação de ondas em zonas costeiras extensas;· DREAMS: modelo linear para simular a propagação de ondas em zonas costeiras abrigadas;· BOUSS-WMH: modelo não linear para simular a propagação de ondas regulares e irregulares em zonas costeiras abrigadas;· GMALHA: gerador de malhas de elementos finitos triangulares;· WAMIT®: que resolve, no domínio da frequência, os problemas de radiação e de difração na interação de um corpo flutuante livre com as ondas;· NPP: gerador de painéis planos sobre o casco do navio;· HYDRO: que transforma os resultados da interação do corpo flutuante livre com as ondas do domínio da frequência para o domínio do tempo;· BAS: que resolve as equações de movimento de um navio amarrado.O desempenho do modelo BOUSS-WMH, em situações de grande não linearidade, foi analisado através da simulação de quatro casos de teste também simulados em laboratório. O BOUSS-WMH provou ser um modelo robusto capaz de representar a deformação da superfície livre de forma adequada, embora, em zonas com forte não linearidade, tenha apresentado algumas limitações, nomeadamente na formaçãode harmónicas de ordem superior à terceira. Verificou-se que os resultados numéricos tendiam a sobrestimar os valores da altura de onda em zonas pouco profundas.Para minorar estas limitações e tornar o modelo mais robusto, fiável e adequado para ser aplicado em casos reais, foram adicionados e implementados novos termos às equações do modelo BOUSS-WMH, nomeadamente, a reflexão parcial, o atrito de fundo e a rebentação (fenómenos físicos inerentes às águas pouco profundas e zonas confinadas).A reflexão parcial foi implementada a partir da adição de um termo de viscosidade na equação de conservação da quantidade de movimento. O modelo com reflexão parcial foi aplicado a dois casos de teste, sendo um deles um caso real (Porto de Vila do Porto, ilha de Santa Maria, Açores), e os resultados foram comparados com medições em ensaios em modelo físico. Verificou-se que a implementação da reflexão parcial diminui em 44% o erro médio quadrático face à simulação com fronteiras totalmente refletoras.O atrito de fundo foi implementado a partir da adição de um termo representativo da tensão tangencial de fundo à equação de conservação da quantidade de movimento. O modelo com atrito de fundo foi aplicado a dois casos de teste, sendo um deles um caso real (Praia de Faro). Verificou-se que a introdução do atrito de fundo dissipa energia como pretendido e reduz as instabilidades numéricas.A rebentação foi implementada com a introdução de um termo de viscosidade turbulenta. Dois critérios distintos foram introduzidos para definir o início e o final da rebentação: a) FSA, baseado na aceleração da elevação da superfície livre; b) RTFN, baseado no número Froude relativo da cava. O modelo com rebentação foi aplicado a três casos de teste estudados em laboratório. A comparação dos resultados do modelo com as medições dos ensaios físicos permite dizer que, de um modo geral, o modelo numérico simula muito bem a localização da rebentação, bem como a dissipação de energia a ela associada. O critério de rebentação baseado no método de RTFN tende a sobrestimar a altura das ondas, mas apresenta índices de concordância maiores do que os obtidos com o critério FSA. A vantagem deste método reside no menor número de parâmetros a calibrar. Os resultados alcançados foram muito bons, dada a complexidade do fenómeno de rebentação das ondas e a simplicidade dos termos utilizados na sua modelação.O gerador GMALHA também foi alvo de melhoramentos nos algoritmos de geração, de forma a possibilitar a existência de ilhas no interior do domínio de cálculo.O modelo numérico WAMIT® foi desenvolvido para aplicações em mar aberto sem muitas estruturas refletoras próximas do navio. Todas as estruturas devem ser discretizadas, o que aumenta exponencialmente o esforço computacional exigido, e são obrigatoriamente totalmente reflectoras. Ora esta é uma limitação importante quando se está no interior de um porto, primeiro devido às dimensões e depois devido às fronteiras parcialmente refletoras. Assim, o modelo WAMIT® foi modificado de modo a incluir as características das ondas modeladas previamente com um dos modelos de propagação. Foram desenvolvidas duas novas rotinas HASK_D e HASK_B, respetivamente para o modelo linear DREAMS e para o modelo não linear BOUSS-WMH. Estas rotinas e as modificações ao modelo WAMIT® foram validadas com a aplicação a casos simples onde o modelo WAMIT® original também podia ser aplicado. A presença de um campo de ondas difratadas pela presença de uma estrutura fixa foi simulada com sucesso. Existem diferenças pouco significativas, principalmente nas altas frequências. Uma explicação para estas diferenças relaciona-se com o facto de que as ondas com períodos baixos apresentam maiores instabilidades numéricas à medida que diminui o número de pontos por comprimento de onda na malha de elementos finitos.Com o intuito de facilitar o tratamento da informação, que, num estudo deste género, tem um volume e complexidade consideráveis, foram desenvolvidas interfaces para cada um dos modelos e estes foram acoplados numa ferramenta integrada
Ano: 2015
Autor(es): Pinheiro, L.
Keywords: Elementos finitos; Navios amarrados; Propagação de ondas; SWAMS; Forças de difração; Modelos numéricos; Portos
Info
Modelação Numérica da Rebentação das Ondas. Sistema Integrado de Modelação da Agitação Marítima
Esta dissertação divide-se em duas partes principais, uma dedicada ao desenvolvimento e validação de um novo modelo numérico para a propagação de ondas fracamente não lineares com a inclusão da dissipação de energia por rebentação, e uma dedicada ao desenvolvimento de um sistema integrado de modelação da agitação marítima em zonas costeiras e portuárias, GUIOMAR, para a utilização simplificada de modelos de propagação de ondas e para o apoio à tomada de decisão em estudos de engenharia costeira (quer correntes, quer situações de emergência).Na primeira parte, as equações de Beji e Nadaoka (1997a,b), são extendidas de forma a incluírem dissipação de energia por rebentação, através da consideração das tensões de Reynolds na equação da conservação da quantidade de movimento e da dedução de uma formulação alternativa da equação não linear de declive suave dependente do tempo. O modelo desenvolvido neste trabalho é aplicado a três testes de ensaios em modelo físico, que abrangem praias de declive constante e um perfil barra-fossa. Os testes efectuados permitiram seleccionar o critério de rebentação e definir regras para a determinação dos valores para as condições críticas de início e de fim da rebentação, calibrar o valor a atribuir ao parâmetro dos termos dissipativos e estabelecer regras para a sua definição, assim como para avaliar o desempenho do modelo. Foram obtidos bons resultados, embora o empolamento da onda tenha sido subestimado. Verificou-se uma estreita relação entre o valor a utilizar para a condição crítica de início da rebentação e a declividade da onda à entrada do domínio de cálculo, e entre o valor óptimo para o parâmetro dos termos dissipativos e o número de Ursell à entrada do domínio.Na segunda parte, apresenta-se a estrutura do sistema GUIOMAR, assim como as interfaces suas constituintes. O sistema GUIOMAR permite a utilização simplificada de modelos de agitação marítima, incluindo a construção, manipulação e visualização de dados e resultados. Estes factores contribuem para a sua afirmação como uma ferramenta de apoio à tomada de decisão em estudos de engenharia costeira e portuária. Foca-se também a avaliação das variabilidades de origem geográfica nos resultados dos modelos, nomeadamente as associadas à escala e resolução dos dados batimétricos utilizados, aos métodos de interpolação mais comuns e à subjectividade inerente à definição das linhas de costa. Da análise do comportamento dos modelos de agitação marítima constatou-se a necessidade de ter em atenção no respectivo estudo, aos dados de batimetria utilizados e às interpolações a que se recorre. São delineadas sugestões/recomendações para um futuro utilizador do sistema GUIOMAR aquando da utilização de métodos de interpolação.
Ano: 2010
Autor(es): Zózimo, A. C.
Keywords: Zonas portuárias e costeiras; Variabilidades de origem geográfica; Sistema integrado; Sistema de informação geográfica; Modelação numérica da rebentação das ondas; Equação não linear de declive suave dependente do tempo
Info
Contribution to the design of Multi-Functional Artificial Reefs
Portugal is one of many countries in the world that suffers from coastal erosion.Conventional ways of protecting a coastline appear to entail some disadvantages. Aninnovative and interesting way of protecting a local coastal zone by means of multifunctionalartificial reefs avoids some of them. A multi-functional artificial reef is asubmerged breakwater which, besides the helping to protect the local coastline, canoffer other purposes; in particular, it may enhance the surfing possibilities and theenvironmental value of the local area. The structure has several positive side-effects:first, it provides an unimpaired visual amenity; second, it offers tourist and economicbenefits by improving the surfing conditions; third, it can enhance the environmentalvalue of the area where it is built, and fourth, with a proper design the down drifterosion can be minimal.An optimal reef design for the Leirosa beach, located to the south of Figueira da Foz,midway along the Portuguese West Atlantic coast, is investigated. In order to achievethis optimal design several steps were conducted. A preliminary design, as a result of atheoretical study, achieved step by step, is proposed for a multi-functional artificial reefmaking use of the theory and of the state of the art. The proposed reef geometry is usedas initial design in the numerical and physical tests, which were executed to analyze thecapacity of a multi-functional reef breakwater to protect the local coastline of Leirosa,Portugal, and to increase the local surfing possibilities. Further, a physical and anumerical 2DV study have been executed to get insight in the influence of the length ofthe reef and in the submergence on the Iribarren number, which is an indication for thebreaker type. In this study, also the optimum height for the reef has been defined.Subsequently, in order to be able to investigate the three design parameters: the reefangle, the geometry of the reef (with or without a platform) and the horizontaldimensions in 2DH for the optimal geometry, the Boussinesq-type COULWAVE modelwas validated; finally, a 2DH numerical study has been conducted with COULWAVEto define the best values for these parameters.In terms of surfability, and for the conditions of the local coastline of Leirosa, thefollowing values for the main parameters were found: a reef angle of 66°, a structureheight of 3.20 m, a reef geometry composed by a delta without platform, a reefsubmergence of 1.50 m and a structure seaward slope of 1:10.The calibration of the numerical model COULWAVE has shown that it can be usedvery well to simulate the hydrodynamics around an MFAR. However, improvementsare recommended on the actual version of COULWAVE. These improvements includemainly the introduction of bottom porous structures. Even though the mean velocityfields give an indication about the sediment transport, they are not sufficient to get a deep insight about the influence of an MFAR on the morphology around the structure.For a proper design of the reef, particularly for the protection of a local coastline, amorphological study has to be executed. / Portugal é um dos muitos países no mundo que sofre os efeitos adversos da erosãocosteira. É sabido que as formas convencionais de protecção costeira apresentamalgumas desvantagens. Uma forma particularmente interessante de proteger oucontribuir para a protecção local de um troço da zona costeira, minimizando algumasdas desvantagens das estruturas tradicionais, é através da construção de recifesartificiais multi-funcionais. Um recife artificial multi-funcional é um quebra-marsubmerso construído com objectivos de ajudar a proteger localmente a zona costeira eoferecer outras valências, nomeadamente, aumentar as possibilidades locais para aprática de surf e contribuir para a valorização ambiental da região.Um recife artificial com estas características apresenta diversos aspectos positivos. Comefeito, esta estrutura: (1) não tem qualquer impacto visual negativo; (2) oferece osbenefícios económicos resultantes da atracção turística, ao melhorar as condições locaispara práticas de surf e/ou outros desportos náuticos (mergulho e pesca); (3) permiteaumentar o valor ambiental da região, e (4) com um projecto adequado permite reduzirao mínimo os efeitos adversos da erosão local ou prolongar mesmo a praia por efeito detômbolo.Neste trabalho são investigados os parâmetros característicos de um recife artificialmulti-funcional, tendo como caso de estudo um troço da zona costeira da Leirosa,situado a cerca de 15 km a sul da Figueira da Foz, Portugal. Para obter as característicasóptimas da estrutura a construir foram realizados vários estudos ao longo de diversasetapas. Como resultado do estado da arte e de um estudo teórico, realizado passo apasso, é proposta uma geometria preliminar de recife artificial multi-funcional. Estageometria preliminar do recife é em seguida utilizada como projecto inicial de diversostestes em modelo numérico e em modelo físico para analisar a sua capacidade naprotecção da praia e do cordão dunar da Leirosa e, em simultâneo, analisar aspossibilidades locais para a prática de surf.Além dos testes realizados em modelo físico, foi efectuado um estudo numéricoutilizando o modelo 2DV COBRAS-UC para obter informações sobre a influência docomprimento do recife e da submersão do mesmo no número Iribarren, sendo este umparâmetro que fornece uma boa indicação do tipo de rebentação a esperar. Neste estudofoi igualmente investigada a altura ideal para o recife. Posteriormente, procedeu-se àvalidação de um modelo numérico 2DH do tipo Boussinesq, o modelo COULWAVE,tendo como objectivo analisar a influência de três parâmetros de projecto com particulardesenvolvimento no plano horizontal: o ângulo do recife, a geometria do recife (com ousem plataforma) e as suas dimensões horizontais. Uma vez validado, foi este modelo numérico utilizado para analisar a influência daqueles parâmetros e identificar os seusvalores mais adequados para a definição de uma geometria optimizada do recife.Em termos de características para a prática de surf, e para as condições locais do litoralda Leirosa, foram identificados os seguintes valores para os principais parâmetros: umaestrutura de recife com ângulo de 66º, altura de 3.20 m e uma geometria composta porum delta sem plataforma, uma submersão do recife de 1.50 m e um declive de barlamarda estrutura de 1:10.A calibração do modelo numérico COULWAVE mostrou que este pode ser usado comsucesso para simular a hidrodinâmica em torno de um recife artificial multi-funcional.No entanto, recomendam-se melhorias na actual versão deste modelo. Tais melhoriasincluem principalmente a possibilidade de modelação com fundos porosos. Embora oscampos da velocidade média forneçam uma primeira indicação sobre o transporte desedimentos, esta informação não é por si só suficientemente profunda para obterconclusões definitivas quanto à influência de uma estrutura com estas característicassobre a morfologia local. Para uma concepção mais adequada da estrutura final dorecife, em particular para efeitos de protecção costeira, deverá ser elaborado um estudomorfodinâmico.
Ano: 2009
Autor(es): Voorde, M.
Info
Transformações Não-Lineares De Ondas Marítimas Em Zonas Portuárias. Análise pelo Método dos Elementos Finitos. Teses e programas de investigação do LNEC
N/A
Ano: 2003
Número Páginas:
400 p.
Autor(es): Fortes, C. J. E. M.
Editor: ******
Info
Modelação Estocástica Numérica e Física da Agitação Marítima
N/A
Ano: 2002
Número Páginas:
434 p.
Autor(es): Capitão, R.
Editor: ******
Info
Transformações Não-Lineares De Ondas Marítimas Em Zonas Portuárias. Análise pelo Método dos Elementos Finitos
N/A
Ano: 2002
Número Páginas:
400 p.
Autor(es): Fortes, C. J. E. M.
Editor: ******
Info
Modelação Estocástica Numérica e Física da Agitação Marítima
N/A
Ano: 2001
Número Páginas:
434 p.
Autor(es): Capitão, R.
Editor: ******
Info
Transformações Não-Lineares De Ondas Marítimas Em Zonas Portuárias. Análise pelo Método dos Elementos Finitos
N/A
Ano: 2001
Número Páginas:
400 p.
Autor(es): Fortes, C. J. E. M.
Editor: ******
Info
Interacción de ondas y grupos de ondas con estructuras de protección de costas
N/A
Ano: 1999
Número Páginas:
257 p.
Autor(es): Neves, M. G.
Editor: ******
|
|
Relatórios
Info
Estudo da estabilidade de molhes de talude.
N/A
Ano: 1962
Autor(es): Castanho, J.
Info
Estudo sobre o prolongamento do molhe de Ponta Delgada.
N/A
Ano: 1962
Autor(es): Abecasis, F.; Vera-Cruz, D.
Info
Estudo sobre o prolongamento do molhe de Ponta Delgada. Estabilidade de caixões durante a construção
N/A
Ano: 1962
Número Páginas:
4.
Autor(es): Abecasis, F.; Vera-Cruz, D.
Info
Licitacion publica nacional e internacional n.º 3611. Estudio del puerto de Comodoro Rivadavia. Plan de trabajos
N/A
Ano: 1962
Número Páginas:
7.
Autor(es): ROCHA, M.
Info
Propagação e rebentação de ondas. Concurso para obtenção do certificado de estágio
N/A
Ano: 1962
Número Páginas:
65 pp.
Autor(es): Vera-Cruz, D.
Info
Assoreamentos no porto da Póvoa de Varzim. Estudo das condições de agitação no porto. 2.º Relatório.
N/A
Ano: 1961
Autor(es): Abecasis, F.; Vera-Cruz, D.
Info
Ensaio em modelo reduzido do porto de Angra do Heroísmo.
N/A
Ano: 1961
Autor(es): Abecasis, F.; Reis de Carvalho, J.
Info
Ensaio em modelo reduzido do porto de Angra do Heroísmo.
N/A
Ano: 1961
Autor(es): Reis de Carvalho, J.
Info
Ensaio em modelo reduzido do porto novo de Santo Antão.
N/A
Ano: 1961
Autor(es): Vera-Cruz, D.
Info
Ensaio em modelo reduzido do Porto Novo de Santo Antão. Ensaios de agitação (II Parte). Relatório
N/A
Ano: 1961
Autor(es): Abecasis, F.
|
|
|
Outros
Info
Risk Perception of Coastal Flooding in Praia da Vitória, Azores
The project LIFEGARACHICO (LIFE20 CCA/ES/001641) proposes the creation of an effective flexible adaptation framework for the coastal unicipalities of Macaronesia, making specific local interventions to increase their resilience against extreme coastal events resulting from climate change. Social surveys were applied in Praia da Vitória, Terceira, Azores, to understand the civilians perception of coastal flooding events.
Ano: 2024
Autor(es): Lima, L.; Fortes, C. J. E. M.; Zózimo, A. C.; Pinheiro, L.
Info
Social Surveys Indicating the Resilience of Urban Coastal Localities in Macaronesia to Coastal Flooding
The project LIFE-GARACHICO (LIFE20 CCA/ES/001641) proposes the creation of an effective flexible adaptation framework for the coastal municipalities of Macaronesia, making specific local interventions to increase their resilience against extreme coastal events resulting from climate change. Social surveys were applied in both localities in order to understand the civilians perception of coastal flooding events
Ano: 2024
Autor(es): Fortes, C. J. E. M.; Lima, L.; Zózimo, A. C.; Pinheiro, L.; Modino, R.; Barroso, J. G.; Cruz, C. G.; Comes, L.
Info
Avaliação da Estabilidade do Quebra-Mar Destacado do Porto Inglês (Ilha do Maio
Realização de ensaios em modelo físico reduzido 3D para a avaliação da estabilidade do quebra-mar destacado do Porto Inlgês, Maio, Cabo-Verde
Ano: 2023
Autor(es): Soares, R. E.; Fortes, C. J. E. M.; Lemos, R.; Costa, C.E.B.C.; Silva, L. G.; Pinheiro, L.; Neves, M. G.
Info
Physical modelling of wave-induced plant drag coefficient
N/A
Ano: 2019
Autor(es): Reis, R.; Suzuki , T.; LUBBERS, B.; Fortes, C. J. E. M.; Pires Silva, A.
Keywords: Flume experiments; Wave propagation; Plant flexibility; Vegetation
Info
Evaluation of wave overtopping consequences for risk assessment using ANP methodology
Adverse sea conditions can cause emergency situations associated to wave overtopping, which endanger the safety of people and goods, with negative impacts for society, the economy and the environment. Therefore, a methodology to assess the overtopping risk in port and coastal areas is essential for a proper planning and management of these areas.The Portuguese National Laboratory for Civil Engineering (LNEC) has been developing the HIDRALERTA system, an integrated decision-support tool for port and coastal management, which focus in preventing and supporting the management of emergency situations and the long-term planning of interventions in the study areas. It enables the user to calculate the risk for various port and coastal activities, starting with the characterization of sea-waves, wind field and tide levels at the study regions.HIDRALERTA was already tested at different places in Portugal: Costa da Caparica beach, in Almada, and port of Praia da Vitória, in the Azores. However, in both cases, the evaluation of wave overtopping / flooding consequences was carried out using a quite simple approach. In fact, such consequences were assessed only at a global and qualitative level, based on information provided by the responsible authorities.This paper aims at testing a quantitative methodology to evaluate overtopping consequences: the Analytic Network Process (ANP). This methodology was proposed to solve complex decision-making problems and its objective is to calculate the priorities among decision elements which will define their final weights on the consequences. The ANP approach can produce interdependencies between criteria and compute the respective weight of each criterion. Such a network model with dependence and feedback improves the priorities derived from judgements and makes prediction more accurate. So, the ANP allows comparison among clusters of elements. Furthermore, in this methodology we take into account the impacts of alternatives on the importance of criteria, by normalizing the comparison matrix. Thus, the given alternatives can influence the ranking of criteria.Application of the ANP methodology to the case study of Praia da Vitória port is presented and a consequences level map is obtained.
Ano: 2016
Autor(es): Lourenço, I.; Santos, J. A.; Fortes, C. J. E. M.; Reis, M. T. L. G. V.; Poseiro, P.; Craveiro, J.
Keywords: Praia da Vitória port; Analytic Network Process (ANP); Risk; Wave overtopping
Info
Sistema de alerta para a navegação em portos
Apresenta-se um novo sistema de alerta para navios amarrados em portos que se encontra em desenvolvimento no Laboratório Nacional de Engenharia Civil. O sistema SWAMS_ALERTA é um sistema de previsão e alerta baseado na avaliação do risco associado ao comportamento de navios amarrados em zonas portuárias.Os riscos associados a navios amarrados decorrem de movimentos excessivos, quer no plano horizontal, quer no plano vertical. Estes movimentos podem condicionar as atividades portuárias, nomeadamente cargas e descargas, mas em casos extremos podem levar à ocorrência de situações de emergência, tais como rotura de cabos, rotura de cabeços de amarração ou até colisões com o cais. As consequências deste tipo de situações envolvem sempre grandes prejuízos materiais e por vezes humanos.O sistema SWAMS_ALERTA utiliza as medições e previsões de agitação marítima para a determinação dos seus efeitos em termos de movimentos nos seis graus de liberdade e forças nas amarras e defensas quando o navio se encontra estacionado no cais. Para tal, o sistema recorre a uma série de modelos numéricos, interligados entre si. A comparação destes valores com valores máximos admissíveis pré-estabelecidos permite a avaliação, em tempo real, de situações de emergência e a emissão de alertas dirigidos às entidades portuárias.O sistema é constituído por 4 módulos:I
Ano: 2016
Autor(es): Pinheiro, L.; Lopes, P.; Fortes, C. J. E. M.; Santos, J. A.
Keywords: Navios amarrados; Riscos; Sistema de previsão e alerta
CONSTRUÇÃO DA EXPANSÃO DO PORTO DE SAL-REI 1ª FASE (ILHA DA BOAVISTA CABO VERDE) Ensaios bidimensionais em modelo reduzido Parecer Técnico
Info
Definição da Geometria para um Recife Artificial Multifuncional em Leirosa (Poster)
N/A
Ano: 2010
Número Páginas:
1.
Autor(es): Neves, M. G.; Antunes do Carmo, J.
Editor: Gestão Costeira: Vulnerabilidades e Riscos na Região Centro
Info
Probabilistic Assessment of Hazards from Overtopping of Seawalls in Urban Environments: Annual Report of pOST-dOCTORAL rESEARCH
N/A
Ano: 2005
Autor(es): Reis, M. T. L. G. V.
Editor: ******
Info
Guidelines and Recommendations for River Information Services - PIANC
N/A
Ano: 2004
Autor(es): J.A Covas
Editor: ******
|