Effects of Hall current on hydromagnetic free-convective flow through a porous medium

An analysis of the effects of Hall current on hydromagnetic free-convective flow through a porous medium bounded by a vertical plate is theoretically investigated when a strong magnetic field is imposed in a direction which is perpendicular to the free stream and makes an angle to the vertical direction. The influence of Hall currents on the flow is studied for various values of .Nomenclature c _p specific heat at constant pressure - e electrical charge - E Eckert number - E electrical field intensity - g acceleration due to gravity - G Grashof number - H ₀ applied magnetic field - H magnetic field intensity - (j _x , j _y , j _z ) components of current densityJ - J current density - K permeability of porous medium - M magnetic parameter - m Hall parameter - n _e electron number density - P Prandtl number - q velocity vector - (T, T _w , T ) temperature - t time - (u, v, w) components of the velocity vectorq - U ₀ uniform velocity - v ₀ suction velocity - (x, y, z) Cartesian coordinates Greek Symbols angle - coefficient of volume expansion - _e cyclotron frequency - frequency - dimensionless temperature - thermal conductivity - coefficient of viscosity - magnetic permeability - kinematic viscosity - mass density of fluid - _e charge density - electrical conductivity - _e electron collision time     

On the A-transform of the exponential integrals

In this paper, a technique of recursive analysis is developed for the integral transform A of the exponential integral functionsE _n which is denoted as _n (). The main result of this analysis enables us to establish a two-term recurrence formula for _n (0) and a three-term recurrence formula for _n (); 0. A computational algorithm based on these formulae is also constructed and its numerical results forn=2(1)25 are presented to 15-digit accuracy.     

Shear angle of magnetic fields

In this paper we introduce a new parameter, the shear angle of vector magnetic fields, , to describe the non-potentiality of magnetic fields in active regions, which is defined as the angle between the observed vector magnetic field and its corresponding current-free field. In the case of highly inclined field configurations, this angle is approximately equal to the angular shear, , defined by Hagyardet al. (1984). The angular shear, , can be considered as the projection of the shear angle, , on the photosphere. For the active region studied, the shear angle, , seems to have a better and neater correspondence with flare activity than does . The shear angle, , gives a clearer explanation of the non-potentiality of magnetic fields. It is a better measure of the deviation of the observed magnetic field from a potential field, and is directly related to the magnetic free energy stored in non-potential fields.     

On gravitational orbits and the virial theorem

We emphasize the sharp distinctions between different one-body gravitational trajectories made by the ratio of time averagesR(t)–E _kin/E _pot.R is calculated as a function of the eccentricity (e) and of the energy (E). Whent, independently ofe andE, R1/2 for closed orbits (this clearly illustrates the fulfillment of the virial theorem in classical mechanics); whereasR1, at any time, for open orbits.     

Redshift cutoff and evolution of QSOs

In this paper I present a new evolution model of QSOs luminosity. The model is based on edges distribution of apparent magnitude-redshift of QSOs. After the quasars were formed, the luminosities were increasing until they attained their maximum value atz=2+a, where –0.1a0.6, then the luminosities were decreasing. If the QSOs originate from superconducting cosmic string of same initial massM _i 10¹² M , the formation epochs are different, most of the quasars start atz _cutoff5.6. The most luminous QSOs start at later epochz _cutoff5.15. The present sky survey echniques may give us the possibility to see the formation of QSOs at apparent magnitudem _V 22.5 by chance of 0.3%.     

On a transformation of the differential equations of the lunar theory

