Hydrodynamic models as applied to the investigation of magnetic plasma stability

In the present paper magnetohydrodynamic models are employed to investigate the stability of an inhomogeneous magnetic plasma with respect to perturbations in which the electric field may be regarded as a potential field (rot E 0). A hydrodynamic model, actually an extension of the well-known Chew-Goldberg er-Low model [1], is used to investigate motions transverse to a strong magnetic field in a collisionless plasma. The total viscous stress tensor is given; this includes, together with magnetic viscosity, the so-called inertial viscosity.Ordinary two-fluid hydrodynamics is used in the case of strong collisions=. It is shown that the collisional viscosity leads to flute-type instability in the case when, collisions being neglected, the flute mode is stabilized by a finite Larmor radius. A treatment is also given of the case when epithermal high-frequency oscillations (not leading immediately to anomalous diffusion) cause instability in the low-frequency (drift) oscillations in a manner similar to the collisional electron viscosity, leading to anomalous diffusion.Notation f particle distribution function - E electric field component - H₀ magnetic field - density - V particle velocity - e charge - m, M electron and ion mass - _i, _e ion and electron cyclotron frequencies - viscous stress tensor - P pressure - r_i Larmor radius - P pressure tensor - t time - frequency - T temperature - collision frequency - collision time - j current density - _i, _e ion and electron drift frequencies - k_x, k_y, k_z wave-vector components - n₀ particle density - g acceleration due to gravity. The authors are grateful to A. A. Galeev for valuable discussion.     

Hypersonic flow past blunt-edged delta wings

A method is proposed for calculating hypersonic ideal-gas flow past blunt-edged delta wings with aspect ratios = 100–200. Systematic wing flow calculations are carried out on the intervals 6 M 20, 0 20, 60 80; the results are analyzed in terms of hypersonic similarity parameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 175–179, September–October, 1990.     

The mirror relations and nonlinear viscoelasticity of polymer melts

An attempt is made to incorporate into a quasilinear viscoelastic constitutive equation of the Boltzmann superposition type the two mirror relations of Gleissle, as well as his relation between the steady-state first normal-stress difference and the shear viscosity curve. It is shown that the three relations can hold separately within this constitutive model, but not simultaneously, because they require a different nonlinear strain measure, namelyS ₁₂ () = – a ( – 1) (a = 0 for 1,a = 1 for 1) for the mirroring of the viscosities,S ₁₂ () = – a (–k ²/) (a = 0 for k, a = 1 for k) for the mirroring of the first normal-stress coefficients, and for the third relation. Here denotes the shear strain and erf the error function. Experimental data on melts of a low-density polyethylene, a high-density polyethylene and a polypropylene show that the mirror relations are passable approximations, but that the third relation meets reality surprisingly close if the right value ofk is used.     

Instability of a plane axisymmetric flow with a shear discontinuity in the velocity profile

It is proposed to investigate the stability of a plane axisymmetric flow with an angular velocity profile (r) such that the angular velocity is constant when r < r_O – L and r > r_O + L but varies monotonically from ₁ to ₂ near the point r_O, the thickness of the transition zone being small L r_O, whereas the change in velocity is not small ¦₂ – ₁¦ ₂, ₁. Obviously, as L O short-wave disturbances with respect to the azimuthal coordinate (k=m/r_O 1/r_O) will be unstable with a growth rate-close to the Kelvin—Helmholtz growth rate. In the case L=O (i.e., for a profile with a shear-discontinuity) we find the instability growth rate _O and show that where the thickness of the discontinuity L is finite (but small) the growth rate does not differ from _O up to terms proportional to kL 1 and 1/m 1. Using this example it is possible to investigate the effect of rotation on the flow stability. It is important to note that stabilization (or destabilization) of the flow in question by rotation occurs only for three-dimensional or axisymmetric perturbations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–114, January–February, 1985.     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

The onset of Taylor-like vortices in the flow induced by an impulsively started rotating cylinder

The onset of instability in the flow by an impulsively started rotating cylinder is analyzed under linear theory. It is well-known that at the critical Taylor number T_c=1695 the secondary flow in form of Taylor vortices sets in under the narrow-gap approximation. Here the dimensionless critical time _c to mark the onset of instability for TT_c is presented as a function of the Taylor number T. Available experimental data of water indicate that deviation of the velocity profiles from the primary flow occurs starting from a certain time 4_c. It seems evident that during _c4_c the secondary flow is very weak and the primary state of time-dependent annular Couette flow is maintained.     

