Accounting for adsorption and desorption in Lattice Boltzmann simulations

We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. With the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g., in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic and also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a lattice Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account in the time-dependent diffusion coefficient in a more complex porous medium.

Estimating the value of the metal-ligand bond dissociation enthalpy<D> (M-L) for adducts using empirical equations supported by TG data

In this work is presented and tested (for 106 adducts, mainly of the zinc group halides) two empirical equations supported in TG data to estimate the value of the metal-ligand bond dissociation enthalpy for adducts: <D> (M-O) = t i / g if t i < 420 K and <D> (M-O) = (t i / g ) - 7,75 . 10-2 . t i if t i > 420 K. In this empirical equations, t i is the thermodynamic temperature of the beginning of the thermal decomposition of the adduct, as determined by thermogravimetry, andg is a constant factor that is function of the metal halide considered and of the number of ligands, but is not dependant of the ligand itself. To half of the tested adducts the difference between experimental and calculated values was less than 5%. To about 80% of the tested adducts, the difference between the experimental (calorimetric) and the calculated (using the proposed equations) values are less than 15%.

Modeling reaction kinetics and mass transfer in ozonation in water solutions

This dissertation is based on four articles dealing with modeling of ozonation. The literature part of this considers some models for hydrodynamics in bubble column simulation. A literature review of methods for obtaining mass transfer coefficients is presented. The methods presented to obtain mass transfer are general models and can be applied to any gas-liquid system. Ozonation reaction models and methods for obtaining stoichiometric coefficients and reaction rate coefficients for ozonation reactions are discussed in the final section of the literature part. In the first article, ozone gas-liquid mass transfer into water in a bubble column was investigated for different pH values. A more general method for estimation of mass transfer and Henry’s coefficient was developed from the Beltrán method. The ozone volumetric mass transfer coefficient and the Henry’s coefficient were determined simultaneously by parameter estimation using a nonlinear optimization method. A minor dependence of the Henry’s law constant on pH was detected at the pH range 4 - 9. In the second article, a new method using the axial dispersion model for estimation of ozone self-decomposition kinetics in a semi-batch bubble column reactor was developed. The reaction rate coefficients for literature equations of ozone decomposition and the gas phase dispersion coefficient were estimated and compared with the literature data. The reaction order in the pH range 7-10 with respect to ozone 1.12 and 0.51 the hydroxyl ion were obtained, which is in good agreement with literature. The model parameters were determined by parameter estimation using a nonlinear optimization method. Sensitivity analysis was conducted using object function method to obtain information about the reliability and identifiability of the estimated parameters. In the third article, the reaction rate coefficients and the stoichiometric coefficients in the reaction of ozone with the model component p-nitrophenol were estimated at low pH of water using nonlinear optimization. A novel method for estimation of multireaction model parameters in ozonation was developed. In this method the concentration of unknown intermediate compounds is presented as a residual COD (chemical oxygen demand) calculated from the measured COD and the theoretical COD for the known species. The decomposition rate of p-nitrophenol on the pathway producing hydroquinone was found to be about two times faster than the p-nitrophenol decomposition rate on the pathway producing 4- nitrocatechol. In the fourth article, the reaction kinetics of p-nitrophenol ozonation was studied in a bubble column at pH 2. Using the new reaction kinetic model presented in the previous article, the reaction kinetic parameters, rate coefficients, and stoichiometric coefficients as well as the mass transfer coefficient were estimated with nonlinear estimation. The decomposition rate of pnitrophenol was found to be equal both on the pathway producing hydroquinone and on the path way producing 4-nitrocathecol. Comparison of the rate coefficients with the case at initial pH 5 indicates that the p-nitrophenol degradation producing 4- nitrocathecol is more selective towards molecular ozone than the reaction producing hydroquinone. The identifiability and reliability of the estimated parameters were analyzed with the Marcov chain Monte Carlo (MCMC) method. @All rights reserved. No part of the publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior permission of the author.

Desenvolvimento e aplicação de células eletroquímicas em camada delgada

