Unsteady hydromagnetic thermal boundary layer flow

With viscous dissipation and Joule heating taken into account, solution of the energy equation is obtained for unsteady hydromagnetic thermal boundary layer flow past a porous wall (e.g., surface of a star) in presence of a transverse magnetic field, under the condition of zero heat transfer between the fluid and the boundary — the so-called plate thermometer problem in MHD. Solution of the problem, in the form of power series, is obtained under certain valid simplifying assumptions, when (i) the wall is subjected to a normal velocity of suction/injection which is proportional tot ^–1/2, and (ii) the wall has a velocity given by t ⁿ . The variation of temperature is shown graphically and is followed by a quantitative discussion therein also signifying the importance of retaining the Joule heating term in the energy equation.Part I is the article inAstrophysics and Space Science, Vol.45, No. 2, 1976, pp. 397–410.On study-leave from Defence Science Laboratory, Delhi, India.     

The magnetic field in a protostellar cloud

General forms of theB-p relation are investigated in both the isothermal and the non-isothermal regions. The magnetic flux dissipation either by ambipolar diffusion or by Ohmic dissipation has been studied. The rates of heating due to the magnetic dissipation processes have been calculated in comparison with the rate of compressional heating.The magnetic field strength is derived as a function of flux/mass ratio, mass, density, and geometry of the isothermal cloud. In the non-isothermal region, the temperature is added to the above-mentioned variables.It has been found that the magnetic flux starts to dissipate via ambipolar diffusion at neutral density ofn>3×10⁹ cm^–3. Ambipolar diffusion continues effective until reaching densities ofn>10¹¹ cm^–3, where Ohmic dissipation dominates. Under some conditions, the electrons evaporate from the grain surface atn>10¹³ cm^–3, while the ions are still adsorbed on the grain surfce. In this case, the magnetic flux loss returns to be influenced by ambipolar diffusion.The rates of heating by both Ohmic dissipation _OD and ambipolar diffusion _AD are found to be smaller than the rate of compressional heating _C in case of magnetic dissipation. Assuming that the magnetic field is frozen in the medium, then _C is smaller than both _OD and _AD . The above results of heating were found in the non-isothermal region.     

The evolution of the moon: A finite element approach

Thermal convection has considerable influence on the thermal evolution of terrestrial planets. Previous numerical models of planetary convection have solved the system of partial differential equations by finite difference methods, or have approximated it by parametrized methods. We have evaluated the applicability of a finite element solution of these equations. Our model analyses the thermal history of a self-gravitating spherical planetary body; it includes the effects of viscous dissipation, internal melting, adiabatic gradient, core formation, variable viscosity, decay of radioactive nucleides, and a depth dependent initial temperature profile. Reflecting current interest, physical parameters corresponding to the Moon were selected for the model.Although no initial basalt ocean is assumed for the Moon, partial melting is observed very early in its history; this is presumably related to the formation of the basalt maria. The convection pattern appears to be dominated by an L-2 mode. The present-day lithospheric thickness in the model is 600 km, with core-mantle temperatures close to 1600 K. Surface heat flux is 25.3 mW m^–2, higher than the steady state-value by about 16%.The finite element method is clearly applicable to the problem of planetary evolution, but much faster solution algorithms will be necessary if a sufficient number of models are to be examined by this method.Notation coefficient of thermal expansion - _ij Kronecker delta - absolute or dynamic viscosity - perturbation in temperature - thermal diffusivity - kinematic viscosity - density - stress tensor - B.P. before present - c specific heat at constant pressure or volume (Boussinesq approximation) - d depth of convection - E ^* activation energy for creep - g gravity - Ga billions of years - H(t) heat generation per unit mass per unit time at timet - k Boltzmann's constant - K mean thermal conductivity - Ma millions of years - p pressure - q heat flux - q _ss steady-state heat flux - Ra Rayleigh number - S volumetric heat sources, includes radioactivity and viscous dissipation - t time - T temperature - u verocity vector - V ^* activation volume for creep     

On the dissipation rates for Alfvén waves in the solar transition region

