Astrocomp: a web service for the use of high performance computers in Astrophysics

Astrocomp is a joint project, developed by the INAF-Astrophysical Observatory of Catania, University of Roma La Sapienza and Enea. The project has the goal of providing the scientific community of a web-based user-friendly interface which allows running parallel codes on a set of high-performance computing (HPC) resources, without any need for specific knowledge about parallel programming and Operating Systems commands. Astrocomp provides, also, computing time on a set of parallel computing systems, available to the authorized user. At present, the portal makes a few codes available, among which: FLY, a cosmological code for studying three-dimensional collisionless self-gravitating systems with periodic boundary conditions; ATD, a parallel tree-code for the simulation of the dynamics of boundary-free collisional and collisionless self-gravitating systems and MARA, a code for stellar light curves analysis. Other codes are going to be added to the portal.Comment: LaTeX with elsart.cls and harvard.sty (included). 7 pages. To be submitted to a specific journa

Development of a decision support system framework for cultural heritage management

Decision support systems (DSSs) have been traditionally identified as useful information technology tools in a variety of fields, including the context of cultural heritage. However, to the best of our knowledge, no prior study has developed a DSS framework that incorporates all the main decision areas simultaneously in the context of cultural heritage. We fill this gap by focusing on design-science research and specifically by developing a DSS framework whose features support all the main decision areas for the sustainable management of cultural assets in a comprehensive manner. The main decision-making areas considered in our study encompass demand manage-ment, segmentation and communication, pricing, space management, and services management. For these areas, we select appropriate decision-making supporting techniques and data management solutions. The development of our framework, in the form of a web-based system, results in an architectural solution that is able to satisfy critical requirements such as ease of use and response time. We present an application of the innovative DSS framework to a museum and discuss the main managerial implications and future improvements

Timescales and radial lengthscales of quasi-periodic density structures observed by the Helios probes

Quasi-periodic density structures are important mesoscale structures that constitute the slow solar wind. The associated periodicities play a fundamental role in the understanding of the release processes of solar coronal plasma and provide important constraints on models trying to understand the origin of the slow solar wind. However, observational restrictions have limited their study to coronagraph images near the Sun and in situ measurements at 1 AU, preventing a clear connection between in situ structures to their coronal sources. A first step toward a better understanding of this problem is the analysis of Helios measurements to probe the solar wind between 0.3 and 0.6 AU. Three instances of quasi-periodic density structures have been identified showing characteristic timescales of ≈31-33 minutes and ≈112-120 minutes. Using the mean radial bulk speed corresponding to the time intervals of the structures, the related radial lengthscales are derived showing a characteristic value of ≈120-150 RE

A computer simulation protocol to assess the accuracy of a Radio Stereometric Analysis (RSA) image processor according to the ISO-5725

Radio-Stereometric-Analysis and x-ray fluoroscopy are radiological techniques that require dedicated software to process data. The accurate calibration of these software is therefore critical. The aim of this work is to produce a protocol for evaluating the softwares' accuracy according to the ISO-5725. A series of computer simulations of the radiological setup and images were employed. The noise level of the images was also changed to evaluate the accuracy with different image qualities. The protocol was tested on a custom software developed by the authors. Radiological scene reconstruction accuracy was of (0.092 +- 0.14) mm for tube position, and (0.38 +- 0.31) mm / (2.09 +- 1.39) deg for detectors oriented in a direction other than the source-detector direction. In the source-detector direction the accuracy was of (2.68 +- 3.08) mm for tube position, and of (0.16 +- 0.27) mm / (0.075 +- 1.16) deg for the detectors. These disparate results are widely discussed in the literature. Model positioning and orientation was also highly accurate: (0.22 +- 0.46) mm / (0.26 +- 0.22) deg. Accuracy was not affected by the noise level. The protocol was able to assess the accuracy of the RSA system. It was also useful to detect and fix hidden bugs. It was also useful to detect and resolve hidden bugs in the software, and in optimizing the algorithms

Population Balance Models for Particulate Flows in Porous Media: Breakage and Shear-Induced Events

