Development of procedures for the simulation of atmospheric flows over complex terrain, using OpenFOAM

The main purpose of this work was the development of procedures for the simulation of atmospheric ows over complex terrain, using OpenFOAM. For this aim, tools and procedures were developed apart from this code for the preprocessing and data extraction, which were thereafter applied in the simulation of a real case. For the generation of the computational domain, a systematic method able to translate the terrain elevation model to a native OpenFOAM format (blockMeshDict) was developed. The outcome was a structured mesh, in which the user has the ability to de ne the number of control volumes and its dimensions. With this procedure, the di culties of case set up and the high computation computational e ort reported in literature associated to the use of snappyHexMesh, the OpenFOAM resource explored until then for the accomplishment of this task, were considered to be overwhelmed. Developed procedures for the generation of boundary conditions allowed for the automatic creation of idealized inlet vertical pro les, de nition of wall functions boundary conditions and the calculation of internal eld rst guesses for the iterative solution process, having as input experimental data supplied by the user. The applicability of the generated boundary conditions was limited to the simulation of turbulent, steady-state, incompressible and neutrally strati ed atmospheric ows, always recurring to RaNS (Reynolds-averaged Navier-Stokes) models. For the modelling of terrain roughness, the developed procedure allowed to the user the de nition of idealized conditions, like an uniform aerodynamic roughness length or making its value variable as a function of topography characteristic values, or the using of real site data, and it was complemented by the development of techniques for the visual inspection of generated roughness maps. The absence and the non inclusion of a forest canopy model limited the applicability of this procedure to low aerodynamic roughness lengths. The developed tools and procedures were then applied in the simulation of a neutrally strati ed atmospheric ow over the Askervein hill. In the performed simulations was evaluated the solution sensibility to di erent convection schemes, mesh dimensions, ground roughness and formulations of the k - ε and k - ω models. When compared to experimental data, calculated values showed a good agreement of speed-up in hill top and lee side, with a relative error of less than 10% at a height of 10 m above ground level. Turbulent kinetic energy was considered to be well simulated in the hill windward and hill top, and grossly predicted in the lee side, where a zone of ow separation was also identi ed. Despite the need of more work to evaluate the importance of the downstream recirculation zone in the quality of gathered results, the agreement between the calculated and experimental values and the OpenFOAM sensibility to the tested parameters were considered to be generally in line with the simulations presented in the reviewed bibliographic sources.

A qualitative description of soil parameters variation due to a prescribed fire in Portuguese northwestern forests using Fuzzy Boolean Nets — The case study of Cabreira mountain

In this paper we study the modifications that occurred in some forest soil properties after a prescribed fire. The research focused on the alterations of soil pH, soil moisture and soil organic matter content during a two-year span, from 2008 to 2009. The study site is located in Anjos, Vieira do Minho municipality, a forest site that has suffered from recurrent wildfires for several decades. Furze (Ulex, sp.), broom (Cytisus, sp.), gorse (Chamaespartum tridentatum) and a very few disperse adult pine (Pinus sylvestris) are the predominant vegetation type in the study area. The average height of this shrub vegetation is around 1.5 m. The prescribed fire was conducted by the National Forestry Authority (AFN) in November 2008. Fuzzy Boolean Nets (FBN) were used to evaluate the alteration in soil parameters when compared with adjacent spots where: i) no fire occurrence was registered since 1998; ii) fire occurrence was registered in 2008; and iii) vegetation pruning by mechanical cut was done in Spring six months prior to the prescribed fire event. Results suggest that in the particular case of the studied site, Anjos, the observed soil properties alterations cannot be related with the prescribed fire.

Measurements of environmental background radiation in the surrounding area of a coal-fired power plant

Certain materials used and produced in a wide range of non-nuclear industries contain enhanced activity concentrations of natural radionuclides. In particular, electricity production from coal is one of the major sources of increased exposure to man from enhanced naturally occurring materials. Over the past decades there has been some discussion about the elevated natural background radiation in the area near coal-fired power plants due to high uranium and thorium content present in coal. This work describes the methodology developed to assess the radiological impact due to natural radiation background increasing levels, potentially originated by a coal-fired power plant’s operation. Gamma radiation measurements have been done with two different instruments: a scintillometer (SPP2 NF, Saphymo) and a gamma ray spectrometer with energy discrimination (Falcon 5000, Canberra). A total of 40 relevant sampling points were established at locations within 20 km from the power plant: 15 urban and 25 suburban measured stations. The highest values were measured at the sampling points near to the power plant and those located in the area within the 6 and 20 km from the stacks. This may be explained by the presence of a huge coal pile (1.3 million tons) located near the stacks contributing to the dispersion of unburned coal and, on the other hand, the height of the stacks (225 m) which may influence ash’s dispersion up to a distance of 20 km. In situ gamma radiation measurements with energy discrimination identified natural emitting nuclides as well as their decay products (212Pb, 214Pb, 226Ra 232Th, 228Ac, 234Th 234Pa, 235U, etc.). This work has been primarily done to in order to assess the impact of a coal-fired power plant operation on the background radiation level in the surrounding area. According to the results, an increase or at least an influence has been identified both qualitatively and quantitatively.

