Simulation Techniques for Cosmological Simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.     

The Radio Plasma Imager investigation on the IMAGE spacecraft

Radio plasma imaging uses total reflection of electromagnetic waves from plasmas whose plasma frequencies equal the radio sounding frequency and whose electron density gradients are parallel to the wave normals. The Radio Plasma Imager (RPI) has two orthogonal 500-m long dipole antennas in the spin plane for near omni-directional transmission. The third antenna is a 20-m dipole along the spin axis. Echoes from the magnetopause, plasmasphere and cusp will be received with the three orthogonal antennas, allowing the determination of their angle-of-arrival. Thus it will be possible to create image fragments of the reflecting density structures. The instrument can execute a large variety of programmable measuring options at frequencies between 3 kHz and 3 MHz. Tuning of the transmit antennas provides optimum power transfer from the 10 W transmitter to the antennas. The instrument can operate in three active sounding modes: (1) remote sounding to probe magnetospheric boundaries, (2) local (relaxation) sounding to probe the local plasma frequency and scalar magnetic field, and (3) whistler stimulation sounding. In addition, there is a passive mode to record natural emissions, and to determine the local electron density, the scalar magnetic field, and temperature by using a thermal noise spectroscopy technique.     

Radio plasma imager simulations and measurements

The Radio Plasma Imager (RPI) will be the first-of-its kind instrument designed to use radio wave sounding techniques to perform repetitive remote sensing measurements of electron number density (N _e) structures and the dynamics of the magnetosphere and plasmasphere. RPI will fly on the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission to be launched early in the year 2000. The design of the RPI is based on recent advances in radio transmitter and receiver design and modern digital processing techniques perfected for ground-based ionospheric sounding over the last two decades. Free-space electromagnetic waves transmitted by the RPI located in the low-density magnetospheric cavity will be reflected at distant plasma cutoffs. The location and characteristics of the plasma at those remote reflection points can then be derived from measurements of the echo amplitude, phase, delay time, frequency, polarization, Doppler shift, and echo direction. The 500 m tip-to-tip X and Y (spin plane) antennas and 20 m Z axis antenna on RPI will be used to measures echoes coming from distances of several R _E. RPI will operate at frequencies between 3 kHz to 3 MHz and will provide quantitative N _e values from 10^–1 to 10⁵ cm^–3. Ray tracing calculations, combined with specific radio imager instrument characteristics, enables simulations of RPI measurements. These simulations have been performed throughout an IMAGE orbit and under different model magnetospheric conditions. They dramatically show that radio sounding can be used quite successfully to measure a wealth of magnetospheric phenomena such as magnetopause boundary motions and plasmapause dynamics. The radio imaging technique will provide a truly exciting opportunity to study global magnetospheric dynamics in a way that was never before possible.     

The Time Structure of Ground Level Enhancements in?Solar Cycle 23

In a recent paper McCracken et al. (J. Geophys. Res. 113:A12101, 2008) proposed that the Ground Level Enhancement (GLE) of 20 January 2005 may have been produced by more than one acceleration mechanism, with the first acceleration due to the solar flare and the second one due to the CME associated with that event. They also noted several other GLEs with similar multiple pulse structures. This paper systematically investigates all the GLEs of solar cycle 23, from GLE 55 on 6 November 1997 to GLE 70 on 13 December 2006, to study their morphology and pulse structure, and to determine whether the multiple structures that may be found in these events are qualitatively similar to that of the GLE of 20 January 2005. We use all the data of all NMs that saw each event, to have as much directional and spectral information as possible. It is shown that three of these 16 events do contain such double-pulse structures, and the properties of these three are discussed in some detail.     

Telerobotics: problems and research needs

With major emphasis on simulation, a university laboratory telerobotics facility permits problems to be approached by groups of graduate students. Helmet-mounded displays provide realism; the slaving of the display to the human operator's viewpoint gives a sense of `telepresence' that may be useful for prolonged tasks. Using top-down 3-D model control of distant images allows distant images to be reduced to a few parameters to update the model used for display to the human operator in a preview model to circumvent, in part, the communication delay. Also, the model can be used as a format for supervisory control and permit short-term local autonomous operations. Image processing algorithms can be made simpler and faster without trying to construct sensible images from the bottom. Control studies of telerobots lead to preferential manual control modes and, in this university environment, to basic paradigms for human motion and thence, perhaps, to redesign of robotic control, trajectory path planning, and rehabilitation prosthetics. Speculation as to future industrial drives for this telerobotic field suggests efficient roles for government agencies such as NASA     

The Cassini Visual And Infrared Mapping Spectrometer (Vims) Investigation

The Cassini visual and infrared mapping spectrometer (VIMS) investigation is a multidisciplinary study of the Saturnian system. Visual and near-infrared imaging spectroscopy and high-speed spectrophotometry are the observational techniques. The scope of the investigation includes the rings, the surfaces of the icy satellites and Titan, and the atmospheres of Saturn and Titan. In this paper, we will elucidate the major scientific and measurement goals of the investigation, the major characteristics of the Cassini VIMS instrument, the instrument calibration, and operation, and the results of the recent Cassini flybys of Venus and the Earth–Moon system.This revised version was published online in July 2005 with a corrected cover date.     