Theoretical and practical aspects concerning construction and use of a thin layer electrochemical cell are presented. Construction was realized by a simple technic and geometry was optimized in order to minimize the problems of electrical resistance. A well known redox system was studied ((Fe(III)(CN)6(3-)/Fe(II)(CN) 6(4-)) using two experimental methods, cyclic voltammetry and pulse chronopotentiometry. A numerical integration based-program was developed to calculate the voltammetric current in case of nernstian and non-nernstian behaviours and a diffusional model was used to treat the chronopotentiometric data. Thermodynamic (potential, concentration) and kinetic parameters (diffusion coefficient) were successfully determined for the redox system studied in this work.

Offdesign operation of kraft recovery boiler

A model to solve heat and mass balances during the offdesign load calculations was created. These equations are complex and nonlinear. The main new ideas used in the created offdesign model of a kraft recovery boiler are the use of heat flows as torn iteration variables instead of the current practice of using the mass flows, vectorizing equation solving, thus speeding up the process, using non dimensional variables for solving the multiple heat transfer surface problem and using a new procedure for calculating pressure losses. Recovery boiler heat and mass balances are reduced to vector form. It is shown that these vectorized equations can be solved virtually without iteration. The iteration speed is enhanced by the use of the derived method of calculating multiple heat transfer surfaces simultaneously. To achieve this quick convergence the heat flows were used as the torn iteration parameters. A new method to handle pressure loss calculations with linearization was presented. This method enabled less time to be spent calculating pressure losses. The derived vector representation of the steam generator was used to calculate offdesign operation parameters for a 3000 tds/d example recovery boiler. The model was used to study recovery boiler part load operation and the effect of the black liquor dry solids increase on recovery boiler dimensioning. Heat flows to surface elements for part load calculations can be closely approximated with a previously defined exponent function. The exponential method can be used for the prediction of fouling in kraft recovery boilers. For similar furnaces the firing of 80 % dry solids liquor produces lower hearth heat release rate than the 65 % dry solids liquor if we fire at constant steam flow. The furnace outlet temperatures show that capacity increase with firing rate increase produces higher loadings than capacity increase with dry solids increase. The economizers, boiler banks and furnaces can be dimensioned smaller if we increase the black liquor dry solids content. The main problem with increased black liquor dry solids content is the decrease in the heat available to superheat. Whenever possible the furnace exit temperature should be increased by decreasing the furnace height. The increase in the furnace exit temperature is usually opposed because of fear of increased corrosion.

Evanescent field of vectorial highly non-paraxial beams

In terms of the Fourier spectrum, a simple but general analytical expression is given for the evanescent field associated to a certain kind of non-paraxial exact solutions of the Maxwell equations. This expression enables one to compare the relative weight of the evanescent wave with regard to the propagating field. In addition, in those cases in which the evanescent term is significant, the magnitude of the field components across the transverse profile (including the evanescent features) can be determined. These results are applied to some illustrative examples.

The influence of a certain volume of fluid in the lower tear meniscus on intraocular pressure measurements using the non-contact tonometer

Background: Intraocular pressure (IOP) is the pressure inside the eye that helps to maintain the integrity and the suitable form of the ocular globe. Precise and accurate measures of IOP are needed for the diagnosis as well as follow-up of glaucoma. In daily clinical practice, Goldmann applanation tonometer (GAT) and Non-contact tonometer (NCT) are the most common devices for measuring IOP. A close agreement between these methods has been showed, particularly in normotensive patients and a poor agreement, especially when IOP levels are above the normal range. Ophthalmologists have noticed a poor agreement between NCT and GAT, observing that by using NCT and after comparing with GAT, there is an overestimation of IOP readings, and particularly it occurs when the eyes are tearful. Previous studies investigate the effect of tears in Non-contact tonometer readings by the instillation of artificial tears, concluding in one of the studies that the variation was less than 1mmHg and not clinically significant, in contrast with another study which the increases were sadistically significant. Tear menisci are a thin strip of tear fluid located between the bulbar conjunctiva and the eyelid margins. We think that the overestimation of IOP readings using NCT could be due to the presence of a higher volume of tear in the lower tear meniscus which might cause an optical interference in the optoelectronic applanation monitoring system of this deviceObjectives: To research the influence of a certain volume of fluid in the lower tear meniscus on IOP measurements using the NCT in healthy eyes. Moreover, to investigate the agreement between IOP readings obtained by NCT and GAT in the presence and absence of this volume of fluidMethods: The study design will be transversal for diagnostic tests of repeated measures. We will study patients with no ocular pathology and IOP<21mmHg. It will consist in the measurement of IOP using NCT before and after the instillation of COLIRCUSÍ FLUOTEST, used as a volume of fluid in the lower tear meniscus, to observe if there will be differences using the paired t-test. Moreover, we will take IOP measures by GAT in order to know the agreement between these methods after and before the application of these eyedrops, using the ICC (intraclass correlation coefficient) and the Bland-Altmann method

On the Application of Beam on Elastic Foundation Theory to the Analysis of Stiffened Plate Strips

The main objective of this thesis is to show that plate strips subjected to transverse line loads can be analysed by using the beam on elastic foundation (BEF) approach. It is shown that the elastic behaviour of both the centre line section of a semi infinite plate supported along two edges, and the free edge of a cantilever plate strip can be accurately predicted by calculations based on the two parameter BEF theory. The transverse bending stiffness of the plate strip forms the foundation. The foundation modulus is shown, mathematically and physically, to be the zero order term of the fourth order differential equation governing the behaviour of BEF, whereas the torsion rigidity of the plate acts like pre tension in the second order term. Direct equivalence is obtained for harmonic line loading by comparing the differential equations of Levy's method (a simply supported plate) with the BEF method. By equating the second and zero order terms of the semi infinite BEF model for each harmonic component, two parameters are obtained for a simply supported plate of width B: the characteristic length, 1/ λ, and the normalized sum, n, being the effect of axial loading and stiffening resulting from the torsion stiffness, nlin. This procedure gives the following result for the first mode when a uniaxial stress field was assumed (ν = 0): 1/λ = √2B/π and nlin = 1. For constant line loading, which is the superimposition of harmonic components, slightly differing foundation parameters are obtained when the maximum deflection and bending moment values of the theoretical plate, with v = 0, and BEF analysis solutions are equated: 1 /λ= 1.47B/π and nlin. = 0.59 for a simply supported plate; and 1/λ = 0.99B/π and nlin = 0.25 for a fixed plate. The BEF parameters of the plate strip with a free edge are determined based solely on finite element analysis (FEA) results: 1/λ = 1.29B/π and nlin. = 0.65, where B is the double width of the cantilever plate strip. The stress biaxial, v > 0, is shown not to affect the values of the BEF parameters significantly the result of the geometric nonlinearity caused by in plane, axial and biaxial loading is studied theoretically by comparing the differential equations of Levy's method with the BEF approach. The BEF model is generalised to take into account the elastic rotation stiffness of the longitudinal edges. Finally, formulae are presented that take into account the effect of Poisson's ratio, and geometric non linearity, on bending behaviour resulting from axial and transverse inplane loading. It is also shown that the BEF parameters of the semi infinite model are valid for linear elastic analysis of a plate strip of finite length. The BEF model was verified by applying it to the analysis of bending stresses caused by misalignments in a laboratory test panel. In summary, it can be concluded that the advantages of the BEF theory are that it is a simple tool, and that it is accurate enough for specific stress analysis of semi infinite and finite plate bending problems.

Correlação entre parâmetros da cinética de intumescimento com características estruturais e hidrofílicas de hidrogéis de poliacrilamida e metilcelulose

In this paper, the effects of acrylamide (AAm), methylcellulose (MC) contents, pH and ionic strength on kinetic, network and hydrophilic properties of polyacrylamide and methylcellulose hydrogels were investigated. The hydrogels were characterized by evaluating of network [average molecular weight between crosslinks (M C), crosslink density (q) and the number of elastically effective chains (Ve)], and kinetic parameters [diffusional exponent (n), diffusion constant (k) and diffusion coefficient (D)]. Such properties were controlled by adjusting of the AAm, MC contents, pH and ionic strength factors. Due to high hydrophilicity and fast water-uptake, the PAAm-MC hydrogels can be considered as materials for potential applications in agricultural fields, mainly in controlled release of water or pesticides.

Dehydration and volatilization non isothermic kinetic of the solid state aluminium 8-hydroxyquinolinate

Al(C9H6ON)3.2.5H2O was precipitated from the mixture of an aqueous solution of aluminium ion and an acid solution of 8-hydroxyquinoline, by increasing the pH value to 9.5 with ammonia aqueous solution. The TG curves in nitrogen atmosphere present mass losses due to dehydration, partial volatilisation (sublimation plus vaporisation) of the anhydrous compound followed by thermal decomposition with the formation of a mixture of carbonaceous and residues. The relation between sublimation and vaporisation depends on the heating rate used. The non isothermic integral isoconventional methods as linear equations of Ozawa-Flynn-Wall and Kissinger-Akahira-Sunose (KAS) were used to obtain the kinetic parameters from TG and DTA curves, respectively. Despite the fact that both dehydration and volatilisation reactions follow the linearity by using both methods, only for the volatilisation reaction the validity condition, 20<= E/RT<= 50, was verified.

Learning Robot Motions from Visual Sensing

Learning from demonstration becomes increasingly popular as an efficient way of robot programming. Not only a scientific interest acts as an inspiration in this case but also the possibility of producing the machines that would find application in different areas of life: robots helping with daily routine at home, high performance automata in industries or friendly toys for children. One way to teach a robot to fulfill complex tasks is to start with simple training exercises, combining them to form more difficult behavior. The objective of the Master’s thesis work was to study robot programming with visual input. Dynamic movement primitives (DMPs) were chosen as a tool for motion learning and generation. Assuming a movement to be a spring system influenced by an external force, making this system move, DMPs represent the motion as a set of non-linear differential equations. During the experiments the properties of DMP, such as temporal and spacial invariance, were examined. The effect of the DMP parameters, including spring coefficient, damping factor, temporal scaling, on the trajectory generated were studied.

A generalized 1D particle transport method for convection diffusion reaction model

This work is devoted to the development of numerical method to deal with convection diffusion dominated problem with reaction term, non - stiff chemical reaction and stiff chemical reaction. The technique is based on the unifying Eulerian - Lagrangian schemes (particle transport method) under the framework of operator splitting method. In the computational domain, the particle set is assigned to solve the convection reaction subproblem along the characteristic curves created by convective velocity. At each time step, convection, diffusion and reaction terms are solved separately by assuming that, each phenomenon occurs separately in a sequential fashion. Moreover, adaptivities and projection techniques are used to add particles in the regions of high gradients (steep fronts) and discontinuities and transfer a solution from particle set onto grid point respectively. The numerical results show that, the particle transport method has improved the solutions of CDR problems. Nevertheless, the method is time consumer when compared with other classical technique e.g., method of lines. Apart from this advantage, the particle transport method can be used to simulate problems that involve movingsteep/smooth fronts such as separation of two or more elements in the system.

Electrochemical study and dependence of 'transition state' in Co(II) and Ni(II) complexes with some antibiotics and cephalothin

Electrode kinetics and study of 'transition state' with applied potential in case of [M - antibiotics - cephalothin] system were reported at pH = 7.30 ± 0.01 at suitable supporting electrolyte at 25.0ºC. The M = Co or Ni and antibiotics were doxycycline, chlortetracycline, oxytetracycline, tetracycline, minocycline, amoxicillin and chloramphenicol used as primary ligands and cephalothin as secondary ligand. Kinetic parameters viz. transfer coefficient (a), degree of irreversibility (l), diffusion coefficient (D) and rate constant (k) were determined. The values of a and k varied from 0.41 to 0.59 and 2.60 X 10-3 cm s-1 to 9.67 X 10-3 cm s-1 in case of [Co - antibiotics - cephalothin] system. In case of [Ni - antibiotics - cephalothin], a and k varied from 0.41 to 0.58 and 2.34 X 10-3 cm s-1 to 9.19 X 10-3 cm s-1 respectively confirmed that transition state behaves between oxidant and reductant response to applied potential and it adjusts it self in such a way that the same is located midway between dropping mercury electrode and solution interface. The values of rate constant confirmed the quasireversible nature of electrode processes. The stability constants (logb) of complexes were also determined.

Non-Markovianity and Quantum Correlations in Qubit-Systems

In this Thesis I discuss the exact dynamics of simple non-Markovian systems. I focus on fundamental questions at the core of non-Markovian theory and investigate the dynamics of quantum correlations under non-Markovian decoherence. In the first context I present the connection between two different non-Markovian approaches, and compare two distinct definitions of non-Markovianity. The general aim is to characterize in exemplary cases which part of the environment is responsible for the feedback of information typical of non- Markovian dynamics. I also show how such a feedback of information is not always described by certain types of master equations commonly used to tackle non-Markovian dynamics. In the second context I characterize the dynamics of two qubits in a common non-Markovian reservoir, and introduce a new dynamical effect in a wellknown model, i.e., two qubits under depolarizing channels. In the first model the exact solution of the dynamics is found, and the entanglement behavior is extensively studied. The non-Markovianity of the reservoir and reservoirmediated-interaction between the qubits cause non-trivial dynamical features. The dynamical interplay between different types of correlations is also investigated. In the second model the study of quantum and classical correlations demonstrates the existence of a new effect: the sudden transition between classical and quantum decoherence. This phenomenon involves the complete preservation of the initial quantum correlations for long intervals of time of the order of the relaxation time of the system.

Modeling of electrolyte sorption – from phase equilibria to dynamic separation systems

In this thesis, general approach is devised to model electrolyte sorption from aqueous solutions on solid materials. Electrolyte sorption is often considered as unwanted phenomenon in ion exchange and its potential as an independent separation method has not been fully explored. The solid sorbents studied here are porous and non-porous organic or inorganic materials with or without specific functional groups attached on the solid matrix. Accordingly, the sorption mechanisms include physical adsorption, chemisorption on the functional groups and partition restricted by electrostatic or steric factors. The model is tested in four Cases Studies dealing with chelating adsorption of transition metal mixtures, physical adsorption of metal and metalloid complexes from chloride solutions, size exclusion of electrolytes in nano-porous materials and electrolyte exclusion of electrolyte/non-electrolyte mixtures. The model parameters are estimated using experimental data from equilibrium and batch kinetic measurements, and they are used to simulate actual single-column fixed-bed separations. Phase equilibrium between the solution and solid phases is described using thermodynamic Gibbs-Donnan model and various adsorption models depending on the properties of the sorbent. The 3-dimensional thermodynamic approach is used for volume sorption in gel-type ion exchangers and in nano-porous adsorbents, and satisfactory correlation is obtained provided that both mixing and exclusion effects are adequately taken into account. 2-Dimensional surface adsorption models are successfully applied to physical adsorption of complex species and to chelating adsorption of transition metal salts. In the latter case, comparison is also made with complex formation models. Results of the mass transport studies show that uptake rates even in a competitive high-affinity system can be described by constant diffusion coefficients, when the adsorbent structure and the phase equilibrium conditions are adequately included in the model. Furthermore, a simplified solution based on the linear driving force approximation and the shrinking-core model is developed for very non-linear adsorption systems. In each Case Study, the actual separation is carried out batch-wise in fixed-beds and the experimental data are simulated/correlated using the parameters derived from equilibrium and kinetic data. Good agreement between the calculated and experimental break-through curves is usually obtained indicating that the proposed approach is useful in systems, which at first sight are very different. For example, the important improvement in copper separation from concentrated zinc sulfate solution at elevated temperatures can be correctly predicted by the model. In some cases, however, re-adjustment of model parameters is needed due to e.g. high solution viscosity.