About fractional quantization and fractional variational principles

In this paper, a new method of finding the fractional Euler–Lagrange equations within Caputo derivative is proposed by making use of the fractional generalization of the classical Faá di Bruno formula. The fractional Euler–Lagrange and the fractional Hamilton equations are obtained within the 1 + 1 field formalism. One illustrative example is analyzed.     

Numerical methods for a class of jump–diffusion systems with random magnitudes

Stochastic differential delay equations with Poisson driven jumps of random magnitude are popular as models in mathematical finance. In this paper, we shall deal with convergence of the semi-implicit Euler method for nonlinear stochastic differential delay equations with random jump magnitudes and show that the approximate solutions strongly converge to the exact solutions with the order 1 − 1/q (q > 1). This result is more general than what they deal with the jump of deterministic magnitude.     

Bifurcation and chaos in discrete-time predator–prey system

The discrete-time predator–prey system obtained by Euler method is investigated. The conditions of existence for flip bifurcation and Hopf bifurcation are derived by using center manifold theorem and bifurcation theory. And numerical simulation results not only show the consistence with the theoretical analysis but also display the new and interesting dynamical behaviors, including period-3, 5, 6, 7, 8, 9, 10, 12, 18, 20, 22, 30, 39-orbits in different chaotic regions, attracting invariant circle, period-doubling bifurcation from period-10 leading to chaos, inverse period-doubling bifurcation from period-5 leading to chaos, interior crisis and boundary crisis, intermittency mechanic, onset of chaos suddenly and sudden disappearance of the chaotic dynamics, attracting chaotic set, and non-attracting chaotic set. In particular, we observe that when the prey is in chaotic dynamic, the predator can tend to extinction or to a stable equilibrium. The computations of Lyapunov exponents confirm the dynamical behaviors. The analysis and results in this paper are interesting in mathematics and biology.     

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.     

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.     

On mathematical models and numerical simulation of the fluidization of polydisperse suspensions

The well-known Masliyah–Lockett–Bassoon (MLB) model for sedimentation of small particles is extended to fluidization of polydisperse suspensions. For N particle species that differ in size and density, this model leads to a first-order system of N conservation laws, which in general is of mixed (in the case N = 2, hyperbolic–elliptic) type. By a simple algebraic steady-state analysis, we derive necessary compatibility conditions on the size and density parameters that admit the formation of stationary fluidized beds. We then proceed to determine the composition of polydisperse fluidized beds of given compatible species by varying the fluidization velocity and the initial composition of the suspensions, and prove that, within the framework of the MLB model combined with the Richardson–Zaki formula, the constructed bidisperse beds always cause the equations to be hyperbolic. This means that these states are always predicted to be stable. The transient behaviour of the MLB model applied to fluidization is illustrated by three numerical examples, in which the system of conservation laws is solved for N = 2, N = 3 and N = 5, respectively. These examples illustrate the effects of bed expansion and layer inversion caused by successively increasing the applied fluidization velocity and show that the predicted fluidized states are indeed attained.     

RBFs approximation method for time fractional partial differential equations

In this paper, radial basis functions (RBFs) approximation method is implemented for time fractional advection–diffusion equation on a bounded domain. In this method the first order time derivative is replaced by the Caputo fractional derivative of order α ∈ (0, 1], and spatial derivatives are approximated by the derivative of interpolation in the Kansa method. Stability and convergence of the method is discussed. Several numerical examples are include to demonstrate effectiveness and accuracy of the method.     

Construction of analytic solution to chaotic dynamical systems using the Homotopy analysis method

Based on a new kind of analytic method, namely the Homotopy analysis method, an analytic approach to solve non-linear, chaotic system of ordinary differential equations is presented. The method is applied to Lorenz system; this system depends on the three parameters: σ, b and the so-called bifurcation parameter R are real constants. Two cases are considered. The first case is when R = 20.5 which corresponds to the transition region and the second case corresponds to R = 23.5 which corresponds to the chaotic region.The validity of the method is verified by comparing the approximation series solution with the results obtained using the standard numerical techniques such as Runge-Kutta method.     

On the stability of fractional neutron point kinetics (FNPK)

The aim of this work is investigate the stability of fractional neutron point kinetics (FNPK). The method applied in this work considers the stability of FNPK as a linear fractional differential equation by transforming the s − plane to the W − plane. The FNPK equations is an approximation of the dynamics of the reactor that includes three new terms related to fractional derivatives, which are explored in this work with an aim to understand their effect in the system stability. Theoretical study of reactor dynamical systems plays a significant role in understanding the behavior of neutron density, which is important in the analysis of reactor safety. The fractional relaxation time (τ^α) for values of fractional-order derivative (α) were analyzed, and the minimum absolute phase was obtained in order to establish the stability of the system. The results show that nuclear reactor stability with FNPK is a function of the fractional relaxation time.     

