Diffusion coefficients in dense and supercritical fluids

Os coeficientes de difusão (D 12) são propriedades fundamentais na investigação e na indústria, mas a falta de dados experimentais e a inexistência de equações que os estimem com precisão e confiança em fases comprimidas ou condensadas constituem limitações importantes. Os objetivos principais deste trabalho compreendem: i) a compilação de uma grande base de dados para valores de D 12 de sistemas gasosos, líquidos e supercríticos; ii) o desenvolvimento e validação de novos modelos de coeficientes de difusão a diluição infinita, aplicáveis em amplas gamas de temperatura e densidade, para sistemas contendo componentes muito distintos em termos de polaridade, tamanho e simetria; iii) a montagem e teste de uma instalação experimental para medir coeficientes de difusão em líquidos e fluidos supercríticos. Relativamente à modelação, uma nova expressão para coeficientes de difusão a diluição infinita de esferas rígidas foi desenvolvida e validada usando dados de dinâmica molecular (desvio relativo absoluto médio, AARD = 4.44%) Foram também estudados os coeficientes de difusão binários de sistemas reais. Para tal, foi compilada uma extensa base de dados de difusividades de sistemas reais em gases e solventes densos (622 sistemas binários num total de 9407 pontos experimentais e 358 moléculas) e a mesma foi usada na validação dos novos modelos desenvolvidos nesta tese. Um conjunto de novos modelos foi proposto para o cálculo de coeficientes de difusão a diluição infinita usando diferentes abordagens: i) dois modelos de base molecular com um parâmetro específico para cada sistema, aplicáveis em sistemas gasosos, líquidos e supercríticos, em que natureza do solvente se encontra limitada a apolar ou fracamente polar (AARDs globais na gama 4.26-4.40%); ii) dois modelos de base molecular biparamétricos, aplicáveis em todos os estados físicos, para qualquer tipo de soluto diluído em qualquer solvente (apolar, fracamente polar e polar). Ambos os modelos dão origem a erros globais entre 2.74% e 3.65%; iii) uma correlação com um parâmetro, específica para coeficientes de difusão em dióxido de carbono supercrítico (SC-CO2) e água líquida (AARD = 3.56%); iv) nove correlações empíricas e semi-empíricas que envolvem dois parâmetros, dependentes apenas da temperatura e/ou densidade do solvente e/ou viscosidade do solvente. Estes últimos modelos são muito simples e exibem excelentes resultados (AARDs entre 2.78% e 4.44%) em sistemas líquidos e supercríticos; e v) duas equações preditivas para difusividades de solutos em SC-CO2, em que os erros globais de ambas são inferiores a 6.80%. No global, deve realçar-se o facto de os novos modelos abrangerem a grande variedade de sistemas e moléculas geralmente encontrados. Os resultados obtidos são consistentemente melhores do que os obtidos com os modelos e abordagens encontrados na literatura. No caso das correlações com um ou dois parâmetros, mostrou-se que estes mesmos parâmetros podem ser ajustados usando um conjunto muito pequeno de dados, e posteriormente serem utilizados na previsão de valores de D 12 longe do conjunto original de pontos. Uma nova instalação experimental para medir coeficientes de difusão binários por técnicas cromatográficas foi montada e testada. O equipamento, o procedimento experimental e os cálculos analíticos necessários à obtenção dos valores de D 12 pelo método de abertura do pico cromatográfico, foram avaliados através da medição de difusividades de tolueno e acetona em SC-CO2. Seguidamente, foram medidos coeficientes de difusão de eucaliptol em SC-CO2 nas gamas de 202 – 252 bar e 313.15 – 333.15 K. Os resultados experimentais foram analisados através de correlações e modelos preditivos para D12.

Novel critical phenomena in optimization driven processes on networks

