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.     

Bifurcations of travelling wave solutions in a new integrable equation with peakon and compactons

Degasperis and Procesi applied the method of asymptotic integrability and obtain Degasperis–Procesi equation. They showed that it has peakon solutions, which has a discontinuous first derivative at the wave peak, but they did not explain the reason that the peakon solution arises. In this paper, we study these non-smooth solutions of the generalized Degasperis–Procesi equation u_t − u_txx + (b + 1)uu_x = bu_xu_xx + uu_xxx, show the reason that the non-smooth travelling wave arise and investigate global dynamical behavior and obtain the parameter condition under which peakon, compacton and another travelling wave solutions engender. Under some parameter condition, this equation has infinitely many compacton solutions. Finally, we give some explicit expression of peakon and compacton solutions.     

Bright solitons of the variants of the Novikov–Veselov equation with constant and variable coefficients

In this paper, the existence of the bright soliton solution of four variants of the Novikov–Veselov equation with constant and time varying coefficients will be studied. We analyze the solitary wave solutions of the Novikov–Veselov equation in the cases of constant coefficients, time-dependent coefficients and damping term, generalized form, and in 1 + N dimensions with variable coefficients and forcing term. We use the solitary wave ansatz method to derive these solutions. The physical parameters in the soliton solutions are obtained as functions of the dependent coefficients. Parametric conditions for the existence of the exact solutions are given. The solitary wave ansatz method presents a wider applicability for handling nonlinear wave equations.     

Existence and asymptotics of traveling waves for nonlocal diffusion systems

In this paper, we deal with the existence and asymptotic behavior of traveling waves for nonlocal diffusion systems with delayed monostable reaction terms. We obtain the existence of traveling wave front by using upper-lower solutions method and Schauder’s fixed point theorem for c > c₁(τ) and using a limiting argument for c = c₁(τ). Moreover, we find a priori asymptotic behavior of traveling waves with the help of Ikehara’s Theorem by constructing a Laplace transform representation of a solution. Especially, the delay can slow the minimal wave speed for ?₂f(0, 0) > 0 and the delay is independent of the minimal wave speed for ?₂f(0, 0) = 0.     

Numerical modelling of continuous concrete box girder bridges considering construction stages

The aim of this study is to analyse the concrete continuous box girder bridges by considering segmentally construction stages through balanced cantilever method. Time-dependent material properties of concrete and steel are also taken into account. Budan Bridge is selected as a numerical example. The Bridge constructed with balanced cantilever method and located on Artvin–Erzurum highway, Turkey, at 55 + 729.00–56 + 079.00 km is modelled using SAP2000 program. Geometric nonlinearities are taken into consideration in the analysis using P-Delta and large displacement criterion. Time-dependent material properties are considered as compressive strength, aging, shrinkage and creep for concrete, and relaxation for steel. The structural behaviour of the bridge at different construction stages is examined. Variation of internal forces such as bending moment, shear forces and axial forces, and displacements for bridge deck and pier are given with detail.Analyses show that, to obtain real behaviour of concrete bridges, segmentally construction stage analysis using time dependent material properties and geometric nonlinearity should be considered, because construction period continue along time and loads may change during this period and after.     

Soliton-like and periodic form solutions to (2 + 1)-dimensional Toda equation

In this paper, we present a further extended tanh method for constructing exact solutions to nonlinear difference-differential equation(s) (NDDEs) and Lattice equations. By using this method via symbolic computation system MAPLE, we obtain abundant soliton-like and period-form solutions to the (2 + 1)-dimensional Toda equation. Solitary wave solutions are merely a special case in one family. This method can also be used to other nonlinear difference differential equations.     

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.     

Second-order approximation of angle-ply composite laminated thin plate under combined excitations

The influence of the quadratic and cubic terms on non-linear dynamic characteristics of the angle-ply composite laminated rectangular plate with parametric and external excitations is investigated. The method of multiple time scale perturbation is applied to solve the non-linear differential equations describing the system up to and including the second-order approximation. All possible resonance cases are extracted and investigated at this approximation order. Two cases of the sub-harmonic resonances cases (Ω₂ ? 2ω₁ and Ω₂ ? 2ω₂) in the presence of 1:2 internal resonance ω₂ ? 2ω₁ are considered. The stability of the system is investigated using both frequency response equations and phase-plane method. It is quite clear that some of the simultaneous resonance cases are undesirable in the design of such system as they represent some of the worst behavior of the system. Such cases should be avoided as working conditions for the system. Some recommendations regarding the different parameters of the system are reported. Comparison with the available published work is reported.     

