Combining Legendre's Polynomials and Genetic Algorithm in the Solution of Nonlinear Initial-Value Problems

This work develops a method for solving ordinary differential equations, that is, initial-value problems, with solutions approximated by using Legendre's polynomials. An iterative procedure for the adjustment of the polynomial coefficients is developed, based on the genetic algorithm. This procedure is applied to several examples providing comparisons between its results and the best polynomial fitting when numerical solutions by the traditional Runge-Kutta or Adams methods are available. The resulting algorithm provides reliable solutions even if the numerical solutions are not available, that is, when the mass matrix is singular or the equation produces unstable running processes.

Theory of small-signal ac response of a dielectric liquid containing two groups of ions

The analysis of Macdonald for electrolytes is generalized to the case in which two groups of ions are present. We assume that the electrolyte can be considered as a dispersion of ions in a dielectric liquid, and that the ionic recombination can be neglected. We present the differential equations governing the ionic redistribution when the liquid is subjected to an external electric field, describing the simultaneous diffusion of the two groups of ions in the presence of their own space charge fields. We investigate the influence of the ions on the impedance spectroscopy of an electrolytic cell. In the analysis, we assume that each group of ions have equal mobility, the electrodes perfectly block and that the adsorption phenomena can be neglected. In this framework, it is shown that the real part of the electrical impedance of the cell has a frequency dependence presenting two plateaux, related to a type of ambipolar and free diffusion coefficients. The importance of the considered problem on the ionic characterization performed by means of the impedance spectroscopy technique was discussed. (c) 2008 American Institute of Physics.

Nonlinear waves in a quark gluon plasma

We study the propagation of perturbations in the quark gluon plasma. This subject has been addressed in other works and in most of the theoretical descriptions of this phenomenon the hydrodynamic equations have been linearized for simplicity. We propose an alternative approach, also based on hydrodynamics but taking into account the nonlinear terms of the equations. We show that these terms may lead to localized waves or even solitons. We use a simple equation of state for the QGP and expand the hydrodynamic equations around equilibrium configurations. The resulting differential equations describe the propagation of perturbations in the energy density. We solve them numerically and find that localized perturbations can propagate for long distances in the plasma. Under certain conditions our solutions mimic the propagation of Korteweg-de Vries solitons.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

POSITIVITY PROPERTIES OF THE FOURIER TRANSFORM AND THE STABILITY OF PERIODIC TRAVELLING-WAVE SOLUTIONS

In this paper we establish a method to obtain the stability of periodic travelling-wave solutions for equations of Korteweg-de Vries-type u(t) + u(p)u(x) - Mu(x) = 0, with M being a general pseudodifferential operator and where p >= 1 is an integer. Our approach uses the theory of totally positive operators, the Poisson summation theorem, and the theory of Jacobi elliptic functions. In particular we obtain the stability of a family of periodic travelling waves solutions for the Benjamin Ono equation. The present technique gives a new way to obtain the existence and stability of cnoidal and dnoidal waves solutions associated with the Korteweg-de Vries and modified Korteweg-de Vries equations, respectively. The theory has prospects for the study of periodic travelling-wave solutions of other partial differential equations.

Self-organized distribution of periodicity and chaos in an electrochemical oscillator

We report a detailed numerical investigation of a prototype electrochemical oscillator, in terms of high-resolution phase diagrams for an experimentally relevant section of the control (parameter) space. The prototype model consists of a set of three autonomous ordinary differential equations which captures the general features of electrochemical oscillators characterized by a partially hidden negative differential resistance in an N-shaped current-voltage stationary curve. By computing Lyapunov exponents, we provide a detailed discrimination between chaotic and periodic phases of the electrochemical oscillator. Such phases reveal the existence of an intricate structure of domains of periodicity self-organized into a chaotic background. Shrimp-like periodic regions previously observed in other discrete and continuous systems were also observed here, which corroborate the universal nature of the occurrence of such structures. In addition, we have also found a structured period distribution within the order region. Finally we discuss the possible experimental realization of comparable phase diagrams.

Mechanism and model of the oscillatory electro-oxidation of methanol

