The residue harmonic balance for fractional order van der Pol like oscillators

When seeking a solution in series form, the number of terms needed to satisfy some preset requirements is unknown in the beginning. An iterative formulation is proposed so that when an approximation is available, the number of effective terms can be doubled in one iteration by solving a set of linear equations. This is a new extension of the Newton iteration in solving nonlinear algebraic equations to solving nonlinear differential equations by series. When Fourier series is employed, the method is called the residue harmonic balance. In this paper, the fractional order van der Pol oscillator with fractional restoring and damping forces is considered. The residue harmonic balance method is used for generating the higher-order approximations to the angular frequency and the period solutions of above mentioned fractional oscillator. The highly accurate solutions to angular frequency and limit cycle of the fractional order van der Pol equations are obtained analytically. The results that are obtained reveal that the proposed method is very effective for obtaining asymptotic solutions of autonomous nonlinear oscillation systems containing fractional derivatives. The influence of the fractional order on the geometry of the limit cycle is investigated for the first time.     

Approximate analytical solutions to a conservative oscillator using global residue harmonic balance method

In the current research paper, a conservative system comprising of a mass grounded by linear and nonlinear springs in series connection is studied. The equation of motion for the aforementioned system has been derived as a nonlinear ordinary differential equation with inertia and static–type cubic nonlinearities. The global residue harmonic balance method is applied to obtain an approximate analytical frequency and periodic solution of the problem. Using the obtained analytical expressions, the influences of the hardening and softening nonlinear spring on the non–dimensional frequency are investigated. The results show that developing the system nonlinearity leads the displacement of the mass and the deflection of linear spring to approach each other. Moreover, comparison of the results obtained using the proposed procedure with those achieved by other methods such as numerical method, variational iteration method and harmonic balance approach demonstrates the accuracy and advantages of the current approach.     

Harmonic balance approach to the periodic solutions of the (an)harmonic relativistic oscillator

The first-order harmonic balance method via the first Fourier coefficient is used to construct two approximate frequency-amplitude relations for the relativistic oscillator for which the nonlinearity (anharmonicity) is a relativistic effect due to the time line dilation along the world line. Making a change of variable, a new nonlinear differential equation is obtained and two procedures are used to approximately solve this differential equation. In the first the differential equation is rewritten in a form that does not contain a square-root expression, while in the second the differential equation is solved directly. The approximate frequency obtained using the second procedure is more accurate than the frequency obtained with the first due to the fact that, in the second procedure, application of the harmonic balance method produces an infinite set of harmonics, while in the first procedure only two harmonics are produced. Both approximate frequencies are valid for the complete range of oscillation amplitudes, and excellent agreement of the approximate frequencies with the exact one are demonstrated and discussed. The discrepancy between the first-order approximate frequency obtained by means of the second procedure and the exact frequency never exceeds 1.6%. We also obtained the approximate frequency by applying the second-order harmonic balance method and in this case the relative error is as low 0.31% for all the range of values of amplitude of oscillation A.     

A new formulation for the existence and calculation of nonlinear normal modes

A new formulation is presented here for the existence and calculation of nonlinear normal modes in undamped nonlinear autonomous mechanical systems. As in the linear case an expression is developed for the mode in terms of the amplitude, mode shape and frequency, with the distinctive feature that the last two quantities are amplitude and total phase dependent. The dynamic of the periodic response is defined by a one-dimensional nonlinear differential equation governing the total phase motion. The period of the oscillations, depending only on the amplitude, is easily deduced. It is established that the frequency and the mode shape provide the solution to a 2π-periodic nonlinear eigenvalue problem, from which a numerical Galerkin procedure is developed for approximating the nonlinear modes. The procedure is applied to various mechanical systems with two degrees of freedom.     

Analytic solution of transversal oscillation of quintic non-linear beam with energy balance method and global residue harmonic balance method