Explicit series solution of travelling waves with a front of Fisher equation

In this paper, an analytic technique, namely the homotopy analysis method, is employed to solve the Fisher equation, which describes a family of travelling waves with a front. The explicit series solution for all possible wave speeds 0 < c < +∞ is given. Such kind of explicit series solution has never been reported, to the best of author’s knowledge. Our series solution indicates that the solution contains an oscillation part when 0 < c < 2. The proposed analytic approach is general, and can be applied to solve other similar nonlinear travelling wave problems.     

On “P” property and the column-W property

P-matrices play an important role in the well-posedness of a linear complementarity problem (LCP). Similarly, the well-posedness of a horizontal linear complementarity problem (HLCP) is closely related to the column-W property of a matrix k-tuple.In this paper we first consider the problem of generating P-matrices from a given pair of matrices. Given a matrix pair (D, F) where D is a square matrix of order m and matrix F has m rows, “what are the conditions under which there exists a matrix G such that (D + FG) is a P-matrix?”. We obtain necessary and sufficient conditions for the special case when the column rank of F is m ? 1. A decision algorithm of complexity O(m²) to check whether the given pair of matrices (D, F) is P-matrisable is obtained. We also obtain a necessary and an independent sufficient condition for the general case when rank(F) is less than m ? 1.We then generalise the P-matrix generating problem to the generation of matrix k-tuples satisfying the column-W property from a given matrix (k + 1)-tuple. That is, given a matrix (k + 1)-tuple (D₁, … ,D_k, F), where D_js are square matrices of order m and F is a matrix having m rows, we determine the conditions under which the matrix k-tuple (D₁ + FG₁, … ,D_k + FG_k) satisfies the column-W property. As in the case of P-matrices we obtain necessary and sufficient conditions for the case when rank(F) = m ? 1. Using these conditions a decision algorithm of complexity O(km²) to check whether the given matrix (k + 1)-tuple is column-W matrisable is obtained. Then for the case when rank(F) is less than m ? 1, we obtain a necessary and an independent sufficient condition.For a special sub-class of P-matrices we give a polynomial time decision algorithm for P-matrisability. Finally, we obtain a geometric characterisation of column-W property by generalising the well known separation theorem for P-matrices.     

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.     

Wave propagation properties in oscillatory chains with cubic nonlinearities via nonlinear map approach

Free wave propagation properties in one-dimensional chains of nonlinear oscillators are investigated by means of nonlinear maps. In this realm, the governing difference equations are regarded as symplectic nonlinear transformations relating the amplitudes in adjacent chain sites (n, n + 1) thereby considering a dynamical system where the location index n plays the role of the discrete time. Thus, wave propagation becomes synonymous of stability: finding regions of propagating wave solutions is equivalent to finding regions of linearly stable map solutions. Mechanical models of chains of linearly coupled nonlinear oscillators are investigated. Pass- and stop-band regions of the mono-coupled periodic system are analytically determined for period-q orbits as they are governed by the eigenvalues of the linearized 2D map arising from linear stability analysis of periodic orbits. Then, equivalent chains of nonlinear oscillators in complex domain are tackled. Also in this case, where a 4D real map governs the wave transmission, the nonlinear pass- and stop-bands for periodic orbits are analytically determined by extending the 2D map analysis. The analytical findings concerning the propagation properties are then compared with numerical results obtained through nonlinear map iteration.     

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.     

The effect of a small background inhomogeneity on the asymptotic properties of linear perturbations

General regularities in the evolution of one-dimensional unstable linear perturbations on a weakly inhomogeneous background are studied when, at the initial instant, the perturbations are concentrated in the δ-neighbourhood of a certain point. Times are considered when these perturbations do not fall outside the limits of a certain domain of size l such that δ ? l ? L, where L is the large characteristic size of the background inhomogeneity. With contain assumptions, the effect of the background inhomogeneity on the asymptotic behaviour of the perturbations at long times is taken into account in a general form. The first corrections to the well known asymptotic relation for the evolution of perturbations on a homogeneous background, that arise because of background inhomogeneity, are obtained using Hamilton's method. An example of the use of the proposed approximate method is considered and the error in the approximation is estimated.     

Study of wave–wind interaction at a seawall using a numerical wave channel