The Cassini Cosmic Dust Analyzer

The Cassini-Huygens Cosmic Dust Analyzer (CDA) is intended to provide direct observations of dust grains with masses between 10⁻¹⁹ and 10⁻⁹ kg in interplanetary space and in the jovian and saturnian systems, to investigate their physical, chemical and dynamical properties as functions of the distances to the Sun, to Jupiter and to Saturn and its satellites and rings, to study their interaction with the saturnian rings, satellites and magnetosphere. Chemical composition of interplanetary meteoroids will be compared with asteroidal and cometary dust, as well as with Saturn dust, ejecta from rings and satellites. Ring and satellites phenomena which might be effects of meteoroid impacts will be compared with the interplanetary dust environment. Electrical charges of particulate matter in the magnetosphere and its consequences will be studied, e.g. the effects of the ambient plasma and the magnetic field on the trajectories of dust particles as well as fragmentation of particles due to electrostatic disruption.The investigation will be performed with an instrument that measures the mass, composition, electric charge, speed, and flight direction of individual dust particles. It is a highly reliable and versatile instrument with a mass sensitivity 10⁶ times higher than that of the Pioneer 10 and 11 dust detectors which measured dust in the saturnian system. The Cosmic Dust Analyzer has significant inheritance from former space instrumentation developed for the VEGA, Giotto, Galileo, and Ulysses missions. It will reliably measure impacts from as low as 1 impact per month up to 10⁴ impacts per second. The instrument weighs 17 kg and consumes 12 W, the integrated time-of-flight mass spectrometer has a mass resolution of up to 50. The nominal data transmission rate is 524 bits/s and varies between 50 and 4192 bps.This revised version was published online in July 2005 with a corrected cover date.     

Fiber optic fence sensor developments

Many detection technologies have been employed for perimeter detection, to sensitize a barrier. These outdoor perimeter fence detection sensors must reliably detect intruders attempting to cut or climb the barrier, while ignoring the effects of environmental noise, including nearby activity. In recent conference proceedings, the new IntelliFIBER/spl trade/ fiber optic based product was introduced and compared with previous technologies. IntelliFIBER is designed to provide the advantages of a nonconductive "dielectric" cable sensor, e.g., resistance to electromagnetic interference and the ability to provide longer cable zones. It utilizes the proven processor of the Intelli-FLEX/spl trade/ triboelectric cable fence sensor, which provides relay or bidirectional networked serial communications and is compatible with its calibration module. IntelliFIBER also capitalizes on the Intelli-FLEX's adaptive detection algorithms that were developed through extensive field-testing, to provide exceptional immunity to environmental alarms while still detecting the skilled intruder.     

Limits in tracking with extended Kalman filters

The classical linearized conversion of measurements from polar or spherical coordinates to Cartesian ones generates a bias restricting the use of this conversion to cases where the bias can be neglected. In this work, the validity limits for the classical 2D transformation from polar to Cartesian coordinates, as derived in previous work, are shown to be too restrictive and the limits for the 3D transformation from spherical to Cartesian coordinates are introduced. Furthermore, quantitative measures for the performance degradation of the commonly used extended Kalman filter (EKF) in comparison with the best linear unbiased estimation (BLUE) filter are obtained by simulating typical tracking scenarios.     

Why space physics needs to go beyond the MHD box

Magnetohydrodynamic (MHD) theory has been used in space physics for more than forty years, yet many important questions about space plasmas remain unanswered. We still do not understand how the solar wind is accelerated, how mass, momentum and energy are transported into the magnetosphere and what mechanisms initiate substorms. Questions have been raised from the beginning of the space era whether MHD theory can describe correctly space plasmas that are collisionless and rarely in thermal equilibrium. Ideal MHD fluids do not induce electromotive force, hence they lose the capability to interact electromagnetically. No currents and magnetic fields are generated, rendering ideal MHD theory not very useful for space plasmas. Observations from the plasma sheet are used as examples to show how collisionless plasmas behave. Interpreting these observations using MHD and ideal MHD concepts can lead to misleading conclusions. Notably, the bursty bulk flows (BBF) with large mean velocities left(〈 v 〉≥400 km s right) that have been interpreted previously as E×B flows are shown to involve much more complicated physics. The sources of these nonvanishing 〈 v 〉 events, while still not known, are intimately related to mechanisms that create large phase space gradients that include beams and acceleration of ions to MeV energies. The distributions of these nonvanishing 〈 v 〉 events are associated with large amplitude variations of the magnetic field at frequencies up to and exceeding the local Larmor frequency where MHD theory is not valid. Understanding collisionless plasma dynamics such as substorms in the plasma sheet requires the self-consistency that only kinetic theory can provide. Kinetic modeling is still undergoing continual development with many studies limited to one and two dimensions, but there is urgent need to improve these models as more and more data show kinetic physics is fundamentally important. Only then will we be able to make progress and obtain a correct picture of how collisionless plasmas work in space.