Eigenfunctions and fundamental solutions of the fractional Laplace and Dirac operators: the Riemann-Liouville case

In this paper we study eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator $\Delta_+^{(\alpha,\beta,\gamma)}:= D_{x_0^+}^{1+\alpha} +D_{y_0^+}^{1+\beta} +D_{z_0^+}^{1+\gamma},$ where $(\alpha, \beta, \gamma) \in \,]0,1]^3$, and the fractional derivatives $D_{x_0^+}^{1+\alpha}$, $D_{y_0^+}^{1+\beta}$, $D_{z_0^+}^{1+\gamma}$ are in the Riemann-Liouville sense. Applying operational techniques via two-dimensional Laplace transform we describe a complete family of eigenfunctions and fundamental solutions of the operator $\Delta_+^{(\alpha,\beta,\gamma)}$ in classes of functions admitting a summable fractional derivative. Making use of the Mittag-Leffler function, a symbolic operational form of the solutions is presented. From the obtained family of fundamental solutions we deduce a family of fundamental solutions of the fractional Dirac operator, which factorizes the fractional Laplace operator. We apply also the method of separation of variables to obtain eigenfunctions and fundamental solutions.

Characterizing transport through a crowded environment with different obstacle sizes

Transport through crowded environments is often classified as anomalous, rather than classical, Fickian diffusion. Several studies have sought to describe such transport processes using either a continuous time random walk or fractional order differential equation. For both these models the transport is characterized by a parameter α, where α = 1 is associated with Fickian diffusion and α < 1 is associated with anomalous subdiffusion. Here, we simulate a single agent migrating through a crowded environment populated by impenetrable, immobile obstacles and estimate α from mean squared displacement data. We also simulate the transport of a population of such agents through a similar crowded environment and match averaged agent density profiles to the solution of a related fractional order differential equation to obtain an alternative estimate of α. We examine the relationship between our estimate of α and the properties of the obstacle field for both a single agent and a population of agents; we show that in both cases, α decreases as the obstacle density increases, and that the rate of decrease is greater for smaller obstacles. Our work suggests that it may be inappropriate to model transport through a crowded environment using widely reported approaches including power laws to describe the mean squared displacement and fractional order differential equations to represent the averaged agent density profiles.

Modeling transport through an environment crowded by a mixture of obstacles of different shapes and sizes

Many biological environments are crowded by macromolecules, organelles and cells which can impede the transport of other cells and molecules. Previous studies have sought to describe these effects using either random walk models or fractional order diffusion equations. Here we examine the transport of both a single agent and a population of agents through an environment containing obstacles of varying size and shape, whose relative densities are drawn from a specified distribution. Our simulation results for a single agent indicate that smaller obstacles are more effective at retarding transport than larger obstacles; these findings are consistent with our simulations of the collective motion of populations of agents. In an attempt to explore whether these kinds of stochastic random walk simulations can be described using a fractional order diffusion equation framework, we calibrate the solution of such a differential equation to our averaged agent density information. Our approach suggests that these kinds of commonly used differential equation models ought to be used with care since we are unable to match the solution of a fractional order diffusion equation to our data in a consistent fashion over a finite time period.

Modeling, control and simulation of full-power converter wind turbines equipped with permanent magnet synchronous generator

In this paper, two wind turbines equipped with a permanent magnet synchronous generator (PMSG) and respectively with a two-level or a multilevel converter are simulated in order to access the malfunction transient performance. Three different drive train mass models, respectively, one, two and three mass models, are considered in order to model the bending flexibility of the blades. Moreover, a fractional-order control strategy is studied comparatively to a classical integer-order control strategy. Computer simulations are carried out, and conclusions about the total harmonic distortion (THD) of the electric current injected into the electric grid are in favor of the fractional-order control strategy.

Comparative study of power converter topologies and control strategies for the harmonic performance of variable-speed wind turbine generator systems

Power converters play a vital role in the integration of wind power into the electrical grid. Variable-speed wind turbine generator systems have a considerable interest of application for grid connection at constant frequency. In this paper, comprehensive simulation studies are carried out with three power converter topologies: matrix, two-level and multilevel. A fractional-order control strategy is studied for the variable-speed operation of wind turbine generator systems. The studies are in order to compare power converter topologies and control strategies. The studies reveal that the multilevel converter and the proposed fractional-order control strategy enable an improvement in the power quality, in comparison with the other power converters using a classical integer-order control strategy. (C) 2010 Elsevier Ltd. All rights reserved.

Transient analysis of variable-speed wind turbines at wind speed disturbances and a pitch control malfunction

