“Thermal layer” effect in MHD: Interaction of a parallel shock with a layer of decreased density

Interaction of a parallel fast MHD shock with a layer of decreased density is discussed using ideal MHD approach. This is an extrapolation of gas dynamic thermal layer effect on ideal MHD. Computer simulations show that a magnetic field of a moderate intensity ( 1) may change the character of the flow for intermediate Mach numbers (M 5) and a new raking regime may occur which is not observed in the absence of a magnetic field. Self similar precursor analogous to that in gas dynamics may develop in the case of highM and low density in the layer but magnetic forces essentially decrease its growth rate. This problem appears in connection with cosmical shock propagation where planetary magnetic tails play the role of the thermal layer, and it may also be observed in the laboratory when the shock is strong enough to heat the walls ahead of it.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Shock heating of plasma in a rapidly increasing magnetic field

The experimental excitation of intense collisionless shock waves (M 5) with subsequent plasma compression by the magnetic field of a shock coil is described. A magnetic plug > 20 kOe is produced in 100 × 10^–9 sec by a current generator, a long line with 250-kV water insulation and a characteristic impedance of l At an initial deuterium-plasma density of 2 × 10¹⁴ cm^–3, shock waves with a front width of 20c/₀₃and a velocity of 5 × 10₇ cm/sec are recorded. The ion energy after the accumulation, determined from the neutron yield, turns out to be 2 ke V. Axial shock waves excited by the plasma flow beneath the shock coil are observed.Translated from Zhurnal Prikladnoi Mekhaniki i Teknicheskoi Fiziki, Vol. 11, No. 2, pp. 28–38, March–April, 1970.The authors thank G. I. Budker and R. Z. Sagdeev for formulating the problem, R. I. Soloukhin for interest in the study, and S. P. Shalamov for construction of the apparatus.     

Waves of finite amplitude in a hot plasma

A study is made of a plane shock wave of arbitrary strength propagating in a hot rarefied plasma across the magnetic field. The question of the propagation of nonstationary waves of finite but small amplitude under these conditions is examined.Fairly detailed studies have been made of waves of finite amplitude in a cold rarefied plasma. The profile of such waves is formed as the result of nonlinear and dispersion effects, the dispersion effects being caused by electron inertia and plasma anisotropy. If the gas-kinetic pressure of the plasma is taken into account, then dispersion effects appear which are associated with the fact that the Larmor radius of the ions is finite. Stationary waves of small but finite amplitude propagating across the magnetic field in a hot plasma (when the gas-kinetic pressure p is comparable with the magnetic pressure H²/87) have been treated in [1, 2]. In [1] an isolated rarefaction wave was found in a hot plasma, instead of the compression wave characteristic of a cold plasma, and a qualitative picture of the shock wave structure was given. In [2] a study was made of a small-amplitude shock wave with the finite size of the ion Larmor radius taken into account. The present paper investigates the structure of shock waves of arbitrary strength which propagate across the magnetic field in a fairly hot rarefied plasma, and also examines nonstationary waves of finite but small amplitude excited in a plasma by a magnetic piston acting over a limited time interval.Notation p gas-kinetic pressure - H magnetic field - u, v macroscopic velocities along the x and y axes - density - m_e(m_i) mass of electron (ion) - plasma conductivity - _H ion-cyclotron frequency - V_A Alfvèn velocity - c velocity of light - adiabatic exponent - V specific volume - _0e(_0i) electron (ion) plasma frequency - S₀ velocity of sound. In conclusion the author thanks R. Z. Sagdeev and N. N. Yanenko for discussing the paper, and also R. N, Makarov for helping with the numerical computations.     

Motion of a conductive piston in a channel with variable inductance

The motion of a conductive piston in the channel of a magnetohydrodynamic (MHD) generator of the conduction type with compound electrodes is considered. Formulas are obtained for calculation of the energy characteristics of the pulse MHD generator for various operational regimes. It is shown that in an MHD generator at magnetic Reynolds number values Re_m = ₀u₀ 1 (where ₀ is the permeability of a vacuum, is the electrical conductivity of the piston, u₀ is the initial velocity, and is the characteristic dimension), the energy transferred to an ohmic load may significantly exceed the values obtained in [1, 2]. Conditions for high-efficiency transformation of piston kinetic energy to electrical energy are considered for limiting values of the ratio of the latter to initial magnetic field energy in the generator channel.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 41–46, November–December, 1973.The authors thank V. I. Yakovlev for his helpful evaluation.     

