Compressible axially symmetric laminar boundary layer flow in the stagnation region of a blunt body in the presence of magnetic field

In this paper, the steady laminar viscous hypersonic flow of an electrically conducting fluid in the region of the stagnation point of an insulating blunt body in the presence of a radial magnetic field is studied by similarity solution approach, taking into account the variation of the product of density and viscosity across the boundary layer. The two coupled non-linear ordinary differential equations are solved simultaneously using Runge-Kutta-Gill method. It has been found that the effect of the variation of the product of density and viscosity on skin friction coefficient and Nusselt number is appreciable. The skin friction coefficient increases but Nusselt number decreases as the magnetic field or the total enthalpy at the wall increases

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.

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.

Viscous flow computations using a meshless solver, LSFD-U

The Upwind-Least Squares Finite Difference (LSFD-U) scheme has been successfully applied for inviscid flow computations. In the present work, we extend the procedure for computing viscous flows. Different ways of discretizing the viscous fluxes are analysed for the positivity, which determines the robustness of the solution procedure. The scheme which is found to be more positive is employed for viscous flux computation. The numerical results for validating the procedure are presented.

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.

Three-Dimensional Unsteady Flow With Heat and Mass Transfer Over a Continuous Stretching Surface

A numerical solution of the unsteady boundary layer equations under similarity assumptions is obtained. The solution represents the three-dimensional unsteady fluid motion caused by the time-dependent stretching of a flat boundary. It has been shown that a self-similar solution exists when either the rate of stretching is decreasing with time or it is constant. Three different numerical techniques are applied and a comparison is made among them as well as with earlier results. Analysis is made for various situations like deceleration in stretching of the boundary, mass transfer at the surface, saddle and nodal point flows, and the effect of a magnetic field. Both the constant temperature and constant heat flux conditions at the wall have been studied.

Second-order boundary-layer effects for the unsteady laminar incompressible three-dimensional stagnation-point flow

All the second-order boundary-layer effects on the unsteady laminar incompressible flow at the stagnation-point of a three-dimensional body for both nodal and saddle point regions have been studied. It has been assumed that the free-stream velocity, wall temperature and mass transfer vary arbitrarily with time. The effect of the Prandtl number has been taken into account. The partial differential equations governing the flow have been derived for the first time and then solved numerically unsteady free-stream velocity distributions, the nature of the using an implicit finite-difference scheme. It is found that the stagnation point and the mass transfer strongly affect the skin friction and heat transfer whereas the effects of the Prandtl number and the variation of the wall temperature with time are only on the heat transfer. The skin friction due to the combined effects of first- and second-order boundary layers is less than the skin friction due to, the first-order boundary layers whereas the heat transfer has the opposite behaviour. Suction increases the skin friction and heat transfer but injection does the opposite

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 mathematical analysis of nonlinear waves in a fluid filled visco-elastic tube

Our investigations in this paper are centred around the mathematical analysis of a ldquomodal waverdquo problem. We have considered the axisymmetric flow of an inviscid liquid in a thinwalled viscoelastic tube under certain simplifying assumptions. We have first derived the propagation space equations in the long wave limit and also given a general procedure to derive these equations for arbitrary wave length, when the flow is irrotational. We have used the method of operators of multiple scales to derive the nonlinear Schrödinger equation governing the modulation of periodic waves and we have elaborated on the ldquolong modulated wavesrdquo and the ldquomodulated long wavesrdquo. We have also examined the existence and stability of Stokes waves in this system. This is followed by a discussion of the progressive wave solutions of the long wave equations. One of the most important results of our paper is that the propagation space equations are no longer partial differential equations but they are in terms of pseudo-differential operators.Die vorliegenden Untersuchungen beziehen sich auf die mathematische Behandlung des ldquorModalwellenrdquo-Problems. Die achsensymmetrische Strömung einer nichtviskosen Flüssigkeit in einem dünnwandigen viskoelastischen Rohr, unter bestimmten vereinfachenden Annahmen, wird betrachtet. Zuerst werden die Gleichungen des Ausbreitungsraumes im Langwellenbereich abgeleitet und eine allgemeine Methode zur Herleitung dieser Gleichungen für beliebige Wellenlängen bei nichtrotierender Strömung angegeben. Eine Operatorenmethode mit multiplem Maßstab wird verwendet zur Herleitung der nichtlinearen Schrödinger-Gleichung für die Modulation der periodischen Wellen, und die ldquorlangmodulierten Wellenrdquo sowie die ldquormodulierten Langwellenrdquo werden aufgezeigt. Weiters wird die Existenz und die Stabilität der Stokes-Wellen im System untersucht. Anschließend werden die progressiven Wellenlösungen der Langwellengleichungen diskutiert. Eines der wichtigsten Ergebnisse dieser Arbeit ist, daß die Gleichungen des Ausbreitungsraumes keine partiellen Differentialgleichungen mehr sind, sondern Ausdrücke von Pseudo-Differentialoperatoren.