Beams are very important element in structures because of its widespread usage in steel construction such as buildings, bridges, aircraft arm, aerospace vehicles and many other industrial usages. To attain a proper design of the beam structures, it is essential to realize how the beam vibrates in its transverse mode which in turn yields the natural frequency of the system. The main objective of present study is to obtain highly accurate analytical solutions for vibrations of a beam with quintic nonlinearity. Two new analytical methods are the improved energy balance method, and the global residue harmonic balance method. Subsequently, the analytical results are compared with the numerical solution of the exact equation in order to evaluate the correctness of the applied approaches.     

Stochastic Quantization of Time-Dependent Systems by the Haba and Kleinert Method

The stochastic quantization method recently developed by Haba and Kleinert is extended to non-autonomous mechanical systems, in the case of the time-dependent harmonic oscillator. In comparison with the autonomous case, the quantization procedure involves the solution of a nonlinear, auxiliary equation. Using a rescaling transformation, the Schrö-din-ger equation for the time-dependent harmonic oscillator is obtained after averaging of a classical stochastic differential equation.     

Stochastic Quantization of Time-Dependent Systems by the Haba and Kleinert Method

The stochastic quantization method recently developed by Haba and Kleinert is extended to non-autonomous mechanical systems, in the case of the time-dependent harmonic oscillator. In comparison with the autonomous case, the quantization procedure involves the solution of a nonlinear, auxiliary equation. Using a rescaling transformation, the Schrödinger equation for the time-dependent harmonic oscillator is obtained after averaging of a classical stochastic differential equation.     

Path-Integral Approach to the Statistical Physics of One-Dimensional Random Systems

A path-integral method is extended and developed to investigate the statistical physics of one-dimensional random systems. Evaluation of the one-particle partition function and density matrix is simplified to finding a solution for a second-order ordinary differential equation. This makes it possible to obtain analytic solutions or conduct accurate numerical calculations for the random systems. With this approach, an analytical solution for the Gaussian model is obtained and the statistical physics of the Frisch–Lloyd model is studied.     

非线性声波方程二次谐波高阶近似解及实验验证

非线性声波方程的传统二次谐波二阶近似解不够精确,对测量材料的非线性系数造成较大误差.利用摄动法展开非线性声波方程得到一系列非齐次偏微分方程,根据低阶解的性质拟定高阶特解的形式,通过符号计算求出高次谐波特解,对所有二次谐波成分求和最终获得二次谐波的高阶近似解.在水中开展非线性声学实验验证高阶近似解,结果表明:二次谐波相对...     

Numerical investigation of the translational motion of bubbles: The comparison of capabilities of the time-resolved and the time-averaged methods

In the present study, the accuracies of two different numerical approaches used to model the translational motion of acoustic cavitational bubble in a standing acoustic field are compared. The less accurate but less computational demanding approach is to decouple the equation of translational motion from the radial oscillation, and solve it by calculating the time-averaged forces exerted on the bubble for one acoustic cycle. The second approach is to solve the coupled ordinary differential equations directly, which provides more accurate results with higher computational effort. The investigations are carried out in the parameter space of the driving frequency, pressure amplitude and equilibrium radius. Results showed that both models are capable to reveal stable equilibrium positions; however, in the case of oscillatory solutions, the difference in terms of translational frequency may be more than three fold, and the amplitude of translational motion obtained by the time-averaged method is roughly 1.5 times higher compared to the time-resolved solution at particular sets of parameters. This observation implies that where the transient behaviour is important, the time-resolved approach is the proper choice for reliable results.     

Boundary conditions for differential approximations

Low order spherical harmonic (P-N) approximations are applied to a radiative transfer Marshak wave problem. A modified Milne boundary condition is developed for the P-2 approximation, similar to one suggested earlier for the P-1 approximation. Comparison with exact Monte Carlo results suggests that this modified P-2 method may be an accurate and generally applicable differential approximation to the equation of transfer. The Monte Carlo results presented should be useful for testing other approximate formulations of radiative transfer and validating time dependent numerical solution methods for the equation of transfer.     

