Modelling and optimization of the ion exchange membrane bioreactor for removal of anionic pollutants from drinking water streams

Dissertação para obtenção do Grau de Doutor em Engenharia Química, especialidade de Engenharia Bioquímica

Experimental and model-based analysis of the steady-state and dynamic operating behaviour of the direct methanol fuel cell (DMFC)

Direct methanol fuel cell, DMFC, model, mass transport, Maxwell-Stefan, Flory-Huggins, crossover, polymer electrolyte membrane, Nafion

Contractive probability metrics and asymptotic behavior of dissipative kinetic equations

The present notes are intended to present a detailed review of the existing results in dissipative kinetic theory which make use of the contraction properties of two main families of probability metrics: optimal mass transport and Fourier-based metrics. The first part of the notes is devoted to a self-consistent summary and presentation of the properties of both probability metrics, including new aspects on the relationships between them and other metrics of wide use in probability theory. These results are of independent interest with potential use in other contexts in Partial Differential Equations and Probability Theory. The second part of the notes makes a different presentation of the asymptotic behavior of Inelastic Maxwell Models than the one presented in the literature and it shows a new example of application: particle's bath heating. We show how starting from the contraction properties in probability metrics, one can deduce the existence, uniqueness and asymptotic stability in classical spaces. A global strategy with this aim is set up and applied in two dissipative models.

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

State-of-the-art of 3D cultures (organs-on-a-chip) in safety testing and pathophysiology.

Integrated approaches using different in vitro methods in combination with bioinformatics can (i) increase the success rate and speed of drug development; (ii) improve the accuracy of toxicological risk assessment; and (iii) increase our understanding of disease. Three-dimensional (3D) cell culture models are important building blocks of this strategy which has emerged during the last years. The majority of these models are organotypic, i.e., they aim to reproduce major functions of an organ or organ system. This implies in many cases that more than one cell type forms the 3D structure, and often matrix elements play an important role. This review summarizes the state of the art concerning commonalities of the different models. For instance, the theory of mass transport/metabolite exchange in 3D systems and the special analytical requirements for test endpoints in organotypic cultures are discussed in detail. In the next part, 3D model systems for selected organs--liver, lung, skin, brain--are presented and characterized in dedicated chapters. Also, 3D approaches to the modeling of tumors are presented and discussed. All chapters give a historical background, illustrate the large variety of approaches, and highlight up- and downsides as well as specific requirements. Moreover, they refer to the application in disease modeling, drug discovery and safety assessment. Finally, consensus recommendations indicate a roadmap for the successful implementation of 3D models in routine screening. It is expected that the use of such models will accelerate progress by reducing error rates and wrong predictions from compound testing.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Impaired uptake of glutathione by hepatic mitochondria from chronic ethanol-fed rats. Tracer kinetic studies in vitro and in vivo and susceptibility to oxidant stress.

Isolated hepatocytes incubated with [35S]-methionine were examined for the time-dependent accumulation of [35S]-glutathione (GSH) in cytosol and mitochondria, the latter confirmed by density gradient purification. In GSH-depleted and -repleted hepatocytes, the increase of specific activity of mitochondrial GSH lagged behind cytosol, reaching nearly the same specific activity by 1-2 h. However, in hepatocytes from ethanol-fed rats, the rate of increase of total GSH specific radioactivity in mitochondria was markedly suppressed. In in vivo steady-state experiments, the mass transport of GSH from cytosol to mitochondria and vice versa was 18 nmol/min per g liver, indicating that the half-life of mitochondrial GSH was approximately 18 min in controls. The fractional transport rate of GSH from cytosol to mitochondria, but not mitochondria to cytosol, was significantly reduced in the livers of ethanol-fed rats. Thus, ethanol-fed rats exhibit a decreased mitochondrial GSH pool size due to an impaired entry of cytosol GSH into mitochondria. Hepatocytes from ethanol-fed rats exhibited a greater susceptibility to the oxidant stress-induced cell death from tert-butylhydroperoxide. Incubation with glutathione monoethyl ester normalized the mitochondrial GSH and protected against the increased susceptibility to t-butylhydroperoxide, which was directly related to the lowered mitochondrial GSH pool size in ethanol-fed cells.

Ultramicroeletrodos. Parte II: construção e aplicações

This work describes the techniques of construction and several applications of ultramicroelectrodes in electrochemistry and electroanalytical chemistry. Disc shaped UME are produced by embedding metal wires on insulating materials such as glass or epoxy resin. In the field of electrochemistry, UME have been applied in studies of the hydrogen evolution reaction and the electrocrystallization of metals. The negligible values of sensibility for ohmic drop and the enhanced mass transport rate by spherical diffusion are the main advantages of UME in these applications. New important conclusions regarding the phenomena under study were drawn from the experimental results. The applications in electroanalytical chemistry involved the determination of contaminants such as heavy metals and nitrites in natural waters and food products. The use of UME requires little sample manipulation and, in general, no need for oxygen removal or the addition of supporting electrolytes.

