磁流体冲激波在非均匀磁场中的传播

本文应用Chisnell- no方法,求解了在理想介貭中,垂直磁流体冲激波在非均勻磁場中的传播問題。这种方法,把非均勻介貭分解成无限小的弱間断面,根据气体动力学中波与間断面相互作用的原理,算出激波通过弱間断面时的强度变化,然后用积分求得激波通过整个非均勻区时的强度变化。作者引入了激波的特征速度（它是激波在波前后介貭中传播速度的几何平均值）作为輔助参量,得到形式上比較簡单的激波传播方程。然后考虑了磁压力远大于气体压力的强磁介貭中的激波传播問題,并进行了数值积分。采用的介貭密度模型有三种:（1）阿尔芬波速为常数;（2）密度不变;（3）密度与磁場强度成正比。計算結果表明:当激波由弱磁場向强磁場传播时,激波的强度逐漸变弱。其中,在阿尔芬波速为常数的介貭中,激波强度的衰减最为緩慢;在密度不变的介貭中,激波强度的衰減最为迅速;而在密度与磁場成正比的介貭中,激波强度的衰減則介乎上述两种密度分布之間。作者联系磁流体冲激波在地球外层空間的传播問題进行了討論,密度的模型采取大气啃昔的观測結果（卽上述第三种密度分布）,并进行了适当的外推,估計了在十个地球半径处的磁流体冲激波传到地面时的强度,求出了激波在地面引起的磁場变化与激波初始速度之間的关系。根据上述簡化模型,計算結果表明,在十个地球半径处初始速度为10⁸厘米/秒的激波,传到地面引起的磁場变化大約为60伽（亻馬）,这个数值的量級恰好与中低緯度强磁暴的急始变幅相符。     

Radiative magnetogasdynamic shock in a self-gravitating system

Summary By taking radiative parameters independent of magnetic field, the aim of this paper is to investigate the propagation of spherical shock in self-gravitating conducting gas spheres, such as stars, caused by an instantaneous central explosion of finite amount of energy. Ionisation, dissociation have not been taken into account and the medium is assumed to be a fully perfect gas.     

Hydromagnetic waves in a thin rotating spherical shell

We consider an electrically conducting fluid confined to a thin rotating spherical shell in which the Elsasser and magnetic Reynolds numbers are assumed to be large while the Rossby number is assumed to vanish in an appropriate limit. This may be taken as a simple model for a possible stable layer at the top of the Earth's outer core. It may also be a model for the thin shells which are thought to be a source of the magnetic fields of some planets such as Mercury or Uranus. Linear hydromagnetic waves are studied using a multiple scale asymptotic scheme in which boundary layers and the associated boundary conditions determine the structure of the waves. These waves are assumed to be of the form of an asymptotic series expanded about an ambient magnetic field which vanishes on the equatorial plane and velocity and pressure fields which do not. They take the form of short wave, slowly varying wave trains. The results are compared to the author's previous work on such waves in cylindrical geometry in which the boundary conditions play no role. The approximation obtained is significantly different from that obtained in the previous work in that an essential singularity appears at the equator and nonequatorial wave regions appear.     

On the disturbance due to a point source in an inhomogeneous medium

Summary The problem of a point source in an isotropic, inhomogeneous fluid medium is discussed. It is assumed that the density of the fluid is constant and the acoustic velocity varies with depth asc=c ₀(1 +m z) wherem is a constant andc ₀ is, the velocity at the level of the origin. An approximate expression for the field due to a point source in such a medium is obtained when the medium is infinite as well as when it is semi-infinite. It is found that the results obtained agree with the WKB solution of the problem.     

Electron velocity distribution function in a plasma with temperature gradient and in the presence of suprathermal electrons: application to incoherent-scatter plasma lines

