Unsteady Incompressible Three-Dimensional Asymmetric Stagnation-Point Boundary Layers

The unsteady laminar incompressible boundary-layer flow near the three-dimensional asymmetric stagnation point has been studied under the assumptions that the free-stream velocity, wall temperature, and surface mass transfer vary arbitrarily with time. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. It is found that in contrast with the symmetric flow, the maximum heat transfer occurs away from the stagnation point due to the decrease in the boundary-layer thickness. The effect of the variation of the wall temperature with time on heat transfer is strong. The skin friction and heat transfer due to asymmetric flow only are comparatively less affected by the mass transfer as compared to those of symmetric flow.

Unsteady Natural Convection Flow over a Heated Cylinder Buried in a Fluid Saturated Porous Medium

Transient natural convection flow on a heated cylinder buried in a semi-infinite liquid-saturated porous medium has been studied. The unsteadiness in the problem arises due to the cylinder which is heated (cooled) suddenly and then maintained at that temperature. The coupled partial differential equations governing the flow and heat transfer are cast into stream function-temperature formulation, and the solutions are obtained from the initial time to the time when steady state is reached. The heat transfer is found to change significantly with increasing time in a small time interval immediately after the start of the impulsive change, and steady state is reached after some time. The average Nusselt number is found to increase with Rayleigh number When the surface of the cylinder is suddenly cooled, there is a change in the direction of the heat transfer in a small time interval immediately after the start of the impulsive change in the surface temperature;however when the surface temperature is suddenly increased, no such phenomenon is observed.

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.

Global lopsided instability in a purely stellar galactic disc

It is shown that pure exponential discs in spiral galaxies are capable of supporting slowly varying discrete global lopsided modes, which can explain the observed features of lopsidedness in the stellar discs. Using linearized fluid dynamical equations with the softened self-gravity and pressure of the perturbation as the collective effect, we derive self-consistently a quadratic eigenvalue equation for the lopsided perturbation in the galactic disc. On solving this, we find that the ground-state mode shows the observed characteristics of the lopsidedness in a galactic disc, namely the fractional Fourier amplitude A(1), increases smoothly with the radius. These lopsided patterns precess in the disc with a very slow pattern speed with no preferred sense of precession. We show that the lopsided modes in the stellar disc are long-lived because of a substantial reduction (approximately a factor of 10 compared to the local free precession rate) in the differential precession. The numerical solution of the equations shows that the groundstate lopsided modes are either very slowly precessing stationary normal mode oscillations of the disc or growing modes with a slow growth rate depending on the relative importance of the collective effect of the self-gravity. N-body simulations are performed to test the spontaneous growth of lopsidedness in a pure stellar disc. Both approaches are then compared and interpreted in terms of long-lived global m = 1 instabilities, with almost zero pattern speed.

Gelatin hydrogel electrolytes and their application to electrochemical supercapacitors

Gelatin hydrogel electrolytes (GHEs) with varying NaCl concentrations have been prepared by cross-linking an aqueous solution of gelatin with aqueous glutaraldehyde and characterized by scanning electron microscopy, differential scanning calorimetry, cyclic voltammetry, electrochemical impedance spectroscopy, and galvanostatic chronopotentiometry. Glass transition temperatures for GHEs range between 339.6 and 376.9 K depending on the dopant concentration. Ionic conductivity behavior of GHEs was studied with varying concentrations of gelatin, glutaraldehyde, and NaCl, and found to vary between 10(-3) and 10(-1) S cm(-1). GHEs have a potential window of about 1 V. Undoped and 0.25 N NaCl-doped GHEs follow Arrhenius equations with activation energy values of 1.94 and 1.88 x 10(-4) eV, respectively. Electrochemical supercapacitors (ESs) employing these GHEs in conjunction with Black Pearl Carbon electrodes are assembled and studied. Optimal values for capacitance, phase angle, and relaxation time constant of 81 F g(-1), 75 degrees, and 0.03 s are obtained for 3 N NaCl-doped GHE, respectively. ES with pristine GHE exhibits a cycle life of 4.3 h vs 4.7 h for the ES with 3 N NaCl-doped GHE. (c) 2007 The Electrochemical Society.

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.

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.

The buoyant plume above a heat source

The boundary-layer type conservation equations of mass, momentum and energy for the steady free turbulent flow in gravitational convection over heat sources are set up for both two-dimensional and axisymmetric cases. These are reduced to ordinary differential equations in a similarity parameter by suitable transformations. The three classical hypotheses of turbulent diffusion-the Constant Exchange Coefficient hypothesis, Prandtl's Momentum Transfer theory and Taylor's Vorticity Transfer theory-are then incorporated into these equations in succession. The resulting equations are solved numerically and the results compared with some experimental results on gravitational convection over heat sources reported by Rouse et al.

