Avaliação do desempenho energético e da qualidade do ar de uma piscina coberta

Este trabalho teve como propósito fazer uma avaliação do desempenho energético e da qualidade do ar no interior das instalações de uma Piscina Municipal Coberta, localizada na zona norte de Portugal, sendo estabelecidos os seguintes objetivos: caracterização geral da piscina, no que respeita aos seus diferentes espaços e equipamentos, cálculo dos consumos térmicos e elétricos bem como o registo das concentrações de elementos poluentes para controlo da qualidade do ar no interior da piscina, tendo como base a legislação atualmente em vigor. A caracterização geral da piscina permitiu verificar algumas inconformidades como a temperatura da água nos tanques de natação que tem valores superiores aos recomendados e a sala de primeiros socorros que não possui acesso direto ao exterior. Acrescente-se que o pavimento nos chuveiros da casa de banho feminina e os valores de pH para água do tanque grande e pequeno não estão sempre dentro da gama de recomendação. O caudal da renovação de ar está a ser operado manualmente e quando está a funcionar a 50% da sua capacidade máxima, que acontece numa parte do dia, apenas consegue renovar 77,5% do caudal recomendado pelo RSECE. Para se obter o valor recomendado é necessário ter pelo menos 7 horas com o caudal a 100% da capacidade máxima. A avaria na UTA2 originou que 40% dos registos diários da humidade relativa interior estivessem fora da gama de valores recomendados e que esta é fortemente dependente da humidade no exterior e pode ser agravada quando as portas dos envidraçados da nave são abertas. Analisando ainda a quantidade de água removida na desumidificação do ar com a água evaporada em condições de Outono-Inverno ou Primavera-Verão, este estudo permitiu concluir que todas as combinações demonstraram a necessidade de desumidificação salvo a combinação Outono-Inverno e UTA2 a funcionar a 100% da sua capacidade máxima. Os isolamentos das tubagens na sala das caldeiras foram observados e comparados com as soluções recomendadas pelas empresas especialistas e verificou-se que alguns estão mal colocados com parcial ou total degradação, promovendo perdas térmicas. No caso das perdas calorificas por evaporação, estas representaram cerca de 67,78% das perdas totais. Como tal, estudou-se a aplicação de uma cobertura sobre o plano de água durante o período de inatividade da piscina (8 horas) e verificou-se que o resultado seria uma poupança de 654,8 kWh/dia, na ausência de evaporação da água, mais 88,00 kWh/dia do período da UTA2 a funcionar a 50% da sua capacidade, perfazendo um total de 742,8 kWh/dia. A aplicação da cobertura permite obter um VAL de valor positivo, uma TIR de 22,77% e sendo este valor superior ao WACC (Weight Average Cost of Capital), o projeto torna-se viável com um Pay-Back de 3,17 anos. Caracterizou-se também o consumo total diário em eletricidade, e verificou-se que as unidades de climatização, as bombas de circulação de água, a iluminação, e outros equipamentos representam, respetivamente, cerca de 67,81, 25,26, 2,68 e 3,91% da energia elétrica total consumida. Por fim, a análise à qualidade do ar no interior da nave em Maio e Setembro identificou que as concentrações de ozono apresentavam valores no limite do aceitável em Maio e superiores ao valor de emissão em Setembro. Os compostos orgânicos voláteis também apresentavam valores em Maio 4,98 vezes superior e em Setembro 6,87 vezes superior aos valores máximos exigidos pelo D.L. nº 79/2006. Houve ainda altas concentrações de radão registadas na casa dos filtros, em Maio com um valor 11,49 vezes superior, no entanto esse valor desceu em Setembro para 1,08 vezes, mesmo assim superior ao exigido pelo D.L. nº 79/2006.

Transferência de calor em fluidos Não-Newtonianos

