Flow and heat transfer in a stagnation-point flow over a stretching sheet with a magnetic field

The flow and heat transfer problem in the boundary layer induced by a continuous moving surface is important in many manufacturing processes in industry such as the boundary layer along material handling conveyers, the aerodynamic extrusion of plastic sheet, the cooling of an infinite metalic plate in a cooling bath (which may also be electrolyte). Glass blowing, continuous casting and spinning of fibres also involve the flow due to a stretching surface. Sakiadis [1] was the first to study the flow induced by a semi-infinite moving wall in an ambient fluid. On the other hand, Crane [2] first studied the flow over a linearly stretching sheet in an ambient fluid. Subsequently, Crane [3] also investigated the corresponding heat transfer problem. Since then several authors [4-8] have studied various aspects of this problem such as the effects of mass transfer, variable wall temperature, constant heat flux, magnetic field etc. Recently, Andersson [9] has obtained an exact solution of the Navier-Stokes equations for the MHD flow over a linearly stretching sheet in an ambient fluid. Also Chiam [10] has studied the heat transfer with variable thermal conductivity on a stretching sheet when the velocities of the sheet and the free stream are equal.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

A numerical study of thermocapillary flow in a rectangular cavity during laser melting

The surface tension gradient driven flow that occurs during laser melting has been studied. The vorticity-streamfunction form of the Navier-Stokes equations and the energy equation has been solved by the ‘Alternative Direction Implicit’ method. It has been shown that the inertia forces in the melt strongly influence the flow pattern in the melt. The convection in the melt modifies the isotherms in the melt at high surface tension Reynolds number and high Prandtl number. The buoyancy driven flow has been shown to be negligible compared to the surface tension gradient driven flow in laser melting.

Effect of inclination on laminar film condensation

The effect of inclination on laminar film condensation over and under isothermal flat plates is investigated analytically. The complete set of Navier Stokes equations in two dimensions is considered. Analysed as a perturbation problem, the zero-order perturbation represents the boundary layer equations. First and second order perturbations are solved to bring about the leading edge effects. Corresponding velocity and temperature profiles are presented. The results show decrease in heat transfer with larger ∥inclinations∥ from the vertical. Comparison with experimental data of Gerstmann and Griffith indicates a closer agreement of the present results than the analytical results of the same authors.

Boundary-layer effects on internal flow choking in dual-thrust solid rocket motors

Theoretical studies have been carried out to examine internal flow choking in the inert simulators of a dual-thrust motor. Using a two-dimensional k-omega turbulence model, detailed parametric studies have been carried out to examine aerodynamic choking and the existence of a fluid throat at the transition region during the startup transient of dual-thrust motors. This code solves standard k-omega turbulence equations with shear flow corrections using a coupled second-order-implicit unsteady formulation. In the numerical study, a fully implicit finite volume scheme of the compressible, Reynolds-averaged, Navier-Stokes equations is employed. It was observed that, at the subsonic inflow conditions, there is a possibility of the occurrence of internal flow choking in dual-thrust motors due to the formation of a fluid throat at the beginning of the transition region induced by area blockage caused by boundary-layer-displacement thickness. It has been observed that a 55% increase in the upstream port area of the dual-thrust motor contributes to a 25% reduction in blockage factor at the transition region, which could negate the internal How choking and supplement with an early choking of the dual-thrust motor nozzle. If the height of the upstream port relative to the motor length is too small, the developing boundary layers from either side of the port can interact, leading to a choked,flow. On the other hand, if the developing boundary layers are far enough apart, then choking does not occur. The blockage factor is greater in magnitude for the choked case than for the unchoked case. More tangible explanations are presented in this paper for the boundary-layer blockage and the internal flow choking in dual-thrust motors, which hitherto has been unexplored.

An unsteady flow due to eccentrically rotating porous disk and a fluid at infinity

An exact solution of the unsteady Navier-Stokes equations is obtained for the flow due to non-coaxial rotations of a porous disk, executing non-torsional oscillations in its own plane, and a fluid at infinity. It is shown that the infinite number of solutions existing for a flow confined between two disks reduce to a single unique solution in the case of a single disk. The adjustment of the unsteady flow near the rotating disk to the flow at infinity rotating about a different axis is explained.

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.

Nonlinear effects on rotating disk flow in a cylindrical casing

The flow due to a finite disk rotating in an incompressible viscous fluid has been studied. A modified Newton-gradient finite difference scheme is used to obtain the solution of full Navier-Stokes equations numerically for different disk and cylinder sizes for a wide range of Reynolds numbers. The introduction of the aspect ratio and the disk-shroud gap, significantly alters the flow characteristics in the region under consideration, The frictional torque calculated from the flow data reveals that the contribution due to nonlinear terms is not negligible even at a low Reynolds number. For large Reynolds numbers, the flow structure reveals a strong boundary layer character.

Injection cooling of blunt bodies flying at high Mach numbers

A study of transpiration cooling of blunt bodies such as a hemicylinder is made by solving Navier-Stokes equations. An upwind, implicit time-marching code is developed for this purpose. The study is conducted for both perfect-gas and real-gas (chemical equilibrium) flows. Investigations are carried out for a special wall condition that is referred to as no heat flow into the wall condition. The effects of air injection on wall temperature are analyzed. Analyses are carried out for Mach numbers ranging between 6-10 and Reynolds numbers ranging between 10(6)-10(7). Studies are made for spatially constant as well as spatially varying mass injection rate distributions, White cold air injection reduces the wall temperature substantially, transpiration cooling is relatively less effective when the gas is in chemical equilibrium.

