Unsteady helical flows of Oldroyd-B fluids

The unsteady helical flow of an Oldroyd-B fluid, in an infinite circular cylinder, is studied by using finite Hankel transforms. The motion is produced by the cylinder that, at time t = 0⁺, is subject to torsional and longitudinal time-dependent shear stresses. The solutions that have been obtained, presented under series form, satisfy all imposed initial and boundary conditions. The corresponding solutions for Maxwell, second grade and Newtonian fluids are obtained as limiting cases of general solutions. Finally, the influence of the pertinent parameters on the fluid motion is underlined by graphical illustrations.     

Exact solutions for the rotational flow of a generalized Maxwell fluid between two circular cylinders

Unsteady flow of an incompressible generalized Maxwell fluid between two coaxial circular cylinders is studied by means of the Laplace and finite Hankel transforms. The motion of the fluid is produced by the rotation of cylinders around their common axis. The solutions that have been obtained, written in integral and series form in terms of the generalized G_a,b,c(·, t)-functions, are presented as a sum of the Newtonian solutions and the corresponding non-Newtonian contributions. They satisfy all imposed initial and boundary conditions and for λ → 0 reduce to the solutions corresponding to the Newtonian fluids performing the same solution. Furthermore, the corresponding solutions for ordinary Maxwell fluids are also obtained for β = 1. Finally, in order to reveal some relevant physical aspects of the obtained results, the diagrams of the velocity field ω(r, t) have been depicted against r and t for different values of the material and fractional parameters.     

Unsteady longitudinal flow of a generalized Oldroyd-B fluid in cylindrical domains

Considering a fractional derivative model the unsteady flow of an Oldroyd-B fluid between two infinite coaxial circular cylinders is studied by using finite Hankel and Laplace transforms. The motion is produced by the inner cylinder which is subject to a time dependent longitudinal shear stress at time t = 0⁺. The solution obtained under series form in terms of generalized G and R functions, satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Oldroyd-B, generalized and ordinary Maxwell, and Newtonian fluids are obtained as limiting cases of our general solutions. The influence of pertinent parameters on the fluid motion as well as a comparison between models is illustrated graphically.     

Fractional dynamics of systems with long-range interaction

We consider one-dimensional chain of coupled linear and nonlinear oscillators with long-range powerwise interaction defined by a term proportional to 1/∣n − m∣^α+1. Continuous medium equation for this system can be obtained in the so-called infrared limit when the wave number tends to zero. We construct a transform operator that maps the system of large number of ordinary differential equations of motion of the particles into a partial differential equation with the Riesz fractional derivative of order α, when 0 < α < 2. Few models of coupled oscillators are considered and their synchronized states and localized structures are discussed in details. Particularly, we discuss some solutions of time-dependent fractional Ginzburg–Landau (or nonlinear Schrodinger) equation.     

Analysis of projectile motion in view of fractional calculus

The projectile motion is examined by means of the fractional calculus. The fractional differential equations of the projectile motion are introduced by generalizing Newton’s second law and Caputo’s fractional derivative is considered to use the physical initial conditions. In the absence of air resistance it is found that at certain conditions, the range and the maximum height of the projectile obtained by using the fractional calculus give the same results of the classical calculus. It is also found that, unlike the classical projectile motion, the launching angle that maximizes the horizontal range is a function of the arbitrary order of the fractional derivative α. Moreover, in a resistant medium, the parametric equations are expressed in terms of Mittag-Leffler function and the results agree with those of the classical projectile as α → 2. Moreover, the trajectories of the projectile are discussed in graphs and compared with those of the classical calculus. In order to explore the validity of modelling the projectile motion by the fractional approach, we compared our results with the experimental data of mortar.     

The dynamics of systems with large gyroscopic forces and the realization of constraints

Lagrangian systems with a large multiplier N on the gyroscopic terms are considered. Simplified equations of motion of general form with holonomic constraints are obtained in the first approximation with respect to the small parameter ɛ = 1/N. The structure of the solutions of the precessional equations is examined.     

Some new exact analytical solutions for helical flows of second grade fluids