Steady flow of a liquid metal in a rectangular MHD channel

The known experimental studies of steady flows of a liquid metal in magnetohydrodynamic (MHD) channels of rectangular section [1–4] were performed only for a few values of the Reynolds number, which does not permit a clear delineation of the fundamental governing laws of the flow in the zone of transition from laminar to turbulent flow. In addition, the study of turbulent MHD flows has been limited to two-dimensional channels.Below we present some results of experimental studies of the effect of a transverse magnetic field on the resistance coefficient for mercury flow in an MHD channel with side ratio 1 to 2.5. The choice of a channel with this side ratio was dictated by the need for studies of the intermediate case between flows in two-dimensional and square channels, which differ significantly from one another because of the different effect of the walls parallel to the magnetic field. In our studies, for each value of the Hartmann number the investigations were made for 30–50 values of the Reynolds number.Notation B₀ flux density of the applied magnetic field - M Hartmann number - R Reynolds number - _tm resistance factor of turbulent MHD flow - _* critical value of the resistance factor - geometric parameter of channel - the component of resistance factor in ordinary hydrodynamics due to pulsations - normed function - electric conductivity of metal - viscosity of metal - R₀ shydraulic radius - N smagnetic field parameter     

Supersonic nonequilibrium flow past segmental bodies of a mixture simulating the atmosphere of venus

Calculations of the flow of the mixture 0.94 CO₂+0.05 N₂+0.01 Ar past the forward portion of segmentai bodies are presented. The temperature, pressure, and concentration distributions are given as a function of the pressure ahead of the shock wave and the body velocity. Analysis of the concentration distribution makes it possible to formulate a simplified model for the chemical reaction kinetics in the shock layer that reflects the primary flow characteristics. The density distributions are used to verify the validity of the binary similarity law throughout the shock layer region calculated.The flow of a CO₂+N₂+Ar gas mixture of varying composition past a spherical nose was examined in [1]. The basic flow properties in the shock layer were studied, particularly flow dependence on the free-stream CO₂ and N₂ concentration.New revised data on the properties of the Venusian atmosphere have appeared in the literature [2, 3] One is the dominant CO₂ concentration. This finding permits more rigorous formulation of the problem of blunt body motion in the Venus atmosphere, and attention can be concentrated on revising the CO₂ thermodynamic and kinetic properties that must be used in the calculation.The problem of supersonic nonequilibrium flow past a blunt body is solved within the framework of the problem formulation of [4].Notation V body velocity - shock wave standoff - universal gas constant - ratio of frozen specific heats - hRt/m enthalpy per unit mass undisturbed stream P pressure - density - T temperature - m molecular weight - c_p specific heat at constant pressure - (X) concentration of component X (number of particles in unit mass) - R body radius of curvature at the stagnation point - _j rate of j-th chemical reaction shock layer P V ² pressure - density - TT temperature - mm molecular weight Translated from Izv. AN SSSR. Mekhanika Zhidkosti i Gaza, Vol. 5, No. 2, pp. 67–72, March–April, 1970.The author thanks V. P. Stulov for guidance in this study.     

Instability of rarefaction shocks in systems of conservation laws

It is well-known that rarefaction shocks are unstable solutions of nonlinear hyperbolic conservation laws. Indeed, for scalar equations rarefaction shocks are unstable in the class of smooth solutions, but for systems one can only say in general that rarefaction shocks are unstable in the larger class of weak solutions. (Here unstable refers to a lack of continuous dependence upon perturbations of the initial data.) Since stability in the class of weak solutions is not well understood, ([T, TE]), entropy considerations have played a leading role in ruling out shocks that violate the laws of physics. However, for non-strictly hyperbolic systems the analogy with the equations of gas dynamics breaks down, and general entropy or admissibility criteria for the variety of shocks which appear, (see, e.g., [IMPT]), are not known. In this paper we address the question of when the instability of a shock can be demonstrated within the class of smooth solutions alone. We show by elementary constructions that this occurs whenever there exists an alternative solution to the Riemann problem with the same shock data which consists entirely of rarefaction waves and contact discontinuities with at least one non-zero rarefaction wave. We show that for 2×2 strictly hyperbolic, genuinely nonlinear systems the condition is both necessary and sufficient. We show too that for the full 3×3 (Euler) equations of gas dynamics with polytropic equations of state, rarefaction shocks of moderate strength are unstable in the class of smooth solutions if and only if the adiabatic gas constant satisfies 1 < < 5/3 (see Theorem 8). More precisely, there is a constant y _*, 0 < y _* < 1, depending only on , such that if y _* p _lp _rp _l for 1-shocks, and if y _* p _rP _lp _r for 3-shocks (where p _r and p _l denote the pressures on both sides of the rarefaction shock), then the shock is unstable if and only if 1 < < 5/3. Thus for such shocks, the theory of the Riemann problem for polytropic gases in the range 1 < < 5/3 can be rigorously developed with a knowledge of the smooth solutions alone by using stability under smoothing as an admissibility criterion, rather than by using the classical entropy inequalities.     

