scholarly journals Parametric Excitation and Evolutionary Dynamics

2013 ◽  
Vol 80 (5) ◽  
Author(s):  
Rocio E. Ruelas ◽  
David G. Rand ◽  
Richard H. Rand

Parametric excitation refers to dynamics problems in which the forcing function enters into the governing differential equation as a variable coefficient. Evolutionary dynamics refers to a mathematical model of natural selection (the “replicator” equation) which involves a combination of game theory and differential equations. In this paper we apply perturbation theory to investigate parametric resonance in a replicator equation having periodic coefficients. In particular, we study evolution in the Rock-Paper-Scissors game, which has biological and social applications. Here periodic coefficients could represent seasonal variation. We show that 2:1 subharmonic resonance can destabilize the usual “Rock-Paper-Scissors” equilibrium for parameters located in a resonant tongue in parameter space. However, we also show that the tongue may be absent or very small if the forcing parameters are chosen appropriately.

Author(s):  
Rocio E Ruelas ◽  
David G Rand ◽  
Richard H Rand

Nonlinear parametric excitation refers to the nonlinear analysis of a system of ordinary differential equations with periodic coefficients. In contrast to linear parametric excitation, which offers determinations of the stability of equilibria, nonlinear parametric excitation has as its goal the structure of the phase space, as given by a portrait of the Poincare map. In this article, perturbation methods and numerical integration are applied to the replicator equation with periodic coefficients, being a model from evolutionary game theory where evolutionary dynamics are added to classical game theory using differential equations. In particular, we study evolution in the Rock–Paper–Scissors game, which has biological and social applications. Here, periodic coefficients could represent seasonal variation.


2014 ◽  
Vol 23 (3-4) ◽  
pp. 109-128 ◽  
Author(s):  
Ivo Senjanović ◽  
Marko Tomić ◽  
Neven Hadžić

AbstractThe Timoshenko beam theory, which deals with the deflection and rotation in two partial differential equations of motion, is transformed into a single partial differential equation with pure bending deflection as a potential function for the determination of the total deflection, rotation angle, and sectional forces. Inclusion of a nonlocal stress parameter results in the extension of the governing differential equation from the 4th to the 6th order. A simply supported nanotube is considered, and the governing differential equation is decomposed into a system of ordinary differential equations by employing the modal superposition method, separation of variables, and the Galerkin method. Both moving nanoparticle gravity and inertia force are consistently taken into account, resulting in ordinary and parametric excitation, respectively. As a novelty, the parameters are split into a constant and a time-dependent part. The former is added to the ordinary system of equations, which is solved analytically in the frequency domain by the harmonic balance method, while the system with variable coefficients is solved in the time domain by the perturbation method. The effects of slenderness ratio, nonlocal parameter, stiffness of elastic medium, nanoparticle gravity and inertia force, and velocity on the free and forced nanotube response are also investigated. Special attention is paid to the influence of damping on resonance. Performed parametric analysis is physically transparent due to the obtained semi-analytical solution. Some analytical results of illustrative examples are compared with numerical ones from the relevant literature, and notable differences are discussed.


2006 ◽  
Vol 11 (1) ◽  
pp. 13-32 ◽  
Author(s):  
B. Bandyrskii ◽  
I. Lazurchak ◽  
V. Makarov ◽  
M. Sapagovas

The paper deals with numerical methods for eigenvalue problem for the second order ordinary differential operator with variable coefficient subject to nonlocal integral condition. FD-method (functional-discrete method) is derived and analyzed for calculating of eigenvalues, particulary complex eigenvalues. The convergence of FD-method is proved. Finally numerical procedures are suggested and computational results are schown.


