ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation

Resolving numerically Vlasov-Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm consisting in representing the phase-space sheet with a conforming, self-adaptive simplicial tessellation of which the vertices follow the Lagrangian equations of motion. The algorithm is implemented both in six- and four-dimensional phase-space. Refinement of the tessellation mesh is performed using the bisection method and a local representation of the phase-space sheet at second order relying on additional tracers created when needed at runtime. In order to preserve in the best way the Hamiltonian nature of the system, refinement is anisotropic and constrained by measurements of local Poincar\'e invariants. Resolution of Poisson equation is performed using the fast Fourier method on a regular rectangular grid, similarly to particle in cells codes. To compute the density projected onto this grid, the intersection of the tessellation and the grid is calculated using the method of Franklin and Kankanhalli (1993) generalised to linear order. As preliminary tests of the code, we study in four dimensional phase-space the evolution of an initially small patch in a chaotic potential and the cosmological collapse of a fluctuation composed of two sinusoidal waves. We also perform a "warm" dark matter simulation in six-dimensional phase-space that we use to check the parallel scaling of the code.Comment: Code and illustration movies available at: http://www.vlasix.org/index.php?n=Main.ColDICE - Article submitted to Journal of Computational Physic

Dynamical criticality in the collective activity of a population of retinal neurons

Recent experimental results based on multi-electrode and imaging techniques have reinvigorated the idea that large neural networks operate near a critical point, between order and disorder. However, evidence for criticality has relied on the definition of arbitrary order parameters, or on models that do not address the dynamical nature of network activity. Here we introduce a novel approach to assess criticality that overcomes these limitations, while encompassing and generalizing previous criteria. We find a simple model to describe the global activity of large populations of ganglion cells in the rat retina, and show that their statistics are poised near a critical point. Taking into account the temporal dynamics of the activity greatly enhances the evidence for criticality, revealing it where previous methods would not. The approach is general and could be used in other biological networks

The benefits of cooperation under uncertainty: the case of climate change

This article presents an analysis of the behavior of countries defining their climate policies in an uncertain context. The analysis is made using the S-CWS model, a stochastic version of an integrated assessment growth model. The model includes a stochastic definition of the climate sensitivity parameter. We show that the impact of uncertainty on policy design critically depends on the shape of the damage function. We also examine the benefits of cooperation in the context of uncertainty: we highlight the existence of an additional benefit of cooperation, namely risk reduction.cooperation, uncertainty, climate change, integrated assessment model

Phase equilibrium of the CO2/glycerol system: Experimental data by in situ FT-IR spectroscopy and thermodynamic modeling

Phase equilibrium experimental data for the CO2/glycerol system are reported in this paper. The measurements were performed using an in situ FT-IR method for temperatures ranging from 40 ◦C to 200 ◦C and pressures up to 35.0 MPa, allowing determination of the mutual solubility of both compounds. Concerning the CO2 rich phase, it was observed that the glycerol solubility in CO2 was extremely low (in the range of 10−5 in mole fraction) in the pressure and temperature domains investigated here. Conversely, the glycerol rich phase dissolved CO2 at mole fractions up to 0.13. Negligible swelling of the glycerol rich phase has been observed. Modeling of the phase equilibrium has been performed using the Peng–Robinson equation of state (PR EoS) with classical van der Waals one fluid and EoS/GE based mixing rules (PSRK and MHV2). Satisfactory agreement was observed between modeling results and experimental measurements when PSRK mixing rules are used in combination with UNIQUAC model, although UNIFAC predictive approach gives unsatisfactory representation of experimental behavior

Logical segmentation for article extraction in digitized old newspapers

Newspapers are documents made of news item and informative articles. They are not meant to be red iteratively: the reader can pick his items in any order he fancies. Ignoring this structural property, most digitized newspaper archives only offer access by issue or at best by page to their content. We have built a digitization workflow that automatically extracts newspaper articles from images, which allows indexing and retrieval of information at the article level. Our back-end system extracts the logical structure of the page to produce the informative units: the articles. Each image is labelled at the pixel level, through a machine learning based method, then the page logical structure is constructed up from there by the detection of structuring entities such as horizontal and vertical separators, titles and text lines. This logical structure is stored in a METS wrapper associated to the ALTO file produced by the system including the OCRed text. Our front-end system provides a web high definition visualisation of images, textual indexing and retrieval facilities, searching and reading at the article level. Articles transcriptions can be collaboratively corrected, which as a consequence allows for better indexing. We are currently testing our system on the archives of the Journal de Rouen, one of France eldest local newspaper. These 250 years of publication amount to 300 000 pages of very variable image quality and layout complexity. Test year 1808 can be consulted at plair.univ-rouen.fr.Comment: ACM Document Engineering, France (2012

Numerical simulation of a stroke: numerical problems and methodology

The numerical simulation of an ischemic stroke is a challenging problem: a complicated geometry, and a very stiff large system of reaction-diffusion equations. This paper, intended for mathematicians as well as for biologist, gives a survey and an introduction to the numerical methods used and some results