Electrode layers in flows of an inviscid plasma with good conductivity

The structure of the electromagnetic electrode layers that are produced in flows across a magnetic field by a completely ionized and inviscid plasma with good conductivity and a high magnetic Reynolds number is examined in a linear approximation. Flow past a corrugated wall and flow in a plane channel of slowly varying cross section with segmented electrodes are taken as specific examples. The possibility is demonstrated of the formation of nondissipative electrode layers with thicknesses on the order of the Debye distance or electron Larmor radius and of dissipative layers with thicknesses on the order of the skin thickness, as calculated from the diffusion rate in a magnetic field [2].In plasma flow in a transverse magnetic field, near the walls, along with the gasdynamie boundary layers, which owe their formation to viscosity, thermal conductivity, etc. (because of the presence of electromagnetic fields, their structures may vary considerably from that of ordinary gasdynamic layers), proper electromagnetic boundary layers may also be produced. An example of such layers is the Debye layer in which the quasi-neutrality of the plasma is upset. No less important, in a number of cases, is the quasi-neutral electromagnetic boundary layer, in which there is an abrupt change in the frozen-in parameter k=B/p (B is the magnetic field and p is the density of the medium). This layer plays a special role when we must explicitly allow for the Hall effect and the related formation of a longitudinal electric field (in the direction of the veloeiryv of the medium). We will call this the magnetic layer. The magnetic boundary layer can be dissipative as well as noudissipative (see below). The dissipative magnetic layer has been examined in a number of papers: for an incompressible medium with a given motion law in [1], for a compressible medium with good conductivity in [2], and with poor conductivity in [3]. In the present paper, particular attention will be devoted to nondissipative magnetic boundary layers.     

Evolution of the Diffusion Mixing Layer of Two Gases upon Interaction with Shock Waves

A mathematical model of mechanics of a twovelocity twotemperature mixture of gases is developed. Based on this model, evolution of the mixing layer of two gases with different densities under the action of shock and compression waves is considered by methods of mathematical simulation in the onedimensional unsteady approximation. In the asymptotic approximation of the full model, a solution of an initialboundary problem is obtained, which describes the formation of a diffusion layer between two gases. Problems of interaction of shock and compression waves with the diffusion layer are solved numerically in the full formulation. It is shown that the layer is compressed as the shock wave traverses it; the magnitude of compression depends on shockwave intensity. As the shock wave passes from the heavy gas to the light gas, the mixing layer becomes overcompressed and expands after shockwave transition. The wave pattern of the flow is described in detail. The calculated evolution of the mixinglayer width is in good agreement with experimental data.     

Construction of exact numerical solutions of the stationary traveling wave type for viscous thin films

An effective numerical procedure, based on the Galerkin method, for finding solutions of the stationary traveling wave type in the complete formulation is proposed for the case of viscous liquid films. Examples of a viscous film flowing freely down a vertical surface have been calculated. The calculations have been made for various values of the dimensionless surface tension , including =0. The method makes it possible to predict a number of bifurcations that occur as decreases. The existence of numerous families of stationary traveling waves when 1 was demonstrated in [6]. The present study shows that as 1 all but one of these families of wave solutions disappear. The shape of the periodic and solitary waves and the pressure distribution in the film are found for various . When =0 and the wave number is fairly small, the periodic solution has a singularity, as predicted in [14]: at the crest of the wave a corner point appears; the angle between the tangents at this point =140–150. The method proposed can be used to calculate other wavy film flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 94–100, May–June, 1990.     

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

