On order reduction of non-linear equations of mechanics and mathematical physics,new integrable equations and exact solutions

Some classes of non-linear equations of mechanics and mathematical physics are described that admit order reduction through the use of a hydrodynamic-type transformation, where a first-order partial derivative is taken as a new independent variable and a second-order partial derivative is taken as the new dependent variable. The results obtained are used for order reduction of hydrodynamic equations (Navier–Stokes, Euler, and boundary layer) and deriving exact solutions to these equations. Associated Bäcklund transformations are constructed for evolution equations of general form (special cases include Burgers, Korteweg-de Vries, and many other non-linear equations of mathematical physics). A number of new integrable non-linear equations, inclusive of the generalized Calogero equation, are considered.     

Group properties and invariant solutions for infinitesimal transformations of a non-linear wave equation

This paper concerns the infinitesimal group analysis for a second order non-linear wave equation involving non-homogeneous processes. In the first part we characterize the most general expression for the generator of the Lie group. Then we calculate the invariant surfaces in some special cases, obtaining the corresponding ordinary differential equations whose integration allow us to get classes of solutions for the original equation.     

Exact analytic solutions of the porous media and the gas pressure diffusion ODEs in non-linear mechanics

Two kinds of second-order non-linear ordinary differential equations (ODEs) appearing in mathematical physics and non-linear mechanics are analyzed in this paper. The one concerns the Kidder equation in porous media and the second the gas pressure diffusion equation. Both these equations are strongly non-linear including quadratic first-order derivatives (damping terms). By a series of admissible functional transformations we reduce the prescribed equations to Abel's equations of the second kind of the normal form that they do not admit exact analytic solutions in terms of known (tabulated) functions. According to a mathematical methodology recently developed concerning the construction of exact analytic solutions of the above class of Abel's equations, we succeed in performing the exact analytic solutions of both Kidder's and gas pressure diffusion equations. The boundary and initial data being used in the above constructions are in accordance with each specific problem under considerations.     

Group classification of one-dimensional equations of fluids with internal energy depending on the density and the gradient of the density

Group analysis provides a regular procedure for mathematical modeling by classifying differential equations with respect to arbitrary elements. This article presents the group classification of one-dimensional equations of fluids, where the internal energy is a function of the density and the gradient of the density. The equivalence Lie group and the admitted Lie group are provided. The group classification separates all models into 21 different classes according to the admitted Lie group. Invariant solutions of one particular model are obtained.      

Rotationally symmetric internal gravity waves

Many mathematical models formulated in terms of non-linear differential equations can successfully be treated and solved by Lie group methods. Lie group analysis is especially valuable in investigating non-linear differential equations, for its algorithms act here as reliably as for linear cases. The aim of this article is to provide the group theoretical modeling of internal waves in the ocean. The approach is based on a new concept of conservation laws that is utilized to systematically derive the conservation laws of non-linear equations describing propagation of internal waves in the ocean. It was shown in our previous publication that uni-directional internal wave beams can be obtained as invariant solutions of non-linear equations of motion. The main goal of the present publication is to thoroughly analyze another physically significant exact solution, namely the rotationally symmetric solution and the energy carried by this solution. It is shown that the rotationally symmetric solution and its energy are presented by means of a bounded oscillating function.     

Post-buckling behaviour of slender structures with a bi-linear bending moment-curvature relationship