The plasma dispersion function and the reduced velocity distribution function are calculated numerically for any arbitrary velocity distribution function with cylindrical symmetry along the magnetic field. The electron velocity distribution is separated into two distributions representing the distribution of the ambient electrons and the suprathermal electrons. The velocity distribution function of the ambient electrons is modelled by a near-Maxwellian distribution function in presence of a temperature gradient and a potential electric field. The velocity distribution function of the suprathermal electrons is derived from a numerical model of the angular energy flux spectrum obtained by solving the transport equation of electrons. The numerical method used to calculate the plasma dispersion function and the reduced velocity distribution is described. The numerical code is used with simulated data to evaluate the Doppler frequency asymmetry between the up- and downshifted plasma lines of the incoherent-scatter plasma lines at different wave vectors. It is shown that the observed Doppler asymmetry is more dependent on deviation from the Maxwellian through the thermal part for high-frequency radars, while for low-frequency radars the Doppler asymmetry depends more on the presence of a suprathermal population. It is also seen that the full evaluation of the plasma dispersion function gives larger Doppler asymmetry than the heat flow approximation for Langmuir waves with phase velocity about three to six times the mean thermal velocity. For such waves the moment expansion of the dispersion function is not fully valid and the full calculation of the dispersion function is needed.     

On the propagation of a reflected radiative shock in a transverse magnetic field

Summary The paper considers a shock wave to be reflected from a wall in the presence of a transverse magnetic field when the radiation pressure and radiation energy are taken into account. It has been shown that under certain physical conditions, the velocity of the reflected radiative shock would be greater than the corresponding velocity when the magnetic field and radiation effects are neglected.     

地球主磁场电磁转矩导致的地球差异旋转

为探索驱动地球系统差异旋转的力源,选择整个地球作为研究对象,应用经典电磁学理论,分析地球固体介质中的电荷在地球主磁场中的运动特点,发现存在一个与地球自转方向相反的切向洛伦兹力.通过电荷与介质间的相互作用,切向洛伦兹力传递至介质迫使介质西漂.为探索地球差异旋转的规律,建立了地球薄圆筒圈层模型.应用经典物理学理论和方法推导出了地球主磁场电磁转矩及其产生的角加速度公式.研究得到四点主要结论:1)作用于半径不同的地球薄圆筒圈层的地球主磁场电磁转矩及其产生的角加速度绝对值不同:地轴及赤道附近圈层的小,其自转相对较快;半径等于3～(1/2)倍地球半径的薄圆筒圈层及其相邻圈层的大,其自转相对较慢.2)同一薄圆筒圈层中的差异旋转缘于介质的介电常数、阻力系数及质量密度的差异.3)地球差异旋转缘于地球的自转、正负电荷的非对等分布及介质的介电常数、阻力系数、质量密度的差异.4)地球差异旋转导致地壳运动,孕育地震,地球主磁场是地球差异旋转和地震孕育的敏感因子.     

Galactic dynamo models without sharp boundaries

Abstract

A new numerical approach is introduced which allows investigation into the conditions for dynamogeneration of axisymmetric and non-axisymmetric large-scale magnetic field modes in galaxy models which are defined by axisymmetric distributions of the α-parameter, the angular velocity and the electrical conductivity. The velocity field is assumed to be localized, however, the common assumption of a sharp boundary of the conducting region is dropped.

The possible anisotropy of the α-tensor is taken into account. The critical dynamo numbers (excitation conditions) for different modes are obtained by a direct method. The required steady states are attained by the use of an artificial non-linearity.

Initial test calculations demonstrate the efficacy of this new concept.     

TRM in low magnetic fields: a minimum field that can be recorded by large multidomain grains

Thermally acquired remanent magnetization is important for the estimation of the past magnetic field present at the time of cooling. Rocks that cool slowly commonly contain magnetic grains of millimeter scale. This study investigated 1-mm-sized magnetic minerals of iron, iron–nickel, magnetite, and hematite and concluded that the thermoremanent magnetization (TRM) acquired by these grains did not accurately record the ambient magnetic fields less than 1 μT. Instead, the TRM of these grains fluctuated around a constant value. Consequently, the magnetic grain ability to record the ambient field accurately is reduced. Above the critical field, TRM acquisition is governed by an empirical law and is proportional to saturation magnetization (M_s). The efficiency of TRM is inversely proportional to the mineral's saturation magnetization M_s and is related to the number of domains in the magnetic grains. The absolute field for which we have an onset of TRM sensitivity is inversely proportional to the size of the magnetic grain. These results have implications for previous reports of random directions in meteorites during alternating field demagnetization, or thermal demagnetization of TRM. Extraterrestrial magnetic fields in our solar system are weaker than the geomagnetic field by several orders of magnitude. Extraterrestrial rocks commonly contain large iron-based magnetic minerals as a common part of their composition, and therefore ignoring this behavior of multidomain grains can result in erroneous paleofield estimates.     

