Publications

Selected interpolation and integration finite element Eulerian-Lagrangian methods for the solution of the transport equation are systematically compared, through a combination of formal analysis and numerical experimentation. Methods that approach conditions at the feet of the characteristic lines with the perspective of "integration" rather than "interpolation" tend to have superior accuracy. However, integration methods that track quadrature points are only conditionally stable. Stability criteria are provided for these methods, and are only mildly restrictive for most applications. A method that performs an exact integration is proposed and shown to have the best accuracy among the studied methods. Key Words: Transport equation. Eulerian-Lagrangian methods. Stability analysis. Accuracy analysis. Numerical experimentation.

N/A

N/A

N/A

N/A

N/A

N/A

N/A

A baía da cidade da Praia, na ilha de Santiago, Cabo Verde, acolhe um dos principais portos comerciais e de cruzeiros de Cabo Verde, tendo também uma relação importante, do ponto de vista ambiental e paisagístico, com a frente ribeirinha da cidade. Procedeu-se a uma análise da evolução morfológica do sistema, com recurso a levantamentos hidrográficos (de 1882 a 2016) e imagens aéreas (de 1968 a 2018), para determinar tendências evolutivas históricas. Realizou-se também a aplicação de um sistema de modelos numéricos, envolvendo a simulação da hidrodinâmica, do transporte sedimentar e da evolução do fundo, para identificar os processos litorais dominantes. Calculou-se uma tendência média de avanço da praia Grande ou da Gamboa de 0,6 a 1,0 m/ano, correspondendo este último valor à média dos últimos 50 anos. As condições de maior hidrodinamismo no interior da baía estão associadas aos períodos de agitação marítima de maior altura. As intervenções mais recentes, nomeadamente na zona portuária, tendem a diminuir esse hidrodinamismo.

This study investigates the roleof infragravity waves (IGW) on the hydrodynamics of a wave-dominated inlet (Albufeira lagoon, Portugal) during storm Leslie (October 2018) through the analysis of field observations. Measured significant IGW heights reached up to 0.9 m in the surf zone, at adjacent ocean beaches, and 0.4 m at the flood-delta. At the flood-delta, IGW were more energetic than wind-generated short-waves by up to a factor of 2. These large IGW are expected to contribute significantly to sediment transport during storm conditions.

Existing knowledge on the evolution of the Tagus estuary inlet morphology derives mostlyfrom analyses of bathymetric surveys and aerial images. The limitations of this approachprevent a detailed understanding of the processes responsible for this evolution, as well as theprediction of the future evolution. A 2DH morhodynamic model was therefore applied to shednew light into the sediment dynamics of this tidal inlet. Description of the method

Due to their exposure to various forcings such as tides and waves, estuary mouths constitute dynamic coastal environments where sediment transport can result in abrupt and significant morphological changes. They concentrate relevant socio-economic and ecological values, which will become even more challenging in the climate change context. Despite their importance worldwide, the dynamics of these environments remains partly understood only and fosters many studies from the scientific community. This work is part of a citizen-science project (named Nouveaux Commanditaires Sciences) which aims at studying the hydro-sedimentary dynamics of the Tagus estuary mouth (bordered by a metropolitan area of 2 million inhabitants), especially during storms, by combining in-situ field observations, remote sensing data and numerical modelling. This research program has been set up following the request for a better understanding of the coastal erosion and marine flooding during storms expressed by severely affected inhabitants from an informal neighbourhood, Segundo Torrão (Southern margin of the inlet). In addition, local populations provide support for instrumental deployments and data acquisitions. In return, new knowledge derived from scientific studies will be used to improve coastal management over these threatened areas. In this presentation, we will present the main actions that are being carried out in the scope of this project.

To increase the response capability to coastal flooding events in an accelerating sea level rise context, the ability to predict floods has to improve and be integrated with other components of risk, such as vulnerability and exposure. The MOSAIC.pt project aims at developing an advanced flood risk management framework for coastal zones, integrating robust predictive tools and multi-source data with the different dimensions of the risk, to adequately contribute to emergency planning and response

Estuarine uses and activities may be negatively affected by climatic variability andclimate change

Os observatórios costeiros podem apoiar a gestão diária e a longo prazo dos sistemascosteiros, permitindo antecipar eventos de contaminação, contribuir para a vigilânciacontínua destas zonas e implementar e adaptar planos de gestão. Ao integrarem previsões,análises de cenários e indicadores, para além de observações históricas e em tempo real, estasferramentas têm aumentado de complexidade ao longo dos anos.

