Publications

The generation of tidal asymmetries in estuaries is analyzed by examining the differences in celerity at high and low tide. An analytical expression for the difference between ebb and flood durations is introduced and applied to a cross-section composed by a channel flanked by tidal flats. The expression depends on the cross-sectional average of the celerity, for which two approximations are proposed and verified with a numerical model. Results confirm that tidal flats enhance ebb dominance, while large tidal amplitudes promote flood dominance. Maximum ebb dominance occurs for tidal flats at, or above, mean water level, depending on tidal amplitude and the extent of the tidal flats. Finally, qualitative arguments indicate that friction also enhances flood dominance.

N/A

A modeling system to simulate tidally driven morphodynamics is developed and assessed. The system couples an existing shallow water model with a new morphological updating model. Four innovative features to improve the accuracy and efficiency of the system are implemented and tested. First, two alternatives are proposed to update the velocity field without running the hydrodynamic model. Both alternatives, based on an assumption on the behavior of friction, are more accurate than the traditional "continuity correction". Second, hydrodynamic results are provided in the frequency domain. The harmonic synthesis of the velocity time series simplifies interpolation and extrapolation of the hydrodynamic model results and the use of different time steps for the various components of the modeling system. Third, the time-integration of the sediment fluxes is performed with a Runge-Kutta method so the time step adapts in space and time to the flow characteristics. Finally, a criterion to determine the need for a new hydrodynamic simulation is derived, based on the importance of the errors introduced by the outdated flow field relative to those introduced by the sand-transport formulae.

N/A

N/A

N/A

N/A

N/A

Descreve-se um modelo de transporte para sistemas costeiros, particularmente apropriado em situações de fortes gradientes de concentração. Este modelo baseia-se numa formulação Euleriana-Lagrangiana para malhas não-estruturadas e apresenta as seguintes características: 1) discretização por volumes finitos e tracking de elevada precisão; 2) malhas distintas para o escoamento e transporte; 3) use da técnica de flux-corrected transporte ou de filtro não-linear para remoção de oscilações. Analisam-se as propriedades numéricas através da análise do erro de truncatura e de um teste numérico. A combinação de várias técnicas numéricas permite obter estabilidade incondicional, conservação de massa e elevada precisão, evitando simultaneamente oscilações numéricas.Palavras-Chave: Volumes finitos; métodos Eulerianos-Lagrangianos; oscilações numéricas.

The Barbarie spit (Langue de Barbarie) is a 100

ShorelineS is a recent coastline evolution model allowing for extreme coastline curvatures and spit formation at high-angle wave incidence. This paper investigates ShorelineS

The XBeach modelling system was used to hindcast the morphological evolution of an ephemeral tidal inlet during Hurricane Leslie (2018). Simulations of morphological evolution were compared against daily observations from topographic surveys conducted at the Albufeira Lagoon during this event. On the one hand, numerical simulations successfully reproduced the observed accumulation pattern at the flood delta. On the other hand, numerical simulations displayed a slight tendency to overestimate the sediment erosion that took place at the transitional channel.

During storm conditions, infragravity (IG) waves with periods between 25 s and 300 s can dominate the seasurfaceoscillations near the shoreline. Through processes like wave run-up, overwash and wave overtopping,IG waves can hence become major drivers of coastal hazards. Although IG dynamics are relatively wellunderstood in open beaches, their contribution of the hydro-sedimentary dynamics of more complex coastalsettings, such as tidal inlets characterized by their ebb and flood deltas connected through a transitional channel, is not well understood. Here, we assess the performance of the widely used XBeach numerical model to predict IG wave dynamics through a comparison against field observations collected on an ephemeral tidal inlet (Albufeira Lagoon) during storm Leslie (2018). XBeach accurately predicted the observed IG significant wave heights reaching up to 1 m over a 2-day period. This performance is assessed at both the nearby beaches and the flood-delta.

The inlet of the Óbidos Lagoon is a wave-dominated tidal inlet located on the Portuguese coast. It is known for its unstable behaviour, which often leads to erosion issues on its margin or to its complete closure. Over twenty years of satellite data were compiled to document this behaviour. The reconstructed behaviour was compared to the contemporary offshore wave climate. The comparison showed that the inlet