On the Temperature of a Rarefied Gas in Non‐Equilibrium

This paper addresses the problem of the proper definition of temperature of a gas in nonequilibrium. It shows that the mean kinetic energy of the atoms of a rarefied gas is not a good measure for thethermodynamic temperature, because in general it jumps at a wall, and because it is nonmonotone in a onedimensional process of stationary heat conduction. The jump of the kinetic temperature is calculated and found to be about 5K in a rarefied gas. The basis for the calculations is provided by the arguments of extended thermodynamics of 14 moments. An essential tool is the minimax principle of entropy production recently postulated by Struchtrup Weiss [1], because it furnishes one important boundary condition.Sommario. Il lavoro riguarda la corretta definizione della temperatura di un gas in condizioni di nonequilibrio. Si mostra come lenergia cinetica media degli atomi di un gas rarefatto non sia una buona misura della temperatura termodinamica poiché in generale, essa risulta discontinua su una parete e nonmonotona in un processo unidimensionale di conduzione stazionaria del calore. Viene calcolato il salto della temperatura cinetica che risulta pari a circa 5K in un gas rarefatto. La base per il calcolo è fornita dal contesto della termodinamica estesa di 14 momenti. Uno strumento essenziale è rappresentato dal principio di minimax di produzione di entropia recentemente postulato da Struchtrup and Weiss [1], che fornisce unimportante condizione alcontorno.     

Shear Instability at the “Explosion Product–Metal” Interface for Sliding Detonation of an Explosive Charge

Periodic perturbations at the explosion product–metal interface were studied experimentally. Experiments were performed for both spherical and plane geometry. Critical conditions of wave formation (detonation velocity of an explosive charge D 6.9 mm/sec) are determined, and an explanation of this effect is given. It is found experimentally that a dynamic pulse causes intense plastic strains at the explosion products–metal interface, leading to thermal softening of the steel boundary layer. In this layer, Kelvin–Helmholtz instability occurs. Calculationanalytical estimates of the critical boundary unstable wavelength agree satisfactorily with experimental results.     

Calculation of heat transfer accompanying hypersonic three-dimensional flow of air over slender blunt-nosed cones

