Modelagem Matemática e Análise de Modelos Matemáticos na Educação Matemática

In this paper we present two studies, the first one completed and the second one in development, which are based in teaching approaches that propose the qualitative study of mathematical models as a strategy for the teaching and learning of mathematical concepts. These teaching approaches focus on subjects from Higher Education such as Introduction to Ordinary Differential Equations and Topics of Differential and Integral Calculus. We denominate this common aspect of the teaching approaches as Model Analysis and in a preliminary level we relate it with Mathematical Modeling. Furthermore, we discuss some questions related with the choice of the theme and the role of Digital Technologies when Model Analysis is applied.

On an energy exchange process and appearance of chaos in a non-ideal portal frame dynamical system

In this paper, we consider non-ideal excitation devices such as DC motors with restrictenergy output capacity. When such motors are attached to structures which needexcitation power levels similar to the source power capacity, jump phenomena and theincrease in power required near resonance characterize the Sommerfeld Effect, actingas a sort of an energy sink. One of the problems often faced by designers of suchstructures is how to drive the system through resonance and avoid this energy sink.Our basic structural model is a simple portal frame driven by a num-ideal powersource-(NIPF). We also investigate the absorption of resonant vibrations (nonlinearand chaotic) by means of a nonlinear sub-structure known as a Nonlinear Energy Sink(NES). An energy exchange process between the NIPF and NES in the passagethrough resonance is investigated, as well the suppression of chaos.

Estudos de técnicas de reinicialização em métodos level-set e suas aplicações em escoamentos bifásicos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Equações diferenciais ordinárias não suaves autônomas e não autônomas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Equações diferenciais ordinárias não suaves autônomas e não autônomas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

On estimating the basic reproduction number in distinct stages of a contagious disease spreading

In epidemiology, the basic reproduction number R-0 is usually defined as the average number of new infections caused by a single infective individual introduced into a completely susceptible population. According to this definition. R-0 is related to the initial stage of the spreading of a contagious disease. However, from epidemiological models based on ordinary differential equations (ODE), R-0 is commonly derived from a linear stability analysis and interpreted as a bifurcation parameter: typically, when R-0 >1, the contagious disease tends to persist in the population because the endemic stationary solution is asymptotically stable: when R-0 <1, the corresponding pathogen tends to naturally disappear because the disease-free stationary solution is asymptotically stable. Here we intend to answer the following question: Do these two different approaches for calculating R-0 give the same numerical values? In other words, is the number of secondary infections caused by a unique sick individual equal to the threshold obtained from stability analysis of steady states of ODE? For finding the answer, we use a susceptibleinfective-recovered (SIR) model described in terms of ODE and also in terms of a probabilistic cellular automaton (PCA), where each individual (corresponding to a cell of the PCA lattice) is connected to others by a random network favoring local contacts. The values of R-0 obtained from both approaches are compared, showing good agreement. (C) 2012 Elsevier B.V. All rights reserved.

Continuity of Dynamical Structures for Nonautonomous Evolution Equations Under Singular Perturbations

In this paper we study the continuity of invariant sets for nonautonomous infinite-dimensional dynamical systems under singular perturbations. We extend the existing results on lower-semicontinuity of attractors of autonomous and nonautonomous dynamical systems. This is accomplished through a detailed analysis of the structure of the invariant sets and its behavior under perturbation. We prove that a bounded hyperbolic global solutions persists under singular perturbations and that their nonlinear unstable manifold behave continuously. To accomplish this, we need to establish results on roughness of exponential dichotomies under these singular perturbations. Our results imply that, if the limiting pullback attractor of a nonautonomous dynamical system is the closure of a countable union of unstable manifolds of global bounded hyperbolic solutions, then it behaves continuously (upper and lower) under singular perturbations.

Leptogenesis via hypermagnetic fields and baryon asymmetry

