Symmetry and ordinary differential constraints

The representative generalized symmetries of any ordinary differential equation are described in terms of its invariants. This identifies the evolution equations compatible with a given constraint. The restriction of the flow of a compatible equation to the solution space of the constraint is generated by the corresponding internal symmetry. This reduces the evolution equation to a finite dimensional system of first-order ordinary differential equations. The Euler–Lagrange equation of any conserved density of a given evolution equation yields such a reduction. Other examples include the generalized method of separation of variables, the characterization of separable evolution equations, and the characterization of equations with complete families of wave solutions. A Newton equation is compatible with an ordinary differential constraint if and only if the constraint is affine, with force field symmetry, in which case the equation reduces to a finite-dimensional dynamical system. Newton equations with complete families of characteristic solutions reduce to central force problems on solution spaces of linear constraints.     

Lie group analysis for the effect of temperature-dependent fluid viscosity and thermophoresis particle deposition on free convective heat and mass transfer under variable stream conditions

This paper examines a steady two-dimensional flow of incompressible fluid over a vertical stretching sheet. The fluid viscosity is assumed to vary as a linear function of temperature. A scaling group of transformations is applied to the governing equa- tions. The system remains invariant due to some relations among the transformation parameters. After finding three absolute invariants, a third-order ordinary differential equation corresponding to the momentum equation and two second-order ordinary differential equations corresponding to energy and diffusion equations are derived. The equations along with the boundary conditions are solved numerically. It is found that the decrease in the temperature-dependent fluid viscosity makes the velocity decrease with the increasing distance of the stretching sheet. At a particular point of the sheet, the fluid velocity decreases but the temperature increases with the decreasing viscosity. The impact of the thermophoresis particle deposition plays an important role in the concentration boundary layer. The obtained results are presented graphically and discussed.     

Effects of chemical reaction on mixed convection flow of a polar fluid through a porous medium in the presence of internal heat generation

This paper reports a detailed numerical investigation on mixed convection flow of a polar fluid through a porous medium due to the combined effects of thermal and mass diffusion. The energy equation accounts for heat generation or absorption, while the nth order homogeneous chemical reaction between the fluid and the diffusing species is included in the mass diffusion equation. The governing equations of the linear momentum, angular momentum, energy and concentration are obtained in a non-similar form by introducing a suitable group of transformations. The final set of non-similar coupled non-linear partial differential equations is solved using an implicit finite-difference scheme in combination with quasi-linearization technique. The effects of various parameters on the velocity, angular velocity, temperature and concentration fields are investigated. Numerical results for the skin friction coefficient, wall stress of angular velocity, Nusselt number and Sherwood number are also presented.     

Solitary wave solutions of the nonlinear generalized Pochhammer–Chree and regularized long wave equations

In this paper, we implement some fast and high accuracy numerical algorithms to obtain the solitary wave solutions of generalized Pochhammer?CChree (PC) and regularized long wave (RLW) equations. We employ the discrete Fourier transform to discretize the original partial differential equations (PDEs) in space and obtain a system of ordinary differential equations (ODEs) in Fourier space which will be solved with fourth order time-stepping methods. The proposed methods are fast and accurate due to the use of the fast Fourier transform in combination with explicit fourth-order time stepping methods. For RLW equation we investigate the propagation of a single solitary and interaction of two and three solitary waves. Moreover, three invariants of motion (mass, energy, and momentum) are evaluated to determine the conservation properties of the problem, and the numerical schemes lead to accurate results. The numerical results are compared with analytical solutions and with those of other recently published methods to confirm the accuracy and efficiency of the presented schemes.     

Bounded solutions of linear differential equations in a banach space

For a linear inhomogeneous differential equation in a Banach space, we find a criterion for the existence of solutions that are bounded on the entire real axis under the assumption that the homogeneous equation admits an exponential dichotomy on the semiaxes. This result is a generalization of the Palmer lemma to the case of infinite-dimensional spaces. We consider examples of countable systems of ordinary differential equations that have bounded solutions. __________ Translated from Neliniini Kolyvannya, Vol. 9, No. 1, pp. 3–14, January–March, 2006.     

