The ambivalent effect of lattice structure on a spatial game

The evolution of cooperation is studied in lattice-structured populations, in which each individual who adopts one of the following strategies ‘always defect’ (ALLD), ‘tit-for-tat’ (TFT), and ‘always cooperate’ (ALLC) plays the repeated Prisoner’s Dilemma game with its neighbors according to an asynchronous update rule. Computer simulations are applied to analyse the dynamics depending on major parameters. Mathematical analyses based on invasion probability analysis, mean-field approximation, as well as pair approximation are also used. We find that the lattice structure promotes the evolution of cooperation compared with a non-spatial population, this is also confirmed by invasion probability analysis in one dimension. Meanwhile, it also inhibits the evolution of cooperation due to the advantage of being spiteful, which indicates the key role of specific life-history assumptions. Mean-field approximation fails to predict the outcome of computer simulations. Pair approximation is accurate in two dimensions but fails in one dimension.     

Efficient and stable perfectly matched layer for CEM

An efficient unsplit perfectly matched layer for numerical simulation of electromagnetic waves in unbounded domains is derived via a complex change of variables. In order to surround a Cartesian grid with the PML, the time-dependent PML requires only one (scalar) auxiliary variable in two space dimensions and six (scalar) auxiliary variables in three space dimensions. It is therefore cheap and straightforward to implement. We use Fourier and energy methods to prove the stability of the PML. We extend the stability result to a semi-discrete PML approximated by central finite differences of arbitrary order of accuracy and to a fully discrete problem for the ‘Leap-Frog’ schemes. This makes precise the usefulness of the derived PML model for longtime simulations. Numerical experiments are presented, illustrating the accuracy and stability of the PML.     

Spatial‐skin effect for eigenvibrations of a thick cascade junction with ‘heavy’ concentrated masses

A spectral problem for the Laplace operator in a thick cascade junction with concentrated masses is considered. This cascade junction consists of the junction's body and a great number of ?‐alternating thin rods belonging to two classes. One class consists of rods of finite length, and the second one consists of rods of small length of order . The density of the junction is of order on the rods from the second class and outside of them. The asymptotic behavior of eigenvalues and eigenfunctions of this problem is studied as ? → 0. There exist five qualitatively different cases in the asymptotic behavior of eigenmagnitudes as ? → 0, namely the case of ‘light’ concentrated (α ∈ (0,1)), ‘middle’ concentrated (α = 1), and ‘heavy’ concentrated masses (α ∈ (1, + ∞ )) that we divide into ‘slightly heavy’ concentrated (α ∈ (1,2)), ‘intermediate heavy’ concentrated (α = 2), and ‘very heavy’ concentrated masses (α > 2). In the paper, we study in detail the influence of the concentrated masses on the asymptotic behavior if α ∈ (1,2). We construct the leading terms of asymptotic expansions both for the eigenvalues and eigenfunctions and prove the corresponding asymptotic estimates. Copyright © 2013 John Wiley & Sons, Ltd.     

Spectral asymptotics and basis properties of fourth order differential operators with regular boundary conditions

In this paper, we consider the problem where λ is a spectral parameter; q(x) ∈ L₁(0,1) is complex‐valued function; α_s, s = 1,2,3, are arbitrary complex constants that satisfy α₂ = α₁ + α₃ and σ = 0,1. The boundary conditions of this problem are regular, but not strongly regular. Asymptotic formulae for eigenvalues and eigenfunctions of the considered boundary value problem are established. It is proved that all the eigenvalues, except for finite number, are simple and the system of root functions of this spectral problem forms a basis in the space L_p(0,1), 1 < p < ∞ , when ; moreover, this basis is unconditional for p = 2. We note that the considered problem was previously investigated in the condition of α₂ ≠ α₁ + α₃. Copyright © 2013 John Wiley & Sons, Ltd.     

Steady flow in a deformable porous medium