We study baryon asymmetry generation originated from the leptogenesis in the presence of hypermagnetic fields in the early Universe plasma before the electroweak phase I ransition (EWPT). For the simplest Chern-Simons (CS) wave configuration of hypermagnetic field we find the baryon asymmetry growth when the hypermagnetic field value changes due to alpha(2)-dynamo and the lepton asymmetry rises due to the Abelian anomaly. We solve the corresponding integro-differential equations for the lepton asymmetries describing such selfconsistent dynamics for lepto- and baryogenesis in the two scenarios: (i) when a primordial lepton asymmetry sits in right electrons e(R); and (ii) when, in addition to e(R), a left lepton asyninwtty for e(L) and v(eL) at due to chirality flip reactions provided by in Iiigg,s decays at the temperatures, T < T-RL similar to 10 TeV. We find that the baryon asymmetry of the Universe (BAU) rises very fast through such leptogenesis, especially, in strong hypermagnetic fields. Varying (decreasing) the CS wave number parameter k(0) < 10(-7) T-EW one can recover the observable value of BAU, eta(B) similar to 10(-9), where k(0) = 10(-7) T-EW corresponds to the ataxinittat value for CS wave number surviving ohmic dissipation of hypermagnetic field. In the scenario (ii) one predicts the essential difference of the lepton numbers of right- and left electrons at EWPT time, L-eR - L-eL similar to (mu(eR) / mu(eL))/T-EW = Delta mu/T-EW similar or equal to 10(-5) that can be used as an initial condition for chiral asymmetry after EWPT.

Period adding cascades: Experiment and modeling in air bubbling

Period adding cascades have been observed experimentally/numerically in the dynamics of neurons and pancreatic cells, lasers, electric circuits, chemical reactions, oceanic internal waves, and also in air bubbling. We show that the period adding cascades appearing in bubbling from a nozzle submerged in a viscous liquid can be reproduced by a simple model, based on some hydrodynamical principles, dealing with the time evolution of two variables, bubble position and pressure of the air chamber, through a system of differential equations with a rule of detachment based on force balance. The model further reduces to an iterating one-dimensional map giving the pressures at the detachments, where time between bubbles come out as an observable of the dynamics. The model has not only good agreement with experimental data, but is also able to predict the influence of the main parameters involved, like the length of the hose connecting the air supplier with the needle, the needle radius and the needle length. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3695345]

STRUCTURE AND BIFURCATION OF PULLBACK ATTRACTORS IN A NON-AUTONOMOUS CHAFEE-INFANTE EQUATION

The Chafee-Infante equation is one of the canonical infinite-dimensional dynamical systems for which a complete description of the global attractor is available. In this paper we study the structure of the pullback attractor for a non-autonomous version of this equation, u(t) = u(xx) + lambda(xx) - lambda u beta(t)u(3), and investigate the bifurcations that this attractor undergoes as A is varied. We are able to describe these in some detail, despite the fact that our model is truly non-autonomous; i.e., we do not restrict to 'small perturbations' of the autonomous case.

SUFFICIENT CONDITIONS ON DIFFUSIVITY FOR THE EXISTENCE AND NONEXISTENCE OF STABLE EQUILIBRIA WITH NONLINEAR FLUX ON THE BOUNDARY

A reaction-diffusion equation with variable diffusivity and non-linear flux boundary condition is considered. The goal is to give sufficient conditions on the diffusivity function for nonexistence and also for existence of nonconstant stable stationary solutions. Applications are given for the main result of nonexistence.

The tangential variation of a localized flux-type eigenvalue problem

In this work the differentiability of the principal eigenvalue lambda = lambda(1)(Gamma) to the localized Steklov problem -Delta u + qu = 0 in Omega, partial derivative u/partial derivative nu = lambda chi(Gamma)(x)u on partial derivative Omega, where Gamma subset of partial derivative Omega is a smooth subdomain of partial derivative Omega and chi(Gamma) is its characteristic function relative to partial derivative Omega, is shown. As a key point, the flux subdomain Gamma is regarded here as the variable with respect to which such differentiation is performed. An explicit formula for the derivative of lambda(1) (Gamma) with respect to Gamma is obtained. The lack of regularity up to the boundary of the first derivative of the principal eigenfunctions is a further intrinsic feature of the problem. Therefore, the whole analysis must be done in the weak sense of H(1)(Omega). The study is of interest in mathematical models in morphogenesis. (C) 2011 Elsevier Inc. All rights reserved.