The work presented in this Ph.D thesis was developed in the context of complex network theory, from a statistical physics standpoint. We examine two distinct problems in this research field, taking a special interest in their respective critical properties. In both cases, the emergence of criticality is driven by a local optimization dynamics. Firstly, a recently introduced class of percolation problems that attracted a significant amount of attention from the scientific community, and was quickly followed up by an abundance of other works. Percolation transitions were believed to be continuous, until, recently, an 'explosive' percolation problem was reported to undergo a discontinuous transition, in [93]. The system's evolution is driven by a metropolis-like algorithm, apparently producing a discontinuous jump on the giant component's size at the percolation threshold. This finding was subsequently supported by number of other experimental studies [96, 97, 98, 99, 100, 101]. However, in [1] we have proved that the explosive percolation transition is actually continuous. The discontinuity which was observed in the evolution of the giant component's relative size is explained by the unusual smallness of the corresponding critical exponent, combined with the finiteness of the systems considered in experiments. Therefore, the size of the jump vanishes as the system's size goes to infinity. Additionally, we provide the complete theoretical description of the critical properties for a generalized version of the explosive percolation model [2], as well as a method [3] for a precise calculation of percolation's critical properties from numerical data (useful when exact results are not available). Secondly, we study a network flow optimization model, where the dynamics consists of consecutive mergings and splittings of currents flowing in the network. The current conservation constraint does not impose any particular criterion for the split of current among channels outgoing nodes, allowing us to introduce an asymmetrical rule, observed in several real systems. We solved analytically the dynamic equations describing this model in the high and low current regimes. The solutions found are compared with numerical results, for the two regimes, showing an excellent agreement. Surprisingly, in the low current regime, this model exhibits some features usually associated with continuous phase transitions.

Multiplicity results for some classes of Schrödinger-Poisson systems