The helical flow of a second grade fluid, between two infinite coaxial circular cylinders, is studied using Laplace and finite Hankel transforms. The motion of the fluid is due to the inner cylinder that, at time t = 0⁺ begins to rotate around its axis, and to slide along the same axis due to hyperbolic sine or cosine shear stresses. The components of the velocity field and the resulting shear stresses are presented in series form in terms of Bessel functions J₀(•), Y₀(•), J₁(•), Y₁(•), J₂(•) and Y₂(•). The solutions that have been obtained satisfy all imposed initial and boundary conditions and are presented as a sum of large-time and transient solutions. Furthermore, the solutions for Newtonian fluids performing the same motion are also obtained as special cases of general solutions. Finally, the solutions that have been obtained are compared and the influence of pertinent parameters on the fluid motion is discussed. A comparison between second grade and Newtonian fluids is analyzed by graphical illustrations.     

A cartographic study of the phase space of the elliptic restricted three body problem. Application to the Sun–Jupiter–Asteroid system

This work deals with numerical investigations of the phase space of the planar elliptic restricted three body model. The Sun–Jupiter–Asteroid system is considered and the fast Lyapunov indicator (FLI) is used as a tool to examine various types of orbits on which the infinitesimal mass can undergo. The FLI is computed on given grids of initial conditions regularly spaced in the domain 1.5 AU ? a ? 6 AU and 0 ? e ? 0.5 and for various choices of initial angles: the argument of perihelion ω and mean anomaly M. On the obtained charts the stability regions, the chaotic zones and the geography of resonances are clearly displayed. Moreover, the ‘V’ shaped layers associated with the mean motion resonances of low order with its chaotic zones due to separatrix splitting and libration regions are clearly distinguished. Their size is discussed as a function of the resonance order and the parameters entering into the perturbing function. The results are discussed and compared with analytical studies concerning the subject.     

An approximate solution of nonlinear fractional reaction–diffusion equation

The article presents a mathematical model of nonlinear reaction diffusion equation with fractional time derivative α (0 < α ? 1) in the form of a rapidly convergent series with easily computable components. Fractional reaction diffusion equation is used for modeling of merging travel solutions in nonlinear system for popular dynamics. The fractional derivatives are described in the Caputo sense. The anomalous behaviors of the nonlinear problems in the form of sub- and super-diffusion due to the presence of reaction term are shown graphically for different particular cases.     

Wave propagation in nonlocal elastic continua modelled by a fractional calculus approach

In the present paper, the wave propagation in one-dimensional elastic continua, characterized by nonlocal interactions modeled by fractional calculus, is investigated. Spatial derivatives of non-integer order 1 < α < 2 are involved in the governing equation, which is solved by fractional finite differences. The influence of long-range interactions is then analyzed as α varies: the resonant frequencies and the standing waves of a nonlocal bar are evaluated and the deviations from the classical (local) ones, recovered by imposing α = 2, are discussed.     

Dynamics of landslide model with time delay and periodic parameter perturbations

In present paper, we analyze the dynamics of a single-block model on an inclined slope with Dieterich–Ruina friction law under the variation of two new introduced parameters: time delay T_d and initial shear stress μ. It is assumed that this phenomenological model qualitatively simulates the motion along the infinite creeping slope. The introduction of time delay is proposed to mimic the memory effect of the sliding surface and it is generally considered as a function of history of sliding. On the other hand, periodic perturbation of initial shear stress emulates external triggering effect of long-distant earthquakes or some non-natural vibration source. The effects of variation of a single observed parameter, T_d or μ, as well as their co-action, are estimated for three different sliding regimes: β < 1, β = 1 and β > 1, where β stands for the ratio of long-term to short-term stress changes. The results of standard local bifurcation analysis indicate the onset of complex dynamics for very low values of time delay. On the other side, numerical approach confirms an additional complexity that was not observed by local analysis, due to the possible effect of global bifurcations. The most complex dynamics is detected for β < 1, with a complete Ruelle–Takens–Newhouse route to chaos under the variation of T_d, or the co-action of both parameters T_d and μ. These results correspond well with the previous experimental observations on clay and siltstone with low clay fraction. In the same regime, the perturbation of only a single parameter, μ, renders the oscillatory motion of the block. Within the velocity-independent regime, β = 1, the inclusion and variation of T_d generates a transition to equilibrium state, whereas the small oscillations of μ induce oscillatory motion with decreasing amplitude. The co-action of both parameters, in the same regime, causes the decrease of block’s velocity. As for β > 1, highly-frequent, limit-amplitude oscillations of initial stress give rise to oscillatory motion. Also for β > 1, in case of perturbing only the initial shear stress, with smaller amplitude, velocity of the block changes exponentially fast. If the time delay is introduced, besides the stress perturbation, within the same regime, the co-action of T_d (T_d < 0.1) and small oscillations of μ induce the onset of deterministic chaos.     

