Investigation of gas dynamics numerical schemes based on the solution of one-dimensional test problems

The application of the optimized dispersion relation preserving scheme (DRP-scheme) in combination with the explicit optimized two-layer Runge-Kutta scheme is presented to solve a system of one-dimensional and quasi-one-dimensional Euler equations using as an example the solution of four test problems, namely, discontinuity disintegration in a tube (Sod’s problem); transfer of the lowamplitude Gaussian pulse; acoustic wave propagation through the transonic nozzle; acoustic wave-shock interaction. Also given are the comparison of the calculation results using different schemes: DRP, CABARET, CE-SE and the standard Lax-Wendroff schemes as well as the solutions obtained with the use of software packages.     

Computer modelling of a penetrator thermal sensor

The Philae lander is part of the Rosetta mission to investigate comet 67P/Churyumov-Gerasimenko. It will use a harpoon like device to anchor itself onto the surface. The anchor will perhaps reach depths of 1–2 m. In the anchor is a temperature sensor that will measure the boundary temperature as part of the MUPUS experiment. As the anchor attains thermal equilibrium with the comet ice it may be possible to extract the thermal properties of the surrounding ice, such as the thermal diffusivity, by using the temperature sensor data. The anchor is not an optimal shape for a thermal probe and application of analytical solutions to the heat equation is inappropriate. We prepare a numerical model to fit temperature sensor data and extract the thermal diffusivity. Penetrator probes mechanically compact the material immediately surrounding them as they enter the target. If the thermal properties, composition and dimensions of the penetrator are known, then the thermal properties of this pristine material may be recovered although this will be a challenging measurement. We report on investigations, using a numerical thermal model, to simulate a variety of scenarios that the anchor may encounter and how they will affect the measurement.     

Wolf–Rayet optically thick winds with Alfvén waves

The Wolf–Rayet (WR) stars are hot luminous objects which are suffering an extreme mass loss via a continuous stellar wind. The high values of mass loss rates and high terminal velocities of the WR stellar winds constitute a challenge to the theories of radiation driven winds. Several authors incorporated magnetic forces to the line driven mechanism in order to explain these characteristics of the wind. Observations indicate that the WR stellar winds may reach, at the photosphere, velocities of the order of the terminal values, which means that an important part of the wind acceleration occurs at the optically thick region. The aim of this study is to analyze a model in which the wind in a WR star begins to be accelerated in the optically thick part of the wind. We used as initial conditions stellar parameters taken from the literature and solved the energy, mass and momentum equations. We demonstrate that the acceleration only by radiative forces is prevented by the general behavior of the opacities. Combining radiative forces plus a flux of Alfvén waves, we found in the simulations a fast drop in the wind density profile which strongly reduces the extension of the optically thick region and the wind becomes optically thin too close its base. The understanding how the WR wind initiate is still an open issue.     

Solar rhythms in the characteristics of the Arctic frontal zone in the North Atlantic

Long-term changes of the Arctic frontal zone characteristics near the south-eastern coasts of Greenland were considered, the NCEP/NCAR reanalysis data being used. It was found that in the cold half of the year the temperature gradients in the layer 1000–500 hPa in the region under study reveal strong ∼10-yr and ∼22-yr periodicities that seem to be related to solar activity cycles. The results obtained suggest the influence of solar activity and cosmic ray variations on the structure of the temperature field of the troposphere resulting in the changes of the temperature contrasts in the Arctic frontal zone that, in turn, may affect the intensity of cyclogenesis at middle latitudes. The detected effects seem to indicate an important part of frontal zones in the mechanism of solar activity and cosmic ray variation influence on the development of extratropical baric systems. It is suggested that the variations of the temperature gradients revealed in the Arctic frontal zone are due to the radiative forcing of cloudiness changes which may be associated with geomagnetic activity and cosmic ray variations.     

Studies of the geocenter motion using 16-years DORIS data

An accuracy of geocenter motion estimation is strongly dependent on the geodetic network size and stations distribution over the Earth’s surface. From this point of view DORIS system has an advantage, as its ground network of beacons consists of more than 50 sites, equally distributed over the Earth’s surface. Aiming to study variations of the geocenter movements, the results of DORIS data analysis for the time span 1993.0–2009.0 (inawd06.snx series), performed at the Analysis Centre of the Institute of astronomy of the Russian Academy of Sciences, have been used. DORIS data processing was made with GIPSY/OASIS II software, developed by Jet Propulsion Laboratory and modified for DORIS data processing by Institute Géographique National. Standard deviations of stations coordinates are estimated at the level 0.5–4.0 cm (internal consistency), depending on the number of satellites used in the solution. RMS of estimated components of the DORIS satellites orbits, compared with the solutions of other IDS analysis centres, do not exceed 1–2 cm. Weekly solutions for coordinates have been transformed from free network solutions (inawd06.snx series) to a well defined terrestrial reference frame ITRF2005 with the use of seven parameters of Helmert transformation, which were examined with a view to study variations of the geocenter movements (ina05wd01.geoc time series). In order to estimate linear trend, amplitudes, periods and phases of geocenter variation a method of linear regression was applied. The evaluated amplitudes of annual variations are of the order of 5–7 mm for X and Y components and 27–29 mm for Z component. Semi-annual amplitudes are also noticeable in all components (1–34 mm for X, Y and Z components). Secular trends in the DORIS geocenter coordinates are: −1.2, −0.1 and −0.3 mm/year for X, Y and Z directions respectively.     

