The effect of Poisson noise on SPH calculations

A simple algorithm is presented which generates a lattice-like, regular spacing of smoothed particle hydrodynamics (SPH) particles in discs, with any prescribed density gradient. Using this for comparison we demonstrate the effect of Poisson noise on SPH estimates of density, pressure and viscous forces when calculated using randomly distributed particles. The standard deviation of density and pressure is typically increased by greater than an order of magnitude. In a disc with a Keplerian velocity profile, the effectiveness of the Balsara switch in reducing the shear component of SPH artificial viscosity is greatly enhanced when the particles are properly spaced, reducing the magnitude of viscosity by two orders of magnitude. Noise problems are exacerbated, not removed, by increasing the numbers of SPH particles, if the number of neighbours used is kept constant. However, comparison of the evolution of a disc created using perfectly spaced particles and a disc with identical parameters but randomly placed particles, reveals very similar results. Although there are subtle differences in the evolution, and the smooth disc takes longer to begin developing structure, in both cases the identical number of objects is created by gravitational collapse. SPH disc simulations do not depend on initial density disturbances to evolve objects by gravitational collapse, which gives added credence to the validity of the results. It also appears that complicated disc settling procedures are unnecessary.     

Non-singular time domain BEM with applications to 3D inertial soil–structure interaction

A direct time domain boundary element method is presented based on the Stokes fundamental solutions, discretized in both time and space, and an efficient time step-by-step solution that minimizes the accumulation of errors. A non-singular numerical integration procedure, in the Cauchy sense, is proposed for the generation of the associated influence matrices. This methodology is shown to be efficient for the solution of a number of computationally intensive problems in the area of soil–structure interaction. In addition, an algorithm for the direct calculation of the response of massive foundations to externally applied forces and/or obliquely incident seismic waves is introduced. The accuracy and computational efficiency of the proposed methodologies is established through a number of comparison studies.     

Dynamic 3-D soil–railway track interaction by BEM–FEM

A study on the dynamic response of a railway track is presented via a 3-D formulation based on the frequency domain Boundary Element Method (BEM) and the Finite Element Method (FEM). The railway track consists of a group of surface, massive, rigid footings resting on a viscoelastic half-space and connected by an overlaying rail structure. The BEM, employing the full-space fundamental solutions and quadrilateral elements, is used for the simulation of the elastic half-space while the FEM is used to model the rigid footings and the rail superstructure. The loading function consists of a set of externally applied, harmonic or transient loads. Frequency as well as transient, by way of FFT, results are presented for various modes of vibration. Various numerical studies assess the through-the-soil interaction of the adjacent footings, the influence of soil damping, the effect of the overlaying structure on the frequency content of the system, and the effective simulation of an infinitely long railway track by a truncated one.     

Ionospheric specification and forecasting based on observations from European ionosondes participating in DIAS project

There are two main objectives of the DIAS (European Digital Upper Atmosphere Server) project. First, it establishes a pan-European repository of raw and derived digital data describing the state of ionospheric part of the upper atmosphere, which is capable of ingesting real-time information and maintaining historical data collections provided by most operating ionospheric stations in Europe. Second, the DIAS system produces and distributes, based on the raw data collection, several operational products required by various user groups for ionospheric nowcasting, prediction and forecasting purposes. The project completed on May 2006 and the DIAS server operates since then continuously. The basic products that are delivered are real-time and historical ionograms from all DIAS-affiliated ionospheric stations, frequency plots and maps of the ionosphere over Europe based on the foF2, M(3000)F2, MUF and electron density parameters, as well as long term prediction and short term forecasting up to 24 hour ahead. The paper describes use of the ionospheric measurements in the DIAS modelling techniques for specification, predict-tion and forecasting of the ionosphere over the European region, and details the final products available to the DIAS user community.     

The diffuse-to-global and diffuse-to-direct-beam spectral irradiance ratios as turbidity indexes in an urban environment

Continuous measurements of solar spectral radiation using the Multi-filter Rotating Shadow Band Radiometer (MFRSR) are performed at the Actinometric Station of the National Observatory of Athens (ASNOA). The present study utilizes three clear-sky days of continuous observations, from local sunrise to local sunset, in order to investigate the daily variation of the radiation components (diffuse, global and direct-beam) as well as their ratios (diffuse-to-global, DGR, and diffuse-to-direct-beam, DDR) under different atmospheric conditions. Both ratios have received great scientific interest, especially for investigating solar irradiance modifications under various atmospheric conditions, aerosol load and optical properties. Apart from this, the present study shows that the DDR can constitute a measure of atmospheric turbidity when it is determined at longer wavelengths, while the DGR cannot. The effect of the solar zenith angle (SZA) on both ratios is significant at the shorter wavelengths with varying sensitivities depending on the aerosol field and sun elevation. The present study confirms the results obtained by previous solar irradiance measurements in Athens and also those computed via radiative transfer codes and sheds light on the scientific knowledge of the use of spectral DDR as an atmospheric turbidity index.     

Numerical simulation of dynamic soil–structure interaction in shaking table testing