As wind power generation undergoes rapid growth, new technical challenges emerge: dynamic stability and power quality. The influence of wind speed disturbances and a pitch control malfunction on the quality of the energy injected into the electric grid is studied for variable-speed wind turbines with different power-electronic converter topologies. Additionally, a new control strategy is proposed for the variable-speed operation of wind turbines with permanent magnet synchronous generators. The performance of disturbance attenuation and system robustness is ascertained. Simulation results are presented and conclusions are duly drawn. (C) 2010 Elsevier Ltd. All rights reserved.

Analysis of temperature time-series: embedding dynamics into the MDS method

Global warming and the associated climate changes are being the subject of intensive research due to their major impact on social, economic and health aspects of the human life. Surface temperature time-series characterise Earth as a slow dynamics spatiotemporal system, evidencing long memory behaviour, typical of fractional order systems. Such phenomena are difficult to model and analyse, demanding for alternative approaches. This paper studies the complex correlations between global temperature time-series using the Multidimensional scaling (MDS) approach. MDS provides a graphical representation of the pattern of climatic similarities between regions around the globe. The similarities are quantified through two mathematical indices that correlate the monthly average temperatures observed in meteorological stations, over a given period of time. Furthermore, time dynamics is analysed by performing the MDS analysis over slices sampling the time series. MDS generates maps describing the stations’ locus in the perspective that, if they are perceived to be similar to each other, then they are placed on the map forming clusters. We show that MDS provides an intuitive and useful visual representation of the complex relationships that are present among temperature time-series, which are not perceived on traditional geographic maps. Moreover, MDS avoids sensitivity to the irregular distribution density of the meteorological stations.

Dynamical modelling of a genetic algorithm

This work addresses the signal propagation and the fractional-order dynamics during the evolution of a genetic algorithm (GA). In order to investigate the phenomena involved in the GA population evolution, the mutation is exposed to excitation perturbations during some generations and the corresponding fitness variations are evaluated. Three distinct fitness functions are used to study their influence in the GA dynamics. The input and output signals are studied revealing a fractional-order dynamic evolution, characteristic of a long-term system memory.

New findings on the dynamics of HIV and TB coinfection models

In this paper we study a model for HIV and TB coinfection. We consider the integer order and the fractional order versions of the model. Let α∈[0.78,1.0] be the order of the fractional derivative, then the integer order model is obtained for α=1.0. The model includes vertical transmission for HIV and treatment for both diseases. We compute the reproduction number of the integer order model and HIV and TB submodels, and the stability of the disease free equilibrium. We sketch the bifurcation diagrams of the integer order model, for variation of the average number of sexual partners per person and per unit time, and the tuberculosis transmission rate. We analyze numerical results of the fractional order model for different values of α, including α=1. The results show distinct types of transients, for variation of α. Moreover, we speculate, from observation of the numerical results, that the order of the fractional derivative may behave as a bifurcation parameter for the model. We conclude that the dynamics of the integer and the fractional order versions of the model are very rich and that together these versions may provide a better understanding of the dynamics of HIV and TB coinfection.

Analysis of Forest Fires by means of Pseudo Phase Plane and Multidimensional Scaling Methods

Forest fires dynamics is often characterized by the absence of a characteristic length-scale, long range correlations in space and time, and long memory, which are features also associated with fractional order systems. In this paper a public domain forest fires catalogue, containing information of events for Portugal, covering the period from 1980 up to 2012, is tackled. The events are modelled as time series of Dirac impulses with amplitude proportional to the burnt area. The time series are viewed as the system output and are interpreted as a manifestation of the system dynamics. In the first phase we use the pseudo phase plane (PPP) technique to describe forest fires dynamics. In the second phase we use multidimensional scaling (MDS) visualization tools. The PPP allows the representation of forest fires dynamics in two-dimensional space, by taking time series representative of the phenomena. The MDS approach generates maps where objects that are perceived to be similar to each other are placed on the map forming clusters. The results are analysed in order to extract relationships among the data and to better understand forest fires behaviour.

The Persistence of Memory

This paper analyzes several natural and man-made complex phenomena in the perspective of dynamical systems. Such phenomena are often characterized by the absence of a characteristic length-scale, long range correlations and persistent memory, which are features also associated to fractional order systems. For each system, the output, interpreted as a manifestation of the system dynamics, is analyzed by means of the Fourier transform. The amplitude spectrum is approximated by a power law function and the parameters are interpreted as an underlying signature of the system dynamics. The complex systems under analysis are then compared in a global perspective in order to unveil and visualize hidden relationships among them.