An index 0 Differential-Algebraic equation formulation for multibody dynamics: Holonomic constraints

The Lagrange multiplier form of index 3 differential-algebraic equations of motion for holonomically constrained multibody systems is transformed using tangent space generalized coordinates to an index 0 form that is equivalent to an ordinary differential equation. The index 0 formulation includes embedded tolerances that assure satisfaction of position, velocity, and acceleration constraints and is solved using established explicit and implicit numerical integration methods. Numerical experiments with two spatial applications show that the formulation accurately satisfies constraints, preserves invariants due to conservation laws, and behaves as if applied to an ordinary differential equation.     

Scaling Transformation for the Effect of Temperature-Dependent Fluid Viscosity with Thermophoresis Particle Deposition on MHD-Free Convective Heat and Mass Transfer Over a Porous Stretching Surface

This article concerns with a steady two-dimensional flow of an electrically conducting incompressible fluid over a vertical stretching sheet. The flow is permeated by a uniform transverse magnetic field. The fluid viscosity is assumed to vary as a linear function of temperature. A scaling group of transformations is applied to the governing equations. The system remains invariant due to some relations among the parameters of the transformations. After finding three absolute invariants, a third-order ordinary differential equation corresponding to the momentum equation, and two second-order ordinary differential equations corresponding to energy and diffusion equations are derived. The equations along with the boundary conditions are solved numerically. It is found that the decrease in the temperature-dependent fluid viscosity makes the velocity to decrease with the increasing distance of the stretching sheet. At a particular point of the sheet, the fluid velocity decreases with the decreasing viscosity but the temperature increases in this case. Impact of thermophoresis particle deposition in the presence of temperature-dependent fluid viscosity plays an important role on the concentration boundary layer. The results, thus, obtained are presented graphically and discussed.     

Finite strain solutions in compressible isotropic elasticity

Three classes of compressible isotropic elastic solids are introduced, for each of which the strain energy, expressed as a function of the three principal invariants of the stretch tensors, is linear in two of its arguments and nonlinear in the third argument. One of these is the class of harmonic materials. Several deformation fields are examined, for which the governing equations, for general compressible isotropic elastic response, reduce to a nonlinear ordinary differential equation. For the three special classes of materials, this differential equation may be solved in closed form, giving a deformation field which is independent of the material response function, or its solution may be reduced to a single quadrature, involving the nonlinear material response function.     

Axisymmetric spreading of a thin power-law fluid under gravity on a horizontal plane

Separable solutions admitted by a nonlinear partial differential equation modeling the axisymmetric spreading under gravity of a thin power-law fluid on a horizontal plane are investigated. The model equation is reduced to a highly nonlinear second-order ordinary differential equation for the spatial variable. Using the techniques of Lie group analysis, the nonlinear ordinary differential equation is linearized and solved. As a consequence of this linearization, new results are obtained.     

Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

Transient analysis of Cosserat rod with inextensibility and unshearability constraints using the least-squares finite element model

In this study, time-dependent fully discretized least-squares finite element model is developed for the transient response of Cosserat rod having inextensibility and unshearability constraints to simulate a surgical thread in space. Starting from the kinematics of the rod for large deformation, the linear and angular momentum equations along with constraint conditions for the sake of completeness are derived. Then, the α-family of time derivarive approximation is used to reduce the governing equations of motion to obtain a semi-discretized system of equations, which are then fully discretized using the least-squares approach to obtain the non-linear finite element equations. Newton׳s method is utilized to solve the non-linear finite element equations. Dynamic response due to impulse force and time-dependent follower force at the free end of the rod is presented as numerical examples.     

Effects of slip condition on MHD stagnation-point flow over a power-law stretching sheet