The present paper concerns Alfvén waves, in a resistive and viscous atmosphere, under a steep temperature gradient (Section 1). The dissipative Alfvén wave equation is deduced assuming uniform vertical background magnetic field, and allowing for arbitrary profiles of Alfvén speed, and viscous and resistive diffusivities as functions of altitude (Section 2). A three-parameter family of temperature profiles, allowing for independent choice of initial and asymptotic temperature, and of initial temperature gradient, is used to re-write the wave equation, with the temperature as the independent variable, instead of altitude (Section 3). It is shown that, for the conditions prevailing in the solar transition region between the chromosphere and corona, two approximations of the dissipative wave equations may be considered, the simplest leading to solution in terms of Gaussian hypergeometric functions (Section 4). The exact analytical solution allows calculation of the (i) velocity and (ii) magnetic field perturbations, (iii) kinetic, (iv) magnetic and (v) total energy density, (vi) energy flux, (vii) rate-of-strain and (viii) electric current, and (ix) viscous, (x) resistive and (xi) total rate of dissipation (Section 5). These are plotted versus temperature, across the transition region from the chromosphere to the corona, for the quiet and active Sun (Section 6). The feasibility of heating of the transition region by dissipation of Alfvén waves is discussed (Section 7), by comparing empirical heating rates, with theoretical values for a range of physical conditions, including initial velocity perturbations 5 to 15 km s ^–1, background magnetic field 12 to 120 G, wave periods 60 to 300 s, thickness of the transition region 100 to 300 km, resistive and anomalous diffusivities to 100 and viscous and turbulent diffusivities to 100 . The conclusion is that dissipation of Alfvén waves is not an effective heating mechanism for the transition region and corona, although it may be for the chromosphere (see Campos and Mendes, 1995, and references therein).     

Plasma flow around coronal loops

Current dissipation models of coronal loop heating are studied. Turbulent current dissipation is shown to lead to a time dependent process because of an enormous mass motion induced in the current layer. A stationary heating process involves only ohmic heating, which requires a large current layer. To insure MHD stability, the loop must be composed of many elements with the oppositely directed currents. A stationary current dissipation process induces the plasma motion across the magnetic field into the loop and down the loop with the speeds v 10⁴ cm s^–1 and v 10⁴ cm s^–1, respectively. The pressure of the loop is also estimated to be proportional to the current density: p/J=6.3 × 10^-8dyn/statamp.     

Numerical simulation of coronal heating by resonant absorption of Alfvén waves

The heating of coronal loops by resonant absorption of Alfvén waves is studied in compressible, resistive magnetohydrodynamics. The loops are approximated by straight cylindrical, axisymmetric plasma columns and the incident waves which excite the coronal loops are modelled by a periodic external driver. The stationary state of this system is determined with a numerical code based on the finite element method. Since the power spectrum of the incident waves is not well known, the intrinsic dissipation is computed. The intrinsic dissipation spectrum is independent of the external driver and reflects the intrinsic ability of the coronal loops to extract energy from incident waves by the mechanism of resonant absorption.The numerical results show that resonant absorption is very efficient for typical parameter values occurring in the loops of the solar corona. A considerable part of the energy supplied by the external driver, is actually dissipated Ohmically and converted into heat. The heating of the plasma is localized in a narrow resonant layer with a width proportional to ^1/3. The energy dissipation rate is almost independent of the resistivity for the relevant values of this parameter. The efficiency of the heating mechanism and the localization of the heating strongly depend on the frequency of the external driver. Resonant absorption is extremely efficient when the plasma is excited with a frequency near the frequency of a so-called collective mode.     

Observations of steady anomalous magnetic heating in thin current sheets

A faint steadily emitting loop-like structure has been observed by HXIS in its low energy channels (3.5–8.0 keV) on November 5/6, 1980. These HXIS observations have permitted us to follow the thermal evolution of this loop for a period of about 15 hr and from this study we conclude that only a fraction of 0.1% of the volume of the loop is steadily heated at the rather large rate of 0.6 erg cm^-3 s^-1. We interpret this heating as the dissipation of magnetic fields in thin current sheets and we find that the dissipation with classical resistivity is very unlikely, while ion-kinetic tearing, as proposed by Galeev et al. (1981), suits the observations very well. The enhancement of the resistivity over the classical resistivity then turns out to be a factor 4 × 10⁴. Dissipation in extremely thin sheets via the ion-acoustic instability (Duijveman et al., 1981) cannot be completely excluded when the cross-field heat conductivity is anomalously enhanced by a factor 400.We identify the source of the X-ray emission in this paper with the H filament in the same region. The hot X-ray emitting plasma and the cool plasma radiating in H are thermally separated by the strong magnetic field.The main conclusion of the paper is that for the first time direct evidence is found for the steady dissipation of coronal magnetic fields via enhanced resistivity in thin current sheets.     