Unsteady hydromagnetic flow in an annulus with porous Walls when the inner boundary starts moving from rest

Summary The problem of unsteady hydromagnetic flow in a porous annulus when the inner boundary starts moving from rest with a constant longitudinal velocity is considered. For the case of equal Reynolds number and magnetic Reynolds number Laplace transforms of the solutions for the axial components of velocity and magnetic field are obtained in terms of modified Bessel functions. The inversions are effected for the cases of very large and small Reynolds number when the pressure gradient vanishes.     

The self-similar,non-linear evolution of rotating magnetic flux ropes

We study, in the ideal MHD approximation, the non-linear evolution of cylindrical magnetic flux tubes differentially rotating about their symmetry axis. Our force balance consists of inertial terms, which include the centrifugal force, the gradient of the axial magnetic pressure, the magnetic pinch force and the gradient of the gas pressure. We employ the “separable” class of self-similar magnetic fields, defined recently. Taking the gas to be a polytrope, we reduce the problem to a single, ordinary differential equation for the evolution function. In general, two regimes of evolution are possible; expansion and oscillation. We investigate the specific effect rotation has on these two modes of evolution. We focus on critical values of the flux rope parameters and show that rotation can suppress the oscillatory mode. We estimate the critical value of the angular velocity _crit, above which the magnetic flux rope always expands, regardless of the value of the initial energy. Studying small-amplitude oscillations of the rope, we find that torsional oscillations are superimposed on the rotation and that they have a frequency equal to that of the radial oscillations. By setting the axial component of the magnetic field to zero, we study small-amplitude oscillations of a rigidly rotating pinch. We find that the frequency of oscillation is inversely proportional to the angular velocity of rotation ; the product being proportional to the inverse square of the Alfvén time. The period of large-amplitude oscillations of a rotating flux rope of low beta increases exponentially with the energy of the equivalent 1D oscillator. With respect to large-amplitude oscillations of a non-rotating flux rope, the only change brought about by rotation is to introduce a multiplicative factor greater than unity, which further increases the period. This multiplicative factor depends on the ratio of the azimuthal speed to the Alfvén speed. Finally, considering interplanetary magnetic clouds as cylindrical flux ropes, we inquire whether they rotate. We find that at 1 AU only a minority do. We discuss data on two magnetic clouds where we interpret the presence in each of vortical plasma motion about the symmetry axis as a sign of rotation. Our estimates for the angular velocities suggest that the parameters of the two magnetic clouds are below critical values. The two clouds differ in many respects (such as age, bulk flow speed, size, handedness of the magnetic field, etc.), and we find that their rotational parameters reflect some of these differences, particularly the difference in age. In both clouds, a rough estimate of the radial electric field in the rigidly rotating core, calculated in a non-rotating frame, yields values of the order mV m⁻¹.     

Nonlinear waves in rotating magnetohydrodynamics

We consider an electrically conducting fluid in rotating cylindrical coordinates in which the Elsasser and magnetic Reynolds numbers are assumed to be large while the Rossby number is assumed to vanish in an appropriate limit. This may be taken as a simple model for the Earth's outer core. Fully nonlinear waves dominated by the nonlinear Lorentz forces are studied using the method of geometric optics (essentially WKB). These waves are assumed to be of the form of an asymptotic series expanded about ambient magnetic and velocity fields which vanish on the equatorial plane. They take the form of short wave, slowly varying wave trains. The first-order approximation is sinusoidal and basically the same as in the linear problem, with a dispersion relation modified by the appearance of mean terms. These mean terms, as well the undetermined amplitude functions, are found by suppressing secular terms in a “fast” variable in the second-order approximation. The interaction of the mean terms with the dispersion relation is the primary cause of behaviors which differ from the linear case. In particular, new singularities appear in the wave amplitude functions and an initial value problem results in a singularity in one of the mean terms which propagates through the fluid. The singularities corresponding to the linear ones are shown to develop when the corresponding waves propagate toward the equatorial plane.     