Certain classes of slender structures of complex cross-section or fabricated from specialised materials can exhibit a bi-linear bending moment-curvature relationship that has a strong influence on their global structural behaviour. This condition may be encountered, for instance, in (a) non-linear elastic or inelastic post-buckling problems if the cross-section stiffness may be well approximated by a bi-linear model; (b) multi-layered structures such as stranded cables, power transmission lines, umbilical cables and flexible pipes where the drop in the bending stiffness is associated with an internal friction mechanism. This paper presents a mathematical formulation and an analytical solution for such slender structures with a bi-linear bending moment versus curvature constitutive behaviour and subjected to axial terminal forces. A set of five first-order non-linear ordinary differential equations are derived from considering geometrical compatibility, equilibrium of forces and moments and constitutive equations, with hinged boundary conditions prescribed at both ends, resulting a complex two-point boundary value problem. The variables are non-dimensionalised and solutions are developed for monotonic and unloading conditions. The results are presented in non-dimensional graphs for a range of critical curvatures and reductions in bending stiffness, and it is shown how these parameters affect the structure's post-buckling behaviour.     

Lie symmetry solutions for two-phase mass flows

We apply Lie symmetry method to a set of non-linear partial differential equations, which describes a two-phase rapid gravity mass flow as a mixture of solid particles and viscous fluid down a slope (Pudasaini, J. Geophys. Res. 117 (2012) F03010, 28 pp [1]). In order to systematically explore the mathematical structure and underlying physics of the two-phase mixture flow, we generate several similarity forms in general form and construct self-similar solutions. Our analysis generalizes the results, obtained by applying the Lie symmetry method to relatively simple single-phase pressure-driven gravity mass flows, to the two-phase mass flows that include several dominant driving forces and strong phase-interactions. Analytical and numerical solutions are presented for the symmetry-reduced homogeneous and non-homogeneous systems of equations. Analytical and numerical results show that the new models presented here can adequately describe the dynamics of two-phase debris flows, and produce observable phenomena that are consistent with the physics of the flow. The solutions are strongly dependent on the choice of the symmetry-reduced model, as characterized by the group parameters, and the physical parameters of the flows. These solutions reveal strong non-linear and distinct dynamic evolutions, and phase-interactions between the solid and fluid phases, namely the phase-heights and phase-velocities.     

Intermediate asymptotics,scaling laws and renormalization group in continuum mechanics

Scaling laws and self-similar solutions are very popular concepts in modern continuum mechanics. In the present paper these concepts are analyzed both from the viewpoint of intermediate asymptotics, known in classical mathematical physics and fluid mechanics, and from the viewpoint of the renormalization group technique, known in modern theoretical physics. The definition of the renormalization group is proposed, related to the intermediate asymptotics with incomplete similarity. The general presentation is illustrated by examples of essentially non-linear problems where all analytical properties of the solutions and their asymptotics are rigorously proved, as well by an example from turbulence, where the rigorous problem statement is missing. General lecture delivered at the 11th Italian National Congress of Theoretical and Applied Mechanics (AIMETA), Trento, Sept./Oct. 1992.     

Steady flows of non-Newtonian fluids past a porous plate with suction or injection

The problem of the steady flow of three classes of non-linear fluids of the differential type past a porous plate with uniform suction or injection is studied. The flow which is studied is the counterpart of the classical ‘asymptotic suction’ problem, within the context of the non-Newtonian fluid models. The non-linear differential equations resulting from the balance of momentum and mass, coupled with suitable boundary conditions, are solved numerically either by a finite difference method or by a collocation method with a B-spline function basis. The manner in which the various material parameters affect the structure of the boundary layer is delineated. The issue of paucity of boundary conditions for general non-linear fluids of the differential type, and a method for augmenting the boundary conditions for a certain class of flow problems, is illustrated. A comparison is made of the numerical solutions with the solutions from a regular perturbation approach, as well as a singular perturbation.     

On a Korteweg–de Vries-like equation with higher degree non-linearity