Topics in geometry, analysis and inverse problems

The thesis consists of three independent parts. Part I: Polynomial amoebas We study the amoeba of a polynomial, as de ned by Gelfand, Kapranov and Zelevinsky. A central role in the treatment is played by a certain convex function which is linear in each complement component of the amoeba, which we call the Ronkin function. This function is used in two di erent ways. First, we use it to construct a polyhedral complex, which we call a spine, approximating the amoeba. Second, the Monge-Ampere measure of the Ronkin function has interesting properties which we explore. This measure can be used to derive an upper bound on the area of an amoeba in two dimensions. We also obtain results on the number of complement components of an amoeba, and consider possible extensions of the theory to varieties of codimension higher than 1. Part II: Differential equations in the complex plane We consider polynomials in one complex variable arising as eigenfunctions of certain differential operators, and obtain results on the distribution of their zeros. We show that in the limit when the degree of the polynomial approaches innity, its zeros are distributed according to a certain probability measure. This measure has its support on the union of nitely many curve segments, and can be characterized by a simple condition on its Cauchy transform. Part III: Radon transforms and tomography This part is concerned with different weighted Radon transforms in two dimensions, in particular the problem of inverting such transforms. We obtain stability results of this inverse problem for rather general classes of weights, including weights of attenuation type with data acquisition limited to a 180 degrees range of angles. We also derive an inversion formula for the exponential Radon transform, with the same restriction on the angle.

Partial melting within the lower crust – Constraints from bulk rock geochemical data and trace element mineral analyses of granulites and migmatites from Central Finland

This PhD thesis concerns geochemical constraints on recycling and partial melting of Archean continental crust. A natural example of such processes was found in the Iisalmi area of Central Finland. The rocks from this area are Middle to Late Archean in age and experienced metamorphism and partial melting between 2.7-2.63 Ga. The work is based on extensive field work. It is furthermore founded on bulk rock geochemical data as well as in-situ analyses of minerals. All geochemical data were obtained at the Institute of Geosciences, University of Mainz using X-ray fluorescence, solution ICP-MS and laser ablation-ICP-MS for bulk rock geochemical analyses. Mineral analyses were accomplished by electron microprobe and laser ablation ICP-MS. Fluid inclusions were studied by microscope on a heating-freezing-stage at the Geoscience Center, University Göttingen. Part I focuses on the development of a new analytical method for bulk rock trace element determination by laser ablation-ICP-MS using homogeneous glasses fused from rock powder on an Iridium strip heater. This method is applicable for mafic rock samples whose melts have low viscosities and homogenize quickly at temperatures of ~1200°C. Highly viscous melts of felsic samples prevent melting and homogenization at comparable temperatures. Fusion of felsic samples can be enabled by addition of MgO to the rock powder and adjustment of melting temperature and melting duration to the rock composition. Advantages of the fusion method are low detection limits compared to XRF analyses and avoidance of wet-chemical processing and use of strong acids as in solution ICP-MS as well as smaller sample volumes compared to the other methods. Part II of the thesis uses bulk rock geochemical data and results from fluid inclusion studies for discrimination of melting processes observed in different rock types. Fluid inclusion studies demonstrate a major change in fluid composition from CO2-dominated fluids in granulites to aqueous fluids in TTG gneisses and amphibolites. Partial melts were generated in the dry, CO2-rich environment by dehydration melting reactions of amphibole which in addition to tonalitic melts produced the anhydrous mineral assemblages of granulites (grt + cpx + pl ± amph or opx + cpx + pl + amph). Trace element modeling showed that mafic granulites are residues of 10-30 % melt extraction from amphibolitic precursor rocks. The maximum degree of melting in intermediate granulites was ~10 % as inferred from modal abundances of amphibole, clinopyroxene and orthopyroxene. Carbonic inclusions are absent in upper-amphibolite facies migmatites whereas aqueous inclusion with up to 20 wt% NaCl are abundant. This suggests that melting within TTG gneisses and amphibolites took place in the presence of an aqueous fluid phase that enabled melting at the wet solidus at temperatures of 700-750°C. The strong disruption of pre-metamorphic structures in some outcrops suggests that the maximum amount of melt in TTG gneisses was ~25 vol%. The presence of leucosomes in all rock types is taken as the principle evidence for melt formation. However, mineralogical appearance as well as major and trace element composition of many leucosomes imply that leucosomes seldom represent frozen in-situ melts. They are better considered as remnants of the melt channel network, e.g. ways on which melts escaped from the system. Part III of the thesis describes how analyses of minerals from a specific rock type (granulite) can be used to determine partition coefficients between different minerals and between minerals and melt suitable for lower crustal conditions. The trace element analyses by laser ablation-ICP-MS show coherent distribution among the principal mineral phases independent of rock composition. REE contents in amphibole are about 3 times higher than REE contents in clinopyroxene from the same sample. This consistency has to be taken into consideration in models of lower crustal melting where amphibole is replaced by clinopyroxene in the course of melting. A lack of equilibrium is observed between matrix clinopyroxene / amphibole and garnet porphyroblasts which suggests a late stage growth of garnet and slow diffusion and equilibration of the REE during metamorphism. The data provide a first set of distribution coefficients of the transition metals (Sc, V, Cr, Ni) in the lower crust. In addition, analyses of ilmenite and apatite demonstrate the strong influence of accessory phases on trace element distribution. Apatite contains high amounts of REE and Sr while ilmenite incorporates about 20-30 times higher amounts of Nb and Ta than amphibole. Furthermore, trace element mineral analyses provide evidence for magmatic processes such as melt depletion, melt segregation, accumulation and fractionation as well as metasomatism having operated in this high-grade anatectic area.