The stochastic soliton-like solutions of stochastic KdV equations

By means of a generalized method and symbolic computation, we consider a stochastic KdV equation U_t + f(t)U ♢ U_x + g(t)U_xxx = W(t) ♢ R^♢(t, U, U_x, U_xxx). We construct new and more general formal solutions. At the same time, we recover all the solutions found by Xie [Phys. Lett. A 310 (2003) 161]. The solutions obtained include the nontravelling wave and coefficient function’s stochastic soliton-like solutions, singular stochastic soliton-like solutions, stochastic triangular functions solutions.     

New exact solution profiles of the nonlinear dispersion Drinfel’d–Sokolov (D(m, n)) system

In this paper, to understand the role of nonlinear dispersion in the coupled systems, the nonlinear dispersion Drinfel’d–Sokolov system (called D(m, n) system) is investigated. As a consequence, many types of compacton and solitary pattern solutions are obtained. Moreover, some solitary wave solutions are also deduced for differential parameters m, n. When n = 1, the D(m, 1) system with linear dispersion is shown to possess also compacton and solitary pattern solutions, which contain the known results. Moreover, some rational solutions of D(m, n) system are also deduced.     

Free convection heat and mass transfer in a power law fluid past an inclined surface with thermophoresis

The problem of heat and mass transfer in a power law, two-dimensional, laminar, boundary layer flow of a viscous incompressible fluid over an inclined plate with heat generation and thermophoresis is investigated by the characteristic function method. The governing non-linear partial differential equations describing the flow and heat transfer problem are transformed into a set of coupled non-linear ordinary differential equation which was solved using Runge–Kutta shooting method. Exact solutions for the dimensionless temperature and concentration profiles, are presented graphically for different physical parameters and for the different power law exponents 0 < n < 0.5 and for n > 0.5.     

Backbone fractal dimension and fractal hybrid orbital of protein structure

Fractal geometry analysis provides a useful and desirable tool to characterize the configuration and structure of proteins. In this paper we examined the fractal properties of 750 folded proteins from four different structural classes, namely (1) the α-class (dominated by α-helices), (2) the β-class (dominated by β-pleated sheets), (3) the (α/β)-class (α-helices and β-sheets alternately mixed) and (4) the (α + β)-class (α-helices and β-sheets largely segregated) by using two fractal dimension methods, i.e. “the local fractal dimension” and “the backbone fractal dimension” (a new and useful quantitative parameter). The results showed that the protein molecules exhibit a fractal behavior in the range of 1 ? N ? 15 (N is the number of the interval between two adjacent amino acid residues), and the value of backbone fractal dimension is distinctly greater than that of local fractal dimension for the same protein. The average value of two fractal dimensions decreased in order of α > α/β > α + β > β. Moreover, the mathematical formula for the hybrid orbital model of protein based on the concept of backbone fractal dimension is in good coincidence with that of the similarity dimension. So it is a very accurate and simple method to analyze the hybrid orbital model of protein by using the backbone fractal dimension.     

Homotopy analysis of unsteady boundary layer flow adjacent to permeable stretching surface in a porous medium