The geometric characteristics of sand waves in the main channel of the outer Tagus estuary are investigated, using high-resolution surveys conducted in 2019 and 2020. The relationship between bedform height and spacing over a 1400 m transect, at ca. 20 m depth, is in good agreement with the global mean trend equation proposed by Flemming (1988). The consistent ebb-oriented asymmetry suggests that the direction of net sediment transport is seaward. The apparent direction of sand wave migration obtained from the comparison between surveys is more variable, but on average also ebb-oriented.

O aprofundamento dos canais de acesso por dragagem, associado ou não à construção de molhes e quebra-mares, tem sido uma característica da evolução dos portos nacionais nos últimos 100 anos. Nesta comunicação, analisa-se a evolução da cota de dragagem dos acessos marítimos aos portos do Continente com base em fontes documentais e levantamentos hidrográficos. No início do século XX verificava-se uma diferença apreciável entre as condições oferecidas pelos portos situados em estuários de maiores dimensões (Lisboa e Setúbal) e artificiais (Leixões), por um lado, e os portos situados em estuários e sistemas lagunares com menor prisma de maré, por outro. Entre 1920 e a atualidade, principalmente desde 1950, todos os portos comerciais em atividade (excetuando Sines) registaram um aumento da profundidade dos canais de acesso entre cerca de 5,5 m e 11,2 m. Considerando os portos no seu conjunto, estima-se que em média, entre 1950 e 2020, se tenha verificado um aumento da profundidade de dragagem na zona da barra de 1,0 m por década.

