Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

Parâmetro de regularização em problemas inversos: estudo numérico com a transformada de Radon

In general, an inverse problem corresponds to find a value of an element x in a suitable vector space, given a vector y measuring it, in some sense. When we discretize the problem, it usually boils down to solve an equation system f(x) = y, where f : U Rm ! Rn represents the step function in any domain U of the appropriate Rm. As a general rule, we arrive to an ill-posed problem. The resolution of inverse problems has been widely researched along the last decades, because many problems in science and industry consist in determining unknowns that we try to know, by observing its effects under certain indirect measures. Our general subject of this dissertation is the choice of Tykhonov´s regulaziration parameter of a poorly conditioned linear problem, as we are going to discuss on chapter 1 of this dissertation, focusing on the three most popular methods in nowadays literature of the area. Our more specific focus in this dissertation consists in the simulations reported on chapter 2, aiming to compare the performance of the three methods in the recuperation of images measured with the Radon transform, perturbed by the addition of gaussian i.i.d. noise. We choosed a difference operator as regularizer of the problem. The contribution we try to make, in this dissertation, mainly consists on the discussion of numerical simulations we execute, as is exposed in Chapter 2. We understand that the meaning of this dissertation lays much more on the questions which it raises than on saying something definitive about the subject. Partly, for beeing based on numerical experiments with no new mathematical results associated to it, partly for being about numerical experiments made with a single operator. On the other hand, we got some observations which seemed to us interesting on the simulations performed, considered the literature of the area. In special, we highlight observations we resume, at the conclusion of this work, about the different vocations of methods like GCV and L-curve and, also, about the optimal parameters tendency observed in the L-curve method of grouping themselves in a small gap, strongly correlated with the behavior of the generalized singular value decomposition curve of the involved operators, under reasonably broad regularity conditions in the images to be recovered

Modelling the lethargic crab disease.

The lethargic crab disease (LCD) is an emergent infirmity that has decimated native populations of the mangrove land crab (Ucides cordatus, Decapoda: Ocypodidae) along the Brazilian coast. Several potential etiological agents have been linked with LCD, but only in 2005 was it proved that it is caused by an ascomycete fungus. This is the first attempt to develop a mathematical model to describe the epidemiological dynamics of LCD. The model presents four possible scenarios, namely, the trivial equilibrium, the disease-free equilibrium, endemic equilibrium, and limit cycles arising from a Hopf bifurcation. The threshold values depend on the basic reproductive number of crabs and fungi, and on the infection rate. These scenarios depend on both the biological assumptions and the temporal evolution of the disease. Numerical simulations corroborate the analytical results and illustrate the different temporal dynamics of the crab and fungus populations.

Population dynamics, life stage and ecological modeling in Chrysomya albiceps (Wiedemann) (Diptera: Calliphoridae)

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

Optimisation of a continuous flash fermentation for butanol production using the response surface methodology

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Time series analysis for minority game simulations of financial markets

The minority game (MG) model introduced recently provides promising insights into the understanding of the evolution of prices, indices and rates in the financial markets. In this paper we perform a time series analysis of the model employing tools from statistics, dynamical systems theory and stochastic processes. Using benchmark systems and a financial index for comparison, several conclusions are obtained about the generating mechanism for this kind of evolution. The motion is deterministic, driven by occasional random external perturbation. When the interval between two successive perturbations is sufficiently large, one can find low dimensional chaos in this regime. However, the full motion of the MG model is found to be similar to that of the first differences of the SP500 index: stochastic, nonlinear and (unit root) stationary. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Vortices in nonlocal Gross-Pitaevskii equation

We consider vortices in the nonlocal two-dimensional Gross-Pitaevskii equation with the interaction potential having Lorentz-shaped dependence on the relative momentum. It is shown that in the Fourier series expansion with respect to the polar angle, the unstable modes of the axial n-fold vortex have orbital numbers l satisfying 0 < \l\ < 2\n\, as in the local model. Numerical simulations show that nonlocality slightly decreases the threshold rotation frequency above which the nonvortex state ceases to be the global energy minimum and decreases the frequency of the anomalous mode of the 1-vortex. In the case of higher axial vortices, nonlocality leads to instability against splitting with the creation of antivortices and gives rise to additional anomalous modes with higher orbital numbers. Despite new instability channels with the creation of antivortices, for a stationary solution comprised of vortices and antivortices there always exists another vortex solution, composed solely of vortices, with the same total vorticity but with a lower energy.

Mathematical glimpse on the Y chromosome degeneration