Sonogashira coupling on an extended gold surface in vacuo: reaction of phenylacetylene with iodobenzene on Au(111)

Temperature-programmed reaction measurements supported by scanning tunneling microscopy have shown that phenylacetylene and iodobenzene react on smooth Au(111) under vacuum conditions to yield biphenyl and diphenyldiacetylene, the result of homocoupling of the reactant molecules. They also produce diphenylacetylene, the result of Sonogashira cross-coupling, prototypical of a class of reactions that are of paramount importance in synthetic organic chemistry and whose mechanism remains controversial. Roughened Au(111) is completely inert toward all three reactions, indicating that the availability of crystallographically well-defined adsorption sites is crucially important. High-resolution X-ray photoelectron spectroscopy and near-edge X-ray absorption fine structure spectroscopy show that the reactants are initially present as intact, essentially flat-lying molecules and that the temperature threshold for Sonogashira coupling coincides with that for C−I bond scission in the iodobenzene reactant. The fractional-order kinetics and low temperature associated with desorption of the Sonogashira product suggest that the reaction occurs at the boundaries of islands of adsorbed reactants and that its appearance in the gas phase is rate-limited by the surface reaction. These findings demonstrate unambiguously and for the first time that this heterogeneous cross-coupling chemistry is an intrinsic property of extended, metallic pure gold surfaces: no other species, including solvent molecules, basic or charged (ionic) species are necessary to mediate the process.

System identification algorithms for the analysis of dielectric responses from broadband spectroscopies

We discuss the modeling of dielectric responses for an electromagnetically excited network of capacitors and resistors using a systems identification framework. Standard models that assume integral order dynamics are augmented to incorporate fractional order dynamics. This enables us to relate more faithfully the modeled responses to those reported in the Dielectrics literature.

Estratégias de controle de ordem fracionária aplicadas ao amortecimento de oscilações eletromecânicas em sistemas elétricos de potência

Neste trabalho é proposta uma nova metodologia de projeto de estabilizadores de sistemas de potência baseada em teoria de sistemas de ordem fracionária (ESP-OF). A estratégia é baseada em uma generalização do projeto de compensadores do tipo rede avançoatraso (lead-lag) para o domínio de funções de transferência de ordem fracionária. Uma nova variável de projeto, a qual define a ordem da dinâmica fracionária do controlador, é sintonizada para se obter um compromisso entre um bom desempenho no amortecimento do modo eletromecânico dominante e uma robustez ampliada do ESP-OF. O desempenho do ESP-OF foi avaliado experimentalmente, em um sistema de potência em escala reduzida, localizado no Laboratório de Sistemas de Potência da Universidade Federal do Pará. A referida planta de teste apresenta uma estrutura típica do tipo gerador síncrono conectado a um barramento infinito e exibe um modo dominante de oscilação eletromecânica, de amortecimento extremamente reduzido, cujo valor da frequência natural é em torno de 1,2 Hz. O ESP-OF foi então projetado para ampliar o amortecimento relativo desse modo alvo, para toda a faixa de operação admissível. Para fins de implementação prática, primeiramente foram realizados testes experimentais para a identificação de um modelo nominal da planta, sob a forma de uma função de transferência pulsada, para uso na fase de projeto. O modelo obtido experimentalmente foi então validado e posteriormente utilizado tanto para o projeto do ESP-OF quanto para o projeto de um ESP convencional (utilizado para fins de comparação de desempenho). As leis de controle amortecedor do ESP-OF foram calculadas, convertidas para a forma de equações a diferenças e, subsequentemente, embarcadas em sistema digital baseado em microcontrolador DSPIC. Diversos testes de resposta ao impulso foram realizadas sob diferentes condições operacionais. As respectivas respostas dinâmicas dos sinais de saída da planta (desvio de potencia ativa) e do esforço de controle foram registradas para fins de análise. Os resultados experimentais mostraram que o ESP fracionário apresentou um desemprenho dinâmico e robustez similar em comparação com o desempenho obtido por um ESP convencional, para toda a faixa de operação investigada.

Soluções numéricas de Equações Diferenciais Fracionária e Ordinária para o cálculo da distribuição da concentração em um sistema de compartimento

The purpose of this work was the study of numerical methods for differential equations of fractional order and ordinary. These methods were applied to the problem of calculating the distribution of the concentration of a given substance over time in a given physical system. The two compartment model was used for representation of this system. Comparison between numerical solutions obtained were performed and, in particular, also compared with the analytical solution of this problem. Finally, estimates for the error between the solutions were calculated