As inundações costeiras afetam anualmente milhões de pessoas em todo o mundo, com grandes perdas socioeconómicas que irão aumentar face à subida do nível médio do mar e ao aumento dos níveis extremos de água. A extensa linha de costa continental portuguesa, caracterizada por grande diversidade natural e de ocupação, apresenta elevada exposição a fenómenos de galgamento e inundação costeira. Estima-se que o impacto da tempestade Hércules, no início do ano 2014, tenha sido cerca de 16 milhões de euros (Santos et al., 2015). As consequências do impacto de níveis extremos na zona costeira dependem das características do território, pelo que a capacidade de obter e integrar esse conhecimento em ferramentas de previsão que suportem adequadamente as necessidades da gestão do risco deve ser uma prioridade. No âmbito do projeto Mosaic.pt (http://mosaic.lnec.pt) desenvolveram-se ferramentas de avaliação do risco de inundação costeira, integrando a dimensão territorial. A complexidade e a diversidade da costa continental recomendam uma abordagem holística na identificação e diferenciação das características territoriais, que dependem da escala espacial considerada. Assim, a abordagem metodológica do projeto seguiu duas escalas espaciais de análise: a costa continental portuguesa e três observatórios representativos do trecho costeiro Ovar - Marinha Grande, onde se concentra o maior número de ocorrências de inundação costeira (Tavares et al., 2021). Estes observatórios são o setor Barra

The continental coastal zone of Portugal is characterized by its morpho-sedimentary diversity. The high population density and range of land uses, occupation and activities make it an area of great national strategic value, and a multi-hazard zone. In this study, a multidimensional methodology to analyse, evaluate and interpret local vulnerability (Coastal Territorial Vulnerability Index - CTVI) has been developed and applied in three study areas.

Continental Portugal presents an extensive and diversified coastal zone. This diversity gives rise to the existence of different coastal typologies such as sand beaches, rocky, cliffed and artificialized coasts, which support distinct levels of human occupancy. To identify and distinguish the different coastal typologies, cluster analysis was used based on a set of variables, expressing the territorial complexity, as well as the occurrences and impacts resulting from coastal flooding.

N/A

The hydrodynamics and sediment dynamics of a low energy estuarine beach was analysed. The beach is subjected to local wind waves generated in an area of restricted fetch, wake waves generated by catamarans and a semi diurnal meso-tidal regime. Only storm events of short duration modify the beach profile, which, once the normal hydrodynamic conditions are restored, naturally recovers its initial shape. The beach is characterised by having a steep upper slope until mean sea level, followed, seaward, by a low gradient terrace, and a bimodal sediment distribution. To investigate the hydro sedimentologic beach behaviour, wind, topo-hydrographic and sedimentologic data were used as input and validation of process based numerical modelling. Wave generation and transformation in an area of restricted fetch, nearshore circulation, sediment transport and morphological evolution were simulated for average annual conditions and storm conditions. The statistical analysis of a six-year wind data series allowed to derive the average wind regime, based on which, the average annual wave regime and the average annual longshore sediment budget, 14.5x103 m3year-1, were calculated. The contribution of the individual components of the representative wave regime, discretised by directional sector and height class of incidence, was evaluated as well as the spatial distribution of the longshore sediment transport in the active part of the beach profile. The characteristics of the wave groups generated by the passage of catamarans at different speeds were estimated and their action on the beach morphology was simulated. Although when the catamarans travel at 20 knots speed the average annual wake wave energy dissipated at the beach is 2.5 times higher than the average annual wind wave energy, the erosion effect on the beach profile is still not relevant for the present traffic. In opposition, short duration storm events generate the formation of an erosion scarp at the upper part of the beach face. The numerical modelling of this phenomenon allowed to acknowledge on the protective effect that the low gradient terrace has on the beach face.

O objectivo do presente trabalho foi o estabelecimento de um modelo matemático do fenómeno da refracção das ondas de gravidade, apoiado num pro¬grama cálculo automático bastante aperfeiçoado e enquadrado por uma técnica de utilização prática que permitisse a resolução rápida dos problemasinerentes ao cálculo de um conjunto de diagramas de refracção, independentemente do número de ondas monocromáticas a analisar.Numa primeira parte, apresentam-se algumas considerações sumárias so¬bre o método de estabelecimento do programa de cálculo e sobre a necessidade de se dispor de um modelo matemático de refracção. A segunda parte é dedicada ao estudo do fenómeno da refracção das on¬das de gravidade e seu tratamento numérico. No Cap.3 indicam-se, de forma Muito sintética, os fundamentos da teoria das ondas de pequena amplitude sendo deduzidas as principais expressões utilizadas no estudo da refracção. É apresentada a teoria da refracção de uma onda sinusoidal, sendo abordado um vasto leque de problemas tais como os aspectos energético e geométrico da refracção, as soluções analíticas para fundos algébricos, as zonas de convergência e cáusticas, a refracção por efeito de correntes, e as aproximações da teoria da refracção. Por último criticam-se os métodos gráficos habituais de cálculo de planos de ondulação. No Cap. 4 são deduzidas as equações de base do tratamento numérico da refracção das ondas sinusoidais como um problema de valores iniciais. No Cap. 5 são apresentadas as dife¬rentes técnicas de cálculo e de interpolação das principais grandezas di¬rectamente ligadas a obtenção de uma ortogonal e dos sucessivos valores da amplitude relativa da onda. São indicados os métodos numéricos de integração das equações diferenciais do raio de onda e do coeficiente de afasta¬mento de ortogonais. E feita referencia a técnicas de alisamento da batimetria. Em todos os casos são apresentados os métodos utilizados no LNEC as¬sim como outros referidos na bibliografia. No Cap. 6 são indicadas, para os vários programas preparados no LNEC, as técnicas de cálculo e de inter¬polação e os métodos numéricos de integração usados. São apresentadas, de uma forma exaustiva, todas as verificações dos programas levadas a cabo após o seu estabelecimento, designadamente a verificação analítica e a comparação dos resultados obtidos com os métodos numéricos e gráficos para fundos ir¬regulares. No âmbito da verificação dos programas é ainda feita referência a resultados obtidos por outros autores e a ensaios em modelo reduzido. Por último tiram-se conclusões sobre a escolha dos programas a lançar em exploração.A terceira e última parte é consagrada à análise do actual programa de cálculo de diagramas de refracção. Assim, no Cap. 7 indicam-se, em li¬nhas muito gerais, o tipo de dados a fornecer, as hipóteses de exploração do programa e o seu modo de funcionamento; é ainda feita uma análise da estrutura do programa, relativamente pormenorizada. 0 Cap. 8 é consagrado à discrição da forma de utilização prática do programa, sendo discutidas as perspectivas de custo da sua utilização cobrindo todas as hipóteses de exploração do programa, ilustrados com os dados e resultados respectivos. No Cap. 9 apresentam-se as conclusões do trabalho e indicam-se as linhas de orientação a seguir no que respeita ao aperfeiçoamento do modelo matemático da refracção e ao estabelecimento de outros modelos matemáticos dependentes ou não daquele.0 texto principal de dissertação encerra-se com 96 referencias bibliográficas, correspondentes aos principais trabalhos consultados pelo autor, na realização do presente estudo; esta bibliografia refere-se a Abril de 1975. Em Anexo no Apêndice I, apresentam-se os elementos referentes aos programas inicialmente estabelecidos no LNEC e cujo aperfeiçoamento posterior conduziu ao programa actualmente em exploração. São indicados, em linhas gerais, as formas possíveis de exploração daqueles programas, o seu modo de funcionamento e os dados a fornecer. É feita uma análise de pormenor da es¬trutura dos programas, sendo descrita a sua forma de utilização prática, e apresentados vários exemplos de aplicação ilustrados com os dados e resul¬tados respectivos. São indicadas algumas verificações efectuadas a estes programas.

0 caminhamento das areias ao longo das costas marítimas é regulado por leis ainda incompletamente apreendidas, mas cujo conhecimento cada vez mais perfeito se torna indispensável no projecto e ensaio em modelo reduzido de grande parte das obras destinadas a solucionar problemas de hidráulica marítima.Começa-se neste trabalho por caracterizar, sob os pontos de vista hidrodinâmico e de movimentação aluvionar, a orla costeira, zona em que o transporte litoral se processa, abordando-se a seguir o aspecto prático mais importante do fenómeno, ou seja, a previsão dos caudais sólidos a partir das características das ondas e do material móvel que forma a praia.Com este último objectivo analisa-se o conjunto dos dados experimentais disponíveis, segundo uma perspectiva baseada na análise dimensional do fenómeno, para a seguir se proceder a um estudo experimental focado sobre a influência do peso específico e do diâmetro do material móvel.Apresentam-se ainda, em anexo, os conceitos básicos respeitantes às correntes longitudinais, fenómeno que teve também de ser objecto de um estudo aprofundado, por estar intimamente relacionado com o transporte sólido litoral.

Este trabalho constitui um estudo em extensão sobre as característi¬cas mais importantes dos perfis transversais das praias e sobre a evolução resultante da actuação da ondulação. Pretendeu-se assim obter um conhecimento geral sobre este assunto, de modo a permitir, por um lado, o melhor planeamento dos trabalhos marítimos e, por outro, uma melhor interpretação dos modelos reduzidos de fundo móvel.Utiliza-se neste trabalho uma nomenclatura original em língua portuguesa, na tentativa de estabelecer uma equivalência com os termos consagrados na literatura siglo-saxónica.Dentro do esquema geral acima enunciado, são estudadas com maior atenção: o mecanismo da movimentação dos materiais nos fundos; a evolução so¬frida pelos perfis em resultado da actuação das ondas e destas em conjugação com as marés; as características da faixa de rebentação das praias; e as variações da granulometria do material móvel.Quanto à movimentação dos materiais em fundos de maior profundidade, apresenta-se uma interpretação sobre as profundidades a que se inicia o transporte dos sedimentos e o sentido desse transporte.Em relação à evolução dos perfis descrevem-se as suas fases mais significativas, quer para o caso de ondas de período e altura crescentes como decrescentes. Introduzem-se os conceitos de perfis de máxima erosão e máximo enchimento.

O presente trabalho baseia-se num estudo experimental sobre o funcionamento hidráulico de esporões marítimos no caso em que estes limitam troços de praia que se comportam, na sua evolução, como unidades fisiográficamente independentes. Pretende-se, com este estudo, contribuir para o conhecimento do mecanismo de evolução de praias em que são construídos campos de esporões, com vista ao projecto de dispositivos funcionais deste tipo de defesa costeira. Acessoriamente, tirando partido dos esporões como meio de medição do caudal sólido litoral, fa¬zem-se algumas considerações sobre as relações existentes entre o transporte litoral e as características da onda actuante.O relatório é constituído pelas seguintes partes : a Introdução, em que se define o objectivo do estudo e as suas características; o capítulo I, em que se apresentam noções gerais sobre esporões; o Capi¬tulo II, em que se dá conhecimento do estudo experimental realizado, definindo o seu planeamento, a tecnologia experimental, as características dos ensaios realizados, a interpretarão dos resultados obtidos e as conclusões; o Capitulo III, que diz respeito a casos reais de emprego de esporões, em que, após umas generalidades, se apresentam as principais características das obras de defesa da Cova do Vapor - Costa da Caparica, do campo de esporões da praia de Espinho e do esporão da praia do Estoril, finalizando com algumas conclusões; o Capítulo IV em que se apresentam considerações sumárias sobre caudal sólido litoral, as quais compreendem algumas generalidades, considerações teóricas, resultados experimen¬tais, observações na natureza e conclusões; finalmente, as Conclusões Finais, onde se agrupam todas as considerações do estudo que se julgam de primordial interesse.

The final technical report of the project

O presente estudo tem como objetivo caracterizar o transporte sedimentar longitudinal potencial (capacidade de transporte) entre as embocaduras dos rios Mondego e Lis no período 1952-2010. Enquadra-se no âmbito de um estudo mais alargado, cujo objetivo é o estudo da dinâmica sedimentar costeira e a previsão da evolução futura da linha de costa do mesmo trecho.O estudo baseia-se em modelação numérica dos processos físicos costeiros para determinar a série temporal (com frequência de 6 horas) de transporte sedimentar longitudinal potencial induzido pela agitação marítima e variação do nível do mar devida à maré astronómica e meteorológica, no período 1952-2010. Concluiu-se sobre a elevada variabilidade interanual e sazonal do transporte sedimentar longitudinal potencial, Qs, e confirmou-se a direção dominante da deriva para sul (em 46 dos 59 anos em estudo). Obtiveram-se volumes anuais médios de Qs Total e Resultante para sul de 918 e 270 x103 m3 e volumes anuais extremos de Qs Total 578 e 1 500 x103 m3. Da análise do diagrama de balanço sedimentar do ano de Qs Total médio e da distribuição de Qs ao longo do perfil transversal, concluiu se que grande parte do volume mobilizado se deve aos eventos com Hrms 2.00-5.00 m e Dir 295-305ºN e identificaram-se as zonas onde essa mobilização é mais intensa. Da análise à sazonalidade concluiu-se que no verão marítimo apenas existe potencial para mobilizar cerca de 60% do volume sedimentar mobilizado no inverno, sendo evidente nas duas estações marítimas a deriva resultante na direção sul. Da análise adicional aos anos de Qs Sul e Qs Norte máximos identificaram-se extensões da zona ativa do perfil de praia de 1 500, 2 270 e 1 970 m nos três anos referidos, respetivamente.

Apresenta-se uma metodologia de estudo de desequilíbrios morfológicos costeiros de grande envergadura. Esta metodologia, baseada na modelação numérica, foi utilizada para: analisar as evoluções sofridas pela linha de costa em toda a extensão espacial e temporal, do Douro ao Cabo Mondego, desde 1950 até ao presente; e antever a evolução nas próximas décadas face a intervenções alternativas de mitigação e controlo do processo erosivo. Usaram-se os dados disponíveis, que, embora apresentando algumas falhas, permitiram confirmar a viabilidade e potencial da metodologia. Em síntese, a metodologia que se propõe baseia-se numa interpretação global dos fenómenos em curso, quantificada através de balanços sedimentares e da aplicação de um modelo numérico de evolução de linha de costa. Recomenda-se a análise de possíveis intervenções considerando um comprimento de costa e um período de tempo adequados aos processos de evolução morfológica em curso.

Apresentam-se e discutem-se os resultados do Programa de Monitorização na fase pré-obra, referente ao estudo Monitorização da Hidrodinâmica da Abertura e Aprofundamento dos Canais da Zona Inferior da Lagoa de Óbidos, com vista à caracterização da situação inicial. Os dados adquiridos incluem um levantamento topo-hidrográfico e um ortofotomapa, dados de agitação marítima, níveis do mar na zona da Lagoa e níveis no seu interior. Os resultados evidenciam forte assoreamento da zona inferior da Lagoa, traduzido num claro predomínio das áreas situadas acima do Zero Hidrográfico. O canal de embocadura, situado no terço central da Lagoa, reflete a intervenção recente com vista à sua abertura, apresentando-se pouco desenvolvido. Este assoreamento traduz-se em amplitudes de maré muito reduzidas no interior da laguna (cerca de 10% da amplitude no mar), em que as principais variações de níveis ocorrem numa escala de tempo quinzenal. Simultaneamente, a Lagoa de Óbidos apresenta-se muito dominada pela enchente.

This report presents the main activities and results of the project RealQual: Towards real-time highresolutionmonitoring and prediction of water quality in estuaries and coastal areas, a collaborative project of the FLAD / NSF 2013, Portugal - U.S. Research Networks Program. RealQual aimed at developing the foundations of a quality-controlled network for real-time, high-resolution monitoring and prediction of the water quality in estuaries and coastal areas. Regarding the real-time monitoring component, an online monitoring station for water quality was maintained operational in the Tagus estuary. Several procedures were tested for data quality control using the data from the SATURN OHSU-CMOP (Oregon Health & Science University

A Ria Formosa (sul de Portugal) é uma laguna muito produtiva, de elevada importância ecológica eeconómica. Assim, perceber os processos que afetam a produtividade deste ecossistema éfundamental para a sua gestão sustentável. O projeto COALA visa estudar a dinâmica de trocas denutrientes, clorofila e material particulado entre a Ria Formosa e o oceano adjacente, assim como ainfluência dos principais forçamentos (e.g. maré) nesta dinâmica, utilizando uma abordagemmultidisciplinar que combina modelação numérica e campanhas de campo. O modelo hidrodinâmicoSELFE foi aplicado na Ria Formosa, representando adequadamente os principais padrões decirculação. Este modelo foi combinado com um modelo de transporte lagrangeano para estudar ainfluência dos forçamentos ambientais na dinâmica de trocas lagoa

N/A

A granulometria das partículas em suspensão foi analisada por difracção laser em ensaiosem canal anular com sedimento do Esteiro dos Frades, Ria de Aveiro. Para velocidadesdecrescentes, verifica-se uma diminuição do diâmetro mediano (D50), a par de umadiminuição da concentração de sedimento em suspensão. Para velocidades crescentes,verifica-se, pelo contrário, um aumento de D50, a par de um aumento da concentração.Conclui-se que, embora o sedimento seja predominantemente silto-argiloso, o efeito dagranulometria nos processos de transporte é relevante.

A descrição topo-hidrográfica dos sistemas lagunares constitui informação fundamental para acompreensão da dinâmica hidrológica, morfológica e sedimentar aí operante, bem como paraa respetiva modelação, numérica ou física. No queconhecidos levantamentos topoQuintino (1988) apresente uma figura de localização das estações de amostragem que incluialgumas linhas de contorno batimétrico (profundidade, estimadas em relação ao plano de água existente no momento de realizaçãodas campanhas de campo.Figura 1. Esboço da batimetria da Lagoa de Albufeira (adaptado de Quintino, 1988,Em comparação com estruturas similares do litoral nacional, a laguna de Albufeira éinvulgarmente profunda: a profundidade máximaLagoa Grande e 2 m na zona central da Lagoa Pequena, tendo o canal de ligação entre ambasmenos de 2 m. No entanto, a espessura da coluna de águapelo que, com a barra fechada e em anos de pluviosidade elevada,excepcionalmente, como no inverno de 1990Grande e 6 m na Lagoa Pequena.bem como na zona terminal (interna) dos depósitos interiores,lagunar.zOcano Atlântico3 1013BocaVelhaBico D.FernandoRibªCRIAÇÃO E IMPLEMENTAÇÃO DE UM SISTEMA DE MONITORIZAÇÃO NO LITORAL ABRANGIDO PELADE JURISDIÇÃO DA ADMINISTRAÇÃO DA REGIÃO HIDROGRÁFICA DOrespeita à Lagoa de Albufeira, não sãotopo-batimétricos da totalidade do espaço lagunar, emboraFigura 1), obtidas por interpolação de sondas dein Freitas e Ferreira,2004).em 1988 rondava 13 m no corpo oeste davaria consideravelmenteesta pode exceder (embora

O presente documento constitui o relatório final do projeto MorFeed, financiado pela Fundação para a Ciência e a Tecnologia. Listam-se abaixo os principais elementos identificadores do projeto:Referência do Projeto: PTDC/AAC-AMB/100092/2008Área Científica: Ambiente e Alterações ClimáticasInvestigadora Responsável: Paula Maria de Santos FreireInstituição Proponente: Laboratório Nacional de Engenharia Civil (LNEC)Parceiro: Fundação da Faculdade de Ciências da Universidade de Lisboa (FFCUL)Data de Início: 01-04-2010Data de Fim: 30-09-2013Financiamento Concedido:

N/A