Coastal forecast systems have become reliable tools to support daily coastal management and strategic long-term planning. However, the availability of these systems for most coastal areas remains scarce, in particular for small estuaries and coastal bays, which limits the efficiency of management for these areas from both economic and environmental perspectives. Building a new forecast system or maintaining in operation an existing one requires considerable effort, from both financial and personnel points of view. Accurate forecasts require a team comprising experts in coastal processes and computer science and the availability of robust and large computational resources to maintain the predictions available on a daily basis. Herein, a new concept in forecast tools is presented, building the foundations for forecast systems to be built and made available for all coastal systems worldwide. Denoted OPENCoastS, this freely available web service builds forecasts systems for the coastal area selected by the user and maintains them in operation using the European Open Science Cloud infrastructure. Forecasts are based on the SCHISM community model

The influence of freshwater discharge and sea level rise in the salinity in the upper Tagusestuary is assessed using a numerical model. Results suggest that the river flow is the main driver of the salinity in the upper estuary and that salinity intrusion increases with the duration of the droughts. For the analyzed scenarios, salinity reaches concentrations that are inadequate for irrigation during some periods, which can harm the agricultural activities.

The Tagus estuary (Portugal) supports diverse uses and activities and itsecological value is well recognized. This study aimed to characterize the water quality inthe Tagus estuary during one year, as part of the Tagus estuary observatory. Three fieldcampaigns were performed in 2018 covering distinct seasons (Spring

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.

N/A

O desenvolvimento de actividades balneares, associadas ao turismo e lazer,tem conduzido a inúmeras intervenções para melhoramento e expansão depraias existentes. Mais raramente tem sido encarada a criação de praiasartificiais, em costas rochosas, sendo escassa a informação existente sobrea elaboração deste tipo de projectos.No sentido de contribuir para preencher esta lacuna, desenvolve-se umametodologia estruturada que contempla a análise dos aspectos e problemasdeterminantes da estabilidade de praias artificiais.

N/A

O presente relatório tem como objectivo a descrição da campanha de campo designada porMORFEED IX, realizada no âmbito do Projecto MORFEED (Mophodynamic feedback of estuarinemargins to climate change), financiado pela Fundação para a Ciência e Tecnologia (contratoPTDC/AAC-AMB/100092/2008).O principal objectivo do projecto MORFEED consiste no desenvolvimento de um modeloconceptual que traduza a resposta das margens estuarinas a diferentes cenários de alteraçõesclimáticas. Com vista à prossecução deste objectivo foi seleccionado o estuário do Tejo como localde estudo.Neste projecto participam o Laboratório Nacional de Engenharia Civil (LNEC), que coordena, ea Fundação da Faculdade de Ciências da Universidade de Lisboa (FFCUL) através dos centros deinvestigação Instituto Dom Luís (IDL), Centro de Geologia (CGUL) e Centro de Oceanografia (CO).A campanha de campo MORFEED IX realizou-se na Póvoa de Santa Íria (imediações doparque industrial da Solvay), margem direita do Estuário do Tejo, durante a manhã do dia 23 deJulho de 2012.

O trecho compreendido entre a foz do Douro e o cabo Mondego apresenta o maisacentuado processo de evolução da costa portuguesa, com erosões a sul deEspinho e a sul da embocadura da Ria de Aveiro e acumulação a norte destaembocadura. A comparação de linhas de água de diferentes datas, complementadacom a estimativa dos balanços de areia resultantes de perdas e ganhos ocorridosnas fronteiras norte e sul e no interior das três zonas referidas permitiu: quantificaras áreas e volumes perdidos pelo trecho de costa; caracterizar a evolução edistribuição espacial dessas perdas ao longo tempo; e interpretar de formaquantificada os processos de erosão e acumulação ocorridos entre a foz do Douro eo cabo Mondego.

Este estudo enquadra-se na área da investigação científica da dinâmica litoral e visa melhorar aavaliação e previsão da vulnerabilidade/resiliência de sistemas dunares, que constituem a maisimportante defesa natural costeira na interface terra-mar, à acção das ondas em condições detempestade marítima caracterizadas por fortes ventos, baixas pressões atmosféricas e curta duração(horas a dias). O objectivo do estudo foi avaliar o desempenho do modelo XBeach, com elevadopotencial para erosão de praias e dunas, relativamente a um dos modelos mais eficazes nesta área,vastamente aplicado a nível mundial em projectos de investigação e consultoria, o modelo Litprof (DHI,2008) do sistema de modelação Litpack. A avaliação do desempenho dos modelos foi feita com base nacomparação da evolução morfológica e em indicadores de desempenho e erro. Apresentam-se ediscutem-se os parâmetros de calibração por defeito (standard set of parameter settings) em ambos osmodelos com os quais se reproduziu o desenvolvimento de um perfil de erosão do tipo barra-fossaobtido em canal de grande escala de laboratório. Com base nas conclusões da avaliação dedesempenho dos dois modelos para estas condições (parâmetros por defeito), recomenda-se que,apesar do indicador de erro utilizado (BSS) apontar para um melhor desempenho do modelo XBeach,seja o modelo Litprof a aplicar num problema de engenharia uma vez que a geometria do perfil é a quemais se aproxima do perfil observado e a erosão da face de praia é sobrestimada, estando assim dolado da segurança. Realizou-se um conjunto de testes aos parâmetros de calibração para os doismodelos e concluiu-se que os parâmetros mais influentes na evolução morfológica foram break,gammax, facua e lws, no caso do modelo XBeach, e Gamma2 e scale parameter , no caso do modelo Litprof. Osparâmetros facua e lws do modelo XBeach foram aqueles que conferiram ao perfil final uma geometriamais próxima da configuração barra-fossa observada. Dos testes aos parâmetros de calibraçãorealizados concluiu-se que ambos os modelos reproduzem o desenvolvimento da barra submersacaracterístico do perfil de erosão, mas foi o modelo Litprof que apresentou o melhor desempenho.

Este estudo enquadra-se na área da investigação científica da dinâmica litoral e visa melhorar aavaliação e previsão da vulnerabilidade/resiliência de sistemas dunares à acção das ondas em condiçõesde tempestade marítima. O objectivo do estudo foi avaliar o desempenho dos modelos XBeach (comelevado potencial para erosão de praias e dunas) e Litprof do sistema de modelação Litpack (um dosmodelos mais eficazes nesta área, vastamente aplicado a nível mundial em projectos de investigação econsultoria) na erosão de dunas. Testaram-se os modelos para um caso de verificação de erosão deduna realizado em canal de grande escala de laboratório. A avaliação do desempenho dos modelos foifeita com base na análise da evolução morfológica e em indicadores de impacto (volume de erosão,recuo da duna e recuo do topo da duna) e de erro (Brier Skill Score). Apresentam-se e testam-se osparâmetros por defeito (standard set of parameter settings) para ambos os modelos. Da avaliação dedesempenho dos dois modelos com os parâmetros por defeito concluiu-se que o desempenho do modeloXBeach é razoável e do modelo Litprof é fraco. Recomenda-se que numa aplicação de engenharia paraprevisão da erosão de dunas, em que importa não só a precisão como também estar do lado dasegurança, o modelo XBeach com os parâmetros por defeito deve ser usado com precaução, uma vezque o volume de erosão foi subestimado. Apresentam-se, discutem-se e testam-se os parâmetros decalibração para ambos os modelos. Concluiu-se que os parâmetros mais influentes na evoluçãomorfológica foram: beta, break, facua, gammax, hswitch, lws e wetslp, para o modelo XBeach eMaximum Angle of Bed Slope, para o modelo Litprof. Os parâmetros lws e wetslp do modelo XBeachforam aqueles que conferiram ao perfil final uma geometria mais próxima da configuração observada. Osresultados obtidos para o indicador de erro BSS evidenciam que o melhor desempenho foi obtido comalteração do parâmetro lws e que o segundo melhor desempenho foi obtido com a alteração doparâmetro wetslp, sendo a ambos atribuída a classificação de excelente. Uma vez que a alteraçãomorfológica alcançada com a modificação do parâmetro lws permitiu melhorar a previsão das duascaracterísticas fundamentais sob o ponto de vista da engenharia que são o recuo do topo da duna e olimite da extensão da zona activa (onde se observou a formação da barra submersa durante aexperiência), considera-se que este parâmetro é de grande relevância na simulação da evolução daerosão de dunas com o modelo XBeach. No seu melhor desempenho o modelo Litprof calibradoreproduziu correctamente o recuo do topo da duna, simulou a formação de uma barra submersa naposição observada experimentalmente, simulou incorrectamente o volume de erosão (cerca de metadedo observado) e o declive da duna, e consequentemente gerou um avanço da duna ao nível do mar emvez de recuo. Por este motivo recomenda-se precaução na aplicação do modelo Litprof para previsão daerosão de dunas. Verificou-se também que a teoria de onda não foi relevante na evolução do perfil nestecaso de estudo. Da comparação dos modelos morfodinâmicos XBeach e Litprof nas duas fases, i.e., nafase de teste com os parâmetros por defeito e na fase de calibração, concluiu-se que foi o modeloXbeach que apresentou o melhor desempenho neste caso de estudo. A execução deste estudo permitiutestar e ficar a conhecer a elevada capacidade do modelo XBeach e a razoável capacidade do modeloLitprof na previsão da erosão de dunas.

N/A

N/A

N/A

N/A