Boundary value problems for third-order nonlinear ordinary differential equations

In this paper, we describe how to analyze boundary value problems for third-order nonlinear ordinary differential equations over an infinite interval. Several physical problems of interest are governed by such systems. The seminumerical schemes described here offer some advantages over solutions obtained by using traditional methods such as finite differences, shooting method, etc. These techniques also reveal the analytic structure of the solution function. For illustrative purposes, several physical problems, mainly drawn from fluid mechanics, are considered; they clearly demonstrate the efficiency of the techniques presented here.

Secondary flows induced by the rotation of a sphere or coaxial cones in micropolar fluids

The micropolar fluids like Newtonian and Non-Newtonian fluids cannot sustain a simple shearing motion, wherein only one component of velocity is present. They exhibit both primary and secondary motions when the boundaries are subject to slow rotations. The primary motion, as in Non-Newtonian fluids, characterized by the equation due to Rivlin-Ericksen, Oldroyd, Walters etc., resembles that of Newtonian fluid for slow steady rotation. We further notice that the micro-rotation becomes identically equal to the vorticity present in the fluid and the condition b) of "Wall vorticity" can alone be satisfied at the boundaries. As regards, the secondary motion, we notice that it can be determined by the above procedure for a special class of fluids, namely that for which j0(n2-n3)=4 n3/l2. Moreover for this class of fluids, the micro-rotation is identical with the vorticity of the fluid everywhere. Also the stream function for the secondary flow is identical with that for the Newtonian fluid with a suitable definition of the Reynolds number. In contrast with the Non-Newtonian fluids, characterized by the equation due to Rivlin-Ericksen, Oldroyd, Walters etc., this class of micropolar fluids does not show separation. This is in conformity with the statement of Condiff and Dahler (3) that in any steady flow, internal spin matches the vorticity everywhere provided that (i) spin boundary conditions are satisfied, (ii) body torques and non-conservative body forces are absent, and (iii) inertial and spin-inertial terms are either negligible or vanish identically.

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.

Large deflections of rectangular plates

Approximate solutions for the non-linear bending of thin rectangular plates are presented considering large deflections for various boundary conditions. In the case of stress-free edges, solutions are given for von Kármán's equations in terms of the stress function and the deflection of the plate. In the case of immovable edges, equations are constructed in terms of the three displacements and these are solved. The solution is given by using double series consisting of the appropriate Beam Functions which satisfy the boundary conditions. The differential equations are satisfied by using the orthogonality properties of the series. Numerical results for square plates with uniform lateral load indicate good convergence of the series solution presented here.

Effect of impulsive change of wall velocity and wall temperature in viscous heat-conducting gas

A study of compression waves produced in a viscous heat-conducting gas by the impulsive start of a one-dimensional piston and by the inpulsive change of piston wall temperature is made using Laplace Transform Technique for Prandt1 number unity. Expressions for velocity, temperature and density have also been obtained using small-time expansion procedure in this case. For arbitrary Prandt1 number solutions have been developed using large-time expansion procedure. A number of graphs exhibiting the distribution of the fluid velocity, temperature and density have been drawn.

Extended self-similarity works for the Burgers equation and why

Extended self-similarity (ESS), a procedure that remarkably extends the range of scaling for structure functions in Navier-Stokes turbulence and thus allows improved determination of intermittency exponents, has never been fully explained. We show that ESS applies to Burgers turbulence at high Reynolds numbers and we give the theoretical explanation of the numerically observed improved scaling at both the IR and UV end, in total a gain of about three quarters of a decade: there is a reduction of subdominant contributions to scaling when going from the standard structure function representation to the ESS representation. We conjecture that a similar situation holds for three-dimensional incompressible turbulence and suggest ways of capturing subdominant contributions to scaling.

The effect of viscosity and surface tension on the Taylor instability of the surface of a cavitation bubble*