Foi realizado um estudo sobre o efeito do tipo de fluidos na transferência de calor. Pretende-se determinar a influência da concentração da solução de Goma de Xantano, do número de Reynolds, do número de Weissenberg, da temperatura e do tempo de escoamento no coeficiente de transferência de calor, jH. O estudo da transferência de calor foi feito num permutador de tubos duplos concêntricos. Já o estudo da reologia foi realizado num reómetro. Na caracterização reológica das soluções de XG, a viscosidade aumenta com a concentração das soluções, diminui para taxas de deformação crescentes e com o aumento da temperatura para ambas as soluções. Os dados mostram um aumento da intensidade da pseudoplasticidade com a concentração do polímero, sendo os valores representados pelo modelo de Sisco. A degradação da solução de 0,20% de goma de xantano a 25 ºC, com o escoamento, é muito acelerada. Os resultados dos ensaios apresentam uma diminuição da viscosidade de 9,4% a 22,9%, para tempos de escoamento de 12 a 47 horas, respectivamente. Num escoamento turbulento em conduta de secção circular constante os resultados mostram uma redução de arrasto total de 18 para 33%. Para a solução de 0,10 % de XG, verifica-se um aumento do calor transferido de 115% e de 130%, quando a temperatura aumenta de 25 ºC para 36 ºC, respectivamente. A água apresenta valores de calor transferido superiores, cerca de 170%, aos da solução de 0,1 %XG. O factor de correlação empírico de Colbourn (jH), utilizado neste trabalho apresenta valores de acordo com a relação de Cho and Hartnett (1985): jH<2/f. Quando o caudal do fluido quente aumenta verifica-se uma diminuição do factor jH. Em relação ao tempo de escoamento verifica-se uma diminuição de cerca de 70% do coeficiente de transferência de calor ao fim de 47 horas. Finalmente verificamos uma diminuição do factor de transferência de calor com o aumento da temperatura do fluido quente, para ambas as concentrações de goma de xantano. Para as soluções de 0,10 e 0,20% de XG essa diminuição variou entre 38 e 15% e entre 34 e 3%, respectivamente.

Modelação do transporte e dispersão de um marcador corado no Rio Febros

O rio Febros é um pequeno curso de água, situado no concelho de Vila Nova de Gaia, com cerca de 15 km de extensão, cuja bacia hidrográfica ocupa uma área de aproximadamente 35,4 km2. Nasce em Seixezelo e desagua na margem esquerda do Rio Douro no Cais do Esteiro, em Avintes. Em Maio de 2008, um acidente de viação teve como consequência o derrame de cerca de quatro toneladas de ácido clorídrico que rapidamente convergiu às águas do rio. Apenas um dia depois, o pH desceu para três e muitos foram os peixes que morreram. A solução adoptada para evitar o desaire foi introduzir milhares de litros de água de modo a diluir o ácido presente, ao longo de todo o curso de água. Tal facto não evitou a destruição de parte de um ecossistema, que ainda nos dias de hoje se encontra em recuperação. De forma a avaliar-se o impacto destas possíveis perturbações sejam estas de origem antropogénica ou natural é necessário possuir conhecimentos dos processos químicos tais como a advecção, a mistura devida à dispersão e a transferência de massa ar/água. Estes processos irão determinar o movimento e destino das substâncias que poderão ser descarregadas no rio. Para tal, recorrer-se-á ao estudo hidrogeométrico do curso de água assim como ao estudo do comportamento de um marcador, simulando uma possível descarga. A rodamina WT será o marcador a ser utilizado devido à panóplia de características ambientalmente favoráveis. Os estudos de campo com este corante, realizados em sequência de descarga previamente estudada, fornecem uma das melhores fontes de informação para verificação e validação de modelos hidráulicos utilizados em estudos de qualidade de águas e protecção ambiental. Escolheram-se dois pontos de descarga no Febros, um em Casal Drijo e outro no Parque Biológico de Gaia, possuindo cada um deles, a jusante, duas estações de monitorização. Pelo modelo ADE os valores obtidos para o coeficiente de dispersão longitudinal para as estações Pontão d’ Alheira, Pinheiral, Menesas e Giestas foram, respectivamente, 0,3622; 0,5468; 1,6832 e 1,7504 m2/s. Para a mesma sequência de estações, os valores da velocidade de escoamento obtidos neste trabalho experimental foram de 0,0633; 0,0684; 0,1548 e 0,1645 m/s. Quanto ao modelo TS, os valores obtidos para o coeficiente de dispersão longitudinal para as estações Pontão d’ Alheira, Pinheiral, Menesas e Giestas foram, respectivamente, 0,2339; 0,1618; 0,5057e 1,1320 m2/s. Para a mesma sequência de estações, os valores da velocidade de escoamento obtidos neste trabalho experimental foram de 0,0652; 0,0775; 0,1891 e 0,1676 m/s. Os resultados foram ajustados por um método directo, o método dos momentos, e por dois métodos indirectos, os modelos ADE e TS. O melhor ajuste corresponde ao modelo TS onde os valores do coeficiente de dispersão longitudinal e da velocidade de escoamento são aqueles que melhor se aproximam da realidade. Quanto ao método dos momentos, o valor estimado para a velocidade é de 0,162 m/s e para o coeficiente de dispersão longitudinal de 9,769 m2/s. Não obstante, a compreensão da hidrodinâmica do rio e das suas características, bem como a adequação de modelos matemáticos no tratamento de resultados formam uma estratégia de protecção ambiental inerente a futuros impactos que possam suceder.