A mechanism for the kinetic instabilities observed in the galvanostatic electro-oxidation of methanol is suggested and a model developed. The model is investigated using stoichiometric network analysis as well as concepts from algebraic geometry (polynomial rings and ideal theory) revealing the occurrence of a Hopf and a saddle-node bifurcation. These analytical solutions are confirmed by numerical integration of the system of differential equations. (C) 2010 American Institute of Physics

Exercise intensity optimization for men with high cardiorespiratory fitness

Exercise intensity is a key parameter for exercise prescription but the optimal range for individuals with high cardiorespiratory fitness is unknown. The aims of this study were (1) to determine optimal heart rate ranges for men with high cardiorespiratory fitness based on percentages of maximal oxygen consumption (%VO(2max)) and reserve oxygen consumption (%VO(2reserve)) corresponding to the ventilatory threshold and respiratory compensation point, and ( 2) to verify the effect of advancing age on the exercise intensities. Maximal cardiorespiratory testing was performed on 210 trained men. Linear regression equations were calculated using paired data points between percentage of maximal heart rate (%HR(max)) and %VO(2max) and between percentage of heart rate reserve (%HRR) and %VO(2reserve) attained at each minute during the test. Values of %VO(2max) and %VO(2reserve) at the ventilatory threshold and respiratory compensation point were used to calculate the corresponding values of %HRmax and %HRR, respectively. The ranges of exercise intensity in relation to the ventilatory threshold and respiratory compensation point were achieved at 78-93% of HR(max) and 70-93% of HRR, respectively. Although absolute heart rate decreased with advancing age, there were no age-related differences in %HR(max) and %HRR at the ventilatory thresholds. Thus, in men with high cardiorespiratory fitness, the ranges of exercise intensity based on %HR(max) and %HRR regarding ventilatory threshold were 78-93% and 70-93% respectively, and were not influenced by advancing age.

Bending of stochastic Kirchhoff plates on Winkler foundations via the Galerkin method and the Askey-Wiener scheme

In this paper, the method of Galerkin and the Askey-Wiener scheme are used to obtain approximate solutions to the stochastic displacement response of Kirchhoff plates with uncertain parameters. Theoretical and numerical results are presented. The Lax-Milgram lemma is used to express the conditions for existence and uniqueness of the solution. Uncertainties in plate and foundation stiffness are modeled by respecting these conditions, hence using Legendre polynomials indexed in uniform random variables. The space of approximate solutions is built using results of density between the space of continuous functions and Sobolev spaces. Approximate Galerkin solutions are compared with results of Monte Carlo simulation, in terms of first and second order moments and in terms of histograms of the displacement response. Numerical results for two example problems show very fast convergence to the exact solution, at excellent accuracies. The Askey-Wiener Galerkin scheme developed herein is able to reproduce the histogram of the displacement response. The scheme is shown to be a theoretically sound and efficient method for the solution of stochastic problems in engineering. (C) 2009 Elsevier Ltd. All rights reserved.

Chaos-Galerkin solution of stochastic Timoshenko bending problems

This paper presents an accurate and efficient solution for the random transverse and angular displacement fields of uncertain Timoshenko beams. Approximate, numerical solutions are obtained using the Galerkin method and chaos polynomials. The Chaos-Galerkin scheme is constructed by respecting the theoretical conditions for existence and uniqueness of the solution. Numerical results show fast convergence to the exact solution, at excellent accuracies. The developed Chaos-Galerkin scheme accurately approximates the complete cumulative distribution function of the displacement responses. The Chaos-Galerkin scheme developed herein is a theoretically sound and efficient method for the solution of stochastic problems in engineering. (C) 2011 Elsevier Ltd. All rights reserved.

Galerkin Solution of Stochastic Beam Bending on Winkler Foundations

In this paper, the Askey-Wiener scheme and the Galerkin method are used to obtain approximate solutions to stochastic beam bending on Winkler foundation. The study addresses Euler-Bernoulli beams with uncertainty in the bending stiffness modulus and in the stiffness of the foundation. Uncertainties are represented by parameterized stochastic processes. The random behavior of beam response is modeled using the Askey-Wiener scheme. One contribution of the paper is a sketch of proof of existence and uniqueness of the solution to problems involving fourth order operators applied to random fields. From the approximate Galerkin solution, expected value and variance of beam displacement responses are derived, and compared with corresponding estimates obtained via Monte Carlo simulation. Results show very fast convergence and excellent accuracies in comparison to Monte Carlo simulation. The Askey-Wiener Galerkin scheme presented herein is shown to be a theoretically solid and numerically efficient method for the solution of stochastic problems in engineering.