Dynamic response of a partially debonded pipeline embedded in a saturated poroelastic medium to harmonic plane waves

A practical model of partially debonded pipeline embedded in a saturated poroelastic medium is proposed, and the dynamic response of this model to harmonic plane waves is theoretical investigated. Biot׳s poroelastic theory is introduced to describe the dynamic equations of the saturated poroelastic medium, and the potentials obtained from Helmholtz decomposition theorem are expressed by wave function expansion method. The debonding areas around the pipeline are assumed to be filled with water. The disturbed solutions of basic field equations in these areas are expressed in terms of a scalar velocity potential. Different boundary conditions of bonded and debonded areas are adopted, and the expanded coefficients are obtained. An example of one partially debonded area is presented and analyzed. It is found that the stresses in the perfectly bonded and debonded areas show great difference, and the jump of dynamic stress at the connection points between these two areas is great in the case of low frequency. The effect of debonded areas on the dynamic stress under different thicknesses of lining is also examined.     

Surfatron particle acceleration in the near-Earth shock wave

A theoretical model, within the scope of which the process of surfatron particle acceleration (surfing) at the shock front was studied, is proposed to describe the near-Earth shock wave structure. The most complex shock front geometry, when it is assumed that the magnetic field vectors and the velocity of the plasma flow incident on the front are oriented arbitrarily with respect to the front plane, has been considered in detail. The theoretical ultimate energies of protons, accelerated due to surfing at the shock front, agree with the observational data.     

On the dissipation of atmospheric alfvén waves in uniform and non-uniform magnetic fields

Abstract

We deduce the dissipative Alfvén wave equation in a medium stratified in one direction, with a transverse magnetic field, in the presence of dissipation by fluid viscosity and electrical resistance; the dissipative Alfvén wave equation generalizes earlier results for homogeneous (Cowling, 1960) and inhomogeneous (Campos, 1983a) media, and corrects an error in the literature (Heyvaerts and Priest, 1983). The wave equation is solved exactly in two cases: a uniform magnetic field, and a magnetic field decreasing with height. In both cases the mean state is assumed to be isothermal, with a constant rate of ionization, so that the magnetic diffusivity is constant, but the dynamic viscosity increases with height. There are therefore two regions, a low- (high-) altitude region where electrical resistance dominates fluid viscosity (or vice versa), and an asymptotic regime relevant to the uppermost (lowermost) layers. The two regions are separated by a transition layer, across which the wave field is continuous and whose structure is expressible by hypergeometric functions, with different arguments in the low- and high-altitude regions, and over the whole altitude range. These exact solutions allow the amplitude and phase of the wave field to be plotted as a function of height for a variety of magnetoatmospheric mean states. They show that wave dissipation is more localized and intense when the magnetic field decreases with height than when it is uniform.     

Modelling of GPR waves for lossy media obeying a complex power law of frequency for dielectric permittivity