Soil remediation time to achieve clean-up goals I: influence of soil water content

The current models are not simple enough to allow a quick estimation of the remediation time. This work reports the development of an easy and relatively rapid procedure for the forecasting of the remediation time using vapour extraction. Sandy soils contaminated with cyclohexane and prepared with different water contents were studied. The remediation times estimated through the mathematical fitting of experimental results were compared with those of real soils. The main objectives were: (i) to predict, through a simple mathematical fitting, the remediation time of soils with water contents different from those used in the experiments; (ii) to analyse the influence of soil water content on the: (ii1) remediation time; (ii2) remediation efficiency; and (ii3) distribution of contaminants in the different phases present into the soil matrix after the remediation process. For sandy soils with negligible contents of clay and natural organic matter, artificially contaminated with cyclohexane before vapour extraction, it was concluded that (i) if the soil water content belonged to the range considered in the experiments with the prepared soils, then the remediation time of real soils of similar characteristics could be successfully predicted, with relative differences not higher than 10%, through a simple mathematical fitting of experimental results; (ii) increasing soil water content from 0% to 6% had the following consequences: (ii1) increased remediation time (1.8–4.9 h, respectively); (ii2) decreased remediation efficiency (99–97%, respectively); and (ii3) decreased the amount of contaminant adsorbed onto the soil and in the non-aqueous liquid phase, thus increasing the amount of contaminant in the aqueous and gaseous phases.

Soil remediation time to achieve clean-up goals II: influence of natural organic matter and water contents

This work reports a relatively rapid procedure for the forecasting of the remediation time (RT) of sandy soils contaminated with cyclohexane using vapour extraction. The RT estimated through the mathematical fitting of experimental results was compared with that of real soils. The main objectives were: (i) to predict the RT of soils with natural organic matter (NOM) and water contents different from those used in experiments; and (ii) to analyse the time and efficiency of remediation, and the distribution of contaminants into the soil matrix after the remediation process, according to the soil contents of: (ii1) NOM; and (ii2) water. For sandy soils with negligible clay contents, artificially contaminated with cyclohexane before vapour extraction, it was concluded that: (i) if the NOM and water contents belonged to the range of the prepared soils, the RT of real soils could be predicted with relative differences not higher than 12%; (ii1) the increase of NOM content from 0% to 7.5% increased the RT (1.8–13 h) and decreased the remediation efficiency (RE) (99–90%) and (ii2) the increase of soil water content from 0% to 6% increased the RT (1.8–4.9 h) and decreased the RE (99–97%). NOM increases the monolayer capacity leading to a higher sorption into the solid phase. Increasing of soil water content reduces the mass transfer coefficient between phases. Concluding, NOM and water contents influence negatively the remediation process, turning it less efficient and more time consuming, and consequently more expensive.

Study of vanadium doped ZnO films prepared by dc reactive magnetron sputtering at different substrate temperatures

ZnO films doped with vanadium (ZnO:V) have been prepared by dc reactive magnetron sputtering technique at different substrate temperatures (RT–500 C). The effects of the substrate temperature on ZnO:V films properties have been studied. XRD measurements show that only ZnO polycrystalline structure has been obtained, no V2O5 or VO2 crystal phase can be observed. It has been found that the film prepared at low substrate temperature has a preferred orientation along the (002) direction. As the substrate temperature is increased, the (002) peak intensity decreases. When the substrate temperature reaches the 500 C, the film shows a random orientation. SEM measurements show a clear formation of the nano-grains in the sample surface when the substrate temperature is higher than 400 C. The optical properties of the films have been studied by measuring the specular transmittance. The refractive index has been calculated by fitting the transmittance spectra using OJL model combined with harmonic oscillator.