Who should wear mask against airborne infections? Altering the contact network for controlling the spread of contagious diseases

There are several ways of controlling the propagation of a contagious disease. For instance, to reduce the spreading of an airborne infection, individuals can be encouraged to remain in their homes and/or to wear face masks outside their domiciles. However, when a limited amount of masks is available, who should use them: the susceptible subjects, the infective persons or both populations? Here we employ susceptible-infective-recovered (SIR) models described in terms of ordinary differential equations and probabilistic cellular automata in order to investigate how the deletion of links in the random complex network representing the social contacts among individuals affects the dynamics of a contagious disease. The inspiration for this study comes from recent discussions about the impact of measures usually recommended by health public organizations for preventing the propagation of the swine influenza A (H1N1) virus. Our answer to this question can be valid for other eco-epidemiological systems. (C) 2010 Elsevier BM. All rights reserved.

Network of phase-locking oscillators and a possible model for neural synchronization

In order to model the synchronization of brain signals, a three-node fully-connected network is presented. The nodes are considered to be voltage control oscillator neurons (VCON) allowing to conjecture about how the whole process depends on synaptic gains, free-running frequencies and delays. The VCON, represented by phase-locked loops (PLL), are fully-connected and, as a consequence, an asymptotically stable synchronous state appears. Here, an expression for the synchronous state frequency is derived and the parameter dependence of its stability is discussed. Numerical simulations are performed providing conditions for the use of the derived formulae. Model differential equations are hard to be analytically treated, but some simplifying assumptions combined with simulations provide an alternative formulation for the long-term behavior of the fully-connected VCON network. Regarding this kind of network as models for brain frequency signal processing, with each PLL representing a neuron (VCON), conditions for their synchronization are proposed, considering the different bands of brain activity signals and relating them to synaptic gains, delays and free-running frequencies. For the delta waves, the synchronous state depends strongly on the delays. However, for alpha, beta and theta waves, the free-running individual frequencies determine the synchronous state. (C) 2011 Elsevier B.V. All rights reserved.

On state representations of nonlinear implicit systems

This work considers a semi-implicit system A, that is, a pair (S, y), where S is an explicit system described by a state representation (x)over dot(t) = f(t, x(t), u(t)), where x(t) is an element of R(n) and u(t) is an element of R(m), which is subject to a set of algebraic constraints y(t) = h(t, x(t), u(t)) = 0, where y(t) is an element of R(l). An input candidate is a set of functions v = (v(1),.... v(s)), which may depend on time t, on x, and on u and its derivatives up to a Finite order. The problem of finding a (local) proper state representation (z)over dot = g(t, z, v) with input v for the implicit system Delta is studied in this article. The main result shows necessary and sufficient conditions for the solution of this problem, under mild assumptions on the class of admissible state representations of Delta. These solvability conditions rely on an integrability test that is computed from the explicit system S. The approach of this article is the infinite-dimensional differential geometric setting of Fliess, Levine, Martin, and Rouchon (1999) (`A Lie-Backlund Approach to Equivalence and Flatness of Nonlinear Systems`, IEEE Transactions on Automatic Control, 44(5), (922-937)).

On the basic reproduction number and the topological properties of the contact network: An epidemiological study in mainly locally connected cellular automata

We study the spreading of contagious diseases in a population of constant size using susceptible-infective-recovered (SIR) models described in terms of ordinary differential equations (ODEs) and probabilistic cellular automata (PCA). In the PCA model, each individual (represented by a cell in the lattice) is mainly locally connected to others. We investigate how the topological properties of the random network representing contacts among individuals influence the transient behavior and the permanent regime of the epidemiological system described by ODE and PCA. Our main conclusions are: (1) the basic reproduction number (commonly called R(0)) related to a disease propagation in a population cannot be uniquely determined from some features of transient behavior of the infective group; (2) R(0) cannot be associated to a unique combination of clustering coefficient and average shortest path length characterizing the contact network. We discuss how these results can embarrass the specification of control strategies for combating disease propagations. (C) 2009 Elsevier B.V. All rights reserved.