This work contains a transformation of Hill-Brown differential equations for the coordinates of the satellite to a type which can be integrated in a literal form using an analytical programming language. The differential equation for the parallax of the satellite is also established. Its use facilitates the computation of Hill's periodic intermediary orbit of the satellite and provides a good check for the expansion of the coordinates and frequencies. The knowledge of the expansion of the parallax facilitates the formation of differential equations for terms with a given characteristic. These differential equations are put into a form which favors the solution by means of iteration on the computer. As in the classical theory we obtain the expansions of the coordinates and of the parallax in the form of trigonometric series in four arguments and in powers of the constants of integration. We expand the differential operators into series in squares of the constants of integration. Only the terms of order zero in these expansions are employed in the integration of the differential equations. The remaining terms are responsible for producing the cross-effects between the perturbations of different order. By applying the averaging operator to the right sides of the differential equations we deduce the expansion of the frequencies in powers of squares of the constants of integration.Basic Notations f the gravitational constant - E the mass of the planet - M the mass of the satellite - t dynamical time - x, y, z planetocentric coordinates of the satellite - u x+y–1 - s x–y–1 - the planetocentric distance of the satellite - w 1/ - ₀ the variational part of - w ₀ the variational part ofw, - n the mean daily sidereal motion of the satellite - a the mean semi-major axis of the satellite defined by means of the Kepler relation:a ³ n ²=f(E+M) - a the mean semi-major axis defined as the constant factor attached to the variational solution - e the constant of the eccentricity of the satellite - the sine of one half the orbital inclination of the satellite relative to the orbit of the sun - c(n–n) the anomalistic frequency of the satellite - c ₀ the part ofc independent frome,e, and - g(n–n) the draconitic frequency of the satellite, - g ₀ the part ofg independent frome,e, and - exp (n–n)t–1 - D d/d - e the eccentricity of the solar planetocentric orbit - a the semi-major axis of the solar orbit - n the mean daily motion of the sun in its orbit around the planet - m n/(n–n) - a/a-the parallactic factor - the disturbing function     

Inverse compton effect and the colorimetric characteristics of flare stars

It is shown that the observed color diagrams(U-B) _f (B-V) _f for pure flare emission of UV Cet type flare stars may be explained within the framework of a fast electron hypothesis. We point out the essential influence on these color indices of the two following factors: (a) the deviations of the normal radiation capability of the star in the infrared region of spectra (on 3.6 m, 4.4 m, and 5.5 m) from the Planckian distribution; (b) the location of the cloud (source) of fast electrons around the star (flare geometry effect). Under the real conditions of the generation of flares around the star the frequency transformation law at the photon-electron interaction has a view =n²₀, wheren may take the different values-from 0.15 up to 4; it depends on the cloud-star-observer geometry. By the observed colors of the flare emission may be understood, in principle, the location of flare source around the star. A possible role of reflection effect at the generation of stellar flares is outlined.     

Effects of metal abundances on the evolutionary period changes of classical Cepheids

Some peculiarities in the behaviour of a model self-gravitating system described by hydrodynamical equations and isothermal equation of state connected with the presence of thermodynamical fluctuations in real systems were investigated in numerical experiment. The values of density and velocity , , respectively, were computed by numerical code perturbed on each time-step and in each computational cell by random values , for modeling such fluctuations. Perturbed values _i = _i + _i ,v _i = _i + v _i were used to initiate the next step of computations. This procedure is equivalent to an introduction into original hydrodynamical equations of Langevin sources which are random functions. It is shown that these small fluctuations (= v =0,² =v ² = 10^–8) grow many times in marginally-stable state.     

Expanding and superdense astrophysical systems in general relativistic hierarchical cosmology

