On the maximum time step in weakly compressible SPH

International audienceIn the SPH method for viscous fluids, the time step is subject to empirical stability criteria. We proceed to a stability analysis of the Weakly Compressible SPH equations using the von Neumann approach in arbitrary space dimension for unbounded flow. Considering the continuous SPH interpolant based on integrals, we obtain a theoretical stability criterion for the time step, depending on the kernel standard deviation, the speed of sound and the viscosity. The stability domain appears to be almost independent of the kernel choice for a given space discretisation. Numerical tests show that the theory is very accurate, despite the approximations made. We then extend the theory in order to study the influence of the method used to compute the density, of the gradient and divergence SPH operators, of background pressure, of the model used for viscous forces and of a constant velocity gradient. The influence of time integration scheme is also studied, and proved to be prominent. All of the above theoretical developments give excellent agreement against numerical results. It is found that velocity gradients almost do not affect stability, provided some background pressure is used. Finally, the case of bounded flows is briefly addressed from numerical tests in three cases: a laminar Poiseuille flow in a pipe, a lid-driven cavity and the collapse of a water column on a wedge

Onzain – Les Basses Terres Noires

Lien Atlas (MCC) : http://atlas.patrimoines.culture.fr/atlas/trunk/index.php?ap_theme=DOM_2.01.02&ap_bbox=1.115;47.479;1.215;47.554 Le dépôt des Basses Terres Noires, sur la commune d’Onzain, a été mis au jour clandestinement en novembre 1997 par un utilisateur de détecteur de métaux qui avait extrait du sol un important lot d’objets en bronze. Une fouille en urgence absolue a été menée en janvier 1998 sur le lieu de la découverte, situé sur la basse terrasse de la Loire, en rive droite du fl..

Onzain – Les Basses Terres Noires

Lien Atlas (MCC) : http://atlas.patrimoines.culture.fr/atlas/trunk/index.php?ap_theme=DOM_2.01.02&ap_bbox=1.115;47.479;1.215;47.554 Le dépôt des Basses Terres Noires, sur la commune d’Onzain, a été mis au jour clandestinement en novembre 1997 par un utilisateur de détecteur de métaux qui avait extrait du sol un important lot d’objets en bronze. Une fouille en urgence absolue a été menée en janvier 1998 sur le lieu de la découverte, situé sur la basse terrasse de la Loire, en rive droite du fl..

Résultats significatifs en Martinique pour l’année 2016

19 autorisations d’opérations archéologiques ont été délivrées au cours de l’année 2016. 17 d’entre-elles se rapportent à des diagnostics d’archéologie préventive, tandis que les deux autres concernent des recherches programmées. Le présent bilan ne pourra cependant rendre compte que des seules opérations préventives. Les recherches programmées, en raison d’impondérables, n’ont pu être mises en œuvre sur le terrain au cours de l’année et se sont trouvées reportées. Il s’agit tout d’abord d’un..

Résultats significatifs en Martinique pour les années 2008, 2009 et 2010

Des prospections diachroniques ont été mises en œuvre pour l’enrichissement de la carte archéologique tantôt dans des secteurs où des projets d’aménagement, de gestion ou de mise en valeur invitent à une prise en compte du patrimoine archéologique, tantôt pour mesurer les effets des phénomènes cycloniques et des houles sur la conservation de certains sites du littoral. Archéologie précolombienne La recherche en archéologie précolombienne s’est nourrie des résultats de plusieurs opérations. À ..

Résultats significatifs en Martinique pour les années 2011, 2012 et 2013

Un programme de prospection des sites littoraux amérindiens de Martinique a permis une actualisation de la carte archéologique, de détecter de nouveaux indices d’occupations précolombiennes et d’identifier les portions du littoral où l’érosion marine se fait la plus menaçante pour la conservation de ces sites. Une prospection inventaire entre Anse Couleuvre, au Prêcheur, et le bourg de Grand’Rivière a également permis le repérage de quelques structures d’époque coloniale. Archéologie précolom..

Résultats significatifs en Martinique pour l’année 2017

24 opérations archéologiques de terrain ont été réalisées sur le territoire martiniquais au cours de l’année 2017. On ne sera guère étonné de constater que l’essentiel de l’activité relève de l’archéologie préventive, pour laquelle on dénombre 13 diagnostics et 6 fouilles préventives. Ce dernier chiffre doit d’ailleurs être remarqué en comparaison de l’année 2016 où aucune opération de ce type n’avait été mise en œuvre par les aménageurs. La recherche programmée connaît également une progress..

Unified semi-analytical wall boundary conditions applied to 2-D incompressible SPH

International audienceThis work aims at improving the 2-D incompressible SPH model (ISPH) by adapting it to the unified semi-analytical wall boundary conditions proposed by Ferrand et al. [10]. The ISPH algorithm considered is as proposed by Lind et al. [25], based on the projection method with a divergence-free velocity field and using a stabilising procedure based on particle shifting. However, we consider an extension of this model to Reynolds-Averaged Navier-Stokes equations based on the k- turbulent closure model, as done in [10]. The discrete SPH operators are modified by the new description of the wall boundary conditions. In particular, a boundary term appears in the Laplacian operator, which makes it possible to accurately impose a von Neumann pressure wall boundary condition that corresponds to impermeability. The shifting and free-surface detection algorithms have also been adapted to the new boundary conditions. Moreover, a new way to compute the wall renormalisation factor in the frame of the unified semi-analytical boundary conditions is proposed in order to decrease the computational time. We present several verifications to the present approach, including a lid-driven cavity, a water column collapsing on a wedge and a periodic schematic fish-pass. Our results are compared to Finite Volumes methods, using Volume of Fluids in the case of free-surface flows. We briefly investigate the convergence of the method and prove its ability to model complex free-surface and turbulent flows. The results are generally improved when compared to a weakly compressible SPH model with the same boundary conditions, especially in terms of pressure prediction

Modeling the evolution of spectral induced polarization during calcite precipitation on glass beads

International audienceWhen pH and alkalinity increase, calcite frequently precipitates and hence modifies the petrophysical properties of porous media. The complex conductivity method can be used to directly monitor calcite precipitation in porous media because it is very sensitive to the evolution of the pore structure and its connectivity. We have developed a mechanistic grain polarization model considering the electrochemical polarization of the Stern layer surrounding calcite particles. This model depends on the surface charge density and mobility of the counter-ions in the Stern layer. Our induced polarization model predicts the evolution of the size of calcite particles, of the pore structure and connectivity during spectral induced polarization experiments of calcite precipitation on glass beads pack. Model predictions are in very good agreement with the complex conductivity measurements. During the first phase of calcite precipitation experiment, calcite crystals growth, and the inverted particle size distribution moves towards larger calcite particles. When calcite continues to precipitate and during pore clogging, inverted particle size distribution moves towards smaller particles because large particles do not polarize sufficiently. The pore clogging is also responsible for the decrease of the connectivity of the pores, which is observed through the increasing electrical formation factor of the porous medium

A concurrent, generational garbage collector for a multithreaded implementation of ML

International audienceThis paper presents the design and implementation of a "quasi real-time" garbage collector for Concurrent Caml Light, an implementation of ML with threads. This two-generation system combines a fast, asynchronous copying collector on the young generation with a non-disruptive concurrent marking collector on the old generation. This design crucially relies on the ML compile-time distinction between mutable and immutable objects