TWO-FLUID MODEL FOR HEATING OF THE SOLAR CORONA AND ACCELERATION OF THE SOLAR WIND BY HIGH-FREQUENCY ALFVÉN WAVES

A model of the solar corona and wind is developed which includes for the first time the heating and acceleration effects of high-frequency Alfvén waves in the frequency range between 1 Hz and 1 kHz. The waves are assumed to be generated by the small-scale magnetic activity in the chromospheric network. The wave dissipation near the gyro-frequency, which decreases with increasing solar distance, leads to strong coronal heating. The resulting heating function is different from other artificial heating functions used in previous model calculations. The associated thermal pressure-gradient force and wave pressure-gradient force together can accelerate the wind to high velocities, such as those observed by Helios and Ulysses. Classical Coulomb heat conduction is also considered and turns out to play a role in shaping the temperature profiles of the heated protons. The time-dependent two-fluid (electrons and protons) model equations and the time-dependent wave-spectrum equation are numerically integrated versus solar distance out to about 0.3 AU. The solutions finally converge and settle on time-stationary profiles which are discussed in detail. The model computations can be made to fit the observed density profiles of a polar coronal hole and polar plume with the sonic point occurring at 2.4 R and 3.2 R , respectively. The solar wind speeds obtained at 63 R are 740 km s^-1 and 540 km s^-1; the mass flux is 2.1 and 2.2 × 10⁸ cm^-2 s^-1 (normalized to 1 AU), respectively. The proton temperature increases from a value of 4 × 10⁵ K at the lower boundary to 2 × 10⁶ K in the corona near 2 R .     

Hydromagnetic free convection and mass transfer flow of non-Newtonian fluid through a porous medium bounded by an infinite vertical limiting surface with constant suction

An analysis of a two-dimensional steady-free convection and mass transfer flow of an incompressible, viscous, and electrically conductive non-Newtonian fluid through a porous medium bounded by a vertical infinite limiting surface (plane wall) has been presented in the presence of a transverse magnetic field. Approximate solutions to the coupled nonlinear equations governing the flow are derived and expression for the velocity, temperature, concentration, the rate of heat transfer, and the skin-friction are derived. Effects of Gr (Grashof number), Gm (modified Grashof number),M ^* (non-Newtonian parameter),N (magnetic parameter), and permeabilityK of the porous medium on the velocity, the skin-friction and the rate of heat transfer are discussed when the surface is subjected to a constant suction velocity.     

Grain formation in the expanding gas flow around cool luminous stars

We studied grain formation process and flow structure around cool luminous mass-loss stars. The nucleation and growth theory of Yamamoto and Hasegawa was extended to the case of expanding gas flow.The envelope was assumed to be steady, spherically symmetric, in thermal and radiative equilibrium, optically thin, and driven by radiation pressure on grains. For oxygen rich stars, Mg-silicate was found to be the first condensate which can drive the gas effectively. The following stellar parameters were chosen; stellar massM _*=1M , effective temperatureT _*=3000K, stellar luminosityL _* from 7.5×10³ to 2.0×10⁴ L , and mass-loss rate |M| from 1.0×10^–6 to 1.0×10^–4 M yr^–1.Main results of calculations are as follows; (1) grain condensation temperatureT _c9801080 K; (2) total gas pressure at the condensation pointP _t6×10^–116×10^–9 atm; (3) scale parameterA _c10³6×10⁴; and (4) final grain sizer _f=400Å1m. For the smaller |M| or the largerL _*, these values are the smaller. We recognized two types of flow solutions (1) Dust driven flow for large |M|, which reaches the sonic point near the condensation point; and (2) Modified Parker flow for small |M| for which grain sizer _f is almost independent of |M|.A comparison with observational results ofL _* and gas terminal velocityV suggests that Mg-silicate grains are of submicron size, which are effective for interstellar extinction in visible and infrared. Fe-grains condense in the rarefied outflow with a size probably smaller than 100Å, which may contribute for interstellar ultraviolet extinction. The envelope has three-layer structure inner dense region with small outflow velocity, grain formation layer and outer rarefied region with fast outflow velocity.No flow solutions exist forM _* greater than a critical stellar massM _*cr for a given stellar luminosityL _* and mass-loss rate |M|.For example, critical stellar massM _*cr=1.8M forL _*=10⁴ L ,T _*=3000 K, and |M|=10^-5 M yr^-1.     

Is a Power-Law Form of the Correlation Gamma Function Sufficient to State that a Distribution Has Fractal Properties?