The steady two-dimensional magnetohydrodynamic stagnation flow towards a nonlinear stretching surface is studied. The no-slip condition on the solid boundary is replaced with a partial slip condition. A scaling group transformation is used to get the invariants. Using the invariants, a third-order ordinary differential equation corresponding to the momentum is obtained. An analytical solution is obtained in a series form using a homotopy analysis method. Reliability and efficiency of series solutions are shown by the good agreement with numerical results presented in the literature. The effects of the slip parameter, the magnetic field parameter, the velocity ratio parameter, the suction velocity parameter, and the power law exponent on the flow are investigated. The results show that the velocity and shear stress profiles are greatly influenced by these parameters.     

Construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body

We consider the problem of construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body so as to ensure the transition of the rigid body from an arbitrary initial angular position to the required final angular position. For the functionals to be minimized, we use combined performance functionals, one of which characterizes the expenditure of time and of the squared modulus of the angular momentum vector in a given proportion, while the other characterizes the expenditure of time and momentum of the modulus of the angular momentum vector necessary to change the rigid body orientation. The control (the vector of the rigid body angular momentum) is assumed to be bounded in the modulus. The problem is solved by using Pontryagin’s maximum principle and the quaternion differential equation [1, 2] relating the vector of the dynamically symmetric rigid body angular momentum to the quaternion of orientation of the coordinate system rotating with respect to the rigid body about its dynamical symmetry axis at an angular velocity proportional to the angular momentum vector projection on the axis. The use of such a model of rotational motion leads to the problem of optimal control with the moving right end of the trajectory and significantly simplifies the analytic study of the problem of construction of optimal laws of variation in the angular momentum vector, because this model explicitly exploits the body angular momentum quaternion (control) instead of the rigid body absolute angular velocity quaternion. We construct general analytic solutions of the differential equations for the boundary-value problems which form systems of nine nonlinear differential equations. It is shown that the process of solving the differential boundary-value problems is reduced to solving two scalar algebraic transcendental equations.     

Global Weak Solutions to Equations of Motion for Magnetic Fluids

We study the differential system governing the flow of an incompressible ferrofluid under the action of a magnetic field. The system consists of the Navier–Stokes equations, the angular momentum equation, the magnetization equation, and the magnetostatic equations. We prove, by using the Galerkin method, a global in time existence of weak solutions with finite energy of an initial boundary-value problem and establish the long-time behavior of such solutions. The main difficulty is due to the singularity of the gradient magnetic force.      

Invariants of a Remarkable Family of Nonlinear Equations

In classical literature, invariants of families of differentialequations were considered for linear equations only, e.g. the renownedLaplace invariants for linear hyperbolic partial differential equationsand invariants of linear ordinary differential equations with variablecoefficients. The restriction to linear equations was essential inpioneering works of Cockle, Laguerre, Halphen, andForsyth for tackling the problem of invariants of differentialequations. Lie regretted that these authors did not use advantagesprovided by his theory of infinite continuous groups, but he himself didnot undertake further developments in this direction.Recently, the present author considered the possibility hinted byLie's remark and introduced the infinitesimal technique in thetheory of invariants of families of differential equations thatwas lacking in old methods. In consequence, a simple unifiedapproach was developed for calculation of invariants of algebraicand differential equations independent on the assumption oflinearity of equations. It was employed recently for calculationof Laplace type invariants for parabolic equations. Here, themethod is applied to calculation of invariants for the family ofnonlinear equations appearing in the problem on linearization ofnonlinear ordinary differential equations.     

Evolving Graphs by Singular Weighted Curvature

A new notion of solutions is introduced to study degenerate nonlinear parabolic equations in one space dimension whose diffusion effect is so strong at particular slopes of the unknowns that the equation is no longer a partial differential equation. By extending the theory of viscosity solutions, a comparison principle is established. For periodic continuous initial data a unique global continuous solution (periodic in space) is constructed. The theory applies to motion of interfacial curves by crystalline energy or more generally by anisotropic interfacial energy with corners when the curves are the graphs of functions. Even if the driving force term (homogeneous in space) exists, the initial-value problem is solvable for general nonadmissible continuous (periodic) initial data. (Accepted July 5, 1996)     