Transport and particulate processes are ubiquitous in environmental, industrial and biological applications, often involving complex geometries and porous media. In this work we present a general population balance model for particle transport at the pore-scale, including aggregation, breakage and surface deposition. The various terms in the equations are analysed with a dimensional analysis, including a novel collision-induced breakage mechanism, and split into one- and two-particles processes. While the first are linear processes, they might both depend on local flow properties (e.g. shear). This means that the upscaling (via volume averaging and homogenisation) to a macroscopic (Darcy-scale) description requires closures assumptions. We discuss this problem and derive an effective macroscopic term for the shear-induced events, such as breakage caused by shear forces on the transported particles. We focus on breakage events as prototype for linear shear-induced events and derive upscaled breakage frequencies in periodic geometries, starting from nonlinear power-law dependence on the local fluid shear rate. Results are presented for a two-dimensional channel flow and a three dimensional regular arrangement of spheres, for arbitrarily fast (mixing-limited) events. Implications for linearised shear-induced collisions are also discussed. This work lays the foundations of a new general framework for multiscale modelling of particulate flows

Hierarchical information clustering by means of topologically embedded graphs

We introduce a graph-theoretic approach to extract clusters and hierarchies in complex data-sets in an unsupervised and deterministic manner, without the use of any prior information. This is achieved by building topologically embedded networks containing the subset of most significant links and analyzing the network structure. For a planar embedding, this method provides both the intra-cluster hierarchy, which describes the way clusters are composed, and the inter-cluster hierarchy which describes how clusters gather together. We discuss performance, robustness and reliability of this method by first investigating several artificial data-sets, finding that it can outperform significantly other established approaches. Then we show that our method can successfully differentiate meaningful clusters and hierarchies in a variety of real data-sets. In particular, we find that the application to gene expression patterns of lymphoma samples uncovers biologically significant groups of genes which play key-roles in diagnosis, prognosis and treatment of some of the most relevant human lymphoid malignancies.Comment: 33 Pages, 18 Figures, 5 Table

Stochastic Cellular Automata Model for Stock Market Dynamics

In the present work we introduce a stochastic cellular automata model in order to simulate the dynamics of the stock market. A direct percolation method is used to create a hierarchy of clusters of active traders on a two dimensional grid. Active traders are characterised by the decision to buy, (+1), or sell, (-1), a stock at a certain discrete time step. The remaining cells are inactive,(0). The trading dynamics is then determined by the stochastic interaction between traders belonging to the same cluster. Most of the stylized aspects of the financial market time series are reproduced by the model.Comment: 17 pages and 7 figure

Ensembles of climate change models for risk assessment of nuclear power plants

Climate change affects technical Systems, Structures and Infrastructures (SSIs), changing the environmental context for which SSI were originally designed. In order to prevent any risk growth beyond acceptable levels, the climate change effects must be accounted for into risk assessment models. Climate models can provide future climate data, such as air temperature and pressure. However, the reliability of climate models is a major concern due to the uncertainty in the temperature and pressure future projections. In this work, we consider five climate change models (individually unable to accurately provide historical recorded temperatures and, thus, also future projections), and ensemble their projections for integration in a probabilistic safety assessment, conditional on climate projections. As case study, we consider the Passive Containment Cooling System (PCCS) of two AP1000 Nuclear Power Plants (NPPs). Results provided by the different ensembles are compared. Finally, a risk-based classification approach is performed to identify critical future temperatures, which may lead to PCCS risks beyond acceptable levels

Asymptotic Energy Dependence of Hadronic Total Cross Sections from Lattice QCD

The nonperturbative approach to soft high-energy hadron-hadron scattering, based on the analytic continuation of Wilson-loop correlation functions from Euclidean to Minkowskian theory, allows to investigate the asymptotic energy dependence of hadron-hadron total cross sections in lattice QCD. In this paper we will show, using best fits of the lattice data with proper functional forms satisfying unitarity and other physical constraints, how indications emerge in favor of a universal asymptotic high-energy behavior of the kind B log^2 s for hadronic total cross sections.Comment: Revised and extended version; 29 pages, 4 figure