Um problema de dispersão de poluentes em rios e canais por meio do método de separação de variáveis

Neste trabalho obtém-se uma solução analítica para a equação de advecção-difusão aplicada a problemas de dispersão de poluentes em rios e canais. Para tanto, consideram-se os casos unidimensionais e bidimensionais em regime transiente com coeficientes de difusividade e velocidades constantes. A abordagem utilizada para a resolução deste problema é o método de Separação de Variáveis. Os modelos resolvidos foram simulados utilizando o MatLab. Apresentam-se os resultados das simulações numéricas em formato gráfico. Os resultados de algumas simulações numéricas existem na literatura e puderam ser comparados. O modelo proposto mostrou-se coerente em relação aos dados considerados. Para outras simulações não foram encontrados comparativos na literatura, todavia esses problemas governados por equações diferenciais parciais, mesmo lineares, não são de fácil solução analítica. Sendo que, muitas delas representam importantes problemas de matemática e física, com diversas aplicações na engenharia. Dessa forma, é de grande importância a disponibilidade de um maior número de problemas-teste para avaliação de desempenho de formulações numéricas, cada vez mais eficazes, já que soluções analíticas oferecem uma base mais segura para comparação de resultados.

Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems on Anisotropically Refined Meshes

In this paper we consider the a posteriori and a priori error analysis of discontinuous Galerkin interior penalty methods for second-order partial differential equations with nonnegative characteristic form on anisotropically refined computational meshes. In particular, we discuss the question of error estimation for linear target functionals, such as the outflow flux and the local average of the solution. Based on our a posteriori error bound we design and implement the corresponding adaptive algorithm to ensure reliable and efficient control of the error in the prescribed functional to within a given tolerance. This involves exploiting both local isotropic and anisotropic mesh refinement. The theoretical results are illustrated by a series of numerical experiments.

Discontinuous Galerkin Methods on hp-Anisotropic Meshes I: A Priori Error Analysis

We consider the a priori error analysis of hp-version interior penalty discontinuous Galerkin methods for second-order partial differential equations with nonnegative characteristic form under weak assumptions on the mesh design and the local finite element spaces employed. In particular, we prove a priori hp-error bounds for linear target functionals of the solution, on (possibly) anisotropic computational meshes with anisotropic tensor-product polynomial basis functions. The theoretical results are illustrated by a numerical experiment.

A bistable reaction-diffusion system in a stretching flow

We examine the evolution of a bistable reaction in a one-dimensional stretching flow, as a model for chaotic advection. We derive two reduced systems of ordinary differential equations (ODEs) for the dynamics of the governing advection-reaction-diffusion partial differential equations (PDE), for pulse-like and for plateau-like solutions, based on a non-perturbative approach. This reduction allows us to study the dynamics in two cases: first, close to a saddle-node bifurcation at which a pair of nontrivial steady states are born as the dimensionless reaction rate (Damkoehler number) is increased, and, second, for large Damkoehler number, far away from the bifurcation. The main aim is to investigate the initial-value problem and to determine when an initial condition subject to chaotic stirring will decay to zero and when it will give rise to a nonzero final state. Comparisons with full PDE simulations show that the reduced pulse model accurately predicts the threshold amplitude for a pulse initial condition to give rise to a nontrivial final steady state, and that the reduced plateau model gives an accurate picture of the dynamics of the system at large Damkoehler number. Published in Physica D (2006)

Discontinuous Galerkin Methods on hp-Anisotropic Meshes II: A Posteriori Error Analysis and Adaptivity

We consider the a posteriori error analysis and hp-adaptation strategies for hp-version interior penalty discontinuous Galerkin methods for second-order partial differential equations with nonnegative characteristic form on anisotropically refined computational meshes with anisotropically enriched elemental polynomial degrees. In particular, we exploit duality based hp-error estimates for linear target functionals of the solution and design and implement the corresponding adaptive algorithms to ensure reliable and efficient control of the error in the prescribed functional to within a given tolerance. This involves exploiting both local isotropic and anisotropic mesh refinement and isotropic and anisotropic polynomial degree enrichment. The superiority of the proposed algorithm in comparison with standard hp-isotropic mesh refinement algorithms and an h-anisotropic/p-isotropic adaptive procedure is illustrated by a series of numerical experiments.