Estimates of azimuthal numbers associated with elementary elliptic cylinder wave functions

The paper deals with issues related to the construction of solutions, 2 π-periodic in the angular variable, of the Mathieu differential equation for the circular elliptic cylinder harmonics, the associated characteristic values, and the azimuthal numbers needed to form the elementary elliptic cylinder wave functions. A superposition of the latter is one possible form for representing the analytic solution of the thermoelastic wave propagation problem in long waveguides with elliptic cross-section contour. The classical Sturm-Liouville problem for the Mathieu equation is reduced to a spectral problem for a linear self-adjoint operator in the Hilbert space of infinite square summable two-sided sequences. An approach is proposed that permits one to derive rather simple algorithms for computing the characteristic values of the angular Mathieu equation with real parameters and the corresponding eigenfunctions. Priority is given to the application of the most symmetric forms and equations that have not yet been used in the theory of the Mathieu equation. These algorithms amount to constructing a matrix diagonalizing an infinite symmetric pentadiagonal matrix. The problem of generalizing the notion of azimuthal number of a wave propagating in a cylindrical waveguide to the case of elliptic geometry is considered. Two-sided mutually refining estimates are constructed for the spectral values of the Mathieu differential operator with periodic and half-periodic (antiperiodic) boundary conditions.     

Heat transfer characteristics of thin power-law liquid films over horizontal stretching sheet with internal heating and variable thermal coefficient

The effect of internal heating source on the film momentum and thermal transport characteristic of thin finite power-law liquids over an accelerating unsteady horizontal stretched interface is studied. Unlike most classical works in this field, a general surface temperature distribution of the liquid film and the generalized Fourier’s law for varying thermal conductivity are taken into consideration. Appropriate similarity transformations are used to convert the strongly nonlinear governing partial differential equations (PDEs) into a boundary value problem with a group of two-point ordinary differential equations (ODEs). The correspondence between the liquid film thickness and the unsteadiness parameter is derived with the BVP4C program in MATLAB. Numerical solutions to the self-similarity ODEs are obtained using the shooting technique combined with a Runge-Kutta iteration program and Newton’s scheme. The effects of the involved physical parameters on the fluid’s horizontal velocity and temperature distribution are presented and discussed.     

Stability analysis of a fractional viscoelastic plate strip in supersonic flow under axial loading

The stability of a viscoelastic plate strip, subjected to an axial load with the Kelvin–Voigt fractional order constitutive relationship is studied. Based on the classical plate theory, the structural formulation of the plate is obtained by using the Newton’s second law and the aerodynamic force due to the fluid flow is evaluated by piston theory. The Galerkin method is employed to discretize the equation of motion into a set of ordinary differential equations. To determine the stability margin of plate the obtained set of ordinary differential equations are solved using the Laplace transform method. The effects of variation of the governing parameters such as axial force, retardation time, fractional order and boundary conditions on the stability margin of fractional viscoelastic panel are investigated and finally some conclusions are outlined.     

Geometrical methods in the analysis of ordinary differential equations

Summary In several fields of engineering research, particularly in the study of vibrations, electrical circuits and in some problems of fluid mechanics, approximations which lead to linear differential equations are proving inadequate. This circumstance is focussing the attention of research workers and engineers on non-linear problems. This article gives an account, without proofs, but with literature references, of methods for the qualitative integration of non-linear ordinary differential equations of the first order, i.e. for the determination of the pattern of the integral curves of such equations. The use of such geometrical methods becomes necessary in cases when the equation cannot be integrated in closed form. Simple and complex patterns associated with singular points are discussed, and criteria for their classification are given. A method of determining the asymptotic behaviour of the family of solutions is given, and criteria for the existence of closed curves in the family of solutions, as well as the occurrence of limit cycles, are discussed. A brief discussion of the Kronecker index and of the mutual relation between several singular points is added. The text is illustrated with several examples selected from the fields of vibration, compressible fluid flow and electrical circuits.