Analysis of a weakly non-linear autonomous oscillator by means of the field method

In this paper the field method is extended to the study of oscillatory systems with two degrees of freedom and weak quadratic non-linearity. The basic field method concept is combined with the technique of multiple time scales and the solution for both non-resonant case and the case out of the first resonance are found. The qualitative analysis of behavior in the resonant area is done by determining the values of the “adelphic” integral.     

Conservation laws by means of a new mixed method

In this paper, by using a mixed approach, recently introduced by the authors, some conservation laws of partial differential equations are derived. The method merges the Ibragimov’s method and the one by Anco and Bluman. In particular, by applying this new mixed method, we determine all zero-th order conservation laws of Chaplygin and Shallow Water equations, as well as new conservation laws for a second order partial differential equation involving an arbitrary function.     

Approximate symmetries and conservation laws of the geodesic equations for the Schwarzschild metric

Approximate symmetries have been defined in the context of differential equations and systems of differential equations. They give approximately, conserved quantities for Lagrangian systems. In this paper, the exact and the approximate symmetries of the system of geodesic equations for the Schwarzschild metric, and in particular for the radial equation of motion, are studied. It is noted that there is an ambiguity in the formulation of approximate symmetries that needs to be clarified by consideration of the Lagrangian for the system of equations. The significance of approximate symmetries in this context is discussed.     

Conservation laws and solutions of a quantum drift-diffusion model for semiconductors

A non-linear system of partial differential equations describing a quantum drift-diffusion model for semiconductor devices is investigated by methods of group analysis. An infinite number of conservation laws associated with symmetries of the model are found. These conservation laws are used for representing the system of equations under consideration in the conservation form. Exact solutions provided by the method of conservation laws are discussed. These solutions are different from invariant solutions.     

Lagrangian dynamics and possible isochronous behavior in several classes of non-linear second order oscillators via the use of Jacobi last multiplier

In this paper, we employ the technique of Jacobi Last Multiplier (JLM) to derive Lagrangians for several important and topical classes of non-linear second-order oscillators, including systems with variable and parametric dissipation, a generalized anharmonic oscillator, and a generalized Lane–Emden equation. For several of these systems, it is very difficult to obtain the Lagrangians directly, i.e., by solving the inverse problem of matching the Euler–Lagrange equations to the actual oscillator equation. In order to facilitate the derivation of exact solutions, and also investigate possible isochronous behavior in the analyzed systems, we next invoke some recent theoretical results and attempt to map the potential term to either the simple harmonic oscillator or the isotonic potential for specific values of the coefficient parameters of each non-linear oscillator. We find non-trivial parameter sets corresponding to isochronous dynamics in some of the considered systems, but none in others. Finally, the Lagrangians obtained here are coupled to Noether׳s theorem, leading to non-trivial conservation laws for several of the oscillators.     

A stable gradient reconstruction for the MUSCL schemes applied to systems of conservation laws

The numerical simulations of compressible flows need more and more accuracy. Several approaches can be used to increase the order of the scheme. One of the most popular has been introduced by [11]: the so-called MUSCL scheme. In the present work, we focus our attention on the robustness of this method. We propose a relevant CFL restriction and a new gradient reconstruction strategy to enforce stability of the method. The novelty stays in the fact that non conservation argument is involved in the gradient reconstruction. Numerical results are performed using this new variant of the MUSCL scheme.     

Free vibration of two coupled nonlinear oscillators

An analytical investigation is carried out on the free vibration of a two degree of freedom weakly nonlinear oscillator. Namely, the method of multiple time scales is first applied in deriving modulation equations for a van der Pol oscillator coupled with a Duffing oscillator. For the case of non-resonant oscillations, these equations are in standard normal form of a codimension two (Hopf-Hopf) bifurcation, which permits a complete analysis to be performed. Three different types of asymptotic states-corresponding to trivial, periodic and quasiperiodic motions of the original system-are obtained and their stability is analyzed. Transitions between these different solutions are also identified and analyzed in terms of two appropriate parameters. Then, effects of a coupling, a detuning, a nonlinear stiffness and a damping parameter are investigated numerically in a systematic manner. The results are interpreted in terms of classical engineering terminology and are related to some relatively new findings in the area of nonlinear dynamical systems.     

On symmetries, conservation laws and invariant solutions of the foam-drainage equation