The effective length method [1, 2] has been used to make systematic calculations of the heat transfer for laminar and turbulent boundary layers on slender blunt-nosed cones at small angles of attack ( + 5° in a separationless hypersonic air stream dissociating in equilibrium (half-angles of the cones 0 20°, angles of attack 0 15°, Mach numbers 5 M 25). The parameters of the gas at the outer edge of the boundary layer were taken equal to the inviscid parameters on the surface of the cones. Analysis of the results leads to simple approximate dependences for the heat transfer coefficients.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1981.     

Some problems concerning nonstationary currents in dense plasma systems confined by walls

Nonstationary currents are examined in a dense magnetized plasma with 1, in which energy release and heat loss by thermal conduction and radiation are possible. Solutions are found in two limiting cases: ¦f¦ ¦ div (T)¦ and ¦f¦ ¦ div(T)¦ (f is the radiation intensity, is the coefficient of heat conduction, and T is the temperature). In the first case a solution was obtained of some problems of the cooling and heating of a plasma illustrated in part by the evolution in time of the temperature profile in the boundary layer. In the second case an isomorphic solution was found for an arbitrary dependence of the coefficient of heat conduction on the temperature, pressure, and magnetic field.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 3–8, January–February, 1972.The author is grateful to G. I. Budker for formulating the problem.     

Unsteady Radiative–Convective Heat Transfer in a Flow Emitting‐Absorbing and Scattering Medium around an Ablating Plate

A conjugate problem of radiative–convective heat transfer in a turbulent hightemperature gasdisperse flow around a thermally thin ablating plate is considered. The plate experiences intense radiative heating by an external source, which is a blackbody. The temperature fields and the distributions of heat fluxes along the plate under unsteady conditions are calculated. The data gained make it possible to examine the effect of the Stark number and phasetransition heat in the plate material on the time evolution of the thermal state of the boundarylayer medium and the plate itself being heated by a hightemperature radiation source.     

Asymptotic analysis of equilibrium states for rotating turbulent flows

The equilibrium states of homogeneous turbulence simultaneously subjected to a mean velocity gradient and a rotation are examined by using asymptotic analysis. The present work is concerned with the asymptotic behavior of quantities such as the turbulent kinetic energy and its dissipation rate associated with the fixed point (/kS)=0, whereS is the shear rate. The classical form of the model transport equation for (Hanjalic and Launder, 1972) is used. The present analysis shows that, asymptotically, the turbulent kinetic energy (a) undergoes a power-law decay with time for (P/)<1, (b) is independent of time for (P/)=1, (c) undergoes a power-law growth with time for 1<(P/)<(C ₂–1), and (d) is represented by an exponential law versus time for (P/)=(C ₂–1)/(C ₁–1) and (/kS)>0 whereP is the production rate. For the commonly used second-order models the equilibrium solutions forP/,II, andIII (whereII andIII are respectively the second and third invariants of the anisotropy tensor) depend on the rotation number when (P/kS)=(/kS)=0. The variation of (P/kS) andII versusR given by the second-order model of Yakhot and Orzag are compared with results of Rapid Distortion Theory corrected for decay (Townsend, 1970).     

Investigation of the flow past wings of infinite span with a blunt leading edge in the vorticity interaction regime

An asymptotic analysis of the Navier-Stokes equations is carried out for the case of hypersonic flow past wings of infinite span with a blunt leading edge when 0, Re , and M . Analytic solutions are obtained for an inviscid shock layer and inviscid boundary layer. The results of a numerical solution of the problems of vorticity interaction at the blunt edge and on the lateral surface of the wing are presented. These solutions are compared with the solution of the equations of a thin viscous shock layer and on the basis of this comparison the boundaries of the asymptotic regions are estimated.deceasedTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 120–127, November–December, 1987.     

Numerical analysis of axisymmetric buckling of conical shells

Nonlinear boundaryvalue problems of axisymmetric buckling of conical shells under a uniformly distributed normal pressure are solved by the shooting method. The problems are formulated for a system of six firstorder ordinary differential equations with independent rotation and displacement fields. Simply supported and clamped cases are considered. Branching solutions of the boundaryvalue problems are studied for different pressures and geometrical parameters of the shells. The nonmonotonic and discontinuous curves of equilibrium states obtained show that collapse, i.e., snapthrough instability is possible. For a simply supported shell, multivalued solutions are obtained for both external and internal pressure. For a clamped thinwalled shell, theoretical results are compared with experimental data.     

Mixed convection flow over a vertical plate with localized heating (cooling), magnetic field and suction (injection)

An analysis is carried out to study the effects of localized heating (cooling), suction (injection), buoyancy forces and magnetic field for the mixed convection flow on a heated vertical plate. The localized heating or cooling introduces a finite discontinuity in the mathematical formulation of the problem and increases its complexity. In order to overcome this difficulty, a non-uniform distribution of wall temperature is taken at finite sections of the plate. The nonlinear coupled parabolic partial differential equations governing the flow have been solved by using an implicit finite-difference scheme. The effect of the localized heating or cooling is found to be very significant on the heat transfer, but its effect on the skin friction is comparatively small. The buoyancy, magnetic and suction parameters increase the skin friction and heat transfer. The positive buoyancy force (beyond a certain value) causes an overshoot in the velocity profiles.A mass transfer constant - B magnetic field - C_fx skin friction coefficient in the x-direction - C_p specific heat at constant pressure, kJ.kg^–1.K - C_v specific heat at constant volume, kJ.kg^–1.K^–1 - E electric field - g acceleration due to gravity, 9.81 m.s^–2 - Gr Grashof number - h heat transfer coefficient, W.m².K^–1 - Ha Hartmann number - k thermal conductivity, W.m^–1.K - L characteristic length, m - M magnetic parameter - Nu_x local Nusselt number - p pressure, Pa, N.m^–2 - Pr Prandtl number - q heat flux, W.m^–2 - Re Reynolds number - Re_m magnetic Reynolds number - T temperature, K - T_o constant plate temperature, K - u,v velocity components, m.s^–1 - V characteristic velocity, m.s^–1 - x,y Cartesian coordinates - thermal diffusivity, m².s^–1 - coefficient of thermal expansion, K^–1 - , transformed similarity variables - dynamic viscosity, kg.m^–1.s^–1 - ₀ magnetic permeability - kinematic viscosity, m².s^–1 - density, kg.m^–3 - buoyancy parameter - electrical conductivity - stream function, m².s^–1 - dimensionless constant - dimensionless temperature, K - w, conditions at the wall and at infinity     

Uniqueness and the maximum principle for quasilinear elliptic equations with quadratic growth conditions

In this paper, we show that the maximum principle holds for quasilinear elliptic equations with quadratic growth under general structure conditions.Two typical particular cases of our results are the following. On one hand, we prove that the equation (1) {ie77-01} where {ie77-02} and {ie77-03} satisfies the maximum principle for solutions in H ¹()L(), i.e., that two solutions u ₁, u ₂H¹() L() of (1) such that u ₁u₂ on , satisfy u ₁u₂ in . This implies in particular the uniqueness of the solution of (1) in H ₀ ¹ ()L().On the other hand, we prove that the equation (2) {ie77-04} where fH^–1() and g(u)>0, g(0)=0, satisfies the maximum principle for solutions uH¹() such that g(u)¦Du|{²L¹(). Again this implies the uniqueness of the solution of (2) in the class uH ₀ ¹ () with g(u)¦Du|{²L¹().In both cases, the method of proof consists in making a certain change of function u=(v) in equation (1) or (2), and in proving that the transformed equation, which is of the form (3) {ie77-05}satisfies a certain structure condition, which using ((v₁ -v ₂)⁺)ⁿ for some n>0 as a test function, allows us to prove the maximum principle.     

Qualitative behavior of dissipative wave equations on bounded domains

The qualitative behavior of solutions of the mixed problem u_tt = u-a(x)u_t in IR x , u=0 on IR x , is studied in the case when a>0 and IRⁿ is bounded. Roughly speaking, if aa_min>0, then solutions decay at least as fast as exp t( –1/2a_min), with the possible exception of a finite dimensional set of smooth solutions whose existence is associated with a phenomenon of overdamping. If a_max is sufficiently small, depending on , then no overdamping occurs.Partially supported by NSF grant NSF GP 34260.This work was partially supported by the National Science Foundation under Grant No. GP 34260     

Effect of the specific heat ratio on the stability and laminar-turbulent transition of a supersonic boundary layer

At high supersonic velocities the specific heat ratio of the flow impinging on a body may change in connection with the considerable increase in stagnation temperature and chemical reactions (strictly speaking, in this case the flow can be described by means of the perfect gas model with an effective specific heat ratio _ef; as a rule, _gf<). This may entail a change in the point of laminar-tubulent boundary layer transition on the model. This paper is concerned with the determination of the effect of the specific heat ratio on boundary layer transition.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 179–183, March–April, 1989.The author wishes to thank M. I. Yaroslavtsev and V. A. Dmitriev for assisting with the experiments.     

Direct derivation of the two-shaft system balance conditions for the second order reciprocal inertia forces on the V-type eight cylinders internal combustion engines with a plane crankshaft

By employing the four shafts balance concept paper [1] has reported a balance regime for the second order reciprocal inertia forces on the V-type eight cylinder internal combustion engines with a plane crankshaft. Thereafter, paper [2] has acquired a two-shafts balance regime, but through a rather tedious roudabout degenerating manipulation. The present article has, but starting out directly from the two-shafts balance concept, successfully acquired the same results as those in paper [2]. In addition, we propose, herein, a third balance system which might be, in general, called the slipper balance regime.     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters A interfacial area of the - interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the - interface contained within the averaging volume, m² - A ^* interfacial area of the - interface contained within a unit cell, m² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - p pressure in the -phase, N/m² - p intrinsic phase average pressure for the -phase, N/m² - p –p , spatial deviation of the pressure in the -phase, N/m² - r ₀ radius of the averaging volume and radius of a capillary tube, m - v velocity vector for the -phase, m/s - v phase average velocity vector for the -phase, m/s - v intrinsic phase average velocity vector for the -phase, m/s - v –v , spatial deviation of the velocity vector for the -phase, m/s - V averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ Greek Letters V/V, volume fraction of the -phase - mass density of the -phase, kg/m³ - viscosity of the -phase, Nt/m² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1