Análise de tecnologias de impermeabilização e isolamento em reabilitação da envolvente de edifícios

Neste trabalho é abordado o estágio efectuado por um período de sete meses na E.A. (Engenheiros Associados) e um estudo de caso sobre as tecnologias aplicadas no isolamento térmico de paredes exteriores. Na primeira parte deste relatório é efectuada uma breve caracterização da empresa e da sua actividade no mercado. A Engenheiros e Associados tem a sua estrutura assente nos seguintes sectores: Técnico-Comercial, de Aprovisionamento e de Gestão Administrativa; da qual foi efectuada uma descrição do trabalho desenvolvido em cada um desses sectores, nomeadamente a permanência na obra e as visitas e diligências às obras. A segunda parte deste relatório, que é a parte fulcral do trabalho desenvolvido, assenta sobre o sistema de ETICS (External Thermal Insulation Composite System), ou seja, o sistema de reboco delgado armado sobre isolamento térmico. Foi utilizado o sistema Cappotto na Obra de Requalificação da Urbanização de Vila d’Este – 1ª Fase. Após uma breve apresentação e evolução da história do ETICS ao longo do tempo, é exposta a caracterização da obra onde vai ser aplicado o sistema, que recai sobre a tipologia dos edifícios. Refira-se que, neste relatório, os edifícios por serem de construção túnel se dividem em três tipos: tipo plana ou corrente, tipo ângulo e tipo topo. Na Analise de Patologias efectuada evidenciam-se as fissuras e a humidade presentes em quase toda a extensão dos referidos edifícios antes do tratamento, tornando-as assim as principais anomalias destes edifícios. Foram elemento de foco as anomalias existentes nos edifícios, como: Deterioração do fibrocimento; Insuficiência das caleiras; Deficiências das impermeabilizações; Ausência de rufos; Deficiência das ligações; Degradação do revestimento e pintura; Fissuração do reboco; Degradação dos forros exteriores; Deterioração das padieiras; Deterioração dos peitoris; Deterioração das juntas de dilatação; Infiltrações e condensações; Ruptura das canalizações; Deslocamentos; Deficiências de ventilação; Deficiências de estanqueidade; Inexistência de Sistema de combate a incêndios. Para estas patologias são apresentadas as propostas de solução de forma a eliminar as mesmas. Sendo o ETICS escolhido por ser o sistema que elimina a maior parte destas patologias, tendo em conta uma relação de qualidade/preço, é abordada de uma forma detalhada e extensiva da aplicação do sistema antes, durante e depois da obra em si. Assim como, é feita uma descrição pormenorizada do material utilizado para a implementação do sistema. A análise dos pontos críticos refere-se a zonas sensíveis onde há a necessidade de reforço do sistema com vista a eliminar o aparecimento posterior de patologias, como por exemplo as características de suporte. Após a aplicação do sistema podem aparecer algumas patologias das quais se destaca o facto dos produtos serem preparados em obra, o erro humano nas dosagens do fabricante e no acrescento de água sem necessidade e a par das condições climatéricas, são as causas mais comuns do aparecimento de anomalias no sistema de ETICS, provocando fissurações e infiltrações, que são descritas neste relatório. É também abordada a manutenção e reparação do sistema, onde a manutenção refere-se à lavagem e à remoção de microorganismos das paredes e posterior pintura. A reparação dos danos divide-se em dois tipos, áreas até 2 cm2 e áreas maiores que 2 cm2. Por fim, são apresentados alguns rendimentos, que foram possíveis obter ao longo do desenvolvimento do trabalho, dos materiais e da mão-de-obra, dando origem aos custos directos. Sendo também abordadas as vantagens e desvantagens do sistema desde o seu início até à sua conclusão.

Power law analysis of financial index dynamics

Power law PL and fractional calculus are two faces of phenomena with long memory behavior. This paper applies PL description to analyze different periods of the business cycle. With such purpose the evolution of ten important stock market indices DAX, Dow Jones, NASDAQ, Nikkei, NYSE, S&P500, SSEC, HSI, TWII, and BSE over time is studied. An evolutionary algorithm is used for the fitting of the PL parameters. It is observed that the PL curve fitting constitutes a good tool for revealing the signal main characteristics leading to the emergence of the global financial dynamic evolution.

