Synthesis of FPGA-based accelerators implementing recursive algorithms

O desenvolvimento de sistemas computacionais é um processo complexo, com múltiplas etapas, que requer uma análise profunda do problema, levando em consideração as limitações e os requisitos aplicáveis. Tal tarefa envolve a exploração de técnicas alternativas e de algoritmos computacionais para optimizar o sistema e satisfazer os requisitos estabelecidos. Neste contexto, uma das mais importantes etapas é a análise e implementação de algoritmos computacionais. Enormes avanços tecnológicos no âmbito das FPGAs (Field-Programmable Gate Arrays) tornaram possível o desenvolvimento de sistemas de engenharia extremamente complexos. Contudo, o número de transístores disponíveis por chip está a crescer mais rapidamente do que a capacidade que temos para desenvolver sistemas que tirem proveito desse crescimento. Esta limitação já bem conhecida, antes de se revelar com FPGAs, já se verificava com ASICs (Application-Specific Integrated Circuits) e tem vindo a aumentar continuamente. O desenvolvimento de sistemas com base em FPGAs de alta capacidade envolve uma grande variedade de ferramentas, incluindo métodos para a implementação eficiente de algoritmos computacionais. Esta tese pretende proporcionar uma contribuição nesta área, tirando partido da reutilização, do aumento do nível de abstracção e de especificações algorítmicas mais automatizadas e claras. Mais especificamente, é apresentado um estudo que foi levado a cabo no sentido de obter critérios relativos à implementação em hardware de algoritmos recursivos versus iterativos. Depois de serem apresentadas algumas das estratégias para implementar recursividade em hardware mais significativas, descreve-se, em pormenor, um conjunto de algoritmos para resolver problemas de pesquisa combinatória (considerados enquanto exemplos de aplicação). Versões recursivas e iterativas destes algoritmos foram implementados e testados em FPGA. Com base nos resultados obtidos, é feita uma cuidada análise comparativa. Novas ferramentas e técnicas de investigação que foram desenvolvidas no âmbito desta tese são também discutidas e demonstradas.

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.

Ecotoxicological tests using Chironomus riparius

Os sistemas aquáticos naturais podem estar sujeitos frequentemente a entrada de tóxicos, quer seja através da lixiviação dos campos agrícolas ou da descarga por parte de unidades industriais. Avaliar o impacto potencial destes contaminantes nos sistemas aquáticos é muito importante, porque pode promover consequências sérias no balanço ecológico dos ecossistemas. Os efeitos de níveis sub-letais destes tóxicos nas populações aquáticas são detectados, em muitos casos, somente após diversas gerações, dependendo da espécie e do contaminante. O comportamento animal é considerado como sendo a primeira linha de defesa perante estímulos ambientais, e pode ser uma representação de alterações fisiológicas no organismo, sendo portanto um indicador excelente de alterações ambientais. O desenvolvimento dos sistemas de aviso prévio que integram parâmetros comportamentais pode ajudar a prever mais rapidamente possíveis alterações ao nível das populações naturais, do que a utilização de testes ecotoxicológicos padrão com a mesma finalidade. O conhecimento acerca de possíveis implicações devido a alterações comportamentais, em organismos bentónicos e em populações do campo sujeitas a tóxicos, é ainda escasso. Sabendo isto, neste estudo pretendeu-se investigar como o comportamento de Chironomus riparius – usando um biomonitor em tempo real – e outros parâmetros tais como crescimento, emergência de adultos, bioacumulação e biomarcadores, são afectados pela exposição a imidacloprid e ao mercúrio, que foram seleccionados como contaminantes. Os resultados demonstraram que a exposição às concentrações sub-letais de imidacloprid afecta o crescimento e o comportamento dos quironomídeos e que estes organismos podem recuperar de uma exposição curta ao insecticida. O comportamento que corresponde à ventilação de C. riparius revelou-se como um parâmetro mais sensível do que a locomoção e do que as respostas bioquímicas, quando as larvas foram sujeitas ao imidacloprid. Larvas de C. riparius expostas a concentrações sub-letais de mercúrio apresentaram uma tendência de diminuição de actividade comportamental, em testes com concentrações crescentes do tóxico; o crescimento das larvas foi também prejudicado, e as taxas de emergência de adultos e o tempo de desenvolvimento apresentaram retardamento. Estes organismos podem bioacumular rapidamente o mercúrio em condições de não alimentação e apresentam uma lenta depuração deste metal. Estes efeitos podem, em último caso, conduzir a prováveis repercussões ao nível da população e das comunidades. As reduções em actividades comportamentais, mesmo em concentrações baixas, podem diminuir a quantidade de tempo gasta na procura de alimento, produzindo efeitos aos níveis morfo-fisiológicos, e assim afectar severamente o desempenho dos quironomídeos no ambiente. O uso destes factores comportamentais como um parâmetro ecotoxicológico sub-letal relevante ao nível da toxicologia aumentará a versatilidade dos testes, permitindo uma resposta comportamental mensurável e quantitativa ao nível do organismo, utilizando uma avaliação não destrutiva, e assim certificando que esta aproximação pode ser usada em testes ecotoxicológicos futuros.