Avaliação isocinética da força muscular dos isquiotibiais e quadricípite em atletas da equipa nacional portuguesa de taekwondo

Objectivos: Avaliar a força muscular bilateralmente ao nível dos músculos quadricípite e isquiotibiais em atletas da equipa nacional Portuguesa de Taekwondo. Amostra: Foi constituída por 6 dos 10 elementos que constituem a população de Taekwondistas do sexo masculino inscritos na Federação Portuguesa de Taekwondo, com presença regular nas selecções nacionais todos com um mínimo de 6 anos de prática. Os atletas apresentaram uma idade média 17,5 (+ 1,9) anos, com uma altura de 181,2 (+ 2,8) cm e com uma massa corporal total de 74,1 (+ 11,6) kg Metodologia: Procedeu-se à avaliação da força muscular dos participantes ao nível da musculatura flexora e extensora do joelho, no dinamómetro Biodex System 4, a uma velocidade de execução de 60º/s (4 repetições) e de 180º/s (10 repetições) com 60 segundos de intervalo, numa amplitude de movimento compreendida entre os 90 e os 0 graus. Todos os dados foram tratados no programa SPSS, versão 18.0, com um nível de significância de 0,05. Resultados: Verificou-se a existência de diferenças estatísticas significativas na análise do Peak Torque (p=0,023) aos 180º/seg para os flexores do joelho e do Peak Torque 30º (p=0,023) aos 180º/seg na acção dos extensores e flexores (p=0,037) do joelho, entre o membro dominante e não dominante. Constatou-se ainda a existência de um rácio Isquiotibiais/Quadricípite (55%) dentro dos valores normais do equilíbrio muscular do joelho. Discussão/Conclusão: Concluiu-se que foram encontradas diferenças nos níveis de força obtidos entre o membro dominante e não dominante, no entanto não foram encontrados desequilíbrios musculares clinicamente significativos. Assim como na relação Isquiotibiais/Quadricípite do próprio membro, não apontando por isso risco de lesão articular do joelho.

Dimensionamento de um secador em leito fluidizado para secagem de cereais

