Dissipative cylindrical fast magnetoacoustic waves in planetary magnetospheres

Nonlinear cylindrical fast magnetoacoustic waves are investigated in a dissipative magnetoplasma comprising of electrons, positrons, and ions. In this regard, cylindrical Kadomtsev-Petviashvili-Burgers (CKPB) equation is derived using the small amplitude perturbation expansion method. Furthermore, cylindrical Burgers-Kadomtsev-Petviashvili (Cyl Burgers-KP) for a fast magnetoacoustic wave is derived, for the first time, for spatial scales larger than the electron/positron skin depths, c/ω _p(e,p). Using the tangent hyperbolic method, the solutions of both planar KPB and Burgers-KP equations are obtained and then subsequently used as an initial profile to solve their respective counterparts in the cylindrical geometry. The effect of positron concentration, kinematic viscosity, and plasma β are explored both for the KPB and the Burgers-KP shock waves and the differences between the two are highlighted. The temporal evolution of the cylindrical fast magnetoacoustic wave is also numerically investigated. The present study may be beneficial to study the propagation characteristics of nonlinear electromagnetic shock waves in planetary magnetospheres.     

High-resolution spectroscopy of QY Sge: an obscured RV Tauri variable?

We describe and demonstrate a method for increasing the resolution locally in a Smoothed Particle Hydrodynamic (SPH) simulation, by splitting particles. We show that in simulations of self-gravitating collapse (of the sort that are presumed to occur in star formation) the method is stable, and affords great savings in computer time and memory. When applied to the standard Boss & Bodenheimer test – which has been shown to depend critically on fulfilment of the Jeans Condition – the results are comparable both with those obtained using Adaptive Mesh Refinement, and with those obtained using a standard high-resolution SPH simulation, but they are achieved with considerably less computational resource. Further development and testing is required before the method can safely be applied to more general flows.     

A simulation study of the solar wind including the solar rotation effect

An axisymmetric solar wind structure including the solar rotation effect is studied by the method of MHD computer simulation. For the case of the radial magnetic field configuration, the simulation result is fairly well coincident with the steady-state solution. For the case of the dipole magnetic field configuration, the properties of the solution depend on the ratio of the gas pressure to the magnetic pressure-ratio) in the model. If the-ratio is small, a clearly defined stagnation region appears in the wind, in which the flow speed is very small and the azimuthal magnetic field is very weak because of the corotation of the plasma. If the-ratio is greater than 1, the plasma is not effectively trapped by the magnetic field so that the stagnation region is not clearly defined in the solution.     

A special perturbation method in orbital dynamics

The special perturbation method considered in this paper combines simplicity of computer implementation, speed and precision, and can propagate the orbit of any material particle. The paper describes the evolution of some orbital elements based in Euler parameters, which are constants in the unperturbed problem, but which evolve in the time scale imposed by the perturbation. The variation of parameters technique is used to develop expressions for the derivatives of seven elements for the general case, which includes any type of perturbation. These basic differential equations are slightly modified by introducing one additional equation for the time, reaching a total order of eight. The method was developed in the Grupo de Dinámica de Tethers (GDT) of the UPM, as a tool for dynamic simulations of tethers. However, it can be used in any other field and with any kind of orbit and perturbation. It is free of singularities related to small inclination and/or eccentricity. The use of Euler parameters makes it robust. The perturbation forces are handled in a very simple way: the method requires their components in the orbital frame or in an inertial frame. A comparison with other schemes is performed in the paper to show the good performance of the method.     

Implications of Non-cylindrical Flux Ropes for Magnetic Cloud Reconstruction Techniques and the Interpretation of Double Flux Rope Events