Bottom up strategies for the morphology control of BaTiO3 particulates

As propriedades funcionais dos materiais ferroeléctricos tais como a polarização reversível, piroelectricidade, piezoelectricidade, elevada actividade óptica não linear e comportamento dieléctrico não linear são fundamentais para a sua aplicação em sensores, microactuadores, detectores de infravermelhos, filtros de fase de microondas e memórias não-voláteis. Nos últimos anos, motivado pelas necessidades industriais de redução do tamanho dos dispositivos microelectrónicos, aumentando a eficiência volumétrica, tem sido feito um grande esforço ao nível da investigação para desenvolver estruturas ferroeléctricas à escala micro- e nano- métrica. É sabido que a redução de tamanho em materiais ferroeléctricos afecta significamente as suas propriedades. Neste sentido e considerando que foi previsto teoreticamente por cálculos ab initio que estruturas do tipo nanocilindros e nanodiscos apresentariam um novo tipo de ordem ferroeléctrica e, na expectativa de alcançar conhecimento para o desenvolvimento de uma nova geração de dispositivos microelectróncos, existe um grande interesse em desenvolver métodos de fabrico de nanoestruturas ferroeléctricas unidimensionais (1D) tais como nanocilindros e nanotubos. As estratégias de fabrico de nanoestruturas 1D até agora descritas na literatura indicam claramente as dificuldades inerentes à sua preparação. Existem duas grandes vias de síntese destas nanoestruturas: i) o método “topdown” que consiste na redução de tamanho de um dado material até à obtenção duma estrutura 1D; e ii) o método “bottom-up” em que átomos, iões e moléculas são agrupados para formar um material 1D. O método “top down” envolve em geral técnicas de desgaste, como o uso do feixe de electrões, que apesar de permitirem elevada precisão no posicionamento e no controlo do tamanho, falham em termos de resolução, exigem muito tempo e causam facilmente defeitos que deterioram as propriedades físicas destes materiais. Na metodologia “bottom up” a utilização de moléculas ou estruturas “molde” tem sido a mais explorada. As estructuras 1D podem também ser preparadas sem recorrer a “moldes”. Neste caso a agregação orientada é promovida pelo recurso a aditivos que controlam o crescimento dos cristais em direcções preferenciais. Neste contexto, neste trabalho utilizaram-se duas estratégias “bottom up” de baixo custo para a preparação de nanopartículas de titanato de bário (BaTiO3) com morfologia controlada: 1) síntese química (em solução e em fase vapor) com utilização de nanotubos de titanato TiNTs) como “moldes” e precursores de titânio 2) síntese química em solução com presença de aditivos. Os nanotubos de titanato de sódio foram preparados por síntese hidrotermal. Como existiam muitas dúvidas acerca da natureza estrutural e do mecanismo de formação dos NTs, a parte inicial do trabalho foi dedicada à realização de um estudo sistemático dos parâmetros intervenientes na síntese e à caracterização da sua estrutura e microestrutura. Foi demonstrado que os NTs têm a fórmula geral A2Ti2O5 (A = H+ or Na+), e não TiO2 (anátase) com defendido por vários autores na literatura, e podem ser preparados por método hidrotermal em meio fortemente alcalino usando como fonte de titânio TiO2 comercial na forma de anátase ou rútilo. A menor reactividade do rútilo exige temperaturas de síntese superiores ou tempos de reacção mais longos. A forma tubular resulta do tratamento hidrotermal e não de processos de lavagem e neutralização subsequentes. Se os NTs forem tratados após a síntese hidrotérmica em água a 200 ºC, transformam-se em nanocilindros. Uma das partes principais desta tese consistiu na investigação do papel dos NTs de titanato no crescimento anisotrópico de BaTiO3. O potencial funcionamento dos NTs como “moldes” para além de precursores foi testado em reacção com hidróxido de bário em síntese em solução e por reacção com um precursor orgânico de bário em fase vapor. Tendo por base os estudos cinéticos realizados, bem como as alterações estruturais e morfológicas das amostras, é possível concluir que a formação do BaTiO3 a partir de NTs de titanato de sódio, ocorre por dois mecanismos dependendo da temperatura e tempo de reacção. Assim, a baixa temperatura e curto tempo de reacção verifica-se que se formam partículas dendríticas de BaTiO3 cuja superfície é bastante irregular (“wild”) e que apresentam estrutura pseudo-cúbica. Estas partículas formam-se por reacção topotáctica na fronteira dos nanotubos de titanato de sódio. A temperaturas mais altas e/ou reacções mais longas, a reacção é controlada por um mecanismo de dissolução e precipitação com formação de dendrites de BaTiO3 tetragonais com superfície mais regular (“seaweed”). A microscopia de força piezoeléctrica mostrou que as dendrites “seaweeds“ possuem actividade piezoeléctrica superior à das dendrites “wild”, o que confirma o papel desempenhado pela estrutura e pela concentração de defeitos na rede na coerência e ordem ferroeléctrica de nanoestruturas. Os nossos resultados confirmam que os NTs de titanato não actuam facilmente como “moldes” na síntese em solução de BaTiO3 já que a velocidade de dissolução dos NTs em condições alcalinas é superior à velocidade de formação do BaTiO3. Assumindo que a velocidade de reacção dos NTs com o precursor de bário é superior em fase vapor, efectuou-se a deposição de um precursor orgânico de bário por deposição química de vapor sobre um filme de NTs de titnato de sódio depositados por deposição electroforética. Estudou-se a estabilidade dos NTs nas diferentes condições do reactor. Quando os NTs são tratados a temperaturas superiores a 700 ºC, ocorre a transformação dos NTs em nanocilindros de anatase por um mecanismo de agregação orientada. Quando se faz a deposição do precursor de bário, seguida de calcinação a 700 ºC em atmosfera oxidante de O2, verifica-se que a superficie dos NTs fica coberta com nanocristais de BaTiO3 independentemente da concentração de bário. O papel dos NTs de titanato no crescimento anisotrópico de BaTiO3 em fase vapor é assim descrito pela primeira vez. Em relação à metodologias de crescimento de partículas na ausência de “moldes” mas com aditivos fez-se um estudo sistemático utilizando 5 aditivos de natureza differente. As diferenças entre aditivos foram sistematizadas tendo em conta as diferenças estruturais e morfológicas verificadas. Está provado que os aditivos podem funcionar como modificadores de crescimento cristalino por alteração do seu padrão de crescimento ou por alteração da cinética de crescimento das faces cristalográficas do cristal. Entre os aditivos testados verificou-se que o ácido poliacrilíco adsorve em faces específicas do BaTiO3 alterando a cinética de crescimento e induzindo a agregação orientada das partículas. O polivinilpirrolidona, o docecilsulfato de sódio e hidroxipropilmetilcelulose actuam mais como inibidores de crescimento do que como modificadores do tipo de crescimento. A D-frutose aumenta a energia de activação da etapa de nucleação não ocorrendo formação de BaTiO3 para as mesmas condições dos outros aditivos. Esta tese clarifica o papel dos NTs de titanato de sódio enquanto precursores e “moldes” no crescimento anisotrópico de BaTiO3 em solução e em fase vapor. É feita também a abordagem do controlo morfológico do BaTiO3 através do uso de aditivos. As estratégias de preparação de BaTiO3 propostas são de baixo custo, reprodutíveis e fáceis de efectuar. Os resultados contribuem para uma melhor compreensão da relação tamanho – morfologia – propriedade em materiais ferroeléctricos nanométricos com vista à sua potencial aplicação.