This paper is devoted to a consideration of the following problem: A spherical mass of fluid of density varrho1, viscosity μ1 and external radius R is surrounded by a fluid of density varrho2 and viscosity μ2.The fluids are immiscible and incompressible. The interface is accelerated radially by g1: to study the effect of viscosity and surface tension on the stability of the interface. By analyzing the problem in spherical harmonics the mathematical problem is reduced to one of solution of the characteristic determinant equation. The particular case of a cavity bubble, where the viscosity μ1 of the fluid inside the bubble is negligible in comparison with the viscosity μ2 of the fluid outside the bubble, is considered in some detail. It is shown that viscosity has a stabilizing role on the interface; and when g1 > T(n − 1) (n + 2)/R2(varrho2 − varrho1) the stabilizing role of both viscosity and surface tension is more pronounced than would result when either of them is taken individually.

Large-scale structures in turbulent and reverse-transitional sink flow boundary layers

Aspects of large-scale organized structures in sink flow turbulent and reverse-transitional boundary layers are studied experimentally using hot-wire anemometry. Each of the present sink flow boundary layers is in a state of 'perfect equilibrium' or 'exact self-preservation' in the sense of Townsend (The Structure of Turbulent Shear Flow, 1st and 2nd edns, 1956, 1976, Cambridge University Press) and Rotta (Progr. Aeronaut. Sci., vol. 2, 1962, pp. 1-220) and conforms to the notion of 'pure wall-flow' (Coles, J. Aerosp. Sci., vol. 24, 1957, pp. 495-506), at least for the turbulent cases. It is found that the characteristic inclination angle of the structure undergoes a systematic decrease with the increase in strength of the streamwise favourable pressure gradient. Detectable wall-normal extent of the structure is found to be typically half of the boundary layer thickness. Streamwise extent of the structure shows marked increase as the favourable pressure gradient is made progressively severe. Proposals for the typical eddy forms in sink flow turbulent and reverse-transitional flows are presented, and the possibility of structural self-organization (i.e. individual hairpin vortices forming streamwise coherent hairpin packets) in these flows is also discussed. It is further indicated that these structural ideas may be used to explain, from a structural viewpoint, the phenomenon of soft relaminarization or reverse transition of turbulent boundary layers when subjected to strong streamwise favourable pressure gradients. Taylor's 'frozen turbulence' hypothesis is experimentally shown to be valid for flows in the present study even though large streamwise accelerations are involved, the flow being even reverse transitional in some cases. Possible conditions, which are required to be satisfied for the safe use of Taylor's hypothesis in pressure-gradient-driven flows, are also outlined. Measured convection velocities are found to be fairly close to the local mean velocities (typically 90% or more) suggesting that the structure gets convected downstream almost along with the mean flow.

Natural modes of a coupled non-linear system

The natural modes of a non-linear system with two degrees of freedom are investigated. The system, which may contain either hard or soft springs, is shown to possess three modes of vibration one of which does not have any counterpart in the linear theory. The stability analysis indicates the existence of seven different modal stability patterns depending on the values of two parameters of non-linearity.

Stability and transition in the flow of polymer solutions

The linear stability analysis of a plane Couette flow of viscoelastic fluid have been studied with the emphasis on two dimensional disturbances with wave number k similar to Re-1/2, where Re is Reynolds number based on maximum velocity and channel width. We employ three models to represent the dilute polymer solution: the classical Oldroyd-B model, the Oldroyd-B model with artificial diffusivity and the non-homogeneous polymer model. The result of the linear stability analysis is found to be sensitive to the polymer model used. While the plane Couette flow is found to be stable to infinitesimal disturbances for the first two models, the last one exhibits a linear instability.

Saturated liquid densities of cryogenic liquids and refrigerants

