Hydraulics of stepped chutes: experimental-numerical-theoretical study

O estudo de escoamentos turbulentos em descarregadores em degraus tem sido um desafio para os investigadores. A macro-rugosidade do leito, a ondulação da superfície livre, a intermitência da localização da secção inicial de entrada de ar e o escoamento bi-fásico a jusante da secção inicial de entrada de ar fazem com que a caracterização do escoamento deslizante sobre turbilhões em descarregadores em degraus não seja simples. Actualmente, é possível combinar técnicas de medição fiáveis com simulações numéricas e análise teórica. Nesta dissertação, o estudo experimental baseia-se em resultados experimentais obtidos em duas instalações experimentais: a instalação A, do Laboratório Nacional de Engenharia Civil (LNEC), dotada de um descarregador em degraus com declive de 1V:0.75H e a instalação B, do Instituto Superior Técnico (IST), dotada de um descarregador em degraus com declive de 1V:2H. Uma sonda de ar, um tubo de Pitot modificado e vários hidrómetros permitiram o estudo do escoamento deslizante sobre turbilhões mono- e bifásico ao longo do descarregador em degraus da instalação A. Na bacia de dissipação de energia, a observação visual do escoamento foi auxiliada pelas leituras efectuadas em tomadas de pressão localizadas na soleira da bacia e pelas leituras da altura do escoamento efectuadas em réguas graduadas localizadas nas paredes da bacia. Na instalação B foram utilizados um tubo de Pitot e vários hidrómetros para estudar a região não arejada do escoamento deslizante sobre turbilhões no descarregador em degraus. O estudo numérico é baseado em simulações numéricas da região não arejada do escoamento deslizante sobre turbilhões sobre descarregadores em degraus com o código comercial de CFD FLOW-3D® de modo a reproduzir as condições ensaiadas experimentalmente. As simulações beneficiaram da técnica de blocos múltiplos (multi-block) num sistema de coordenadas cartesianas, da determinação da superfície livre pelo método TruVOF e da utilização de dois modelos de turbulência: os modelos k- e RNG k-. Por último, o estudo teórico consistiu em desenvolver um modelo simplificado 1D para determinar as características hidráulicas principais do trecho não arejado do escoamento deslizante sobre turbilhões em descarregadores em degraus. O modelo foi desenvolvido a partir das equações de Navier-Stokes, conjuntamente com resultados experimentais e numéricos. Os resultados apresentados nesta dissertação contribuem para o conhecimento do escoamento deslizante sobre turbilhões em descarregadores em degraus, nomeadamente na região não arejada, na secção inicial de entrada de ar e na região arejada. A hidráulica dos dissipadores de energia a jusante de descarregadores em degraus, em particular das bacias tipo III do USBR, é também objecto de estudo. Em relação à região não arejada do escoamento em descarregadores com declive acentuado, e com base em resultados experimentais e numéricos, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, da espessura da camada limite, da concentração média de ar, do coeficiente de energia cinética, da dissipação de energia, do factor de resistência e do coeficiente n da fórmula de Manning. São ainda propostas expressões adimensionais para a energia cinética turbulenta e sua dissipação. Para declives moderados, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, do coeficiente de energia cinética e da energia específica residual. São ainda apresentados valores do expoente 1/N da expressão adimensional da distribuição de velocidades, quer para descarregadores com declive acentuado quer com declive moderado. Em conformidade com outros estudos centrados em escoamentos de parede e com derivações teóricas, para a região não arejada do escoamento em descarregadores em degraus de acentuado declive, observa-se que o factor de resistência depende da macro-rugosidade criada pelos degraus e da geometria da secção transversal e que o coeficiente n da fórmula de Manning aumenta com a rugosidade. A descrição estatística da turbulência do escoamento é igualmente explorada, contribuindo para o conhecimento da estrutura do escoamento. Observou-se que para números de Reynolds rugoso não superiores a 6.8x104 a energia cinética turbulência e a sua dissipação cumprem leis de semelhança. Estas expressões adimensionais estão de acordo com os resultados obtidos por outros autores para escoamentos completamente desenvolvidos em canais abertos e no escoamento em rios com leito de gravilha. Em acréscimo, a taxa de dissipação de energia, quer para descarregadores de declive acentuado quer de moderado declive, é baixa. Por último, observa-se que os valores da média temporal da concentração de ar entre 0 e 1 medidos na região não arejada do escoamento dizem respeito não só ao ar capturado entre ondas de água, na zona de ondulação da superfície livre, mas também ao ar emulsionado no escoamento, i.e., sob a forma de bolhas de ar, quando perto da secção média inicial de entrada de ar, devido à diferença entre localizações instantânea e média temporal. Foram revistas metodologias e fórmulas para estimar a localização da secção inicial de entrada de ar e apresentadas expressões para estimar a concentração média de ar e a altura equivalente de água nessa secção. Relativamente à região de escoamento arejado em descarregadores em degraus com declive acentuado, os resultados experimentais apresentados nesta dissertação permitiram estimar a influência da definição da superfície livre nos parâmetros hidráulicos da região do escoamento arejado e estimar a máxima elevação do escoamento nesta região do escoamento. Com base nos resultados experimentais obtidos na bacia de dissipação de energia do tipo III do USBR localizada a jusante do descarregador em degraus da instalação A, observou-se que os perfis da altura piezométrica e da altura do escoamento tendem a seguir o perfil recomendado pelo USBR para bacias tipo III. A excepção ocorre à entrada da bacia, onde as alturas piezométricas apresentadas nesta dissertação excedem largamente as apresentadas pelo USBR. É ainda observado que, tal como entre as bacias tipo I e tipo III do USBR, o ressalto hidráulico estabiliza muito mais rapidamente numa bacia tipo III a jusante de um descarregador em degraus do que uma bacia tipo I a jusante do mesmo descarregador em degraus. Finalmente, observa-se que os blocos de amortecimento a colocar no descarregador não têm influência visível nos resultados da altura piezométrica nem da altura do escoamento ao longo da bacia. Relativamente às simulações numéricas do escoamento não arejado, a proximidade entre resultados experimentais e numéricos permite validar o modelo teórico e a integração numérica usados no FLOW-3D®. As simulações desenvolvidas também mostraram que o modelo de turbulência k- permite representar as características do escoamento não arejado em descarregadores em degraus, uma vez que não foram observadas diferenças significativas entre as simulações com este modelo e com o modelo RNG k-. Finalmente, observou-se que o modelo de entrada de ar usado no FLOW-3D® é válido para estimar a localização da secção inicial de entrada de ar. Por último, a proximidade entre os resultados obtidos da aplicação do modelo teórico desenvolvido no âmbito desta dissertação e os resultados experimentais indica que as hipóteses e simplificações consideradas no desenvolvimento do modelo são adequadas.