Magnetic clouds (MCs) are a subset of interplanetary coronal mass ejections (ICMEs) which exhibit signatures consistent with a magnetic flux rope structure. Techniques for reconstructing flux rope orientation from single-point in situ observations typically assume the flux rope is locally cylindrical, e.g., minimum variance analysis (MVA) and force-free flux rope (FFFR) fitting. In this study, we outline a non-cylindrical magnetic flux rope model, in which the flux rope radius and axial curvature can both vary along the length of the axis. This model is not necessarily intended to represent the global structure of MCs, but it can be used to quantify the error in MC reconstruction resulting from the cylindrical approximation. When the local flux rope axis is approximately perpendicular to the heliocentric radial direction, which is also the effective spacecraft trajectory through a magnetic cloud, the error in using cylindrical reconstruction methods is relatively small (≈ 10^∘). However, as the local axis orientation becomes increasingly aligned with the radial direction, the spacecraft trajectory may pass close to the axis at two separate locations. This results in a magnetic field time series which deviates significantly from encounters with a force-free flux rope, and consequently the error in the axis orientation derived from cylindrical reconstructions can be as much as 90^∘. Such two-axis encounters can result in an apparent ‘double flux rope’ signature in the magnetic field time series, sometimes observed in spacecraft data. Analysing each axis encounter independently produces reasonably accurate axis orientations with MVA, but larger errors with FFFR fitting.     

An implementation ofN-body chain regularization

The chain regularization method (Mikkola and Aarseth 1990) for high accuracy computation of particle motions in smallN-body systems has been reformulated. We discuss the transformation formulae, equations of motion and selection of a chain of interparticle vectors such that the critical interactions requiring regularization are included in the chain. The Kustaanheimo-Stiefel (KS) coordinate transformation and a time transformation is used to regularize the dominant terms of the equations of motion. The method has been implemented for an arbitrary number of bodies, with the option of external perturbations. This formulation has been succesfully tested in a generalN-body program for strongly interacting subsystems. An easy to use computer program, written inFortran, is available on request.     

Galaxies in the connection machine

Many problems in galactic dynamics require computer speeds orders of magnitude larger than what can be achieved on current single-processor machines. In the near future such speeds are likely to become available through computer architectures based on large-scale, fine-grained parallelism. An example of a highly parallel computer is the Connection Machine, with up to 65,636 processors. We have benchmarked gravitationalN-body algorithms on the Connection Machine, and compared those with similar benchmarks which we have obtained on more traditional vector supercomputers. Our conclusions are: (1) The direct summation algorithm, with of orderN ² interactions forN particles, can be made to run with high efficiency on either type of computer. As a result, the Connection Machine clearly wins in speed over all supercomuters tested, with the exception of an 8-processor ETA, which shows a comparable performance. (2) A more efficient tree algorithm reduces the growth of the number of interactions fromN ² toN logN. However, the greater complexity of this algorithm causes a considerable degradation of efficiency, by a factor which is larger on the Connection Machine than on vector supercomputers. As a result, out tree code runs at comparable speeds on both types of machines, with the notable exception of the 8-processor ETA, which has an extrapolated speed for running our tree code which is higher than any of the other machines we have tested.     

Three-Dimensional Particle Acceleration in Electromagnetic Cylinder and Torus

Particle acceleration via Poynting vector with toroidal magnetic field is studied in 3D PIC simulation of electron-positron plasma. We choose two different initial magnetic field configurations to compare how the particle acceleration is affected by the expansion of electromagnetic wave. In the cylindrical case, the electromagnetic field strength decays as (ct)⁻², and particles are accelerated in the radial direction as well as the axial direction. Rayleigh-Taylor instability is also observed at the center of the cylinder. In the torus case, the field strength decays as (ct)⁻³, making the acceleration less efficient. Particles accelerated in the axial direction by E × B force creates strong charge separation.     

Practical considerations in the numerical solution of theN-body problem by Taylor series methods

It is shown that in the numerical integration ofN-body problems, as much importance should be given to considerations of the computer programming language to be used as to questions of the accummulation of round-off and truncation error, the stability of the method chosen and the problem being treated. By careful programming processing time may be cut by a factor of 2 or 3 which is an important consideration in extended numerical investigations. The relative usefulness of differing strategies for determining the step size is discussed and in addition the usefulness is shown of treatingN-body problems by a Taylor series method.     

Sweeping-magnetic-twist mechanism for the acceleration of jets in the solar atmosphere

