Equações diferenciais implícitas com descontinuidade

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

Sobre a estabilidade estrutural e regularizações de campos de vetores descontínuos

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

Soluções quase automórficas para equações diferenciais abstratas de segunda ordem

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

Cálculo fracionário aplicado em dinâmica tumoral: método da Transformada Diferencial generalizada

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)

The effect of flex fuel vehicles in the Brazilian light road transportation

The retaking of the ethanol program in the year 2003 as a fuel for light road transportation in Brazil through the introduction of flex fuel vehicles fleet was a good strategy to overcome the difficulties of the ethanol production sector and did work to increase its market share relative to gasoline. This process, however, may cause a future disequilibrium on the food production and on the refining oil derivates structure. In order to analyze the substitution process resultant of the competition between two opponents fighting for the same market, in this case the gasoline/ethanol substitution process, a method derived from the biomathematics based on the non-linear differential equations (NLDE) system is utilized. A brief description of the method is presented. Numerical adherence of the method to explain several substitution phenomena that occurred in the past is presented in the previous author`s paper, in which the urban gas pipeline system substitution of bottled LPG in the dwelling sector and the substitution of the urban diesel transportation fleet by compressed natural gas (CNG) buses is presented. The proposed method is particularly suitable for prospective analysis and scenarios assessment. (c) 2008 Elsevier Ltd. All rights reserved.

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.

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.

Influence of the magnetization damping on dynamic hysteresis loops in single domain particles

This article reports on the influence of the magnetization damping on dynamic hysteresis loops in single-domain particles with uniaxial anisotropy. The approach is based on the Neel-Brown theory and the hierarchy of differential recurrence relations, which follow from averaging over the realizations of the stochastic Landau-Lifshitz equation. A new method of solution is proposed, where the resulting system of differential equations is solved directly using optimized algorithms to explore its sparsity. All parameters involved in uniaxial systems are treated in detail, with particular attention given to the frequency dependence. It is shown that in the ferromagnetic resonance region, novel phenomena are observed for even moderately low values of the damping. The hysteresis loops assume remarkably unusual shapes, which are also followed by a pronounced reduction of their heights. Also demonstrated is that these features remain for randomly oriented ensembles and, moreover, are approximately independent of temperature and particle size. (C) 2012 American Institute of Physics. [doi:10.1063/1.3684629]

Delay nonlinear boundary conditions as limit of reactions concentrating in the boundary

In this work, we are interested in the dynamic behavior of a parabolic problem with nonlinear boundary conditions and delay in the boundary. We construct a reaction-diffusion problem with delay in the interior, where the reaction term is concentrated in a neighborhood of the boundary and this neighborhood shrinks to boundary, as a parameter epsilon goes to zero. We analyze the limit of the solutions of this concentrated problem and prove that these solutions converge in certain continuous function spaces to the unique solution of the parabolic problem with delay in the boundary. This convergence result allows us to approximate the solution of equations with delay acting on the boundary by solutions of equations with delay acting in the interior and it may contribute to analyze the dynamic behavior of delay equations when the delay is at the boundary. (C) 2012 Elsevier Inc. All rights reserved.

Globally solvable systems of complex vector fields

We consider a class of involutive systems of n smooth vector fields on the n + 1 dimensional torus. We obtain a complete characterization for the global solvability of this class in terms of Liouville forms and of the connectedness of all sublevel and superlevel sets of the primitive of a certain 1-form in the minimal covering space.

Modeling of continuous thermal processing of a non-Newtonian liquid food under diffusive laminar flow in a tubular system

The classic conservative approach for thermal process design can lead to over-processing, especially for laminar flow, when a significant distribution of temperature and of residence time occurs. In order to optimize quality retention, a more comprehensive model is required. A model comprising differential equations for mass and heat transfer is proposed for the simulation of the continuous thermal processing of a non-Newtonian food in a tubular system. The model takes into account the contribution from heating and cooling sections, the heat exchange with the ambient air and effective diffusion associated with non-ideal laminar flow. The study case of soursop juice processing was used to test the model. Various simulations were performed to evaluate the effect of the model assumptions. An expressive difference in the predicted lethality was observed between the classic approach and the proposed model. The main advantage of the model is its flexibility to represent different aspects with a small computational time, making it suitable for process evaluation and design. (C) 2012 Elsevier Ltd. All rights reserved.

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]

STABLE PIECEWISE POLYNOMIAL VECTOR FIELDS

Let N = {y > 0} and S = {y < 0} be the semi-planes of R-2 having as common boundary the line D = {y = 0}. Let X and Y be polynomial vector fields defined in N and S, respectively, leading to a discontinuous piecewise polynomial vector field Z = (X, Y). This work pursues the stability and the transition analysis of solutions of Z between N and S, started by Filippov (1988) and Kozlova (1984) and reformulated by Sotomayor-Teixeira (1995) in terms of the regularization method. This method consists in analyzing a one parameter family of continuous vector fields Z(epsilon), defined by averaging X and Y. This family approaches Z when the parameter goes to zero. The results of Sotomayor-Teixeira and Sotomayor-Machado (2002) providing conditions on (X, Y) for the regularized vector fields to be structurally stable on planar compact connected regions are extended to discontinuous piecewise polynomial vector fields on R-2. Pertinent genericity results for vector fields satisfying the above stability conditions are also extended to the present case. A procedure for the study of discontinuous piecewise vector fields at infinity through a compactification is proposed here.