Filter matrix based on interpolation wavelets for solving Fredholm integral equations

The interpolation wavelet is used to solve the Fredholm integral equation of the second kind in this study. Hence, by the extension of interpolation wavelets that [−1, 1] is divided to 2^N+1 (N ⩾  − 1) subinterval, we have polynomials with a degree less than M + 1 in each new interval. Therefore, by considering the two-scale relation the filter coefficients and filter matrix are used as the proof of theorems. The important point is interpolation wavelets lead to more sparse matrix when we try to solve integral equation by an approximate kernel decomposed to a lower and upper resolution. Using n-time, where (n ⩾ 2), two-scale relation in this method errors of approximate solution as O((2^−(N+1))ⁿ⁺¹). Also, the filter coefficient simplifies the proof of some theorems and the order of convergence is estimated by numerical errors.     

Simple regression in view of elliptical models

For the simple linear model Y = θ1 + βx + ? where the error vector follows the elliptically contoured distribution, we consider the unrestricted, restricted, preliminary test and shrinkage estimators for the intercept parameter, θ when it is suspected that the slope parameter β may be β_o. The exact bias and MSE expressions are derived and the mean-square relative efficiency is taken to determine the relative dominance properties of the proposed estimators in comparison. In the continuation, the optimal level of significance of the preliminary test estimator is tabulated and some graphical result are also displayed.     

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.     

Entropy estimation of the Hénon attractor

The topological entropy of the Hénon attractor is estimated using a function that describes the stable and unstable manifolds of the Hénon map. This function provides an accurate estimate of the length of curves in the attractor. The estimation method presented here can be applied to cases in which the invariant set is not hyperbolic. From the result of the length calculation, we have estimated the topological entropy h as h ～ 0.49703 for the original parameters a = 1.4 and b = 0.3 adopted by Hénon.     

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.     

Estimating the ultimate bound and positively invariant set for a synchronous motor and its application in chaos synchronization

To estimate the ultimate bound and positively invariant set for a dynamic system is an important but quite challenging task in general. This paper has investigated the ultimate bound and positively invariant set of a permanent magnet synchronous motor system. We combine the Lyapunov stability theory with the comparison principle method. For this system, we derive a three-dimensional ellipsoidal ultimate bound and positively invariant set for all the positive values of its parameters σ, γ. In addition, the two-dimensional bound with respect to x ? y are established. Then, it is the two-dimensional estimation about x ? z. Finally, the result is applied to the study of completely chaos synchronization. Numerical simulations are presented to show the effectiveness of the proposed chaos synchronization scheme. At the same time, one numerical example illustrating a localization of a chaotic attractor is presented as well. Numerical simulation is consistent with the results of theoretical calculation.     

线性中立型随机延迟微分方程Euler方法的均方稳定性

本文讨论Euler方法用于求解线性中立型随机延迟微分方程初值问题时数值解的稳定性,利用了一种不同于以往文献中的证明技巧,给出了Euler方法均方稳定的一个充分条件.文末的数值试验证实了本文所获理论结果的正确性.     

Discrepancy principle for DSM II

Let Ay = f, A is a linear operator in a Hilbert space H, y ⊥ N(A) ≔ {u : Au = 0}, R(A) ≔ {h : h = Au, u ∈ D(A)} is not closed, ∥f_δ − f∥ ⩽ δ. Given f_δ, one wants to construct u_δ such that lim_δ→0∥u_δ − y∥ = 0. Two versions of discrepancy principles for the DSM (dynamical systems method) for finding the stopping time and calculating the stable solution u_δ to the original equation Ay = f are formulated and mathematically justified.     

Conditional contractivity of Runge–Kutta methods for nonlinear differential equations with many variable delays

This paper deals with conditional contractivity properties of Runge–Kutta (RK) methods with variable step-size applied to nonlinear differential equations with many variable delays (MDDEs). The concepts of CRN_m(ω, H)- and BN_f(μ, ?)-stability are introduced. It is shown that the numerical solution produced by a BN_f(μ, ?)-stable Runge–Kutta method with an appropriate interpolation is contractive. In particular, these results are also novel for nonlinear differential equations with many constant delays or single variable delay. To obtain BN_f(μ, ?)-stable methods, (k, l)-algebraically stable Runge–Kutta methods are also investigated.     

Characterizations of PG(n − 1, q)\PG(k − 1, q) by Numerical and Polynomial Invariants

We show that the simple matroid PG(n − 1, q)\PG(k − 1, q), for n ≥ 4 and 1 ≤ k ≤ n − 2, is characterized by a variety of numerical and polynomial invariants. In particular, any matroid that has the same Tutte polynomial as PG(n − 1, q)\PG(k − 1, q) is isomorphic to PG(n − 1, q)\PG(k − 1, q).