A magnetodynamic mechanism for the acceleration of jets in the solar atmosphere (surges, Brueckner's EUV jets, and so on) is proposed, and a 2.5-dimensional MHD simulation is performed to show how this mechanism operates in the situation of the chromosphere-corona region of the solar atmosphere. It is seen from the result of simulation that together with the release of the magnetic twist, e.g., into a reconnected open flux tube, the mass in the high density twisted loop is driven out into the open flux tube due both to the pinch effect progressing with the packet of the magnetic twist into the open flux tube, and to the j × B force at the front of the packet of the unwinding twist in the off-axis part of the tube. The former, the progressing pinch, is accompanied by an accelerated hot blob, while the latter, the unwinding front of the magnetic twist, drives a cool cylindrical flow, both with velocities of the order of the local Alfvén velocity. One of the characteristic properties of the jet in our model is that the jet, consisting of hot core and cool sheath, has a helical velocity field in it, explaining the thus-far unexplained observed feature.The sudden release of the magnetic twist into an open flux tube is most likely to be due to the reconnection between a twisted loop and the open flux tube. The mass is driven out in the relaxation process of the magnetic twist from the twisted loop to the open flux tube.     

Numerical simulation of reconnection in an emerging magnetic flux region

The resistive MHD equations are numerically solved in two dimensions for an initial-boundary-value problem which simulates reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. The emerging region is modelled by a cylindrical flux tube with a poloidal magnetic field lying in the same plane as the external, coronal field. The plasma betas of the emerging and coronal regions are 1.0 and 0.1, respectively, and the magnetic Reynolds number for the system is 2 × 10³. At the beginning of the simulation the tube starts to emerge through the base of the rectangular computational domain, and, when the tube is halfway into the computational domain, its position is held fixed so that no more flux of plasma enters through the base. Because the time-scale of the emergence is slower than the Alfvén time-scale, but faster than the reconnection time-scale, a region of closed loops forms at the base. These loops are gradually opened and reconnected with the overlying, external magnetic field as time proceeds.The evolution of the plasma can be divided into four phases as follows: First, an initial, quasi-steady phase during which most of the emergence is completed. During this phase, reconnection initially occurs at the slow rate predicted by the Sweet model of diffusive reconnection, but increases steadily until the fast rate predicted by the Petschek model of slow-shock reconnection is approached. Second, an impulsive phase with large-scale, super-magnetosonic flows. This phase appears to be triggered when the internal mechanical equilibrium inside the emerging flux tube is upset by reconnection acting on the outer layers of the flux tube. During the impulsive phase most of the flux tube pinches off from the base to form a cylindrical magnetic island, and temporarily the reconnection rate exceeds the steady-state Petschek rate. (At the time of the peak reconnection rate, the diffusion region at the X-line is not fully resolved, and so this may be a numerical artifact.) Third, a second quasi-steady phase during which the magnetic island created in the impulsive phase is slowly dissipated by continuing, but low-level, reconnection. And fourth, a static, non-evolving phase containing a potential, current-free field and virtually no flow.During the short time in the impulsive phase when the reconnection rate exceeds the steady-state Petschek rate, a pile-up of magnetic flux at the neutral line occurs. At the same time the existing Petschek-slow-mode shocks are shed and replaced by new ones; and, for a while, both new and old sets of slow shocks coexist.     

Time evolution of nonplanar electron acoustic shock waves in a plasma with superthermal electrons

The propagation of cylindrical and spherical electron acoustic (EA) shock waves in unmagnetized plasmas consisting of cold fluid electrons, hot electrons obeying a superthermal distribution and stationary ions, has been investigated. The standard reductive perturbation method (RPM) has been employed to derive the cylindrical/spherical Korteweg-de-Vries-Burger (KdVB) equation which governs the dynamics of the EA shock structures. The effects of nonplanar geometry, plasma kinematic viscosity and electron suprathermality on the temporal evolution of the cylindrical and spherical EA shock waves are numerically examined.     

Computational approach to the radio-emission of pulsars

Radio-emission of pulsars is investigated through computer simulation. We assume that electronpositron plasma is thermal and beam particles have _b in the plasma rest frame. The main purpose of this study is to search for unstable electromagnetic waves which may be caused by beam-plasma interactions. We obtain the dominant waves in the range of radio frequencies in the observer frame. The results show that transverse electromagnetic waves grow rapidly, while electrostatic waves have very slow growth.     

A large-eddy simulation of turbulent compressible convection: differential rotation in the solar convection zone

We present the results of two simulations of the convection zone, obtained by solving the full hydrodynamic equations in a section of a spherical shell. The first simulation has cylindrical rotation contours (parallel to the rotation axis) and a strong meridional circulation, which traverses the entire depth. The second simulation has isorotation contours about mid-way between cylinders and cones, and a weak meridional circulation, concentrated in the uppermost part of the shell.
We show that the solar differential rotation is directly related to a latitudinal entropy gradient, which pervades into the deep layers of the convection zone. We also offer an explanation of the angular velocity shear found at low latitudes near the top. A non-zero correlation between radial and zonal velocity fluctuations produces a significant Reynolds stress in that region. This constitutes a net transport of angular momentum inwards, which causes a slight modification of the overall structure of the differential rotation near the top. In essence, the thermodynamics controls the dynamics through the Taylor–Proudman momentum balance . The Reynolds stresses only become significant in the surface layers, where they generate a weak meridional circulation and an angular velocity 'bump'.     

Photoelectric drift scans II: time pulse evaluation,limb profiles,and the solar diameter

The method of computer controlled photoelectric drift scans has been improved by virtue of a new timing technique, allowing for an accuracy of 4 ms or equivalently 50 km on the Sun. With this technique, our previous result for the solar semidiameter has been confirmed: R = 960.277″. The extreme solar limb intensity profile at 5012 Å has been derived from drift scans. A comparison with the computed profile has been made, and parameters for the base of the chromosphere have been derived.

     

Computer simulation of chemical nucleation

The problem of nucleation at chemical instabilities is investigated by means of microscopic computer simulation. The first-order transition of interest involves a new kind of nucleation arising from chemical transformations rather than physical forces. Here it is the chemical state of matter, and not matter itself, which is spatially localized to form the nucleus for transition between different chemical states. First, the concepts of chemical instability, nonequilibrium phase transition, and dissipative structure are reviewed briefly. Then recently developed methods of reactive molecular dynamics are used to study chemical nucleation in a simple model chemical reaction. Finally, the connection of these studies to nucleation and condensation processes involving physical and chemical interactions is explored.Invited contribution to the Proceedings of a Workshop onThermodynamics and Kinetics of Dust Formation in the Space Medium held at the Lunar and Planetary Institute, Houston, 6–8 September, 1978.     

基于CPLD的高精度可调脉冲信号发生器研制

为满足精密时间间隔测量设备的测试需要,研制了一种时间间隔可调的高精度脉冲信号发生器。利用计算机串口控制的方式,结合复杂可编程逻辑器件（CPLD）集成度高、可靠性好及工作速度快的优点,采用A1tera公司的设计软件QuartusII进行设计仿真及实现。仿真与实测实验表明,该脉冲信号发生器不仅可以产生单路可调脉冲信号,而且能产生多路可调脉冲信号,产生的单路秒脉冲信号的1s取样Allan方差为1．84×10^-11;产生的时间间隔为100ns的多路脉冲信号的1s取样Allan方差为2．36×10^-11,2路信号之间的时间间隔数据系列的峰一峰值为101ps,可以满足多通道时间间隔测量设备测试要求的稳定度与准确度。     

Analytical solutions of cylindrical shock waves in magnetogasdynamics,I

The C.C.W. method has been used to investigate the propagation of converging and diverging cylindrical shock waves in a non-uniform medium under the influence of a magnetic field of constant strength. A comparison has also been made between the two types of cylindrical shock waves, simultaneously for both weak and strong cases of the magnetic field. Density distribution is assumed to be _o = r ^–, where is the density at the axis of symmetry and a constant. The analytical expressions for shock velocity and shock strength as well as the pressure, the density, and the particle velocity just behind the shock front have been derived for both the cases.     

Large mirror figuring and testing

We describe a method for figuring and testing large aspheric mirrors using a rectangular, flexible lap (the so-called Membrane Tool) and a vibration stabilized interferometer. The rear side of the lap is covered with computer controlled dynamic pressure actuators which determine the amount of material to be removed for surface error correction. This method has been developed in the laboratory and tested to some extent by figuring thef/2.2 primary of the 3.6 m ESO-NTT. We describe the ongoing developments and the manufacturing plan for 8 m-class mirrors.Paper presented at the Symposium on the JNLT and Related Engineering Developments, Tokyo, November 29–December 2, 1988.     

天体物理学中动力学研究的计算机模拟

简要介绍了计算机模拟的基本原理和方法 ,以及数值模拟计算在天体物理学中的主要应用结果和前景