Rational Chebyshev pseudospectral approach for solving Thomas-Fermi equation

In this Letter we propose a pseudospectral method for solving Thomas-Fermi equation which is a nonlinear ordinary differential equation on semi-infinite interval. This approach is based on rational Chebyshev pseudospectral method. This method reduces the solution of this problem to the solution of a system of algebraic equations. Comparison with some numerical solutions shows that the present solution is highly accurate.     

Analytical periodic motions in a parametrically excited, nonlinear rotating blade

The stability and bifurcation analyses of periodic motions in a rotating blade subject to a torsional excitation are investigated. For high speed rotations, cubic geometric nonlinearity and gyroscopic effects of the rotating blade are considered. From the Galerkin method, the partial differential equation of the nonlinear rotating blade is simplified to the ordinary differential equations, and periodic motions and stability of the rotating blade are studied by the generalized harmonic balance method. The analytical and numerical results of periodic solutions are compared. The rich dynamics and co-existing periodic solutions of the nonlinear rotating blades are investigated.     

Spline-Galerkin Solution of Dynamic Equations for Particle Comminution and Collection

The population balance equation is solved for particles undergoing a combination of growth, comminution, and collection. The approximation method is to use a weighted Galerkin technique with cubic B-splines and an implicit scheme for solving the system of ordinary differential equations. The cubic splines are defined on a graded mesh. The performance of the method is investigated by solving a model problem with simple but nonsmooth kernels. The weight function is chosen so that singularities in the equation can be easily treated. A self-similar solution for comminuted particles is shown to be a useful representation for the solution of the population balance equation provided that this equation is solved over a sufficiently long time interval. Stationary solutions of the equation are obtained for a model that describes both particle comminution and collection.     

LIMIT CYCLE FLUTTER OF AIRFOILS IN STEADY AND UNSTEADY FLOWS

The limit cycle flutter of a two-dimensional wing with non-linear pitching stiffness is investigated. For modelling the aerodynamic forces of the wing steady linear and non-linear models as well as an unsteady model were used. The flutter speed was calculated using the harmonic balance method and by predicting Hopf bifurcation. Analytical solutions based on the centre manifold theory and normal forms were obtained as were results given by the harmonic balance method. The analytical solutions were compared with those obtained by numerical integration. The results show that the harmonic balance method can forecast flutter speed with a good accuracy while analytical solutions based on centre manifold theorem are accurate only in a small neighbourhood of the bifurcation point. The oscillation of the airfoil after flutter for two different models, linear and non-linear pitching stiffness were compared with each other and the flutter speeds for two linear steady and an unsteady aerodynamic model calculated. The obtained results show that flutter analysis based on the linear steady model is conservative only for the ratios of plunge frequency to pitch frequency lower than 1.     

Formulation of Zeeman modulation as a signal filter

The Bloch equation containing a Zeeman modulation field is solved analytically by treating the Zeeman modulation frequency as a perturbation. The absorption and dispersion signals at both 0 degrees and 90 degrees modulation phase are obtained. The solutions are valid to first order in the modulation frequency, but are otherwise valid for any value of modulation amplitude or microwave amplitude. A first order treatment of modulation frequency is shown to be a valid approximation over a wide range of typical experimental EPR conditions. The solutions derived from the Bloch equation suggest that the effect of over-modulation on first and second harmonic EPR spectra can be formulated as a mathematical filter that smoothes and broadens the under-modulated signal. The only adjustable filter parameter is a width that is equivalent to the applied peak-to-peak modulation amplitude. The true spin-spin and spin-lattice relaxation rates are completely determined from the under-modulated spectrum. The filters derived from the analytic solutions of the Bloch equation in the linear limit of modulation frequency are tested against numerical solutions of the Bloch equation that are valid for any modulation frequency to show their applicability. The filters are further tested using experimental EPR spectra. Experimental under-modulated spectra are mathematically filtered and compared with the experimental over-modulated spectra. The application of modulation filters to STEPR spectra is explored and limitations are discussed.     

