The three dimensional incompressible laminar boundary layer on a spinning cone

Summary The motion of an incompressible viscous fluid induced by a spinning cone is analytically studied and similar solutions of the relevant steady state boundary equations are obtained. Some of the numerical results are shown to be obtainable from the Karman-Cochran solution for the infinite disc.Symbols and Notation p Pressure - p Pressure at infinity - p ₀ Pressure at the wall - Density - Transverse component of velocity - Normal component of velocity - Radial component of velocity - Angular velocity - Semi-vertex angle - Re Reynolds number with respect to _o - _o Transverse component of velocity at the cone surface - Kinematic viscosity This research is sponsored by the Air Force Office of Scientific Research, Fluid Mechanics Division, under Contract Number AF 18(600)-498.     

Flow in vicinity of blunt body stagnation point in diverging hypersonic stream

In the hypersonic thin shock layer approximation for a small ratio k of the densities before and after the normal shock wave the solution of [1] for the vicinity of the stagnation point of a smooth blunt body is extended to the case of nonuniform outer flow. It is shown that the effect of this nonuniformity can be taken into account with the aid of the effective shock wave radius of curvature R_*, whose introduction makes it possible to reduce to universal relations the data for different nonuniform outer flows with practically the same similarity criterion k. The results of the study are compared with numerical calculations of highly underexpanded jet flow past a sphere.Notations x, y a curvilinear coordinate system with axes directed respectively along and normal to the body surface with origin at the forward stagnation point - R radius of curvature of the meridional plane of the body surface - uV, vV., , p V ² respectively the velocity projections on the x, y axes, density, and pressure - and V freestream density and velocity The indices =0 and=1 apply to plane and axisymmetric flows Izv. AN SSSR, Mekhanika Zhidkosti i Gaza, Vol. 5, No. 3, pp. 102–105, 1970.     

Robustness of Exponential Dichotomies in Infinite-Dimensional Dynamical Systems

In this paper we examine the issue of the robustness, or stability, of an exponential dichotomy, or an exponential trichotomy, in a dynamical system on an Banach space W. These two hyperbolic structures describe long-time dynamical properties of the associated time-varying linearized equation _t +A=B(t) , where the linear operator B(t) is the evaluation of a suitable Fréchet derivative along a given solution in the set K in W. Our main objective is to show, under reasonable conditions, that if B(t)=B(, t) depends continuously on a parameter and there is an exponential dichotomy, or exponential trichotomy, at a value ₀, then there is an exponential dichotomy, or exponential trichotomy, for all near ₀.We present several illustrations indicating the significance of this robustness property.     

A direct formulation of correlations for wall-proximity corrections of hot-wire measurements

Correlations for corrections to hot-wire data for the effects of wall proximity within the viscous sublayer are usually presented in the form u/u = F (y u /). The application of such correlations requires a prior knowledge of the wall shear stress; alternatively, the correlation must be used in an iterative fashion. It is shown in the present note that any such correlation may be recast with no loss of generality in the explicit form u/u _m = f (y u _m/), which is more convenient for use.List of symbols u difference between measured and actual velocities, u _m – u - u _m measured velocity - u shear velocity, - u ⁺ on-dimensional velocity, u/u - y distance from wall - y ⁺ non-dimensional distance from wall, y u / - fluid density - fluid kinematic viscosity - _s wall shear stress     

Thermal radiation effects on MHD-Rayleigh flow with constant surface heat flux