In this thesis, we study the existence and multiplicity of solutions of the following class of Schr odinger-Poisson systems: u + u + l(x) u = (x; u) in R3; = l(x)u2 in R3; where l 2 L2(R3) or l 2 L1(R3). And we consider that the nonlinearity satis es the following three kinds of cases: (i) a subcritical exponent with (x; u) = k(x)jujp 2u + h(x)u (4 p < 2 ) under an inde nite case; (ii) a general inde nite nonlinearity with (x; u) = k(x)g(u) + h(x)u; (iii) a critical growth exponent with (x; u) = k(x)juj2 2u + h(x)jujq 2u (2 q < 2 ). It is worth mentioning that the thesis contains three main innovations except overcoming several di culties, which are generated by the systems themselves. First, as an unknown referee said in his report, we are the rst authors concerning the existence of multiple positive solutions for Schr odinger- Poisson systems with an inde nite nonlinearity. Second, we nd an interesting phenomenon in Chapter 2 and Chapter 3 that we do not need the condition R R3 k(x)ep 1dx < 0 with an inde nite noncoercive case, where e1 is the rst eigenfunction of +id in H1(R3) with weight function h. A similar condition has been shown to be a su cient and necessary condition to the existence of positive solutions for semilinear elliptic equations with inde nite nonlinearity for a bounded domain (see e.g. Alama-Tarantello, Calc. Var. PDE 1 (1993), 439{475), or to be a su cient condition to the existence of positive solutions for semilinear elliptic equations with inde nite nonlinearity in RN (see e.g. Costa-Tehrani, Calc. Var. PDE 13 (2001), 159{189). Moreover, the process used in this case can be applied to study other aspects of the Schr odinger-Poisson systems and it gives a way to study the Kirchho system and quasilinear Schr odinger system. Finally, to get sign changing solutions in Chapter 5, we follow the spirit of Hirano-Shioji, Proc. Roy. Soc. Edinburgh Sect. A 137 (2007), 333, but the procedure is simpler than that they have proposed in their paper.

Calculus of variations on time scales and applications to economics

We consider some problems of the calculus of variations on time scales. On the beginning our attention is paid on two inverse extremal problems on arbitrary time scales. Firstly, using the Euler-Lagrange equation and the strengthened Legendre condition, we derive a general form for a variation functional that attains a local minimum at a given point of the vector space. Furthermore, we prove a necessary condition for a dynamic integro-differential equation to be an Euler-Lagrange equation. New and interesting results for the discrete and quantum calculus are obtained as particular cases. Afterwards, we prove Euler-Lagrange type equations and transversality conditions for generalized infinite horizon problems. Next we investigate the composition of a certain scalar function with delta and nabla integrals of a vector valued field. Euler-Lagrange equations in integral form, transversality conditions, and necessary optimality conditions for isoperimetric problems, on an arbitrary time scale, are proved. In the end, two main issues of application of time scales in economic, with interesting results, are presented. In the former case we consider a firm that wants to program its production and investment policies to reach a given production rate and to maximize its future market competitiveness. The model which describes firm activities is studied in two different ways: using classical discretizations; and applying discrete versions of our result on time scales. In the end we compare the cost functional values obtained from those two approaches. The latter problem is more complex and relates to rate of inflation, p, and rate of unemployment, u, which inflict a social loss. Using known relations between p, u, and the expected rate of inflation π, we rewrite the social loss function as a function of π. We present this model in the time scale framework and find an optimal path π that minimizes the total social loss over a given time interval.

Um modelo de medição e gestão de desempenho para as universidades públicas

Nas últimas décadas as instituições de ensino superior têm sido alvo de uma pressão crescente para aumentar a sua eficiência e a sua eficácia. Fatores como a diversificação da oferta formativa, a massificação, o aumento da internacionalização, entre outros, conduziram a uma maior monitorização das instituições e, por sua vez, geraram o desenvolvimento de novos sistemas de informação. A constante solicitação de informação, quer pelo Estado, quer pelo mercado constitui a base para a definição do objetivo desta investigação: construir um modelo integrado de medição e gestão de desempenho para as universidades públicas e testá-lo no universo português. Para a construção conceptual deste modelo foi realizada uma revisão da literatura baseada em diferentes contextos: organizacional, europeu e nacional. Este modelo foi depois aplicado às universidades públicas portuguesas recorrendo a dados disponibilizados por entidades públicas, tendo em conta a área de educação e formação (CNAEF). Uma análise descritiva aos dados constituiu uma contribuição desta investigação para a prática, no sentido de que permitiu fornecer um conjunto de recomendações às universidades e às entidades oficiais responsáveis pela recolha de dados a nível nacional, relativamente aos sistemas de informação e processos de recolha de dados. O modelo proposto constitui a contribuição teórica desta investigação, no sentido de que integra, no mesmo modelo, as diferentes dimensões de desempenho, a visão dos diferentes stakeholders, quer no contexto do ensino, quer no da investigação, quer no da própria gestão da universidade. A vertente analítica deste modelo, representada pelas diferentes relações entre os grupos de indicadores, foi testada recorrendo à técnica de análise de equações estruturais baseada na variância (nomeadamente Partial Least Squares) em quatro áreas CNAEF. Os resultados demonstraram que o comportamento em termos de medição e gestão de desempenho difere consoante a CNAEF, identificando um maior ajustamento às áreas hard e evidenciando que as áreas soft necessitam de indicadores mais ajustados às suas especificidades.

Characterization of RSOA modulators in radio over fiber links

In this work physical and behavioral models for a bulk Reflective Semiconductor Optical Amplifier (RSOA) modulator in Radio over Fiber (RoF) links are proposed. The transmission performance of the RSOA modulator is predicted under broadband signal drive. At first, the simplified physical model for the RSOA modulator in RoF links is proposed, which is based on the rate equation and traveling-wave equations with several assumptions. The model is implemented with the Symbolically Defined Devices (SDD) in Advanced Design System (ADS) and validated with experimental results. Detailed analysis regarding optical gain, harmonic and intermodulation distortions, and transmission performance is performed. The distribution of the carrier and Amplified Spontaneous Emission (ASE) is also demonstrated. Behavioral modeling of the RSOA modulator is to enable us to investigate the nonlinear distortion of the RSOA modulator from another perspective in system level. The Amplitude-to-Amplitude Conversion (AM-AM) and Amplitude-to-Phase Conversion (AM-PM) distortions of the RSOA modulator are demonstrated based on an Artificial Neural Network (ANN) and a generalized polynomial model. Another behavioral model based on Xparameters was obtained from the physical model. Compensation of the nonlinearity of the RSOA modulator is carried out based on a memory polynomial model. The nonlinear distortion of the RSOA modulator is reduced successfully. The improvement of the 3rd order intermodulation distortion is up to 17 dB. The Error Vector Magnitude (EVM) is improved from 6.1% to 2.0%. In the last part of this work, the performance of Fibre Optic Networks for Distributed and Extendible Heterogeneous Radio Architectures and Service Provisioning (FUTON) systems, which is the four-channel virtual Multiple Input Multiple Output (MIMO), is predicted by using the developed physical model. Based on Subcarrier Multiplexing (SCM) techniques, four-channel signals with 100 MHz bandwidth per channel are generated and used to drive the RSOA modulator. The transmission performance of the RSOA modulator under the broadband multi channels is depicted with the figure of merit, EVM under di erent adrature Amplitude Modulation (QAM) level of 64 and 254 for various number of Orthogonal Frequency Division Multiplexing (OFDM) subcarriers of 64, 512, 1024 and 2048.

Electromechanical properties of engineered lead free potassium sodium niobate based materials

K0.5Na0.5NbO3 (KNN), is the most promising lead free material for substituting lead zirconate titanate (PZT) which is still the market leader used for sensors and actuators. To make KNN a real competitor, it is necessary to understand and to improve its properties. This goal is pursued in the present work via different approaches aiming to study KNN intrinsic properties and then to identify appropriate strategies like doping and texturing for designing better KNN materials for an intended application. Hence, polycrystalline KNN ceramics (undoped, non-stoichiometric; NST and doped), high-quality KNN single crystals and textured KNN based ceramics were successfully synthesized and characterized in this work. Polycrystalline undoped, non-stoichiometric (NST) and Mn doped KNN ceramics were prepared by conventional ceramic processing. Structure, microstructure and electrical properties were measured. It was observed that the window for mono-phasic compositions was very narrow for both NST ceramics and Mn doped ceramics. For NST ceramics the variation of A/B ratio influenced the polarization (P-E) hysteresis loop and better piezoelectric and dielectric responses could be found for small stoichiometry deviations (A/B = 0.97). Regarding Mn doping, as compared to undoped KNN which showed leaky polarization (P-E) hysteresis loops, B-site Mn doped ceramics showed a well saturated, less-leaky hysteresis loop and a significant properties improvement. Impedance spectroscopy was used to assess the role of Mn and a relation between charge transport – defects and ferroelectric response in K0.5Na0.5NbO3 (KNN) and Mn doped KNN ceramics could be established. At room temperature the conduction in KNN which is associated with holes transport is suppressed by Mn doping. Hence Mn addition increases the resistivity of the ceramic, which proved to be very helpful for improving the saturation of the P-E loop. At high temperatures the conduction is dominated by the motion of ionized oxygen vacancies whose concentration increases with Mn doping. Single crystals of potassium sodium niobate (KNN) were grown by a modified high temperature flux method. A boron-modified flux was used to obtain the crystals at a relatively low temperature. XRD, EDS and ICP analysis proved the chemical and crystallographic quality of the crystals. The grown KNN crystals exhibit higher dielectric permittivity (29,100) at the tetragonal-to-cubic phase transition temperature, higher remnant polarization (19.4 μC/cm2) and piezoelectric coefficient (160 pC/N) when compared with the standard KNN ceramics. KNN single crystals domain structure was characterized for the first time by piezoforce response microscopy. It could be observed that <001> - oriented potassium sodium niobate (KNN) single crystals reveal a long range ordered domain pattern of parallel 180° domains with zig-zag 90° domains. From the comparison of KNN Single crystals to ceramics, It is argued that the presence in KNN single crystal (and absence in KNN ceramics) of such a long range order specific domain pattern that is its fingerprint accounts for the improved properties of single crystals. These results have broad implications for the expanded use of KNN materials, by establishing a relation between the domain patterns and the dielectric and ferroelectric response of single crystals and ceramics and by indicating ways of achieving maximised properties in KNN materials. Polarized Raman analysis of ferroelectric potassium sodium niobate (K0.5Na0.5)NbO3 (KNN) single crystals was performed. For the first time, an evidence is provided that supports the assignment of KNN single crystals structure to the monoclinic symmetry at room temperature. Intensities of A′, A″ and mixed A′+A″ phonons have been theoretically calculated and compared with the experimental data in dependence of crystal rotation, which allowed the precise determination of the Raman tensor coefficients for (non-leaking) modes in monoclinic KNN. In relation to the previous literature, this study clarifies that assigning monoclinic phase is more suitable than the orthorhombic one. In addition, this study is the basis for non-destructive assessments of domain distribution by Raman spectroscopy in KNN-based lead-free ferroelectrics with complex structures. Searching a deeper understanding of the electrical behaviour of both KNN single crystal and polycrystalline materials for the sake of designing optimized KNN materials, a comparative study at the level of charge transport and point defects was carried out by impedance spectroscopy. KNN single crystals showed lower conductivity than polycrystals from room temperature up to 200 ºC, but above this temperature polycrystalline KNN displays lower conductivity. The low temperature (T < 200 ºC) behaviour reflects the different processing conditions of both ceramics and single crystals, which account for less defects prone to charge transport in the case of single crystals. As temperature increases (T > 200 ºC) single crystals become more conductive than polycrystalline samples, in which grain boundaries act as barriers to charge transport. For even higher temperatures the conductivity difference between both is increased due to the contribution of ionic conduction in single crystals. Indeed the values of activation energy calculated to the high temperature range (T > 300 ºC) were 1.60 and 0.97 eV, confirming the charge transport due to ionic conduction and ionized oxygen vacancies in single crystals and polycrystalline KNN, respectively. It is suggested that single crystals with low defects content and improved electromechanical properties could be a better choice for room temperature applications, though at high temperatures less conductive ceramics may be the choice, depending on the targeted use. Aiming at engineering the properties of KNN polycrystals towards the performance of single crystals, the preparation and properties study of (001) – oriented (K0.5Na0.5)0.98Li0.02NbO3 (KNNL) ceramics obtained by templated grain growth (TGG) using KNN single crystals as templates was undertaken. The choice of KNN single crystals templates is related with their better properties and to their unique domain structure which were envisaged as a tool for templating better properties in KNN ceramics too. X-ray diffraction analysis revealed for the templated ceramics a monoclinic structure at room temperature and a Lotgering factor (f) of 40% which confirmed texture development. These textured ceramics exhibit a long range ordered domain pattern consisting of 90º and 180º domains, similar to the one observed in the single crystals. Enhanced dielectric (13017 at TC), ferroelectric (2Pr = 42.8 μC/cm2) and piezoelectric (d33 = 280 pC/N) properties are observed for textured KNNL ceramics as compared to the randomly oriented ones. This behaviour is suggested to be due to the long range ordered domain patterns observed in the textured ceramics. The obtained results as compared with the data previously reported on texture KNN based ceramics confirm that superior properties were found due to ordered repeated domain pattern. This study provides an useful approach towards properties improvement of KNN-based piezoelectric ceramics. Overall, the present results bring a significant contribution to the pool of knowledge on the properties of sodium potassium niobate materials: a relation between the domain patterns and di-, ferro-, and piezo-electric response of single crystals and ceramics was demonstrated and ways of engineering maximised properties in KNN materials, for example by texturing were established. This contribution is envisaged to have broad implications for the expanded use of KNN over the alternative lead-based materials.

Controle activo de estruturas adaptativas com sensores e actuadores piezoeléctricos

Nos últimos anos o estudo de estruturas inteligentes tem atraído vários investigadores devido às suas potenciais vantagens numa larga gama de aplicações, tais como controle de forma, supressão de vibrações, atenuação de ruído e detecção de dano. O uso de materiais “inteligentes” tal como os materiais piezoeléctricos na forma de lâminas ou “patches”, embebidas ou coladas na superfície de estruturas construídas de materiais compósitos, permite assim obter estruturas que por um lado são adaptativas e por outro revelam excelentes propriedades mecânicas, aumentando assim bastante o desempenho e a fiabilidade de sistemas estruturais. Os materiais piezoeléctricos têm a propriedade de gerar uma carga eléctrica sob a acção duma carga mecânica e o reverso, isto é, aplicando um campo eléctrico nos elementos piezoeléctricos da estrutura, esta deforma-se. Neste trabalho, é apresentado um modelo de elementos finitos, baseado na teoria clássica de placas, desenvolvido para a análise do controle activo em estática e dinâmica lineares de estruturas integrando sensores e actuadores piezoeléctrico na forma de lâminas, os quais introduzem um grau de liberdade referente ao potencial eléctrico, por cada camada piezoeléctrica do elemento finito. É utilizado método de Newmark para a solução iterativa das equações de equilíbrio. Apresentam-se os resultados obtidos em três exemplos ilustrativos.

Análise do risco associado à rotura da barragem de Fernandilho. Modelação hidrodinâmica da onda de cheia

As cheias provocadas por roturas de barragens constituem eventos catastróficos que podem conduzir a perdas consideráveis de vidas humanas e de bens materiais. Para melhorar a segurança, em relação a este risco potencial, torna-se necessário avaliar a segurança da população e das estruturas, existentes ou a construir, e estimar a extensão dos danos em caso de acidente. Neste estudo, após a indicação da legislação em vigor, apresenta-se uma descrição sumária do modelo DamBreak. Descrevem-se as equações básicas que definem o modelo de formação da brecha e que caracterizam o processo de formação da cheia. Também se referem as equações de Saint-Venant, as quais caracterizam o processo de propagação da cheia. De seguida é feita uma descrição da barragem de Fernandilho, do leito do Barranco d

Análise geometricamente não-linear de estruturas adaptativas piezolaminadas

Neste trabalho apresenta-se um modelo de elementos finitos, baseado na teoria clássica de placas, para a análise linear e não-linear de estruturas do tipo placa/casca integrando sensores e actuadores piezoeléctricos. É usado um simples e eficiente elemento placa/casca triangular plano de 3 nós, e em cuja formulação se introduz um grau de liberdade referente ao potencial eléctrico, por cada camada piezoeléctrica do elemento finito. É utilizada a formulação Lagrangeana actualizada associada à tecnica de Newton - Raphson para a solução iterativa das equações de equilibrio .O modelo pode ser aplicado a cascas piezolaminadas com geometria e carregamento arbitrários. Apresentam-se vários exemplos ilustrativos cujos resultados mostram a eficiencia do modelo proposto.

Placas elasto-electro-magnéticas

Neste trabalho apresenta-se um modelo de elementos finitos, baseado na teoria clássica de placas, para a análise linear e não-linear de estruturas do tipo placa/casca integrando sensores e actuadores piezoeléctricos. É usado um simples e eficiente elemento placa/casca triangular plano de 3 nós, e em cuja formulação se introduz um grau de liberdade referente ao potencial eléctrico, por cada camada piezoeléctrica do elemento finito. É utilizada a formulação Lagrangeana actualizada associada à tecnica de Newton - Raphson para a solução iterativa das equações de equilibrio .O modelo pode ser aplicado a cascas piezolaminadas com geometria e carregamento arbitrários. Apresentam-se vários exemplos ilustrativos cujos resultados mostram a eficiencia do modelo proposto.

Métodos dos momentos

Os métodos utilizados na resolução de problemas, nos vários ramos da Engenharia ou ciências aplicadas, baseiam-se, actualmente, em uma de duas categorias: métodos analíticos e métodos numéricos. É preferível a utilização dos métodos analíticos, na resolução de equações ou outros modelos matemáticos, sempre que possível, uma vez que formamões gerais em vez de particularizadas, para além de uma maior informação.

Modelação do impacte produzido pelo aumento de nutrientes na produção primária do estuário do Sado

Os modelos matemáticos são ferramentas poderosas para avaliar o impacto humano em sistemas costeiros e estuarinos de estrutura complexa. Neste artigo utiliza-se um sistema de modelos para caracterizar o estado trófico do estuário do Sado e para prever a sua evolução baseada em diferentes cenários de intervenção humana. Analisa-se em pormenor o efeito do aumento de nutrientes introduzidos no estuário através do rio. Este estudo permite tirar conclusões sobre o impacte que esse aumento terá sobre o estuário. O sistema de modelação é composto por um modelo hidrodinâmico tridimensional às equações primitivas e modelos de transporte acoplados a modelos de transporte de sedimentos, modelos de qualidade da água e modelos ecológicos. O modelo hidrodinâmico encontra-se calibrado e validado para o sistema. Utilizam-se dados de campo de nutrientes e de fitoplâncton para comparação com os resultados do modelo. Analisam-se as características do sistema e efectua-se uma análise de sensibilidade do comportamento biológico do estuário face ao aumento da carga de nutrientes. As conclusões principais são: Confirmação da adequabilidade da metodologia utilizada para a simulação da situação actual do estuário; Pequena sensibilidade do estuário a aumentos moderados da carga de nutrientes; Sensibilidade moderada do estuário para elevadas cargas de nutrientes.

Qualidade da água na Ria Formosa: combinando técnicas de amostragem e de modelação para uma gestão integrada das águas residuais

Aplica-se uma técnica combinada de amostragem clássica, perfilagem in situ com sonda multiparamétrica e modelação matemática para avaliar a qualidade da água da Ria Formosa como meio receptor dos efluentes das ETAR’s. O sistema de modelação desenvolvido baseia-se num modelo hidrodinâmico que resolve as equações de águas pouco profundas pelo método dos volumes finitos acoplado a um modelo de transporte Lagrangiano para a simulação das plumas microbiológicas. Este sistema é usado como ferramenta de gestão para ensaiar diferentes cenários de descarga com variação da localização do ponto de descarga, analisando o impacto das alternativas no meio receptor.

Compreensão e interpretação em matemática: dificuldades de alunos do 9º ano na resolução de problemas

Dissertação mest., Matemática, Universidade do Algarve, 2009