A semi-continuous hierarchy, (i.e., one in which there are galaxies outside clusters, clusters outside superclusters etc.), is examined using an expression of the field equations of general relativity in a form due to Podurets, Misner and Sharp. It is shown (a) that for a sufficiently populous hierarchy, the thinning factor( _i+1/ _i [r _i /r _i+1] is approximately equal to the exponentN in a continuous density law (=aR ^–N) provided (r _i /r _i+1)^3-1; (b) that a hierarchical Universe will not look decidedly asymmetric to an observer like a human being because such salient observers live close to the densest elements of the hierarchy (viz stars), the probability of the Universe looking spherically symmetric (dipole anisotropy0.1 to such an observer being of order unity; (c) the existence of a semi-continuous or continuous hierarchy (Peebles) requires that 2 if galaxies, not presently bound to clusters were once members of such systems; (d) there are now in existence no less than ten arguments for believing 2, though recent number counts by Sandageet al. seem to be in contradiction to such a value; (e) Hubble's law, withH independent of distance, can be proved approximately in a relativistic hierarchy provided (i)N=2, (ii)2GM(R)/c ² R1; (iii)Rc (iv)M0 in a system of massM, sizeR (f) Hubble's law holds also in a hierarchy with density jumps; (g)H100 km s^–1 Mpc^–1; (h) objects forming the stellar level of the hierarchy (in a cosmology of the Wilson type) must once have had 2GM/c ² R1; (i) there is a finite pressurep=2Ga in all astrophysical systems (a=R ^N ,N2); (j) for the Galaxy, theory predictsp _G7×10^–12 dyn cm^–2, observation givesp _G5×10^–12 dyn cm^–2; (k) if the mass-defect (or excess binding energy) hypothesis is taken as a postulate, all non-collapsed astrophysical systems must be non-static, and any non-static, p0 systems must in any case be losing mass; (1) the predicted mass-loss rate from the Sun is 10¹² g s^–1, compared to 10¹¹ g s^–1 in the observed solar wind; (m) the mass-loss rates known by observation imply timescales of 5×10⁹ years for the Sun and 10¹⁰ years for other astrophysical systems; (n) degenerate superdense objects composed of fermions must haveN-2 if they were ever at their Schwarzschild radii and comprised a finite numberN _B of baryons; (o)N _B10⁵⁷N for degenerate fermion and boson systems; (p)285-4; (q) the metric coefficients for superdense bodies give equations of motion that imply equal maximum luminosities for all evolving superdense bodies (L _max10⁵⁹ erg s^–1); (r) larger bodies have longer time-scales of energy radiation atL _max (10^–5 s for stars,1 h for QSO's) (s) expansion velocities are c soon after the initial loss of equilibrium in a superdense object; (t) if the density parametera(t) in aR ^–N isa=a (non-atomic constants of physicsc, G, A), andA, thenN=2; (u) N2 is necessary to giveMM at the stellar level of the hierarchy;(v) systems larger than, and including, galaxies must have formed by clumping of smaller systems and not (as advocated by Wertz and others) in a multiple big bang.     

The scattering of light by dielectric particles

We consider homogeneous particles with a frequency-dependent index of refractionn =n +in = 1 +(), ()1, in an outer electromagnetic field. From general considerations we derive an expression for the total scattering cross-section of these particles. We calculate the efficiency factor for the sphere, the infinite circular cylinder, and for an ellipsoidal particle with the main axesa=1,b=2,c=1. In contrast with the ray approximation, the present method can be used to calculate the polarization of the scattered light from ellipsoids.     

The radiative relaxation time in the chromosphere

The cooling effect of emission in the spectral lines, which dominates over continuous emission in the chromosphere and becomes important first around the temperature minimum, modifies greatly the radiative relaxation timet _r in the solar atmosphere. This rises from low photospheric values to a maximum of 8 min just aboveT _min, falls in the low chromosphere to 1.5 min because of line emission, but rises again to 6 min atT 7000–8400 K in the chromosphere where hydrogen ionization increases the specific heat.     

Rigorous estimates for the series expansions of Hamiltonian perturbation theory

In the present paper we prove a theorem giving rigorous estimates in the problem of bringing to normal form a nearly integrable Hamiltonian system, using methods of classical perturbation theory, i.e. series expansions in the small parameter . For any order of normalization, we give a lower bound _* ^r for the convergence radius of the normalized Hamiltonian, and a greater bound for the remainder, i.e. the non normalized part of the Hamiltonian. As an application, we consider the case of weakly coupled harmonic oscillators with highly nonresonant frequencies and show how, by optimizing, for fixed , the orderr of normalization, one gets for the remainder a greater bound of the formAe ^–(*₁/) ^a , with positive constantsA,a and ₁ ^* exponential estimate of Nekhoroshev's type.     

Search for the matter clumps in scalesZ ~ 1

The diagramV - log(1 +z _e ) as function of (, ) is considered for the quasars. HereV is the apparent visual magnitude,z _e is the emission line redshift, and are the equatorial coordinates. Two opposite extreme spots NE and SE are observed on the sky, where the inclination of the straight line fitting the dependenceV - log(1 +z _e ) is maximum and minimum. The coordinates of the centres of these extreme spots are ( _NE, _NE) = (282°, +42°) and ( _SE, _SE) = (70°, -38°) with errors 5°. A hypothesis of the Superattractor (SA) is proposed to explain such an effect. Two independent tests of this hypothesis are realized. First, the dependence or the frequency a of the absorbers in QSO spectra on (, ) is investigated. A region of the larger a is found. The coordinates of its centre are (, ) = (82°, - 10°) with error 5°. Second, the cases ofz _a >z _e are plotted in the Mercatorial projection (, ). The most of the casesz -z _e > 0.02 are concentrated within the circle with radiusR = 34° and centre (, ) = (50°, - 15°). The both anomalous regions overlap the Southern extreme spot around SE. The SA direction is (, ) = (67°, -21°) with errors about 12°. The redshift of SA isz _SA = 1.7 ± 0.3 that corresponds to the distancer _SA = (3100 ± 300)h ^–1 Mpc for the Hubble constantH ₀ = 75h kms^–1 Mpc^–1. The SA mass isM _SA ~ 10¹⁸-10²⁰ M . The orientation of the normal to the quasiperiodical large-scale sheet structure on the sky occurs near SA.     

Non-restricted double-averaged three body problem in Hill's case

A limiting case of the problem of three bodies (m ₀,m ₁,m ₂) is considered. The distance between the bodiesm ₀ andm ₁ is assumed to be much less than that between their barycenter and the bodym ₂ so that one may use Hill's approximation for the potential of interaction between the bodiesm ₁ andm ₂. In the absence of resonant relations the potential, double-averaged by the mean longitudes ofm ₁ andm ₂, describes the secular evolution of the orbits in the first approximation of the perturbation theory.As Harrington has shown, this problem is integrable. In the present paper a qualitative investigation of the evolution of the orbits and comparison with the analogous case in the restricted problem are carried out.The set of initial data is found, for which a collision between the bodiesm ₀ andm ₁ takes place.The region of the parameters of the problem is determined, for which plane retrograde motion is unstable.In a special example the results of approximate analysis are compared with those of numerical integration of the exact equations of the three body problem.

---

m ₀,m ₁,m ₂. , m ₀ m ₁. m ₂, m ₁ m ₂ m ₁ m ₂ . , . . , m ₀ m ₁. , . .

     

Solar particle propagation from 1 to 5 AU

A statistical analysis of solar particle events, observed by the GSFC-UNH charged particle detector on board Pioneer 10 and Pioneer 11 from March 1972 to December 1974 (from 1 to 5 AU for each spacecraft), is carried out with the goal of experimentally determining the statistical average interplanetary propagation conditions from 3 to 30 MeV. A numerical propagation model is developed that includes diffusion with a diffusion coefficient of the form k _r =k _o r , convection, adiabatic deceleration, and a variable coronal injection profile. The statistical analysis is carried out by individually analyzing each of five parameters (t _max, (t_max), t ₅, ) that are uniquely defined in a solar particle event. Each of the five parameter data sets were analyzed in terms of both a spacecraft-solar flare connection longitude 50°, and a numerical model that employed a variable exponential decaying coronal injection profile.The five individual parameter analyses are combined with the results that the statistical average radial interplanetary diffusion coefficient from 1 to 5 AU is given by k _r = (1.2 ± 0.4) × 10²¹ cm² s^-1 with = 0.0± 0.3 for 3.4 to 5.2 MeV protons and k _r = (2.6 ± 0.6) × 10²¹ cm² s^-1 with () = 0.0± 0.3 for 24 to 30 MeV protons. Using the classical relationship for the radial scattering mean free path _r, i.e. k _r = _r/3, we obtain _r = 0.09 ± 0.03 AU and 0.075 ± 0.020 AU for the low and high energy data, respectively. These results show, from 1 to 5 AU and from 3 to 30 MeV, that _r is both independent of radial distance and approximately independent of rigidity (for _r~P , where P = rigidity, = -0.15 ± 0.20).The above diffusion coefficients are inconsistent With both the predictions of the diffusion coefficient from present theoretical transport models and with the diffusion coefficient used in modulation studies at low energies.     

A momentum distribution of high energy particles around Crab pulsar

We determine the momentum distribution of the relativistic particles near the Crab pulsar from the observed X- and -ray spectra (10³10⁹ eV), provided that the curvature radiation is responsible for it. The power law spectrum for the relativistic electrons,f() ^–5, reproduces a close fit to the observed high-energy photon spectrum. The theoretically determined upper limit to the momentum (due to radiation damping), _M 8×10⁶, corresponds to the upper cut-off energy of the -ray spectrum, 10⁹ eV. The lower limit to the momentum, _m 1.8×10⁵, is chosen such that flattening of the X-ray spectrum below 10 keV is simulated. The number density of these electrons is found to be much higher than the Goldreich-Julian density. We also discuss pulse shape and polarization of high-energy photons. The extremely high density of particles and the steep momentum spectrum are difficult to understand. This may imply that another, more efficient, mechanism is in operation.     

Cyclotron waves in the solar wind

Cyclotron waves in the solar wind near 1 AU with frequencies well below the electron cyclotron frequency and wavelengths much larger than the electron cyclotron radius but less than the proton cyclotron radius are considered. The cyclotron radii are defined from parallel thermal velocity of electron component and proton component with respect to the interplanetary magnetic field. No LH cyclotron waves are found to propagate for _p < 0, where _p 1 –T _p/T _p is the temperature anisotropy of the proton component with respect to the interplanetary magnetic field. The damping or growth of RH cyclotron waves is found to depend on the frequency range and the temperature anisotropy of the proton component. The RH cyclotron waves are damped in the frequency range _r | _p | _p for _p < 0, where _p is the proton cyclotron frequency. RH cyclotron instabilities occur in the frequency range | _p | _p > _r > | _p | _p /(1– _r ) for _p < 0. The marginal state is at _r =| _p | _p .Abstract presented at theInternational Symposium on Solar-Terrestrial, São Paulo, Brazil, 17–22 June, 1974     

Stability of outer planetary orbits around binary stars: A comparison of Hill's and Laplace's stability criteria

A comparison is made between the stability criteria of Hill and that of Laplace to determine the stability of outer planetary orbits encircling binary stars. The restricted, analytically determined results of Hill's method by Szebehely and co-workers and the general, numerically integrated results of Laplace's method by Graziani and Black are compared for varying values of the mass parameter =m ₂/(m ₁+m ₂). For 00.15, the closest orbit (lower limit of radius) an outer planet in a binary system can have and still remain stable is determined by Hill's stability criterion. For >0.15, the critical radius is determined by Laplace's stability criterion. It appears that the Graziani-Black stability criterion describes the critical orbit within a few percent for all values of .     

Rapidly rotating polytropes in the post-Newtonian approximation to general relativity

The study of uniformly polytropes with axial symmetry is extended to include all rotational terms of order ⁴, where is the angular velocity, consistently within the first post-Newtonian approximation to general relativity. The equilibrium structure is determined by treating the effects of rotation and post-Newtonian gravitation as independent perturbations on the classical polytropic structure. The perturbation effects are characterized by a rotation parameter = ²/2G _c and a relativity parameter, =p _c / _c C ² , wherep _c and _c are the central pressure and density respectively. The solution to the structural problem is obtained by following Chandrasekhar's series expansion technique and is complete to the post-Newtonian rotation terms of order ². The critical rotation parameterv _c , which characterizes the configuration with maximum uniform rotation, is accurately evaluated as a function of . Numerical values for all the structural parameters needed to determine the equilibrium configurations are presented for polytropes with indicesn=1, 1.5, 2, 2 5, 3, and 3.5.     

Free-convection effects on MHD flow past a porous plate with step function change in suction velocity

Free convection effects on MHD flow past a semi infinite porous flat plate is studied when the time dependent suction velocity changes in step function form. The solution of the problem is obtained in closed form for the fluid with unit Prandtl number. It is observed that for both cooling and heating of the plate the suction velocity enhances the velocity field. The heat transfer is higher with increase in suction velocity.Notations B intensity of magnetic field - G Grashof number - H magnetic field parameter,H=(M+1/4) ^1/2–1/2 - M magnetic field parameter - N _u Nusselt number - P Prandtl number of the fluid - r suction parameter - T temperature of the fluid - T _w temperature of the plate - T temperature of the fluid at infinity - t time - t non-dimensional time - u velocity of the fluid parallel to the plate - u non-dimensional velocity - U velocity of the free stream - suction velocity - ₁ suction velocity att0 - ₂ suction velocity att>0 - x,y coordinate axes parallel and normal to the plate, respectively - y non-dimensional distance normal to the plate - coefficient of volume expansion - thermal diffusivity - kinematic viscosity - electric conductivity of the fluid - density of the fluid - non-dimensional temperature of the fluid - shear stress at the plate - non dimensional shear stress - erf error function - erfc complementary error function