Extent of force indeterminacy in packings of frictional rigid disks

Static packings of frictional rigid particles are investigated by means of discrete element simulations. We explore the ensemble of allowed force realizations in the space of contact forces for a given packing structure. We estimate the extent of force indeterminacy with different methods. The indeterminacy exhibits a nonmonotonic dependence on the interparticle friction coefficient. We verify directly that larger force-indeterminacy is accompanied by a more robust behavior against local perturbations. We also investigate the local indeterminacy of individual contact forces. The probability distribution of local indeterminacy changes its shape depending on friction. We find that local indeterminacy tends to be larger on force chains for intermediate friction. This correlation disappears in the large friction limit.Comment: 5 pages, 6 figure

Stress-strain behavior and geometrical properties of packings of elongated particles

We present a numerical analysis of the effect of particle elongation on the quasistatic behavior of sheared granular media by means of the Contact Dynamics method. The particle shapes are rounded-cap rectangles characterized by their elongation. The macroscopic and microstructural properties of several packings subjected to biaxial compression are analyzed as a function of particle elongation. We find that the shear strength is an increasing linear function of elongation. Performing an additive decomposition of the stress tensor based on a harmonic approximation of the angular dependence of branch vectors, contact normals and forces, we show that the increasing mobilization of friction force and the associated anisotropy are key effects of particle elongation. These effects are correlated with partial nematic ordering of the particles which tend to be oriented perpendicular to the major principal stress direction and form side-to-side contacts. However, the force transmission is found to be mainly guided by cap-to-side contacts, which represent the largest fraction of contacts for the most elongated particles. Another interesting finding is that, in contrast to shear strength, the solid fraction first increases with particle elongation, but declines as the particles become more elongated. It is also remarkable that the coordination number does not follow this trend so that the packings of more elongated particles are looser but more strongly connected.Comment: Submited to Physical Review

The scenario of two-dimensional instabilities of the cylinder wake under EHD forcing: A linear stability analysis

We propose to study the stability properties of an air flow wake forced by a dielectric barrier discharge (DBD) actuator, which is a type of electrohydrodynamic (EHD) actuator. These actuators add momentum to the flow around a cylinder in regions close to the wall and, in our case, are symmetrically disposed near the boundary layer separation point. Since the forcing frequencies, typical of DBD, are much higher than the natural shedding frequency of the flow, we will be considering the forcing actuation as stationary. In the first part, the flow around a circular cylinder modified by EHD actuators will be experimentally studied by means of particle image velocimetry (PIV). In the second part, the EHD actuators have been numerically implemented as a boundary condition on the cylinder surface. Using this boundary condition, the computationally obtained base flow is then compared with the experimental one in order to relate the control parameters from both methodologies. After validating the obtained agreement, we study the Hopf bifurcation that appears once the flow starts the vortex shedding through experimental and computational approaches. For the base flow derived from experimentally obtained snapshots, we monitor the evolution of the velocity amplitude oscillations. As to the computationally obtained base flow, its stability is analyzed by solving a global eigenvalue problem obtained from the linearized Navier–Stokes equations. Finally, the critical parameters obtained from both approaches are compared

The post-Higgs MSSM scenario: Habemus MSSM?

We analyze the Minimal Supersymmetric extension of the Standard Model that we have after the discovery of the Higgs boson at the LHC, the hMSSM (habemus MSSM?), i.e. a model in which the lighter $h$ boson has a mass of approximately 125 GeV which, together with the non-observation of superparticles at the LHC, indicates that the SUSY-breaking scale $M_S$ is rather high, $M_S > 1$ TeV. We first demonstrate that the value $M_h \approx 125$ GeV fixes the dominant radiative corrections that enter the MSSM Higgs boson masses, leading to a Higgs sector that can be described, to a good approximation, by only two free parameters. In a second step, we consider the direct supersymmetric radiative corrections and show that, to a good approximation, the phenomenology of the lighter Higgs state can be described by its mass and three couplings: those to massive gauge bosons and to top and bottom quarks. We perform a fit of these couplings using the latest LHC data on the production and decay rates of the light $h$ boson and combine it with the limits from the negative search of the heavier $H,A$ and $H^\pm$ states, taking into account the current uncertainties.Comment: 1+12 pages, pdflatex, 7 figure

Spatial fluctuations of a surviving particle in the trapping reaction