The propagation properties of linear wave motions in magnetic and/or velocity shears which vary in one coordinate z (say) are usually governed by a second order linear ordinary differential equation in the independent variable z. It is proved that associated with any such differential equation there always exists a quantity A which is independent of z. By employing A a measure of the intensity of the wave, this result is used to investigate the general propagation properties of hydromagnetic-gravity waves (e.g. critical level absorption, valve effects and wave amplification) in magnetic and/or velocity shears, using a full wave treatment. When variations in the basic state are included, the governing differential equation usually has more singularities than it has in the W.K.B.J. approximation, which neglects all variations in the background state. The study of a wide variety of models shows that critical level behaviour occurs only at the singularities predicted by the W.K.B.J. approximation. Although the solutions of the differential equation are necessarily singular at the irregularities whose presence is solely due to the inclusion of variations in the basic state, the intensity of the wave (as measured by A) is continuous there. Also the valve effect is found to persist whatever the relation between the wavelength of the wave and the scale of variations of the background state. In addition, it is shown that a hydromagnetic-gravity wave incident upon a finite magnetic and/or velocity shear can be amplified (or over-reflected) in the absence of any critical levels within the shear layer. In a Boussinesq fluid rotating uniformly about the vertical, wave amplification can occur if the horizontal vertically sheared flow and magnetic field are perpendicular. In a compressible isothermal fluid, on the other hand, wave amplification not only occurs in both magnetic-velocity and velocity shears but also in a magnetic shear acting alone.


2016 ◽  
Vol 6 (1) ◽  
pp. 19 ◽  
Author(s):  
Ahmad Salah Edeen Nassef ◽  
Mohammed A. Dahim

<p class="1Body">This paper was investigating the buckling problem of reinforced concrete columns considering the reinforced concrete as bi – modular material. Governing differential equations was driven. The relation between the non-dimensional transverse deflection and non-dimensional distance between centroid axis and the neutral axis "eccentricity" was drawn to enable the solution of the governing differential equation. The new approach was verified with different experimental results and different codes of practice.<strong></strong></p>


2012 ◽  
Vol 2012 ◽  
pp. 1-13
Author(s):  
P. G. Siddheshwar ◽  
A. Chan ◽  
U. S. Mahabaleswar

The paper discusses the boundary layer flow of Walters’ liquid B over a stretching sheet. The stretching is assumed to be a quadratic function of the coordinate along the direction of stretching. The study encompasses within its realm both Walters’ liquid B and second order liquid. The velocity distribution is obtained by solving the nonlinear governing differential equation. Analytical expressions are obtained for stream function and velocity components as functions of the viscoelastic and stretching related parameters. It is shown that the viscoelasticity goes hand in hand with quadratic stretching in enhancing the lifting of the liquid as we go along the sheet.


2016 ◽  
Vol 78 (6-5) ◽  
Author(s):  
Nurul Akmal Mohamed ◽  
Nur Fadhilah Ibrahim ◽  
Mohd Rozni Md Yusof ◽  
Nurul Farihan Mohamed ◽  
Nurul Huda Mohamed

In this paper, we present the numerical results of the Boundary-Domain Integro-Differential Equation (BDIDE) associated to Dirichlet problem for an elliptic type Partial Differential Equation (PDE) with a variable coefficient. The numerical constructions are based on discretizing the boundary of the problem region by utilizing continuous linear iso-parametric elements while the domain of the problem region is meshed by using iso-parametric quadrilateral bilinear domain elements. We also use a semi-analytic method to handle the integration that exhibits logarithmic singularity instead of using Gauss-Laguare quadrature formula. The numerical results that employed the semi-analytic method give better accuracy as compared to those when we use Gauss-Laguerre quadrature formula. The system of equations that obtained by the discretized BDIDE is solved by an iterative method (Neumann series expansion) as well as a direct method (LU decomposition method). From our numerical experiments on all test domains, the relative errors of the solutions when applying semi-analytic method are smaller than when we use Gauss-Laguerre quadrature formula for the integration with logarithmic singularity. Unlike Dirichlet Boundary Integral Equation (BIE), the spectral properties of the Dirichlet BDIDE is not known. The Neumann iterations will converge to the solution if and only if the spectral radius of matrix operator is less than 1. In our numerical experiment on all the test domains, the Neumann series does converge. It gives some conclusions for the spectral properties of the Dirichlet BDIDE even though more experiments on the general Dirichlet problems need to be carried out.


Sign in / Sign up

Export Citation Format

Share Document