Transporte de massa em polímeros intrinsecamente condutores: importância, técnicas e modelos teóricos

In this work we discuss the aspects related to the phenomenon of mass transport in thin electroactive polymer films. Such phenomenon must be considered because the properties and consequent applications of these materials largely depend on the movement of charge carriers, i.e. ions, electrons or holes. The most recent majority of the techniques, methods and theoretical models used in this type of study are gathered and discussed, providing an easy and critical way for choosing the methodology for an investigation.

Princípio e utilização da técnica de efeito miragem de concentração

The aim of this work is to present the theoretical and experimental aspects of the mirage effect technique. We are especially interested in the concentration mirage effect, which is a powerful tool in the study of electrochemical reactions that produce ionic movements close to electrodes and to get some fundamental information on mass transport and charge transfer during electrochemical processes. Limitations of this technique are discussed as well as the recent attempts to overcome them.

Aplicação do ultra-som em sistemas eletroquímicos: considerações teóricas e experimentais

The aim of this review is to present and discuss the applications of ultrasound in electrochemical systems such as in sonoelectroanalysis and sonoelectrolysis for the electrochemical combustion of organic compounds. Initially, theoretical and experimental aspects are discussed, particularly those related to the enhancement of mass transport and the surface cleaning effects. Some results are included to illustrate alternative geometries for the experimental measurements and the working electrodes used in these systems. In the sequence, the available publications are presented and discussed to demonstrate that ultrasound combined with electrochemical techniques is a powerful set-up for the detection of analytes such as metals and/or organic compounds in hostile media and for the effective destruction of toxic organic substances. At the end, a table summarizes the results already published in the literature.

Chelating Adsorbents in Purification of Hydrometallurgical Solutions

In this thesis, equilibrium and dynamic sorption properties of weakly basic chelating adsorbents were studied to explain removal of copper, nickel from a concentrated zinc sulfate solution in a hydrometallurgical process. Silica-supported chelating composites containing either branched poly(ethyleneimine) (BPEI) or 2-(aminomethyl)pyridine (AMP) as a functional group were used. The adsorbents are commercially available from Purity Systems Inc, USA as WP-1® and CuWRAM®, respectively. The fundamental interactions between the adsorbents, sulfuric acid and metal sulfates were studied in detail and the results were used to find the best conditions for removal of copper and nickel from an authentic ZnSO4 process solution. In particular, the effect of acid concentration and temperature on the separation efficiency was considered. Both experimental and modeling aspectswere covered in all cases. Metal sorption is considerably affected by the chemical properties of the studied adsorbents and by the separation conditions. In the case of WP-1, acid affinity is so high that column separation of copper, nickel and zinc has to be done using the adsorbent in base-form. On the other hand, the basicity of CuWRAM is significantly lower and protonated adsorbent can be used. Increasing temperature decreases the basicity and the metals affinity of both adsorbents, but the uptake capacities remain practically unchanged. Moreover, increasing temperature substantially enhances intra-particle mass transport and decreases viscosities thus allowing significantly higher feed flow rates in the fixed-bed separation. The copper selectivity of both adsorbents is very high even in the presence of a 250-fold excess of zinc. However, because of the basicity of WP-1, metal precipitation is a serious problem and therefore only CuWRAM is suitable for the practical industrial application. The optimum temperature for copper removal appears to be around 60 oC and an alternative solution purification method is proposed. The Ni/Zn selectivity of both WP-1 and CuWRAM is insufficient for removal of the very small amounts of nickel present in the concentrated ZnSO4 solution.

Eletrodo compósito à base de grafite-araldite®: aplicações didáticas

A composite electrode prepared by mixing a commercial epoxy resin Araldite® and graphite powder is proposed to be used in didactic experiments. The electrode is prepared by the students and applied in simple experiments to demonstrate the effect of the composite composition on the conductivity and the voltammetric response of the resulting electrode, as well as the response in relation to the scan rate dependence on mass transport. The possibility of using the composite electrode in quantitative analysis is also demonstrated.

Potencialidades do uso de ultrassom na determinação do pesticida carbaril empregando eletrodos de diamante

The potentiality of the use of ultrasound radiation in association with a boron-doped diamond electrode was evaluated on the voltammetric determination of the pesticide carbaryl. Improvements in the sensitivity, limit of detection and reproducibility of the measurements were observed due to both, the enhancement of mass transport and the cleaning of the electrode surface provided by ultrasound. Satisfactory recovery levels for carbaryl in pure water (96-98%) and pineapple juice (89-92%) for quiescent and sonovoltammetric methodologies were obtained. These methodologies can be alternative tools for the analyses of pesticides in fruit samples, mainly the insonated condition that improve the analytical performance and dispense intermediary cleanings of the electrode surface.