This paper provides an insight into the numerical simulation of soil–structure interaction (SSI) phenomena studied in a shaking table facility. The shaking table test is purposely designed to confirm the ability of the numerical substructure technique to simulate the SSI phenomenon. A model foundation–structure system with strong SSI potential is embedded in a dry bed of sand deposited within a purpose designed shaking-table soil container. The experimental system is subjected to a strong ground motion. The numerical simulation of the complete soil–foundation–structure system is conducted in the linear viscoelastic domain using the substructure approach. The matching of the experimental and numerical responses in both frequency and in time domain is satisfying. Many important aspects of SSI that are apparent in the experiment are captured by the numerical simulation. Furthermore, the numerical modelling is shown to be adequate for practical engineering design purposes.     

Meta-analysis of electron cyclotron resonance absorption features detected in high-mass X-ray binaries

Using recent compilations of detailed X-ray observations and spectral models of exceptional quality, we record the electron cyclotron resonance absorption(ECRA) features that have been detected in 45 pulsating high-mass X-ray binaries(HMXBs) and ultraluminous X-ray(ULX) sources harboring neutron stars, although seven of these detections are still questionable and another 21 are single and/or not independently confirmed. From the comprehensive catalogs of Jaisawal Naik and Staubert et al.and from several additional recent observations, we produce two lists of HMXB ECRA sources: a list of 17 sources in which multiple ECRA lines or single very low-energy lines are seen, in which we can reasonably assume that the lowest energy reveals the fundamental cyclotron level for each source; and a"contaminated" list of 38 sources including the 21 detections of single ECRA lines that may(not) be higherlevel harmonics. Both lists confirm a previous result that we have obtained independently by modeling the propeller lines of Magellanic HMXB pulsars: the surface dipolar magnetic fields B*of HMXB neutron stars are segregated around five distinct values with B*= 0.28 ± 0.08, 0.55 ± 0.11, 1.3 ± 0.37, 3.0 ± 0.68 and 7.9 ± 3.1, in units of TG. However, an explanation of this phenomenon is currently lacking. We have found no correlation between these B*values and the corresponding observed spin periods, spin period derivatives, orbital periods, maximum X-ray luminosities, neutron star masses or companion star masses.     

A striking confluence between theory and observations of high-mass X-ray binary pulsars

We analyze the most powerful X-ray outbursts from neutron stars in eleven Magellanic high-mass X-ray binaries and three pulsating ultraluminous X-ray sources. Most of the outbursts rise to L_(max) which is about the level of the Eddington luminosity, while the remaining more powerful outbursts also appear to recognize that limit when their emissions are assumed to be anisotropic and beamed toward our direction. We use the measurements of pulsar spin periods P_S and their derivatives P_S to calculate the X-ray luminosities L_p in their faintest accreting("propeller-line") states. In five cases with unknown P_S, we use the lowest observed X-ray luminosities, which only adds to the heterogeneity of the sample. Then we calculate the ratios L_p/L_(max) and we obtain an outstanding confluence of theory and observations from which we conclude that work done on both fronts is accurate and the results are trustworthy: sources known to reside on the lowest Magellanic propeller line are all located on/near that line, whereas other sources jump higher and reach higher-lying propeller lines. These jumps can be interpreted in only one way, higher-lying pulsars have stronger surface magnetic fields in agreement with previous empirical results in whichP_S and L_p values were not used.     

Burned area mapping and post-fire impact assessment in the Kassandra peninsula (Greece) using Landsat TM and Quickbird data

The purpose of this study was to assess the environmental impacts of forest fires on part of the Mediterranean basin. The study area is on the Kassandra peninsula, prefecture of Halkidiki, Greece. A maximum likelihood supervised classification was applied to a post-fire Landsat TM image for mapping the exact burned area. Land-cover types that had been affected by fire were identified with the aid of a CORINE land-cover type layer. Results showed an overall classification accuracy of 95%, and 83% of the total burned area was ‘forest areas’. A normalized difference vegetation index threshold technique was applied to a post-fire Quickbird image which had been recorded six years after the fire event to assess the vegetation recovery and to identify the vegetation species that were dominant in burned areas. Four classes were identified: ‘bare soil’, ‘sparse shrubs’, ‘dense shrubs’ and ‘tree and shrub communities’. Results showed that ‘shrublands’ is the main vegetation type which has prevailed (65%) and that vegetation recovery is homogeneous in burned areas.     

Recursive anisotropy: a spatial taphonomic study of the Early Pleistocene vertebrate assemblage of Tsiotra Vryssi,Mygdonia Basin,Greece

By applying advanced spatial statistical methods, spatial taphonomy complements the traditional taphonomic approach and enhances our understanding of biostratinomic and diagenetic processes. In this study, we elaborate on a specific aspect – spatial anisotropy – of taphonomic processes. We aim to unravel the taphonomic history of the Early Pleistocene vertebrate assemblage of Tsiotra Vryssi (Mygdonia Basin, Greece). Circular statistics are used for the fabric analysis of elongated elements; geostatistics (directional variograms), wavelet and point pattern analyses are applied for detecting anisotropy at the assemblage level. The anisotropy of magnetic susceptibility (AMS) of sedimentary magnetic minerals is also investigated. The results, integrated with preliminary remarks about the differential preservation of skeletal elements, sedimentological and micromorphological observations, suggest multiple dispersion events and recurrent spatial re‐arrangement of a lag, (peri)autochthonous assemblage, consistent with the cyclical lateral switching of a braided fluvial system. Furthermore, this study offers a contribution to the building of a spatial taphonomic referential framework for the interpretation of other fossil vertebrate assemblages, including archaeo‐palaeontological ones.