This study deals with symmetry group properties and conservation laws of the foam-drainage equation. Firstly, we study the classical Lie symmetries, optimal systems, similarity reductions and similarity solutions of the foam-drainage equation which are obtained through the Lie group method of infinitesimal transformations. Secondly, using the new general theorem on non-local conservation laws and partial Lagrangian approach, local and non-local conservation laws are also studied and, finally, non-classical symmetries are derived.     

Large deflections of cantilever beams of non-linear elastic material under a combined loading

Large deflection of cantilever beams made of Ludwick type material subjected to a combined loading consisting of a uniformly distributed load and one vertical concentrated load at the free end was investigated. Governing equation was derived by using the shearing force formulation instead of the bending moment formulation because in the case of large deflected member, the shearing force formulation possesses some computational advantages over the bending moment formulation. Since the problem involves both geometrical and material non-linearities, the governing equation is complicated non-linear differential equation, which would in general require numerical solutions to determine the large deflection for a given loading. Numerical solution was obtained by using Butcher's fifth order Runge-Kutta method and are presented in a tabulated form.     

A comparison of numerical methods applied to non-linear adsorption columns

Eight numerical schemes (first-order upstream finite difference, MacCormack, explicit Taylor–Galerkin, random choice, flux-corrected transport, ENO, TVD, and Euler–Lagrange methods) are compared on the basis of their computational efficiency for one-dimensional non-linear convection–diffusion problems. For the ideal chromatographic equation for which an exact solution exists, errors plotted against computational times show that the best methods are the random choice, Euler–Lagrange and flux-corrected MacCormack methods. Even when significant diffusion is added to the model, steep gradients are possible because of non-linearities. In such an instance, the random choice and flux-corrected transport methods give the best performance. One can now tackle more complicated problems and refer to this comparative study in order to choose an adequate numerical method which will provide sufficiently accurate results at a reasonable cost.     

Dynamics of two delay coupled van der Pol oscillators

In this paper, the dynamics of a system of two van der Pol oscillators with delayed position and velocity coupling is studied by the method of averaging together with truncation of Taylor expansions. According to the slow-flow equations, the dynamics of 1:1 internal resonance is more complex than that of non-1:1 internal resonance. For 1:1 internal resonance, the stability and the number of periodic solutions vary with different time delay for given coupling coefficients. The condition necessary for saddle-node and Hopf bifurcations for symmetric modes, namely in-phase and out-of-phase modes, are determined. The numerical results, obtained from direct integration of the original equation, are found to be in good agreement with analytical predictions.     

Cohesive zone laws for fatigue crack growth: Numerical field projection of the micromechanical damage process in an elasto-plastic medium

Cohesive zone failure models are widely used to simulate fatigue crack propagation under cyclic loading, but the model parameters are phenomenological and are not closely tied to the underlying micromechanics of the problem. In this paper, we will inversely extract the cohesive zone laws for fatigue crack growth in an elasto-plastic ductile solid using a field projection method (FPM), which projects the equivalent tractions and separations at the cohesive crack-tip from field information outside the process zone. In our small-scale yielding model, a single row of discrete voids is deployed directly ahead of a crack in an elasto-plastic medium subjected to cyclic mode I K-field loading. Damage accumulation under cyclic loading is captured by the growth of voids within the micro-voiding zone ahead of the crack, while the evolution of the cohesive zone law representing the micro-voiding zone is inversely extracted via the FPM. We show that the field-projected cohesive zone law captures the essential micromechanisms of fatigue crack growth in the ductile medium: from loading and unloading hysteresis caused by void growth and plastic hardening, to the softening damage locus associated with crack propagation via a void by void growth mechanism. The results demonstrate the effectiveness of the FPM in obtaining a micromechanics-based cohesive zone law in-place of phenomenological models, which opens the way for a unified treatment of fatigue crack problems.     

Response and stability of strongly non-linear oscillators under wide-band random excitation

A new stochastic averaging procedure for single-degree-of-freedom strongly non-linear oscillators with lightly linear and (or) non-linear dampings subject to weakly external and (or) parametric excitations of wide-band random processes is developed by using the so-called generalized harmonic functions. The procedure is applied to predict the response of Duffing–van der Pol oscillator under both external and parametric excitations of wide-band stationary random processes. The analytical stationary probability density is verified by digital simulation and the factors affecting the accuracy of the procedure are analyzed. The proposed procedure is also applied to study the asymptotic stability in probability and stochastic Hopf bifurcation of Duffing–van der Pol oscillator under parametric excitations of wide-band stationary random processes in both stiffness and damping terms. The stability conditions and bifurcation parameter are simply determined by examining the asymptotic behaviors of averaged square-root of total energy and averaged total energy, respectively, at its boundaries. It is shown that the stability analysis using linearized equation is correct only if the linear stiffness term does not vanish.     

