Decoupling in non-linear systems with two degrees of freedom

The problem of decoupling a class of non-linear two degrees of freedom systems is studied. The coupled non-linear differential equations of motion of the system are shown to be equivalent to a pair of uncoupled equations. This equivalence is established through transformation techniques involving the transformation of both the dependent and independent variables. The sufficient conditions on the form of the non-linearity, for the case wherein the transformed equations are linear, are presented. Several particular cases of interest are also illustrated.

A multistep transversal linearization (MTL) method in non-linear dynamics through a Magnus characterization

A new form of a multi-step transversal linearization (MTL) method is developed and numerically explored in this study for a numeric-analytical integration of non-linear dynamical systems under deterministic excitations. As with other transversal linearization methods, the present version also requires that the linearized solution manifold transversally intersects the non-linear solution manifold at a chosen set of points or cross-section in the state space. However, a major point of departure of the present method is that it has the flexibility of treating non-linear damping and stiffness terms of the original system as damping and stiffness terms in the transversally linearized system, even though these linearized terms become explicit functions of time. From this perspective, the present development is closely related to the popular practice of tangent-space linearization adopted in finite element (FE) based solutions of non-linear problems in structural dynamics. The only difference is that the MTL method would require construction of transversal system matrices in lieu of the tangent system matrices needed within an FE framework. The resulting time-varying linearized system matrix is then treated as a Lie element using Magnus’ characterization [W. Magnus, On the exponential solution of differential equations for a linear operator, Commun. Pure Appl. Math., VII (1954) 649–673] and the associated fundamental solution matrix (FSM) is obtained through repeated Lie-bracket operations (or nested commutators). An advantage of this approach is that the underlying exponential transformation could preserve certain intrinsic structural properties of the solution of the non-linear problem. Yet another advantage of the transversal linearization lies in the non-unique representation of the linearized vector field – an aspect that has been specifically exploited in this study to enhance the spectral stability of the proposed family of methods and thus contain the temporal propagation of local errors. A simple analysis of the formal orders of accuracy is provided within a finite dimensional framework. Only a limited numerical exploration of the method is presently provided for a couple of popularly known non-linear oscillators, viz. a hardening Duffing oscillator, which has a non-linear stiffness term, and the van der Pol oscillator, which is self-excited and has a non-linear damping term.

Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

Unsteady MHD mixed convection flow over an impulsively stretched permeable vertical surface in a quiescent fluid

The unsteady mixed convection flow of an incompressible laminar electrically conducting fluid over an impulsively stretched permeable vertical surface in an unbounded quiescent fluid in the presence of a transverse magnetic field has been investigated. At the same time, the surface temperature is suddenly increased from the surrounding fluid temperature or a constant heat flux is suddenly imposed on the surface. The problem is formulated in such a way that for small time it is governed by Rayleigh type of equation and for large time by Crane type of equation. The non-linear coupled parabolic partial differential equations governing the unsteady mixed convection flow under boundary layer approximations have been solved analytically by using the homotopy analysis method as well as numerically by an implicit finite difference scheme. The local skin friction coefficient and the local Nusselt number are found to decrease rapidly with time in a small time interval and they tend to steady-state values for t* >= 5. They also increase with the buoyancy force and suction, but decrease with injection rate. The local skin friction coefficient increases with the magnetic field, but the local Nusselt number decreases. There is a smooth transition from the unsteady state to the steady state. (C) 2010 Elsevier Ltd. All rights reserved.

A new approach to the study of non-linear non-autonomous systems

In this paper, a new approach to the study of non-linear, non-autonomous systems is presented. The method outlined is based on the idea of solving the governing differential equations of order n by a process of successive reduction of their order. This is achieved by the use of “differential transformation functions”. The value of the technique presented in the study of problems arising in the field of non-linear mechanics and the like, is illustrated by means of suitable examples drawn from different fields such as vibrations, rigid body dynamics, etc.

Application of ultraspherical polynomials to non-linear autonomous systems

This paper deals with the approximate solutions of non-linear autonomous systems by the application of ultraspherical polynomials. From the differential equations for amplitude and phase, set up by the method of variation of parameters, the approximate solutions are obtained by a generalized averaging technique based on the ultraspherical polynomial expansions. The method is illustrated with examples and the results are compared with the digital and analog computer solutions. There is a close agreement between the analytical and exact results.

An approximate analysis of non-linear non-conservative systems subjected to step function excitation

This paper deals with the approximate analysis of the step response of non-linear nonconservative systems by the application of ultraspherical polynomials. From the differential equations for amplitude and phase, set up by the method of variation of parameters, the approximate solutions are obtained by a generalized averaging technique based on ultraspherical polynomial expansions. The Krylov-Bogoliubov results are given by a particular set of these polynomials. The method has been applied to study the step response of a cubic spring mass system in presence of viscous, material, quadratic, and mixed types of damping. The approximate results are compared with the digital and analogue computer solutions and a close agreement has been found between the analytical and the exact results.

Natural Convection on a Horizontal Cone in a Porous Medium with Non-Uniform Wall Temperature/Concentration or Heat/Mass Flux and Suction/Injection