An “ideal equivalent gas method” for the study of shock waves in supersonic real gas flows

Summary A method developed by the author for the systematic study of the thermodynamic and dynamic properties of the gas behind a shock wave is reported.The method is applicable to supersonic flow regimes for which the excitation, dissociation and ionization effects invalidate the usually adopted hypothesis of ideal gas.An Ideal Equivalent Gas, having the ratio of the specific heats _s dependent on Mach number and altitude of flight is postulated.On the basis of the mass, momentum and energy conservation equations, valid through the shock wave, the relations defining the thermodynamic and dynamic state of the gas behind the shock wave are derived. These relations establish an extension of the classic relations valid for the ideal gas and reduce to them identically for _s=.The dependence of the ratio of specific heats _s of the Ideal Equivalent Gas on Mach number and altitude has been established, over a wide range, on the basis of the real gas solutions derived by Huber.

---

Sommario Nella presente nota viene esposto un metodo sviluppato dall'autore per lo studio sistematico dello stato termodinamico e dinamico del gas a valle di un'onda d'urto in regime supersonico, allorchè cioè gli effetti dell'eccitazione dei gradi di libertà vibrazionali delle molecole e della loro dissociazione e successiva ionizzazione invalidano l'ipotesi di gas ideale generalmente adottata.Viene definito un gas ideale equivalente avente rapporto dei calori specifici _s funzione del numero di Mach e della quota di volo ed in base alle equazioni di conservazione della massa, della quantità di moto e dell'energia, valide attraverso all'onda d'urto, vengono derivate delle relazioni definenti lo stato termodinamico e dinamico del gas a valle dell'onda d'urto. Tali relazioni costituiscono una estensione delle classiche relazioni dell'urto valide per il gas ideale alle quali si riducono per _s=.La dipendenza del rapporto dei calori specifici _s del gas ideale equivalente, dal numero di Mach e dalla quota è stata stabilita sulla base delle soluzioni ottenute da Huber per il gas reale.

---

     

Electrical parameters of a channel with finite electrodes with allowance for the electrode potential drop

Spatial problems involving the electric field in an MHD channel were formulated in [1] with allowance for the electrode potential drop. It was assumed that the electrode layer had a small thickness, so that relationships on the boundary of the layer could be applied to the surface of the electrode. It was assumed that the electrode potential drop ° could be represented as a function of the current density jn at the electrode in the form of a known function ° =f (j_n) determined experimentally or deduced from the appropriate electrode-layer theory. An approximate method was then put forward for solving such problems by reducing them to the determination of the electric field from a known distribution of the magnetic field and the gas-dynamic parameters. It was shown that when =°/ E is small (E is the characteristic induced or applied potential difference), the solution can be sought in the form of series in powers of . In the zero-order approximation, the electric field is determined without taking into account the electrode processes. The first approximation gives a correction of the order of . The quantity °, which is present in the boundary conditions on the electrode in the first-order approximation, is determined from the current density calculated in the zero-order approximation.One of the problems discussed in [1] was concerned with the electric current in a channel with one pair of symmetric electrodes. Its solution was found in the first approximation in the form of the integral Keldysh-Sedov formula. In this paper we report an analysis of the solution for ° taken in the form of a step function.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

A Numerical Study of the Critical Gas Saturation in a Porous Medium

We use porenetwork simulations to study the dependence of the critical gas saturation in solutiongas drive processes on the geometric parameters of the porous medium. We show that for a variety of growth regimes (including global and local percolation, instantaneous and sequential nucleation, and masstransfer driven processes), the critical gas saturation, Sgc, follows a powerlaw scaling with the final nucleation fraction (fraction of sites activated), fq. For 3D processes, this relation reads Sgcfq0.16, indicating a sensitive dependence of Sgc to fq at very small values of fq.     

A third order scheme in the vlasov theory and some consequences in the evolution of few body systems

Summary A set of implemented evolution equations, describing the coherent nonlinear interaction of plasma waves, based on the perturbation method, has been derived, taking into account initial value effects and third order nonlinearities in the modal amplitudes.The equations reduce, in the appropriate limit, to well known stochastic triplets of hydrodynamic type.It is argued that, the stochastization mechanism of the decay instability in strongly damped regime, may be interpreted as a Duffing-type behavior.

---