Self-similarity principle: the reduced description of randomness

A new general fitting method based on the Self-Similar (SS) organization of random sequences is presented. The proposed analytical function helps to fit the response of many complex systems when their recorded data form a self-similar curve. The verified SS principle opens new possibilities for the fitting of economical, meteorological and other complex data when the mathematical model is absent but the reduced description in terms of some universal set of the fitting parameters is necessary. This fitting function is verified on economical (price of a commodity versus time) and weather (the Earth’s mean temperature surface data versus time) and for these nontrivial cases it becomes possible to receive a very good fit of initial data set. The general conditions of application of this fitting method describing the response of many complex systems and the forecast possibilities are discussed.

Complex-order forced van der Pol oscillator

In this paper we consider a complex-order forced van der Pol oscillator. The complex derivative Dα1jβ, with α, β ∈ ℝ+, is a generalization of the concept of an integer derivative, where α = 1, β = 0. The Fourier transforms of the periodic solutions of the complex-order forced van der Pol oscillator are computed for various values of parameters such as frequency ω and amplitude b of the external forcing, the damping μ, and parameters α and β. Moreover, we consider two cases: (i) b = 1, μ = {1.0, 5.0, 10.0}, and ω = {0.5, 2.46, 5.0, 20.0}; (ii) ω = 20.0, μ = {1.0, 5.0, 10.0}, and b = {1.0, 5.0, 10.0}. We verified that most of the signal energy is concentrated in the fundamental harmonic ω0. We also observed that the fundamental frequency of the oscillations ω0 varies with α and μ. For the range of tested values, the numerical fitting led to logarithmic approximations for system (7) in the two cases (i) and (ii). In conclusion, we verify that by varying the parameter values α and β of the complex-order derivative in expression (7), we accomplished a very effective way of perturbing the dynamical behavior of the forced van der Pol oscillator, which is no longer limited to parameters b and ω.

Sistema de detecção e localização de veículos subaquáticos

Mestrado em Engenharia Electrotécnica e de Computadores.Área de Especialização de Sistemas Autónomos

Modelação e identificação de motor DC

O modelo matemático de um sistema real permite o conhecimento do seu comportamento dinâmico e é geralmente utilizado em problemas de engenharia. Por vezes os parâmetros utilizados pelo modelo são desconhecidos ou imprecisos. O envelhecimento e o desgaste do material são fatores a ter em conta pois podem causar alterações no comportamento do sistema real, podendo ser necessário efetuar uma nova estimação dos seus parâmetros. Para resolver este problema é utilizado o software desenvolvido pela empresa MathWorks, nomeadamente, o Matlab e o Simulink, em conjunto com a plataforma Arduíno cujo Hardware é open-source. A partir de dados obtidos do sistema real será aplicado um Ajuste de curvas (Curve Fitting) pelo Método dos Mínimos Quadrados de forma a aproximar o modelo simulado ao modelo do sistema real. O sistema desenvolvido permite a obtenção de novos valores dos parâmetros, de uma forma simples e eficaz, com vista a uma melhor aproximação do sistema real em estudo. A solução encontrada é validada com recurso a diferentes sinais de entrada aplicados ao sistema e os seus resultados comparados com os resultados do novo modelo obtido. O desempenho da solução encontrada é avaliado através do método das somas quadráticas dos erros entre resultados obtidos através de simulação e resultados obtidos experimentalmente do sistema real.

Estudo da turvação em rede de abastecimento de água a alta pressão no grande porto