Microstructure design of Titanate-based electroceramics

Electrocerâmicos são uma classe de materiais avançados com propriedades eléctricas valiosas para aplicações. Estas propriedades são geralmente muito dependentes da microestrutura dos materiais. Portanto, o objectivo geral deste trabalho é investigar o desenho da resposta dieléctrica de filmes espessos obtidos por Deposição Electroforética (EPD) e cerâmicos monolíticos, através do controlo da evolução da microestrutura durante a sinterização de electrocerâmicos à base de titanatos. Aplicações sem fios na indústria microelectrónica e de comunicações, em rápido crescimento, tornaram-se um importante mercado para os fabricantes de semicondutores. Devido à constante necessidade de miniaturização, redução de custos e maior funcionalidade e integração, a tecnologia de filmes espessos está a tornar-se uma abordagem de processamento de materiais funcionais cada vez mais importante. Uma técnica adequada neste contexto é EPD. Os filmes espessos resultantes necessitam de um passo subsequente de sinterização que é afectada pelo substrato subjacente, tendo este um forte efeito sobre a evolução da microestrutura. Relacionado com a miniaturização e a discriminação do sinal, materiais dieléctricos usados como componentes operando a frequências das microondas em aplicações na industria microelectrónica de comunicações devem apresentar baixas perdas dieléctricas e elevadas permitividade dieléctrica e estabilidade com a temperatura. Materiais do sistema BaO-Ln2O3- TiO2 (BLnT: Ln = La ou Nd), como BaLa4Ti4O15 (BLT) e Ba4.5Nd9Ti18O54 (BNT), cumprem esses requisitos e são interessantes para aplicações, por exemplo, em estações de base para comunicações móveis ou em ressonadores para telefones móveis, onde a miniaturização dos dispositivos é muito importante. Por sua vez, o titanato de estrôncio (SrTiO3, STO) é um ferroeléctrico incipiente com constante dieléctrica elevada e baixas perdas, que encontra aplicação em, por exemplo, condensadores de camada interna, tirando partido de fronteiras de grão altamente resistivas. A dependência da permitividade dieléctrica do campo eléctrico aplicado torna este material muito interessante para aplicações em dispositivos de microondas sintonizáveis. Materiais à base de STO são também interessantes para aplicações termoeléctricas, que podem contribuir para a redução da actual dependência de combustíveis fósseis por meio da geração de energia a partir de calor desaproveitado. No entanto, as mesmas fronteiras de grão resistivas são um obstáculo relativamente à eficiência do STO para aplicações termoeléctricas. Para além do efeito do substrato durante a sinterização constrangida, outros factores, como a presença de fase líquida, a não-estequiometria ou a temperatura de sinterização, afectam significativamente não apenas a microestrutura dos materiais funcionais, mas também a sua resposta dieléctrica. Se adequadamente compreendidos, estes factores podem ser intencionalmente usados para desenhar a microestrutura dos electrocerâmicos e, desta forma, as suas propriedades dieléctricas. O efeito da não-estequiometria (razão Sr/Ti 0.995-1.02) no crescimento de grão e resposta dieléctrica de cerâmicos de STO foi investigado neste trabalho. A mobilidade das fronteiras de grão aumenta com a diminuição da razão Sr/Ti. A resistividade do interior dos grãos e das fronteiras de grão é sistematicamente diminuída em amostras não-estequiométricas de STO, em comparação com o material estequiométrico. O efeito é muito mais forte para as fronteiras de grão do que para o seu interior. Dependências sistemáticas da não-estequiometria foram também observadas relativamente à dependência da condutividade da temperatura (muito mais afectada no caso da contribuição das fronteiras de grão), à capacitância do interior e fronteiras de grão e à espessura das fronteiras de grão. Uma anomalia no crescimento de grão em cerâmicos de STO ricos em Ti foi também observada e sistematicamente analisada. Foram detectadas três descontinuidades na dependência do tipo Arrhenius do crescimento de grão relativamente à temperatura com diminuições no tamanho de grão a temperaturas em torno de 1500, 1550 e 1605 °C. Além disso, descontinuidades semelhantes foram também observadas na dependência da energia de activação relativamente à condutividade das fronteiras de grão e na espessura das fronteiras de grão, avaliadas por Espectroscopia de Impedância. Estas notáveis coincidências suportam fortemente a formação de diferentes complexos de fronteira de grão com transições entre os regimes de crescimento de grão observados, que podem ser correlacionados com diferentes mobilidades de fronteira de grão e propriedades dieléctricas. Um modelo é sugerido, que se baseia na diminuição da fase líquida localizada nas fronteiras de grão, como o aumento da temperatura de sinterização, um cenário compatível com um fenómeno de solubilidade retrógrada, observado anteriormente em metais e semicondutores, mas não em cerâmicos. A EPD de filmes espessos de STO em substratos de folha de Pt e a sinterização constrangida dos filmes fabricados foram também preliminarmente tratadas. Filmes espessos de STO foram depositados com êxito por EPD sobre substratos de Pt e, depois de sinterizados, atingiram densidades elevadas. Um aumento da densificação e do tamanho de grão assim como o alargamento da distribuição de tamanho do grão foram observados com a diminuição da razão Sr/Ti, tal como anteriormente observado em amostras cerâmicas. Grãos equiaxiados foram observados para todas as composições, mas um certo grau de anisotropia na orientação dos poros foi detectado: os poros revelaram uma orientação vertical preferencial. Este trabalho focou-se também na sinterização constrangida do sistema BLnT (Ln = La ou Nd), nomeadamente de filmes espessos de BLT e BNT sobre substratos de folha de platina, e na relação do desenvolvimento de anisotropia microestrutural com as propriedades dieléctricas. As observações durante a sinterização constrangida foram comparadas com cerâmicos monolíticos equivalentes sinterizados livremente. Filmes espessos de BLnT (Ln = La ou Nd) com elevada densidade foram obtidos por EPD e subsequente sinterização constrangida. A anisometria cristalográfica do material em conjunto com um passo de sinterização constrangida resultou em grãos alongados e microestruturas anisotrópicas. O efeito do stress do substrato durante a sinterização constrangida originou graus mais elevados de anisotropia (grãos e poros alongados e orientação preferencial, bem como textura cristalográfica) nos filmes sinterizados relativamente aos cerâmicos equivalentes sinterizados livremente, não obstante o estado equivalente das amostras em verde. A densificação dos filmes de BLnT (Ln = La ou Nd) é retardada em comparação com os cerâmicos, mas depois de longos tempos de sinterização densidades semelhantes são obtidas. No entanto, em oposição a observações na sinterização constrangida de outros sistemas, o crescimento do grão em filmes de BLnT (Ln = La ou Nd) é favorecido pelo constrangimento causado pelo substrato. Além disso, grãos e poros alongados orientados paralelamente ao substrato foram desenvolvidos durante a sinterização constrangida de filmes espessos. Verificou-se uma forte correlação entre a evolução de grãos e poros, que começou assim que o crescimento do grão se iniciou. Um efeito da tensão do substrato no aumento do crescimento de grão, bem como um forte “Zener pinning”, origina microestruturas altamente texturizadas, o que também é observado a nível cristalográfico. Efeitos marcantes da anisotropia microestrutural foram também detectados nas propriedades dieléctricas dos filmes de BLnT (Ln = La ou Nd). Juntamente com o aumento da razão de aspecto dos grãos, do factor de orientação e do grau de textura cristalográfica, a permitividade relativa é ligeiramente diminuída e o coeficiente de temperatura da permitividade evolui de negativo para positivo com o aumento do tempo isotérmico de sinterização. Este trabalho mostra que a não-estequiometria pode ser usada para controlar a mobilidade das fronteiras de grão e, portanto, desenhar a microestrutura e as propriedades dieléctricas de electrocerâmicos à base de STO, com ênfase nas propriedades das fronteiras de grão. O papel da não-estequiometria no STO e dos complexos de fronteira de grão no desenvolvimento microestrutural é discutido e novas oportunidades para desenhar as propriedades de materiais funcionais são abertas. As observações relativamente à sinterização constrangida apontam para o efeito de tensões mecânicas desenvolvidas devido ao substrato subjacente no desenvolvimento da microestrutura de materiais funcionais. É assim esperado que a escolha adequada de substrato permitia desenhar a microestrutura de filmes espessos funcionais com desempenho optimizado. “Stress Assisted Grain Growth” (SAGG) é então proposto como uma técnica potencial para desenhar a microestrutura de materiais funcionais, originando microestruturas anisotrópicas texturizadas com propriedades desejadas.