High-speed streams of the solar wind near the earths’s orbit and their sources on the sun according to stereoscopic observations in the minimum of the 23rd cycle

Based on the results of plasma and magnetic measurements at three different points of the heliosphere and telescopic observations of the Sun from these points we study simultaneously high-speed streams (HSS) of the solar wind (SW) near the Earths’s orbit and coronal holes (CH) that have generated them. The data from spacecraft STEREO-A, STEREO-B, ACE, and SOHO are used together with ground-based observations from March 2007 to May 2008. In this period there existed HSS whose sources represented CH of various polarity, geometry, and location relative to the heliographic and heliomagnetic equators. Dependence of SW parameters on mutual positions of spacecraft with respect to CH and heliospheric current sheet, and also on heliolatitude and geometry of the CH is revealed. A difference of more than 5° in locations of spacecraft with respect to the heliospheric current sheet in November 2007 allowed us to discover a heliolatitude velocity gradient of the SW streams between the STEREO-A and STEREO-B spacecraft. On the average this gradient at that time was equal to 20 km/s per degree. Substantial variations in SW streams associated with variations of the HSS SW sources during a few hours or days were also observed. This variability makes it difficult to use the data of spacecraft STEREO-B for sufficiently accurate prediction of SW properties in the near-Earth space by the method of simple advanced ti me shift due to heliolongitude difference between a spacecraft and the Earth even in solar activity minimum.     

Electron flux variations at altitudes of 600–800 km in the second half of 2014. preliminary results of an experiment using RELEC equipment onboard the satellite <Emphasis Type="Italic">VERNOV</Emphasis>

The results of measurements of fluxes and spectra carried out using the RELEC (relativistic electrons) equipment onboard the VERNOV satellite in the second half of 2014 are presented. The VERNOV satellite was launched on July 8, 2014 in a sun-synchronous orbit with an altitude from 640 to 830 km and an inclination of 98.4°. Scientific information from the satellite was first received on July 20, 2014. The comparative analysis of electron fluxes using data from RELEC and using experimental data on the electron detection by satellites Elektro-L (positioned at a geostationary orbit) and Meteor-M no. 2 (positioned at a circular polar orbit at an altitude of about 800 km as the VERNOV satellite) will make it possible to study the spatial distribution pattern of energetic electrons in near-Earth space in more detail.     

Variations in cosmic ray cutoff rigidities during the great geomagnetic storm of November 2004

A very strong interplanetary and magnetospheric disturbance observed on 7–13 November 2004 can be regarded as one of the strongest events during the entire period of space observations. In this paper we report on the studies of cosmic ray cutoff rigidity variations during 7–13 November 2004 showing how storm conditions can affect the direct cosmic ray access to the inner magnetosphere. Effective cutoff rigidities have been calculated for selected points on the ground by tracing trajectories of cosmic ray particles through the magnetospheric magnetic field of the “storm-oriented” Tsyganenko 2003 model. Cutoff rigidity variations have also been determined by the spectrographic global survey method on the basis of experimental data of the neutron monitor network. Relations between the calculated and experimental cutoff rigidities and the geomagnetic Dst-index and interplanetary parameters have been investigated. Correlation coefficients between the cutoff rigidities obtained by the trajectory tracing method and the spectrographic global survey method have been found to be in the limits 0.76–0.89 for all stations except the low-latitude station Tokyo (0.35). The most pronounced correlation has been revealed between the cutoff rigidities that exhibited a very large variation of ∼1–1.5 GV during the magnetic storm and the Dst index.     

Height of shock formation in the solar corona inferred from observations of type II radio bursts and coronal mass ejections

Employing coronagraphic and EUV observations close to the solar surface made by the Solar Terrestrial Relations Observatory (STEREO) mission, we determined the heliocentric distance of coronal mass ejections (CMEs) at the starting time of associated metric type II bursts. We used the wave diameter and leading edge methods and measured the CME heights for a set of 32 metric type II bursts from solar cycle 24. We minimized the projection effects by making the measurements from a view that is roughly orthogonal to the direction of the ejection. We also chose image frames close to the onset times of the type II bursts, so no extrapolation was necessary. We found that the CMEs were located in the heliocentric distance range from 1.20 to 1.93 solar radii (Rs), with mean and median values of 1.43 and 1.38 Rs, respectively. We conclusively find that the shock formation can occur at heights substantially below 1.5 Rs. In a few cases, the CME height at type II onset was close to 2 Rs. In these cases, the starting frequency of the type II bursts was very low, in the range 25–40 MHz, which confirms that the shock can also form at larger heights. The starting frequencies of metric type II bursts have a weak correlation with the measured CME/shock heights and are consistent with the rapid decline of density with height in the inner corona.