We consider the trapping reaction, $A+B\to B$, where $A$ and $B$ particles have a diffusive dynamics characterized by diffusion constants $D_A$ and $D_B$. The interaction with $B$ particles can be formally incorporated in an effective dynamics for one $A$ particle as was recently shown by Bray {\it et al}. [Phys. Rev. E {\bf 67}, 060102 (2003)]. We use this method to compute, in space dimension $d=1$, the asymptotic behaviour of the spatial fluctuation, $^{1/2}$, for a surviving $A$ particle in the perturbative regime, $D_A/D_B\ll 1$, for the case of an initially uniform distribution of $B$ particles. We show that, for $t\gg 1$, $^{1/2} \propto t^{\phi}$ with $\phi=1/4$. By contrast, the fluctuations of paths constrained to return to their starting point at time $t$ grow with the larger exponent 1/3. Numerical tests are consistent with these predictions.Comment: 10 pages, 5 figure

Force transmission in a packing of pentagonal particles

We perform a detailed analysis of the contact force network in a dense confined packing of pentagonal particles simulated by means of the contact dynamics method. The effect of particle shape is evidenced by comparing the data from pentagon packing and from a packing with identical characteristics except for the circular shape of the particles. A counterintuitive finding of this work is that, under steady shearing, the pentagon packing develops a lower structural anisotropy than the disk packing. We show that this weakness is compensated by a higher force anisotropy, leading to enhanced shear strength of the pentagon packing. We revisit "strong" and "weak" force networks in the pentagon packing, but our simulation data provide also evidence for a large class of "very weak" forces carried mainly by vertex-to-edge contacts. The strong force chains are mostly composed of edge-to-edge contacts with a marked zig-zag aspect and a decreasing exponential probability distribution as in a disk packing

Dynamic Provenance for SPARQL Update

While the Semantic Web currently can exhibit provenance information by using the W3C PROV standards, there is a "missing link" in connecting PROV to storing and querying for dynamic changes to RDF graphs using SPARQL. Solving this problem would be required for such clear use-cases as the creation of version control systems for RDF. While some provenance models and annotation techniques for storing and querying provenance data originally developed with databases or workflows in mind transfer readily to RDF and SPARQL, these techniques do not readily adapt to describing changes in dynamic RDF datasets over time. In this paper we explore how to adapt the dynamic copy-paste provenance model of Buneman et al. [2] to RDF datasets that change over time in response to SPARQL updates, how to represent the resulting provenance records themselves as RDF in a manner compatible with W3C PROV, and how the provenance information can be defined by reinterpreting SPARQL updates. The primary contribution of this paper is a semantic framework that enables the semantics of SPARQL Update to be used as the basis for a 'cut-and-paste' provenance model in a principled manner.Comment: Pre-publication version of ISWC 2014 pape

Pore Stabilization in Cohesive Granular Systems

Cohesive powders tend to form porous aggregates which can be compacted by applying an external pressure. This process is modelled using the Contact Dynamics method supplemented with a cohesion law and rolling friction. Starting with ballistic deposits of varying density, we investigate how the porosity of the compacted sample depends on the cohesion strength and the friction coefficients. This allows to explain different pore stabilization mechanisms. The final porosity depends on the cohesion force scaled by the external pressure and on the lateral distance between branches of the ballistic deposit r_capt. Even if cohesion is switched off, pores can be stabilized by Coulomb friction alone. This effect is weak for round particles, as long as the friction coefficient is smaller than 1. However, for nonspherical particles the effect is much stronger.Comment: 10 pages, 15 figure

Generalized Forward-Backward Splitting

This paper introduces the generalized forward-backward splitting algorithm for minimizing convex functions of the form $F + \sum_{i=1}^n G_i$, where $F$ has a Lipschitz-continuous gradient and the $G_i$'s are simple in the sense that their Moreau proximity operators are easy to compute. While the forward-backward algorithm cannot deal with more than $n = 1$ non-smooth function, our method generalizes it to the case of arbitrary $n$. Our method makes an explicit use of the regularity of $F$ in the forward step, and the proximity operators of the $G_i$'s are applied in parallel in the backward step. This allows the generalized forward backward to efficiently address an important class of convex problems. We prove its convergence in infinite dimension, and its robustness to errors on the computation of the proximity operators and of the gradient of $F$. Examples on inverse problems in imaging demonstrate the advantage of the proposed methods in comparison to other splitting algorithms.Comment: 24 pages, 4 figure

An interactive medical image segmentation system based on the optimal management of regions of interest using topological medical knowledge

This paper presents an original interactive system for efficient medical image segmentation in computer aided diagnosis. The main originality concerns the method used to manage, according to an a priori topological-based structural model, regions of interest (ROIs) within which computations can be constrained. The goal is then to avoid the processing of irrelevant image points, therefore improving and accelerating segmentations. In the case of a hierarchical modeling procedure, our ROI management method enables, for delineating a given medical structure, to optimally determine image points of interest by taking previously segmented structures into account. We propose a mathematical formulation of the method as well as a possible implementation within an interactive system. We also detail an experience report focussing on the segmentation of several abdominal structures from a CT image. It illustrates the behavior and the potential of our method