Fractional differential equations with dependence on the Caputo-Katugampola derivative

In this paper we present a new type of fractional operator, the Caputo–Katugampola derivative. The Caputo and the Caputo–Hadamard fractional derivatives are special cases of this new operator. An existence and uniqueness theorem for a fractional Cauchy type problem, with dependence on the Caputo–Katugampola derivative, is proven. A decomposition formula for the Caputo–Katugampola derivative is obtained. This formula allows us to provide a simple numerical procedure to solve the fractional differential equation.

Modeling some real phenomena by fractional differential equations

This paper deals with fractional differential equations, with dependence on a Caputo fractional derivative of real order. The goal is to show, based on concrete examples and experimental data from several experiments, that fractional differential equations may model more efficiently certain problems than ordinary differential equations. A numerical optimization approach based on least squares approximation is used to determine the order of the fractional operator that better describes real data, as well as other related parameters.

Similarity solutions of a Becker-Döring system with time-dependent monomer input

We formulate the Becker-Döring equations for cluster growth in the presence of a time-dependent source of monomer input. In the case of size-independent aggregation and ragmentation rate coefficients we find similarity solutions which are approached in the large time limit. The form of the solutions depends on the rate of monomer input and whether fragmentation is present in the model; four distinct types of solution are found.

Robust solutions of bi-blend recipe optimization with quadratic constraints

Production companies use raw materials to compose end-products. They often make different products with the same raw materials. In this research, the focus lies on the production of two end-products consisting of (partly) the same raw materials as cheap as possible. Each of the products has its own demand and quality requirements consisting of quadratic constraints. The minimization of the costs, given the quadratic constraints is a global optimization problem, which can be difficult because of possible local optima. Therefore, the multi modal character of the (bi-) blend problem is investigated. Standard optimization packages (solvers) in Matlab and GAMS were tested on their ability to solve the problem. In total 20 test cases were generated and taken from literature to test solvers on their effectiveness and efficiency to solve the problem. The research also gives insight in adjusting the quadratic constraints of the problem in order to make a robust problem formulation of the bi-blend problem.

Determination of local density of states in amorphous oxide semiconductors with Kelvin Probe Force Microscopy

Amorphous semiconductors are important materials as they can be deposited by physical deposition techniques on large areas and even on plastic substrates. Therefore, they are crucial for transistors in large active matrices for imaging and transparent wearable electronics. The most widely applied candidate for amorphous thin film transistors production is Indium Gallium Zinc Oxide (IGZO). It is attracting much interest because of its optical transparency, facile processing by sputtering deposition and notable improved charge carrier mobility with respect to hydrogenated amorphous silicon a-Si:H. Degradation of the device and long-term performance issues have been observed if IGZO thin film transistors are subjected to electrical stress, leading to a modification of IGZO channel properties and subthreshold slope. Therefore, it is of great interest to have a reliable and precise method to study the conduction band tail, and the density of states in amorphous semiconductors. The aim of this thesis is to develop a local technique using Kelvin Probe Force Microscopy to study the evolution of IGZO DOS properties. The work is divided into three main parts. First, solutions to the non-linear Poisson-Boltzmann equation of a metal-insulator-semiconductor junction describing the charge accumulation and its relation to DOS properties are elaborated. Second macroscopic techniques such as capacitance voltage (CV) measurements and photocurrent spectroscopy are applied to obtain a non-local estimate of band-tail DOS properties in thin film transistor samples. The third part of my my thesis is dedicated to the KPFM measurements. By fitting the data to the developed numerical model, important parameters describing the amorphous conduction band tail are obtained. The results are in excellent agreement with the macroscopic characterizations. KPFM result is comparable also with non-local optoelectronic characterizations, such as photocurrent spectroscopy.

Local and nonlocal integro-differential reaction-diffusion models for protein dynamics in Alzheimer's disease