Wind energy resource modelling in Portugal and its future large-scale alteration due to anthropogenic induced climate changes

The high dependence of Portugal from foreign energy sources (mainly fossil fuels), together with the international commitments assumed by Portugal and the national strategy in terms of energy policy, as well as resources sustainability and climate change issues, inevitably force Portugal to invest in its energetic self-sufficiency. The 20/20/20 Strategy defined by the European Union defines that in 2020 60% of the total electricity consumption must come from renewable energy sources. Wind energy is currently a major source of electricity generation in Portugal, producing about 23% of the national total electricity consumption in 2013. The National Energy Strategy 2020 (ENE2020), which aims to ensure the national compliance of the European Strategy 20/20/20, states that about half of this 60% target will be provided by wind energy. This work aims to implement and optimise a numerical weather prediction model in the simulation and modelling of the wind energy resource in Portugal, both in offshore and onshore areas. The numerical model optimisation consisted in the determination of which initial and boundary conditions and planetary boundary layer physical parameterizations options provide wind power flux (or energy density), wind speed and direction simulations closest to in situ measured wind data. Specifically for offshore areas, it is also intended to evaluate if the numerical model, once optimised, is able to produce power flux, wind speed and direction simulations more consistent with in situ measured data than wind measurements collected by satellites. This work also aims to study and analyse possible impacts that anthropogenic climate changes may have on the future wind energetic resource in Europe. The results show that the ECMWF reanalysis ERA-Interim are those that, among all the forcing databases currently available to drive numerical weather prediction models, allow wind power flux, wind speed and direction simulations more consistent with in situ wind measurements. It was also found that the Pleim-Xiu and ACM2 planetary boundary layer parameterizations are the ones that showed the best performance in terms of wind power flux, wind speed and direction simulations. This model optimisation allowed a significant reduction of the wind power flux, wind speed and direction simulations errors and, specifically for offshore areas, wind power flux, wind speed and direction simulations more consistent with in situ wind measurements than data obtained from satellites, which is a very valuable and interesting achievement. This work also revealed that future anthropogenic climate changes can negatively impact future European wind energy resource, due to tendencies towards a reduction in future wind speeds especially by the end of the current century and under stronger radiative forcing conditions.