This investigation indicates an ambiguity in interpreting the results of applying the apparatus of the correlation gamma function [(r) and ^*(r)] to analyze the spatial distribution of objects from some sample. The presence of a linear section in the dependences of log() on log(r) and of log(^*) on log(r) proves to be insufficient to state that the distribution has fractal properties (self-similarity). It is shown that a change in the geometrical boundaries of the sample may influence the form of the dependences of log() on log(r) and of log(^*) on log(r).     

Thermal properties of the solar wind plasma

In the solar wind, electrostatic ion cyclotron waves can be excited by electrons when the flow velocity becomes supersonic. The waves reduce the proton temperature anisotropy and heat the protons effectively. Temperature equations for T _e ,T _p and T _p are solved numerically in the region from 1 AU to the Sun, with the non-thermal proton heating rate included as a parameter. Distributions of T _e ,T _p , T _p and the proton heating rate are determined and found to be in good agreement with the proton heating rate expected from the linear growth rate of electrostatic ion cyclotron waves. The electron thermal conductivity is reduced approximately 2–3 times smaller than the usual collisional one due to the plasma wave instabilities. Effective energy exchange rates from proton-proton and electron-proton interactions are 1–10 and 10–100 times larger than the Coulomb collision rates v ppand v ep,respectively.     

Correlation Properties of a Distribution of Objects from the PSCz Survey

The properties of a sample of IR galaxies from the Point Source Catalog z (PSCz) Survey were investigated using the correlation gamma function [(r) and ^*(r)]. The results with different volume-limited subsamples indicate that the regions of strong correlations (a power-law decrease in density with distance with an exponent ₁ 1) are limited to a scale of 10-15 Mpc. A break is present at this scale in the dependences of log() on log(r) and of log(^*) on log(r). Such a break is also observed for various other samples of galaxies and clusters (at different scales). After the break the density dependence changes to another regime, corresponding to a fully uniform distribution for bright galaxies in the northern Galactic hemisphere. For some subsamples of the southern and northern hemispheres the latter regime corresponds to some decrease in density with distance. Indications of significant differences between the distributions of objects in the southern and northern hemispheres are obtained. It is shown that the section of the gamma function after the break, even when its extent is small, is a significant indicator of actual correlation properties of the distributions. The results of the analysis are in good agreement, on the whole, with preceding studies of the PSCz survey.     

Study of alpha component dynamics in the solar wind using the Prognoz satellite

The analysis of solar wind He⁺⁺ and H⁺ ion distribution functions, collected over five months by the satellite Prognoz 1, shows that these are in general maxwellian but that often tails appear at higher speeds. The existing relation V-T, the observation of ratios of T/Tp 3.83 and V/Vp 1.035 give evidence of preferential He⁺⁺ ion heating and acceleration. The criteria for heating by dissipation of hydromagnetic waves proposed by Barnes and Hung (1973) are tested experimentally. Finally, multifluid models are likely to predict certain observations such as dependence of the velocity ratio V/Vp on the solar wind flux.     

Tidally Induced Volcanism

The dissipation of tidal energy causes the ongoing silicate volcanism on Jupiter's satellite, Io, and cryovolcanism almost certainly has resurfaced parts of Saturn's satellite, Enceladus, at various epochs distributed over the latter's history. The maintenance of tidal dissipation in Io and the occurrence of the same on Enceladus depends crucially on the maintenance of the respective orbital eccentricities by the existence of mean motion resonances with nearby satellites. A formation of the resonances among the Galilean satellites by differential expansion of the satellite orbits from tides raised on Jupiter by the satellites means the onset of the volcanism on Io could be relatively recent. If, on the other hand, the resonances formed by differential migration from resonant interactions of the satellites with the disk of gas and particles from which they formed, Io would have been at least intermittently volcanically active throughout its history. Either means of assembling the Galilean satellite resonances lead to the same constraint on the dissipation function of Jupiter Q _J 10⁶, where the currently high heat flux from Io seems to favor episodic heating as Io's eccentricity periodically increases and decreases. Either of the two models might account for sufficient tidal dissipation in the icy satellite Enceladus to cause at least occasional cryovolcanism over much of its history. However, both models are assumption-dependent and not secure, so uncertainty remains on how tidal dissipation resurfaced Enceladus.     

Can gravitational effects damp Alfvén waves?

We show that Alfvén-gravity waves propagating in a gravitationally stratified atmosphere do not suffer damping as a result of the rate of working of the gravity drift current on the electric field of the waves. A self-consistent treatment involving conservation of total energy, Poynting's theorem, and the rate of working of the various drift currents on the electric field demonstrates that dissipation only arises from real dissipative processes such as Ohmic heating or viscous effects, otherwise the system is adiabatic.     