Vibrational spectroscopic and electrochemical investigations of multi-centered heme proteins in biomimetic membrane architectures

Membrane proteins play a major role in every living cell. They are the key factors in the cell’s metabolism and in other functions, for example in cell-cell interaction, signal transduction, and transport of ions and nutrients. Cytochrome c oxidase (CcO), as one of the membrane proteins of the respiratory chain, plays a significant role in the energy transformation of higher organisms. CcO is a multi centered heme protein, utilizing redox energy to actively transport protons across the mitochondrial membrane. One aim of this dissertation is to investigate single steps in the mechanism of the ion transfer process coupled to electron transfer, which are not fully understood. The protein-tethered bilayer lipid membrane is a general approach to immobilize membrane proteins in an oriented fashion on a planar electrode embedded in a biomimetic membrane. This system enables the combination of electrochemical techniques with surface enhanced resonance Raman (SERRS), surface enhanced reflection absorption infrared (SEIRAS), and surface plasmon spectroscopy to study protein mediated electron and ion transport processes. The orientation of the enzymes within the surface confined architecture can be controlled by specific site-mutations, i.e. the insertion of a poly-histidine tag to different subunits of the enzyme. CcO can, thus, be oriented uniformly with its natural electron pathway entry pointing either towards or away from the electrode surface. The first orientation allows an ultra-fast direct electron transfer(ET) into the protein, not provided by conventional systems, which can be leveraged to study intrinsic charge transfer processes. The second orientation permits to study the interaction with its natural electron donor cytochrome c. Electrochemical and SERR measurements show conclusively that the redox site structure and the activity of the surface confined enzyme are preserved. Therefore, this biomimetic system offers a unique platform to study the kinetics of the ET processes in order to clarify mechanistic properties of the enzyme. Highly sensitive and ultra fast electrochemical techniques allow the separation of ET steps between all four redox centres including the determination of ET rates. Furthermore, proton transfer coupled to ET could be directly measured and discriminated from other ion transfer processes, revealing novel mechanistic information of the proton transfer mechanism of cytochrome c oxidase. In order to study the kinetics of the ET inside the protein, including the catalytic center, time resolved SEIRAS and SERRS measurements were performed to gain more insight into the structural and coordination changes of the heme environment. The electrical behaviour of tethered membrane systems and membrane intrinsic proteins as well as related charge transfer processes were simulated by solving the respective sets of differential equations, utilizing a software package called SPICE. This helps to understand charge transfer processes across membranes and to develop models that can help to elucidate mechanisms of complex enzymatic processes.