Scaling behavior of coupled conservative weakly nonlinear oscillators

The scaling of the solution of coupled conservative weakly nonlinear oscillators is demonstrated and analyzed through evaluating the normal modes and their bifurcation with an equivalent linearization technique and calculating the general solutions with a method of multiple seales. The scaling law for coupled Duffing oscillators is that the coupling intensity should be proportional to the total energy of the system.Present address: Department of Chemistry and Chemistry and Chemical Engineering, Stevens Institute of Technology, Hoboken, New Jersey 07030, U.S.A.     

Dynamical analysis of two coupled parametrically excited van der Pol oscillators

The dynamical behavior of two coupled parametrically excited van der pol oscillators is investigated in this paper. Based on the averaged equations, the transition boundaries are sought to divide the parameter space into a set of regions, which correspond to different types of solutions. Two types of periodic solutions may bifurcate from the initial equilibrium. The periodic solutions may lose their stabilities via a generalized static bifurcation, which leads to stable quasi-periodic solutions, or via a generalized Hopf bifurcation, which leads to stable 3D tori. The instabilities of both the quasi-periodic solutions and the 3D tori may directly lead to chaos with the variation of the parameters. Two symmetric chaotic attractors are observed and for certain values of the parameters, the two attractors may interact with each other to form another enlarged chaotic attractor.     

Hybrid dynamics of two coupled oscillators that can impact a fixed stop

We consider two linearly coupled masses, where one mass can have inelastic impacts with a fixed, rigid stop. This leads to the study of a two degree of freedom, piecewise linear, frictionless, unforced, constrained mechanical system. The system is governed by three types of dynamics: coupled harmonic oscillation, simple harmonic motion and discrete rebounds. Energy is dissipated discontinuously in discrete amounts, through impacts with the stop. We prove the existence of a non-zero measure set of orbits that lead to infinite impacts with the stop in a finite time. We show how to modify the mathematical model so that forward existence and uniqueness of solutions for all time is guaranteed. Existence of hybrid periodic orbits is shown. A geometrical interpretation of the dynamics based on action coordinates is used to visualize numerical simulation results for the asymptotic dynamics.     

Slow Eigenvalues of Self-similar Solutions of the Dafermos Regularization of a System of Conservation Laws: an Analytic Approach

The Dafermos regularization of a system of n hyperbolic conservation laws in one space dimension has, near a Riemann solution consisting of n Lax shock waves, a self-similar solution u = u _ε(X/T). In Lin and Schecter (2003, SIAM J. Math. Anal. 35, 884–921) it is shown that the linearized Dafermos operator at such a solution may have two kinds of eigenvalues: fast eigenvalues of order 1/ε and slow eigenvalues of order one. The fast eigenvalues represent motion in an initial time layer, where near the shock waves solutions quickly converge to traveling-wave-like motion. The slow eigenvalues represent motion after the initial time layer, where motion between the shock waves is dominant. In this paper we use tools from dynamical systems and singular perturbation theory to study the slow eigenvalues. We show how to construct asymptotic expansions of eigenvalue-eigenfunction pairs to any order in ε. We also prove the existence of true eigenvalue-eigenfunction pairs near the asymptotic expansions.     

Application of the Krylov–Bogoliubov–Mitropolsky method to weakly damped strongly non-linear planar Hamiltonian systems

In this paper, an analytical approximation of damped oscillations of some strongly non-linear, planar Hamiltonian systems is considered. To apply the Krylov–Bogoliubov–Mitropolsky method in this strongly non-linear case, we mainly provide the formal and exact solutions of the homogeneous part of the variational equations with periodic coefficients resulting from the Hamiltonian systems. It is shown that these are simply expressed in terms of the partial derivatives of the solutions, written in action-angle variables, of the Hamiltonian systems. Two examples, including a non-linear harmonic oscillator and the Morse oscillator, are presented to illustrate this extension of the method. The approximate first order solution obtained in each case is observed to be quite satisfactory.     

Primary resonance of multiple degree-of-freedom dynamic systems with strong non-linearity using the homotopy analysis method

A homotopy analysis method（HAM）is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existence of small parameters in the considered equation.The HAM provides a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter.Two examples are presented to show that the HAM solutions agree well with the results of the modified Linstedt-Poincar＇e method and the incremental harmonic balance method.