We deal with a non-linear partial differential equation which has been widely investigated owing to its applications in quantum field theory, as well as plasma and solid-state physics. It is the matter of a third order KdV-like equation with higher degree non-linearity in the coefficient of the transport term; it can be derived from a Lagrangian or an Hamiltonian density. In the current literature specific attention has been devoted to the search for traveling-wave solutions, depending upon a positive parameter v, which assesses the speed of the solitary wave. The velocity v is always assumed to be constant, as its dependence on the wave-amplitude is neglected in the mathematical model. In this context, Coffey [On series expansions giving closed-form solutions of Korteweg–de Vries-like equations, SIAM J. Appl. Math. 50 (6) (1990) 1580–1592] exploits an algebraic recursive technique to obtain these solutions in closed form for particular values of v. The aim of this paper is to extend these results by showing that closed-form solutions are achievable for every value of v: to this purpose we supply a proper mathematical framework for these issues by taking into account a suitable special function, namely an elliptic function in the sense of Weierstraß. Furthermore we obtain two classes of the so-called kink solutions, see [M.W. Coffey, On series expansions giving closed-form solutions of Korteweg–de Vries-like equations, SIAM J. Appl. Math. 50 (6) (1990) 1580–1592; B. Dey, Domain wall solutions of KdV-like equations with higher order non-linearity, J. Phys. A 19 (1) (1986) L9–L12], as well as an exponential development of the general solution, for which we prove the convergence. Eventually we show how to implement the resulting functions by means of a symbolic manipulation program.     

A study of the interior ballistic equations

Summary A system of non-linear differential equations, typical of the interior ballistics equations for orthodox guns, is studied. Characteristics of the mathematical pressure curves are obtained for the entire spectrums of the ballistic parameter, burning exponent, and geometric parameter (characteristic of the form function). The mathematical solutions are found to distinguish sharply between regressive and progressive burning geometries, and a constant pressure solution is obtained (a piezometric efficiency of one) for what is defined as the ideally progressive geometry. The approximate validity of isothermal solutions, established by Clemmow for a linear burning law, is extended to the general case.The work on this paper was sponsored by the Office of Ordnance Research, U. S. Army, Contract DA-36-034 ORD-2733 RD.     

Discontinuity formation in solutions of homogeneous non-linear hyperbolic equations possessing smooth initial data

Large classes of non-linear equations, at which previous breakdown theories have been aimed, are obtainable by differentiation of first order equations. The general solutions of these first order equations are also solutions of the corresponding second order equations. These are displayed and employed to calculate the critical time for singularity occurrence. Examples are discussed from gas dynamics, shallow water waves, wave propagation in solids, and electrical transmission lines. This method, when applicable, is simple and yields results which agree with those obtainable from the Ludford and Lax-Jeffrey theories.     

Axisymmetric Vortex Fluid Flow in a Long Elastic Tube

A mathematical model for axisymmetric eddy motion of a perfect incompressible fluid in a long tube with thin elastic walls is proposed. Necessary and sufficient conditions for hyperbolicity of the system of equations of motion for flows with monotonic radial velocity profiles are formulated. The propagation velocities of the characteristics of the system under study and the characteristic shape of this system are calculated. The existence of simple waves continuously attached to a given steady shear flow is proved. The group of transformations admitted by the system is found, and submodels that determine invariant solutions are given. By integrating factorsystems, new classes of exact solutions of equations of motion are found.     

Nonlinear three-dimensional stretched flow of an Oldroyd-B fluid with convective condition,thermal radiation,and mixed convection

The effect of non-linear convection in a laminar three-dimensional Oldroyd-B fluid flow is addressed. The heat transfer phenomenon is explored by considering the non-linear thermal radiation and heat generation/absorption. The boundary layer assumptions are taken into account to govern the mathematical model of the flow analysis. Some suitable similarity variables are introduced to transform the partial differential equations into ordinary differential systems. The Runge-Kutta-Fehlberg fourth-and fifth-order techniques with the shooting method are used to obtain the solutions of the dimensionless velocities and temperature. The effects of various physical parameters on the fluid velocities and temperature are plotted and examined. A comparison with the exact and homotopy perturbation solutions is made for the viscous fluid case, and an excellent match is noted. The numerical values of the wall shear stresses and the heat transfer rate at the wall are tabulated and investigated. The enhancement in the values of the Deborah number shows a reverse behavior on the liquid velocities. The results show that the temperature and the thermal boundary layer are reduced when the nonlinear convection parameter increases. The values of the Nusselt number are higher in the non-linear radiation situation than those in the linear radiation situation.     