Free convection heat transfer between vertical parallel plates

A new mathematical model for the solution of the problem of free convection heat transfer between vertical parallel flat isothermal plates under isothermal boundary conditions, has been presented. The set of boundary layer equations used in the model are transformed to nonlinear coupled differential equations by similarity type variables as obtained by Ostrach for vertical flat plates in an infinite fluid medium. By utilising a parameter ηw* to represent the outer boundary, the governing differential equations are solved numerically for parametric values of Pr = 0.733. 2 and 3, and ηw* = 0.1, 0.5, 1, 2, 3, 4, ... and 8.0. The velocity and temperature profiles are presented. Results indicate that ηw* can effectively classify the system into (1) thin layers where conduction predominates, (2) intermediate layers and (3) thick layers whose results can be predicted by the solutions for vertical flat plates in infinite fluid medium. Heat transfer correlations are presented for the 3 categories. Several experimental and analytical results available in the literature agree with the present correlations.

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.

Free convection heat transfer from vertical cylinders part I: Power law surface temperature variation

This paper reports on the investigations of laminar free convection heat transfer from vertical cylinders and wires whose surface temperature varies along the height according to the relation TW - T∞ = Nxn. The set of boundary layer partial differential equations and the boundary conditions are transformed to a more amenable form and solved by the process of successive substitution. Numerical solutions of the first approximated equations (two-point nonlinear boundary value type of ordinary differential equations) bring about the major contribution to the problem (about 95%), as seen from the solutions of higher approximations. The results reduce to those for the isothermal case when n=0. Criteria for classifying the cylinders into three broad categories, viz., short cylinders, long cylinders and wires, have been developed. For all values of n the same criteria hold. Heat transfer correlations obtained for short cylinders (which coincide with those of flat plates) are checked with those available in the literature. Heat transfer and fluid flow correlations are developed for all the regimes.

A weighted mean square method of linearization in non-linear oscillations

The paper deals with a linearization technique in non-linear oscillations for systems which are governed by second-order non-linear ordinary differential equations. The method is based on approximation of the non-linear function by a linear function such that the error is least in the weighted mean square sense. The method has been applied to cubic, sine, hyperbolic sine, and odd polynomial types of non-linearities and the results obtained are more accurate than those given by existing linearization methods.

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.

Reliability models for existing structures based on dynamic state estimation and data based asymptotic extreme value analysis

The problem of time variant reliability analysis of existing structures subjected to stationary random dynamic excitations is considered. The study assumes that samples of dynamic response of the structure, under the action of external excitations, have been measured at a set of sparse points on the structure. The utilization of these measurements m in updating reliability models, postulated prior to making any measurements, is considered. This is achieved by using dynamic state estimation methods which combine results from Markov process theory and Bayes' theorem. The uncertainties present in measurements as well as in the postulated model for the structural behaviour are accounted for. The samples of external excitations are taken to emanate from known stochastic models and allowance is made for ability (or lack of it) to measure the applied excitations. The future reliability of the structure is modeled using expected structural response conditioned on all the measurements made. This expected response is shown to have a time varying mean and a random component that can be treated as being weakly stationary. For linear systems, an approximate analytical solution for the problem of reliability model updating is obtained by combining theories of discrete Kalman filter and level crossing statistics. For the case of nonlinear systems, the problem is tackled by combining particle filtering strategies with data based extreme value analysis. In all these studies, the governing stochastic differential equations are discretized using the strong forms of Ito-Taylor's discretization schemes. The possibility of using conditional simulation strategies, when applied external actions are measured, is also considered. The proposed procedures are exemplifiedmby considering the reliability analysis of a few low-dimensional dynamical systems based on synthetically generated measurement data. The performance of the procedures developed is also assessed based on a limited amount of pertinent Monte Carlo simulations. (C) 2010 Elsevier Ltd. All rights reserved.

Unsteady MHD forced flow due to a point sink

An analysis is performed to study the unsteady laminar incompressible boundary-layer flow of an electrically conducting fluid in a cone due to a point sink with an applied magnetic field. The unsteadiness in the flow is considered for two types of motion, viz. the motion arising due to the free stream velocity varying continuously with time and the transient motion occurring due to an impulsive change either in the strength of the point sink or in the wall temperature. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The magnetic field increases the skin friction but reduces heat transfer. The heat transfer and temperature field are strongly influenced by the viscous dissipation and Prandtl number. The velocity field is more affected at the early stage of the transient motion, caused by an impulsive change in the strength of the point sink, as compared to the temperature field. When the transient motion is caused by a sudden change in the wall temperature, both skin friction and heat transfer take more time to reach a new steady state. The transient nature of the flow and heat transfer is active for a short time in the case of suction and for a long time in the case of injection. The viscous dissipation prolongs the transient behavior of the flow.