O objectivo desta tese é dimensionar um secador em leito fluidizado para secagem de cereais, nomeadamente, secagem de sementes de trigo. Inicialmente determinaram-se as condições de hidrodinâmica (velocidade de fluidização, TDH, condições mínimas de “slugging”, expansão do leito, dimensionamento do distribuidor e queda de pressão). Com as condições de hidrodinâmica definidas, foi possível estimar as dimensões físicas do secador. Neste ponto, foram realizados estudos relativamente à cinética da secagem e à própria secagem. Foi também estudado o transporte pneumático das sementes. Deste modo, determinaram-se as velocidades necessárias ao transporte pneumático e respectivas quedas de pressão. Por fim, foi realizada uma análise custos para que se soubesse o custo deste sistema de secagem. O estudo da secagem foi feito para uma temperatura de operação de 50ºC, tendo a ressalva que no limite se poderia trabalhar com 60ºC. A velocidade de operação é de 2,43 m/s, a altura do leito fixo é de 0,4 m, a qual sofre uma expansão durante a fluidização, assumindo o valor de 0,79 m. O valor do TDH obtido foi de 1,97 m, que somado à expansão do leito permite obter uma altura total da coluna de 2,76 m. A altura do leito fixo permite retirar o valor do diâmetro que é de 0,52 m. Verifica-se que a altura do leito expandido é inferior à altura mínima de “slugging” (1,20 m), no entanto, a velocidade de operação é superior à velocidade mínima de “slugging” (1,13 m/s). Como só uma das condições mínimas é cumprida, existe a possibilidade da ocorrência de “slugging”. Finalmente, foi necessário dimensionar o distribuidor, que com o diâmetro de orifício de 3 mm, valor inferior ao da partícula (3,48 mm), permite a distruibuição do fluido de secagem na coluna através dos seus 3061 orifícios. O inicio do estudo da secagem centrou-se na determinação do tempo de secagem. Além das duas temperaturas atrás referidas, foram igualmente consideradas duas humidades iniciais para os cereais (21,33% e 18,91%). Temperaturas superiores traduzem-se em tempos de secagem inferiores, paralelamente, teores de humidade inicial inferiores indicam tempos menores. Para a temperatura de 50ºC, os tempos de secagem assumiram os valores de 2,8 horas para a 21,33% de humidade e 2,7 horas para 18,91% de humidade. Foram também tidas em conta três alturas do ano para a captação do ar de secagem, Verão e Inverno representando os extremos, e a Meia- Estação. Para estes três casos, foi possível verificar que a humidade específica do ar não apresenta alterações significativas entre a entrada no secador e a corrente de saída do mesmo equipamento, do mesmo modo que a temperatura de saída pouco difere da de entrada. Este desvio de cerca de 1% para as humidades e para as temperaturas é explicado pela ausência de humidade externa nas sementes e na pouca quantidade de humidade interna. Desta forma, estes desvios de 1% permitem a utilização de uma razão de reciclagem na ordem dos 100% sem que o comportamento da secagem se altere significativamente. O uso de 100% de reciclagem permite uma poupança energética de cerca de 98% no Inverno e na Meia-Estação e de cerca de 93% no Verão. Caso não fosse realizada reciclagem, seria necessário fornecer à corrente de ar cerca de 18,81 kW para elevar a sua temperatura de 20ºC para 50ºC (Meia-Estação), cerca de 24,67 kW para elevar a sua temperatura de 10ºC para 50ºC (Inverno) e na ordem dos 8,90 kW para elevar a sua temperatura dos 35ºC para 50ºC (Verão). No caso do transporte pneumático, existem duas linhas, uma horizontal e uma vertical, logo foi necessário estimar o valor da velocidade das partículas para estes dois casos. Na linha vertical, a velocidade da partícula é cerca de 25,03 m/s e cerca de 35,95 m/s na linha horizontal. O menor valor para a linha vertical prende-se com o facto de nesta zona ter que se vencer a força gravítica. Em ambos os circuitos a velocidade do fluido é cerca de 47,17 m/s. No interior da coluna, a velocidade do fluido tem o valor de 10,90 m/s e a velocidade das partículas é de 1,04 m/s. A queda de pressão total no sistema é cerca de 2408 Pa. A análise de custos ao sistema de secagem indicou que este sistema irá acarretar um custo total (fabrico mais transporte) de cerca de 153035€. Este sistema necessita de electricidade para funcionar, e esta irá acarretar um custo anual de cerca de 7951,4€. Embora este sistema de secagem apresente a possibilidade de se realizar uma razão de reciclagem na ordem dos 100% e também seja possível adaptar o mesmo para diferentes tipos de cereais, e até outros tipos de materiais, desde que possam ser fluidizados, o seu custo impede que a realização deste investimento não seja atractiva, especialmente tendo em consideração que se trata de uma instalação à escala piloto com uma capacidade de 45 kgs.

Fate Hazardous elements in soils surrounding a coal-fired power plant complex

Proceedings of the 13th International UFZ-Deltares Conference on Sustainable Use and Management of Soil, Sediment and Water Resources - 9–12 June 2015 • Copenhagen, Denmark

Estudo sobre viabilidade de produção descentralizada de energia à custa do efeito chaminé

A procura por alternativas ao atual paradigma energético, que se caracteriza por uma predominância indiscutível das fontes combustíveis fósseis, é o motivo primário desta investigação. A energia emitida pelo Sol que chega à Terra diariamente ultrapassa em várias ordens de grandeza a energia que a nossa sociedade atual necessita. O efeito chaminé é uma das formas de aproveitar essa energia. Este efeito tem origem no diferencial de temperaturas existente entre o interior e o exterior de uma chaminé, que provoca um gradiente nas massas volúmicas do fluido entre o interior e o exterior da chaminé, induzindo assim um fluxo de ar. Esta diferença de temperaturas radica na exposição da face exterior da chaminé à radiação solar. No sistema que nos propomos estudar, o ar entra na chaminé por pequenos orifícios situados na sua base, e, ao tomar contacto com as paredes internas da chaminé, aquece desde a temperatura ambiente, Ta, até à temperatura interna, Ti . Este aumento de temperatura torna o ar dentro da chaminé mais “leve” em comparação com o ar mais frio do exterior levando-o a ascender ao longo do interior da chaminé. Este escoamento contém energia cinética que pode, por exemplo, ser transformada em energia elétrica por intermédio de turbinas. A eficiência de conversão da energia será tanto maior quanto menor for a velocidade do ar a jusante da turbina. Esta tecnologia poderá ser instalada de forma descentralizada, como acontece com as atuais centrais concentradoras solares térmicas e fotovoltaicas localizadas na periferia de grandes cidades ou, alternativamente, poderá ser inserida no próprio tecido urbanístico. A investigação demonstra que as dimensões da chaminé, a irradiação e a temperatura do ar são os fatores com maior impacto na potência hidráulica gerada.