Unsteady nonsimilar compressible laminar two-dimensional and axisymmetric boundary-layer flows

Unsteady nonsimilar laminar compressibletwo-dimensional and axisymmetric boundarylayer flows have been studied when external velocity varies arbitrarily with time and the flow is nonhomentropic. The governing nonlinear partial differential equations with three independent variables have been solved using an implicit finite difference scheme with quasilinearization technique from the origin to the point of zero skin-friction. The results have been obtained for (i) an accelerating stream and (ii) a fluctuating stream. The skin friction responds to the fluctuations in the free stream more compared to the heat transfer. It is observed that Mach number and hot wall cause the point of zero skin friction to occur earlier whereas cold wall delays it.

Invariance group properties and exact solutions of equations describing time-dependent free surface flows under gravity

Using the method of infinitesimal transformations, a 6-parameter family of exact solutions describing nonlinear sheared flows with a free surface are found. These solutions are a hybrid between the earlier self-propagating simple wave solutions of Freeman, and decaying solutions of Sachdev. Simple wave solutions are also derived via the method of infinitesimal transformations. Incomplete beta functions seem to characterize these (nonlinear) sheared flows in the absence of critical levels.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

Exact, self-similar, time-dependent free surface flows under gravity

A class of exact, self-similar, time-dependent solutions describing free surface flows under gravity is found which extends the self-propagating class of solutions discovered earlier by Freeman (1972) to those which decay with time.

Rayleigh-testi periodisuudesta: muinaisegyptiläinen ennustuskalenteri sarjana aikapisteitä

Työssä perehdytään "päivienvalikoimistekstien" analysointiin Rayleigh'n testillä. Näissä muinaisegyptiläisissä teksteissä päiville on annettu nk. prognoosi, joka on jokaisen päivän kolmelle osalle joko "hyvä" tai "huono". Analyysin tavoitteena on selvittää, olivatko hyvät ja huonot päivät jakautuneet satunnaisesti vai periodisesti vuoden aikana. Prognooseista löytyvä 29.5 päivän periodi, jolle saadaan huomattavan korkea merkitsevyys, viittaa Kuun synodisen jakson läsnäoloon kalenterissa. Löytyvistä periodeista 7.5 ja 30 ilmenee myös muinaisegyptiläisen siviilikalenterin vaikutus ennustuksiin. Aikasarjoista löydetään myös merkitseviä periodeja, joilla on mahdollisesti numerologinen merkitys. Periodi 2.85 viittaa mahdollisesti muinaisegyptiläiseen havaintoon Algol-tähdestä, joka on yksi tähtitieteen tunnetuimmista pimennysmuuttujista. Kyseessä voi siis olla historian vanhin tunnettu muuttuvan tähden periodin määritys.

Solutions of Boundary Layer Equations Governing Free Convective and Magnetohydrodynamic Flows through Parametric Differentiation

In der vorliegenden Arbeit wird die Methode der parametrischen Differentiation angewendet, um ein System nichtlinearer Gleichungen zu lösen, das zwei- und dreidimensionale freie, konvektive Grenzschichströmungen bzw. eine zweidimensionale magnetohydrodynamische Grenzschichtströmung beherrscht. Der Hauptvorteil dieser Methode besteht darin, daß die nichlinearen Gleichungen auf lineare reduziert werden und die Nichtlinearität auf ein System von Gleichungen erster Ordnung beschränkt wird, das, verglichen mit den ursprünglichen Nichtlinearen Gleichungen, viel leichter gelöst werden kann. Ein anderer Vorzug der Methode ist, daß sie es ermöglicht, die Lösung von einer bekannten, zu einem bestimmten Parameterwert gehörigen Lösung aus durch schrittweises Vorgehen die Lösung für den gesamten Parameterbereich zu erhalten. Die mit dieser Methode gewonnenen Ergebnisse stimmen gut mit den entsprechenden, mit anderen numerischen Verfahren erzielten überein.

Drag reduction by a forward facing aerospike for a large angle blunt cone in high enthalpy flows

Drag reduction studies are conducted using a flat disc tipped aerospike for a 120-degree apex angle blunt cone model in high enthalpy flows. Accelerometer based force balance is used for the drag force measurement in the newly established free piston driven shock tunnel, HST3. Drag reduction upto about 58 percent has been achieved for Mach 8 flow of 5 MJ/kg specific enthalpy at zero degree angle of attack.

Drag reduction by counterflow supersonic jet for a blunt cone in high enthalpy flows

Counterflow supersonic jet is used as a drag reduction device during the experiments in free piston driven shock tunnel, HST3. Accelerometer based force balance is employed to measure the drag force experienced by the 60-degree apex angle blunt cone model without and with the supersonic jet opposing the hypersonic flow. It is observed that the drag force decreases with increase in injection pressure ratio until the critical injection pressure is reached. Maximum reduction in drag force of 44 percent is recorded at the critical injection pressure ratio 22.36. Further increase in injection pressure ratio has reduced the percentage drag reduction. Change in nature of the flowfield around the model has also been observed across the critical injection pressure ratio.

Boundary-layer flows with swirl and large suction

An analytical solution is presented for the laminar swirling flow in a tube. Attention is given to a particular type of swirling flow corresponding to a zero longitudinal acceleration parameter, with large suction at the surface. The investigation shows that in the case of very large rates of suction the velocity overshoot can be almost eliminated. This is even possible in flows with swirls which are characterized by a velocity overshoot in the longitudinal direction.

A computer code for enclosed natural convection in the boundary layer regime

A computer code is developed for the numerical prediction of natural convection in rectangular two-dimensional cavities at high Rayleigh numbers. The governing equations are retained in the primitive variable form. The numerical method is based on finite differences and an ADI scheme. Convective terms may be approximated with either central or hybrid differencing for greater stability. A non-uniform grid distribution is possible for greater efficiency. The pressure is dealt with via a SIMPLE type algorithm and the use of a fast elliptic solver for the solenoidal velocity correction field significantly reduces computing times. Preliminary results indicate that the code is reasonably accurate, robust and fast compared with existing benchmarks and finite difference based codes, particularly at high Rayleigh numbers. Extension to three-dimensional problems and turbulence studies in similar geometries is readily possible and indicated.

The laminar boundary layer regime for natural convection of air in a square cavity

Numerical predictions are obtained for laminar natural convection of air in a square two dimensional cavity at high Rayleigh numbers. Proper resolution of the core reveals weak multi-cellular structure which varies in a complex manner as the effects of convection are increased. The end of the steady laminar regime is numerically estimated to occur at Ra=2.2x10^8.

Laminar and turbulent natural convection in rectangular cavities

A computer program has been developed for the prediction of buoyancy-driven laminar and turbulent flow in rectangular air-filled two-dimensional cavities with differentially heated side walls. Laminar flow predictions for a square cavity and Rayleigh numbers from Ra = 10^3 up to the onset of unsteady flow have been obtained. Accurate solutions for Ra = 5 x 10^6, 10^7, 5 x 10^7 and 10^8 are presented and an estimate for the critical Rayleigh number at which the steady laminar flow becomes unsteady is given for this geometry. Numerical predictions of turbulent flow have been obtained for RaH~0(10^9 -10^11 ) and compared with existing experimental data. A previously developed second moment closure model (Behnia et al. 1987) has been used to model the turbulence. Results indicate that a second moment closure model is capable of predicting the observed flow features.