Thermal radiation heat transfer effects on the Rayleigh flow of gray viscous fluids under the effect of a transverse magnetic field are investigated. The free convection heat transfer problem from constant surface heat flux moving plate is selected for study. It is found that the increasing of the magnetic field number M= H₀² / U₀²decreased velocities inside boundary layer, the increasing of the conduction–radiation parameter Rd=k__R/4aT³ decreased both temperatures and heat transfer rates. It is also found that the increasing of the dimensionless surface heat flux parameter q₀^*=q₀ /(kU₀T₎ increased the temperatures inside the boundary layer and increased the heat transfer rates. Comparison with previous works shows excellent agreement. Different transient velocity profiles, temperature profiles and local Nusselt numbers against different dimensionless groups are drawn.     

Bouncing and Coalescence of Bubble Pairs Rising at High Reynolds Number in Pure Water or Aqueous Surfactant Solutions

The encounter of bubble pairs of O(1 mm) in both pure water and aqueous surfactant solutions was studied experimentally. In pure water, two equally sized bubbles were found to coalesce if the Weber number, W = V² R/, based on the velocity of approach, V, was below a critical value, W_cr = 0.18, where and are the density and surface tension of the liquid respectively and R the equivalent radius of the bubbles. After coalescence bubbles perform volume and shape oscillations.When W_cr is exceeded, bubbles bounce. After bouncing, bubbles can either coalesce or separate without coalescing. This was found to depend on the Weber number, based on the rise velocity U, We = U² R/. If this number was below a critical value, bubbles coalesced after bouncing. The relative motion of the bubbles was found to be damped out by acoustic damping due to surface oscillations rather then by viscosity.If We was above a critical value, which was close to that for path instability of a single bubble (We = 3.3), the bubbles separated after bouncing. This is probably caused by shedding of vortices which dominate the relative motion of the bubbles. This mechanism may cause bubbles in bubbly flows not aggregating in horizontal planes, as was found in calculations based on potential flow theory. For modelling bubbly flows it will therefore be essential to incorporate the influence of vorticity.When surfactants are added to the water it was found that bubbles are prevented to coalesce above a critical concentration, which is nearly identical to that of single rising bubbles. Above this critical concentration, bubbles behave as rigid spheres and trajectories cannot be predicted by potential flow theory.     

Transformation and some solutions of the equations of motion of an ideal gas

The equations of one-dimensional and plane steady adiabatic motion of an ideal gas are transformed to a new form in which the role of the independent variables are played by the stream function and the function introduced by Martin [1, 2], It is shown that the function retains a constant value on a strong shock wave (and on a strong shock for plane flows). For one-dimensional isentropic motions the resulting transformation permits new exact solutions to be obtained from the exact solutions of the equations of motion. It is shown also that the one-dimensional motions of an ideal gas with the equation of state p=f(t) and the one-dimensional adiabatic motions of a gas for which p=f() are equivalent (t is time, is the stream function). It is shown that if k=s=–1, m and n are arbitrary (m+n0) and =1, the general solution of the system of equations which is fundamental in the theory of one-dimensional adiabatic self-similar motions [3] is found in parametric form with the aid of quadratures. Plane adiabatic motions of an ideal gas having the property that the pressure depends only on a single geometric coordinate are studied.     

Momentum balance in highly distorted turbulent pipe flows

An in depth study into the development and decay of distorted turbulent pipe flows in incompressible flow has yielded a vast quantity of experimental data covering a wide range of initial conditions. Sufficient detail on the development of both mean flow and turbulence structure in these flows has been obtained to allow an implied radial static pressure distribution to be calculated. The static pressure distributions determined compare well both qualitatively and quantitatively with earlier experimental work. A strong correlation between static pressure coefficient Cp and axial turbulence intensity is demonstrated.List of symbols C _p static pressure coefficient = (pw-p)/1/2 - D pipe diameter - K turbulent kinetic energy - (r, , z) cylindrical polar co-ordinates. / 0 - R, y pipe radius, distance measured from the pipe wall - U, V axial and radial time mean velocity components - mean value of u - u, u/ , / - u, , w fluctuating velocity components - axial, radial turbulence intensity - turbulent shear stress - u friction velocity, (u ² = ₀/p) - ₀ wall shear stress - ^* boundary layer thickness A version of this paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

On Perturbative Expansions to the Stochastic Flow Problem

When analyzing stochastic steady flow, the hydraulic conductivity naturally appears logarithmically. Often the log conductivity is represented as the sum of an average plus a stochastic fluctuation. To make the problem tractable, the log conductivity fluctuation, f, about the mean log conductivity, lnK _G, is assumed to have finite variance, _f ². Historically, perturbation schemes have involved the assumption that _f ²<1. Here it is shown that _f may not be the most judicious choice of perturbation parameters for steady flow. Instead, we posit that the variance of the gradient of the conductivity fluctuation, _f ², is more appropriate hoice. By solving the problem withthis parameter and studying the solution, this conjecture can be refined and an even more appropriate perturbation parameter, , defined. Since the processes f and f can often be considered independent, further assumptions on f are necessary. In particular, when the two point correlation function for the conductivity is assumed to be exponential or Gaussian, it is possible to estimate the magnitude of f in terms of _f and various length scales. The ratio of the integral scale in the main direction of flow ( _x ) to the total domain length (L^*), _x ²=_x/L^*, plays an important role in the convergence of the perturbation scheme. For _x smaller than a critical value _c, _x < _c, the scheme's perturbation parameter is =_f/_x for one- dimensional flow, and =_f/_x ² for two-dimensional flow with mean flow in the x direction. For _x > _c, the parameter =_f/_x ³ may be thought as the perturbation parameter for two-dimensional flow. The shape of the log conductivity fluctuation two point correlation function, and boundary conditions influence the convergence of the perturbation scheme.     

The chaotic response of a fluttering panel: The influence of maneuvering

The influence of maneuvering on the chaotic response of a fluttering buckled plate on an aircraft has been studied. The governing equations, derived using Lagrangian mechanics, include geometric non-linearities associated with the occurrence of tensile stresses, as well as coupling between the angular velocity of the maneuver and the elastic degrees of freedom. Numerical simulation for periodic and chaotic responses are conducted in order to analyze the influence of the pull-up maneuver on the dynamic behavior of the panel. Long-time histories phase-plane plots, and power spectra of the responses are presented. As the maneuver (load factor) increases, the system exhibits complicated dynamic behavior including a direct and inverse cascade of subharmonic bifurcations, intermittency, and chaos. Beside these classical routes of transition from a periodic state to chaos, our calculations suggest amplitude modulation as a possible new mode of transition to chaos. Consequently this research contributes to the understanding of the mechanisms through which the transition between periodic and strange attractors occurs in, dissipative mechanical systems. In the case of a prescribed time dependent maneuver, a remarkable transition between the different types of limit cycles is presented.Nomenclature a plate length - a _r u _r /h - D plate bending stiffness - E modulus of elasticity - g acceleration due to gravity - h plate thickness - j₁,j₂,j₃ base vectors of the body frame of reference - K spring constant - M Mach number - n 1 + ₀/g - N ₁ applied in-plane force - p–p aerodynamic pressure - P pa ⁴/Dh - q ₀/2 - Q _r generalized Lagrangian forces - R rotation matrix - R ₄ N, a ²/D - t time - kinetic energy - u plate deflection - u displacement of the structure - u _r modal amplitude - v₀ velocity - x coordinates in the inertial frame of reference - z coordinates in the body frame of reference - Ka/(Ka+Eh) - - elastic energy - 2qa ³/D - a/_mh - Poisson's ratio - material coordinates - air density - _m plate density - - _r prescribed functions - _r sin(r z/a) - angular velocity - a/v₀ - skew-symmetric matrix form of the angular velocity     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

Magnetohydrodynamic natural convection heat transfer from radiate vertical porous surfaces

Magnetohydrodynamic natural convection heat transfer from radiate vertical surfaces with fluid suction or injection is considered. The nonsimilarity parameter is found to be the conductive fluid injection or suction along the streamwise coordinate = V{4x/² g(T _w – T )}^1/4. Three dimensionless parameters had been found to describe the problem: the magnetic influence number N = B ² _y /V ², the radiation-conduction parameter R _d = k _R /4aT ³ , and the Gebhart number Ge _x = gx/cp to represent the effect of the viscous dissipation. It is found that increasing the magnetic field strength causes the velocity and the heat transfer rates inside the boundary layer to decrease. Its apparent that increasing the radiation-conduction parameter decreases the velocity and enhances the heat transfer rates. The Gebhart number, i.e, the viscous dissipation had no effect on the present problem.Nomenclature a Stefan-Boltzmann constant - B _y Magnetic field flux density Wb/m² - Cf _x Local skin friction factor - c _p Specific heat capacity - f Dimensionless stream function - Ge _x Gebhart number, gx/cp - g Gravitational acceleration - k Thermal Conductivity - L Length of the plate - N Magnetic influence number, B ² _y /V ² - p Pressure - Pr Prandtl number - q _r Radiative heat flux - q _w (x) Local surface heat flux - Q _w (x) Dimensionless Local surface heat flux - R _d Planck number (Radiation-Conduction parameter), k _R /4aT ³ - T Temperature - T Free stream temperature - T _w Wall temperature - u, v Velocity components in x- and y-directions - V Porous wall suction or injection velocity - V _w Porous wall suction or injection velocity - x, y Axial and normal coordinates - Thermal diffusivity Greek symbols _R Roseland mean absorption coefficient, 4/3R _d - Coefficient of thermal expansion - Nonsimilarity parameter, V{4x/² g(T _w – T )}^1/4 - Peseudo-similarity variable - Dimensionless temperature - _w Ratio of surface temperature to the ambient temperature, T _w /T - Dynamice viscosity - Kinemtic viscosity - Fluid density - Electrical conductivity - _w Local wall shear stress - Dimensional stream function     

An experimental study of the lateral migration of a droplet in a creeping flow

The distribution of droplets in a plane Hagen-Poiseuille flow of dilute suspensions has been measured by a special LDA technique. This method assumes a well defined relation between the velocity of the droplets and their lateral position in the channel. The measurements have shown that the droplet distribution is non-uniform and depends on the viscosity ratio between the droplets and the carrier liquid. The results have been compared with a theory by Chan and Leal describing the lateral migration of suspended droplets.List of symbols a particle radius, m - d half width of the channel, m - Re flow Reynolds number, = 2 _m · d · /µ - flow velocity, m/s - _m flow velocity at the channel axis, m/s - We Weber number, = 2 _m ^Emphasis>/2 · d · / - x distance from center line (x = 0) of the channel, m - non-dimensional distance from the channel center line, _x ^d - y distance along the channel (y = 0 at channel inlet), m - non-dimensional distance along the channel, = y/2d - non-dimensional, normalized distance along the channel, = · _m · µ/ - interfacial tension, N/m - viscosity ratio of dispersed (droplet) phase to viscosity of continuous phase - µ viscosity of continuous phase, Pa · s - density of continuous phase, kg/m³ - phase density difference, kg/m³ Experiments were performed at Max-Planck-Institut, Göttingen     

On the forced motion of elastic solids

Summary This note develops a method for the solution of the elastokinetic boundary value problem for time dependent surface tractions and/or displacements, as well as body forces which are functions of time and space. The method of Williams is extended to resolve three-dimensional problems of elastodynamics by classical mathematical techniques.Nomenclature x _i position vector - t time - u _i displacement vector - _ij stress tensor - F _i vector characterizing body force per unit volume - stress vector acting on surface S with unit outer normal v _i - density - , Lamé's constants - _ij Kronecker delta     

Analytic Approximation of the Homoclinic Orbits of the Lorenz System at σ = 10, b = 8/3 and ρ = 13.926...

We present an iterative technique to analytically approximate the homoclinic loops of the Lorenz system for = 10, b = 8/3 and = _H = 13.926.... First, the local structure of the homoclinic solution for t 0 ± and t ± is analyzed. Then, global approximants are used to match the local expansions. The matching procedure resembles the one used in Padé approximations. The accuracy of the approximation is improved iteratively, with each iteration providing estimates for the initial conditions of the homoclinic orbit, the value of _H, and three undetermined constants in the local expansions. Within three iterations the error in _H falls to the order of 0.1%. Comparisons with numerical integrations are made, and a discussion on ways to extend the technique to other types of homoclinic or heteroclinic orbits, and to improve its accuracy, is given.     

Analyse und Berechnung der Strömungsvorgänge im zylindrischen Scherspalt von Plastifizieraggregaten

Zusammenfassung Dieser Aufsatz zeigt eine Möglichkeit auf, zylindrische Scherteile einer Plastifiziereinheit, auf der strukturviskose Materialien verarbeitet werden, approximativ zu berechnen. Es ist möglich, den Volumenstrom und Druckabfall, die mittlere Schergeschwindigkeit, Scherdeformation und Schubspannung im Scherspalt zu approximieren. Durch diese Gleichungen wird eine Abschätzung der Verteil- und Zerteilvorgänge im Scherelement möglich.

---

A method is described for approximating the flow in cylindrical shearing gaps of plasticating extruder, which is applicable to shear thinning materials. It is possible to calculate the through-put and pressure drop as well as the shear rate, strain and shear stress in the gap. With these equations the distribution and separation process in shearing gaps can be evaluated.

---

D Zylinderdurchmesser - d ₁ Schnecken-Kerndurchmesser der Meteringzone - d _s Durchmesser des zylindrischen Scherteils - K Konstante im Potenzfließgesetz - K _0T Koeffizient des Potenzfließgesetzes - L ₁ Länge der Anlaufschräge - L _s Länge des zylindrischen Scherteils - n Fließindex - n ₀ Drehzahl - p Druckabfall über der Scherteillänge - s Scherspalthöhe - T _M Massetemperatur - ₀ Umfangsgeschwindigkeit - _0x Geschwindigkeitskomponente inx-Richtung - _x, _z Geschwindigkeit inx- bzw.z-Richtung als Funktion der Koordinatey - Volumenstrom - x, z Ortskoordinaten - Exponent des Potenzfließgesetzes - Schergeschwindigkeit - mittlere Schergeschwindigkeit - Viskosität - dimensionslose Höhe - Dichte der Schmelze - Schubspannung - _yx, _yz Schubspannungskomponenten - _xx, _zz Normalspannungskomponenten - _ps dimensionsloser Druckgradient - dimensionsloser Volumenstrom - _x, _z dimensionslose Geschwindigkeit inx- bzw.z-Richtung     

On Full Two-Scale Expansion of the Solutions of Nonlinear Periodic Rapidly Oscillating Problems and Higher-Order Homogenised Variational Problems

We consider a scalar quasilinear equation in the divergence form with periodic rapid oscillations, which may be a model of, e.g., nonlinear conducting, dielectric, or deforming in a restricted way hardening elastic-plastic composites, with outer periodicity conditions of a fixed large period. Under some natural growth assumptions on the stored-energy function, we construct for uniformly elliptic problems a full two-scale asymptotic expansion, which has a precise double-series structure, separating the slow and the fast variables in all orders, so that its slowly varying part solves asymptotically an infinite-order homogenised equation (cf. Bakhvalov, N.S., Panasenko, G.P.: Homogenisation: Averaging Processes in Periodic Media. Nauka, Moscow, 1984 (in Russian); English translation: Kluwer, 1989), and whose higher-order terms depend on the higher gradients of the slowly varying part. We prove the error bound, i.e., that the truncated asymptotic expansion is higher-order close to the actual solution in appropriate norms. The approach is extended to a non-uniformly elliptic case: for two-dimensional power-law potentials we prove the non-degeneracy using topological index methods. Examples and explicit formulae for the higher-order terms are given. In particular, we prove that the first term in the higher-order homogenised equations is related to the first-order corrector to the mean flux, and has in general the form of a fully nonlinear operator which is quadratic with respect to its highest (second) derivative being a linear combination of the second minors of the Hessian with coefficients depending on the first gradient, and in dimension two is of Monge-Ampère type. We show that this term is present at least for some examples (three-phase power-law laminates).In the second part of the paper we extend to this nonlinear context some of the results previously developed by us in the linear case (Smyshlyaev, V.P., Cherednichenko, K.D. J. Mech. Phys. Solids 48, 1325–1357, 2000). In particular, we prove that the slowly varying part of the full asymptotic expansion is the rigorous asymptotics in all orders for the translationally averaged actual solution and flux, or in the sense of a higher-order version of the weak convergence. We then explore to what extent the method of variational truncation of the infinite-order homogenised equation, successfully implemented by us in the linear context in the previous work for constructing explicit higher-order homogenised equations, is extendable to the nonlinear regime. We propose a natural extension and prove that at least under some further natural non-degeneracy assumptions it has a solution (the existence), and that any such solution is close to the actual solution in appropriate norms.Acknowledgement We acknowledge the support of the Department of Mathematical Sciences, University of Bath, U.K. and partial support of the Isaac Newton Institute for Mathematical Sciences, Cambridge, U.K. KDC also acknowledges the support of St. Johns College, Oxford, U.K. and thanks the Mathematical Institute, Oxford for its hospitality.     

Multicomponent turbulent boundary layer on a chemically active surface

When a gaseous mixture flows past chemically active surfaces the boundary layer formed on the wetted body may contain a large number of components with different diffusion properties. This leads to the necessity for studying the diffusion of the components in the multicomponent boundary layer.The use of thebinary boundary layer concept in the general case cannot yield satisfactory results, since replacement of the mutual diffusion coefficients D_ij of the various pairs of components by a single diffusion coefficient D in many cases is a rough approximation.In the general case the number of different diffusion coefficients is equal to N(N–1)/2 (N is the number of components). Usually it is possible to identify groups of components with similar molecular weights. Then the number of different diffusion coefficients may be reduced without large error. However, even in the comparatively simple case when it is possible to divide all the components into two groups with similar molecular weights we must take account of three different diffusion coefficients (one diffusion coefficient in each group and also the diffusion coefficient for the components of one group relative to the components of the other group). Only in particular cases when the gaseous mixture consists of only two components with arbitrary molecular weights, or if all the components of the gaseous mixture have similar molecular weights, can we with justification introduce a single diffusion coefficient (if in this case there are no limitations on the direction of the diffusion).Studies have been published covering the laminar multicomponent boundary layer. An analytic method for solving the equations of the laminar multicomponent boundary layer was developed by Tirskii [1]. There are also studies in which concrete results were obtained by numerical methods with the use of computers (for example, [2, 3]).As far as the author knows, for turbulent flow there are studies (for example, [4, 5]) covering flow with chemical reactions only in the case when all the diffusion coefficients are equal (D_ij=D).The present paper presents a method for calculating the turbulent multicomponent boundary layer with account for several different diffusion coefficients.Notation x, y coordinates - u, v velocity components - density - T temperature - h heat content - H enthalpy - c_i mass concentration of the i-th component - c ₁ ⁽¹⁾ element concentrations in solid body - J_i diffusion flux of the i-th component - m molecular weight - dynamic viscosity coefficient - kinematic viscosity coefficient - heat conduction coefficient - c_p specific heat - adiabatic index - D_ij binary diffusion coefficients - P Prandtl number - S_ij Schmidt number - S_t Stanton number - M Mach number - friction - q radiant thermal flux - boundary layer thickness - D rate of displacement of gas-solid interface - degree of gasification - r_ij weight fraction of element i in component j - _ij stoichiometric coefficients - K_i reaction equilibrium constants - l number of components for which I_i0 Indices i, j component number - w quantities for y=0 - * quantities on the edge of the laminar sublayer - ⁽¹⁾ quantities at the solid body - quantities at the outer edge of the boundary layer - molar transport coefficients     

Heat transfer and equilibrium temperature of a sphere in a supersonic rarefied gas flow

Some results are presented of experimental studies of the equilibrium temperature and heat transfer of a sphere in a supersonic rarefied air flow.The notations D sphere diameter - u, , T,,l, freestream parameters (u is velocity, density, T the thermodynamic temperature,l the molecular mean free path, the viscosity coefficient, the thermal conductivity) - T₀ temperature of the adiabatically stagnated stream - T_e mean equilibrium temperature of the sphere - T_w surface temperature of the cold sphere (T_we) - mean heat transfer coefficient - _e air thermal conductivity at the temperature T_e - P Prandtl number - M Mach number     

Laminar assisting mixed convection heat transfer from a vertical isothermal plate to water with variable physical properties

The steady laminar boundary layer flow, with an external force, along a vertical isothermal plate is studied in this paper. The external force may be produced either by the motion of the plate or by a free stream. The fluid is water whose density-temperature relationship is non-linear at low temperatures and viscosity and thermal conductivity are functions of temperature. The results are obtained with the numerical solution of the boundary layer equations with , k and variable across the boundary layer. Both upward and downward flow is considered. It was found that the variation of , k and with temperature has a strong influence on mixed convection characteristics.Nomenclature c_p water specific heat - f dimensionless stream function - g gravitational acceleration - Gr_x local Grashof number - k thermal conductivity - Nu_x local Nusselt number - Pr Prandtl number - Pr_a ambient Prandtl number - Re_x local Reynolds number - s salinity - T water temperature - T_a ambient water temperature - T_o plate temperature - u vertical velocity - u_a free stream velocity - u_o plate velocity - v horizontal velocity - x vertical coordinate - y horizontal coordinate - pseudo-similarity variable - nondimensional temperature - dynamic viscosity - _f film dynamic viscosity - _o dynamic viscosity at plate surface - kinematic viscosity - buoyancy parameter - water density - _a ambient water density - _f film water density - _o water density at plate surface - physical stream function