This paper presents the study on wind and waves interactions at a seawall using a numerical wave channel. The numerical experiments were conducted for wave overtopping of a 1/4 sloping seawall using several conditions of incident waves and wind speeds. The numerical results were verified against laboratory data in a case for wave overtopping without wind effects. The interaction of waves and wind was analyzed in term of mean wave quantities, overtopping rate and variation of wind velocity at some selected locations. The results showed that the overtopping rate was strongly affected by wind and the wind field was also significantly modified by waves. There exists an effective range of wind speed in comparison with the local shallow wave speed at the breaking location, which gives significant effects to the wave overtopping rates. The maximum of wind adjustment coefficient f_w for wave overtopping rate was strongly related to the mean overtopping rate in the case for no wind. This study also showed that when the mean overtopping rate was greater than 5 × 10⁻⁴ m³/s/m, the maximum of wind adjustment coefficient f_w approached to a specific value of about 1.25.     

Pattern Frequency Sequences and Internal Zeros

Let q be a pattern and let S_n, q(c) be the number of n-permutations having exactly c copies of q. We investigate when the sequence (S_n, q(c))_c ≥ 0 has internal zeros. If q is a monotone pattern it turns out that, except for q = 12 or 21, the nontrivial sequences (those where n is at least the length of q) always have internal zeros. For the pattern q = 1(l + 1)l…2 there are infinitely many sequences which contain internal zeros and when l = 2 there are also infinitely many which do not. In the latter case, the only possible places for internal zeros are the next-to-last or the second-to-last positions. Note that by symmetry this completely determines the existence of internal zeros for all patterns of length at most 3.     

A semi-Markov process with an inverse Gaussian distribution as sojourn time

Let X_n denote the state of a device after n repairs. We assume that the time between two repairs is the time τ taken by a Wiener process {W(t), t ? 0}, starting from w₀ and with drift μ < 0, to reach c ∈ [0, w₀). After the nth repair, the process takes on either the value X_n?1 + 1 or X_n?1 + 2. The probability that X_n = X_n?1 + j, for j = 1, 2, depends on whether τ ? t₀ (a fixed constant) or τ > t₀. The device is considered to be worn out when X_n ? k, where k ∈ {1, 2, …}. This model is based on the ones proposed by Rishel (1991) [1] and Tseng and Peng (2007) [2]. We obtain an explicit expression for the mean lifetime of the device. Numerical methods are used to illustrate the analytical findings.     

Exact chirped solitary-wave solutions for Ginzburg–Landau equation

In this short letter, by applying specially envelope transform and direct ansatz approach to (1 + 1)D Ginzburg–Landau equation the authors obtain a new type of exact solitary wave solution including chirped bright solitary-wave and chirped dark solitary-wave solutions.     

New families of exact soliton-like solutions and dromion solutions to the (2+1)-dimensional higher-order Broer–Kaup equations

Using homogeneous balance method we obtain Bäcklund transformation (BT) and a linear partial differential equation of higher-order Broer–Kaup equations. As a result, new soliton-like solutions and new dromion solution and other exact solutions of (2 + 1)-dimensional higher-order Broer–Kaup equations are given. By analyzing a soliton-like solution, we get some dromions solutions. This method, which can be generalized to some (2 + 1)-dimensional nonlinear evolution equations, is simple and powerful.     

Simulación numérica del flujo turbulento de aire con gotas dispersas de agua a través de separadores de torres de refrigeración

This work presents a numerical study on the turbulent flow of air with dispersed water droplets in separators of mechanical cooling towers. The averaged Navier-Stokes equations are discretised through a finite volume method, using the Fluent and Phoenics codes, and alternatively employing the turbulence models k ? ?, k ? ω and the Reynolds stress model, with low-Re version and wall enhanced treatment refinements. The results obtained are compared with numerical and experimental results taken from the literature. The degree of accuracy obtained with each of the considered models of turbulence is stated. The influence of considering whether or not the simulation of the turbulent dispersion of droplets is analyzed, as well as the effects of other relevant parameters on the collection efficiency and the coefficient of pressure drop. Focusing on four specific eliminators (‘Belgian wave’, ‘H1-V’, ‘L-shaped’ and ‘Zig-zag’), the following ranges of parameters are outlined: 1 ≤ U_e ≤ 5 m/s for the entrance velocity, 2 ≤ D_p ≤ 50 μm for the droplet diameter, 650 ≤ Re ≤ 8.500 for Reynolds number, and 0.05 ≤ P_i ≤ 5 for the inertial parameter. Results reached alternately with Fluent and Phoenics codes are compared. The best results correspond to the simulations performed with Fluent, using the SST k ? ω turbulence model, with values of the dimensionless scaled distance to wall y⁺ in the range 0.2 to 0.5. Finally, correlations are presented to predict the conditions for maximum collection efficiency (100 %), depending on the geometric parameter of removal efficiency of each of the separators, which is introduced in this work.