A model of impulsive loop flares caused by anomalous heat conduction

Numerical simulation is made of the impulsive loop flare caused by transient heat conduction along the loop with an applied axial electric current.It is assumed that a segment near the top of the coronal loop is heated to above 10⁷ K by a heat input that is small compared with the total flare energy, which is given by the magnetic energy of the initial current. Due to the heat conduction, a hump appears in the velocity distribution of electrons, which may excite electron plasma waves with a sufficiently high intensity to cause an anomalous resistivity, as shown theoretically in a previous paper. In that paper, an effect of the plasma waves on the dynamics of electrons was taken into account consistently, but an anomalous heating due to an ohmic dissipation of the initial current under the anomalous resistivity was not taken into account.The aim of the present study is to study the subsequent dynamics of the heated gas caused by the anomalous heating, but in order to avoid an unpractically long computation time, the energy density of the plasma waves is estimated by the energy density of electrons in the velocity hump, without taking into account the effect of the plasma waves consistently in the dynamics of the electrons.The initial current starts to decay gradually by an ohmic dissipation under the anomalous resistivity occurring near the top of the loop to heat this region more. The enhanced heat conduction causes the velocity humps in a wider location. Consequently, the anomalous heating continues and spreads in a self-generating way even after the end of the initial minor heating. Thus the temperature near the loop top becomes above 10⁸ K and the high-temperature region spreads in both directions along the loop with such a high speed as (2–3) × 10⁴ km s^–1, which is nearly equal to the speed of flux-limited heat conduction. On the other hand, induced electric field estimated from the anomalous resistivity is 3.3 × 10⁷ V at the termination of the present simulation, under the modest initial current of 1.5 A m^–2.X-ray emissions expected from the present model loop, show three sources, two footpoints with unequal brightness and a coronal source expanding along the loop in both directions.     

Characteristic actionsħ (s) in the structure of the universe and the problem of the origin of galaxies

From the characteristic actions ^(s), we can derive various relations between the basic characteristic quantities of objects and the fundamental constants in known physical laws. The main physical processes which lead to the formation of objects should be included in some such reletions through the fundamental constants. The problem of the origin of galaxies has been considered on the basis of the theory of actions ^(s). It has been shown that in addition to gravitational effect, the dissipation process of the adiabatic density perturbations arising from the Thomson scattering in the early universe is a crucial process in forming galaxies; and if the Hubble constant has a valueH ₀ 50 km s^–1 Mpc^–1, the protogalaxies might be formed just before recombination.     

Linear hydrodynamical equations coupled with radiative transfer in a non-isothermal atmosphere

A method coupling the hydrodynamical equations and radiative transfer in a realistic solar model atmosphere is described. The influence of the temperature gradient of the model and the radiative dissipation is pointed out.The effect of the large temperature gradient is important in the layers where the optical depth ₅₀₀₀ is greater than 0.5; the ratio between the amplitude of the temperature and the velocity fluctuations decreases with the altitude by a factor 2 between = 1 and = 0.5 and in the case of the acoustic waves, the phase shift between these fluctuations is small.The radiative energy loss in the thick layers ( ₅₀₀₀ = 1) leads to a decrease of the vertical phase velocity of the waves and to a damping of their amplitudes in the layers of intermediate optical depth (10^-2 < ₅₀₀₀ < 0.5). The effect of the dissipation is negligible in the thin layers (₅₀₀₀ < 10^-2).     

The thermal structure of the dayside upper atmosphere of venus above 125 km

A one-dimensional model of the Venus thermosphere has been constructed which includes computation of the heating efficiency of solar ultraviolet radiation, heat loss by radiation to space of infrared-active species, thermal transport by molecular and eddy conduction, and viscous dissipation. By comparing model predictions with results obtained from the Pioneer Venus Orbiter space-craft, the results indicate that energy transport parameterized by eddy heat conduction plays a dominant role in determining thermospheric temperature T. It is suggested that there exists a feedback mechanism linking heating and thermospheric circulation such that eddy cooling maintains an asymptotic temperature T_∞～300°K for both solar-maximum and solar-minimum conditions. We also study the variation in thermospheric temperature with solar zenith angle, atomic oxygen-mixing ratio, rate of vibrational excitation of CO₂ by ground-state O atoms, and the assumed transfer of O(¹D) electronic energy to CO₂ vibrational energy.