The steady natural convection flow on a horizontal cone embedded in a saturated porous medium with non-uniform wall temperature/concentration or heat/mass flux and suction/injection has been investigated. Non-similar solutions have been obtained. The nonlinear couple differential equations under boundary layer approximations governing the flow have been numerically solved. The Nusselt and Sherwood numbers are found to depend on the buoyancy forces, suction/injection rates, variation of wall temperature/concentration or heat/mass flux, Lewis number and the non-Darcy parameter.

Nonsimilar boundary layers for non-Darcy mixed convection flow about a horizontal surface in a saturated porous medium

The nonsimilar non-Darcy mixed convection flow about a heated horizontal surface in a saturated porous medium has been studied when the surface temperature is a power function of distance (Tw = T∞ ± Axλ). The analysis is performed for the cases of parallel and stagnation flows with favourable induced pressure gradient. The partial differential equations governing the flow have been solved numerically using the Keller box method. The heat transfer is enhanced due to the buoyancy parameter and wall temperature, but the non-Darcy parameter reduces it. For non-Darcy flow, the similarity solution exists only for the case of parallel flow.

Non-Darcy mixed convection flow with thermal dispersion on a vertical cylinder in a saturated porous medium

The non-Darcy mixed convection flow on a vertical cylinder embedded in a saturated porous medium has been studied taking into account the effect of thermal dispersion. Both forced flow and buoyancy force dominated cases with constant wall temperature condition have been considered. The governing partial differential equations have been solved numerically using the Keller box method. The results are presented for the buoyancy parameter which cover the entire regime of mixed convection flow ranging from pure forced convection to pure free convection. The effect of thermal dispersion is found to be more pronounced on the heat transfer than on the skin friction and it enhances the heat transfer but reduces the skin friction.

Nonsimilar mixed convection flow of a non-Newtonian fluid past a vertical wedge

The flow and heat transfer characteristics of a second-order fluid over a vertical wedge with buoyancy forces have been analysed. The coupled nonlinear partial differential equations governing the nonsimilar mixed convection flow have been solved numerically using Keller box method. The effects of the buoyancy parameter, viscoelastic parameter, mass transfer parameter, pressure gradient parameter, Prandtl number and viscous dissipation parameter on the skin friction and heat transfer have been examined in detail. Particular cases of the present results match exactly with those available in the literature.

Non-rotating beams isospectral to a given rotating uniform beam

In this paper, we seek to find non-rotating beams with continuous mass and flexural stiffness distributions, that are isospectral to a given uniform rotating beam. The Barcilon-Gottlieb transformation is used to convert the fourth order governing equation of a non-rotating beam, to a canonical fourth order eigenvalue problem. If the coefficients in this canonical equation match with the coefficients of the uniform rotating beam equation, then the non-rotating beam is isospectral to the given rotating beam. The conditions on matching the coefficients leads to a pair of coupled differential equations. We solve these coupled differential equations for a particular case, and thereby obtain a class of non-rotating beams that are isospectral to a uniform rotating beam. However, to obtain isospectral beams, the transformation must leave the boundary conditions invariant. We show that the clamped end boundary condition is always invariant, and for the free end boundary condition to be invariant, we impose certain conditions on the beam characteristics. We also verify numerically that the frequencies of the non-rotating beam obtained using the finite element method (FEM) are the exact frequencies of the uniform rotating beam. Finally, the example of beams having a rectangular cross-section is presented to show the application of our analysis. Since experimental determination of rotating beam frequencies is a difficult task, experiments can be easily conducted on these rectangular non-rotating beams, to calculate the frequencies of the rotating beam. (c) 2012 Elsevier Ltd. All rights reserved.

Analytical test functions for free vibration analysis of rotating non-homogeneous Timoshenko beams

In this paper, the governing equations for free vibration of a non-homogeneous rotating Timoshenko beam, having uniform cross-section, is studied using an inverse problem approach, for both cantilever and pinned-free boundary conditions. The bending displacement and the rotation due to bending are assumed to be simple polynomials which satisfy all four boundary conditions. It is found that for certain polynomial variations of the material mass density, elastic modulus and shear modulus, along the length of the beam, the assumed polynomials serve as simple closed form solutions to the coupled second order governing differential equations with variable coefficients. It is found that there are an infinite number of analytical polynomial functions possible for material mass density, shear modulus and elastic modulus distributions, which share the same frequency and mode shape for a particular mode. The derived results are intended to serve as benchmark solutions for testing approximate or numerical methods used for the vibration analysis of rotating non-homogeneous Timoshenko beams.

Unsteady free convection flow in the stagnation-point region of a rotating sphere

The unsteady free convection boundary-layer flow in the forward stagnation-point region of a sphere, which is rotating with time-dependent angular velocity in an ambient fluid, has been studied. Both constant wall temperature and constant hear flux conditions have been considered. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The skin friction and the heat transfer are enhanced by the buoyancy force. The effect of the buoyancy force is found to be more pronounced for smaller Prandtl numbers than for larger Prandtl numbers. For a given buoyancy force, the heat transfer increases with an increase in Prandtl number, but the skin friction decreases.

Unsteady free convection flow in the stagnation-point region of a three-dimensional body

The unsteady free convection flow in the stagnation-point region of a heated three-dimensional body placed in an ambient fluid is studied under boundary layer approximations. We have considered the case where there is an initial steady state that is perturbed by a step-change in the wall temperature. The non-linear coupled partial differential equations governing the free convection flow are solved numerically using a finite difference scheme. The presented results show the temporal development of the momentum and thermal boundary layer characteristics.