The saturated liquid density, varrholr, data along the liquid vapour coexistence curve published in the literature for several cryogenic liquids, hydrocarbons and halocarbon refrigerants are fitted to a generalized equation of the following form varrholr = 1 + A(1 − Tr + B(1 − Tr)β The values of β, the index in phase density differences power law, have been obtained by means of two approaches namely statistical treatment of saturated fluid phase density difference data and the existence of a maximum in T(varrho1 − varrhov) along the saturation curve. Values of the constants A and B are determined utilizing the fact that Tvarrho1 has a maximum at a characteristic temperature T. Values of A, B and β are tabulated for Ne, Ar, Kr, Xe, N2, O2, methane, ethane, propane, iso-butane, n-butane, propylene, ethylene, CO2, water, ammonia, refrigerants-11, 12, 12B1, 13, 13B1, 14, 21, 22, 23, 32, 40, 113, 114, 115, 142b, 152a, 216, 245 and azeotropes R-500, 502, 503, 504. The average error of prediction is less than 2%.

Stochastic Analysis of Yielding System

The probability distribution of the instantaneous incremental yield of an inelastic system is characterized in terms of a conditional probability and average rate of crossing. The detailed yield statistics of a single degree-of-freedom elasto-plastic system under a Gaussian white noise are obtained for both nonstationary and stationary response. The present analysis indicates that the yield damage is sensitive to viscous damping. The spectra of mean and mean square damage rate are presented.

Unsteady laminar compressible boundary-layer flow at a three-dimensional stagnation point

Unsteady laminar compressible boundary-layer flow with variable properties at a three-dimensional stagnation point for both cold and hot walls has been studied for the case when the velocity of the incident stream varies arbitrarily with time. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. Computations have been carried out for two particular unsteady free-stream velocity distributions: (i) an accelerating stream and (ii) a fluctuating stream. The results indicate that the variation of the density-viscosity product across the boundary layer, the wall temperature and the nature of stagnation point significantly affect the skin friction and heat transfer.

A study on boundary-layer transition induced by free-stream turbulence

Boundary-layer transition at different free-stream turbulence levels has been investigated using the particle-image velocimetry technique. The measurements show organized positive and negative fluctuations of the streamwise fluctuating velocity component, which resemble the forward and backward jet-like structures reported in the direct numerical simulation of bypass transition. These fluctuations are associated with unsteady streaky structures. Large inclined high shear-layer regions are also observed and the organized negative fluctuations are found to appear consistently with these inclined shear layers, along with highly inflectional instantaneous streamwise velocity profiles. These inflectional velocity profiles are similar to those in the ribbon-induced boundary-layer transition. An oscillating-inclined shear layer appears to be the turbulent spot-precursor. The measurements also enabled to compare the actual turbulent spot in bypass transition with the simulated one. A proper orthogonal decomposition analysis of the fluctuating velocity field is carried out. The dominant flow structures of the organized positive and negative fluctuations are captured by the first few eigenfunction modes carrying most of the fluctuating energy. The similarity in the dominant eigenfunctions at different Reynolds numbers suggests that the flow prevails its structural identity even in intermittent flows. This analysis also indicates the possibility of the existence of a spatio-temporal symmetry associated with a travelling wave in the flow.

Rocking response of rectangular rigid blocks under random noise base excitations

The response of a rigid rectangular block resting on a rigid foundation and acted upon simultaneously by a horizontal and a vertical random white-noise excitation is considered. In the equation of motion, the energy dissipation is modeled through a viscous damping term. Under the assumption that the body does not topple, the steady-state joint probability density function of the rotation and the rotational velocity is obtained using the Fokker-Planck equation approach. Closed form solution is obtained for a specific combination of system parameters. A more general but approximate solution to the joint probability density function based on the method of equivalent non-linearization is also presented. Further, the problem of overturning of the block is approached in the framework of the diffusion methods for first passage failure studies. The overturning of the block is deemed incipient when the response trajectories in the phase plane cross the separatrix of the conservative unforced system. Expressions for the moments of first passage time are obtained via a series solution to the governing generalized Pontriagin-Vitt equations. Numerical results illustra- tive of the theoretical solutions are presented and their validity is examined through limited amount of digital simulations.