This communication deals with the unsteady boundary layer flow of a viscous fluid in porous medium started due to the impulsively stretching of the plane wall. The wall is assumed to be porous so that suction or injection is possible. Complete analytic solution which is uniformly valid for all the dimensionless times 0 ⩽ τ < 0 in the whole spatial region 0 ⩽ η < ∞ is obtained by a purely analytic technique, namely the homotopy analysis method. Results are discussed through graphs.     

New exact solitary-wave special solutions for the nonlinear dispersive K(m, n) equations

In this paper, we study the nonlinear dispersive K(m, n) equations: u_t + (u^m)_x − (uⁿ)_xxx = 0 which exhibit solutions with solitary patterns. New exact solitary solutions are found. The two special cases, K(2, 2) and K(3, 3), are chosen to illustrate the concrete features of the decomposition method in K(m, n) equations. The nonlinear equations K(m, n) are studied for two different cases, namely when m = n being odd and even integers. General formulas for the solutions of K(m, n) equations are established.     

Comment on the 3+1 dimensional Kadomtsev–Petviashvili equations

We comment on traveling wave solutions and rational solutions to the 3+1 dimensional Kadomtsev–Petviashvili (KP) equations: (u_t + 6uu_x + u_xxx)_x ± 3u_yy ± 3u_zz = 0. We also show that both of the 3+1 dimensional KP equations do not possess the three-soliton solution. This suggests that none of the 3+1 dimensional KP equations should be integrable, and partially explains why they do not pass the Painlevé test. As by-products, the one-soliton and two-soliton solutions and four classes of specific three-soliton solutions are explicitly presented.     

Principal resonance responses of SDOF systems with small fractional derivative damping under narrow-band random parametric excitation

The principal resonance responses of nonlinear single-degree-of-freedom (SDOF) systems with lightly fractional derivative damping of order α (0 < α < 1) subject to the narrow-band random parametric excitation are investigated. The method of multiple scales is developed to derive two first order stochastic differential equation of amplitude and phase, and then to examine the influences of fractional order and intensity of random excitation on the first-order and second-order moment. As an example, the stochastic Duffing oscillator with fractional derivative damping is considered. The effects of detuning frequency parameter, the intensity of random excitation and the fractional order derivative damping on stability are studied through the largest Lyapunov exponent. The corresponding theoretical results are well verified through direct numerical simulations. In addition, the phenomenon of stochastic jump is analyzed for parametric principal resonance responses via finite differential method. The stochastic jump phenomena indicates that the most probable motion is around the larger non-trivial branch of the amplitude response when the intensity of excitation is very small, and the probable motion of amplitude responses will move from the larger non-trivial branch to trivial branch with the increasing of the intensity of excitation. Such stochastic jump can be considered as bifurcation.     

Free in-plane vibration analysis of a curved beam (arch) with arbitrary various concentrated elements

In the existing literature, the information regarding the exact solutions for free in-plane vibrations of the curved beams (or arches) carrying various concentrated elements is rare, particularly for the case with multiple attachments including eccentricities and mass moments of inertias. For this reason, this paper aims at presenting an effective approach to tackle the title problem. First of all, the un-coupled equation of motion for the circumferential displacement of an arch segment is derived. Next, based on the value of the discriminate parameter for a cubic equation, the exact solutions for the three types of roots of the un-coupled equation are determined and, corresponding to each type of roots, all displacement functions for the arch segment in terms of the real numbers (instead of the complex ones) are obtained. Finally, use of the compatible equations for the displacements and slopes together with the equilibrium equations for the forces and moments at each intermediate node and two ends of the entire curved beam, a frequency equation of the form ∣H(ω)∣ = 0 is obtained. It is found that the conventional approach by using the condition “∣H(ω^t)∣ ? ε” to search for the approximate value of ω^t is difficult even if the convergence tolerance ε is greater than 10⁺³ (i.e., ε > 10⁺³) instead of less than 10^?3 (i.e., ε < 10^?3), however, the half-interval method is one of the effective tools for solving the problem if all coefficients of the determinant ∣H(ω)∣ are the real numbers. In addition to comparing with the existing literature, most of the numerical results obtained from the presented method are compared with those obtained from the conventional finite element method (FEM) and good agreement is achieved.