Response control for the externally excited van der Pol oscillator with non-local feedback

A non-local control force is introduced in such a way to obtain a third-order nonlinear differential equation (jerk dynamics) and to control nonlinear vibrations in an externally excited van der Pol oscillator. Two first-order nonlinear ordinary differential equations governing the modulation of the amplitude and the phase of solutions are derived and subsequently the performance of the control strategy is investigated. Excitation amplitude–response and frequency–response curves are shown. In certain cases when the excitation amplitude is very low an approximate analytic solution corresponding to a modulated two-period quasi-periodic motion can be obtained for the uncontrolled system. Uncontrolled and controlled systems are compared and the appropriate choices for the feedback gains are found in order to reduce the amplitude peak of the response and to exclude the possibility of quasi-periodic motion. Numerical simulation confirms the validity of the new method.     

Signs of life

Programmable electronic calculators provide a speedy means of performing scientific calculations, e.g. the evaluation of a given polynomial for many different x values and the performance of iterative procedures for the solution of polynomial equations. Both of these uses appear in variational calculations of the energy of simple quantum mechanical systems such as the perturbed harmonic oscillator. For this system the traditional perturbation method has some drawbacks, end so it is useful to find its energy levels directly by purely numerical methods. The electronic calculator can do this if the relevant Schrödinger differential equation is transformed into a difference equation. The diffusion equation is a partial differential equation, but can be converted to an ordinary differential equation and thence to a pair of difference equations which can be solved on the calculator. This applies even if the diffusion coefficient depends on the concentration, so that the associated ordinary differential equation is non-linear.     

Accurate analytic presentation of solution of the Schrödinger equation with arbitrary physical potential: Excited states

The quasilinearization method (QLM) is used to approximate analytically, both the ground state and the excited state solutions of the Schrödinger equation for arbitrary potentials. The procedure of approximation was demonstrated on examples of a few often used physical potentials such as the quartic anharmonic oscillator, the Yukawa and the spiked harmonic oscillator potentials. The accurate analytic expressions for the ground and excited state energies and wave functions were presented. These high-precision approximate analytic representations are obtained by first casting the Schrödinger equation into a nonlinear Riccati form and then solving that nonlinear equation analytically in the first QLM iteration. In the QLM the nonlinear differential equation is treated by approximating the nonlinear terms by a sequence of linear expressions. The QLM is iterative but not perturbative and gives stable solutions to nonlinear problems without depending on the existence of a smallness parameter. The method provides final and reasonable results for both small and large values of the coupling constant and is able to handle even super singular potentials for which each term of the perturbation theory is infinite and the perturbation expansion does not exist. The choice of zero iteration is based on general features of solutions near the boundaries. In order to estimate the accuracy of the QLM solutions, the exact numerical solutions were found as well. The first QLM iterate given by analytic expression allows to estimate analytically the role of different parameters and the influence of their variation on different characteristics of the relevant quantum systems.     

A search for the simplest chaotic partial differential equation

A numerical search for the simplest chaotic partial differential equation (PDE) suggests that the Kuramoto-Sivashinsky equation is the simplest chaotic PDE with a quadratic or cubic nonlinearity and periodic boundary conditions. We define the simplicity of an equation, enumerate all autonomous equations with a single quadratic or cubic nonlinearity that are simpler than the Kuramoto-Sivashinsky equation, and then test those equations for chaos, but none appear to be chaotic. However, the search finds several chaotic, ill-posed PDEs; the simplest of these, in the discrete approximation of finitely many, coupled ordinary differential equations (ODEs), is a strikingly simple, chaotic, circulant ODE system.