Critical magnetic elds in superconducting systems with semi-metallic bands

In this work, we study the Zeeman splitting effects in the parallel magnetic field versus temperature phase diagram of two-dimensional superconductors with one graphene-like band and the orbital effects of perpendicular magnetic fields in isotropic two-dimensional semi-metallic superconductors. We show that when parallel magnetic fields are applied to graphene and as the intraband interaction decreases to a critical value, the width of the metastability region present in the phase diagram decreases, vanishing completely at that critical value. In the case of two-band superconductors with one graphene-like band, a new critical interaction, associated primarily with the graphene-like band, is required in order for a second metastability region to be present in the phase diagram. For intermediate values of this interaction, a low-temperature first-order transition line bifurcates at an intermediate temperature into a first-order transition between superconducting phases and a second-order transition line between the normal and the superconducting states. In our study on the upper critical fields in generic semi-metallic superconductors, we find that the pair propagator decays faster than that of a superconductor with a metallic band. As result, the zero field band gap equation does not have solution for weak intraband interactions, meaning that there is a critical intraband interaction value in order for a superconducting phase to be present in semi-metallic superconductors. Finally, we show that the out-of-plane critical magnetic field versus temperature phase diagram displays a positive curvature, contrasting with the parabolic-like behaviour typical of metallic superconductors.