A nonlinear model for a steady flow in a deformable porous medium is considered. The flow is governed by the poroelasticity system consisting of an elasticity equation for the displacement of the porous medium and Darcy's equation for the pressure in the fluid. This poroelasticity system is nonlinear when the permeability in Darcy's equation is assumed to depend on the dilatation of the porous medium. Existence and uniqueness of a weak solution of this poroelasticity system is established under rather weak assumptions on the regularity of the data. Convergence of a finite element approximation is proved and verified through numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.     

Projective synchronization between different fractional‐order hyperchaotic systems with uncertain parameters using proposed modified adaptive projective synchronization technique

In the present article, the authors have proposed a modified projective adaptive synchronization technique for fractional‐order chaotic systems. The adaptive projective synchronization controller and identification parameters law are developed on the basis of Lyapunov direct stability theory. The proposed method is successfully applied for the projective synchronization between fractional‐order hyperchaotic Lü system as drive system and fractional‐order hyperchaotic Lorenz chaotic system as response system. A comparison between the effects on synchronization time due to the presence of fractional‐order time derivatives for modified projective synchronization method and proposed modified adaptive projective synchronization technique is the key feature of the present article. Numerical simulation results, which are carried out using Adams–Boshforth–Moulton method show that the proposed technique is effective, convenient and also faster for projective synchronization of fractional‐order nonlinear dynamical systems. Copyright © 2013 John Wiley & Sons, Ltd.     

A fractional order epidemic model for the simulation of outbreaks of influenza A(H1N1)

In this paper, we propose a nonlinear fractional order model in order to explain and understand the outbreaks of influenza A(H1N1). In the fractional model, the next state depends not only upon its current state but also upon all of its historical states. Thus, the fractional model is more general than the classical epidemic models. In order to deal with the fractional derivatives of the model, we rely on the Caputo operator and on the Grünwald–Letnikov method to numerically approximate the fractional derivatives. We conclude that the nonlinear fractional order epidemic model is well suited to provide numerical results that agree very well with real data of influenza A(H1N1) at the level population. In addition, the proposed model can provide useful information for the understanding, prediction, and control of the transmission of different epidemics worldwide. Copyright © 2013 John Wiley & Sons, Ltd.     

Reflection of plane waves by rough surfaces in the sense of Born approximation

The topic of the present paper is the reflection of electromagnetic plane waves by rough surfaces, that is, by smooth and bounded perturbations of planar faces. Moreover, the contrast between the cover material and the substrate beneath the rough surface is supposed to be low. In this case, a modification of Stearns’ formula based on Born approximation and Fourier techniques is derived for a special class of surfaces. This class contains the graphs of functions where the interface function is a radially modulated almost periodic function. For the Born formula to converge, a sufficient and almost necessary condition is given. A further technical condition is defined, which guarantees the existence of the corresponding far field of the Born approximation. This far field contains plane waves, far‐field terms such as those for bounded scatterers, and, additionally, a new type of terms. The derived formulas can be used for the fast numerical computations of far fields and for the statistics of random rough surfaces. Copyright © 2013 John Wiley & Sons, Ltd.     

Formulação numérica de estruturas compósitas amortecidas utilizando as teorias da Deformação Cisalhante de Primeira Ordem e de Alta Ordem

Composite materials have been used in the design of the aircrafts structures because their low weight and high mechanical strength. However, structures made in composite material are exposed to dynamical and/or static loading environments. Therefore, a major research effort is undertaken in the development of tools numerical for analysis and design of composite structures. This paper presents a numerical formulation of the composite structures using the Finite Element Method (FEM). The damped composite structures, using inserted viscoelastic devices, and undamped composite structures are formulated by FEM. Viscoelastic materials are applied as continuous layers inserted on composite structures. The intrinsic damping of the composite material is included in the studies, too. The First‐order (FSDT) and Higher‐order Shear Deformation (HSDT) theories are formulated. They are distinguished by order of the approximation functions used in the mechanical displacements field. Both theories are computationally implemented using the Serendipity finite element. This is a rectangular finite element with 8 nodes, 5 or 11 degrees of freedom per node. The results are compared with papers predictions. The advantages and disadvantages of using each theory in the modeling of composite (thin or thick) and thick sandwiches structures, including the intrinsic and the viscoelastic damping, are discusses.