The Y chromosomes are genetically degenerate and do not recombine with their matching partners X. Non-recombination of XY pairs has been pointed out as the key factor for the degeneration of the Y chromosome. The aim here is to show that there is a mathematical asymmetry in sex chromosomes which leads to the degeneration of Y chromosomes even in the absence of XX and XY recombination. A model for sex-chromosome evolution in a stationary regime is proposed. The consequences of their asymmetry are analyzed and lead us to a couple of conclusions. First, Y chromosome degeneration shows up v 2 more often than X chromosome degeneration. Second, if nature prohibits female mortalities from beeing exactly 50%, then Y chromosome degeneration is inevitable.

Mathematical models of generalized diffusion

We discuss in this paper equations describing processes involving non-linear and higher-order diffusion. We focus on a particular case (u(t) = 2 lambda (2)(uu(x))(x) + lambda (2)u(xxxx)), which is put into analogy with the KdV equation. A balance of nonlinearity and higher-order diffusion enables the existence of self-similar solutions, describing diffusive shocks. These shocks are continuous solutions with a discontinuous higher-order derivative at the shock front. We argue that they play a role analogous to the soliton solutions in the dispersive case. We also discuss several physical instances where such equations are relevant.

NOISE and ULTRAVIOLET DIVERGENCES IN SIMULATIONS of GINZBURG-LANDAU-LANGEVIN TYPE of EQUATIONS

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

SPH simulations of clumps formation by dissipative collisions of molecular clouds - II. Magnetic case

We performed computer simulations of interstellar cloud-cloud collisions using the three-dimensional smoothed particle magnetohydrodynamics method. In order to study the role of the magnetic field on the process of collision-triggered fragmentation, we focused our attention on head-on supersonic collisions between two identical spherical molecular-clouds. Two extreme configurations of the magnetic field were adopted: parallel and perpendicular to the initial clouds motion. The initial magnetic field strength was approximately 12.0 muG. In the parallel case, much more of the collision debris were retained in the shocking region than in the non-magnetic case where gas escaped freely throughout the symmetry plane. Differently from the non-magnetic case, eddy-like vortices were formed. The regions of highest vorticity and the the regions of highest density are offset. We found clumps formation only in the parallel case, however, they were larger, hotter and less dense than in the analogous non-magnetic case. In the perpendicular case, the compressed field works as a magnetic wall, preventing a stronger compression of the colliding clouds. This last effect inhibits direct contact of the two clouds. In both cases, we found that the field lines show a chaotic aspect in large scales. Also, the field magnitude is considerably amplified in the shock layer. However, the field distribution is almost coherent in the higher density regions.

SPH simulations of clumps formation by dissipative collision of molecular clouds - I. Non magnetic case

Computer experiments of interstellar cloud collisions were performed with a new smoothed-particle-hydrodynamics (SPH) code. The SPH quantities were calculated by using spatially adaptive smoothing lengths and the SPH fluid equations of motion were solved by means of a hierarchical multiple time-scale leapfrog. Such a combination of methods allows the code to deal with a large range of hydrodynamic quantities. A careful treatment of gas cooling by H, H(2), CO and H II, as well as a heating mechanism by cosmic rays and by H(2) production on grains surface, were also included in the code. The gas model reproduces approximately the typical environment of dark molecular clouds. The experiments were performed by impinging two dynamically identical spherical clouds onto each other with a relative velocity of 10 km s(-1) but with a different impact parameter for each case. Each object has an initial density profile obeying an r(-1)-law with a cutoff radius of 10 pc and with an initial temperature of 20 K. As a main result, cloud-cloud collision triggers fragmentation but in expense of a large amount of energy dissipated, which occurred in the head-on case only. Off-center collision did not allow remnants to fragment along the considered time (similar to 6 Myr). However, it dissipated a considerable amount of orbital energy. Structures as small as 0.1 pc, with densities of similar to 10(4) cm(-3), were observed in the more energetic collision.

Short comments on self-synchronization of two non-ideal sources supported by a flexible portal frame structure

A practical problem of synchronization of a non-ideal (i.e. when the excitation is influenced by the response of the system) and non-linear vibrating system was posed and investigated by means of numerical simulations. Two rotating unbalanced motors compose the mathematical model considered here with limited power supply mounted on the horizontal beam of a simple portal frame. As a starting point, the problem is reduced to a four-degrees-of-freedom model and its equations of motion, derived elsewhere via a Lagrangian approach, are presented. The numerical results show the expected phenomena associated with the passage through resonance with limited power. Further, for a two-to-one relationship between the frequencies associated with the first symmetric mode and the sway mode, by using the variation of torque constants, the control of the self-synchronization and synchronization (in the system) are observed at certain levels of excitations.

On non-linear dynamics and control designs applied to the ideal and non-ideal variants of the Fitzhugh-Nagumo (FN) mathematical model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)