Flow induced by an asymmetrically placed disk rotating coaxially inside a cylindrical casing

In this paper, the flow due to a rotating disk non-symmetrically placed with respect to the height of the enclosing stationary cylinder is analyzed numerically. The full Navier-Stokes equations expressed in terms of stream function and vorticity are solved by successive over-relaxation for different disk radii, its distance from the bottom casing and rotational Reynolds numbers. It is observed that the flow pattern is strongly influenced by the size and the position of the disk. When the disk is very close to the top casing and small in radius, there are two regions of different scales and the vortices in the region of small scale are trapped between the disk and the top casing. Further, the variation of the moment coefficient is determined for different positions and sizes of the rotating disk. The calculations shows that the frictional torque increases rapidly, when the disk approaches the top casing. This finding is of importance for the design of vertical rotating disk reactors applied in chemical vapor deposition.

Study of transpiration cooling over a flat plate at hypersonic Mach numbers

Transpiration cooling over a flat plate at hypersonic Mach numbers is analyzed using Navier-Stokes equations, without the assumption of an isothermal wall with a prescribed wall temperature. A new criterion is proposed for determining a relevant range of blowing rates, which is useful in the parametric analysis. The wall temperature is found to decrease with the increasing blowing rate, but this effect is not uniform along the plate. The effect is more pronounced away from the leading edge. The relative change in the wall temperature is affected stronger by blowing at high Reynolds numbers. (AIAA)

Sandwich propellant combustion: Modeling and experimental comparison

This paper reports reacting fluid dynamics calculations for an ammonium percholrate binder sandwich and extracts experimentally observed features including surface profiles and maximum regression rates as a function of pressure and binder thickness. These studies have been carried out by solving the two-dimensional unsteady Navier-Stokes equations with energy and species conservation equations and a kinetic model of three reaction steps (ammonium perchlorate decomposition flame, primary diffusion flame, and final diffusion flame) in the gas phase. The unsteady two-dimensional conduction equation is solved in the condensed phase. The regressing surface is unsteady and two dimensional. Computations have been carried out for a binder thickness range of 25-125 mum and a pressure range of 1.4 to 6.9 MPa. Good comparisons at several levels of detail are used to demonstrate the need for condensed-phase two-dimensional unsteady conduction and three-step gas-phase reactions. The choice of kinetic and thermodynamic parameters is crucial to good comparison with experiments. The choice of activation energy parameters for ammonium percholrate combustion has been made with stability of combustion in addition to experimentally determined values reported in literature. The choice of gas-phase parameters for the diffusion flames are made considering that (a) primary diffusion flame affects the low-pressure behavior and (b) final diffusion flame affects high-pressure behavior. The predictions include the low-pressure deflagration limit of the sandwich apart from others noted above. Finally, this study demonstrates the possibility of making meaningful comparisons with experimental observations on sandwich propellant combustion.

Unsteady mhd flow and heat transfer on a rotating disk in an ambient fluid

An unsteady flow and heat transfer of a viscous incompressible electrically conducting fluid over a rotating infinite disk in an otherwise ambient fluid are studied. The unsteadiness in the flow field is caused by the angular velocity of the disk which varies with time. The magnetic field is applied normal to the disk surface. The new self-similar solution of the Navier-Stokes and energy equations is obtained numerically. The solution obtained here is not only the solution of the Navier-Stokes equations, but also of the boundary layer equations. Also, for a simple scaling factor, it represents the solution of the flow and heat transfer in the forward stagnation-point region of a rotating sphere or over a rotating cone. The asymptotic behaviour of the solution for a large magnetic field or for a large independent variable is also examined. The surface shear stresses in the radial and tangential directions and the surface heat transfer increase as the acceleration parameter increases. Also the surface shear stress in the radial direction and the surface heat transfer decrease with increasing magnetic field, but the surface shear stress in the tangential direction increases. (C) 2002 Editions scientifiques et medicales Elsevier SAS. All rights reserved.

A study of shock wave interaction with a rotating cylinder

Shock wave reflection over a rotating circular cylinder is numerically and experimentally investigated. It is shown that the transition from the regular reflection to the Mach reflection is promoted on the cylinder surface which rotates in the same direction of the incident shock motion, whereas it is retarded on the surface that rotates to the reverse direction. Numerical calculations solving the Navier-Stokes equations using extremely fine grids also reveal that the reflected shock transition from RRdouble right arrowMR is either advanced or retarded depending on whether or not the surface motion favors the incident shock wave. The interpretation of viscous effects on the reflected shock transition is given by the dimensional analysis and from the viewpoint of signal propagation.

Explicit filtering LES for turbulent non-premixed combustion

Two models for large eddy simulation of turbulent reacting flow in homogeneous turbulence were studied. The sub-grid stress arising out of non-linearities of the Navier-Stokes equations were modeled using an explicit filtering approach. A filtered mass density function (FMDF) approach was used for closure of the sub-grid scalar fluctuations. A posteriori calculations, when compared with the results from the direct numerical simulation, indicate that the explicit filtering is adequate in representing the effect of sub-grid stress on the filtered velocity field in the absence of reaction. Discrepancies arise when reactions occur, but the FMDF approach suffices to account for sub-grid scale fluctuations of the reacting scalars, accurately.