A turvação é uma das principais preocupações das empresas de rede de distribuição de água potável. Para se saber a sua origem existem diversos mecanismos, sugeridos em bibliografia, que explicam a formação de depósitos de partículas nas tubagens, que se acumulam em camadas sobrepostas, criando entre si forças de coesão. A turvação pode surgir se o caudal ou a velocidade na tubagem aumentar, por alguma razão, pois um aumento da tensão de corte sobre o depósito, pode superar as forças de coesão levando à sua ressuspensão. Esta explicação é traduzida matematicamente pelo modelo PODDS (Prediction of Discolouration in Drinking Water Distribtuion System). O modelo PODDS foi compilado em EPANET recentemente por uma equipa de investigadores da Universidade de Sheffield e foi utilizado para fazer o estudo num ramal da empresa AdDP - Águas do Douro e Paiva, SA. Os resultados das simulações permitirão prever o caudal mínimo teórico a utilizar no ramal que origina o valor limite de turvação. Foram realizados ensaios de variação do caudal no ramal de Jovim-Ramalde para a calibração do modelo. Os parâmetros obtidos foram: Cmáx igual a 100 NTUm, b e n igual a 1, k é igual -1 NTUm3/N, e por fim, P igual a 0,0001 NTUm3/N.s. Após calibração do modelo realizaram-se simulações com atualização das forças de coesão para conhecer o caudal que provocaria a mobilização das partículas até uma turvação de 4 NTU e concluiu-se que o valor seria de 4,32 m3/s. O tempo necessário para regressar a valores inferiores a 1 NTU também foi determinado e 4 horas e 8 minutos foi o valor alcançado. Por fim realizou-se um estudo da influência da temperatura neste modelo e o resultado mostrou que este parâmetro é pouco influente para a gama de temperaturas considerada.

Previsão do Diagrama de Carga de Subestações da REN utilizando Redes Neuronais

A presente dissertação apresenta o estudo de previsão do diagrama de carga de subestações da Rede Elétrica Nacional (REN) utilizando redes neuronais, com o intuito de verificar a viabilidade do método utilizado, em estudos futuros. Atualmente, a energia elétrica é um bem essencial e desempenha um papel fundamental, tanto a nível económico do país, como a nível de conforto e satisfação individual. Com o desenvolvimento do setor elétrico e o aumento dos produtores torna-se importante a realização da previsão de diagramas de carga, contribuindo para a eficiência das empresas. Esta dissertação tem como objetivo a utilização do modelo das redes neuronais artificiais (RNA) para criar uma rede capaz de realizar a previsão de diagramas de carga, com a finalidade de oferecer a possibilidade de redução de custos e gastos, e a melhoria de qualidade e fiabilidade. Ao longo do trabalho são utilizados dados da carga (em MW), obtidos da REN, da subestação da Prelada e dados como a temperatura, humidade, vento e luminosidade, entre outros. Os dados foram devidamente tratados com a ajuda do software Excel. Com o software MATLAB são realizados treinos com redes neuronais, através da ferramenta Neural Network Fitting Tool, com o objetivo de obter uma rede que forneça os melhores resultados e posteriormente utiliza-la na previsão de novos dados. No processo de previsão, utilizando dados reais das subestações da Prelada e Ermesinde referentes a Março de 2015, comprova-se que com a utilização de RNA é possível obter dados de previsão credíveis, apesar de não ser uma previsão exata. Deste modo, no que diz respeito à previsão de diagramas de carga, as RNA são um bom método a utilizar, uma vez que fornecem, à parte interessada, uma boa previsão do consumo e comportamento das cargas elétricas. Com a finalização deste estudo os resultados obtidos são no mínimo satisfatórios. Consegue-se alcançar através das RNA resultados próximos aos valores que eram esperados, embora não exatamente iguais devido à existência de uma margem de erro na aprendizagem da rede neuronal.

Effects of Soil Compaction and Organic Carbon Content on Preferential Flow in Loamy Field Soils

Abstract: Preferential flow and transport through macropores affect plant water use efficiency and enhance leaching of agrochemicals and the transport of colloids, thereby increasing the risk for contamination of groundwater resources. The effects of soil compaction, expressed in terms of bulk density (BD), and organic carbon (OC) content on preferential flow and transport were investigated using 150 undisturbed soil cores sampled from 15 × 15–m grids on two field sites. Both fields had loamy textures, but one site had significantly higher OC content. Leaching experiments were conducted in each core by applying a constant irrigation rate of 10 mm h−1 with a pulse application of tritium tracer. Five percent tritium mass arrival times and apparent dispersivities were derived from each of the tracer breakthrough curves and correlated with texture, OC content, and BD to assess the spatial distribution of preferential flow and transport across the investigated fields. Soils from both fields showed strong positive correlations between BD and preferential flow. Interestingly, the relationships between BD and tracer transport characteristics were markedly different for the two fields, although the relationship between BD and macroporosity was nearly identical. The difference was likely caused by the higher contents of fines and OC at one of the fields leading to stronger aggregation, smaller matrix permeability, and a more pronounced pipe-like pore system with well-aligned macropores.