Sommario Si presenta un sistema di equazioni di evoluzione descrivente l'interazione coerente nonlineare di onde di plasma, tenendo conto di effetti di condizioni iniziali e nonlinearità del terzo ordine nelle ampiezze modali. Le equazioni si riducono, nel limite appropriato, a ben noti tripletti stocastici di tipo idrodinamico.Si ipotizza che il processo di stocastizzazione della instabilità di decadimento, in regime di forte dissipazione, possa essere interpretato da un modello dinamico del tipo Duffing.

---

Lyst of greek symbols omega (lower case)=frequency - delta (1.c.)=partial derivative - pi (1.c.)=greek pi (=3,14r. units) in (26) - Sigma (Capital case)=Summation symbol - Delta (C.c)=def. as in (14) - epsilon (1.c.) def. as in (10) - sigma (1.c.)=def. as in (36) - gamma (1.c.)=def. as in (37) - beta (1.c.)=def. as in (37) - ro, (1.c.)=def. as in (37) - alfa (1.c.)=def. in connection with _k in (60) - Gamma (C.c.)=def. as in (63) and (66), (67) - delta (1.c.)=def. as in (66) - psi (1.c.)=def. as in (75) - fi (1.c.)=def. in (71) in connection with ₀ - Fi (C.c.)=def. as in (75) - Omega (C.c.)=def. in (78)     

Das erweiterte Korrespondenzgesetz für die dynamische Idealviskosität und die Idealwärmeleitfähigkeit reiner Stoffe

Zusammenfassung Zur Berechnung der dynamischen Idealviskosität Ideal (T) und der Idealwärmeleitfähigkeit ideal (T) benötigt man die kritische TemperaturT _kr, das kritische spezifische Volum _kr, die MolmasseM, den kritischen Parameter _kr und die molare isochore WärmekapazitätC _v(T). Sowohl das theoretisch, als auch das empirisch abgeleitete erweiterte Korrespondenzgesetz ergeben eine für praktische Zwecke ausreichende Genauigkeit für die Meßwertwiedergabe, die bei den assoziierenden Stoffen und den Quantenstoffen jedoch geringer ist als bei den Normalstoffen.

---

The extended correspondence law for the ideal dynamic viscosity and the ideal thermal conductivity of pure substances

For the calculation of the ideal dynamic viscosity Ideal (T) and the ideal thermal conductivity ideal (T) the critical temperatureT _kr, the critical specific volumev _kr, the molecular massM, the critical parameter _kr, and the molar isochoric heat capacityC _v(T) is needed. Not only the theoretically determined but also the empirically determined extended correspondence law gives for practical use a good representation of the measured data, which for the associating substances and the quantum substances is not so good as for the normal substances.

     

Flow Structure in the Separation Zone in the Interaction between Shock Waves and the Boundary Layer on a Plate with Slip

A model which makes it possible to calculate the reverse-flow parameters in the separation zone is constructed on the basis of the results of an integrated experimental study of the characteristics of the separated flow developed in the transition from free to non-free interaction between plane shock waves and the boundary layer on a plate with slip. The effect of the Mach number of the reverse flow in the separation zone on the properties of inner boundary layer separation is analyzed. Features of the interference flow due to boundary layer transition are described. The present study is a continuation of investigations [1–3] devoted to the study of a new steady-state type of interaction between shock waves and the boundary layer on a plate with slip in which the separation line formed would propagate upstream beyond the sharp leading edge if no leading edge was present, i.e., so-called non-free interaction.     

Interaction of a shock wave with a magnetic field in a medium with finite conductivity

We consider the process of the interaction of aplanar shock wave with a magnetic field (impact on a magnetic wall) in a medium having finite conductivity.The problem cannot be solved analytically in the general form. Numerical methods are used to study the problem. A computer is used to calculate the complete system of one-dimensional nonsteady equations of MHD with finite conductivity which depends on temperature in a nonlinear fashion. Results are also presented of particular analytic solutions obtained under simplifying assumptions.We discuss the dependence of the process dynamics on the magnitude of the magnetic field intensity and the law of variation of the medium conductivity with temperature.In the numerical calculations we note the formation of a T-layer, a phenomenon which occurs under definite conditions in unsteady MHD problems [1].In conclusion the authors wish to thank N. G. Basov, A. A. Samarskli, and O. N. Krokhin for posing the problem and for fruitful discussions, and also D. A. Gol'din and A. A. Ivanov for carrying out the numerical calculations.