Low reynolds number flow of a dilute suspension in slowly varying tubes

We have discussed here the flow of a dilute suspension of rigid particles in Newtonian fluid in slowly varying tubes characterized by a small parameter ε. Solutions are presented in the form of asymptotic expansions in powers of ε. The effect of the suspension on the fluid is described by two parameters β and γ which depend on the volume fraction of the particles which we assume to be small. It is found that the presence of the particles accelerate the process of eddy formation near the constriction and shifts the point of separation.

Unsteady mixed convection flow in stagnation region adjacent to a vertical surface

The unsteady two-dimensional laminar mixed convection flow in the stagnation region of a vertical surface has been studied where the buoyancy forces are due to both the temperature and concentration gradients. The unsteadiness in the flow and temperature fields is caused by the time-dependent free stream velocity. Both arbitrary wall temperature and concentration, and arbitrary surface heat and mass flux variations have been considered. The Navier-Stokes equations, the energy equation and the concentration equation, which are coupled nonlinear partial differential equations with three independent variables, have been reduced to a set of nonlinear ordinary differential equations. The analysis has also been done using boundary layer approximations and the difference between the solutions has been discussed. The governing ordinary differential equations for buoyancy assisting and buoyancy opposing regions have been solved numerically using a shooting method. The skin friction, heat transfer and mass transfer coefficients increase with the buoyancy parameter. However, the skin friction coefficient increases with the parameter lambda, which represents the unsteadiness in the free stream velocity, but the heat and mass transfer coefficients decrease. In the case of buoyancy opposed flow, the solution does not exist beyond a certain critical value of the buoyancy parameter. Also, for a certain range of the buoyancy parameter dual solutions exist.

Boundary layer swirling flow through a conical hydrocyclone

The paper deals with the flow and heat-transfer problem of a steady axisymmetric laminar incompressible boundary layer swirling flow of a fluid through a conical hydrocyclone. The implicit finitedifference scheme is used to solve the partial differential equations governing the flow. The effect of swirl is found to be more pronounced on the longitudinal skin friction than on the tangential skin friction and heat transfer. The skin friction and heat transfer increase with swirl or with longitudinal distance. Swirl also gives rise to velocity overshoot in the longitudinal velocity profiles and the magnitude of the velocity overshoot increases as the swirl parameter increases. The results are found to be in good agreement with those of the local nonsimilarity and momentum integral methods but they differ appreciably from those of the local similarity method except for the longitudinal skin friction which is fairly in good agreement with that of the local similarity method.Die Arbeit beschäftigt sich mit der Strömung und dem Wärmeübergang in einem konischen Zyklon unter der Voraussetzung stationärer, achsensymmetrischer, laminarer, inkompressibler Grenzschichtströmung. Ein implizites Differenzenverfahren wird benutzt, um die partiellen Differentialgleichungen zu lösen. Der Einfluß des Dralls ist besonders ausgeprägt auf die longitudinale Komponente der Oberflächenreibung, weniger dagegen bei der tangentialen Komponente und beim Wärmeübergang. Die Oberflächenreibung und der Wärmeübergang nehmen zu mit dem Drall, sowie mit dem longitudinalen Abstand. Der Drall erzeugt ein Überschießen der Geschwindigkeit in der longitudinalen Abstand. Der Drall erzeugt ein Überschießen der Geschwindigkeit in der Längsrichtung. Die Größe des Überschusses nimmt mit wachsendem Drallparameter zu. Die Resultate stimmen gut mit den Ergebnissen der Theorie der lokalen Nichtähnlichkeit und der Impulsintegralmethode überein. Dagegen weichen sie mit Ausnahme der longitudinalen Komponente der Oberflächenreibung beträchtlich von den Resultaten der Theorie der lokalen Ähnlichkeit ab.

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]