In this work, integro-differential reaction-diffusion models are presented for the description of the temporal and spatial evolution of the concentrations of Abeta and tau proteins involved in Alzheimer's disease. Initially, a local model is analysed: this is obtained by coupling with an interaction term two heterodimer models, modified by adding diffusion and Holling functional terms of the second type. We then move on to the presentation of three nonlocal models, which differ according to the type of the growth (exponential, logistic or Gompertzian) considered for healthy proteins. In these models integral terms are introduced to consider the interaction between proteins that are located at different spatial points possibly far apart. For each of the models introduced, the determination of equilibrium points with their stability and a study of the clearance inequalities are carried out. In addition, since the integrals introduced imply a spatial nonlocality in the models exhibited, some general features of nonlocal models are presented. Afterwards, with the aim of developing simulations, it is decided to transfer the nonlocal models to a brain graph called connectome. Therefore, after setting out the construction of such a graph, we move on to the description of Laplacian and convolution operations on a graph. Taking advantage of all these elements, we finally move on to the translation of the continuous models described above into discrete models on the connectome. To conclude, the results of some simulations concerning the discrete models just derived are presented.

Effects Of Partial Inhibition Of Respiratory Complex I On H2o 2 Production By Isolated Brain Mitochondria In Different Respiratory States.

The aim of this work was to characterize the effects of partial inhibition of respiratory complex I by rotenone on H2O2 production by isolated rat brain mitochondria in different respiratory states. Flow cytometric analysis of membrane potential in isolated mitochondria indicated that rotenone leads to uniform respiratory inhibition when added to a suspension of mitochondria. When mitochondria were incubated in the presence of a low concentration of rotenone (10 nm) and NADH-linked substrates, oxygen consumption was reduced from 45.9 ± 1.0 to 26.4 ± 2.6 nmol O2 mg(-1) min(-1) and from 7.8 ± 0.3 to 6.3 ± 0.3 nmol O2 mg(-1) min(-1) in respiratory states 3 (ADP-stimulated respiration) and 4 (resting respiration), respectively. Under these conditions, mitochondrial H2O2 production was stimulated from 12.2 ± 1.1 to 21.0 ± 1.2 pmol H2O2 mg(-1) min(-1) and 56.5 ± 4.7 to 95.0 ± 11.1 pmol H2O2 mg(-1) min(-1) in respiratory states 3 and 4, respectively. Similar results were observed when comparing mitochondrial preparations enriched with synaptic or nonsynaptic mitochondria or when 1-methyl-4-phenylpyridinium ion (MPP(+)) was used as a respiratory complex I inhibitor. Rotenone-stimulated H2O2 production in respiratory states 3 and 4 was associated with a high reduction state of endogenous nicotinamide nucleotides. In succinate-supported mitochondrial respiration, where most of the mitochondrial H2O2 production relies on electron backflow from complex II to complex I, low rotenone concentrations inhibited H2O2 production. Rotenone had no effect on mitochondrial elimination of micromolar concentrations of H2O2. The present results support the conclusion that partial complex I inhibition may result in mitochondrial energy crisis and oxidative stress, the former being predominant under oxidative phosphorylation and the latter under resting respiration conditions.

STM probing of local oscillations of the Fano-Kondo effect: a Doniach-Sunjic approach for the Kondo peak

The use of the scanning tunneling microscope (STM) for the investigation of Kondo adatoms on normal metallic surfaces reveals a Fano-Kondo behavior of the conductance as a function of the tip bias. In this work, the Doniach-Sunjic expression is used to describe the Kondo peak and we analyze the effect of a complex Fano phase, arising from an external magnetic field, on the conductance pattern. It is demonstrated that such phase generates local oscillations of the Fano-Kondo line shape and can lead to the suppression of anti-resonances.

Data collapse, scaling functions, and analytical solutions of generalized growth models

We consider a nontrivial one-species population dynamics model with finite and infinite carrying capacities. Time-dependent intrinsic and extrinsic growth rates are considered in these models. Through the model per capita growth rate we obtain a heuristic general procedure to generate scaling functions to collapse data into a simple linear behavior even if an extrinsic growth rate is included. With this data collapse, all the models studied become independent from the parameters and initial condition. Analytical solutions are found when time-dependent coefficients are considered. These solutions allow us to perceive nontrivial transitions between species extinction and survival and to calculate the transition's critical exponents. Considering an extrinsic growth rate as a cancer treatment, we show that the relevant quantity depends not only on the intensity of the treatment, but also on when the cancerous cell growth is maximum.