The attenuation of ground‐penetrating radar (GPR) energy in the subsurface decreases and shifts the amplitude spectrum of the radar pulse to lower frequencies (absorption) with increasing traveltime and causes also a distortion of wavelet phase (dispersion). The attenuation is often expressed by the quality factor Q. For GPR studies, Q can be estimated from the ratio of the real part to the imaginary part of the dielectric permittivity. We consider a complex power function of frequency for the dielectric permittivity, and show that this dielectric response corresponds to a frequency‐independent‐Q or simply a constant‐Q model. The phase velocity (dispersion relationship) and the absorption coefficient of electromagnetic waves also obey a frequency power law. This approach is easy to use in the frequency domain and the wave propagation can be described by two parameters only, for example Q and the phase velocity at an arbitrary reference frequency. This simplicity makes it practical for any inversion technique. Furthermore, by using the Hilbert transform relating the velocity and the absorption coefficient (which obeys a frequency power law), we find the same dispersion relationship for the phase velocity. Both approaches are valid for a constant value of Q over a restricted frequency‐bandwidth, and are applicable in a material that is assumed to have no instantaneous dielectric response. Many GPR profiles acquired in a dry aeolian environment have shown a strong reflectivity inside dunes. Changes in water content are believed to be the origin of this reflectivity. We model the radar reflections from the bottom of a dry aeolian dune using the 1D wavelet modelling method. We discuss the choice of the reference wavelet in this modelling approach. A trial‐and‐error match of modelled and observed data was performed to estimate the optimum set of parameters characterizing the materials composing the site. Additionally, by combining the complex refractive index method (CRIM) and/or Topp equations for the bulk permittivity (dielectric constant) of moist sandy soils with a frequency power law for the dielectric response, we introduce them into the expression for the reflection coefficient. Using this method, we can estimate the water content and explain its effect on the reflection coefficient and on wavelet modelling.     

On unsteady motion of a rotating fluid bounded by flat porous plates

Summary In this paper, an investigation is made of the unsteady flow generated in a viscous, incompressible and homogeneous fluid bounded by (i) an infinite horizontal porous plate atz=0, or (ii) two infinite horizontal porous plates atz=0 andz=D. The fluid together with the plate(s) is in a state of solid body rotation with a constant angular velocity about the z-axis normal to the plate(s), and additionally, the plate(s) performs non-torsional elliptic harmonic oscillations in its (their) own plane(s). A uniform suction or injection is introduced in the configurations through the porous plate(s) and its influence on the unsteady flow and the associated boundary layers is examined. The unsteady flow field as well as the associated boundary layers is obtained explicitly. In contrast to the unsteady rotating flow without suction, solutions of the present problem with suction exhibit no resonant phenomena. It is shown that the suction is responsible for making the boundary layers thinner and for the elimination of the resonant phenomena. It is confirmed that the velocity field and the associated multiple boundary layers are significantly modified by suction. Physical significances of the mathematical results are discussed. Several limiting cases of interest are recovered from this analysis. The initial value problem for both the configurations is exactly solved by the Heaviside operational calculus combined with the theory of residues.     

Linear and angular momentum conservation in hydraulic jump in diverging channels

This paper addresses the integral conservation of linear and angular momentum in the steady hydraulic jump in a linearly diverging channel.The flow is considered to be divided into a mainstream that conveys the total liquid discharge, and a roller where no average mass transport occurs. It is assumed that no macroscopic rheological relationship holds, so mass, momentum and angular momentum integral balances are independent relationships. Normal stresses are assumed to be hydrostatic on vertical, cylindrical surfaces. Viscous stresses are assumed to be negligible with respect to turbulent stresses. Assuming that the horizontal velocity distribution in the mainstream is uniform and that the horizontal momentum and angular momentum in the roller are negligible with respect to their mainstream counterparts, an analytical solution is obtained for the free surface profile of the flow. This solution is fundamental for finding the sequent depths and their positions. Consequently, it permits solving for the length of the jump, which is assumed to be equal to the length of the roller. Mainstream and roller thicknesses can also be derived from the present solution. This model may also be theoretically used to derive the average shear stresses exerted by the roller on the mainstream and the power losses per unit weight. This second relationship, which returns the well-known classical expression for total power loss in the jump, demonstrates that the strongly idealized mechanical model proposed here is internally consistent.     

Turbulent transport of angular momentum in magnetized fluids and torsional oscillations of the sun

Abstract

Fluxes of angular momentum produced by turbulence in rotating fluids are derived with the effects of a magnetic field included. It is assumed that the rotation is slow but that the magnetic field is of arbitrary strength. A mean magnetic field is shown to produce qualitative changes of the sources of the differential rotation rather than the quenching of differential rotation usually expected. A new equatorward flux of angular momentum arises through the influence of the toroidal magnetic field. The possibility of interpreting the torsional oscillations of the Sun as a consequence of the magnetic perturbations of the turbulent angular momentum fluxes is discussed.