Non-linear in-plane postbuckling of arches with rotational end restraints under uniform radial loading

This paper presents a thorough and comprehensive investigation of non-linear buckling and postbuckling analyses of pin-ended shallow circular arches subjected to a uniform radial load and which have equal elastic rotational end-restraints. The differential equations of equilibrium for non-linear buckling and postbuckling are established based on a virtual work approach. Exact solutions for the non-linear bifurcation, limit point and lowest buckling loads are obtained; in particular, exact solutions for the non-linear postbuckling equilibrium paths are derived. The criteria for switching between fundamental buckling and postbuckling modes are developed in terms of critical values of a geometric parameter for an arch, with exact solutions for these critical values of geometric parameter being obtained. Analytical solutions of non-linear buckling and postbuckling problems for arches with rotational end-restraints are of great interest, since they constitute one of the very few closed-form analyses of buckling and postbuckling behaviour of continuous structural systems. These exact solutions are a contribution to the non-linear structural mechanics of arches, as well as providing useful benchmark solutions for verifying non-linear numerical analyses.     

A unified compatibility method for exact solutions of non-linear flow models of Newtonian and non-Newtonian fluids

Criteria are established for higher order ordinary differential equations to be compatible with lower order ordinary differential equations. Necessary and sufficient compatibility conditions are derived which can be used to construct exact solutions of higher order ordinary differential equations subject to lower order equations. We provide the connection to generalized groups through conditional symmetries. Using this approach of compatibility and generalized groups, new exact solutions of non-linear flow problems arising in the study of Newtonian and non-Newtonian fluids are derived. The ansatz approach for obtaining exact solutions for non-linear flow models of Newtonian and non-Newtonian fluids is unified with the application of the compatibility and generalized group criteria.     

Hypoplasticity theory for granular materials—II: Three-dimensional axially symmetric exact solutions

In part I of this paper, we consider the governing equations of hypoplasticity theory for two-dimensional steady quasi-static plane strain compressible gravity flow and determine some exact analytical solutions applying for certain special cases. Similarly, for the three-dimensional situation considered here in part II, we undertake a similar mathematical investigation to determine some simple solutions of the governing equations for three-dimensional steady quasi-static axially symmetric compressible gravity flow for hypoplastic granular materials. We again find that for certain special cases, we are able to determine some exact solutions for the stress and velocity profiles. We comment that hypoplasticity theory generally gives rise to complicated systems of coupled non-linear differential equations, for which the determination of any analytical solutions is not a trivial matter, and that the solutions determined here might be exploited as benchmarks for full numerical schemes.     

New solutions for surface tension driven spreading of a thin film

The standard fourth-order non-linear PDE modelling the flow of thin fluid film subject to surface tension is studied. The Lie group method is used to reduce the model equation from a fourth-order PDE to a fourth-order ODE. Analytical solutions are obtained for certain cases. Where analytical progress cannot be made, we determine numerical solutions.     

Euler flow solutions for transonic shock wave–boundary layer interaction

Steady 2D Euler flow computations have been performed for a wind tunnel section, designed for research on transonic shock wave–boundary layer interaction. For the discretization of the steady Euler equations, an upwind finite volume technique has been applied. The solution method used is collective, symmetric point Gauss–Seidel relaxation, accelerated by non-linear multigrid. Initial finest grid solutions have been obtained by nested iteration. Automatic grid adaptation has been applied for obtaining sharp shocks. An indication is given of the mathematical quality of four different boundary conditions for the outlet flow. Two transonic flow solutions with shock are presented: a choked and a non-choked flow. Both flow solutions show good shock capturing. A comparison is made with experimental results.